5.
Probleme de elasticitate 3D
5.1. Elementul finit T10
T10 este un element de tip tetraedru cu 10 noduri. Nodurile sunt
situate la col urile tetraedrului i la mijlocul muchiilor drepte.
Coordonatele nodurilor sunt ( )
T
e
x x x
10 2 1
x , ( )
T
e
y y y
10 2 1
y
i ( )
T
e
z z z
10 2 1
z . Caracteristicile materialului izotrop sunt date prin
modulul de elasticitate E i coeficientul de contrac ie transversal .
Gradele de liertate nodale sunt deplasrile u, v i w ale nodurilor,
rezult!nd un numr total de "0 de #rade de liertate pe element. $ectorul
deplasrilor elementului este%
( )
T
e
w w w v v v u u u
10 2 1 10 2 1 10 21 1
a .
&i#ura '.1 (lementul tetraedru cu 10 noduri
'
Formularea elementului
(cua iile elementului rezult printr)o formulare *n deplasri, care
introduce polinoame de apro+imare pentru c!mpul de deplasri u din
interiorul elementului.
P
P 0 0
0 P 0
0 0 P
u
'
1
1
1
]
1
'
"0
2
1
2
2
2
w
v
u
,
'.1
,ici P
2
este matricea termenilor polinomiali de #radul 2 *n x y i z.
[ ] . 1
2 2 2
2
z yz y xz xy x z y x P
'.1
Coeficien ii necunoscu i ( )
T
"0 2 1
treuiesc
determina i *n func ie de deplasrile din nodurile elementului.
-un!nd condi ia ca rela ia '.1 s conduc la deplasrile din
nodurile elementului, *nlocuind coordonalele nodurilor *n '.1, se o ine%
. G a
e '.2
.in rela ia '.2 se pot e+plicita parametrii
e
a G
1
. '."
/nlocuind rela ia '." *n rela ia '.1 se o ine c!mpul de deplasri *n
element%
e e
a a PG u
1
. '.0
,ici este matricea func iilor de interpolare format din func1iile N
i
(x,y),
(i=1,2, ... , 10) corespunztoare nodurilor elementului, care au urmtoarea
proprietate%
2 10 ,..., 2 , 1 , 3
, 0
, 1
2 , 3
'
j i
j i
j i
y x N
j j i .
'.'
&unc iile de interpolare permit determinarea c!mpului de deplasri *n
4
interiorul elementului *n func ie de deplasrile nodurilor. &unc iile de
interpolare sunt polinoame paraolice complete *n x,y i z, iar *n lun#ul
unei muchii deplasrile depind numai de deplasrile nodurilor adiacente.
,cest lucru asi#ur compatiilitatea deplasrilor *ntre elementele vecine de
acela i tip.
-rin aplicarea rela iei #eometrice dintre deplasri i deforma1ii
din teoria elasticit1ii, din c!mpul de deplasri '.1 rezult deforma iile
elementului%
e e
y z
x z
x y
x
x
x
y z
x z
x y
z
y
x
yz
xz
xy
z
y
x
w v
w u
v u
w
v
u
!a a "
P P 0
P 0 P
0 P P
P 0 0
0 P 0
0 0 P
#
1
1
1
1
1
1
1
1
]
1
'
+
+
+
'
1
, 2 , 2
, 2 , 2
, 2 , 2
, 2
, 2
, 2
, ,
, ,
, ,
,
,
,
.
'.4
5atricea ! este matricea deforma ie)deplasare i este o func1ie liniar *n x, y i z,
iar
[ ]
[ ]
[ ] , 2 0 0 0 1 0 0 0
, 0 2 0 0 0 1 0 0
, 0 0 0 2 0 0 1 0
, 2
, 2
, 2
z y x
z y x
z y x
z
y
x
P
P
P
'.4
sunt derivatele matricei P
2
.
Tensiunile din interiorul elementului se o in prin aplicarea rela iei fizice%
# " # $
1
1
2 1 3 2 0 0 0 0 0
0 2 1 3 2 0 0 0 0
0 0 2 1 3 2 0 0 0
0 0 0 1
0 0 0 1
0 0 0 1
1
,
_
1
1
1
1
1
1
1
]
1
+
+
+
'
E
yz
xz
xy
z
y
x
'.
6
,ici " este matricea le#ii lui constitutive pentru elasticitate
tridimensional.
6
(cua iile elementului se deduc pe aza unor rela ii de lucru
mecanic virtual. 7e presupune c deplasrile i deforma iile virtuale
u
i
#
rezult din deplasrile nodale virtuale
a
pe aza acelora i rela ii ca i
cele reale%
e
a u
, e
a ! #
. '.8
-entru e+primarea echilirului sistemului discretizat prin nel elemente,
ac1ionat de for1ele de volum b i *ncrcrile t
9
pe suprafa1a
t
, se poate
scrie urmtoarea rela ie%
+
t
d dV dV
T
V
T
V
T
t u b u $ #
9
'.:
;ela ia 3:2 e+prim condi ia de echiliru su forma L
ef
=L
ext
, care treuie
s fie *ndeplinit pentru orice c!mp virtual deplasri
u
i deforma ii
#
/nlocuind rela iile 302, 342, 362 i 382 *n 3:2, rezult condi iile de
echiliru e+primate *n func ie de vectorul deplasrilor reale i virtuale a i
a
din toate nodurile sistemului discretizat cu elemente finite%
+
nel
T
V
T T
nel
V
T T
el el el
d dV dV 2
9
3 2 3
1
t b a a ! " ! a
. '.10
,v!nd *n vedere c vectorul deplasrilor virtuale
a
treuie s fie
nenul, rela ia 3102 se poate simplifica i rezult%
+
nel
T
V
T
nel
V
T
el el el
d h dV dV 2
9
3 2 3
1
t b a ! " !
, '.11
sau
F %a
, '.12
*n care % este matricea de ri#iditate structural, iar F este vectorul
*ncrcrilor pe structur.
7criind rela iile la nivelul unui element rezult matricea de
ri#iditate a elementului%
8
el
V
T
dV ! " ! &
1
'.1"
i vectorul *ncrcrilor pe element%
+
el el
V
T T
d dV t b f
9
. '.10
/n rela1ia 10, b reprezint vectorul for1elor de volum distriuite, iar t
9
sunt
trac1iunile pe suprafa1.
,naliz!nd rela ia 31"2, se constat c inte#randul are #radul polinomial doi,
astfel rela ia se poate evalua prin cuadratur numeric folosind patru
puncte de inte#rare%
2 , , 3 2 , , 3
0
1
i i i
IV
I i
T
i i i
el
z y x z y x
V
! " ! &
, '.1"
*n care V
el
este volumul elementului
1
1
1
1
]
1
0 0 0
" " "
2 2 2
1 1 1
1
1
1
1
det
4
1
z y x
z y x
z y x
z y x
V
el
, '.1'
iar x
i
, y
i
i z
i
3i=I, II, III, IV2 sunt coordonatele punctelor de inte#rare%
1
1
1
1
]
1
1
1
1
1
]
1
1
1
1
1
]
1
0 0 0
" " "
2 2 2
1 1 1
z y x
z y x
z y x
z y x
a
a
a
a
z y x
z y x
z y x
z y x
IV IV IV
III III III
II II II
I I I
cu a < 0.'8'01020 i < 0.1"81:440.
'.14
5atricea de ri#iditate structural % i vectorul F al *ncrcrilor pe structur
din rela1ia '.12 rezult prin asamlarea matricelor & i vectorilor f ale
elementelor%
:
2. , 1 3
,
, ,
nel e f
e ind ind
e ind ind ind ind
e e
e e e e
+
+
F F
& % %
, '.16
/n rela1ia de mai sus ind
e
reprezint vectorul indicilor de cone+iune
3coresponden1a *ntre cele "0 de deplasri ale elementului i deplasrile
nodurilor structurii2.
;ezolv!nd sistemul de ecua1ii liniare '.12 rezult deplasrile nodurilor%
a<%
)1
F, care permit determinarea eforturilor *n elemente prin rela1ia%
e e
a ! " $
1
. '.18
5.' ()lica*ie+ ,tructuri rutiere multistrat
5odelul de calcul este un semispa1iu *ncrcat cu o presiune
uniform distriuit pe o suprafa1 de contact circular. =in!nd cont de
simetrie se poate lucra pe un sfert de model.
>serva1ie% -rolema are o simetrie a+ial i e+ist i o rezolvare
cu ajutorul unor elemente finite a+ialsimetrice, dar aici se studiaz modelul
". pentru dezvoltri ulterioare, care nu sunt a+ialsimetrice.
10
&i#ura '.2 7tructura rutier multistrat ) modelul de calcul
-entru aplicarea 5(& s)a adoptat o divizare cu ajutorul unei re1ele
re#ulate ca *n fi#ura '.". > unitate de volum astfel rezultat s)a sudivizat *n cinci
elemente de tip teraedru T10. .imensiunile modelului s)au stailit prin *ncercri,
astfel *nc!t diferen1ele fa1 de modelul infinit s fie acceptaile *n calculele
in#inereti. 7uprafa1a circular de raz ! s)a apro+imat cu un octo#on caracterizat
prin dimensiunile d
1
i d
2
, cu d
2
" !.
/n urma calculelor rezult deplasrile nodurilor, deforma1iile i
tensiunile *n elemente. .in punct de vedere practic prezint interes numai c*teva
valori *n puncte caracteristice ale modelului de calcul aa cum sunt prezentate *n
continuare pentru un e+emplu.
-200
0
200
400
600
800
0
200
400
600
-800
-700
-600
-500
-400
-300
-200
-100
0
&i#ura '." 7istem rutier multistrat ) discretizare cu elemente finite
PROGRAMUL S R M 3D - v.06.10
Aplicatia: Exemplu C! "00#
DA$ELE PRO!LEME
% & '(.) *+, -1 & 6..0 mm, -" & 166.' mm, / & 0.6"' MPa
11
St0atul E1Mpa2 3u 4056imea1mm2
1 3600 0.3' )0
" 3000 0.3' '0
3 '000 0.3' #0
) 300 0."( 1'0
' (0 0.)" '00
Divi7a0ea pe -i0ectia x 1mm2
3).0 3'.0 )#.# )#.# 100.0 160.0 "00.0
Divi7a0ea pe -i0ectia 8 1mm2
3).0 3'.0 )#.# )#.# 100.0 160.0 "00.0
Divi7a0ea pe -i0ectia 7 1mm2
)0 '0 #0 (' (' 100 100 1'0 1'0
RE9UL$A$E
Depla6a0ea ve0ticala i3 ax : 0.")( 1mm2
x 8 7 6i4ma x ep6il53 x ep6il53 7
mm mm mm MPa mic05-e: mic05-e:
0 0 )0 -0.'60 -##.30. '#.1(#
0 0 )0 -0.((1 -)'.3"( -1().((#
0 0 .0 -0."01 0..#1 -#3.336
0 0 .0 -0.1#0 ).'03 -'6.#.'
0 0 1(0 1.11) 1)'.)0. -161.6))
0 0 1(0 0.0") 1('.1(( -)#3.)3'
0 0 3"0 0.0'0 1(3.0)) -"#1.300
0 0 3"0 -0.0"1 1#6.#63 -606..01
0
200
400
600
0
200
400
600
-800
-700
-600
-500
-400
-300
-200
-100
12
;i4u0a '.) De:50mata 6i6temului 0utie0 multi6t0at
&i#ura '.' -opularea matricei de ri#iditate 3numrul de necunoscute < 2"222
Taelul '.1 Timpi de rezolvare 3secunde2
.imensiunea
prolemei
3mii de #dl2
5etoda
deplasrilor
rezolvare direct
?@&
5etoda
deplasrilor
rezolvare iterativ
pc#3?,&2
2 0.2" 0.:
10 ".' 10
"" 04 61
'6 1'16 168
60 118' 288
108 "06" 000
118 '8"6 680
14" ) :0"
"20 ) 288'
Aa rezolvrile directe s)folosit al#oritmul implicit din 5atla 3mldivide2.
1"
/n cazul rezolvrii iterative a sistemelor din metoda deplasrilor s)a folosit
al#oritmul #radientului conju#at precondi1ionat prin factorizare incomplet.
Toleran1a la factorizarea incomplet a fost 10
)".
.
5.3 ,tabilirea )arametrilor [Link] )rin metoda /radientului
con0u/at )recondi*ionat
Problema% 5odel ". de sistem rutier multistrat discretizat printr)o re1ea de 2" + 2" +
"" < 160'6 noduri i :480 elemente finite de tip tetraedru paraolic cu 10 noduri.
Numrul de necunoscute% nec<'2"61.
<0&<1i-l,i-l2= ;0&;1i-l2= 3idl% indicii deplasrilor liere2
&i#. '.4. -opularea matricei de ri#iditate
10
Numrul de termeni nenuli *n matricea de ri#iditate a structurii% 3371<02&3(01.#.
37& 30 te0me3i 3e3uli i3 mat0icea -e 0i4i-itate a
6t0uctu0ii
Nr tol
timp
;<cholinc3?r,tol2
3s2
nnz3;2
timp
itera1ii pc#
3s2
nr iteratii
precizie
otinuta
1 1e)2 0' 1"2'226 402 :::6 ".1e)10
2 1e)" 166 08680:" 0: "6: 1e)11
" 1e)0 ''0 1'1'00"0 :2 "02 1e)11
>serva1ie%
Toleran1a tolr<1e)12 folosit la rezolvarea pc# este prea mic. .e fiecare dat
solu1ia a rezultat prin sta#narea itera1iilor.
1ulare cu rezol-itor direct '11 secunde
(+emplu% )pozi1ia 2 din tael
B&A,G,;(A;(7,CT(;D < ".0000 0.0000 "6:.0000
&i#. '.6. -opularea matricei de ri#iditate incomplet factorizate cu toleran1a de 10
)"
1'
0 50 100 150 200 250 300 350 400
-10
-8
-6
-4
-2
0
2
4
iteratii
l
o
g
1
0
(
r
e
z
i
d
u
u
r
e
l
a
t
i
v
)
&i#. '.8. Conver#en1a itera1iilor de #radient conju#at precondi1ionat
Alte exemple
1ulare cu rezol-itor iterati- 2)c/3 4 rezol-itor direct
>> ex03-60m
t10 - eleme3te tet0ae-0u cu 10 35-u0i
+uma0ul -e 35-u0i : )"('
+uma0ul -e eleme3te : ""0'
+uma0ul -e ecuatii : 1013(
$impi -e calcul 16ecu3-e2
Calculul mat0icel50 eleme3tel50 : ".#0
A6am?la0ea mat0icei -e 0i4i-itate 6t0uctu0ala : 1."3
14
;act50i7a0ea C@5le6*8 i3c5mpleta : 1'.(#
Re75lva0ea 6i6temului -e ecuatii li3ia0e : 0.()
+uma0ul -e ite0atii : "3
E05a0ea 0elativa : "..e-011
1ulare cu rezol-itor direct
>> t100-
t10
Elap6e- time i6 ).0))03. 6ec53-6.
Elap6e- time i6 '.....1# 6ec53-6.
>> ex13-
t10 - eleme3te tet0ae-0u cu 10 35-u0i
+uma0ul -e 35-u0i : )3")'
+uma0ul -e eleme3te : ")('0
+uma0ul -e ecuatii : 10#3))
$impi -e calcul 16ecu3-e2
Calculul mat0icel50 eleme3tel50 : 30.6'
A6am?la0ea mat0icei -e 0i4i-itate 6t0uctu0ala : 1'.6)
;act50i7a0ea C@5le6*8 i3c5mpleta : 6#0."#
Re75lva0ea 6i6temului -e ecuatii li3ia0e : "3."0
+uma0ul -e ite0atii : ("
E05a0ea 0elativa : ..3e-011
16
1ulare cu rezol-itor direct
>> ex13-
t10
Elap6e- time i6 )'.(6.1'0 6ec53-6.
Elap6e- time i6 30(3.3630#" 6ec53-6.
>> ex"3-
t10 - eleme3te tet0ae-0u cu 10 35-u0i
+uma0ul -e 35-u0i : )(0#.
+uma0ul -e eleme3te : "(000
+uma0ul -e ecuatii : 11#030
$impi -e calcul 16ecu3-e2
Calculul mat0icel50 eleme3tel50 : 3).)'
A6am?la0ea mat0icei -e 0i4i-itate 6t0uctu0ala : 1#.3#
;act50i7a0ea C@5le6*8 i3c5mpleta : (').)#
Re75lva0ea 6i6temului -e ecuatii li3ia0e : "'.()
+uma0ul -e ite0atii : (3
E05a0ea 0elativa : (.3e-011
1ulare cu rezol-itor direct
>> ex"3-
t10
Elap6e- time i6 ').(3..6' 6ec53-6.
Elap6e- time i6 '#3(.3"0('0 6ec53-6.
18
>> <0&<1i-l,i-l2=
>> 3371<02& 10'3#((3
>> 3371%2 & ((.()"6
>> 3371S2 & 1"..'(10
>> A@56 <0 %
+ame Si7e !8te6 Cla66 Att0i?ute6
<0 11#030x11#030 16.'6)616 -5u?le 6pa06e
% ")x30x"(000 1'''"0000 -5u?le
>> 3371<02 & 10'3#((3
>> 3371R2 & 1"''6()0
>> 3371t0iu1<022 & '3"#3)1
>> 3371t0iu1<022B3371i-l2 & )'.1)
>> 3371R2B3371i-l2& 106.3#
t10 - eleme3te tet0ae-0u cu 10 35-u0i
+uma0ul -e 35-u0i : 1"('('
+uma0ul -e eleme3te : ('0"0
+uma0ul -e ecuatii : 3")3)#
$impi -e calcul 16ecu3-e2
Calculul mat0icel50 eleme3tel50 : 66.'6
A6am?la0ea mat0icei -e 0i4i-itate 6t0uctu0ala : '"1.#)
;act50i7a0ea C@5le6*8 i3c5mpleta : "(#3.11
Re75lva0ea 6i6temului -e ecuatii li3ia0e : 10".'0
+uma0ul -e ite0atii : 1'3
E05a0ea 0elativa : #e-011
3371<02 & ".#1(((0
3371R2 & 36#"').1
3uma0 me-iu -e te0me3i 3e3uli pe li3ie :
3371t0iu1<022B3371i-l2 & )6.)6'0
3371R2B3371i-l2 & 113.'3(0
5.5. !iblio/rafie+
B1D #$n%e&t' and a&&li%ati$n' $f finite ele(ent analy'i', &ourth (d.,
;...CooE, .a., FileG, 2002, p.1:)68.
1: