Leszek Chodor, 26 września 2014
15 czerwca 2025 naprawa awarii i scalenie kilku odrębnych artykułów
W przypadku nieczytelnych treści, proszę powiadomić: leszek@chodor.co
W ciągu ostatnich 24 godzin z artykułu korzystało 13 Czytelników
Artykuł jest naprawiany po poważnej awarii portalu
Zagadnienie skręcania prętów jest ważnym problemem w projektowaniu konstrukcji budowlanych. Współczesne analizy wytrzymałości i stateczności systemów prowadzone są na modelach złożonych z prętów uogólnionych w stosunku do klasycznego modelu Bernoulliego. Pręt uogólniony ma dodatkowym stopień swobody węzła – paczenie (deplanację) przekroju. Uogólniona teoria belkowa Sapountzakis (2013) [1] zbudowana na dorobku kilku pokoleń inżynierów, a zapoczątkowana pracą Marguerre (1940) [2] umożliwia włączenie w proces obliczeniowy analizy wyboczenia giętno-skrętnego (w tym zwichrzenia) prętów. Analiza konstrukcji prętowych z zastosowaniem teorii zginania i skręcania skrępowanego prętów cienkościennych praktycznie w całości zastąpiła tradycyjną analizę z użyciem modelu pręta zwartego i jest stosowana w zasadzie we wszystkich nowoczesnych programach obliczeniowych (program Consteel i szereg innych).
Naprężenia, odkształcenia i przemieszczenia skręcanego pręta istotnie zależą od kształtu przekroju oraz od sposobu podparcia (warunków brzegowych ) końcowych (poprzecznych) ścianek i pośrednich przekrojów, a także od sposobu obciążenie pręta. Ścisłe rozwiązanie zagadnienia skręcania możliwe jest tylko w nielicznych przypadkach, a w zasadzie tylko dla pręta o przekroju okrągłym, ze swobodnymi ściankami poprzecznymi, obciążonymi równomiernymi obciążeniami skręcającymi o takiej samej wartości ale o przeciwnych zwrotach na obu ściankach poprzecznych. W każdym innym przypadku mamy do czynienia z rozwiązaniem przybliżonym, ale akceptowanym w praktycznym wymiarowaniu prętów
Teoria prętów cienkościennych oryginalnie sformułowana przez Własowa (1959) [3] znajduje zastosowanie w analizie nie tylko konstrukcji stalowych, ale również w analizie skręcania prętów i obiektów żelbetowych, mostów skrzynkowych lub budynków wysokich i wszędzie tam, gdzie wystąpi nieswobodne (skrępowane) skręcania prętów. W praktyce skręcanie nieswobodne występuje w każdym pręcie konstrukcji budowlanych, a skręcanie swobodne wystąpi tylko w nielicznych przypadkach, na przykład skręcania walów na łożyskach.
W prętach skręcanych swobodnie wytężenie przekrojów spowodowane jest tylko naprężeniami stycznymi, a podczas skręcania skrępowanego wystąpią również naprężenia normalne, które są spowodowane przez powstanie nieznanej w klasycznej teorii zginania prętów siłę przekrojową – bimoment Bω, stowarzyszony z deplanacją (spaczeniem) przekroju. Powstaną również dodatkowe naprężenia styczne wywołane przez działanie momentu giętno-skrętnego Mω. Rozkład naprężeń ścinających po grubości ścianki zależy od tego, czy przekrój jest otwarty, czy też zamknięty.
Współczesne programy komputerowe, szczególnie ConSteel [4], [5] pozwalają na projektowanie ram złożonych z prętów stalowych z uwzględnieniem niestateczności ogólnej dowolnego typu (wyboczenie, zwichrzenie, wyboczenie giętno-skrętne) bez potrzeby wykorzystywania skomplikowanych formuł analitycznych i normowych. W 2002 roku opublikowano do powszechnego użytku inżynierów program LTBeam [6], który umożliwiał szybkie i łatwe obliczenie momentu krytycznego belek jedno- i wieloprzęsłowych o dowolnych przekrojach bi- i mono-symetrycznych ze stężeniami bocznymi niepodatnymi lub sprężystymi założonymi na dowolnej wysokości przekroju. Od wersji LTBeamN 2.02 za pomocą programu można obliczyć siły krytyczne dla belki zginanej momentem M i ściskanej siłą N. Program został opracowany przez CTICM (Centre Technique Industriel de la Construction Metallique – France) w ramach europejskiego projektu ECSC(European Community for Steel and Coal) – Project No 7210-PR183 : „Lateral torsional buckling of steel and composite beams” – 1999-2002).
Tablice projektowe
Charakterystyki przekrojów skręcanych
Tab.1 Charakterystyki geometryczne skrętne wybranych przekrojów zwartych stosowanych w budownictwie
W przypadku swobodnego skręcania jednostkowy kąt skręcenia Θ wyznacza się z zależności (???).
Tab.2 Współczynniki wskaźnika wytrzymałości (???) i momentu bezwładności (???)na skręcanie dla różnych stosunków boku prostokąta (opracowano na podstawie [7] )
Tab.3 Charakterystyki geometryczne wybranych, zamkniętych przekrojów cienkościennych [8]
A=b⋅t2⋅(t1⋅cos2α+t2⋅cosα)⋅tgα4⋅(t1⋅cosα+t2)⋅(t1+t2⋅cosαB=b38⋅t1⋅t2⋅sinαt1+t2⋅cosα,C=2⋅b2⋅h2⋅t1⋅t2b⋅t2+h⋅t1,,D=b2⋅h276⋅(b⋅t2–h⋅t1)2(b⋅t2+h⋅t1)2⋅(b⋅t1+h⋅t2)
Tab.4 Charakterystyki geometryczne wybranych, zamkniętych przekrojów cienkościennych – wg [8]
Sposób wyznaczania współrzędnych wycinkowych dla dwuteownika, ceownika i zetownika pokazano w przykładzie 1
Tab.5 Bimomenty Bω, momenty giętno-skrętne Mω, momenty swobodnego skręcania Ms oraz kąty skręcenia φ dla prętów cienkościennych [8], tab 1-2
Tab.6 Współczynniki korekcyjne sztywności k na skręcanie kształtowników stalowych

2) Wyniki badań Foplle przedstawiono w poradniku Birger i in (1998) [7], 3) Niezgodziński (2004) [9]
Skręcanie czyste
Skręcanie czyste realizuje się podczas działania na ścianki poprzeczne (czołowe) pręta pryzmatycznego specyficznie dobranego obciążenia o gęstości qv=[0,qvy,qvz] (rys. 1a). Obciążenie qv można statycznie zredukować do pary momentów Ms działającej w płaszczyźnie przekroju poprzecznego początkowego A i końcowego pręta (w miejscu utwierdzenia S). Obciążenie zewnętrzne parą momentów Ms wywołuje moment skręcający Mx w przekrojach pręta.
Mx jest całkowitym momentem skręcającym przekrój wokół osi x, wyznaczonym w drodze statycznej redukcji obciążenia zewnętrznego Ms do rozpatrywanego przekroju pręta w którym jest zrównoważony przez składowe siły przekrojowe: moment czystego skręcania (Saint Venanta) Mv i moment giętno-skrętny Mω:
Mx=Mv+Mω
Moment Mv powstaje przy skręcaniu swobodnym, a przy skręcaniu nieswobodnym dodatkowo działa moment Mω
Można pokazać [10], że gęstość obciążenia qv ścianek czołowych pręta, prowadząca do czystego skręcania musi mieć rozkład
qv={qvx=0qvy=±Θ⋅G⋅(∂φ∂y–z)qvz=±Θ⋅G⋅(∂φ∂z+y)
gdzie Θ jest parametrem, a funkcja φ=φ(y,z) jest funkcją harmoniczną, czyli spełnia warunek zagadnienia brzegowego Neumana
∇φ2=0
i tak dobraną, że spełnia warunki równowagi (statyczne warunki równowagi) na pobocznicy oraz kinematyczne warunki brzegowe , a w tym przypadku zeruje się w punkcie utwierdzenia S (rys. 1a) .

Rys.1. Skręcanie: a) czyste, b) proste, c) swobodne, d) skrępowane – nieswobodne (rysunek złożony z ilustracji z pracy [10] )
Skręcanie proste
Skręcanie proste (rys. 1b) dotyczy problemu technicznego skręcania pręta obciążonego na powierzchniach czołowych parą momentów skręcających Ms (rys. 1b), czyli dowolnym (niekoniecznie spełniającym warunek (2) )rozkładem obciążeń qv, które jest statycznie równoważne momentowi Ms,
Warunek równoważności obciążenia momentami skupionymi Ms oraz obciążenia realizującego czyste skręcanie (2) prowadzi do zależności
Ms=∬A(qvzy−qvyzdA=∬A[ΘG(∂φ∂z+y)y–ΘG(∂φ∂y−z)z]dA=ΘG∬A(∂φ∂zy−∂φ∂yz+y2+z2)=ΘG⋅Iv
gdzie wprowadzono oznaczenie definiowane jako moment bezwładności przekroju na skręcanie (Saint Vennanta).
Ivdef=ΘG∬A(∂φ∂zy−∂φ∂yz+y2+z2)
Skręcanie swobodne
Przy skręcanie swobodnym pręta każdy przekrój poprzeczny może się swobodnie odkształcać w kierunku osi pręta, co wywołuje paczenie deplanację przekroju płaskiego przed obciążeniem (rys. 1c) . W praktyce oznacza to, że w przekroju występują tylko naprężenia styczne, a deplanacja na rozkład naprężeń. Taki przypadek zachodzi, gdy na końcach pręt jest obciążony równymi i przeciwnie skierowanymi momentami skręcającymi Ms (rys. 1b).
Skręcanie nieswobodne
Skręcanie nieswobodne (skrępowane) to takie, w którym odebrano swobodę deplanacji przekroju (rys. 1c) i w rezultacie w pręcie oprócz naprężeń stycznych pojawiają się naprężenia normalne składające się na siłę przekrojową zwaną bimomentem.
W każdym przypadku skręcania należy zwrócić uwagę na swobodę deplanacji przekroju , bo przy jej ograniczeniu, obok naprężeń od swobodnego skręcania (Saint Venanta) szacować naprężenia od momentu giętno-skrętnego, powstających przy nieswobodnym skręcaniu.
Deplanacja przekroju pręta skręcanego
Zasada sztywnego konturu i deplanacja przekroju
Skręcanie pręta znamienne jest tym, że pręt ulega deplanacji (spaczeniu) , to znaczy przekrój płaski przed odkształceniem nie zachowuje płaskości . Zjawisko deplanacji obserwuje się w pręcie skręcanym swobodnie, co pokazano na rys. 2. , ale także w pręcie skręcanym nieswobodnie , w którym ograniczenie swobody deplanacji wywołuje nowe siły przekrojowe.
W rezultacie podczas skręcania nie jest zachowana zasada płaskich przekrojów -fundamentalna zasada stosowana w technicznej teorii zginania,. Założenie teorii zginania ” odcinek prosty i prostopadły do osi pręta przed deformacją pozostaje prosty i prostopadły po deformacji, a pręt nie doznaje odkształceń (od zginania) na linii osi przekrojów” zastępuje zasada sztywnego konturu „linie konturowa przekrojów poprzecznych pręta nie ulega odkształceniu (od skręcania), ale punktu zawarte wewnątrz konturu przemieszczają się wzdłuż osi pręta”.
W sformalizowanym podejściu do analizy skręcania prętów zasada sztywnego konturu oraz deplanacja przekroju jest przedstawiano jako wniosek z podstawowego równania skręcania swobodnego (???).
Kratowy model skręcania
Analizę napreżeń w skręcanym pręcie można przeprowadzić na kratowym modelu skręcania, pokazanym na rys.2 Model taki rozpatruje się przy analizie skręcania przekroju żelbetowego. Na rys.2 pokazano ten model z oznaczeniem kąta nachylenia krzyżulca γ, podczas gdy w analizie konstrukcji żelbetowych kąt ten oznacza się jako Θ i utożsamia a kątem nachylenia krzyżulca betonowego tworzącego się w pręcie betonowym w granicach 1,0≤ctgΘ≤2.
Na kratowym modelu skręcania rys.2. możemy zaobserwować, że w przypadku odebrani swobody skręcania na dalszym końcu pręta w zbrojeniu podłużnym pojawia się rozciąganie w pasach górnych i dolnych, co jest objawem deplanacji przekroju (porównaj z rys.3)
Skręcanie swobodne
Ścisłe rozwiązanie zadania swobodnego skręcania w przypadku dowolnego materiału i kształtu przekroju pręta jest zaskakująco złożone, mimo, że w podstawowych podręcznika tematu zagadnienie swobodnego skręcania jest traktowane jako proste zadanie.
W 1853 roku Saint-Venant (p. np [11]) przedstawił metodę półodwrotną rozwiązania zagadnienia swobodnego skręcania izotropowych sprężystych prętów pryzmatycznych. Zaproponował metodę, w której założył postać przemieszczeń pręta z możliwością spaczenia przekroju (identycznego dla wszystkich przekrojów poprzecznych) i wyprowadził, konsekwentnie w ramach liniowej sprężystości, równania opisujące zadanie swobodnego skręcania prętów pryzmatycznych. Rozkład obciążenia zewnętrznego na końcach pręta jest identyczny z rozkładem naprężeń w każdym przekroju. Wykazał, że sposób przyłożenia statycznie równoważnego obciążenia nie jest istotny w przekrojach dostatecznie odległych od jego końców. Uogólnienie tej idei, którą nazywa się zasadą Saint-Venanta stosowaną w innych przypadkach wytrzymałościowych ciała stałego.
Na rys.2 pokazano postać odkształconą pręta o zwartym przekroju kwadratowym poddanego swobodnemu skręcaniu parą momentów skręcających Ms o przeciwnych zwrotach przyłożonych do powierzchni czołowej.
Saint Venant na podstawie wyników eksperymentów założył następujące pole przemieszczeń przekroju pręta skręcanego swobodnie ( u- przemieszczenie w kierunku osi pręta x; v- przemieszczenie wzdłuż osi y ; w-przemieszczenie pionowe wzdłuż osi z):
u=Θ⋅Φ(y,z);v=–Θ⋅x⋅z;w=Θ⋅x⋅y
gdzie Φ(y,z) – funkcja paczenia zależna od rodzaju przekroju. Dla przekroju okrągłego (walca) Φ(y,z)=1
Jednostkowy kąt skręcenia Θ=dφdx wyznacza się z zależności
Θ=MvGIv
gdzie Mv moment czystego skręcania Saint Venanta. Dla schematu z rys. 2 moment skręcający Saint Venanta przekrój pręta jest stały na długości pręta i równy obciążeniu zewnętrznym momentem skręcająMS : Mv=Ms
Pręt o przekroju zwartym
Z punktu widzenia praktyki istotne znaczenie ma problem czystego skręcania pręta o przekroju okrągłym oraz prostokątnym (rys.2).

Rys.4. Naprężenia w przekroju skręcanym : s) okrągłym, b) prostokątnym [10]
Przekrój okrągły
W przypadku przekroju kołowego rozwiązaniem zagadnienia Neumanna jest funkcja
φ(y,z)=0
W konsekwencji moment bezwładności przy skręcaniu (5) wynosi
Iv=∬Ar2dA=πR42=I0
Maksymalne naprężenia styczne (na obwodzie koła) wynoszą
Przekrój prostokątny
W przypadku przekroju prostokątnego o wymiarach bxh (h-wysokość) (rys. 2b) funkcja paczenia może być przedstawiona w postaci szeregu [10]:
φ(y,z)=y⋅z−∞∑i=0Bnkncosh(knh/2⋅sin(kny)cosh(knz)
gdzie
kn=(2n+1)πb;Bn=(−1)n8b(2b+1)2π2
W konsekwencji moment bezwładności skręcania (5) możemy zapisać w postaci
Iv=[13–66π5⋅(h/b))−1∞∑n=0tgh(knh/2)(2n+1)5]b3h=β(h/b)⋅b3⋅h
Wyrażenie na maksymalne naprężenia styczne można zapisać w postaci:
maxτxz=Mvα(h/b)⋅b2h
skąd wskaźnik wytrzymałości na skręcania
Wv=α(h/b)⋅b2⋅h
Współczynniki α i β są zależne jedynie od stosunku boków prostokąta i zestawiono je w tab.1
W literaturze podawane są jeszcze inne przybliżenia na współczynniki α i β. W pracy [7] i w licznych podręcznikach teorii sprężystości lub wytrzymałości materiałó w zamieszcza się przybliżone formuły:
β≈13[1–0,63(h/b)+0,052(h/b)5];α≈[1–8π2+∞∑n=1,3,…1[n2cosh(nπh)/(2b)]]
Przybliżona formuła Saint-Venanta dla sztywności skrętnej
W praktyce inżynierskiej często potrzebujemy szybkiego szacunku sztywności na skręcanie przekroju o złożonym kształcie w celu podjęcia decyzji o wprowadzeniu przekroju do projektu konstrukcji, W tym przypadku można stosować formułe podaną przez Sait-Venanta w postaci
Iv≈A440⋅I0
gdzie:
A- pole przekroju pręta zwartego lub pole zawarte wewnątrz przekroju zamkniętego.
I0=Iy+Iz – osiowy moment bezwładności równy sumie momentów bezwładności względem osi głównych przekroju.
Formuła (17) powstała po analizie wyrażenia na sztywność elipsy:
Iv=πa3b3a2+b2
gdzie a, b – półosie elipsy, która ma pole przekroju A=πab oraz Io=πab4(a2+b2).
Po uwzględnieniu tego, że 4π2≈40 z (18) uzyskujemy (17) .
Formuła (17) często jest stosowana do przybliżonego wyznaczenia sztywności na skręcanie prętów o dowolnym przekroju. Należy jednak wziąć pod uwagę, że takie przybliżenie może prowadzić do zbyt dużych błędów. Na przykład w przypadku pręta okrągłego z nacięciem błąd (zawyżenie sztywności) może osiągnąć 80% [7]
Przekrój cienkościenny
Definicja pręta cienkościennego
Pręt cienkościenny to cylindryczna lub pryzmatyczna powłoka, której trzy miarodajne wymiary są różnych rzędów: długość jest o rząd większa od szerokości (lub wysokości) przekroju, która z kolei jest o rząd większa od grubości ścianki powłoki [12]. Jak pokazuje doświadczenie, a także analiza statyczno-wytrzymałościowa oparta na hipotezach Wlasowa, pręt można uważać za cienkościenny, jeśli spełnione są warunki [3], [13], [14]:
t≤lNth;l≤LNb
gdzie: Nth=Nb=10,
t – grubość ścianki, l – długość ścianki , L – długość pręta.
Drugi warunek wynika z ogólnej definicji pręta (również o przekroju zwartym). Elementy o mniejszej smukłości w zasadzie powinny być analizowane jako tarcze. W technice często ten warunek osłabia się i wymaga Nb=8, a w konstrukcjach żelbetowych nawet Nb=4.
Przy zmniejszaniu smukłości ścianek l/t lub pręta L/l dokładność rozwiązania zadania jest coraz gorsza.
Pręty cienkościenne dzieli się z uwagi na geometrię przekroju na:
- otwarte, to znaczy takie, w których linia środkowa nie tworzy obwodów, czyli nigdzie nie zamyka się (np. rys. 5, 6),
- zamknięte, toznaczy takie, w których ścianka tworzy obwód zamknięty , rurę, komorę (np. rys. 7) . Przekroje zamknięte mogą być jedno- lub wielokomorowe.
- mieszane, które zawierają części otwarte i części zamknięte
- quasi-zamknięte.(rys. 9 ) i rozdział pręty wielogałęziowe oraz kratownice wielopasowe.
Przekrój otwarty
Na rys. 5 przedstawiono często stosowane w konstrukcjach stalowych

