Metod konačnih elemenata

Prvi ispitni zadatak

1. Razmatramo problem
−u '' ( x) = 7, x ∈ (0,1)
u (0) = 2, u (1) = 3
a) Odrediti slabu formu problema
b) Formulisati metod konačnih elemenata za zadati problem (matrica krutosti, vektor
slobodnih članova)
c) Rešiti sistem jednačina koristeći ravnomernu mrežu sa 4 čvora.

2. Razmatramo problem
−u '' ( x) = f ( x), x ∈ (0,1)
u (0) = u (1) = 0
a) Odrediti slabu formu problema
b) Formulisati metod konačnih elemenata za zadati problem koristeći ravnomernu mrežu
sa 3 čvora ( h = 1/ 2 ).

3. Razmatramo problem
−u '' ( x) + u ( x) = f ( x), x ∈ (0,1)
u (0) = u (1) = 0
a) Odrediti slabu formu problema
b) Formulisati metod konačnih elemenata za zadati problem koristeći ravnomernu mrežu
sa 5 čvora ( h = 1/ 4 ).

4. Razmatramo problem
−u '' ( x) + u ( x) = f ( x), x ∈ (0,1)
u (0) = u (1) = 0
a) Odrediti slabu formu problema
b) Formulisati metod konačnih elemenata za zadati problem koristeći ravnomernu mrežu
sa 5 čvora ( h = 1/ 4 )
c) Ako je u ( x) = sin(π x) odrediti f ( x) i rešititi sistem jednačina datog aproksimiranog
problema
x(1 − x)
d) Ako je u ( x) = odrediti f ( x) i rešititi sistem jednačina datog aproksimiranog
2
problema.
5. Razmatramo problem
 −u '' ( x) = f ( x), x ∈ (0,1)
u (0) = u (1) = 0
a) Odrediti slabu formu problema
b) Formulisati metod konačnih elemenata za zadati problem
c) Ako je u ( x) = 2 x − x 2 odrediti f ( x) i rešiti sistem jednačina koristeći ravnomernu
mrežu sa 3 čvora.

6. Razmatramo problem
−u '' ( x) = f ( x), x ∈ (0,1)
u (0) = u (1) = 0

a) Odrediti slabu formu problema

b) Formulisati metod konačnih elemenata za zadati problem
c) Ako je f ( x) = x − 4 i rešiti sistem jednačina koristeći ravnomernu mrežu sa 3 čvora
d) Ako je f ( x) = 1 i rešiti sistem jednačina koristeći ravnomernu mrežu sa 3 čvora
Drugi ispitni zadatak

1. Za jednodimenzioni konačni element prikazan na slici 1

a) Izvesti izraze za interpolacione funkcije
b) Izvesti elemente matrice krutosti
c) Za nosač prikazan na slici 2 opterećen linijskim paraboličkim opterećenjem
p ( x) = a + bx + cx 2 odrediti nepoznata pomeranja i napone koristeći elemente sa slike
1, ako je E = 30GPa, A = 0.1m 2
u1 u2 u3

E, A

l l
2 2

p0 = 10 kNm

l=4

2. Za jednodimenzioni konačni element sa 3 čvora prikazan na slici

a) Izvesti izraze za interpolacione funkcije
b) Ako je l = 1, u1 = 0, u2 = 5 ⋅10 −6 m, u3 = 15 ⋅10−6 m sračunati pomeranje, dilataciju i
napon
u1 u2 u3

E, A

l
2

3. Za jednodimenzioni konačni element sa 3 čvora prikazan na slici

a) Izvesti izraze za interpolacione funkcije
 b) Izvesti elemente matrice krutosti na glavnoj dijagonali ako je promena poprečnog
x
preseka linearna A = A1 + ( A3 − A1 )
l
u1 u2 u3

E , A1 A3

l
2

4. Za jednodimenzioni konačni element sa 2 čvora prikazan na slici

a) Izvesti izraze za interpolacione funkcije
x
b)Izvesti elemente matrice krutosti ako je E = 30GPa, A = 0.0225m 2 , p ( x) = .
2

x
5 kNm
p( x) = 2

u1 = 0 u2
E , A = const

l = 10m
5. Odrediti nepoznata pomeranja, reakcije oslonca i napone u slučaju:
a) 2 konačna elementa sa 2 stepena slobode
b) 1 konačni element sa 3 stepena slobode
ako je l = 4m, A = 0.6m 2 , p1 = p2 = 5 kN m , E = 30GPa
p1

p2

l
6. Za jednodimenzioni konačni element sa 2 čvora prikazan na slici izvesti izraze za
interpolacione funkcije
u1 u2
ϕ1
ϕ2

l
7. Sračunati pomeranja, deformaciju, napone, reakcije oslonca primenom metode konačnih
elemenata ukoliko diskretizaciju vršimo sa dva konačna elementa sa dva čvora.
p0

Treći ispitni zadatak

1. U pravougaonom konačnom elementu prikazanom na slici vlada konstantno stanje
T T kN
napona čiji je vektor σ x σ y τ xy  = [ 200000 200000 100000] 2 . Sračunati
m
pomeranja čvorova konačnog elementa ako u njemu vlada ravno stanje napona,
E = 210GPa, ν = 0.3 .

v3
v4
u3 u4
3 4

b = 0.7 m

v2
v1 u1 u2
1 2

a = 1m
2. U trougaonom konačnom elementu prikazanom na slici vlada konstantno stanje
T
T 100 −6 60 −6 
deformacije, čiji je vektor ε x ε y γ xy  =  10 10 0  . Sračunati pomeranja
 7 7 
čvorova konačnog elementa.
 y
v3
u3
3(0.2,0.3)

v2
u2
2(0.3, 0.1)
v1

u1
x
1(0,0)
3. Na slici je prikazan trougaoni konačni element. Moduo elastičnosti iznosi E = 30GPa ,
dok je ν = 0.2 . U tabeli su date karakteristike konačnog elementa, za koji je potrebno sračunati
vektor napona i deformacije.

x[m] y[m] u[10 −6 ] v[10 −6 ]

1 3 3 0.5 0.5
2 10 2 -1 -3
3 6 6 2 3

3(6, 6)
u3

v3

v1

u12 1(3,3)

v2
2(10, 2)
u2