A B C D E F G H I K Ł M N O P R S T U W Z

Nieliniowy geometrycznie pręt o przekroju zwartym

Niniejszy artykuł jest Dodatkiem C podręcznika Imperfekcyjna metoda projektowania konstrukcji

Publikacja internetowa w wersji „free” z nieograniczonym prawem cytatu – z powołaniem się na autora i źródło:
Leszek Chodor, (2019), Imperfekcyjna metoda  projektowania konstrukcji, Encyklopedia  πWiki,
[ https://chodor-projekt.net/encyclopedia/imperfekcyjna-metoda-projektowania-konstrukcji/ ]

Historia edycji:
(2019-05-02)  Dodatek C : Neiliniowy geometrycznie pręt o przekroju zwartym Wersja 1.0 (2019-05-02)

Proszę społeczność Inżynierów w internecie o przesyłanie recenzji podręcznika  na adres  wydawnictwa biuro@chodor-projekt.net
Leszek Chodor


Przedstawiono geometrycznie nieliniową metodę elementów skończonych prostych prętów o przekroju zwartym, użyteczną przy analizie wpływu rzędu nieliniowego na efekt P-Delta

Pręt nieliniowy geometrycznie

Wyznaczymy ścisłą, nieliniową macierz sztywności oraz strzałkę ugięcia elementu prętowego przekroju zwartym, ściskanego siłą osiową P. Element [e] z węzłami (1)-(2) pokazano na Rys. C.1.

Rys. C.1. Element prętowy ściskany siłą P

Ścisła macierz sztywności

Nieliniową macierz sztywności ściskanego elementu prętowego wyznaczymy w drodze rozwiązania równania różniczkowego ($\ref{C.1}$) belki-słupa bez obciążenia poprzecznego:

$$\begin{equation} EI w^{iv} +Pw^{i} =0 \label {C.1} \end {equation} $$

Wprowadzamy zmienną  $\xi=\cfrac{x}{L}$ oraz miarę obciążenia $\lambda^2=\cfrac{PL^2}{EI}=\pi^2 \alpha_{cr} \to \lambda=\pi \sqrt{\alpha_{cr}}$ , gdzie $\alpha_{cr}= \cfrac{P}{P_{cr}}$  (1.9) dla obciążenia krytycznego Eulera $P_{cr}=\pi^2 EI/L^2$. W granicznym stanie krytycznym $\alpha_{cr}=1$ i stopień wytężenia siła ściskającą  $\lambda=\pi$.

Równanie pręta ($\ref{C.1}$) można teraz przepisać do postaci

$$\begin{equation} \cfrac{d^4 w}{d \xi^2}+ \lambda^2 \cfrac{d^2 w}{d \xi^2}=0 \label {C.2} \end {equation} $$

Rozwiązanie równania ($\ref{C.2}$) można zapisać w postaci:

$$\begin{equation} w(\xi) = C_1 +C_2 \lambda \xi +C_3 \cos{( \lambda \xi)} + C_4 \sin{( \lambda \xi)} \label {C.3} \end {equation} $$

Po wykonaniu różniczkowania funkcji ugięcia uzyskamy funkcję obrotów

$$ \begin{equation} \cfrac{ d w ( \xi  )} { d \xi } = \lambda [ C_2 + C_3 \sin{ ( \lambda \xi) } + C_4  \cos { ( \lambda \xi )} ]  \label {C.3a} \end {equation} $$

Warunki brzegowe wynoszą:

$$\begin{equation}  w(0)= w_1 ,  \quad w (1)=w_2, \quad \cfrac{d w(0)}{d \xi} = \varphi_1 L, \quad \cfrac{d w(1)}{d \xi} = \varphi_2 L  \label {C.4} \end {equation} $$

Warunki brzegowe  ($\ref{C.4}$) można zapisać w formie macierzowej:

$$\begin{equation}  [D] \cdot \left [ \begin {array}{c} C_1 \\ C_2 \\ C_3 \\ C_4 \end {array} \right ]  =  \left [ \begin {array}{c} w_1 \\ \varphi_1 L \\ w_2 \\ \varphi_2 L \end {array} \right] , \text {gdzie  } [D]=
\left [ \begin {array} {ccc} 1 & 0 &1 &0  \\ 0 & \lambda & 0 & \lambda \\ 1 & \lambda &  \cos {\lambda} & {\sin \lambda} \\ 0 & \lambda & -\lambda \sin{ \lambda} & \lambda \cos{ \lambda}  \end {array} \right ] \label {C.5} \end {equation} $$

Po odwróceniu macierzy [D] , otrzymamy:

$$\begin{equation}  [D]^{-1} =\cfrac{1}{\Theta} \left [ \begin {array}{cccc}
\lambda \beta & \gamma & \lambda (1- \cos{ \lambda})  & – ( \lambda – \sin{ \lambda} ) \\
\lambda \sin{ \lambda} & 1- \cos{\lambda} & -\lambda \sin{\lambda} & 1- \cos { \lambda} \\
\lambda (1- \cos{ \lambda}) & -\gamma & -\lambda(1-\cos{ \lambda}) & \lambda -\sin {\lambda} \\
– \lambda \sin{\lambda}  & \beta & \lambda \sin { \lambda}  & -(1-\cos{ \lambda})\\
\end {array} \right] \label {C.6} \end {equation} $$

gdzie:
$\Theta=2 \lambda (1-\cos {\lambda})- \lambda^2 \sin{\lambda}$,
$ \beta= 1- \cos{\lambda } -\sin{\lambda }$,
$\gamma=\lambda\cos{\lambda }-\sin{\lambda }$.

Stałe całkowania wynoszą

$$\begin{equation}  \left [ \begin {array} {c}  C_1 \\ C_2 \\ C_3 \\ C_4 \end {array} \right ]  = [D]^{-1}
\left [ \begin {array}{c} w_1 \\ \varphi_1 L \\ w_2 \\ \varphi_2 L \end {array} \right ] =
\cfrac{1}{\Theta} \left [ \begin {array} {c}
\lambda (\Sigma w -\varphi_2 L)+ \lambda(\varphi_1 L -\Sigma w ) \cos{\lambda}- (\lambda^2 w_1 -\Delta \varphi L)  \sin{\lambda} \\
\Sigma \varphi L \cdot (1-\cos{\lambda}) -\lambda \Delta w \sin {\lambda} \\
\lambda(\varphi_2L -\Delta w ) -\lambda(\varphi_1 L -\Delta w) \cos{\lambda} -\Delta \varphi L \sin {\lambda}  \\
-\Delta \varphi L \cdot (1-\cos{\lambda}) -\lambda (\varphi_1 L -\Delta w ) \sin {\lambda} \\
\end {array} \right ]  \label {C.7} \end {equation} $$

Siły przekrojowe uzyskamy z zależności:

$$\begin{equation} M=EI w^{ii} = – \cfrac{EI}{L^2} \cfrac{dw^2}{d\xi}= \lambda^2 \cfrac{EI}{L^2} (C_3 \cos{\lambda \xi} +C_4 \sin{\lambda \xi}) \label {C.8} \end {equation} $$

$$\begin{equation} T=M^{i} =  – \cfrac{EI}{L^3} \cfrac{dw^3}{d \xi} = \lambda^3 \cfrac{EI}{L^3} ( – C_3 \sin{\lambda \xi} +C_4 \cos{\lambda \xi}) \label {C.9} \end {equation} $$

Wektor sił węzłowych wynosi

$$\begin{equation}  \left [  \begin {array} {c}  M_1 \\ T_1 \\ M_2 \\ T_2 \end {array} \right ] =
\left [  \begin {array} {c}  M(\xi=0) \\ T(0) \\ M(1)\\ T(1) \end {array} \right ] =
\lambda^2 \cfrac{EI}{L^2} \cdot \left [  \begin {array} {cc} 1 & 0 \\  0 & \cfrac{\lambda}{L}  \\  \cos{\lambda} & \sin {\lambda} \\ \cfrac {- \lambda}{ L \sin{\lambda}} & \cfrac{ \lambda}{ L \cos{\lambda}}
\end {array} \right ]
\cdot \left [  \begin {array} {c}  C_3\\ C_4 \end {array} \right ]
\label {C.10} \end {equation} $$

Po przekształceniach powyższych wyrażeń (pokazano w pracy [1]), uzyskamy macierz sztywności nieliniowego geometrycznie pręta ściskanego w postaci:

$$\begin{equation} [k]^{[e]} =  \left [  \begin {array} {cccccc}
\cfrac{EA}{L} & 0& 0 & \cfrac{-EA}{L} &0 &0 \\
& \cfrac{12 EI}{L^3} \Lambda_1 & \cfrac{6EI}{L^2} \Lambda_2 & 0 & \cfrac{-12 EI}{L^3} \Lambda_1 & \cfrac{6EI}{L^2} \Lambda_2 \\
& & \cfrac{4EI}{L} \Lambda_3 & 0 & \cfrac{-6EI}{L^2} \Lambda _2& \cfrac{2EI}{L} \Lambda_4 \\
& & & \cfrac{EA}{L} & 0 & 0  & \\
& SYM & & &  \cfrac{12EI}{L^3} \Lambda_1 & \cfrac{-6 EI}{L^2} \Lambda_2 \\
& &  & & & \cfrac{4 EI}{L} \Lambda_3 \\
\end {array} \right ]  \label {C.11} \end {equation} $$

Macierz ($\ref{C.11}$) zawiera konwencjonalne współczynniki pręta zginanego zmodyfikowane współczynnikami $\Lambda_i \quad (i=1,..,4)$ [2], [1].
Współczynniki $\Lambda_1$ do $\Lambda_4$ definiują ścisłą geometrycznie nieliniową macierz sztywności pręta ściskanego, zależą od wielkości siły osiowej N w elemencie, wywołanej siłami P na  konstrukcji  na Rys. C.1) i wynoszą:

