e - techmat.vgtu.lt
e - techmat.vgtu.lt
e - techmat.vgtu.lt
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
10. LAUKO UŽDAVINIAI<br />
– šilumos plitimo, elektrostatinių laukų, nesūkurinių srautų, difuzijos, …, uždaviniai. Visi šie<br />
uždaviniai aprašomi panašiomis, diferencialinėmis atžvilgiu lauko kintamojo Φ , lygtimis.<br />
BEMo formuluotes uždaviniams išvesim naudodami variacinius principus, 3D atvejui.<br />
Yra trys uždavinių klasės:<br />
10.1. Pusiausvyros uždaviniai<br />
Šių uždavinių tikslas – nustatyti pastovų, nepriklausantį nuo laiko, Φ pasiskirstymą<br />
erdvėje. Uždaviniai aprašomi kvaziharmonine lygtimi:<br />
∂ ⎛ ∂Φ ⎞ ∂ ⎛ ∂Φ ⎞ ∂ ⎛ ∂Φ ⎞<br />
⎜k<br />
x ⎟ + ⎜k<br />
y<br />
⎟ + ⎜k<br />
z ⎟ =<br />
∂x<br />
⎝ ∂x<br />
⎠ ∂y<br />
⎝ ∂y<br />
⎠ ∂z<br />
⎝ ∂z<br />
⎠<br />
f ( x, y,<br />
z),<br />
(10.1)<br />
čia: Φ – lauko kintamasis,<br />
k x , k y , k z , f – funkcijos, nepriklausančios nuo Φ ,<br />
k x , k y , k z ≠ 0 ir apibrėžti tūryje V,<br />
k x , k y , k z – neizotropinei terpei pateikti svarbiausiose ašyse x, y, z, o izotropinei terpei<br />
k x = k y = k z = k<br />
Kraštinės uždavinio sąlygos – dviejų tipų:<br />
ir<br />
čia:<br />
( x, y,<br />
z) ∈ S1<br />
Φ = Φ<br />
(10.2) – Dirichlet sąlygos<br />
k<br />
x<br />
∂Φ<br />
n<br />
∂x<br />
x<br />
+ k<br />
y<br />
∂Φ<br />
n<br />
∂y<br />
y<br />
+ k<br />
z<br />
∂Φ<br />
n<br />
∂z<br />
z<br />
+ g<br />
( x, y,<br />
z) + h( x,<br />
y,<br />
z) Φ = 0 ∈ S<br />
2<br />
n x , n y , n z - nukreiptų į terpės išorę normalių krypties kosinusai<br />
g, h – žinomos funkcijos<br />
(10.3) – Cauchy;<br />
jeigu g = h = 0 - Neumanno sąlygos<br />
(10.1, 2, 3) sudaro ”eliptinį kraštinių reikšmių uždavinį”. Uždavinys supaprastėja<br />
izotropinei terpei: (1) tampa<br />
2 f<br />
∇ Φ = ( x,<br />
y,<br />
z)<br />
– Poissono lygtis<br />
k<br />
arba, kai f = 0<br />
2<br />
∇ Φ = 0<br />
- Laplace’o lygtis<br />
78
Fizinė lygčių prasmė:<br />
Uždavinys Φ k x , k y , k z f g h<br />
Elektrostatinio<br />
lauko<br />
Elektros jėgų<br />
lauko<br />
intensyvumas<br />
Specifiniai<br />
talpumai<br />
Srovės<br />
ša<strong>lt</strong>iniai<br />
⎯<br />
⎯<br />
Nusistovėjęs<br />
šilumos<br />
plitimo<br />
Temperatūra<br />
Šiluminiai<br />
laidumai<br />
Šilumos<br />
ša<strong>lt</strong>iniai<br />
Šilumos nuotekis<br />
per sienas<br />
Konvekcinio<br />
šilumos<br />
nuotėkio<br />
koeficientas<br />
Nesūkurinio<br />
srauto<br />
Srauto<br />
funkcija<br />
⎯ 0 Kraštiniai<br />
greičiai<br />
0<br />
Fi<strong>lt</strong>racijos Slėgis Laidumai Srauto<br />
ša<strong>lt</strong>inis<br />
⎯<br />
⎯<br />
Tepimo (plonu<br />
skysčio<br />
sluoksniu)<br />
Slėgis<br />
kx, ky ∼ nuo<br />
sluoksnio<br />
storio ir<br />
klampio kz =0<br />
Srautas Nuotėkis ⎯<br />
Kūno sukimo<br />
Įtempimų<br />
funkcija<br />
Atvirkštiniai<br />
šlyties<br />
moduliai<br />
Susukimo<br />
kampas<br />
vienetiniam<br />
ilgiui<br />
⎯<br />
⎯<br />
. . .<br />
79
BEMo priklausomybių išvedimas variaciniu principu<br />
Lygtį (1) pagal Eulerio - Lagrange’o teoremą atitinka funkcionalas<br />
2<br />
2<br />
2<br />
1 ⎛<br />
⎞<br />
1 2<br />
( ) ⎜ ⎛ ∂Φ ⎞ ⎛ ∂Φ ⎞ ⎛ ∂Φ ⎞<br />
⎛<br />
⎞<br />
Ι Φ = k<br />
x<br />
k<br />
y<br />
k<br />
z<br />
2 f dV g h dS<br />
2<br />
2<br />
∫ ⎜ ⎟ + ⎜ ⎟ + ⎜ ⎟ + Φ⎟<br />
+<br />
x y z<br />
∫ ⎜ Φ + Φ ⎟ (10.4)<br />
⎜<br />
⎟<br />
2<br />
V ⎝ ⎝ ∂ ⎠ ⎝ ∂ ⎠ ⎝ ∂ ⎠ ⎠ S ⎝<br />
2<br />
⎠<br />
Ieškoma minimali funkcionalo reikšmė:<br />
δ<br />
Ι<br />
( Φ) = 0<br />
Visą sritį V sudalinam į elementus su n mazgų ir parenkam Φ kitimo dėsnį elemente<br />
Φ<br />
e<br />
=<br />
n<br />
∑<br />
i=<br />
1<br />
N Φ<br />
i<br />
i<br />
=<br />
e<br />
[ N ]{ Φ m },<br />
Φ , i = 1,2,...,<br />
i<br />
n<br />
n - laisvumo laipsnių elemente skaičius<br />
Φ - lauko kintamojo mazginės reikšmės.<br />
i<br />
N i turi tenkinti tik C 0 suderinamumo reikalavimus, nes (4) yra tik pirmosios Φ išvestinės.<br />
Funkcionalo minimumo sąlygos vienam baigtiniam elementui, vienam mazgui i<br />
reikalauja:<br />
e<br />
( Φ )<br />
i n<br />
∂Ι<br />
∂Φ<br />
i<br />
= 0 , = 1,2,...,<br />
e<br />
e<br />
e<br />
⎛ ∂Φ ∂ ⎛ ∂Φ ⎞ ∂Φ ⎞<br />
⎜<br />
⎟<br />
⎜<br />
⎟ +<br />
e<br />
= ∫ k<br />
x<br />
... + f<br />
dV + tik elementui ant pavirsiaus S<br />
2<br />
e⎝<br />
∂x<br />
∂Φ ⎝ ∂x<br />
⎠ ∂Φ<br />
V<br />
i<br />
i ⎠<br />
e<br />
e<br />
⎛ ∂Φ e ∂Φ ⎞<br />
e<br />
+ ∫ ⎜ g + hΦ<br />
⎟dS<br />
2 ⇒<br />
e ⎝ ∂Φ<br />
i<br />
∂Φ<br />
i ⎠<br />
S 2<br />
+<br />
⇒<br />
∫<br />
e<br />
V<br />
Turint omenyje, kad<br />
e<br />
∂Φ ⎡∂N<br />
∂N<br />
i j ⎤<br />
= ⎢ , ,... ⎥{ Φ<br />
m},<br />
∂x<br />
⎣ ∂x<br />
∂x<br />
⎦<br />
e<br />
∂ ⎛ ∂Φ ⎞ ∂N<br />
i<br />
⎜ = ,<br />
i<br />
x<br />
⎟<br />
∂Φ ⎝ ∂ ⎠ ∂x<br />
e<br />
∂Φ<br />
= N<br />
i;<br />
∂Φ<br />
⎛ ⎡∂N<br />
⎤<br />
⎜k<br />
x ⎢<br />
⎝ x ⎥<br />
⎣ ∂ ⎦<br />
i<br />
∂N<br />
i<br />
⎞ e<br />
e<br />
{ Φ<br />
e<br />
m } + ... + fN ( [ ] { e<br />
i<br />
⎟dV<br />
+ gN<br />
i<br />
h N N<br />
i m }) dS 2<br />
∫ + Φ<br />
∂x<br />
⎠<br />
e<br />
S 2<br />
80
Minimumo sąlygos vienam elementui tada<br />
e<br />
( Φ ) ⎫<br />
e<br />
e<br />
= 0 = [ K ]{ Φ } + [ K S ]{ Φ<br />
e<br />
m } + { R<br />
e<br />
2<br />
}<br />
⎧∂Ι<br />
⎨ e ⎬<br />
m<br />
⎩ ∂Φ ⎭<br />
n×n n×n n×1<br />
“Standumo” matricų elementai, atitinkantys laisvumo laipsnius i ir j:<br />
K<br />
e<br />
ij<br />
=<br />
∫<br />
e<br />
V<br />
⎛<br />
⎜k<br />
⎝<br />
∫<br />
x<br />
∂N<br />
i<br />
∂N<br />
∂x<br />
∂x<br />
j<br />
+ k<br />
y<br />
∂N<br />
i<br />
∂N<br />
∂y<br />
∂y<br />
j<br />
+ k<br />
z<br />
∂N<br />
i<br />
∂N<br />
∂z<br />
∂z<br />
j<br />
⎟ ⎞<br />
dV<br />
⎠<br />
e<br />
K e<br />
S ij = hN N dS<br />
tik elementams ant paviršiaus S 2<br />
R<br />
2 i j 2<br />
e<br />
S 2<br />
e<br />
e<br />
i = ∫ fN<br />
idV<br />
+ ∫<br />
e<br />
e<br />
V<br />
S 2<br />
gN dS<br />
i<br />
e<br />
2<br />
e<br />
Taigi, pagrindinė lygtis yra<br />
e e e e<br />
([ K ] + [ K S2 ]){ Φ m } + { R } = 0.<br />
10.2. Tikrinių reikšmių uždaviniai<br />
– mechaninių virpesių (modalinė analizė), mechanikos pradinio pastovumo, rezonanso<br />
akustikoje, stovinčių bangų sekliame vandenyje, elektromagnetinių bangų. Visi šie uždaviniai<br />
aprašomi Helmho<strong>lt</strong>zo lygtimi<br />
∂ ⎛<br />
⎜k<br />
∂x<br />
⎝<br />
x<br />
∂Φ ⎞<br />
⎟ +<br />
∂x<br />
⎠<br />
∂<br />
∂<br />
y<br />
⎛<br />
⎜k<br />
⎝<br />
y<br />
∂Φ ⎞<br />
⎟ +<br />
∂<br />
y ⎠<br />
∂<br />
∂<br />
z<br />
⎛<br />
⎜k<br />
⎝<br />
z<br />
∂Φ ⎞<br />
⎟<br />
∂<br />
z ⎠<br />
+ λ Φ = 0 ( 10.5 )<br />
su Dirichlet ir Neumanno kraštinės sąlygomis. Čia λ – nežinomas parametras.<br />
Eulerio - Lagrange’o funkcionalas Helmho<strong>lt</strong>zo lygčiai yra:<br />
2<br />
1 ⎛<br />
2<br />
2<br />
⎞<br />
2<br />
( ) ⎜ ⎛ ∂Φ ⎞ ⎛ ∂Φ ⎞ ⎛ ∂Φ ⎞<br />
Ι Φ =<br />
⎜ ⎟ + ⎜ ⎟ − Φ ⎟ = 0<br />
2<br />
∫ k<br />
x ⎜ ⎟ + k<br />
y<br />
k<br />
z<br />
λ dV<br />
(10.6)<br />
⎜<br />
⎟<br />
⎝ ⎝ ∂x<br />
⎠ ⎝ ∂y<br />
⎠ ⎝ ∂z<br />
V<br />
⎠ ⎠<br />
81
Panašiai kaip skyriuje 10.1 variaciniu principu gautume:<br />
e<br />
[ ]{ Φ } [ e<br />
m − M<br />
e ]{ Φ<br />
e<br />
m } = 0,<br />
K λ (10.7)<br />
k<br />
m<br />
e<br />
ij<br />
e<br />
ij<br />
=<br />
=<br />
∫<br />
e<br />
V<br />
∫e<br />
V<br />
⎛<br />
⎜k<br />
⎝<br />
x<br />
N N<br />
i<br />
∂N<br />
i<br />
∂N<br />
∂x<br />
∂x<br />
j<br />
dV<br />
e<br />
;<br />
j<br />
+ k<br />
y<br />
∂N<br />
i<br />
∂y<br />
∂N<br />
j<br />
∂y<br />
+ k<br />
z<br />
∂N<br />
i<br />
∂z<br />
∂N<br />
∂z<br />
j<br />
⎟ ⎞<br />
dV<br />
⎠<br />
e<br />
,<br />
N i šiems uždaviniams - taip pat C 0 tipo. (7) lygtis vadinama savųjų (arba tikrinių) reikšmių λ<br />
uždaviniu. Netrivialus sprendinys yra tik tada, kai<br />
[ K] − λ [ M ] = 0.<br />
Determinantą galima išskleisti į n eilės polinomą, turintį n šaknų λ<br />
i<br />
; n – sistemos<br />
laisvumo laipsnių kiekis. Kiekvieną λ<br />
i<br />
atitinka begalinis skaičius savųjų vektorių { Φ m<br />
} i<br />
,<br />
tenkinančių lygtį (7). Dažniausiai savieji vektoriai pateikiami normalizuota forma.<br />
10.3. Plitimo uždaviniai<br />
Tai – lauko kintamojo pasiskirstymo erdvėje priklausomai nuo laiko uždavinys:<br />
difuzijos, bangų sklidimo uždaviniai.<br />
Lygtis:<br />
ir<br />
∂ ⎛<br />
⎜k<br />
∂x<br />
⎝<br />
∂Φ ⎞ ∂ ⎛<br />
⎟ + ⎜k<br />
∂x<br />
⎠ ∂y<br />
⎝<br />
∂Φ ⎞ ∂ ⎛<br />
⎟ + ⎜k<br />
∂y<br />
⎠ ∂z<br />
⎝<br />
∂Φ ⎞<br />
⎟ = f<br />
∂z<br />
⎠<br />
x y<br />
z<br />
,<br />
Kraštinės sąlygos:<br />
( x y,<br />
z,<br />
t) ∈ S ,<br />
,<br />
1<br />
Φ = Φ<br />
t > 0<br />
( x y,<br />
z t) + c Φ & + m Φ &<br />
, (10.8)<br />
k<br />
x<br />
∂Φ<br />
n<br />
∂x<br />
x<br />
+ k<br />
y<br />
∂Φ<br />
n<br />
∂y<br />
y<br />
+ k<br />
z<br />
∂Φ<br />
n<br />
∂z<br />
z<br />
+ g<br />
( x y,<br />
z,<br />
t) + h( x,<br />
y,<br />
z,<br />
t) Φ = 0 ∈ S ,<br />
,<br />
2<br />
t > 0 (10.9)<br />
Pradinės sąlygos:<br />
( x,<br />
y,<br />
z) ∈ ,<br />
( x,<br />
y,<br />
z) ∈ ,<br />
Φ = Φ<br />
0<br />
V t > 0 (10.10)<br />
Φ & = Φ&<br />
t > 0<br />
0<br />
V<br />
82
(8) lygtis vadinama paraboline, jei m = 0 ir c ≠ 0;<br />
hiperboline, jei m ≠ 0.<br />
Parabolinėms lygtims būdingas lauko kintamojo plitimas terpėje begaliniu greičiu<br />
(pvz., paveikus kūną šiluma, jo vidinių taškų temperatūra pakistų akimirksniu), o<br />
hiperbolinėms - baigtiniu greičiu.<br />
BEM formuluotės<br />
Naudojami keli būdai. Dažniausiai uždavinys nagrinėjamas fiksuotu laiko momentu ir<br />
laikoma, kad Φ & ir Φ&<br />
∼ x, y, z. Tada:<br />
Φ<br />
e<br />
n<br />
( x,<br />
y,<br />
z,<br />
t) = ∑ N<br />
i<br />
( x,<br />
y,<br />
z) Φ<br />
i<br />
( t)<br />
i=<br />
1<br />
- tik mazginiai laisvumo laipsniai priklauso nuo laiko. Dabar elemento lygtis šiam fiksuotam<br />
laiko momentui galima išvesti variaciniu principu (jei uždavinys turi funkcionalą) arba<br />
svertinių neatitikčių metodais.<br />
Kitas metodas (J. T. Odenas, 1969): uždavinys nagrinėjamas keturmatėje erdvėje:<br />
Φ<br />
e<br />
n<br />
( x,<br />
y,<br />
z,<br />
t) = ∑ N<br />
i<br />
( x,<br />
y,<br />
z,<br />
t) Φ<br />
i<br />
i=<br />
1<br />
- pridedama dar viena dimensija, o tai labai pabrangina sprendimą.<br />
Pirmuoju metodu galima gauti:<br />
e<br />
[ M ]{ Φ&<br />
e e<br />
} + [ C ]{ Φ&<br />
e<br />
e<br />
m<br />
} + [ ]{ Φ } + [ e<br />
m K<br />
e<br />
m K<br />
e<br />
S 2 ]{ Φ<br />
e<br />
m } + { R ( t)<br />
} = 0<br />
& (10.11)<br />
čia<br />
m<br />
c<br />
k<br />
k<br />
ij<br />
R<br />
ij<br />
ij<br />
=<br />
=<br />
=<br />
e<br />
S2<br />
ij<br />
e<br />
i<br />
∫e<br />
V<br />
∫e<br />
V<br />
e<br />
V<br />
=<br />
mN N<br />
i<br />
i<br />
cN N<br />
j<br />
j<br />
dV<br />
dV<br />
⎛ ∂N<br />
∂N<br />
i j<br />
⎜k<br />
x<br />
⎝ ∂x<br />
∂x<br />
= hN N dS<br />
∫<br />
i j ∫e<br />
S2<br />
e<br />
∫ f N<br />
idV<br />
+<br />
e<br />
V<br />
e<br />
e<br />
+ k<br />
e<br />
2<br />
∫<br />
e<br />
S 2<br />
i<br />
y<br />
g dS<br />
∂N<br />
i<br />
∂y<br />
e<br />
2<br />
∂N<br />
j<br />
∂y<br />
+ k<br />
z<br />
∂N<br />
∂z<br />
i<br />
∂N<br />
∂z<br />
j<br />
⎟ ⎞<br />
dV<br />
⎠<br />
e<br />
Lygtyje (4) [M] dažniausiai įvertina inerciją, [C] - slopinimą arba talpumą, [K] - standumą.<br />
Kai kuriuose uždaviniuose [M]=0 (difuzijos), [C]=0 (neslopinamų svyravimų). Priklausomai<br />
nuo R(t) kitimo skiriami periodiniai, harmoniniai, atsitiktiniai uždaviniai.<br />
Sprendimo metodai (11) lygčiai - dažniausiai rekurentiniai.<br />
83
11. BEM PAKLAIDOS. TESTAI. PROGRAMOS<br />
BEM sprendinio paklaidos (esant korektiškiems elementams, tinkamai naudojant<br />
programą):<br />
1. modeliavimo – atspindi matematinio modelio skirtumą nuo realybės (pvz.,<br />
Kirchhofo b.e. storai plokštelei modeliuoti);<br />
2. diskretizavimo – kyla dėl perėjimo nuo begalinio prie baigtinio laisvumo laipsnių<br />
kiekio modelyje;<br />
3. apvalinimo – kyla dėl baigtinio kompiuterio žodžio ilgio.<br />
Pastarosios paklaidos itin didelę įtaka turi “blogai sąlygotiems” BEM uždaviniams.<br />
Matematinis blogo sąlygotumo pavyzdys:<br />
⎡ 1.00<br />
⎢<br />
⎣−1.00<br />
−1.00⎤⎧x⎫<br />
⎧ 4.00 ⎫ ⎧x⎫<br />
⎧104⎫<br />
⎬ → ⎨ ⎬ = ⎨ ⎬;<br />
1.02<br />
⎥⎨<br />
⎬ = ⎨<br />
⎦⎩y⎭<br />
⎩−<br />
2.00⎭<br />
⎩y⎭<br />
⎩100⎭<br />
pakeitus vieną koeficientą 1%, rezu<strong>lt</strong>atas pakinta ∼ 100%:<br />
⎡ 1.00<br />
⎢<br />
⎣−1.00<br />
−1.00⎤⎧x⎫<br />
⎧ 4.00 ⎫ ⎧x⎫<br />
⎧204⎫<br />
⎬ → ⎨ ⎬ = ⎨ ⎬.<br />
1.01<br />
⎥⎨<br />
⎬ = ⎨<br />
⎦⎩y⎭<br />
⎩−<br />
2.00⎭<br />
⎩y⎭<br />
⎩200⎭<br />
Toks uždavinys 2- jų skaitmenų tikslumų neišsprendžiamas.<br />
Blogai sąlygoto mechanikos uždavinio pavyzdys:<br />
x, u k 2<br />
k 1<br />
p<br />
Jei k 1 >>k 2 :<br />
⎡ k1<br />
⎢<br />
⎣−<br />
k<br />
1<br />
− k1<br />
⎤⎧u1<br />
⎫ ⎧P⎫<br />
⎨ ⎬ = ⎨ ⎬;<br />
k1<br />
+ k<br />
⎥<br />
2 ⎦⎩u<br />
2 ⎭ ⎩0⎭<br />
k<br />
1<br />
+ k<br />
2<br />
→ k ,<br />
1<br />
ir eilutės tampa tiesiškai priklausomos. Tokiu atveju reikia itin didelio skaičiavimo tikslumo:<br />
pvz., jei k 1 =40., k 2 =0.0014, tai, kad gautume sprendinį, turime žinoti k 1 tokiu tikslumu:<br />
k 1 =40.0000. Antraip k 1 paklaida lemia viso skaičiavimo tikslumą.<br />
Šiame sprendinyje modeliavimo paklaida ∼ 0, diskretizavimo = 0.<br />
84
Kitas pavyzdys: rėmas, kuriam ašinis standumas >> lenkimo standumas<br />
turėtų gautis u1 ≅ u 2<br />
, o tai rodo, kad dvi lygtis yra beveik tiesiškai priklausomos.<br />
P 1 2<br />
Kitas pavyzdys: plona plokštelė modeliuota Mindlino b.e.<br />
Blogą sąlygotumą gali ke<strong>lt</strong>i ir netinkamas diskretizavimas:<br />
- iškraipytos formos baigtinis elementas<br />
- elementas su dideliu kraštinių santykiu<br />
- apjungti skirtingų dydžių baigtiniai elementai<br />
Vienintelis “vaistas” nuo blogo sąlygotumo – ilgesnis kompiuterio žodis.<br />
Paklaidų įverčiai<br />
Blogo sąlygotumo koeficientas (condition number), C. Moler, 1967:<br />
C<br />
λ<br />
max<br />
( K ) = , λ − [ K ]<br />
savosios reikšmės.<br />
λmin<br />
≈ lgC ( K )<br />
Pavyzdys: kompiuterio žodžio ilgis – 7 ženklai, C(K)=1E5 – sprendinys turės tik 2<br />
patikimus skaitmenis.<br />
Trūkumai: dažnai pernelyg pesimistinis įvertis; įvertina tik apvalinimo paklaidas,<br />
padarytas iki lygčių sprendimo; neįvertina ir dešinės lygties pusės paklaidų.<br />
Įvertis tinklelio kokybei – diskretizavimo paklaidoms, I. Fried, 1972:<br />
čia<br />
C<br />
2m−1<br />
2m<br />
⎛ h ⎞<br />
= ⎜ h ⎟ ,<br />
⎝ min ⎠<br />
max<br />
n<br />
( K ) b⎜<br />
⎟ N<br />
b1<br />
b = ,<br />
1 − 2υ<br />
b 1 ∼ E, υ ir ∼ nuo uždavinio tipo<br />
h – elementų “diametrai”<br />
N – elementų kiekis<br />
2m – diferencialinių lygčių eilė<br />
n – uždavinio matas<br />
85
Įvertis lygčių sprendimo metu kylančioms paklaidoms, R. McNeal, 1976:<br />
e =<br />
{ U} T<br />
{ ∆P}<br />
{ U} { P} , T<br />
{ ∆ P} = { P} − [ K ]{ U}<br />
- santykis likutinių jėgų darbo su tikrųjų jėgų darbu.<br />
Pavyzdys: jei e ≤ 1E<br />
− 7 – viengubas kompiuterio žodis lygčių sprendimo metu<br />
paklaidų neįneša.<br />
Praktinė įverčių nauda: A. Anon, 1986: sprendinys, matyt, bus priimtinas, jei vienas iš<br />
įverčių rodo blogą padėtį; sprendinys tikrai nepriimtinas, jei visi trys įverčiai yra blogų<br />
reikšmių.<br />
tūriui.<br />
čia<br />
Kiti įverčiai diskretizavimo paklaidoms.<br />
Paprasčiausias įvertis: diskretizuotos struktūros tūris turi būti lygus realios struktūros<br />
Santykinės klaidos įvertis (jo matas - %), R. Cook, 1985:<br />
e r<br />
= ρ ρ h<br />
h =<br />
0<br />
1<br />
1<br />
1<br />
n<br />
q−r<br />
2 0<br />
,<br />
,<br />
N<br />
q = 1 + <br />
r = 0 poslinkių santykinei klaidai<br />
r = 1 deformacijų ir įtempimų klaidai, C 0 tipo uždaviniams<br />
r = 2 deformacijų ir įtempimų klaidai, C 1 tipo uždaviniams<br />
ρ<br />
1<br />
– maksimalus baigtinių elementų matmenų tinklelyje santykis<br />
ρ<br />
2<br />
– maksimalus didžiausio elemento tinklelyje ir mažiausio “diametrų” santykis<br />
N – elementų kiekis<br />
n – matas<br />
Praktinė e r nauda: jei gautas e r ≥ 1 (t.y., 100%) – reikalingas tinklelio smulkinimas; jei<br />
e r ≈ 0.1x – rezu<strong>lt</strong>atai, matyt, yra patikimi. Kartais įvertis yra labai<br />
pesimistiškas.<br />
86
Prisitaikančios (self-adaptive) BEM programos<br />
Galimybę joms rastis suteikė patikimų įverčių nustatymas. N.Kikuchi, 1986: programa<br />
nusprendžia, ar reikia smulkinti ( jei reikia – kur ir iki kokio lygio) elementų tinklelį;<br />
sprendimo etapas; rezu<strong>lt</strong>atų analizė, kol pasiekiama reikiama sprendinio tolerancija arba<br />
leidžia kompiuterio resursai, arba kol neįsivyrauja apvalinimo paklaidų triukšmas.<br />
Tinklelio smulkinimas: h- ir p- strategijos. h strategija: tinklelio smulkinimas;<br />
pageidautina – išsaugant visus ankstesnius laisvumo laipsnius. p -strategija: reikiamuose<br />
mazguose įvedami papildomi laisvumo laipsniai – naudojamas aukštesnės eilės interpoliacinis<br />
polinomas. Tokioje, hierarchinėje sprendimo schemoje nebūtina iš naujo perskaičiuoti [K],<br />
[M], [C]: matricos tik papildomos naujomis eilutėmis ir stulpeliais.<br />
Testai elemento kokybei įvertinti<br />
Tikrinių reikšmių testas.<br />
Idėja: [ K ]{ u} = { P}<br />
; mažų poslinkių atveju, tampriai medžiagai<br />
P = λ u<br />
{ } { };<br />
⎛ ⎞<br />
⎜ K λ ⎟<br />
(11.1)<br />
⎝ ⎠<br />
tada [ ] − [ I ] { u} = 0<br />
T<br />
Tikrinius vektorius normalizavus: { u } i { u} = 1,<br />
o (1) padauginus iš normalizuoto { u } i :<br />
{ u} T<br />
i [ K]{ u} = λ<br />
i<br />
= 2U<br />
i<br />
i<br />
i<br />
– deformacijos energija, esant mazginiams {u} i . Iš čia: kai {u} i atitinka baigtinio elemento<br />
poslinkius kaip kieto nedeformuojamo kūnu poslinkius λ i<br />
turi būti = 0.<br />
Todėl 2D uždaviniui turim gauti 3 nulines tikrines reikšmes, 3D arba kevalui – 6,<br />
simetriniam ašiai uždaviniui – 1. Jei testas rodo per mažai nulinių reikšmių – elementas<br />
netenkina pilnumo reikalavimų; jei per daug – elemente yra mechanizmų.<br />
Testą galima pakartoti elementą perorientavus erdvėje. Jei dabar gaunamos kitos λ –<br />
elementas nėra geometriškai invariantiškas.<br />
“Vieno elemento testas”: tiriamas elemento atsakas į apkrovą, keičiant elemento<br />
parametrus: ilgio/aukščio santykį, kampą ir pan. Teikia informaciją apie leistinas šių dydžių<br />
reikšmes.<br />
“Patch testas”.<br />
Standartiniai testai sulyginimui su analitiniais sprendiniais.<br />
T<br />
Reziume: elemento testavimui nepakanka vieno testo. Be to, testai gali tik parodyti<br />
klaidas, tačiau negali įrodyti, kad klaidų nėra.<br />
87
BEM programos<br />
Siūloma ∼ 100 komercinių paketų. Visos ankstesnės programos – Fortrano kalba, vėlesnės – ir<br />
C++. Programų apimtis: ∼ 100 000 ÷ 1 000 000 operatorių. Galimybės: tiesinė ir netiesinė<br />
analizė, stacionarūs ir nestacionarūs procesai, mechanikos, šilumos, elektromagnetizmo<br />
taikymai, mu<strong>lt</strong>ifizikos uždaviniai. Modifikuotos versijos pasirodo kas 1-2 metai, programos<br />
pritaikytos įvairioms platformoms; yra ir suprastintos versijos PC.<br />
Dokumentacijos apimtis – tūkstančiai puslapių; kursai vartotojams, konsu<strong>lt</strong>acijos<br />
telefonu ir pan.<br />
Per pastaruosius keliolika metų rinkoje nepasirodė nė viena nauja BEM programa: itin<br />
aukšti “starto” kaštai, rinkos inercija.<br />
Kūrimo kaštai (J. Rice, 1983): analizuota 400 su BEM susijusių projektų, kurių trukmė<br />
1-8 metai, 2-200 darbuotojų. Įskaitytas laikas algoritmo, programos kūrimui, testavimui,<br />
dokumentavimui. Produktyvumas matuotas programos operatorių kiekiu<br />
1darbuotojui/1mėnesiui.<br />
Rezu<strong>lt</strong>atai: 5-5000 eilučių, vidurkis 200 eilučių; 2/3 projektų pakliuvo į ribas 75-550<br />
eilučių.<br />
Kaina didelių programų projektų - 2-10 mln. USD. Palaikymo kaštai: kaip taisyklė<br />
viršija kūrimo kaštus.<br />
88