Galīgo elementu metode

Galīgo elementu metode Galīgo elementu metode

04.02.2015 Views

30 NODAĻA 2. VIENDIMENSIJU PROBLĒMAS Tā kā B = 1 [−1 1] l tad K e = AE [ ] −1 [ ] −1 1 (2.41) l 1 jeb sareizinot matricas K e = AE [ ] [ ] 1 −1 K e = 11 K12 e l −1 1 K21 e K22 e (2.42) Redzam, ka elementa stinguma matrica ir simetriska matrica K12 e = K21, e vai vispārīgā veidā Kij e = Kji. e 2.5.2 Elementa ārējie spēki Enerǵētiski ekvivalentos elementa ārējos spēku var noteikt, izmantojot ārējo spēku darba izteiksmi (2.36) un aproksimāciju (2.26) W e = ∫ l 0 u T pdx = ∫ l 0 (N 1 q 1 + N 2 q 2 )pdx (2.43) Šo izteiksmi uzrakstām sekojošā formā ⎡ ∫ l ⎤ W e = [ ] p N 1 dx q 1 q 2 ⎢ 0 ⎣ ∫ l ⎥ ⎦ = [ ] [ ] f q 1 q 1 2 = q T f (2.44) f 2 p N 2 dx 0 kur elementa ārējo spēku vektors ir ⎡ ⎤ [ ] pl f1 ⎢ f = = ⎣ 2 ⎥ f 2 pl ⎦ (2.45) 2 Šī izteiksme iegūta, izskaitļojot šādus integrāļus ∫ l 0 ∫ l 0 N 1 dx = l 2 N 2 dx = l 2 ∫ +1 −1 ∫ +1 −1 1 2 (1 − ξ)dξ = l 2 1 2 (1 + ξ)dξ = l 2 (2.46) Tipiskais galīgais elements ar mezglu spēkiem un mezglu pārvietojumiem ir attēlots Zīm. 2.12

2.6. GALĪGO ELEMENTU MONTĀŽA KONSTRUKCIJĀ 31 1 e 2 ξ q 1 ,f 1 q 2 ,f 2 Zīmējums 2.12: Galīgais elements ar mezglu spēkiem un mezglu pārvietojumiem 2.6 Galīgo elementu montāža konstrukcijā Aplūkosim galīgo elementu montāžu konstrukcijā, kas attēlota Zīm. 2.4. Šai galīgo elementu montāžai var uzrakstīt tā saucamo savienojumu matricu. Tā attēlota Tabulā 2.6. Tabula 2.1: Galīgo elementu savienojumu matrica Elements Mezgli Numerācija e 1 2 Lokālie numuri 1 1 2 Globālie numuri 2 2 3 Globālie numuri Šī savienojumu matrica nosaka, kādā veidā galīgie elementi ir samontēti konstrukcijā. Lielai konstrukcijai šo savienojumu matricu neveido ar roku, bet izmantojot speciālas ǵeometrisko objektu modelēšanas programmas. Ievietojot formulā(2.34) elementa deformācijas enerǵijas (2.37) un ārējo spēku darba (2.44) izteiksmes, iegūstam Π= ∑ Π e = ∑ 1 2 qT K e q − ∑ q T f (2.47) e e e Šeit visas konstrukcijas potenciālā enerǵija ir izteikta ar katra elementa mezglu pārvietojumiem q, kas ir tie paši visas konstrukcijas mezglu pārvietojumi Q. Tikai q tiek definēti katram elementam atsevišķi, bet Q tiek definēti visai konstrukcijai. Sakarību starp elementa pārvietojumiem q un visas konstrukcijas pārvietojumiem Q nosaka Tabulā 2.6definētā galīgo elementu savienojumu matrica (sk. Zīm. 2.5). Tā elementam 1 vektors q T =[Q 1 ,Q 2 ], bet elementam 2 vektors q T =[Q 2 ,Q 3 ]. Tāpēc visas konstrukcijas potenciālo enerǵiju (2.47) var izteikt ar visas konstrukcijas parametriem Π= 1 2 QT KQ − Q T F (2.48)

30 NODAĻA 2. VIENDIMENSIJU PROBLĒMAS<br />

Tā kā<br />

B = 1 [−1 1]<br />

l<br />

tad<br />

K e = AE [ ] −1 [ ]<br />

−1 1 (2.41)<br />

l 1<br />

jeb sareizinot matricas<br />

K e = AE [ ] [ ]<br />

1 −1 K<br />

e<br />

= 11 K12<br />

e<br />

l −1 1 K21 e K22<br />

e (2.42)<br />

Redzam, ka elementa stinguma matrica ir simetriska matrica K12 e = K21, e vai<br />

vispārīgā veidā Kij e = Kji.<br />

e<br />

2.5.2 Elementa ārējie spēki<br />

Enerǵētiski ekvivalentos elementa ārējos spēku var noteikt, izmantojot ārējo<br />

spēku darba izteiksmi (2.36) un aproksimāciju (2.26)<br />

W e =<br />

∫ l<br />

0<br />

u T pdx =<br />

∫ l<br />

0<br />

(N 1 q 1 + N 2 q 2 )pdx (2.43)<br />

Šo izteiksmi uzrakstām sekojošā formā<br />

⎡ ∫ l<br />

⎤<br />

W e = [ ]<br />

p N 1 dx<br />

q 1 q 2<br />

⎢ 0<br />

⎣<br />

∫ l<br />

⎥<br />

⎦ = [ ] [ ]<br />

f<br />

q 1 q 1<br />

2 = q T f (2.44)<br />

f 2<br />

p N 2 dx<br />

0<br />

kur elementa ārējo spēku vektors ir<br />

⎡ ⎤<br />

[ ] pl<br />

f1 ⎢<br />

f = = ⎣ 2 ⎥<br />

f 2<br />

pl ⎦ (2.45)<br />

2<br />

Šī izteiksme iegūta, izskaitļojot šādus integrāļus<br />

∫ l<br />

0<br />

∫ l<br />

0<br />

N 1 dx = l 2<br />

N 2 dx = l 2<br />

∫ +1<br />

−1<br />

∫ +1<br />

−1<br />

1<br />

2 (1 − ξ)dξ = l 2<br />

1<br />

2 (1 + ξ)dξ = l 2<br />

(2.46)<br />

Tipiskais galīgais elements ar mezglu spēkiem un mezglu pārvietojumiem<br />

ir attēlots Zīm. 2.12

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

Saved successfully!

Ooh no, something went wrong!