MES w statyce ośrodka ciągłego
Piotr Pluciński
e-mail: Piotr.Plucinski@pk.edu.pl
Jerzy Pamin
e-mail: Jerzy.Pamin@pk.edu.pl
Katedra Technologii Informatycznych w Inżynierii Wydział Inżynierii Lądowej Politechniki Krakowskiej
Strona domowa: www.CCE.pk.edu.pl
Metody obliczeniowe, 2020 P.Plucińskic
Zakres prezentacji
1 Stan równowagi
2 Dyskretyzacja MES
3 Płaski stan naprężenia
4 Przykład
Stan równowagi
Y Z
X
S
n t
u(x, y, z) ρb(x, y, z)
P
P
0V
Wektor gęstości sił masowych [N/m
3] ρb = ρ
" 0 0
−g
#
Wektor gęstości sił powierzchniowych [N/m
2]
t =
" t
x
t
yt
z#
Przemieszczenie, odkształcenie, naprężenie (notacja Voigta)
u =
u
xu
yu
z, =
xx
xy
xz
xy
yy
yz
xz
yz
zz
→
x
y
z
γ
xyγ
yzγ
zx
, σ =
σ
xxσ
xyσ
xzσ
xyσ
yyσ
yzσ
xzσ
yzσ
zz→
σ
xσ
yσ
zτ
xyτ
yzτ
zx
Metody obliczeniowe, 2020 P.Plucińskic
Stan równowagi
Y Z
X
S
n t
u(x, y, z) ρb(x, y, z)
P
P
0V
Wektor gęstości sił masowych [N/m
3] ρb = ρ
" 0 0
−g
#
Wektor gęstości sił powierzchniowych [N/m
2]
t =
" t
x
t
yt
z#
Przemieszczenie, odkształcenie, naprężenie (notacja Voigta)
u =
u
xu
yu
z, =
xx
xy
xz
xy
yy
yz
xz
yz
zz
→
x
y
z
γ
xyγ
yzγ
zx
, σ =
σ
xxσ
xyσ
xzσ
xyσ
yyσ
yzσ
xzσ
yzσ
zz→
σ
xσ
yσ
zτ
xyτ
yzτ
zx
Stan równowagi
Y Z
X
S
n t
u(x, y, z) ρb(x, y, z)
P
P
0V
Wektor gęstości sił masowych [N/m
3] ρb = ρ
" 0 0
−g
#
Wektor gęstości sił powierzchniowych [N/m
2]
t =
" t
x
t
yt
z#
Przemieszczenie, odkształcenie, naprężenie (notacja Voigta)
u =
u
xu
yu
z, =
xx
xy
xz
xy
yy
yz
xz
yz
zz
→
x
y
z
γ
xyγ
yzγ
zx
, σ =
σ
xxσ
xyσ
xzσ
xyσ
yyσ
yzσ
xzσ
yzσ
zz→
σ
xσ
yσ
zτ
xyτ
yzτ
zx
Metody obliczeniowe, 2020 P.Plucińskic
Stan równowagi
Y Z
X
S
n t
u(x, y, z) ρb(x, y, z)
P
P
0V
Wektor gęstości sił masowych [N/m
3] ρb = ρ
" 0 0
−g
#
Wektor gęstości sił powierzchniowych [N/m
2]
t =
" t
x
t
yt
z#
Przemieszczenie, odkształcenie, naprężenie (notacja Voigta)
x
σ
x
σ
Stan równowagi
Y Z
X
S
n t
u(x, y, z) ρb(x, y, z)
P
P
0V
Równanie równowagi ciała
Z
S
tdS + Z
V
ρbdV = 0
Statyczne warunki brzegowe t = σn gdzie σ – tensor naprężenia
Wykorzystując twierdzenie Greena–Gaussa–Ostrogradzkiego
Z
S
σndS = Z
V
L T σdV gdzie L – macierz operatorów różniczkowych
Metody obliczeniowe, 2020 P.Plucińskic
Stan równowagi
Y Z
X
S
n t
u(x, y, z) ρb(x, y, z)
P
P
0V
Równanie równowagi ciała
Z
S
tdS + Z
V
ρbdV = 0
Statyczne warunki brzegowe t = σn gdzie σ – tensor naprężenia
Wykorzystując twierdzenie Greena–Gaussa–Ostrogradzkiego
Z
S
σndS = Z
V
L T σdV gdzie L – macierz operatorów różniczkowych
Stan równowagi
Y Z
X
S
n t
u(x, y, z) ρb(x, y, z)
P
P
0V
Równanie równowagi ciała
Z
S
tdS + Z
V
ρbdV = 0
Statyczne warunki brzegowe t = σn gdzie σ – tensor naprężenia
Wykorzystując twierdzenie Greena–Gaussa–Ostrogradzkiego
Z
S
σndS = Z
V
L T σdV gdzie L – macierz operatorów różniczkowych
Metody obliczeniowe, 2020 P.Plucińskic
Równania równowagi
Równania Naviera Z
V
L T σ + ρb dV = 0 ⇐⇒L T σ + ρb = 0 ∀P ∈ V
σ ij,j + ρb i = 0
Sformułowanie słabe – funkcja wagowa w ≡ δu – kinematycznie dopuszczalna wariacja przemieszczenia
Z
V
(δu) T L T σ + ρb dV = 0 ∀δu
Równania równowagi
Równania Naviera Z
V
L T σ + ρb dV = 0 ⇐⇒L T σ + ρb = 0 ∀P ∈ V σ ij,j + ρb i = 0
Sformułowanie słabe – funkcja wagowa w ≡ δu – kinematycznie dopuszczalna wariacja przemieszczenia
Z
V
(δu) T L T σ + ρb dV = 0 ∀δu
Metody obliczeniowe, 2020 P.Plucińskic
Równania równowagi
Równania Naviera Z
V
L T σ + ρb dV = 0 ⇐⇒L T σ + ρb = 0 ∀P ∈ V σ ij,j + ρb i = 0
Sformułowanie słabe – funkcja wagowa w ≡ δu – kinematycznie dopuszczalna wariacja przemieszczenia
Z
V
(δu) T L T σ + ρb dV = 0 ∀δu
Równania równowagi
Równania Naviera Z
V
L T σ + ρb dV = 0 ⇐⇒L T σ + ρb = 0 ∀P ∈ V σ ij,j + ρb i = 0
Sformułowanie słabe – funkcja wagowa w ≡ δu – kinematycznie dopuszczalna wariacja przemieszczenia (zgodna z kinematycznymi warunkami brzegowymi)
Z
V
(δu) T L T σ + ρb dV = 0 ∀δu
− Z
V
(Lδu) T σdV + Z
S
(δu) T σndS + Z
V
(δu) T ρbdV = 0
Metody obliczeniowe, 2020 P.Plucińskic
Równania równowagi
Równania Naviera Z
V
L T σ + ρb dV = 0 ⇐⇒L T σ + ρb = 0 ∀P ∈ V σ ij,j + ρb i = 0
Sformułowanie słabe – funkcja wagowa w ≡ δu – kinematycznie dopuszczalna wariacja przemieszczenia (zgodna z kinematycznymi warunkami brzegowymi)
Z
V
(δu) T L T σ + ρb dV = 0 ∀δu
− Z
V
(Lδu) T σdV + Z
S
(δu) T σn t
dS + Z
V
(δu) T ρbdV = 0
Równania równowagi
Równania Naviera Z
V
L T σ + ρb dV = 0 ⇐⇒L T σ + ρb = 0 ∀P ∈ V σ ij,j + ρb i = 0
Sformułowanie słabe – funkcja wagowa w ≡ δu – kinematycznie dopuszczalna wariacja przemieszczenia (zgodna z kinematycznymi warunkami brzegowymi)
Z
V
(δu) T L T σ + ρb dV = 0 ∀δu
− Z
V
(Lδu) T σdV + Z
S
(δu) T tdS + Z
V
(δu) T ρbdV = 0
Metody obliczeniowe, 2020 P.Plucińskic
Równania równowagi
Równania Naviera Z
V
L T σ + ρb dV = 0 ⇐⇒L T σ + ρb = 0 ∀P ∈ V σ ij,j + ρb i = 0
Sformułowanie słabe – funkcja wagowa w ≡ δu – kinematycznie dopuszczalna wariacja przemieszczenia – zasada prac wirtualnych
Z
V
(δu) T L T σ + ρb dV = 0 ∀δu
Z
V
(Lδu) T σdV = Z
S
(δu) T tdS + Z
V
(δu) T ρbdV
Równania równowagi
Równania Naviera Z
V
L T σ + ρb dV = 0 ⇐⇒L T σ + ρb = 0 ∀P ∈ V σ ij,j + ρb i = 0
Sformułowanie słabe – funkcja wagowa w ≡ δu – kinematycznie dopuszczalna wariacja przemieszczenia – zasada prac wirtualnych
Z
V
(δu) T L T σ + ρb dV = 0 ∀δu
Z
V
(Lδu) T σdV = Z
S
(δu) T tdS + Z
V
(δu) T ρbdV
praca sił wewnętrznych praca sił zewnętrznych
Metody obliczeniowe, 2020 P.Plucińskic
Dyskretyzacja MES (n=LWE, N =LSSU, E=LEU)
Aproksymacja pola przemieszczeń u eh =
n
X
i=1
N i e (ξ, η, ζ)d e i = N e d e
N e
[3×3n] =
N 1 e 0 0 . . . N n e 0 0 0 N 1 e 0 . . . 0 N n e 0 0 0 N 1 e . . . 0 0 N n e
d e
[3n×1] =
d e 1 . . . d e n
y z
2 8
9
11 14
17
η ζ
η ζ
i
j l
m n
o p d e
[3n×1]
= T I e
[3n×N ]
d
[N ×1]
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi (ρb e = f e – wektor sił objętościowych)
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
Metody obliczeniowe, 2020 P.Plucińskic
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
E
X
e=1
Z
Ve
(L
eN
eδd
e)
Tσ
edV
e−
Z
Se
(N
eδd
e)
Tt
edS
e−
Z
Ve
(N
eδd
e)
Tf
edV
e= 0
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
E
X
e=1
Z
Ve
(L
eN
eB
eδd
e)
Tσ
edV
e−
Z
Se
(N
eδd
e)
Tt
edS
e−
Z
Ve
(N
eδd
e)
Tf
edV
e= 0
Metody obliczeniowe, 2020 P.Plucińskic
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
E
X
e=1
Z
Ve
(B
eδd
e)
Tσ
edV
e−
Z
Se
(N
eδd
e)
Tt
edS
e−
Z
Ve
(N
eδd
e)
Tf
edV
e= 0
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
E
X
e=1
Z
Ve
(B
eδd
e)
Tσ
edV
e−
Z
Se
(N
eδd
e)
Tt
edS
e−
Z
Ve
(N
eδd
e)
Tf
edV
e= 0
E
X
e=1
(δd
e)
TZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
Metody obliczeniowe, 2020 P.Plucińskic
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
E
X
e=1
Z
Ve
(B
eδd
e)
Tσ
edV
e−
Z
Se
(N
eδd
e)
Tt
edS
e−
Z
Ve
(N
eδd
e)
Tf
edV
e= 0
E
X
e=1
( δd
eI T
eδd
)
TZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
E
X
e=1
Z
Ve
(B
eδd
e)
Tσ
edV
e−
Z
Se
(N
eδd
e)
Tt
edS
e−
Z
Ve
(N
eδd
e)
Tf
edV
e= 0
E
X
e=1
( I T
eδd)
TZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
Metody obliczeniowe, 2020 P.Plucińskic
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
E
X
e=1
Z
Ve
(B
eδd
e)
Tσ
edV
e−
Z
Se
(N
eδd
e)
Tt
edS
e−
Z
Ve
(N
eδd
e)
Tf
edV
e= 0
E
X
e=1
( I T
eδd)
TZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
(δd)
TE
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
E
X
e=1
Z
Ve
(B
eδd
e)
Tσ
edV
e−
Z
Se
(N
eδd
e)
Tt
edS
e−
Z
Ve
(N
eδd
e)
Tf
edV
e= 0
E
X
e=1
( I T
eδd)
TZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
(δd)
T∀δd
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
Metody obliczeniowe, 2020 P.Plucińskic
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
E
X
e=1
Z
Ve
(B
eδd
e)
Tσ
edV
e−
Z
Se
(N
eδd
e)
Tt
edS
e−
Z
Ve
(N
eδd
e)
Tf
edV
e= 0
E
X
e=1
( I T
eδd)
TZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
E
X
e=1
Z
Ve
(B
eδd
e)
Tσ
edV
e−
Z
Se
(N
eδd
e)
Tt
edS
e−
Z
Ve
(N
eδd
e)
Tf
edV
e= 0
E
X
e=1
( I T
eδd)
TZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e=
E
X
e=1
I T
eTZ
Se
N
eTt
edS
e+
Z
Ve
N
eTf
edV
eMetody obliczeniowe, 2020 P.Plucińskic
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
Z
V
e(L e δu e ) T σ e dV e − Z
S
e(δu e ) T t e dS e − Z
V
e(δu e ) T f e dV e
= 0
E
X
e=1
Z
Ve
(B
eδd
e)
Tσ
edV
e−
Z
Se
(N
eδd
e)
Tt
edS
e−
Z
Ve
(N
eδd
e)
Tf
edV
e= 0
E
X
e=1
( I T
eδd)
TZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e− Z
Se
N
eTt
edS
e− Z
Ve
N
eTf
edV
e= 0
E
X Z
EX Z Z
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e=
E
X
e=1
I T
eTZ
Se
N
eTt
edS
e+
Z
Ve
N
eTf
edV
eUwzględnienie związków kinematycznych i konstytutywnych liniowa sprężystość: σ = Dε
liniowy związek kinematyczny: ε = Lu
σ e = D e L e u e = D e L e N e d e = D e B e T I e d
Równanie równowagi
Metody obliczeniowe, 2020 P.Plucińskic
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e=
E
X
e=1
I T
eTZ
Se
N
eTt
edS
e+
Z
Ve
N
eTf
edV
eUwzględnienie związków kinematycznych i konstytutywnych liniowa sprężystość: σ = Dε
liniowy związek kinematyczny: ε = Lu
σ e = D e L e u e = D e L e N e d e = D e B e T I e d
Równanie równowagi
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e=
E
X
e=1
I T
eTZ
Se
N
eTt
edS
e+
Z
Ve
N
eTf
edV
eUwzględnienie związków kinematycznych i konstytutywnych liniowa sprężystość: σ = Dε
liniowy związek kinematyczny: ε = Lu
σ e = D e L e u e = D e L e N e d e = D e B e T I e d
Równanie równowagi
E
X
e=1
I T
eTZ
Ve
B
eTD
eB
eT I
ed dV
e=
E
X
e=1
I T
eTZ
Se
N
eTt
edS
e+
Z
Ve
N
eTf
edV
eMetody obliczeniowe, 2020 P.Plucińskic
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e=
E
X
e=1
I T
eTZ
Se
N
eTt
edS
e+
Z
Ve
N
eTf
edV
eUwzględnienie związków kinematycznych i konstytutywnych liniowa sprężystość: σ = Dε
liniowy związek kinematyczny: ε = Lu
σ e = D e L e u e = D e L e N e d e = D e B e T I e d
Równanie równowagi
E
X T I
eT( Z
B
eTD
eB
edV
e)
I T
ed =
E
X T I
eT( Z
N
eTt
edS
e+
Z
N
eTf
edV
e)
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e=
E
X
e=1
I T
eTZ
Se
N
eTt
edS
e+
Z
Ve
N
eTf
edV
eUwzględnienie związków kinematycznych i konstytutywnych liniowa sprężystość: σ = Dε
liniowy związek kinematyczny: ε = Lu
σ e = D e L e u e = D e L e N e d e = D e B e T I e d
Równanie równowagi
E
X
e=1
I
T eT K ¯ e T I e d =
E
X
e=1
I T eT p ¯ e b +
E
X
e=1
I T eT p ¯ e
Metody obliczeniowe, 2020 P.Plucińskic
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e=
E
X
e=1
I T
eTZ
Se
N
eTt
edS
e+
Z
Ve
N
eTf
edV
eUwzględnienie związków kinematycznych i konstytutywnych liniowa sprężystość: σ = Dε
liniowy związek kinematyczny: ε = Lu
σ e = D e L e u e = D e L e N e d e = D e B e T I e d
Równanie równowagi
E
X eT ¯ e e
E
X eT e
E
X eT e
Równanie równowagi dla układu zdyskretyzowanego
Równanie równowagi
E
X
e=1
I T
eTZ
Ve
B
eTσ
edV
e=
E
X
e=1
I T
eTZ
Se
N
eTt
edS
e+
Z
Ve
N
eTf
edV
eUwzględnienie związków kinematycznych i konstytutywnych liniowa sprężystość: σ = Dε
liniowy związek kinematyczny: ε = Lu
σ e = D e L e u e = D e L e N e d e = D e B e T I e d
Równanie równowagi
Kd = p b + p
Metody obliczeniowe, 2020 P.Plucińskic
Płaski stan naprężenia (σ z = 0)
Wektor funkcji przemieszczeń u = {u(x, y), v(x, y)}
Wektor odkształcenia ε = {ε x , ε y , γ xy }
Wektor naprężenia σ = {σ x , σ y , τ xy }
Wektor intensywności sił
Wektor obciążenia po długości elementu f = {f x , f y }
Macierz związków konstytutywnych
D = E 1 − ν 2
1 ν 0
ν 1 0
0 0 1−ν 2
Macierz operatorów różniczkowych
∂
∂x 0
∂
Płaski stan naprężenia (σ z = 0)
Macierz sztywności
k e = Z
A
eB eT D e B e h e dA e
A e , h e – pole powierzchni i grubość ES
Wektor obciążenia elementu p e =
Z
A
eN eT f e h e dA e
Wektor sił brzegowych
p e b = Z
Γ
eN eT t e h e dΓ e
x e y e
Γ e A e
Metody obliczeniowe, 2020 P.Plucińskic
Elementy skończone dla tarczy
Element trójwęzłowy
u e (x, y) = N e (x, y) d e
N e = N i e 0 N j e 0 N k e 0 0 N i e 0 N j e 0 N k e
, d e =
d 1 d 2 d 3 d 4
d 5
d 6
x e y e
i
k
j d 1 e
d 2
d 3
d 4
d 5
d 6
Elementy skończone dla tarczy
Element trójwęzłowy
u e (x, y) = N e (x, y) d e
N e = N i e 0 N j e 0 N k e 0 0 N i e 0 N j e 0 N k e
, d e =
d 1 d 2 d 3 d 4
d 5
d 6
x e y e
i
k
j d 1 e
d 2
d 3
d 4
d 5
d 6
y
eN
i(x
e, y
e)
x
e1 i
k j
y
eN
j(x
e, y
e)
x
e1
i k
j y
eN
k(x
e, y
e)
x
e1 i
k
j
Metody obliczeniowe, 2020 P.Plucińskic
Tarczowe elementy skończone
Element czterowęzłowy
u e (x, y) = N e (x, y) d e N
e=
N
ie0 N
je0 N
ke0 N
le0 0 N
ie0 N
je0 N
ke0 N
led e = {d 1 , d 2 , d 3 , d 4 , d 5 , d 6 , d 7 , d 8 }
x e y e
i
k
j l
e d 1 d 2
d 3 d 4
d 5 d 6
d 7
d 8
Tarczowe elementy skończone
Element czterowęzłowy
u e (x, y) = N e (x, y) d e N
e=
N
ie0 N
je0 N
ke0 N
le0 0 N
ie0 N
je0 N
ke0 N
led e = {d 1 , d 2 , d 3 , d 4 , d 5 , d 6 , d 7 , d 8 }
x e y e
i
k
j l
e d 1 d 2
d 3 d 4
d 5 d 6
d 7 d 8
y
eN
i(x, y)
x
e1
i
j
k
l y
eN
l(x, y)
x
ei 1 j
k l y
eN
j(x, y)
x
e1 j i
k
l y
eN
k(x, y)
x
e1
i j
k l
Metody obliczeniowe, 2020 P.Plucińskic
Przykład
Statyka tarczy
3 kN/m
7.5 kN/m
4 m 2 m
2 m
E = 18 GPa ν = 0.25 h = 0.2 m
X Y
4 i
5 j
2 k 1
l
Dyskretyzacja
i j
3 k
1 2
d 1
d 2
d 3 d 4
d 5 d 6
d 7 d 8
d 9 d 10
nr elem. nr węzłów
1 4 5 2 1
2 5 3 2
Przykład
Statyka tarczy
3 kN/m
7.5 kN/m
4 m 2 m
2 m
E = 18 GPa ν = 0.25 h = 0.2 m
X Y
4 i
5 j
2 k 1
l
Dyskretyzacja
i j
3 k
1 2
d 1
d 2
d 3 d 4 d 5
d 6
d 7 d 8
d 9 d 10
nr elem. nr węzłów 1 4 5 2 1 2 5 3 2
Metody obliczeniowe, 2020 P.Plucińskic
Przykład
Statyka tarczy
3 kN/m
7.5 kN/m
E = 18 GPa
Y Dyskretyzacja
k 1
l j
3 k 2 d1
d2 d3 d4
d5 d6