Wpływ rozstawu względnego (a/d) na opór aerodynamiczny pionowej drogi kablowej

Wprowadzenie

W konstrukcjach telekomunikacyjnych często występuje konieczność prowadzenia wielu kabli wzdłuż wieży lub masztu. Kable te – ze względu na długość i lokalizację – podlegają znacznym obciążeniom wiatrowym. Wartość tych obciążeń zależy m.in. od średnicy kabla, prędkości wiatru oraz wzajemnego położenia kabli.

W sytuacji, gdy kilka przewodów znajduje się blisko siebie, mogą występować istotne interakcje aerodynamiczne – takie jak ekranowanie lub wzajemne wzmacnianie przepływu – które wpływają na całkowity opór aerodynamiczny układu.

Wprowadzenie

Celem niniejszej analizy jest ilościowa ocena wpływu wzajemnego położenia kabli w wiązce na całkowity współczynnik oporu aerodynamicznego. Skoncentrowano się na typowej konfiguracji napowietrznej trasy kablowej mocowanej do masztu stalowego, w której przewody prowadzone są równolegle w poziomej płaszczyźnie.

Aby umożliwić analizę szeregu wariantów geometrycznych przy zachowaniu rozsądnego czasu obliczeń, przyjęto model dwuwymiarowy (2D), w którym przekrój poprzeczny układu traktowany jest jako nieskończenie długi w kierunku osiowym. Takie podejście stanowi kompromis między dokładnością a nakładem obliczeniowym – pozwala uchwycić kluczowe zjawiska aerodynamiczne związane z oddziaływaniem sąsiednich walców, przy znacznie mniejszym koszcie czasowym i sprzętowym w porównaniu do pełnej analizy 3D.

Wyniki badania mogą być pomocne w projektowaniu układów kablowych na wieżach i masztach, szczególnie tam, gdzie istotne jest ograniczenie obciążeń wiatrowych.

Analiza aerodynamiczna

W ramach oceny wpływu wiatru na elementy konstrukcji opracowano uproszczony model układu 10 równoległych walców o jednakowej średnicy d = 15 mm, reprezentujących kable rozmieszczone w jednej płaszczyźnie. Geometrię analizowano w 2D jako przekrój poprzeczny przez układ – umożliwia to dokładne odwzorowanie zjawisk aerodynamicznych.

Rozstaw między kablami oznaczono jako a, a podstawą analizy był bezwymiarowy stosunek a/d, co pozwala na szersze porównywanie wyników.

Podejście numeryczne

Rys. 1: Schematyczny przekrój poprzeczny analizowanego układu walców w 2D, ze wskazaniem średnicy d, rozstawu a oraz kierunku przepływu powietrza.

Podejście numeryczne

Zastosowano metodę CFD (Computational Fluid Dynamics), polegającą na symulacji przepływu powietrza wokół analizowanego układu z wykorzystaniem podziału domeny obliczeniowej na dyskretną siatkę. Pozwala to na odwzorowanie lokalnych zmian ciśnienia, prędkości i sił działających na obiekt, a tym samym – na wyznaczenie współczynników aerodynamicznych.

Do przeprowadzenia analizy wykorzystano oprogramowanie FloEFD, zintegrowane z systemem CAD Solid Edge. Program korzysta z modelu turbulencji k–ε RNG, który dobrze odwzorowuje przepływy z dużymi gradientami prędkości – takimi jak separacje strugi czy wiry za przeszkodami.

W strefach przyściennych zastosowano funkcje ścienne do odwzorowania gradientów prędkości i efektów lepkości blisko powierzchni kabli. Symulacje przeprowadzono w trybie ustalonym (steady-state), bez analizy dynamiki strukturalnej (brak sprzężenia płyn–konstrukcja, FSI). Oznacza to, że nie uwzględniono zjawisk rezonansowych ani zmiennych w czasie sił wymuszających. W praktyce – dla rzeczywistych konstrukcji – należy stosować odpowiedni współczynnik dynamiczny cscd zgodnie z PN-EN 1993-3-1, jednak nie był on przedmiotem niniejszej analizy.

Na Rysunku 1 przedstawiono schematyczny przekrój 2D analizowanego układu walców wraz z oznaczeniem podstawowych wymiarów (średnicy d oraz odstępu a) oraz kierunku przepływu powietrza.

Parametry środowiskowe

  • Przyjęto następujące parametry środowiskowe:
  • Gęstość powietrza – 1,25 [kg/m3]
  • Ciśnienie podstawowe – 1013 [hPa]
  • Temperatura otoczenia – 20 [°C]
  • Lepkość kinematyczna powietrza – 15 * 10-5 [m2/s]
  • Grawitacja – 9.81 [m/s2]
  • Analizę przeprowadzoną jako przepływ ustalony, czyli niezmienny w czasie.

Wartości te odpowiadają standardowym warunkom referencyjnym dla analiz inżynierskich.

Analizy przeprowadzono przy zadanej prędkości wiatru 35 m/s. Wartość ta została wybrana na podstawie doświadczenia inżynierskiego oraz dostępnych danych referencyjnych i jest typowa dla obliczeń aerodynamicznego dla tego typu obiektów.

Wyznaczanie współczynnika oporu

Dla każdego przypadku (tj. danej wartości a/d) obliczano siłę oporu działającą na zestaw kabli, a następnie określano średni współczynnik oporu aerodynamicznego Cf z relacji:

F– siła oporu działająca na zestaw kabli [N]
A – pole powierzchni czołowej zestawu kabli [m2]
ρ – gęstość powietrza [kg/m3]

Siatka obliczeniowa

Dla przygotowanego modelu CAD utworzono siatkę obliczeniową o strukturze sześciennej, zróżnicowaną pod względem gęstości w zależności od lokalnych wymagań analizy. Zagęszczenie siatki (ang. refinement) stosowano automatycznie w całej domenie – zarówno w obszarach przy powierzchni kabli, jak i w ich otoczeniu, zwłaszcza za przeszkodami, gdzie występują istotne zmiany przepływu. Dodatkowo, w kluczowych rejonach blisko kabli zastosowano ręczne dogęszczenia siatki dla poprawy dokładności wyników.

Na rysunkach w dalszej części przedstawiono przykładowe przekroje siatki obliczeniowej – zarówno ogólne widoki domeny, jak i lokalne zbliżenia na obszary o zwiększonej gęstości siatki.

Siatka obliczeniowa

Rys. 2: Bazowa siatka obliczeniowa – przekrój poziomy

Rys. 3: Bazowa siatka obliczeniowa – przekrój poziomy – widok szczegółowy

Rys. 3: Bazowa siatka obliczeniowa – przekrój poziomy – widok szczegółowy

Rys. 4: Bazowa siatka obliczeniowa – przekrój poziomy – widok pojedynczego kabla.

Rys. 4: Bazowa siatka obliczeniowa – przekrój poziomy – widok pojedynczego kabla.

Prezentacja wyników

W ramach analizy CFD wyznaczono wartości siły aerodynamicznej działającej na model dla różnych wartości odległości a przy stałej wartości wymiaru d=15 mm. Na tej podstawie obliczono bezwymiarowy stosunek a/d oraz odpowiadające mu współczynniki aerodynamiczne Cf. Wyniki te zestawiono w tabeli.

a [mm]F [N]a/d [-]Cd [-]
1.002.800.071.36
1.102.730.071.33
1.212.640.081.28
1.322.670.091.32
1.442.550.101.26
1.562.640.111.31
1.692.450.121.22
1.832.640.131.32
1.982.580.141.31
2.142.750.141.33
2.362.680.161.30
2.592.710.171.31
2.852.630.191.27
3.142.550.211.23
3.452.590.231.25
3.802.610.261.25
4.182.490.281.21
4.592.400.311.16
5.052.210.341.06
5.552.140.371.02
6.092.060.410.97
6.701.960.450.91
7.371.840.500.85
8.101.750.550.80
8.911.630.610.73
9.801.510.670.66
10.781.420.740.62
11.861.340.810.58
13.051.290.890.56
14.351.230.980.53
15.781.171.080.51
17.351.141.190.51
19.061.101.310.48
20.911.061.440.45
22.931.011.580.43
25.120.991.730.42
27.500.961.900.41
30.070.942.070.40
32.870.912.260.39
35.900.892.470.38
39.180.862.690.37
42.730.832.930.36
46.570.823.190.35

Tabela 1. Zestawienie wyników analizy CFD

Na podstawie danych przedstawionych w tabeli sporządzono wykres zależności współczynnika oporu Cf od bezwymiarowego parametru a/d, definiowanego jako stosunek odległości między osiami sąsiednich przewodów do ich średnicy. Wykres ten umożliwia graficzne przedstawienie wpływu zmiany rozstawu kabli na charakterystykę aerodynamiczną wiązki.

Zastosowanie parametru a/d pozwala na uogólnienie wyników – niezależnie od konkretnej średnicy przewodów – co ułatwia ich zastosowanie w różnych przypadkach projektowych. Taka forma prezentacji umożliwia szybkie wychwycenie trendów, takich jak zmniejszanie się lub stabilizacja wartości Cf przy zwiększającym się rozstawie, a także identyfikację zakresów, w których występuje silna interakcja aerodynamiczna między przewodami. W celu uchwycenia ogólnej tendencji zmian współczynnika Cf wraz ze wzrostem parametru a/d, na wykresie naniesiono linię trendu o charakterze logarytmicznym, dopasowaną do wyników obliczeń. Jej równanie ma postać: Cf= -0,252*ln(a/d) + 0,788

Współczynnik determinacji wynosi R2=0,954, co świadczy o bardzo dobrej zgodności modelu z uzyskanymi danymi. Równanie oraz wartość R2 umieszczono również bezpośrednio na wykresie.

Współczynnik determinacji wynosi R2=0,954

Dla zobrazowania charakteru przepływu pokazano pola prędkości w płaszczyźnie pionowej dla trzech wybranych przypadków geometrycznych: a/d<1; 1< a/d <3; oraz a/d >3.

Zestawienie to pozwala prześledzić, jak zmienia się rozkład prędkości w przestrzeni między kablami w zależności od ich wzajemnego położenia.

  • Dla najmniejszych odległości (a/d < 1) przepływ jest silnie zaburzony, a pomiędzy kablami tworzą się wyraźne strefy spowolnionego przepływu.
  • W zakresie pośrednim (1 < a/d < 3) obserwuje się stopniową stabilizację strug i zwiększenie prędkości w przestrzeniach międzyprzewodowych.
  • Dla największych odległości (a/d > 3) profile przepływu stają się bardziej symetryczne i zbliżone do przypadku pojedynczego cylindra.
Rys. 3: Bazowa siatka obliczeniowa – przekrój poziomy – widok szczegółowy

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅0,10

: Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅0,25

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅0,25

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅0,40

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅0,40

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅0,60

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅0,60

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅0,80

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅0,80

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅1,30

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅1,30

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅2,50

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅2,50

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅3,70

Przykładowe pole prędkości wokół układu dla wartości parametru a/d ≅3,70

Podsumowanie

W analizie CFD zbadano wpływ zmiennej odległości a (przy stałym d = 15 mm) na siłę aerodynamiczną działającą na model. Wyniki przedstawiono w formie tabeli oraz wykresu zależności współczynnika Cf od stosunku a/d.

Analiza wykazała różny charakter oddziaływań aerodynamicznych w zależności od stosunku a/d, który można podzielić na trzy zakresy:

  • Dla a/d<1 obserwuje się wyraźne oddziaływanie aerodynamiczne między kablami, co prowadzi do podwyższonych wartości Cf.
  • Dla 1<a/d<3 wpływ sąsiednich kabli jest zauważalny lecz przestaje być znaczący. Przy a/d>3 przekraczającym 3 przepływy wokół kabli zaczynają się zachowywać niezależnie, a wzajemne oddziaływanie staje się pomijalne.
  • Pola prędkości potwierdzają te obserwacje, pokazując zmienność rozkładu przepływu w zależności od rozstawu elementów.

Uwagi końcowe

Uzyskane wyniki analizy CFD mogą być wykorzystane przy projektowaniu zespołów przewodów lub kabli narażonych na działanie wiatru. Dla wartości a/d powyżej 3 przewody można traktować jako osobne walce o niezależnym oddziaływaniu aerodynamicznym. Natomiast dla niższych wartości a/d konieczne jest stosowanie współczynnika bezpieczeństwa lub przeprowadzenie szczegółowej analizy CFD, aby uwzględnić wzajemne oddziaływania kabli.
Należy jednak pamiętać, że analiza została przeprowadzona w warunkach ustalonych i idealizowanych, bez uwzględnienia wpływu pełnoskalowej turbulencji czy efektów drgań własnych kabli, które w rzeczywistych warunkach eksploatacji mogą mieć znaczenie.

Skontaktuj się z nami

Zadzwoń do nas, aby uzyskać więcej informacji lub omówić szczegóły

+48 608 539 911

    Informacja UE

    Carbonpro sp. z o.o. realizuje projekt "Udział w programie akceleracyjnym Concordia Design Accelerator w ramach ścieżki B2B/B2A/Inwestorska" w ramach projektu pn. „Concordia Design Accelerator – We turn ideas into human profit”, w ramach programu Fundusze Europejskie dla Nowoczesnej Gospodarki działanie Startup Booster Poland – Smart Up.

    Cel projektu: Wykonanie walidacji wyników uzyskanych metodą CFD za pomocą badań eksperymentalnych w tunelu aerodynamicznym

    Dofinansowanie: 379.072,36 zł