$$\begin{equation} \begin {cases}
\Lambda_1= & \alpha_c \cdot \Lambda_2 \\
\Lambda_2= & \cfrac{\alpha^2}{3(1-\alpha_c) } \\
\Lambda_3= & \cfrac{3}{4} \Lambda_2+\cfrac {\alpha_c}{4} \\
\Lambda_4= & \cfrac{3}{2} \Lambda_2 -\cfrac{\alpha_c}{2} \\
\end{cases} \label {C.12} \end {equation} $$

gdzie $\alpha_c=\alpha \cdot ctg{\alpha}$,

Współczynnik $\alpha$ wynosi

$$\begin{equation} \alpha= \cfrac{\lambda}{2}= \cfrac{\pi}{2} \sqrt{\cfrac{N}{N_{cr}}} \label {C.13} \end {equation} $$

W przypadku znajomości siły osiowej w elemencie, współczynniki ($\ref{C.12}$), pozwalają wyznaczyć ścisłą sztywność elementu i w rezultacie nieliniową sztywność systemu konstrukcyjnego. W prostych przypadkach metoda jest użyteczna do obliczeń symbolicznych oraz porównawczych.

W przypadku konwencjonalnych obliczeń MES , polegających na rozwiązywaniu liniowego układu równań współczynniki $\Lambda_i$ wywołują nieliniowość i potrzebę iteracyjnego rozwiązywania systemu konstrukcyjnego. Do tego celu użyteczne jest rozwinięcie funkcji  $\Lambda_i \quad (i=1,2,3,4)$ w szereg Taylora podług potęg $\cfrac{NL^2}{EI}$:

$$\begin{equation} \Lambda_i \approx \sum \limits_{i=0}^{i=4} C_{1,i}  \cdot  \left( \cfrac{NL^2}{EI} \right)^i \label {C.14} \end {equation} $$

Współczynniki rozwinięcia zestawiono w Tab. C.1. Rozwinięcie ścisłej macierzy sztywności w szereg potęgowy podług potęg , jest stosowane w klasycznej analizie mechaniki. Człony liniowe są zachowywane w klasycznej, liniowej analizie wyboczenia (LBA), a mnożniki w kolumnie (1) Tab. B.1. określają tak zwaną macierz geometryczną. W teorii 2 rzędu (efekt P-Δ) są zachowywane człony wyższych rzedów  $(\cfrac{NL^2}{EI})^2 $, itd.

Tab. C.1. Współczynniki $C_{1,i}$ rozwinięcia funkcji w szereg Taylora podług potęg $\cfrac{NL^2}{EI}$

W analizach wpływu stopnia nieliniowości teorii na efekty, bardziej użyteczne jest rozwinięcie funkcji szereg podług potęg współczynnika ($\ref{C.10}$) :

$$\begin{equation} \Lambda_i \approx \sum \limits_{i=0}^{i=4} C_{2,i}  \cdot  \alpha^i \label {C.15} \end {equation} $$