Rys.5 Przykłady przekrojów cienkościennych otwartych stosowane w konstrukcjach stalowych : a) dwuteownik walcowany, b) ceownik walcowany, c) teownik walcowany (1/2 I), d) ceownik gięty na zimno z zagiętymi pasami (C+), e) kątownik gięty na zimno, f) przekrój spawany złożony gwiaździsty , g) przekrój prosty gwiaździsty; C- środek ciężkości, S- środek zginania ( ścinania) [15]
Rozważmy przekrój poprzeczny pokazany na rys.6,, który składa się z n elementów prostokątnych o długości li o grubości ti elementu i-tego, spełniających warunek (???)..

Rys. 6 Przekrój otwarty cienkościenny złożony ze ścianek prostych [16]
Przyjmiemy następujące założenia upraszczające dotyczące otwartych profili cienkościennych [16]:
#1 jednostkowy kat skręcenia każdego elementu prostokątnego przekroju poprzecznego jest jednakowy,
#2 suma momentów skręcających poszczególne elementy Mυi jest równa momentowi skręcającemu przyłożonemu do całego profilu Mυ,
Jednostkowy kąt skręcenie występujący w (2) do (5) zgodnie założeniem #1 możemy napisać w postaci
Θ=MvGIv=MviGIvi
Moment bezwładności na skręcanie prostokąta i-tego (p. 13) wynosi:
Iv,i=βi⋅t3i⋅li
gdzie współczynnik βi zależy od proporcji boków prostokąta poddanego czystemu skręcaniu , dla kwadratu ($t_i =l_i$) wynosi 0,141, a dla długiego prostokąta (li/ti>10) wynosi 1/3 , Wytrzymałość materiałów dla wydziałów budowlanych. PWN)), dla innych wartości można je odczytać z tab.1.
Z założenia #2 mamy
Mv=n∑i=1Mv,i
skąd po podstawieniu wyrażeń na momenty skręcające z założenia Z1, otrzymujemy wyrażenie na moment bezwładności czystego skręcania całego przekroju
Iv=n∑i=1Iv,i=n∑i=1βi⋅t3i⋅li≈k3n∑i=1t3i⋅li
W oszacowaniu Iv wprowadzono współczynnik korekcyjny k, uwzględniający zwiększenie momentu bezwładności przekroju otwartego ze względu na wyokrąglenia naroży oraz pochylenie stopek profili wg. tab.6.
W tab.1 zamieszczono charakterystyki kilku otwartych przekrojów cienkościennych [8]
Naprężenie styczne czystego skręcania w złożonym, otwartym przekroju cienkościennym są maksymalne przy krawędzi ścianki (rys. 6) i wynoszą
maxτi=MvIv⋅ti
Przekrój zamknięty
Przekroje zamknięte charakteryzują się dużą odpornością na skręcanie, wielokrotnie większą od profili otwartych. W przypadku obciążenia pręta skręcaniem należy stosować profile o przekroju zamkniętym, a nie otwartym, czyli stosować rury a nie dwuteowniki.
W ściankach przekroju zamkniętego (rys.7) poddanego swobodnemu skręcaniu powstają tylko naprężenia styczne, które są równomiernie rozłożone po grubości t ścianki, czyli inaczej niż w przekroju otwartym (rys. 6), w którym zmieniają się liniowo i są zerowe w środku ścianki. Strumień naprężeń stycznych w ściance profilu zamkniętego wynosi
qv=τv⋅t
i ma stałą wartość, niezależnie od współrzędnej bieżącej przekroju s, określającej położenie punktów na konturze. Strumień naprężeń stycznych (25) oblicza się ze wzoru Bredta
qv=Mv2Ω⋅t
gdzie Ω jest podwojonym polem powierzchni ograniczonej konturem (linią środkową ścianek) – rys.4 w oogólności wyznaczanym z formuły:
Ω=12∮hds

Rys.7 Przekrój zamknięty: S- środek ścinania, s -współrzędna bieżąca obwodu przekroju pręta, A- pole wewnątrz linii środkowej, Ω – podwojone pole środkowe, [17]
Moment bezwładności swobodnego skręcania zamkniętego przekroju cienkościennego Iv, występujący w (23) wyznaczamy z zależności
Iv=4Ω2(∮dst(s))−1=(2Ω)2⋅t0¯s0
gdzie sprowadzona długość obwodu przekroju ¯s0 do grubość ścianki t0 w dowolnie wybranym miejscu przekroju poprzecznego wynosi:
¯s0=∮t0t(s)
Dla t(s)=t=const s0=s (s- rzeczywista długość obwodu przekroju ), moement bezwładności na skręcanie wynosi Iv=ts⋅(2Ω)2,.
W tab.2 zamieszczono charakterystyki kilku zamkniętych przekrojów cienkościennych za pracą [8].
Skręcanie nieswobodne
Skręcanie nieswobodne (skrępowane) przedstawimy według technicznej teorii prętów cienkościennych Własowa (1959) [3].
Istota nieswobodnego skręcania
W większości praktycznych przypadków deplanacja przekrojów pręta skręcanego (rys. 1c i rys.6) nie może rozwijać się swobodnie, czyli skręcanie jest nieswobodne. Przykładem skręcanie skrępowanego jest zwykłe utwierdzenie końca pręta lub specyficzny, symetryczny sposób przyłożenia momentów skręcających, wymuszający płaskość przekroju leżącego w płaszczyźnie symetrii.
W technicznej teorii nieswobodnego skręcania przyjmuje się, że słuszne są zależności wyprowadzone dla czystego skręcania. Jednostkowy kąt skręcenia Θ definiuje się podług zależności (7) w funkcji położenia przekroju na osi pręta (x), proporcjonalnie do zmieniającego się po osi pręta momentu. czystego skręcania Mv(x), który jest częścią cąłkowitego momentu skręcającego przekrój Mx (1). Przemieszczenie podłużne punków przekroju pręta są różne zależnie od położenia punktu na konturze przekroju ( rys.1c); są nazywane przemieszczeniami paczenia. i są proporcjonalne do współrzędnej wycinkowej przekroju ze współczynnikiem proporcjonalności Θ, określonego w zdaniu poprzednim , wg zależności:
u(x.y.x)=−Θ(x)⋅ω
Współrzędna wycinkowa jest specyficzną współrzędną, określającą położenie punktu P (p. rys. 11 dla przekoroju otwartego rys. 13 dla przekroju zamkniętego).
Rozkład przemieszczeń (30) jest odejściem od klasycznego założenia Bernoulliego o płaskim przekroju na rzecz sztywnego konturu, który paczy się (deplanuje). Paczenie przekroju jest stowarzyszone z nową siła przekrojową – bimomentem.
Bimoment – nowa siła przekrojowa o dyskusyjnej naturze matematycznej
Bimoment jest nową silą przekrojową specyficzną da prętów cienkościennych proporcjonalnyądo drugiej pochodnej kąta skręcenia przekroju podług zależności. (???). Bimoment jest analogiem momentu zginajacego stowarzyszoego z krzywizną pręta w zagadneiniu zginania
W tradycyjnej mechanice przyjmuje się, że bimoment jest liczbą (np. [13], [14]), której można przyporządkować nieskończenie wiele bipar (par momentów o przeciwnym zwrocie), a tymi można obciążyć pręt. W pracy [18] przedstawiono twierdzenia dotyczące algebry układów wektorów sił przyłożonych do płaskiej, sztywnej jedynie w swej płaszczyźnie, linii materialnej, to jest przy założeniu sztywnego konturu przekroju pręta cienkościennego. Rozważania prowadzono w ramach klasycznej algebry wektorów, z której można wyciągnąć wniosek , że bimoment jest liczbą (skalarem).
Autor artykułu stoi na stanowisku, że zagadnienie mechaniki prętów cienkościennych oraz algebry sił w tym bimomentu należałoby rozpatrywać w ramach ogólniejszej algebry Clifforda [19]. Zgodnie z tą teorią siła fizyczna jest trójką (skalar, wektor, biwektor), która jest nazywana pseudowektorem (qwaterionem lub auaterionem). Ta trójka obiektów posiada wszystkie cechy tradycyjnego wektora (wartość, kierunek i zwrot), skalara oraz biwekotra. Obiekty takie są analizowane w algebrze geometrycznej Clifforda [20], [21], [22], [23]. Algebra Clifforda jest szczególnie użyteczna w analizie dużych obrotów, a obecnie obserwuje się renesans jej zastosowań.
Bimoment nie jest ani skalarem , ani wektorem – jest obiektem matematycznym o nazwie quaterion (lub pseudowektor).
W podejściu technicznym bimoment często przedstawia ię jako „parą momentów” lub „momentem momentów”., a nigdy jako skalar, wektor lub para wektorów (moment).
W niniejszym artykule nie zajmujemy się sformalizowana, matematyczną teorią bimomentów oraz z ich transformacją, a dla dalszych rozważań wystarczająca jest techniczna definicja bimomentu jako pary momentów.
Przy zginaniu nieswobodnym powstaje jeszcze jedna nowa siła przekrojowa – moment giętno-skrętny Mω (???).
Modelowanie pręta cienkościennego
Kontur przekroju cienkościennego jest sprowadzony do linii środkowych ścianek pręta, które są umieszczone w osiach ścianek a cała grubość odcinka skupia się na linii środkowej. To znaczy przekrój składa się z odcinków, o długości li którym przypisano grubości ti Linię środkowąna rys. 6 oznaczono linią przerywaną. W takim modelu przekroju cienkościennego charakterystyki geometryczne przekroju (A, Iy, Iz, itd) różnią się od wartości wyliczonych z klasycznych definicji (p. przykład 2).
Specyficznym typem podpory w analizie prętów cienkościennych jest podpora widełkowa, czyli taka w której nie jest możliwy kąt obrotu. Na rys. 8 pokazano jedną z możliwych realizacji podpory widełkowej, a obok pokazano stosowany symbol podpory i warunki brzegowe generowane przez podporę widełkową.

