Galīgo elementu metode

Galīgo elementu metode Galīgo elementu metode

04.02.2015 Views

66 NODAĻA 4. DIVDIMENSIJU PROBLĒMAS Tā kā ∫ e dA = A e,kurA e ir elementa laukums, tad jeb U e = 1 2 qT B T DBqtA e (4.40) U e = 1 2 qT K e q (4.41) kur K e ir elementa stinguma matrica K e = B T DBqtA e (4.42) Tā kā matrica D ir simetriska, arī elementa stinguma matrica K e ir simetriska. Visu elementu montāžu apgabalā veicsaskaņā arZīm. 4.2 un Tabulu 4.2. Tādā veidāmēs nosakām visa apgabala deformācijas potenciālo enerǵiju U un globālo stinguma matricu K kur U = ∑ e 1 2 qT K e q = 1 2 QT KQ (4.43) K = ∑ e K e (4.44) Globālā stinguma matrica ir lentveida un retināta matrica. 4.3.4 Elementa mezglu spēki Vispirms nosakām tilpuma spēkus, kas ietilpst funkcionālī (4.34) ∫ ∫ u T ftdA = t (uf x + vf y )dA (4.45) e e Izmantojot aproksimācijas (4.14), iegūstam ∫ ( ∫ ) ∫ u T ftdA = q 1 tf x N 1 dA + q 2 (tf y e e ( ∫ ) ∫ ) +q 3 tf x N 2 dA + q 4 (tf y N 2 dA e e ( ∫ ) ∫ ) +q 5 tf x N 3 dA + q 6 (tf y e N 3 dA e N 1 dA e ) (4.46) Izmantojot elementa formas funkciju definīciju, kas parādīta Zīm. 4.7, redzams, ka ∫ e N 1dA ir piramīdas tilpums. Tās pamats ir A e un augstums ir

4.3. TRĪSSTŪRA GALĪGAIS ELEMENTS 67 1 (bezdimensiju lielums). Tātadpiramīdas tilpums jebkurai formas funkcijai N i ir ∫ N i dA = 1 e 3 A e Tāpēc formulu (4.46) var uzrakstīt ∫ = u T ftdA = q T F e (4.47) e kur F e ir elementa tilpuma spēku, kas reducēti uz mezglu punktiem, vektors F T e = tA e 3 [f x,f y ,f x ,f y ,f x ,f y ] (4.48) Šie mezglu spēki tiek piesummēti globālam slodzes vektoram F , izmantojot elementu savienojumu matricu, kas attēlota tabulā 4.2 F = ∑ e F e Spēku summēšanas procedūra ir līdzīga kā stienim un sijai. ξ =0 ξ =1 h =1 2 η =1 dA dξ A e N 1 1 dη η =0 3 Zīmējums 4.7: Integrālis no elementa formas funkcijas Virsmas spēki T ir izkliedēta slodze uz apgabala kontūru. Šādi spēki darbojas uz elementa malas, kas savieno mezglus uz kontūra (sk. Zīm. 4.8). Šie spēki arī dod papildinājumu globālam mezglu spēku vektoram F .

4.3.<br />

TRĪSSTŪRA GALĪGAIS ELEMENTS 67<br />

1 (bezdimensiju lielums). Tātadpiramīdas tilpums jebkurai formas funkcijai<br />

N i ir<br />

∫<br />

N i dA = 1<br />

e 3 A e<br />

Tāpēc formulu (4.46) var uzrakstīt<br />

∫<br />

= u T ftdA = q T F e (4.47)<br />

e<br />

kur F e ir elementa tilpuma spēku, kas reducēti uz mezglu punktiem, vektors<br />

F T e<br />

= tA e<br />

3 [f x,f y ,f x ,f y ,f x ,f y ] (4.48)<br />

Šie mezglu spēki tiek piesummēti globālam slodzes vektoram F , izmantojot<br />

<strong>elementu</strong> savienojumu matricu, kas attēlota tabulā 4.2<br />

F = ∑ e<br />

F e<br />

Spēku summēšanas procedūra ir līdzīga kā stienim un sijai.<br />

ξ =0 ξ =1<br />

h =1<br />

2<br />

η =1<br />

dA<br />

dξ<br />

A e<br />

N 1<br />

1<br />

dη<br />

η =0<br />

3<br />

Zīmējums 4.7: Integrālis no elementa formas funkcijas<br />

Virsmas spēki T ir izkliedēta slodze uz apgabala kontūru. Šādi spēki<br />

darbojas uz elementa malas, kas savieno mezglus uz kontūra (sk. Zīm. 4.8).<br />

Šie spēki arī dod papildinājumu globālam mezglu spēku vektoram F .

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!