Tab. C.2. Współczynniki $C_{2,i}$ rozwinięcia funkcji w szereg Taylora podług potęg $\alpha$

Strzałka ugięcia prostego pręta ściskanego

Z zależności ($\ref{C.2}$) po podstawieniu uzyskanych wyrażeń na stałe całkowania można wyznaczyć funkcję ugięcia pręta pomiędzy węzłami w zależności od przemieszczeń węzłowych:

$$\begin{equation}  w(\xi)=\cfrac{1}{\Theta}
\left [  \begin {array} {c}
1- \cos{\lambda} +\cos{ \lambda \xi} -\cos{[ \lambda(1-\xi)] } – \lambda \sin { [\lambda(1-\xi)] }  \\
– \lambda \xi – \lambda (1-\xi) \cos{\lambda}+\xi\cos{ [ \lambda(1-\xi)] } +\sin{\lambda} +\sin{[ \lambda\xi] } \sin{[ \lambda(1-\xi)] } \\
1- \cos{\lambda} -\cos{\lambda \xi} +\cos{[ \lambda(1-\xi)] } -\lambda \xi \sin \lambda \\
– \lambda (1- \xi ) – \lambda \xi \cos{\lambda} +\lambda \cos{\lambda \xi } +\sin{\lambda} – \sin{\lambda \xi } -\sin { [\lambda(1-\xi ) ]} \\
\end {array} \right ]
\cdot  \left [  \begin {array} {cc} \lambda w_1 & \varphi_1 L & \lambda w_2 & \varphi_2 L\end {array} \right ]
\label {C.16} \end {equation} $$

Z formuły ($\ref{C.16}$) wynika, że pręt idealnie prosty, ugina się pod wpływem osiowej siły ściskającej w każdym przypadku, gdy węzły doznają przemieszczeń liniowych lub obrotów wstępnych, czyli w każdej konstrukcji rzeczywistej.

W stanie granicznym krytycznym, dla  $\lambda= \pi$ mamy:

$$\begin{equation}  w(\xi)=\cfrac{1}{2 \pi} \left [  \Sigma w +\Sigma \varphi L  (\xi-1/2) + \cfrac{ \pi}{2} [ \Sigma \varphi L  -2 \Delta w ] \cos{ ( \pi \xi)}- \Delta \varphi L sin (\pi \xi) \right ]   \label {C.17} \end {equation} $$

Maksymalne ugięcie pręta, czyli strzałka wygięcia wystąpi  w miejscu pręta o współrzędnej $\xi_{sup}$

$$\begin{equation}  \xi_{sup} = \pm \cfrac{1}{\lambda} \cdot arcos {(\pm \Gamma)} \label {C.18} \end {equation} $$

gdzie $ \Gamma = \cfrac {1} { \sqrt{ (C_3/C_4)^2 +1}} $. Wyrażenie uzyskano w drodze porównania do zera drugiej pochodnej funkcji $\ref{C.3}$), a w rezultacie uzyskano  cztery pierwiastki (przy wszystkich kombinacjach znaków).

Po odrzuceniu pierwiastków ujemnych jako fizycznie niemożliwych i po podstawieniu dodatniej współrzędnej ($\ref{C.18}$) do ($\ref{C.16}$) otrzymamy strzałkę ugięcia pręta w zalezności od warunków brzegowych wyrażonych stałymi $C_3 , C_4$.