Rys. 8 Podpora widełkowa pręta cienkościennego: a ) przykład realizacji podpory, b) symbol i warunki brzegowe ( uzupełniony rys. 4.13 (Bródka J., Broniewicz N., Giżejowski M. , Kształtowniki gięte. Poradnik projektanta, PWT, Rzeszów , 2006))),
Pole przemieszczeń
Pole przemieszczeń liniowych pręta cienkościennego,opisanych formułą (30), zaprezentowano na rys.9. Znamienna dla prętów cienkościennych deplanacja przekroju objawia się różnicami przemieszczeń podłużnych punktów przekroju u1 do u5z konfiguracji początkowej (przed odkształceniem) Ω(0) do konfiguracji po odkształceniu Ω(N). W celu zwiększenia czytelności na rys.9. nie wniesiono oznaczeń kątów obrotu oraz paczenia.
Przemieszczenia poprzeczne do osi pręta: pionowe w oraz poziome v rozpatruje się zgodnie z klasyczną teorią belki Bernoulliego.
Podstawowe zależności różniczkowe
Równanie skręcania nieswobodnego
Równanie skręcania pręta o przekroju cienkościennym (otwartym , zamkniętymi mieszanym), obciążonego rozłożonym wzdłuż jego osi x momentem skręcającym o intensywności ms(x) ma następującą postać
EIω⋅d4φdx4–GIv⋅d2φdx2=ms(x)
gdzie:
φ(x) – kąt skręcenia przekroju pręta wokół osi x,
EIω – sztywność giętno-skrętna przekroju pręta: E- moduł Youmga, Iω – wycinkowy moment bezwładności,
GIv – sztywność skrętna przekroju pręta: G- moduł Kirchoffa, Iv – moment bezwładności czystego skręcania (Saint Venanta)
czyli ma postać różniczkową podobną do klasycznego równania zginania względem osi y belki Timoshenko o zwartym przekroju, sztywności giętnej EIy i sztywności na ścinanie kGA obciążonej rozłożonym, i pionowymi siłami qz
EIy⋅d4w/dx4+kGA⋅(d2wdx2+d2ψdx2)=qz(x) , gdzie w(x) jest ugięciem pręta, wywołanym obciążeniem pręta qz(x).
Ugięcie w(x) belki i kąt obrotu zginania ψ(x)w punkcie belki x zależą od sztywność zginania EIy i, sztywności na ścinanie kGA belki ( k to współczynnik uwzględniający kształt przekroju, G – moduł Kirchhoffa, A – pole przekroju poprzecznego).
Rozwiązanie równania skręcania nieswobodnego
Równaniem stanu pręta cienkościennego jest równanie kąta skręcenia φ (31) . Inne zmienne stanu. w tym siły przekrojowe można wyznaczyć z formuł zależnych tylko od kąta skręcenia i jego pochodnych. Znajdziemy rozwiązanie równania (31) zapisanego w postaci
φ(iv)(x)–k2φ(ii)(x)=ms(x)EIω
gdzie
k2=GIvEIω
Piechnik (2000) [13] pokazał że równanie różniczkowe czwartego rzędu (32) można sprowadzić do równania trzeciego rzędu w postaci
φ(iii)(x)–k2φ(i)(x)=f(x)
gdzie wyraz wolny – funkcja f(x) jest zdefiniowana jako
f(x)=–1¯EIωMx,S
gdzie Mx,S (1) jest momentem sił wokół osi x przekroju o rzędnej x, liczonym względem środka zginania S.
Sprowadzony moduł Younga ¯E wynosi
¯Edef=E1–ν2
i pozwala stosować proste prawe Hooka wiążące odkształcenia i naprężenia wzdłuż osi pręta cienkościennego
σx=¯E⋅εx
Rozwiązanie zwyczajnego równania różniczkowego (34) jest sumą całek: ogólnej równania jednorodnego φo(x) i szczególnej równania niejednorodnego.φs(x)
φ(x)=φo(x)+φs(x)
Zgodnie z metodą podaną przez Eulera poszukujemy rozwiązania równania jednorodnego
φ(iii)(x)–k2φ(i)(x)=0
w postaci erx. Po podstawieniu zmiennych uzyskujemy następujące równanie charakterystyczne
r3−k2r=0
które ma trzy pierwiastki
r1=0,r2=k,r3=−k
Całkę ogólną równania różniczkowego liniowego (39) można zapisać więc jako kombinację liniową funkcji wykładniczych:
Ae0+Bekx+Ce−kx
Po zastąpieniu funkcji wykładniczych ich liniową kombinacją funkcji wykładniczych, zdefiniowanych jako funkcje hiperboliczne
sinh(kx)=ekx−e−kx2 oraz
cosh(kx)=ekx+e−kx2
φo(x)=C1+C2sinh(kx)+C3cosh(kx)
Całka ogólna metody parametrów początkowych
Całkę ogólną (43) można przekształcić do postaci dogodnej do obliczeń metodą parametrów początkowych i korespondującą z równaniem 4 rzędu (32).
φo(x)=C1+C21ksinh(kx)+C31k2¯EIω[1−cosh(kx)]+C41k3¯EIω[kx–sinh(kx)]
Całka szczególna φs(x) zależy od rozkładu momentów skręcających Mx,S(x) po długości pręta i zwykle przyjmuje się ją metodą przewidywania.
Na przykład w przypadku pręta obciążonego momentem skręcającym ms stałym na całej długości pręta (schemat 1.3.5.7 w tab. 4) całkę szczególną można zapisać w postaci
φs(x)=−msk4¯EIω[12(kx)2+1–cosh(kx)]
Rozwiąniem zagadnienia jest suma całki ogólnej i szczególnej (38).
Warunki brzegowe a stałe całkowania
Stałe całkowania C1,C2,C3,C4 wyznacza się ze statycznych i kinematycznych warunków brzegowych.
W tab. 4 zapisano warunki dla kliku często spotykanych schematów belek cienkościennych.
W innych przypadkach dobór warunków brzegowych przebiega z uwzględnieniem zasad:
- C1 jest kątem obrotu φ(x=0) na początku belki (x=0),
- C2 jest pochodną kąta skręcenia φ(x=0),
- C3 jest bimomentem Bφ(x=0) ,
- C4 jest momentem skręcającym Mx(x=0)
przy czym:
- przy braku obrotu φ=0 i swobodzie skręcania Mx=0, warunek brzegowy; φ(iii)–k2⋅φ(i)=0
- dla swobodnego końca obciążonego momentem Ms warunek; ¯EIω∼φ(iii)–k2φ(i)−Mx
- w przypadku swobodnego paczenia Bω=0, czyli φ(ii)=0
- przy braku spaczenia końcowego przekroju u=0, czyli φ(i)=0
Analizy prowadzone w pracy [13] pozwalają postawić warunki brzegowe dla bardziej złożonych przypadków z dokładnością pozwalającą uwzględniać pojedyncze więzi składających się na deplanację przekroju.
Siły przekrojowe w skręcanym pręcie cienkościennym
Skręcanie pręta cienkościennego wywołuje siły przekrojowe, które można wyznaczyć po wyznaczeniu funkcji kąta skrecenia pręta φ(x) ($\ref{38}$) z zależności różniczkowych:
moment czystego skręcania (Saint Venanta)
Mv(x)=GIvφi(x)
bimoment
B(x)=–EIωφii(x)
moment giętno- skrętny
Mω(x)=–EIωφiii(x)
Praktycznie w każdej sytuacji pręty skręcane są jednocześnie zginane i ścinane, Siły przekrojowe od obciążenia zginającego (poprzecznego do pręta w płaszczyżnie i z płaszczyzny zginania określa się w sposób standardowy, znany z podstawowego kursu mechaniki prętów,
Charakterystyki geometryczne przekroju cienkościennego
Przekrój cienkościenny otwarty
Współrzędna wycinkowa
Lokalny krzywoliniowy układ współrzędnych (x,s,n) w dowolnym punkcie powierzchni środkowej pręta (rys. 5) definiujemy w ten sposób, że pierwszą oś układu określamy jako równoległą do osi globalnej x, przyjmując dla niej to samo oznaczenie (rys 5a) o wersorze →ex. Druga oś s jest styczna do linii środkowej, definiuje wersor →es. Oś trzecia n jest wyznaczona wersorem →en=→ex×→es Początek osi spółrzędnej łukowej s przyjmujemy w ustalonym punkcie O na linii środkowej, zwanym początkiem przekroju. Punkt B jest dowolnie przyjętym biegunem z którego prowadzi się promienie ρ(s) do punktu bieżącego M(s) i promień ρ(s+ds) do punktu bieżącego po przyroście współrzędnej łukowej s.
Obok zdefiniowanych wyżej współrzędnych w analizie prętów cienkościennych posługujemy się jeszcze współrzędną wycinkową ω. Przy znajomości linii środkowej, położenia punktu początkowego (O) i bieguna (B) – współrzędna wycinkowa w sposób jednoznaczny wyznacza położenie dowolnego punktu (P). Wartość współrzędnej wycinkowej punktu (P) zależy zarówno od położenia bieguna (B) jak i punktu początkowego (O) – najczęściej punktu zerowej współrzędnej wycinkowej.

Rys. 11 Współrzędne przekroju cienkościennego otwartego (s, n, ω), gdzie ω – współrzędna wycinkowa: a) układ globalny (x,y,z) i lokalny (x,n,s) współrzędnych punktu, b) współrzędna wycinkowa ω zmodyfikowane [24]
Współrzędne wersorów lokalnego układu, w dowolnym punkcie P(s) odniesione do układu globalnego zapiszemy w postaci macierzy [C]:
[C]=[y′(s)→exz′(s)→ex–z′(s)→eny′(s)→en]
gdzie: y′(s)=dy(s)ds, z′(s)=dz(s)ds
Na linii środkowej przekroju poprzecznego pręta przyjmiemy dowolny punkt (P) o wektorze wodzącym →ρ(s) zaczepionym biegunie B o współrzędnych B(0 ,yb ,zb) oraz utworzymy iloczyn skalarny, który oznaczymy przez [dωB](s):
[dωB]=→ρ(s)·→en(s)ds
Geometryczną interpretacją bezwzględnej wartości ||[dω](s)|| jest podwojone pole powierzchni zawartej pomiędzy promieniami wodzącymi →ρ(s) i →ρ(s+ds), bowiem bezwzględna wartość iloczynu skalarnego przedstawia długość rzutu wektora ρ(s) na normalną do krzywej, będącego wysokością wspomnianego trójkąta.
Współrzędną wycinkową względem biegunu (B) definiujemy jako :
ωB(s)=s∫0[dωB](s)=s∫0→ρ(s)⋅→ends
Jeśli współrzędna łukowa przebiega po kilku ściankach(bokach linii środkowej), to wówczas powyższe możemy zapisać jako sumę całek
ωB(s)=∫s0→ρ(s)⋅→ends=j−1∑i=1s′M∫s′O→ρ(s)⋅→ends+s∫sj0→ρ(s)⋅→ends
gdzie: s′O i s′M jest wartością współrzędnej łukowej punktu i-tej krawędzi początkowego (O) i końcowego (P) na odcinku linii środkowej odpowiednio,
Bezwzględna wartość każdego i-tego składnika sumy jest równa podwojonemu polu zawartemu pomiędzy wektorem wodzącym punktu s′O, i-tym bokiem linii środkowej oraz wektorem wodzącym punktu s′M.
Zależność (51) można zapisać we współrzędnych globalnych (y,z) w postaci;
ωB=∫s0(z–zB)dy–(y–yB)dz=∫so(zdy–ydz)–zB∫s0dy+yB∫s0dz=ωO(s)–zB[y(s)–yO]+yB[z(s)−zO]
Ponieważ ωB jest współrzędna wycinkową dla bieguna w punkcie (B), a ωO dla bieguna we wcześniejszym punkcie (O), to wzór (53) określa zależność pomiędzy współrzędnymi wycinkowymi przy zmianie bieguna z dowolnego punktu (O) na dowolny punkt(B). Formuła (53) jest w praktyce stosowana do obliczania współrzędnych wycinkowych przekroju dla znanego puntu początkowego (O) (w którym współrzędna jest znana i zerowa ωO=0).
Środek zginania (ścinania) i główne charakterystyki przekroju
Biegun (B) (rys.5b) bęzie środkiem zginania (S) (w literaturze nazywany również środkiem ścinania lub skręcania), jeśli zerują się wycinkowe momenty odśrodkowe wyznaczone względem tego bieguna:
Iω,ydf=∫cω(s)y(s)t(s)ds=∬Aω(s)y(s)dA=0Iω,zdf=∫cω(s)z(s)t(s)ds=∬Aω(s)z(s)dA=0
Po przekształceniu (54) z zastosowaniem z formuły transformacji współrzędnej wycinkowej- zamiany bieguna (53) można uzyskać formuły do wyznaczania współrzędnych środka zginania (s). Najprostsze formuły na współrzędne środka zginania (S) S(yS,zS) uzyskuje się w układzie głównych centralnych osi bezwładności czyli w układzie osi (y,z) przechodzących przez środek ciężkości (C) S(yC,zC) i w którym moment bezwładności Iy jest największy a moment bezwładności Iz najmniejszy.
yS=yB+1Iy∫Aω′BzdA;zS=zB+1Iz∫Aω′BydA
gdzie ω′B współrzędna wycinkowa punktu (P) wyznaczona z dowolnego (pomocniczego) bieguua (B) zliczana od dowolnego (pomocniczego) punktu początkowego (O’).
Do wyznaczenia głównej współrzędnej wycinkowej oprócz ustalenia głównego bieguna, którym jest środek zginania (S) z zależności (55) należy ustalić jeszcze główny punkt początkowy (O) z zależności
ωS=ω′S–1A⋅∫Aω′SdA
gdzie ω′S jest polem wycinkowym względem środka zginania (S) od pomocniczego punktu początkowego (O’)
Technika obliczania charakterystyk geometrycznych przekroju cienkościennego otwartego
Obliczanie charakterystyk geometrycznych przekroju cienkościennego jest w praktyce prowadzone graficzną c metodą ałkowania iloczynu funkcji (zwaną również metodaąWereszczagina). Jest to technika obliczania całęk iloczynu f⋅g dwóch funkcji f(s) i g(s) poprzez „mnożenie” pola Af pod wykresem funkcji f przez rzędną sg funkcji g pod środkiem pola pola Af
Wyznaczanie charakterystyk dotyczy przekroju cienkościennego sprowadzonego do konturu złożonego z linii środkowych ścianek ( p. modelowanie pręta cienkościennego).
Pole powierzchni
Pole przekroju pręta cienkościennego wyznaczane jest w sposób odmienny od przekroju zwartego. W konsekwencji modelem przekroju #1 pole powierzchni jest sumą iloczynów długości i grubości każdego odcinka
A=∑ili⋅ti⋅1
Wyrażenie (57) ma następującą interpretację całkowania geometrycznego: na odcinku (ściance) „i” funkcja f(s)=li⋅ti, a funkcja g(s)≡1 jest stała na całym odcinku, więc rzędna w środku prostokąta lixti wynosi sg=1.
Momenty statyczne w układzie globalnym
SY=∑ili⋅ti⋅Zi
SZ=∑ili⋅ti⋅Yi
gdzie: Yi=(Yi,1+Yi,2)/2/,;/, Z_i= (Z_{i,1}+ Z_{i,2)/2 $ są współrzędnymi środka ciężkości ścianki „i” w układzie globalnym (X,Y)
Środek ciężkości przekroju (C)
Yc=SZA
Zc=SYA
Centralne momenty bezwładności w układzie globalnym
Pokażemy różnicę w wyznaczaniu momentu bezwładności przekroju zwartego i cienkościennego na przykładzie przekroju dwuteowego z rys 5 o całkowitej wysokości h oraz szerokości półki b. Wysokość środnika wynosi hw=h−2⋅tf
Dla przekroju dwuteownika (hw×tw)/(b×tf) traktowanego jako zwarty moment bezwładności (z pominięciem wyokrąglę naroży) względem osi poziomej Iy oraz pionowej Iz, a także moment dewiacji Iyz wynoszą:
IY,zwarty=2⋅[b⋅t3f12+(b⋅tf)⋅((hw–tf)2)2]+[tw⋅h3w12+(hw⋅tw)⋅02]
IZ,zwarty=2⋅[tf⋅b312+(b⋅tf)⋅02]+[hw⋅t3w12+(hw⋅tw)⋅02]
IYZ,zwarty=2⋅[0+(b⋅tf)⋅hw–tf2⋅0]+[0+(hw⋅tw)⋅0⋅0]=0
Dla przekroju cienkościennego w wyrażeniach na momenty bezwładności składniki, w których występuje grubość ścianki podniesiona do potęgi, o to dlatego, że cała grubość odcinka skupiona jest na jego linii środkowej. W rezultacie dla dwuteownika otrzymujemy
IY,cienko=2⋅[(b⋅tf)⋅((hw−tf)2)2]+[tw⋅h3w12+(hw⋅tw)⋅02]
IZ,cienko=2⋅[tf⋅b312+(b⋅tf)⋅02]+[(hw⋅tw)⋅02]
IYZ,cienko=2⋅[0+(b⋅tf)⋅hw–tf2⋅0]+[0+(hw⋅tw)⋅0⋅0]=0
Kąt nachylenia osi głównych α
tg(2α)(=T)=–IYZIY–IZ→α=arctg(T)2
Główne centralne momenty bezwładności
Podstawowe zależności z kursu wytrzymałości materiałów dla przekrojów zwartych są ważne również dla przekrojów cienkościennych
Im=Iyc+Izc2
IΔ=√(Iyc−Izc2)2+I2yc,zc
Iy=Imax=Im+IΔ=
Iz=Imin=Im–IΔ
Dla przekrojów cienkościennych można też przeprowadzić całkowanie graficzne poprzez „mnożenie” wykresów (p. przykład 2),
Przekrój cienkościenny zamknięty
Współrzędna wycinkowa
Zamknięty przekrój cienkościenny jest niewyznaczalny w schemacie rzeczywistym i w celu wyznaczenia głównych współrzędnych wycinkowych i związanych z nimi charakterystyk należy utworzyć schemat podstawowy – przekrój otwarty) przez przecięcie przekroju w dowolnym punkcie (k) – rys. 13.

Rys. 13. Wyznaczanie głównej współrzędnej wycinkowej przekroju zamkniętego [8]
Uogólniony wycinkowy moment bezwładności dla dowolnego przekroju zamkniętego oblicza się ze wzoru
I¯ω=∮¯ω2dA
w którym ¯ω oznacza uogólnione pole wycinkowe względem środka ścinania S przekroju od głównego punktu początkowego O.
W celu wyznaczenia (73) najpierw oblicza się pomocnicze współrzędne wycinkowe ¯ω′B w podstawowym przekroju otwartym (dla obszaru przeciętego w dowolnym punkcie (K) – rys. 13 ) względem przyjętego (tego samego) bieguna pomocniczego (B) od doolnego pomocniczego pubktu począgtkowego (O.)
Pole wycinkowe względem bieguna B od dowolnego punktu początkowego O’ na konturze, oblicza się ze wzoru
¯ω′B=ω′B−¯s¯s0⋅2Ω
gdzie:
¯s=s∫0t0t(s)ds – sprowadzona współrzędna s (długość konturu) liczona od punktu O’.
Środek zginania (ścinania)
Współrzędne środka ścinania (skręcania) (S) przekroju zamkniętego wyznacza się w układzie centralnych , głównych osi bezwładności, (y, z) z następujących wzorów:
yS=yB+I¯ω,yIy;zS=zB+I¯ω,zIz
gdzie:
yB, zB – współrzędne dowolnie przyjętego bieguna B, który na ogół przyjmuje się na konturze,
$I_{\overline \omega, y}= \oint \ovyerline \omega_B \cdot z \, dA ,I_{\overline \omega, z}= \oint \ovyerline \omega_B \cdot y \, dA ,I_y,I_z$ – główne momenty bezwładności przekroju poprzecznego względem głównej osi poziomej y i pionowej z odpowiednio.
Pręty o przekroju trójkątnym i inne w kształcie wieloboku foremnego, ale o stałej grubości ścianki (t=const) nie ulegają deplanacji i ich charakterystyki wycinkowe są równe zeru.
Stan naprężeń
Naprężenia
Brak swobody deplanacji wywołuje naprężenie normalne do przekroju równe
σω=Edudx=−Ed2φdx2
gdzie E- moduł odkształcalności podłużnej (Younga), a φ jest kątem skręcenia przekroju. Zachodzi
Można pokazać, że naprężenia normalne σω można wyrazić za pomocą nowej siły przekrojowej bimomentu Bω w postaci [16]:
σω=BωIω⋅ω
Moment giętno-skrętny wywołuje naprężenia styczne τω o stałej wartości po grubości ścianki, które można wyznaczyć z zależności:
τω=Mω⋅SωIω⋅t
gdzie: Sω – wycinkowy moment statyczny przekroju, t – grubość ścianki przekroju., przy czym
Mx=Mυ+Mω
Naprężenie styczne od momentu czystego skręcania wyznaczamy ze wzoru[25]:
τυ=MυWυ
gdzie:
Wυ=Iυmaxti
jest wskaźnikiem wytrzymałości przekroju na skręcanie, równym stosunkowi momentu bezwładności na skręcanie i maksymalnej grubości ścianki.
Dla przekroju cienkościennego przekroju cienkościennego otwartego, złożonego moment bezwładności na skręcanie można oszacować z zależności (29)
W przypadku przekroju zamkniętego o powierzchni Ω zamkniętej wewnątrz linii środkowych przekroju, wskaźnik wytrzymałości przekroju na skręcanie, wynosi
Wυ=2⋅Ω⋅min{ti}
gdzie występuje minimalna grubość t ścianki przekroju.
Sztywność na skręcanie przekroju zamkniętego jest wielokrotnie większa od przekroju otwartego i jest proporcjonalna do podwojonego pola zawartego wewnątrz przekroju.
Ze skręcaniem skrępowanym mamy do czynienia w przypadku skrępowania deplanacji przekroju poprzecznego pręta. W wyniku tego w przekroju powstają naprężenia normalne σ¯ω oraz dodatkowe naprężenia styczne τ¯ω . Naprężenia te wyznacza się ze wzorów:
σ¯ω=B¯ω⋅¯ωI¯ω
τ¯ω=–M¯ω⋅¯S¯ωI¯ω⋅t(s)
gdzie:
¯ω – główne, uogólnione pole wycinkowe (współrzędna wycinkowa przekroju);
I¯ω – uogólniony wycinkowy moment bezwładności;
¯S¯ω – uogólniony wycinkowy moment statyczny przekroju.
Dodatnim naprężeniom τ¯ω odpowiada zwrot zgodny z kierunkiem dodatniego przyrostu współrzędnej s.c
W odróżnieniu od prętów o przekrojów otwartym , giętno skrętne naprężenia styczne τ¯ω stanowią samozrównoważone układy sił i tym samym cały moment skręcający Mv przenoszą naprężenia styczne swobodnego skręcania τv.
Całkowite naprężenie styczne w dowolnym punkcie przekroju znajduje się jako sumę τv+τ¯ω.
Charakterystycznym punktem przekroju cienkościennego jest środek zginania (zwany czasami środkiem skręcania lub środkiem ścinania). Jest to punkt położony w przekroju lub poza nim, taki, że obciążenie przechodzące przez ten punkt nie powoduje skręcania pręta, a tylko jego zginanie poprzeczne. W takim przypadku pozostaje słuszna zasada płaskich przekrojów Bernulliego i pręt jest poddany zginaniu bez udziału skręcania.
Pręty wielogałęziowe oraz kratownice wielopasowe
Pręty (kratownice) wielopasowe można przystosować do przenoszenia momentu skręcającego poprzez stosowny dobór zamykającego skratowania lub przewiązek. Przykładem skręcanych kratownic są galerie przenośnikowe, bramki nad autostradami lub słupy energetyczne. W ścianach takich kratownic powstają siły typowe dla swobodnego skręcania prętów o przekroju zamkniętym [8].
Pręty wielogałęziowe
Na rys.14 pokazano cienkościenny przekrój quasi-zamknięty oraz przykłady prętów o dwu gałęziach, które mogą być łączone skratowaniem lub przewiązkami w sposób często stosowany przy konstruowaniu dwugałęziowych slupów. Przekrój pozornie zamknięty skratowaniem lub przewiązkami można zamodelować zastępczym przekrojem cienkościennym i zastosować analogię prętową.