W przypadku symetrycznych warunków brzegowych strzałka ugięcia wystąpi w środku pręta i wyniesie

$$\begin{equation} f=w (\xi =0,5 ) = \cfrac{1}{2} \left ( \Sigma w -\cfrac { tg(\lambda/4) \cdot \Delta \varphi L }{\lambda}\  \right)  \label {C.19} \end {equation} $$

W granicznym stanie krytycznym: $ \lambda=\pi$ i strzałka ugięcia $(\ref{C.19}$) osiągnie wartość:

$$\begin{equation} f_{cr} = \cfrac{\Sigma w}{2} – \cfrac{\Delta \varphi L}{2 \pi}  \label {C.20} \end {equation} $$

Dla $N=0$ jest $\lambda=0$, więc  $f_0 = \lim \limits_{\lambda \to 0}= (\Sigma w) /2 -(\Delta \varphi L) / 8 $, co jest zgodne z rezultatem teorii I rzędu.

Wyrażenia  $(\ref{C.17})$ i $(\ref{C.20})$ są  fundamentalnymi wynikami niniejszego rozdziału. Wskazują ono, że w granicznym stanie  krytycznym (Eulera):

  1. Funkcja ugięcia m postać $(\ref{C.17})$, to znaczy jest mieszaną funkcją liniową z sinusoidą oraz kosinusoidą, można więc ją sprowadzić do szeregu Fouriera
  2. Strzałka ugięcia pręta jest skończona i niezerowa ,  na skutek przemieszczeń i obrotów przekrojów na podporach.

W dalszej części dodatku C rozważymy kilka przypadków szczegółowych.

Nieliniowy geometrycznie pręt wspornikowy

Ściskany pręt wspornikowy obciążony jest poprzecznie siłą H w wierzchołku (1) lub obciążeniem równomiernie rozłożonym h po wysokości. Pręt jest sprężyście utwierdzony w węźle dolnym (1) ze stałą sprężystości $c_\varphi$. -rys. C.2.

Warunki brzegowe można zapisać w postaci:

$$\begin{equation} u_2=0, \quad w_2=0, \quad M_2= -C_{\varphi} \cdot \varphi_2 \label {C.21} \end {equation} $$

Rys. C.2. Pręt wspornikowy

W przypadku konstrukcji jednoelementowej macierz ($ \ref{C.11}$) jest macierzą konstrukcji i można ją bezpośrednio modyfikować w celu uwzględnienia warunków brzegowych ($ \ref{C.21}$). Modyfikację macierzy sztywności dokonamy poprzez wykreślenie wierszy i kolumn , odpowiadających odebranym stopniom swobody i bezpośrednie wpisanie warunku sprężystej podpory do równania kanonicznego. W rezultacie otrzymamy macierz zmodyfikowaną $[\widetilde{K]}$:

$$\begin{equation} [ \widetilde{K]} = \cfrac{EI}{L}
\begin{bmatrix}
\cfrac{EA}{EI} & 0 & 0 & 0 \\
& \cfrac{12 }{L^2} \Lambda_1 & \cfrac{6}{L} \Lambda_2 & \cfrac{6}{L} \Lambda_2 \\
& &  4 \Lambda_3  & 2 \Lambda_4 \\
& SYM & & 4 \Lambda_3+ c_{\varphi} \\
\end{bmatrix}
\begin{matrix}
u_2 \\ w_2 \\ \varphi_2 \\ \varphi_1
\end {matrix}
\label {C.22} \end {equation} $$

Poprzez odwrócenie macierzy ($\refC.22}$)  uzyskujemy macierz podatności nieliniowego geometrycznie wspornika:

$$\begin{equation} [ \widetilde{K]}^{-1} = \cfrac{L}{\Theta}
\begin{bmatrix}
\cfrac{\Theta}{EA} & 0 & 0 & 0 \\
& k_{22} & k_{24}+\Lambda_2 \cdot c_{\varphi}\cdot L/2 & k_{24} \\
& &  3 \Lambda_2^2 -\Lambda_1(4 \Lambda_3+c_{\varphi})  & -3\Lambda_2^2 +2 \Lambda_1 \Lambda_4 \\
& SYM & & 3 \Lambda_2^2 -4 ] \Lambda_1 \Lambda_3 \\
\end{bmatrix}
\begin{matrix}
u_2 \\ w_2 \\ \varphi_2 \\ \varphi_1
\end {matrix}
\label {C.23} \end {equation} $$

gdzie:
$c_{\varphi}=C_{\varphi} \cdot L/EI$,
$k_{24}=\Lambda_2L(2\Lambda_3-\Lambda_4)$,
$k_{22}=-\cfrac{L^3}{3}\cdot(4\Lambda_3^2-\Lambda_4^2+\Lambda_3c_{\varphi} ) $,
$\Theta = – 4EI \cdot [(2\Lambda_3 -\Lambda_4) \cdot \Lambda_1 (2 \Lambda_3 + \Lambda_4) – 3\Lambda_2^2] + c_{\varphi} \cdot L \cdot  ( 3 \Lambda_2^2 –  4 \Lambda_1 \Lambda_3) $.

Przemieszczenia i siły przekrojowe we wsporniku zostały przeanalizowane w treści podręcznika. Przemieszczenia węzłowe uzyskano poprzez przemnożenie macierzy podatności $(\ref{C.23}$) przez wektor równoważników węzłowych, który w przypadku obciążenia rozłożonego h, wynosi

$$\begin{equation} [F]_h = \cfrac{hL}{2} \cdot [\quad 0 \quad  1 \quad  L/6 \quad -L/6 \quad ] \label {C.24} \end {equation} $$

a w przypadku obciążenia skupionego H, wynosi

$$\begin{equation} [F]_H =  [ \quad 0 \quad  H \quad  0  \quad 0 \quad ] \label {C.25} \end {equation} $$

Moment zginający w utwierdzeniu wspornika otrzymano po dodaniu do odpowiedniej współrzędnej wektora siły przywęzłowej, uzyskanej z przemnożenia wektora przemieszczeń węzłowych przez ostatni wiersz pierwotnej macierzy sztywności.

Nieliniowa geometrycznie belka wolnodparta