Rys. r4 Pręty wielogałęziowe: a) przekrój b), c) skratowany, d) z przewiązkami [8]
Sztywność przekroju na czyste skręcanie bez uwzględnienia pozornego zamknięcia wyznacza się standardowo z zależności (\ref23}). Oznaczymy go jako moment bezwładności czystego skręcania przekroju otwartego
Iv,open=k/3∑ili\cdott3i
Zwiększenie sztywności przekroju na skutek zamknięcia go skratowaniem lub przewiązkami można oszacować z zależności [8]:
¯Iv,close={EG(Ωb)2⋅ Ak⋅ sin2αcosα, dla skratowania typu N (bezsłupkowe) – rys. 14c EG(Ωb)2⋅As⋅Ak⋅ sin2α⋅cosαAs+Ak sin3α,dla skratowania typu V – rys. 14b 24EG⋅(Ωb)2⋅Ifa2(1+2b1⋅Ifa⋅Ip+28,8⋅EGIfab1Ap)−1, dla przewiązek – rys. 14d
gdzie:
E,G – moduł Younga i moduł Kirchoffa materiału;
Ω=2⋅b⋅h
b,b1- prześwit i odległość osiowa między gałęziami pręta;
a- odległość osiowa między przewiązkami lub węzłami skratowania.
Ak, As , Af – pole przekroju krzyżulca, słupka, gałęzi;
If- moment bezwładności giętnej pojedynczej gałęzi
Całkowity moment bezwładności czystego skręcania przekroju quasi-zamkniętego jest szacunkowo sumą powyższych części
Iv=Iv,open+IV,close
Równanie różniczkowe pręta cienkościennego rozwiązuje się dla sztywności pręta (87). W wyniku uzyskujemy siły przekrojowe B, M_\omegam_v(47)do(49),przyczymnależyuwzględnić,żeuzyskanymomentM_vodpowiadaswobodnemuskręcaniuprętaqusi−zamkniętegoozwieszonejsztywnościwstosunkudoprętaniezamkniętegoskratowaniemlubprzewiązkami.WceluobliczenianaprężeństycznychswobodnegoskręcanianależymomentM_vpodzielićnaczęśćodpowiadającąprzekrojowiotartemu M_{v,open} = Mv \cdot \cfrac{{I_{v,open}}{I_v} oraz zamkniętemu Mv,close=Mv−Mv,open. Dla tak obliczonych momentów można stosować wyrażenia odpowiednie dla danego typu przekroju ( otwarty lub zamknięty).
Siły w prętach skartowań zamykających przekrój pręta wyznacza się uwzględniając fakt, że
Kratownice wielopasowe
Metoda analogii pręta cienkościennego może by stosowana do wstępnego zaprojektowania skręcanej kratownicy wielopasowej w sposób pokazany na rys. 2.
Do wstępnego wymiarowania prętów kratownic przestrzennych można pomijać zmniejszający wpływ quasi-zamknięcia prea skratowaniem i przyjmować kv (???0 równw jedności. Przy takim założeniu dla kartownicy czteropasowej lub trójpasowej pokazanych na rys. 15, mamy:
a wówczas w przypadku przekrojów pokazanych na rys. 3, mamy:
\begin{equation} \overline Q =\begin {cases}
\cfrac{M_v\cdot a }{2ab}=\cfrac{M_v }{2b} , & \text { dla kraty pionowej kratownicy czteropasowej }\\
\cfrac{M_v\cdot b }{2ab}=\cfrac{M_v }{2a} , & \text { dla kraty poziomej kratownicy czteropasowej }\\
\cfrac{2M_v}{a \cdot \sqrt{3) , & \text { dla karty bocznej kratownicy przekroju trójkąta równobocznego o boku a}\\
\end {cases} \label{89} \end{equation}
Obciążenia poprzeczne Q są przyłożone do wydzielonej kratownicy płaskiej zgodnie z przebiegiem momentu skręcającego Mv(x) po długości pręta x. Siły w prętach kratownicy wyznacza się poprzez rozwiązanie kratownic płaskich pod takim obciążeniem. Należy pamiętać o sumowaniu sił w pasach od poszczególnych kratownic składowych z uwzględnieniem kierunku obciążenia Q. Na przykład dla kratownicy, K1 (rys. 15) :przy działaniu momentu Mv prawoskrętnie:
- obciążenie zewnętrzne: Q pionowo do góry,
- sił w pasach kratownic pas dolny – ściskany, pas górny PG – rozciągany,
- sumaryczne siły osiowe : pas górny lewy – K3(PG)+K4(
Wynika stąd, że w przypadku przekroju kwadratowego a=b siły osiowe w pasach będą zerowe, a całkowity moment skręcający przenosi skratowanie ścian.
Kratownica przestrzenna, przenosząca skręcanie powinna być skratowana na wszystkich płaszczyznach bocznych oraz powinna być wyposażona w poprzeczne przepony w miejscach przyłożenia momentów (w tym w zamocowaniu) oraz nie rzadziej niż co 8-me pole kratownicy. Do wyznaczania przemieszczeń kratownicy od wpływu skręcania należy uwzględniać współczynnik kv<1. Analizując kratownice w modelu zastępczym pręta w złożonym stanie obciążenia (zginanie, ścinanie, skręcanie) należy stosować element Timoshenko.
Przykłady rachunkowe
Przykłady w przygotowaniu
Przykład 1 [Współrzędne wycinkowe dla typowych przekrojów otwartych prętów stalowych ; I, C, Z]
Na podstawie pracy Rykaluk (20126) [15].
Wyznaczyć główne współrzędne wycinkowe \omega, moment bezwładności na skreęanie I_v oraz główny wycinkowy moment bezwładności i_\omega
dla dwuteownika, ceownika i zetownika.
Wykresy głównych współrzędnych wycinkowym podsumowano na rys
Dwuteownik I (H)
Analizowany jest model cienkościennego przekroju ( złożony z linii środkowych ścianek). Wysokość przekroju wynosi h, szerokość stopki b, grubość stopek t_f, a grubość środnika t_w.
Przekrój jest bisymetryczny , więc środek zginania (ścinania) S pokrywa się ze środkiem ciężkości C
Główne współrzędne wycinkowe wyznacza się z bieguna, którym jest środek zginania jako podwojone pole zakreślone przez promień wodzący, jak następuje:
punkt P_0 \to \omega_{P0}= 2\cdot h/2 \cdot 0 = 0
punkt 2 \to : \omega_2 = + 2 \cdot 1/2 \cdot h/2 \cdot b/2 = + \cfrac{b h }{4} (znak + , bo promień wodzący zakreślił w prawo),
punkt 1 \to : \omega_1 = – 2 \cdot 1/2 \cdot h/2 \cdot b/2 = – \cfrac{b h }{4} (znak – , bo promień wodzący zakreślił w lewo),
analogiczne rzędne na stopce dolnej 2ynoszą \omega= +/- bh/4
Moment bezwładności na skręcanie
(\ref{23}) \to I_v = \cfrac{1}{3} \left (2 b t_f^3 + h t_w^3 \right)
Główny wycinkowy moment bezwładności:
I_\omega = \cfrac{b}{6} \cdot \left [ \right ] = \cfrac{h^2}{4} \cdot \cfrac{2 t_f b^3}{12}

Rys P1-1. Współrzędne wycinkowe dla dwuteownika, ceownika i zetownika [15]
Przykład 2 [Charakterystyki geometryczne przekroju otwartego ]
Przykład na podstawie opracowania [26]
Wyznaczyć charakterystyki geometryczne I_y, I_z, I_v, I_\omega , oraz \overline S_y, \overline S_z, \overline S_\omega odciętej części przekroju w punkcie C wg rys P2-1 . Układ globalny osi (Y,Z) jest prostopadły so osi pręta X.

Rys. P2-1 Przekrój cienkościenny, otwarty wg [26]
Dane
Na rys. P2-1 przekrój w liniach środkowych ścianek oznaczono linią grubą
Węzły przekroju:
(1) (y_1 = 20 \, ;\, z_1= 0) cm,
(2) (0\, ;\, 0) cm,
(3) (0\, ; , -30) cm,
Ścianki przekroju:
[1] =(1)-(2) l_1= 20 \, cm, t_1= 2 \, cm
[2]=(2)-(3) l_2= 30 \, cm, t_2= 3 \, cm
[3]=(3)- (4) l_3= 13 \, cm t_3= 1,3 \, cm
Środek ciężkości, osie centralne, osie główne
Na rys. P2-2 osie centralne (przechodzące prze środek ciężkości) naniesiono liniami niebieskimi, a osie główne (osie centrale, ale nachylone pod kątem \alpha, względem których momenty bezładności przyjmują wartości ekstremalne, a moment dewiacji jest zerowy) oznaczono liniami bordowymi.
A= 20 \cdot 2,0 + 30 \cdot 3,0 + 13\cdot 1,3 = 146,9 \, cm^2
Momenty statyczne w układzie głównym (X,Y)
S_Y=20 \cdot 2,0 \cdot 0 + 30 \cdot 3,0 \cdot \cfrac{-30}{2} +13\cdot 1,3\cdot (- 30)= – 1 857 \, cm^3
S_Z = 20 \cdot 2,0 \cdot \cfrac {20}{2} + 30 \cdot 3,0 \cdot 0 +13\cdot 1,3\cdot \cfrac{13}{2}= 509,9 \, cm^3
Współrzędne środka ciężkości
y_c=\cfrac{509,9}{146,9}= 3,47 \, cm
z_c=\cfrac{-1 857}{146,9}= – 12,64 \, cm
Centralne momenty bezwładności
I_{84}c} = 20 \cdot 2,0 \cdot (-12,64)^2 + \cfrac{3,0 \cdot 30^3}{12}+ 30 \cdot 3,0 \cdot (30/2 – 12,64) ^2 +13 \cdot 1,3 \cdot (30 – 12,64)^2 = 18735 cm^4.
I_{zc} =\cfrac{2,0 \cdot 20^3}{12}+ 20\cdot 2.0 \cdot (20/2- 3,47)^2 +30 \cdot 3,0 \cdot (3,47)^2 + \cfrac{1,3 \cdot 13^3}{12}+13 \cdot 1,3 \cdot (12/2 -3,47)^2= 4516 \, cm^4
I_{84}c,zc}= 20\cdot 2,0 \cdot (-12,64)\cdot (20/2-3,47)+ 30\cdot 3,0 \cdot (30/2-12,64)\cdot (-3,47)+13 \cdot 1,3 \cdot (30 – 12,64)\cdot (13/2 -3,47)=-3150 \\, cm^4
Kąt nachylenia osi głównych
tg (2 \alpha) = – \cfrac{-3150}{18 735-4516}= 0,143 \to \alpha = 11,95^0
Główne centralne momenty bezwładności
Stosujemy podstawowe zależności wytrzymałości materiałów dla przekrojów zwartych
I_m= \cfrac{I_{84}c} + I_{zc}}{2} = \cfrac{ 18735+4516 }{2}= 11625,5\,cm^2
I_\Delta = \sqrt{\left( \cfrac{I_{84}c} – I_{zc}}{2}\right)^2 +I_{84}c,zc}^2}= \sqrt{\left( \cfrac{18735 – 4516}{2}\right)^2 +3150^2}= 7776,1 \, cm^4
I_y = I_{max}=I_m+I_\Delta = 11625,5+7776,1 =19402 \, cm^4
I_z = I_{min}=I_m – I_\Delta = 11625,5- 7776,1 =3849 \, cm^4
Współrzędne punktów przekroju
Na rys. P2-3 przedstawiono wykresy współrzędnych głównych y, z oraz współrzędnej wycinkowej względem dowolnego bieguna pomocniczego (B) oraz początkowego punktu zliczania (O). Biegun pomocniczy (B) przyjęto w górnym lewym narożu przekroju, a punku (O) w dolnym lewym narożu przekroju
Współrzędne y, z
sin (11,95^o) =0,2071
cos (11,95^o) =0,9783
y_{(1)} = (20- 3,47) \cdot 0,9783 +(0-12,64) \cdot 0,2071 = 13,55 \, cm,
z_{(1)} = (20- 3,47) \cdot 0,2071 – (0-12,64) \cdot 0,9783 = – 15,79 \, cm,
y_{(2)} = (0- 3,47) \cdot 0,9783 +(0-12,64) \cdot 0,2071 = – 6,01 \, cm,
z_{(2)} = (0- 3,47) \cdot 0,2071 – (0-12,64) \cdot 0,9783 = 11,65 \, cm,
y_{(3)} = (0- 3,47) \cdot 0,9783 +(30-12,64) \cdot 0,2071 = 0,20 \, cm,
z_{(3)} = (0- 3,47) \cdot 0,2071 – (30-12,64) \cdot 0,9783 = = -17,7 \, cm,
y_{(4)} = (13 – 3,47) \cdot 0,9783 +(30-12,64) \cdot 0,2071 = – 12,92 \, cm,
z_{(4)} = (13 – 3,47) \cdot 0,2071 – (30-12,64) \cdot 0,9783 = =-15,01 \, cm,
Pomocnicze współrzędne wycinkowe
\omega_{B, (1)} =0, \omega{B, (2)} = 0, \omega_{B.(3)}=0, \omega_{B, (4)} =- (2 \cdot 1/2 \cdot 30 \cdot 13 =- 380
Momenty bezwładności v z mnożenia wykresów
Ponieważ odcinki konturu przekroju są ukośne w stosunku do układu współrzędnych , więc stosuje się wzory trapezów.
I_y =\\ \cfrac{20}{6} \left [ 2 \cdot \left ( 15,79^2 +1,65^2 \right )+ 2 \cdot 15,79 \cdot 11,65 \right ] \cdot 2,0 + \cfrac{29}{6} \left [ 2 \cdot \left ( 11,65^2 + (-17,70)^2 \right ) + 2 \cdot 11,65 \cdot (-17,70 ) \right ] \cdot 3,0 + \\ \cfrac{13}{6} \left [ 2 \cdot \left ( (-17,10)^2 + (-15,01)^2 \right ) + 2 \cdot (-17,70 ) \cdot (-15,01) \right ] \cdot 1,3 = 10402 \, cm^4
I_y =\\ \cfrac{20}{6} \left [ 2 \cdot \left ( 13,35^2 + (-6.01)^2 \right ) + 2 \cdot 13,35 \cdot (-6,01) \right ] \cdot 2,0 +\\ \cfrac{29}{6} \left [ 2 \cdot \left ( (-6,01) ^2 + (0,20)^2 \right ) + 2 \cdot (-6,01) \cdot 0,20 \right ] \cdot 3,0 + \\ \cfrac{13}{6} \left [ 2 \cdot \left ( 0,20)^2 + 12,92^2 \right ) + 2 \cdot 0,20 \cdot 12,92 \right ] \cdot 1,3 = \\ 3849 \, cm^4
(\ref{23}) \to I_v= \cfrac{1}{3} (20 \cdot 2,0^3+ 30 \cdot 3,0^3 + 13 \cdot 1,4^3 )= 332,9 \, cm^4
Przykłąd 3 [ Charakterystyki geometryczne przekroju zamkniętego]