Warunki brzegowe dla belki wolnopodpartej (belki-słupa, pokazanego na C.3 są następujące

$$\begin{equation} u_1=0, \quad w_1=0, \quad T_2= -C_{\Delta} \cdot u_2 \label {C.26} \end {equation} $$

gdzie: $C_\Delta$ jest stała sprężystości podparcia w węźle prawym (2).

Eys. C.3 Schemat belka-słup

Belkę potraktujemy, jako jeden element, a strzałkę ugięcia wyznaczymy z zależności wyprowadzonych w poprzednim rozdziale. Zmodyfikowana macierz sztywności belki-słupa wynosi:

$$\begin{equation} [ \widetilde K] = \cfrac{EI}{L}
\begin{bmatrix}
4 \Lambda_3  & 0 & \cfrac {-6 \Lambda_2} {L}  & 2 \Lambda_4 \\
&   \cfrac {EA}{EI} & 0 & 0 \\
&& \cfrac{12 \Lambda_1 + c_{\Delta}} {L^2} & \cfrac { – 6\Lambda_2}{L} \\
& SYM  & \cfrac{ -6 \Lambda_2}{L} & 4 \Lambda_3 \\
\end{bmatrix}
\begin{matrix}
\varphi_1 \\u_2 \\ w_2 \\ \varphi_2
\end {matrix}
\label {C.27} \end {equation} $$

Odwrotność macierzy zmodyfikowanej wynosi:

$$\begin{equation} [ \widetilde K]^{-1} = \cfrac{EI}{L}
\begin{bmatrix}
k_{11}  & 0 & 3 L \cdot \Lambda_2  & k_{14} \\
0 &   \cfrac {\Theta} {EA} & 0 & 0 \\
& & L ( 2 \Lambda_3+\Lambda_4 ) & 3 L \cdot \Lambda_2 \\
& SYM  &  & k_{11} \\
\end{bmatrix}
\begin{matrix}
\varphi_1 \\u_2 \\ w_2 \\ \varphi_2
\end {matrix}
\label {C.28} \end {equation} $$

gdzie:
$c_\Delta= C_\Delta \cdot \cfrac{L^3}{EI}$,
$k_{11}= \cfrac{-9\Lambda_2^2 +(c_\Delta + 12 \Lambda_1 ) \Lambda_3}{2 \Lambda_3 -\Lambda_4}$,
$k_{44}=\cfrac{18 \Lambda_2^2 – (c_\Delta + 12 \Lambda_1 ) \Lambda_4}{4 \Lambda_3 -2 \Lambda_4}$,
$\Theta=EI\cdot [36 \Lambda_2^2 -(12\Lambda_1 +c_\Delta)  \cdot (2 \Lambda_3 +\Lambda_4)]$.

Siły węzłowe od obciążenia rozłożonego q

$$\begin{equation} [F]_q =  \cfrac{qL}{2} [ \quad L/6  \quad 0  \quad 1  \quad – L/6 \quad ] \label {C.29} \end {equation} $$

Z formuły ($\ref {C.19}$)   można wyznaczyć ścisłą strzałkę ugięcia pręta od obciążenia ściskającego $N$ bez udziału sił poprzecznych.

$$\begin{equation} f_N =w_N (\xi=0,5) = \cfrac{qL^4}{48EI} \cdot \left ( \cfrac{12 (2 \Lambda_3 + \Lambda_4)}{36 \Lambda_2^2 +(12 c_\Delta  \Lambda_1) \cdot (2 \Lambda_3 +\Lambda_4)} +\cfrac{\alpha \cdot tg(\alpha/2)}{2 \Lambda_3 +\Lambda_4 } \right) \label {C.30} \end {equation} $$

Strzałką ugięcia pochodząca od sił poprzecznych jest całką szczególną zagadnienia, która powinna być dodana do ($\ref{C.30}$). W przypadku obciążenia równomiernie rozłożonego całka szczególna wynosi [2]:

$$\begin{equation}  w_q (\xi) = \cfrac{qL^4}{24 EI}  \to  f_q= w(\xi=0,5)= \cfrac{q L^4}{384 EI}  \label {C.31} \end {equation} $$

Element krzywoliniowy, sinusoidalny

Macierz podatności elementu, którego oś jest  opisana ogólną zależnością parametryczną

$$\begin{equation}  x=x(s) \quad y=y(s)  \label {C.32} \end {equation} $$

podał [2] w postaci

$$\begin{equation} [  K]^{-1} =
\begin{bmatrix}
\int \limits_0^L \left ( \cfrac{(dx/dy)^2}{EA}+\cfrac{y^2}{EI} \right ) ds & \int \limits_0^L \left ( \cfrac{(dx/dy) (dy/ds)} {EA}- \cfrac{xy}{EI} \right ) ds & \int \limits_0^L \cfrac{y}{EI} ds \\
&  \int \limits_0^L \left ( \cfrac{(dy/ds)^2} {EA}- \cfrac{x^2}{EI} \right ) ds & – \int \limits_0^L \cfrac{x}{EI} ds  \\
& SYM  &  \int \limits_0^L \cfrac{1}{EI} ds \\
\end{bmatrix}
\label {C.33} \end {equation} $$

gdzie: $A(s)$ i $I(s)$ – pole przekroju i moment bezwładności sa znanymi funkcjami współrzędnej bieżącej s

Kształt elementu wstępnie wygiętego w łuk sinusoidy opisany w standardowym układzie współrzędnych (x,y)

$$\begin{equation}  y= \eta \sin{(\omega x)} \label {C.34} \end {equation} $$

gdzie:  $\eta$ – amplituda, $\omega= \cfrac{ n \ pi x}{L}$ – częstotliwość dla $n=1,2,3,…$ – liczba naturalna oznaczająca numer postaci własnej elementu,
Równanie krzywej ($\ref{C.34}$) można zapisać w postaci parametrycznej:

$$\begin{equation}  y= \eta \sin {\alpha} ;  \quad x =\cfrac{\alpha}{\omega} \label {C.35} \end {equation}$$

gdzie parametr $\alpha \in  [ \quad \alpha_1= -\omega L/2 ; \quad \alpha_2= \omega L/2 \quad ]$,

Poszczególne wyrazy macierzy podatności $( \ref{C.32}$) wynoszą:

$k_{11}= \int \limits_0^L \left ( \cfrac{(dx/dy)^2}{EA}+\cfrac{y^2}{EI} \right ) ds=\int \limits_{\alpha_1}^{\alpha_2} \left ( \cfrac{(dx/dy)^2}{EA}+\cfrac{y^2}{EI} \right ) d \alpha= \cfrac{ \eta ^2 [ n \pi – \sin {(n \pi) ] } }{2EI}$,

$k_{12}= \int \limits_{\alpha_1}^{\alpha_2} \left ( \cfrac{(dx/dy) (dy/ds)} {EA}- \cfrac{xy}{EI} \right ) d \alpha = \cfrac{ \eta L [ n \pi \cdot \cos {(n \pi /2) } – 2 \sin{(n \pi / 2)]} }{n \pi \cdot EI }$,

$k_{22}= \int \limits_{\alpha_1}^{\alpha_2} \left ( \cfrac{(dy/ds)^2} {EA}- \cfrac{x^2}{EI} \right ) d \alpha = \cfrac{L}{12} \left (  \cfrac{6 \eta^2 [ n \pi +\sin{(n \pi) ] }}{EA}+\cfrac{L^2 \cdot n \pi}{EI}\right )$,

$k_{13}= \int \limits_{\alpha_1}^{\alpha_2} \cfrac{y}{EI} d \alpha  =0$,

$k_{23}= \int \limits_{\alpha_1}^{\alpha_2} \cfrac{-x}{EI} d \alpha  =0$,

$k_{33}= \int \limits_{\alpha_1}^{\alpha_2} \cfrac{1}{EI} d \alpha  = \cfrac{n \pi}{EI}$,

Macierz sztywności uzyskamy przez odwrócenie macierzy podatności. Stosunkowo proste wyrażenie uzyskamy  dla kształtu pręta proporcjonalnego do  pierwszej postaci własnej, czyli dla n=1:

$$\begin{equation} [ K] = \cfrac{1}{\Theta}
\begin{bmatrix}
2 \pi^3 EI ( EA L^2/ \eta +6EI \cdot \eta & 48  \pi L \cdot EA \cdot EI  & 0 \\
&   12 \pi^3 EA \cdot EI  \cdot \eta & 0  \\
& SYM  & EI \Theta / \pi \\
\end{bmatrix}
\label {C.36} \end {equation} $$

gdzie: $\Theta= \eta [ EA L^2 (\pi^2 -96) +6EI \pi^4 \eta^2 ]$

Bibliografia artykułu
  1. Chodor, L. (2013). Belka Timoshenko na sprężystym-podłożu. PiWiki – Inżynierowie i Architekci Chodor-Projekt, [ https://chodor-projekt.net/encyclopedia/belka-timoshenko-sprezyste-podloze/ ]
  2. Livesley R. K. (1975), Matrix methods of structural analysis, Pergamon Press [http://books.google.com/books?id=tYsoAQAAMAAJ ]
Comments : 0
O autorze
* dr inż. Leszek Chodor. Architekt i Inżynier Konstruktor; Rzeczoznawca budowlany. Autor wielu projektów budowli, w tym nagrodzonych w konkursach krajowych i zagranicznych, a między innymi: projektu wykonawczego konstrukcji budynku głównego Centrum "Manufaktura" w Łodzi, projektu budowlanego konstrukcji budynku PSE w Konstancinie Bielawa, projektów konstrukcji "Cersanit" ( Starachowice, Wałbrzych, Nowograd Wołyński-Ukraina), projektu konstrukcji hali widowiskowo-sportowej Arena Szczecin Autor kilkudziesięciu prac naukowych z zakresu teorii konstrukcji budowlanych, architektury oraz platformy BIM w projektowaniu.

Twój komentarz do artykułu

Translate »