Rys.12 Współrzędne wycinkowe oraz wycinkowe momenty statyczne w rurze prostokątnej [8]
Wykresy uogólnionych współrzędnych wycinkowych oraz wycinkowych momentów statycznych dla zamkniętego przekroju prostokątnego o stałej grubości ścianki t=const podano na rys.7.
Dla przekroju prostokątnego o dwóch różnych grubościach ścianek (t1,t2) znajduje się wg wzoru w tabeli 1.
Przykład 4 [Siły przekrojowe w pręcie cienkościennym]
Przykład 1-2 z pracy [8]
Znaleźć kąt skręcenia i siły przekrojowe w pręcie o schemacie statycznym podanym w wrs 6. tab.4, ale dla obciążenia skupionym momentem M_s przyłożonego do dowolnego miejsca po długości pręta, jak pokazano na rys. P2-1

Rys. P2-1 Schemat belki do przykładu [8]
Rozwiązanie belki znajdujemy w postaci (\ref{44}).
Z warunków brzegowych dla x=0 : \varphi =0 ; B_\omega=0 wynika wprost C_1=0 ; C_3=0 , a stąd kąt obrotu w przedziale (I) (0 le x < a) opisuje zależność
$\varphi_I = C_2 \cfrac{1}{k sinh (kx)} +C_4 \cfrac{1}{k^3 \overline E I_\omega} [kx – sinh (kx)]
Literatura
- Sapountzakis, E. J., Bars under Torsional Loading: A Generalized Beam Theory Approach. Review Article. ISRN Civil Engineering, 2013, pp.1–39, [ http://dx.doi.org/10.1155/2013/916581]
- Marguerre M., K., Torsion von voll-und hohlquerschnitten, Der Bauingenieur, vol. 21, pp. 317–322, 1940
- Vlasov, V. Z., Tonkostiennyje uprugije stierzni, Gos. Izd. Fiz Mat. Literat., Moskva 1959 / Thin-Walled Elastic Beams. Israel Program for Scientific Translations, Jerusalem, Israel, 1961
- STRENCO. (2016, 2019). Consteel, Oprogramowanie dla budownictwa, [http://www.strenco.pl/strenco ]
- Consteel Software. (2020), ConSteel 14 Manual , [ http://www.consteelsoftware.com/en/downloads/manuals-documents ]
- CTICM. (2013). LTBeam Software (1.011), [ http://www.cesdb.com/ltbeam.html ]
- Birger I., A., Panovko J., G., Spravochnik. Prochnost, ustojchivost, kolebania, Moschinostrojenie, Moskva 1968
- Niewiadomski Z., Wytrzymałość prętów cienkościennych – w pracy Bogucki W. (Red.), Poradnik projektanta konstrukcji metalowych: praca zbiorowa. Tom I., cz.1.1. Arkady, Warszawa), 1980
- Niezgodziński M., E., Niezgodziński T., Wzory, wykresy i tablice obliczeń wytrzymałościowych , IX Wyd. WNT 2004
- Piechnik S,. Wytrzymałość materiałów dla wydziałów budowlanych, PWN, Warszawa-Kraków 1980
- Timoshenko S., Goodier J.N.: Teoria sprężystości, Arkady, Warszawa 1962
- Birger, I. A., & Panowko, J. G.,Procnost, ustojcivost, kolebanija, Spravocnik (Vol. 1). Mosinostroejnije, 1968
- Piechnik, S. Pręty cienkościenne otwarte. Podręcznik dla studentów wyższych szkół technicznych. Wydawnictwo Politechniki Krakowskiej, Kraków 2000
- Piechnik S. Mechanika techniczna ciała stałego, Wydawnictwo Politechniki Krakowskiej, Kraków, 2007
- Rykaluk K., Zagadnienia stateczności konstrukcji metalowych, Dolnośląskie Wydawnictwo Edukacyjne , Wrocław 2012
- Piechnik S., Wytrzymałość materiałów dla wydziałów budowlanych, PWN, Warszawa-Kraków 1980
- (Niewiadomski Z., Wytrzymałość prętów cienkościennych – w pracy Bogucki W. (Red.), Poradnik projektanta konstrukcji metalowych: praca zbiorowa. Tom I., cz.1.1. Arkady, Warszawa), 1980
- Piechnik, S. (2006). Algebra of systems of forces applied to the flat material line. JTAM, 44(1), 107–125
- Clifford, W. K. (1873). Preliminary sketch of bi-quaternions. Proc. London Math. Soc., 4, 381–395
- Jancewicz, B., & Brzeski, P. (2005). Magnetic field surfaces. European Journal of Physics, 26, 617–634.
- Perwass, C. (2009). Geometric algebra with applications in engineering. Springer
- Gull, S., Lasenby, A., , Doran, C. (1993). Imaginary Numbers are not Real – the Geometric Algebra of Spacetime. Found. Phys., 23(9), 1175–1201
- Jancewicz, B. (2011). Pseudowektory. Foton, 115(Zima 2011), 31–43
- Gawłowski, S. (2006). Podejście statyczne do oceny nośności granicznej prętów cienkościennych otwartych, Praca doktorska, Politechnika Krakowska
- Brzoska Z., Statyka i stateczność konstrukcji prętowych cienkościennych, (Wyd. II, PWN, Warszawa, 1965
- Kasprzyk T., Zadanie 5- rozwiązanie przykładowe, Politechnika wrocławska KMBiIM ZWM, Wrocław maj 2021
________________________________