OSNOVI METODE KONAČNIH

ELEMENATA

Rešeni primeri

2024. godina
 Zadatak 1

Odrediti raspodelu kinematičkih i statičkih veličina u polju konačnih elemenata. Svaki raspon
diskretizovati sa jednim linijskim konačnim elementom (Ojler-Bernulijeva teorija savijanja) i usvojiti
savijanje u X – Y ravni.
Podaci: E = 30 GPa; b/h = 30/70 cm; L1 = 6,0 m; L2 = 8,0 m; q = 30,0 kN/m; P = 50 kN.

x1 ≔ 0 y1 ≔ 0 E ≔ 30 ⋅ 10 6

x2 ≔ 6 y2 ≔ 0 b ≔ 0.3

x3 ≔ 14 y3 ≔ 0 d ≔ 0.7

Nepoznata generalisana pomeranja: 3, 4, 6, 7 i 9

KE 1
čvorovi - mreža: 1, 2
čvorovi - KE: 1, 2
SS GKS: 1, 2, 3, 4, 5, 6
SS KE: 1, 2, 3, 4, 5, 6
Matrica transformacije

x1_1 ≔ x1 = 0.0 x1_2 ≔ x2 = 6.0 y1_1 ≔ y1 = 0.0 y1_2 ≔ y2 = 0.0

L1 ≔ ‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾
((x1_2 - x1_1)) 2 + ((y1_2 - y1_1)) 2 = 6.0

x1_2 - x1_1 y1_2 - y1_1

cos_α_KE1 ≔ ――――= 1.0 sin_α_KE1 ≔ ――――= 0.0
L1 L1
Napomena:
matrica transformacije je jedinična pa su veličine u GKS jednake veličinama u LKS

Matrica krutosti

Napomena:
određeni su samo elementi matrice krutosti KE uz nepoznata generalisana pomeranja

b ⋅ d3
A1 ≔ b ⋅ d = 0.21000 I1 ≔ ――= 0.00858
12

4 ⋅ E ⋅ I1 2 ⋅ E ⋅ I1
k1_33 ≔ ――― = 171500.0 k1_34 ≔ 0.0 k1_36 ≔ ――― = 85750.0
L1 L1

E ⋅ A1
k1_43 ≔ k1_34 = 0.0 k1_44 ≔ ―― = 1050000.0 k1_46 ≔ 0.0
L1

1
 E ⋅ A1
Zadatak
k1_44 ≔ ―― =1 1050000.0
L1

4 ⋅ E ⋅ I1
k1_63 ≔ k1_36 = 85750.0 k1_64 ≔ k1_46 = 0.0 k1_66 ≔ ――― = 171500.0
L1
3 4 6 7 9
⎡ k1_33 k1_34 k1_36 0 0 ⎤ ⎡ 171500.0 0.0 85750.0 0.0 0.0 ⎤ ⎡3⎤
⎢ k1_43 k1_44 k1_46 0 0 ⎥ ⎢ 0.0 1050000.0 0.0 0.0 0.0 ⎥ ⎢4⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥
k1_aa ≔ ⎢ k1_63 k1_64 k1_66 0 0 ⎥ = ⎢ 85750.0 0.0 171500.0 0.0 0.0 ⎥ ⎢6⎥
⎢ 0 0 0 0 0⎥ ⎢ 0.0 0.0 0.0 0.0 0.0 ⎥ ⎢7⎥
⎢⎣ 0 0 0 0 0 ⎥⎦ ⎢⎣ 0.0 0.0 0.0 0.0 0.0 ⎥⎦ ⎢⎣ 9 ⎥⎦

Ekvivalentno opterećenje q ≔ -30.0

⎡ q ⋅ L1 2 ⎤
⎢ ――― ⎥
⎢ 12 ⎥ ⎡ -90.00 ⎤ ⎡3⎤
⎢ 0 ⎥ ⎢ ⎥ ⎢4⎥
⎢ q ⋅ L1 2 ⎥ ⎢ 0.00 ⎥ ⎢ ⎥
Q1_a ≔ ⎢ -――― ⎥ = 90.00 ⎢6⎥
⎢ 12 ⎥ ⎢ 0.00 ⎥
⎢ ⎥ ⎢7⎥
⎢ 0 ⎥ ⎢⎣ 0.00 ⎥⎦ ⎢⎣ 9 ⎥⎦
⎢⎣ 0 ⎥⎦

ili

2 ⋅ x2 x 3 expand x 3 x2 x2 x 3 expand x 3 x2
N1_3 ((x)) ≔ x - ――+ ―― ――― → ― ―+x
- N1_6 ((x)) ≔ -―― + ―― ――― → ― ― -
L1 L1 2 36 3 L1 L1 2 36 6

L1 L1

Q1_3 ≔ ⌠
⌡ N1_3 ((x)) ⋅ q d x = -90.00 Q1_6 ≔ ⌠
⌡ N1_6 ((x)) ⋅ q d x = 90.00
0 0

KE 2
čvorovi - mreža: 2, 3
čvorovi - KE: 1, 2
SS GKS: 4, 5, 6, 7, 8, 9
SS KE: 1, 2, 3, 4, 5, 6

Matrica transformacije

x2_1 ≔ x2 = 6.0 x2_2 ≔ x3 = 14.0 y2_1 ≔ y2 = 0.0 y2_2 ≔ y3 = 0.0

L2 ≔ ‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾
((x2_2 - x2_1)) 2 + ((y2_2 - y2_1)) 2 = 8.0

x2_2 - x2_1 y2_2 - y2_1

cos_α_KE2 ≔ ――――= 1.0 sin_α_KE2 ≔ ――――= 0.0
L2 L2
Napomena:
matrica transformacije je jedinična pa su veličine u GKS jednake veličinama u LKS

2
 Zadatak 1

Matrica krutosti

Napomena:
određeni su samo elementi matrice krutosti KE uz nepoznata generalisana pomeranja

b ⋅ d3
A2 ≔ b ⋅ d = 0.21000 I2 ≔ ――= 0.00858
12

E ⋅ A2 E ⋅ A2
k2_11 ≔ ―― = 787500.0 k2_13 ≔ 0.0 k2_14 ≔ -―― = -787500.0 k2_16 ≔ 0.0
L2 L2

4 ⋅ E ⋅ I2 2 ⋅ E ⋅ I2
k2_31 ≔ k2_13 = 0.0 k2_33 ≔ ――― = 128625.0 k2_34 ≔ 0.0 k2_36 ≔ ――― = 64312.5
L2 L2

E ⋅ A2
k2_41 ≔ k2_14 = -787500.0 k2_43 ≔ k2_34 = 0.0 k2_44 ≔ ―― = 787500.0 k2_46 ≔ 0.0
L2
k2_61 ≔ k2_16 = 0.0 k2_63 ≔ k2_36 = 64312.5
4 ⋅ E ⋅ I2
k2_64 ≔ k2_46 = 0.0 k2_66 ≔ ――― = 128625.0
L2
3 4 6 7 9
⎡0 0 0 0 0 ⎤ ⎡ 0.0 0.0 0.0 0.0 0.0 ⎤ ⎡3⎤
⎢0 k2_11 k2_13 k2_14 k2_16 ⎥ ⎢ 0.0 787500.0 0.0 -787500.0 0.0 ⎥ ⎢4⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥
k2_aa ≔ ⎢ 0 k2_31 k2_33 k2_34 k2_36 ⎥ = ⎢ 0.0 0.0 128625.0 0.0 64312.5 ⎥ ⎢6⎥
⎢0 k2_41 k2_43 k2_44 k2_46 ⎥ ⎢ 0.0 -787500.0 0.0 787500.0 0.0 ⎥ ⎢7⎥
⎢⎣ 0 k2_61 k2_63 k2_64 k2_66 ⎥⎦ ⎢⎣ 0.0 0.0 64312.5 0.0 128625.0 ⎥⎦ ⎢⎣ 9 ⎥⎦

Ekvivalentno opterećenje P ≔ -50

⎡ 0 ⎤
⎢ 0 ⎥ ⎡ 0.00 ⎤ ⎡3⎤
⎢ ⎥ ⎢ 0.00 ⎥ ⎢4⎥
P ⋅ L2 ⎢ ⎥ ⎢ ⎥
Q2_a ≔ ⎢ ―― ⎥ = ⎢ -50.00 ⎥ ⎢6⎥
⎢ 8 ⎥ 0.00 ⎥
⎢ ⎥ ⎢ ⎢7⎥
0
⎢ P ⋅ L2 ⎥ ⎢⎣ 50.00 ⎥⎦ ⎢⎣ 9 ⎥⎦
⎢ -――⎥
⎣ 8 ⎦

ili
2 ⋅ x2 x 3 expand x 3 x2 x2 x 3 expand x 3 x2
N2_3 ((x)) ≔ x - ――+ ―― ――― → ― - ― +x N2_6 ((x)) ≔ -―― + ―― ――― → ― - ―
L2 L2 2 64 4 L2 L2 2 64 8

⎛ L2 ⎞ ⎛ L2 ⎞
Q2_6 ≔ N2_3 ⎜―― ⎟ ⋅ P = -50.00 Q2_9 ≔ N2_6 ⎜―― ⎟ ⋅ P = 50.00
⎝ 2 ⎠ ⎝ 2 ⎠

3
 Zadatak 1

Sistem KE

3 4 6 7 9
⎡ 171500.0 0.0 85750.0 0.0 0.0 ⎤ ⎡3⎤
⎢ 0.0 1837500.0 0.0 -787500.0 0.0 ⎥ ⎢4⎥
⎢ ⎥ ⎢ ⎥
Kaa ≔ k1_aa + k2_aa = ⎢ 85750.0 0.0 300125.0 0.0 64312.5 ⎥ ⎢6⎥
⎢ 0.0 -787500.0 0.0 787500.0 0.0 ⎥ ⎢7⎥
⎢⎣ 0.0 0.0 64312.5 0.0 128625.0 ⎥⎦ ⎢⎣ 9 ⎥⎦

⎡ -90.00 ⎤ ⎡3⎤ ⎡ -0.00065806 ⎤ ⎡3⎤

⎢ 0.00 ⎥ ⎢4⎥ ⎢ 0.00000000 ⎥ ⎢4⎥
⎢ ⎥ ⎢ ⎥ -1 ⎢ ⎥ ⎢ ⎥
Sa ≔ Q1_a + Q2_a = ⎢ 40.00 ⎥ ⎢6⎥ da ≔ Kaa ⋅ Sa = ⎢ 0.00026656 ⎥ ⎢6⎥
⎢ 0.00 ⎥ ⎢7⎥ ⎢ 0.00000000 ⎥ ⎢7⎥
⎢⎣ 50.00 ⎥⎦ ⎢⎣ 9 ⎥⎦ ⎢⎣ 0.00025545 ⎥⎦ ⎢⎣ 9 ⎥⎦

Raspodela kinematičkih i statičkih veličina - KE 1

u1_1 ≔ 0.0 u1_2 ≔ 0.0

v1_1 ≔ 0.0 v1_2 ≔ 0.0

φ1_1 ≔ da = -0.00065806 φ1_2 ≔ da = 0.00026656

0 2

⎡ v1_1 ⎤ ⎡ 0.00000000 ⎤
⎢ φ1_1 ⎥ ⎢ -0.00065806 ⎥
d1 ≔ ⎢ ⎥=⎢ ⎥
⎢ v1_2 ⎥ ⎢ 0.00000000 ⎥
⎣ φ1_2 ⎦ ⎣ 0.00026656 ⎦

v1 ((x)) ≔ N1_3 ((x)) ⋅ φ1_1 + N1_6 ((x)) ⋅ φ1_2

float , 5
v1 ((x)) ――― → -0.000010875 ⋅ x 3 + 0.00017493 ⋅ x 2 - 0.00065806 ⋅ x

d float , 5
φ1 ((x)) ≔ ―― v1 ((x)) ――― → -0.000032625 ⋅ x 2 + 0.00034985 ⋅ x - 0.00065806
dx

⎡ 6 12 ⋅ x 4 6⋅x 6 12 ⋅ x 2 6⋅x ⎤
B1 ((x)) ≔ ⎢ -―― 2
+ ―― -―― + ―― ―― - ―― -―― + ―― ⎥
⎣ L1 L1 3 L1 L1 2 L1 2 L1 3 L1 L1 2 ⎦

⎡ x 1 x 2 x 1 x 1⎤
B1 ((x)) → ⎢ ― - ――- ―-― + ――- ―⎥
⎣ 18 6 6 3 18 6 6 3 ⎦

D1 ≔ E ⋅ I1 = 257250.0

float , 5
M1 ((x)) ≔ D1 ⋅ B1 ((x)) ⋅ d1 ――― → -16.786 ⋅ x + 90.0

4
 Zadatak 1

M1 ((0)) = 90.000 M1 ((L1)) = -10.716

6 ⋅ ((2 ⋅ v1_1 - 2 ⋅ v1_2 + L1 ⋅ ((φ1_1 + φ1_2))))

T1 ≔ -E ⋅ I1 ⋅ ――――――――――――― = 16.786
L1 3

ili

⎛ d ⎞ float , 5 ⎛ d3 ⎞ float , 5
T1 ≔ -⎜―― M1 ((x))⎟ ――― → 16.786 T1 ≔ -E ⋅ I1 ⋅ ⎜―― v1 ((x))⎟ ――― → 16.786
3
⎜⎝ dx ⎟⎠ ⎜⎝ d x ⎟⎠

A ≔ A1 = 0.210 I ≔ I1 = 0.00858 L ≔ L1 = 6.0

⎡ E⋅A E⋅A ⎤
⎢ ―― 0 0 -―― 0 0 ⎥
⎢ L L ⎥
⎢ 0 12 ⋅ E ⋅ I 6⋅E⋅I 12 ⋅ E ⋅ I 6 ⋅ E ⋅ I ⎥
――― ――― 0 -――― ―――
⎢ L3 L2 L3 L2 ⎥
⎢ ⎥
⎢ 0 6⋅E⋅I 4⋅E⋅I 6⋅E⋅I 2⋅E⋅I ⎥
――― ――― 0 -――― ―――
⎢ L2 L L2 L ⎥
k1 ≔ ⎢ ⎥
E⋅A E⋅A
⎢ -―― 0 0 ―― 0 0 ⎥
⎢ L L ⎥
⎢ 12 ⋅ E ⋅ I 6⋅E⋅I 12 ⋅ E ⋅ I 6⋅E⋅I ⎥
⎢ 0 -――― 3
-――― 2
0 ――― 3
-――― ⎥
⎢ L L L L2 ⎥
⎢ 6⋅E⋅I 2⋅E⋅I 6⋅E⋅I 4⋅E⋅I ⎥
0 ――― ――― 0 -――― ―――
⎢ L 2
L L 2
L ⎥
⎣ ⎦

⎡ 10.50000 0.00000 0.00000 -10.50000 0.00000 0.00000 ⎤

⎢ 0.00000 0.14292 0.42875 0.00000 -0.14292 0.42875 ⎥
⎢ ⎥
k1
―― = ⎢ 0.00000 0.42875 1.71500 0.00000 -0.42875 0.85750 ⎥
10 5 ⎢ -10.50000 0.00000 0.00000 10.50000 0.00000 0.00000 ⎥
⎢ 0.00000 -0.14292 -0.42875 0.00000 0.14292 -0.42875 ⎥
⎢⎣ 0.00000 0.42875 0.85750 0.00000 -0.42875 1.71500 ⎥⎦

⎡ 0 ⎤
⎢ q ⋅ L1 ⎥
⎢ ―― ⎥
⎡ u1_1 ⎤ ⎡ 0.00000000 ⎤ ⎢ 2 ⎥ ⎡ 0.0 ⎤ ⎡ 0.000 ⎤
⎢ v1_1 ⎥ ⎢ 0.00000000 ⎥ ⎢ q ⋅ L1 2 ⎥ ⎢ -90.0 ⎥ ⎢ 73.214 ⎥
⎢ ⎥ ⎢ ⎥ ⎢ ――― ⎥ ⎢ ⎥ ⎢ ⎥
φ1_1 ⎥ ⎢ -0.00065806 ⎥ ⎢ 12 ⎥ = ⎢ -90.0 ⎥ 0.000 ⎥
d1 ≔ ⎢ = Q1 ≔ R1 ≔ k1 ⋅ d1 - Q1 = ⎢
⎢ u1_2 ⎥ ⎢ 0.00000000 ⎥ ⎢ 0 ⎥ ⎢ 0.0 ⎥ ⎢ 0.000 ⎥
⎢ v1_2 ⎥ ⎢ 0.00000000 ⎥ ⎢ q ⋅ L1 ⎥ ⎢ -90.0 ⎥ ⎢ 106.786 ⎥
⎢⎣ φ1_2 ⎥⎦ ⎢⎣ 0.00026656 ⎥⎦ ⎢ ―― 2 ⎥ ⎢⎣ 90.0 ⎥⎦ ⎢⎣ -100.714 ⎥⎦
⎢ 2
⎥
⎢ q ⋅ L1 ⎥
-―――
⎢⎣ 12 ⎥⎦

5
 Zadatak 1

Raspodela kinematičkih i statičkih veličina - KE 2

u2_1 ≔ 0.0 u2_2 ≔ 0.0

v2_1 ≔ 0.0 v2_2 ≔ 0.0

φ2_1 ≔ da = 0.00026656 φ2_2 ≔ da = 0.00025545

2 4

⎡ v2_1 ⎤ ⎡ 0.00000000 ⎤
⎢ φ2_1 ⎥ ⎢ 0.00026656 ⎥
d2 ≔ ⎢ ⎥=⎢ ⎥
⎢ v2_2 ⎥ ⎢ 0.00000000 ⎥
⎣ φ2_2 ⎦ ⎣ 0.00025545 ⎦

v2 ((x)) ≔ N2_3 ((x)) ⋅ φ2_1 + N2_6 ((x)) ⋅ φ2_2

float , 5
v2 ((x)) ――― → 0.0000081563 ⋅ x 3 - 0.00009857 ⋅ x 2 + 0.00026656 ⋅ x

d float , 5
φ2 ((x)) ≔ ―― v2 ((x)) ――― → 0.000024469 ⋅ x 2 - 0.00019714 ⋅ x + 0.00026656
dx

⎡ 6 12 ⋅ x 4 6⋅x 6 12 ⋅ x 2 6⋅x ⎤
B2 ((x)) ≔ ⎢ -―― 2
+ ―― 3
-―― + ―― ―― - ―― -―― + ―― ⎥
⎣ L2 L2 L2 L2 L2
2 2
L2 3
L2 L2 2 ⎦

⎡ 3⋅x 3 3⋅x 1 3⋅x 3 3⋅x 1 ⎤

B2 ((x)) → ⎢ ―― - ― ―― - ―-―― + ― ―― - ―⎥
⎣ 128 32 32 2 128 32 32 4⎦

D2 ≔ E ⋅ I2 = 257250.0

float , 5
M2 ((x)) ≔ D2 ⋅ B2 ((x)) ⋅ d2 ――― → 12.589 ⋅ x - 50.714

M2 ((0)) = -50.714 M2 ((L2)) = 49.998

6 ⋅ ((2 ⋅ v2_1 - 2 ⋅ v2_2 + L2 ⋅ ((φ2_1 + φ2_2))))

T2 ≔ -E ⋅ I2 ⋅ ――――――――――――― = -12.589
L2 3

ili

⎛ d ⎞ float , 5 ⎛ d3 ⎞ float , 5
T2 ≔ -⎜―― M2 ((x))⎟ ――― → -12.589 T2 ≔ -E ⋅ I2 ⋅ ⎜―― v2 ((x))⎟ ――― → -12.589
3
⎜⎝ dx ⎟⎠ ⎜⎝ d x ⎟⎠

A ≔ A2 = 0.210 I ≔ I2 = 0.00858 L ≔ L2 = 8.0

6
 Zadatak 1

⎡ E⋅A E⋅A ⎤
⎢ ―― 0 0 -―― 0 0 ⎥
⎢ L L ⎥
⎢ 0 12 ⋅ E ⋅ I 6⋅E⋅I 12 ⋅ E ⋅ I 6 ⋅ E ⋅ I ⎥
――― ――― 0 -――― ―――
⎢ L3 L2 L3 L2 ⎥
⎢ ⎥
⎢ 0 6⋅E⋅I 4⋅E⋅I 6⋅E⋅I 2⋅E⋅I ⎥
――― ――― 0 -――― ―――
⎢ L2 L L2 L ⎥
k2 ≔ ⎢ ⎥
E⋅A E⋅A
⎢ -―― 0 0 ―― 0 0 ⎥
⎢ L L ⎥
⎢ 12 ⋅ E ⋅ I 6⋅E⋅I 12 ⋅ E ⋅ I 6⋅E⋅I ⎥
⎢ 0 -――― 3
-――― 2
0 ――― 3
-――― ⎥
⎢ L L L L2 ⎥
⎢ 6⋅E⋅I 2⋅E⋅I 6⋅E⋅I 4⋅E⋅I ⎥
0 ――― ――― 0 -――― ―――
⎢ L 2
L L 2
L ⎥
⎣ ⎦

⎡ 7.87500 0.00000 0.00000 -7.87500 0.00000 0.00000 ⎤

⎢ 0.00000 0.06029 0.24117 0.00000 -0.06029 0.24117 ⎥
⎢ ⎥
k2 0.00000 0.24117 1.28625 0.00000 -0.24117 0.64313 ⎥
―― = ⎢
10 5 ⎢ -7.87500 0.00000 0.00000 7.87500 0.00000 0.00000 ⎥
⎢ 0.00000 -0.06029 -0.24117 0.00000 0.06029 -0.24117 ⎥
⎢⎣ 0.00000 0.24117 0.64313 0.00000 -0.24117 1.28625 ⎥⎦

⎡ 0 ⎤
⎢ P ⎥
⎢ ― ⎥
⎡ u2_1 ⎤ ⎡ 0.00000000 ⎤ ⎢ 2 ⎥ ⎡ 0.0 ⎤ ⎡ 0.000 ⎤
⎢ v2_1 ⎥ ⎢ 0.00000000 ⎥ ⎢ P ⋅ L2 ⎥ ⎢ -25.0 ⎥ ⎢ 37.589 ⎥
⎢ ⎥ ⎢
φ2_1 ⎥ ⎢ 0.00026656 ⎥
⎥ ⎢ ――8 ⎥ ⎢ -50.0 ⎥ ⎢ ⎥
100.714 ⎥
d2 ≔ ⎢ = Q2 ≔ ⎢ ⎥=⎢ ⎥ R2 ≔ k2 ⋅ d2 - Q2 = ⎢
⎢ u2_2 ⎥ ⎢ 0.00000000 ⎥ ⎢ 0 ⎥ ⎢ 0.0 ⎥ ⎢ 0.000 ⎥
⎢ v2_2 ⎥ ⎢ 0.00000000 ⎥ ⎢ P ⎥ ⎢ -25.0 ⎥ ⎢ 12.411 ⎥
⎢⎣ φ2_2 ⎥⎦ ⎢⎣ 0.00025545 ⎥⎦ ⎢ ― 2 ⎥ ⎢⎣ 50.0 ⎥⎦ ⎢⎣ 0.000 ⎥⎦
⎢ ⎥
P ⋅ L2 ⎥
⎢ -――
⎢⎣ 8 ⎥⎦

7
 Zadatak 2

Odrediti raspodelu pomeranja, deformacije i napona. Numeričku analizu sprovesti sa dva trougaona
konačna elementa prema skici (CST konačni element; ravansko stanje deformacije). Podaci:
E = 30 GPa; ν = 0,2; L1 = 6,0 m; L2 = 6,0 m; L3 = 3,0 m; q1 = 10,0 MPa; q2 = 5,0 MPa.

E ≔ 30 ⋅ 10 6 ν ≔ 0.2 h ≔ 1.0 Nepoznata pomeranja: 5, 6, 7 i 8

x1 ≔ 0 x2 ≔ 6 x3 ≔ 3 x4 ≔ 0

y1 ≔ 0 y2 ≔ 0 y3 ≔ 6 y4 ≔ 6

⎡1-ν ν 0 ⎤ ⎡ 33333333.3 8333333.3 0.0 ⎤

E ⎢ ⎥ ⎢
D ≔ ―――――― ⋅ ν 1-ν 0 = 8333333.3 33333333.3 0.0 ⎥
((1 + ν)) ⋅ ((1 - 2 ⋅ ν)) ⎢ ⎥ ⎢ ⎥
⎢ 0 1-2⋅ν ⎥ ⎣ 0.0 0.0 12500000.0 ⎦
0 ―――
⎢⎣ 2 ⎥⎦

CST KE 1
čvorovi - mreža: 1, 2, 3
čvorovi - KE: 1, 2, 3
SS GKS: 1, 2, 3, 4, 5, 6
SS KE: 1, 2, 3, 4, 5, 6

x1_1 ≔ x1 = 0.0 x1_2 ≔ x2 = 6.0 x1_3 ≔ x3 = 3.0

y1_1 ≔ y1 = 0.0 y1_2 ≔ y2 = 0.0 y1_3 ≔ y3 = 6.0

a1_1 ≔ x1_2 ⋅ y1_3 - x1_3 ⋅ y1_2 = 36.0 a1_2 ≔ x1_3 ⋅ y1_1 - x1_1 ⋅ y1_3 = 0.0
1
a1_3 ≔ x1_1 ⋅ y1_2 - x1_2 ⋅ y1_1 = 0.0 Ae1 ≔ ―⋅ ((a1_1 + a1_2 + a1_3)) = 18.0
2

b1_1 ≔ y1_2 - y1_3 = -6.0 b1_2 ≔ y1_3 - y1_1 = 6.0 b1_3 ≔ y1_1 - y1_2 = 0.0

c1_1 ≔ x1_3 - x1_2 = -3.0 c1_2 ≔ x1_1 - x1_3 = -3.0 c1_3 ≔ x1_2 - x1_1 = 6.0

1
 Zadatak 2

Matrica krutosti

Napomena:
određeni su samo elementi matrice krutosti KE uz nepoznata generalisana pomeranja

h
k1_55 ≔ ――― ⋅ ⎛b1_3 2 ⋅ D + c1_3 2 ⋅ D ⎞ = 6250000.0
4 ⋅ Ae1 ⎜⎝ 0,0 2 , 2⎟
⎠

h
k1_56 ≔ ――― ⋅ ⎛b1_3 ⋅ c1_3 ⋅ ⎛D + D ⎞⎞ = 0.0
4 ⋅ Ae1 ⎝ ⎝ 0,1 2 , 2⎠⎠

k1_65 ≔ k1_56 = 0.0

h
k1_66 ≔ ――― ⋅ ⎛c1_3 2 ⋅ D + b1_3 2 ⋅ D ⎞ = 16666666.7
4 ⋅ Ae1 ⎜⎝ 1,1 2 , 2⎟
⎠
5 6 78 5 6 7 8
⎡ k1_55 k1_56 0 0 ⎤ ⎡ 5 ⎤ ⎡ 62.50000 0.00000 0.00000 0.00000 ⎤ ⎡ 5 ⎤
⎢ k1_65 k1_66 0 0 ⎥ ⎢ 6 ⎥ k_1_aa ⎢ 0.00000 166.66667 0.00000 0.00000 ⎥ ⎢ 6 ⎥
k_1_aa ≔ ⎢ ⎥ ⎢ ⎥ ――― =⎢ ⎥⎢ ⎥
⎢ 0 0 0 0⎥ ⎢7⎥ 10 5 ⎢ 0.00000 0.00000 0.00000 0.00000 ⎥ ⎢ 7 ⎥
⎣ 0 0 0 0⎦ ⎣8⎦ ⎣ 0.00000 0.00000 0.00000 0.00000 ⎦ ⎣ 8 ⎦

ili određivanje cele matrice krutosti

⎡ b1_1 0 b1_2 0 b1_3 0 ⎤ ⎡ -0.167 0.000 0.167 0.000 0.000 0.000 ⎤

1
B1 ≔ ――― ⋅ ⎢ 0 c1_1 0 c1_2 0 c1_3 ⎥ = ⎢ 0.000 -0.083 0.000 -0.083 0.000 0.167 ⎥
2 ⋅ Ae1 ⎢ c1_1 b1_1 c1_2 b1_2 c1_3 b1_3 ⎥ ⎢ -0.083 -0.167 -0.083 0.167 0.167 0.000 ⎥
⎣ ⎦ ⎣ ⎦

k_1 ≔ B1 T ⋅ D ⋅ B1 ⋅ h ⋅ Ae1

1 2 3 4 5 6
⎡ 182.29167 52.08333 -151.04167 -10.41667 -31.25000 -41.66667 ⎤ ⎡ 1 ⎤
⎢ 52.08333 104.16667 10.41667 -20.83333 -62.50000 -83.33333 ⎥ ⎢ 2 ⎥
⎢ ⎥ ⎢ ⎥
k_1
―― = ⎢ -151.04167 10.41667 182.29167 -52.08333 -31.25000 41.66667 ⎥ ⎢ 3 ⎥
10 5
⎢ -10.41667 -20.83333 -52.08333 104.16667 62.50000 -83.33333 ⎥ ⎢ 4 ⎥
⎢ -31.25000 -62.50000 -31.25000 62.50000 62.50000 0.00000 ⎥ ⎢ 5 ⎥
⎢⎣ -41.66667 -83.33333 41.66667 -83.33333 0.00000 166.66667 ⎥⎦ ⎢⎣ 6 ⎥⎦

CST KE 2
čvorovi - mreža: 1, 3, 4
čvorovi - KE: 1, 2, 3
SS GKS: 1, 2, 5, 6, 7, 8
SS KE: 1, 2, 3, 4, 5, 6

x2_1 ≔ x1 = 0.0 x2_2 ≔ x3 = 3.0 x2_3 ≔ x4 = 0.0

y2_1 ≔ y1 = 0.0 y2_2 ≔ y3 = 6.0 y2_3 ≔ y4 = 6.0

a2_1 ≔ x2_2 ⋅ y2_3 - x2_3 ⋅ y2_2 = 18.0 a2_2 ≔ x2_3 ⋅ y2_1 - x2_1 ⋅ y2_3 = 0.0
1
a2_3 ≔ x2_1 ⋅ y2_2 - x2_2 ⋅ y2_1 = 0.0 Ae2 ≔ ―⋅ ((a2_1 + a2_2 + a2_3)) = 9
2

2
 1
Zadatak
Ae2 ≔ ―2⋅ ((a2_1 + a2_2 + a2_3)) = 9
2

b2_1 ≔ y2_2 - y2_3 = 0.0 b2_2 ≔ y2_3 - y2_1 = 6.0 b2_3 ≔ y2_1 - y2_2 = -6.0

c2_1 ≔ x2_3 - x2_2 = -3.0 c2_2 ≔ x2_1 - x2_3 = 0.0 c2_3 ≔ x2_2 - x2_1 = 3.0

Matrica krutosti

Napomena:
određeni su samo elementi matrice krutosti KE uz nepoznata generalisana pomeranja

h
k2_33 ≔ ――― ⋅ ⎛b2_2 2 ⋅ D + c2_2 2 ⋅ D ⎞ = 33333333.3
4 ⋅ Ae2 ⎜⎝ 0,0 2 , 2⎟
⎠

h
k2_34 ≔ ――― ⋅ ⎛b2_2 ⋅ c2_2 ⋅ ⎛D + D ⎞⎞ = 0.0
4 ⋅ Ae2 ⎝ ⎝ 0,1 2 , 2⎠⎠

h
k2_35 ≔ ――― ⋅ ⎛b2_2 ⋅ b2_3 ⋅ D + c2_2 ⋅ c2_3 ⋅ D ⎞ = -33333333.3
4 ⋅ Ae2 ⎝ 0,0 2 , 2⎠

h
k2_36 ≔ ――― ⋅ ⎛b2_2 ⋅ c2_3 ⋅ D + b2_3 ⋅ c2_2 ⋅ D ⎞ = 4166666.7
4 ⋅ Ae2 ⎝ 0,1 2 , 2⎠

k2_43 ≔ k2_34 = 0.0

h
k2_44 ≔ ――― ⋅ ⎛c2_2 2 ⋅ D + b2_2 2 ⋅ D ⎞ = 12500000.0
4 ⋅ Ae2 ⎜⎝ 1,1 2 , 2⎟
⎠

h
k2_45 ≔ ――― ⋅ ⎛b2_3 ⋅ c2_2 ⋅ D + b2_2 ⋅ c2_3 ⋅ D ⎞ = 6250000.0
4 ⋅ Ae2 ⎝ 1,0 2 , 2⎠

h
k2_46 ≔ ――― ⋅ ⎛c2_2 ⋅ c2_3 ⋅ D + b2_2 ⋅ b2_3 ⋅ D ⎞ = -12500000.0
4 ⋅ Ae2 ⎝ 1,1 2 , 2⎠

k2_53 ≔ k2_35 = -33333333.3

k2_54 ≔ k2_45 = 6250000.0

h
k2_55 ≔ ――― ⋅ ⎛b2_3 2 ⋅ D + c2_3 2 ⋅ D ⎞ = 36458333.3
4 ⋅ Ae2 ⎜⎝ 0,0 2 , 2⎟
⎠

h
k2_56 ≔ ――― ⋅ ⎛b2_3 ⋅ c2_3 ⋅ ⎛D + D ⎞⎞ = -10416666.7
4 ⋅ Ae2 ⎝ ⎝ 0,1 2 , 2⎠⎠

k2_63 ≔ k2_36 = 4166666.7

k2_64 ≔ k2_46 = -12500000.0

k2_65 ≔ k2_56 = -10416666.7

h
k2_66 ≔ ――― ⋅ ⎛c2_3 2 ⋅ D + b2_3 2 ⋅ D ⎞ = 20833333.3
4 ⋅ Ae2 ⎜⎝ 1,1 2 , 2⎟
⎠

3
 Zadatak 2

5 6 7 8
⎡ k2_33 k2_34 k2_35 k2_36 ⎤ ⎡ 5 ⎤
⎢ k2_43 k2_44 k2_45 k2_46 ⎥ ⎢ 6 ⎥
k_2_aa ≔ ⎢ ⎥⎢ ⎥
⎢ k2_53 k2_54 k2_55 k2_56 ⎥ ⎢ 7 ⎥
⎣ k2_63 k2_64 k2_65 k2_66 ⎦ ⎣ 8 ⎦
5 6 7 8
⎡ 333.33333 0.00000 -333.33333 41.66667 ⎤ ⎡ 5 ⎤
k_2_aa ⎢ 0.00000 125.00000 62.50000 -125.00000 ⎥ ⎢ 6 ⎥
――― = ⎢ ⎥ ⎢ ⎥
10 5 ⎢ -333.33333 62.50000 364.58333 -104.16667 ⎥ ⎢ 7 ⎥
⎣ 41.66667 -125.00000 -104.16667 208.33333 ⎦ ⎣ 8 ⎦

ili određivanje cele matrice krutosti

⎡ b2_1 0 b2_2 0 b2_3 0 ⎤ ⎡ 0.000 0.000 0.333 0.000 -0.333 0.000 ⎤

1
B2 ≔ ――― ⋅ ⎢ 0 c2_1 0 c2_2 0 c2_3 ⎥ = ⎢ 0.000 -0.167 0.000 0.000 0.000 0.167 ⎥
2 ⋅ Ae2 ⎢ c2_1 b2_1 c2_2 b2_2 c2_3 b2_3 ⎥ ⎢ -0.167 0.000 0.000 0.333 0.167 -0.333 ⎥
⎣ ⎦ ⎣ ⎦

k_2 ≔ B2 T ⋅ D ⋅ B2 ⋅ h ⋅ Ae2

1 2 5 6 7 8
⎡ 31.25000 0.00000 0.00000 -62.50000 -31.25000 62.50000 ⎤ ⎡ 1 ⎤
⎢ 0.00000 83.33333 -41.66667 0.00000 41.66667 -83.33333 ⎥ ⎢ 2 ⎥
⎢ ⎥ ⎢ ⎥
k_2
―― = ⎢ 0.00000 -41.66667 333.33333 0.00000 -333.33333 41.66667 ⎥ ⎢ 5 ⎥
10 5 ⎢ -62.50000 0.00000 0.00000 125.00000 62.50000 -125.00000 ⎥ ⎢ 6 ⎥
⎢ -31.25000 41.66667 -333.33333 62.50000 364.58333 -104.16667 ⎥ ⎢ 7 ⎥
⎢⎣ 62.50000 -83.33333 41.66667 -125.00000 -104.16667 208.33333 ⎥⎦ ⎢⎣ 8 ⎥⎦

Ekvivalentno opterećenje
1 ⎛ q1 ⋅ h ⋅ L2 ⎞
q1 ≔ 10000.0 L2 ≔ 6.0 Q4x ≔ ―⋅ ⎜―――⎟ = 10000.0
3 ⎝ 2 ⎠
1 1
q2 ≔ 5000.0 L3 ≔ 3.0 Q4y ≔ -―⋅ q2 ⋅ h ⋅ L3 = -7500.0 Q3y ≔ -―⋅ q2 ⋅ h ⋅ L3 = -7500.0
2 2
ili

a) trougaono opterećenje
⎛⎛ q1 ⎞ ⎞ s
qx ((s)) ≔ ⎜⎜-―― ⎟ ⋅ s + q1⎟ ⋅ h Ns4 ((s)) ≔ ――
⎝⎝ L2 ⎠ ⎠ L2
L2

Q4x ≔ ⌠
⌡ Ns4 ((s)) ⋅ qx ((s)) d s = 10000.0
0

b) jednakopodeljeno opterećenje
s s
qy ((s)) ≔ -q2 ⋅ h Ns4 ((s)) ≔ 1 - ―― Ns3 ((s)) ≔ ――
L3 L3
L3 L3

Q4y ≔ ⌠
⌡ Ns4 ((s)) ⋅ qy ((s)) d s = -7500.0 Q3y ≔ ⌠
⌡ Ns3 ((s)) ⋅ qy ((s)) d s = -7500.0
0 0

4
 Zadatak 2

⎡ 0 ⎤ ⎡ 0.0 ⎤ ⎡ 5 ⎤
⎢ Q3y ⎥ ⎢ -7500.0 ⎥ ⎢ 6 ⎥
Q_2_a ≔ ⎢ ⎥=⎢ ⎥⎢ ⎥
⎢ Q4x ⎥ ⎢ 10000.0 ⎥ ⎢ 7 ⎥
⎣ Q4y ⎦ ⎣ -7500.0 ⎦ ⎣ 8 ⎦

Sistem KE

5 6 7 8
⎡ 39583333.3 0.0 -33333333.3 4166666.7 ⎤ ⎡ 5 ⎤
⎢ 0.0 29166666.7 6250000.0 -12500000.0 ⎥ ⎢ 6 ⎥
Kaa ≔ k_1_aa + k_2_aa = ⎢ ⎥ ⎢ ⎥
⎢ -33333333.3 6250000.0 36458333.3 -10416666.7 ⎥ ⎢ 7 ⎥
⎣ 4166666.7 -12500000.0 -10416666.7 20833333.3 ⎦ ⎣ 8 ⎦

⎡ 0.0 ⎤ ⎡5⎤ ⎡ 0.0012667 ⎤ ⎡5⎤

⎢ -7500.0 ⎥ ⎢6⎥ ⎢ -0.0007000 ⎥ ⎢6⎥
Sa ≔ Q_2_a = ⎢ ⎥ ⎢ ⎥ da ≔ Kaa -1 ⋅ Sa = ⎢ ⎥ ⎢ ⎥
⎢ 10000.0 ⎥ ⎢7⎥ ⎢ 0.0014667 ⎥ ⎢7⎥
⎣ -7500.0 ⎦ ⎣8⎦ ⎣ -0.0003000 ⎦ ⎣8⎦

Raspodela kinematičkih i statičkih veličina - KE 1

⎡ 0 ⎤ CST KE 1
⎢ 0 ⎥ ⎡ 0.0000000000 ⎤ ⎡1⎤
čvorovi - mreža: 1, 2, 3
⎢ ⎥ ⎢ 0.0000000000 ⎥ ⎢2⎥
čvorovi - KE: 1, 2, 3
⎢ 0 ⎥ ⎢ 0.0000000000 ⎥ ⎢ ⎥
d1 ≔ ⎢ 0 ⎥ = ⎢ ⎥ ⎢3⎥ SS GKS: 1, 2, 3, 4, 5, 6
⎢ da ⎥ ⎢ 0.0000000000 ⎥ ⎢4⎥ SS KE: 1, 2, 3, 4, 5, 6
⎢ 0
⎥ ⎢ 0.0012666667 ⎥ ⎢5⎥
⎢ da ⎥ ⎢⎣ -0.0007000000 ⎥⎦ ⎢⎣ 6 ⎥⎦
⎣ 1⎦

Raspodela pomeranja
1 float , 5
čvor 3 N1_3 ((x , y)) ≔ ――― ⋅ ((a1_3 + b1_3 ⋅ x + c1_3 ⋅ y)) ――― → 0.16667 ⋅ y
2 ⋅ Ae1

float , 5
u1 ((x , y)) ≔ N1_3 ((x , y)) ⋅ d1 ――― → 0.00021112 ⋅ y
4

float , 5
v1 ((x , y)) ≔ N1_3 ((x , y)) ⋅ d1 ――― → -0.00011667 ⋅ y
5

⎡ 0.0000000000 ⎤
Komponente deformacije: ε1 ≔ B1 ⋅ d1 = ⎢ -0.0001166667 ⎥
⎢ ⎥
⎣ 0.0002111111 ⎦

⎡ b1_1 ⋅ D c1_1 ⋅ D b1_2 ⋅ D c1_2 ⋅ D b1_3 ⋅ D c1_3 ⋅ D ⎤

0,0 0,1 0,0 0,1 0,0 0,1
1 ⎢ ⎥
S1 ≔ ――― ⋅ ⎢ b1_1 ⋅ D c1_1 ⋅ D b1_2 ⋅ D c1_2 ⋅ D b1_3 ⋅ D c1_3 ⋅ D ⎥
2 ⋅ Ae1 ⎢ 1,0 1,1 1,0 1,1 1,0 1,1
⎥
⎢ c1_1 ⋅ D b1_1 ⋅ D c1_2 ⋅ D b1_2 ⋅ D c1_3 ⋅ D b1_3 ⋅ D ⎥
⎣ 2,2 2,2 2,2 2,2 2,2 2,2⎦

5
 2 ⋅ Ae1 ⎢ ⎥
⎢ c1_1 ⋅ D b1_1 ⋅ D c1_2 ⋅ D b1_2 ⋅ D c1_3 ⋅ D b1_3 ⋅ D ⎥
⎣ 2,2 2,2 Zadatak
2,2 2 2,2 2,2 2,2⎦

⎡ -5555555.6 -694444.4 5555555.6 -694444.4 0.0 1388888.9 ⎤

⎢
S1 = -1388888.9 -2777777.8 1388888.9 -2777777.8 0.0 5555555.6 ⎥
⎢ ⎥
⎣ -1041666.7 -2083333.3 -1041666.7 2083333.3 2083333.3 0.0 ⎦

⎡ -972.222 ⎤ ⎡ -972.222 ⎤
Komponente napona: σ1 ≔ S1 ⋅ d1 = ⎢ -3888.889 ⎥ ili σ1 ≔ D ⋅ ε1 = ⎢ -3888.889 ⎥
⎢ ⎥ ⎢ ⎥
⎣ 2638.889 ⎦ ⎣ 2638.889 ⎦

σ1z ≔ ν ⋅ ⎛σ1 + σ1 ⎞ = -972.222

⎝ 0 1⎠

Raspodela kinematičkih i statičkih veličina - KE 2

⎡ 0 ⎤
⎢ 0 ⎥ ⎡ CST KE 2
⎤ ⎡1⎤
⎢ ⎥ ⎢ 0.0000000000 ⎥ ⎢2⎥ čvorovi - mreža: 1, 3, 4
da
⎢ 0 ⎥ ⎢ 0.0000000000 ⎥ ⎢ ⎥ čvorovi - KE: 1, 2, 3
⎢ ⎥ 0.0012666667 ⎥ ⎢5⎥ SS GKS: 1, 2, 5, 6, 7, 8
d2 ≔ ⎢ da ⎥ = ⎢
1
⎢ -0.0007000000 ⎥ ⎢6⎥ SS KE: 1, 2, 3, 4, 5, 6
⎢ ⎥ ⎢ ⎥ ⎢7⎥
⎢ 2 ⎥ ⎢ 0.0014666667 ⎥
da
⎢⎣ 8 ⎥⎦
⎢ ⎥ ⎣ -0.0003000000 ⎦
da
⎢⎣ 3 ⎥⎦
Raspodela pomeranja
1 float , 5
čvor 3 N2_2 ((x , y)) ≔ ――― ⋅ ((a2_2 + b2_2 ⋅ x + c2_2 ⋅ y)) ――― → 0.33333 ⋅ x
2 ⋅ Ae2

1 float , 5
čvor 4 N2_3 ((x , y)) ≔ ――― ⋅ ((a2_3 + b2_3 ⋅ x + c2_3 ⋅ y)) ――― → 0.16667 ⋅ y - 0.33333 ⋅ x
2 ⋅ Ae2

float , 5
u2 ((x , y)) ≔ N2_2 ((x , y)) ⋅ d2 + N2_3 ((x , y)) ⋅ d2 ――― → 0.00024445 ⋅ y - 0.000066666 ⋅ x
2 4

float , 5
v2 ((x , y)) ≔ N2_2 ((x , y)) ⋅ d2 + N2_3 ((x , y)) ⋅ d2 ――― → -0.000050001 ⋅ y - 0.00013333 ⋅ x
3 5

⎡ -0.0000666667 ⎤
Komponente deformacije: ε2 ≔ B2 ⋅ d2 = ⎢ -0.0000500000 ⎥
⎢ ⎥
⎣ 0.0001111111 ⎦
⎡ b2_1 ⋅ D c2_1 ⋅ D b2_2 ⋅ D c2_2 ⋅ D b2_3 ⋅ D c2_3 ⋅ D ⎤
0,0 0,1 0,0 0,1 0,0 0,1
1 ⎢ ⎥
S2 ≔ ――― ⋅ ⎢ b2_1 ⋅ D c2_1 ⋅ D b2_2 ⋅ D c2_2 ⋅ D b2_3 ⋅ D c2_3 ⋅ D ⎥
2 ⋅ Ae2 ⎢ 1,0 1,1 1,0 1,1 1,0 1,1
⎥
⎢ c2_1 ⋅ D b2_1 ⋅ D c2_2 ⋅ D b2_2 ⋅ D c2_3 ⋅ D b2_3 ⋅ D ⎥
⎣ 2,2 2,2 2,2 2,2 2,2 2,2⎦

⎡ 0.0 -1388888.9 11111111.1 0.0 -11111111.1 1388888.9 ⎤

S2 = ⎢ 0.0 -5555555.6 2777777.8 0.0 -2777777.8 5555555.6 ⎥
⎢ ⎥
⎣ -2083333.3 0.0 0.0 4166666.7 2083333.3 -4166666.7 ⎦

⎡ -2638.889 ⎤ ⎡ -2638.889 ⎤
Komponente napona: σ2 ≔ S2 ⋅ d2 = ⎢ -2222.222 ⎥ ili σ2 ≔ D ⋅ ε2 = ⎢ -2222.222 ⎥
⎢ ⎥ ⎢ ⎥
⎣ 1388.889 ⎦ ⎣ 1388.889 ⎦

σ2z ≔ ν ⋅ ⎛σ2 + σ2 ⎞ = -972.222

⎝ 0 1⎠

6
 Zadatak 3

Odrediti komponente pomeranja u tački „m” sa prirodnim koordinatama ξm = 0 i ηm = 0, kao i

vrednosti komponenata deformacije i napona u tačkama numeričke integracije. Numeričku analizu
sprovesti sa jednim izoparametarskim četvorougaonim konačnim elementom sa četiri čvora
(ravansko stanje napona). Primeniti red numeričke integracije 2x2. Podaci: E = 30 GPa; ν = 0,2;
h = 0,15 m; L1 = 6,1 m; L2 = 5,2 m; L3 = 4,2 m; P1 = 1000 kN i P2 = 1000 kN.

E ≔ 30 ⋅ 10 6 ν ≔ 0.2 h ≔ 0.15 Nepoznata pomeranja: 5 i 6

x1 ≔ 0 y1 ≔ 0 x2 ≔ 6.1 y2 ≔ 0 x3 ≔ 6.1 y3 ≔ 4.2 x4 ≔ 0 y4 ≔ 5.2

čvorovi - mreža: 1, 2, 3, 4
wa=wb=wc=wd=1.0 w≔1 čvorovi - KE: 1, 2, 3, 4
SS GKS: 1, 2, 3, 4, 5, 6, 7, 8
SS KE: 1, 2, 3, 4, 5, 6, 7, 8

⎡ 1 ν 0 ⎤ ⎡ 31250000.000 6250000.000 0.000 ⎤

E ⎢ ⎥ ⎢
D ≔ ――― ⋅ ν 1
2 ⎢
0 = 6250000.000 31250000.000 0.000 ⎥
1-ν ⎥ ⎢ ⎥
1-ν ⎥ ⎣
⎢ 0 0 ―― 0.000 0.000 12500000.000 ⎦
⎢⎣ 2 ⎥⎦

Matrica krutosti

Tačka numeričke integracije a

1 1
ξ ≔ -―― η ≔ -――
‾‾
3 ‾‾
3
1
J11 ≔ ―⋅ ((x1 ⋅ ((-1 + η)) + x2 ⋅ ((1 - η)) + x3 ⋅ ((1 + η)) + x4 ⋅ ((-1 - η)))) = 3.05000
4

1
J12 ≔ ―⋅ ((y1 ⋅ ((-1 + η)) + y2 ⋅ ((1 - η)) + y3 ⋅ ((1 + η)) + y4 ⋅ ((-1 - η)))) = -0.10566
4

1
J21 ≔ ―⋅ ((x1 ⋅ ((-1 + ξ)) + x2 ⋅ ((-1 - ξ)) + x3 ⋅ ((1 + ξ)) + x4 ⋅ ((1 - ξ)))) = 0.00000
4

1
J22 ≔ ―⋅ ((y1 ⋅ ((-1 + ξ)) + y2 ⋅ ((-1 - ξ)) + y3 ⋅ ((1 + ξ)) + y4 ⋅ ((1 - ξ)))) = 2.49434
4

detJ ≔ J11 ⋅ J22 - J12 ⋅ J21 = 7.60773

1
 Zadatak 3

1
dN1dx ≔ ――― ⋅ ((y2 ⋅ ((1 - η)) + y3 ⋅ ((-ξ + η)) + y4 ⋅ ((-1 + ξ)))) = -0.13477
8 ⋅ detJ

1
dN1dy ≔ ――― ⋅ ((x2 ⋅ ((-1 + η)) + x3 ⋅ ((ξ - η)) + x4 ⋅ ((1 - ξ)))) = -0.15809
8 ⋅ detJ

1
dN2dx ≔ ――― ⋅ ((y1 ⋅ ((-1 + η)) + y3 ⋅ ((1 + ξ)) + y4 ⋅ ((-ξ - η)))) = 0.12782
8 ⋅ detJ

1
dN2dy ≔ ――― ⋅ ((x1 ⋅ ((1 - η)) + x3 ⋅ ((-1 - ξ)) + x4 ⋅ ((ξ + η)))) = -0.04236
8 ⋅ detJ

1
dN3dx ≔ ――― ⋅ ((y1 ⋅ ((ξ - η)) + y2 ⋅ ((-1 - ξ)) + y4 ⋅ ((1 + η)))) = 0.03611
8 ⋅ detJ

1
dN3dy ≔ ――― ⋅ ((x1 ⋅ ((-ξ + η)) + x2 ⋅ ((1 + ξ)) + x4 ⋅ ((-1 - η)))) = 0.04236
8 ⋅ detJ

1
dN4dx ≔ ――― ⋅ ((y1 ⋅ ((1 - ξ)) + y2 ⋅ ((ξ + η)) + y3 ⋅ ((-1 - η)))) = -0.02917
8 ⋅ detJ

1
dN4dy ≔ ――― ⋅ ((x1 ⋅ ((-1 + ξ)) + x2 ⋅ ((-ξ - η)) + x3 ⋅ ((1 + η)))) = 0.15809
8 ⋅ detJ

⎡ dN1dx 0 dN2dx 0 dN3dx 0 dN4dx 0 ⎤

Ba ≔ ⎢ 0 dN1dy 0 dN2dy 0 dN3dy 0 dN4dy ⎥
⎢ ⎥
⎣ dN1dy dN1dx dN2dy dN2dx dN3dy dN3dx dN4dy dN4dx ⎦

⎡ -0.13477 0.00000 0.12782 0.00000 0.03611 0.00000 -0.02917 0.00000 ⎤

Ba = ⎢ 0.00000 -0.15809 0.00000 -0.04236 0.00000 0.04236 0.00000 0.15809 ⎥
⎢ ⎥
⎣ -0.15809 -0.13477 -0.04236 0.12782 0.04236 0.03611 0.15809 -0.02917 ⎦

ka ≔ Ba T ⋅ D ⋅ Ba ⋅ h ⋅ detJ ⋅ w

1 2 3 4 5 6 7 8
⎡ 10.04212 4.55876 -5.18789 -2.47540 -2.69078 -1.22152 -2.16345 -0.86185 ⎤ ⎡ 1 ⎤
⎢ 4.55876 11.50374 -0.62694 -0.06905 -1.22152 -3.08242 -2.71031 -8.35227 ⎥ ⎢ 2 ⎥
⎢ ⎥⎢ ⎥
⎢ -5.18789 -0.62694 6.08260 -1.15857 1.39009 0.16799 -2.28480 1.61753 ⎥ ⎢ 3 ⎥
ka ⎢ -2.47540 -0.06905 -1.15857 2.97058 0.66328 0.01850 2.97069 -2.92003 ⎥ ⎢ 4 ⎥
―― =
10 5 ⎢ -2.69078 -1.22152 1.39009 0.66328 0.72099 0.32730 0.57969 0.23093 ⎥ ⎢ 5 ⎥
⎢ -1.22152 -3.08242 0.16799 0.01850 0.32730 0.82593 0.72622 2.23798 ⎥ ⎢ 6 ⎥
⎢ ⎥⎢ ⎥
-2.16345 -2.71031 -2.28480 2.97069 0.57969 0.72622 3.86855 -0.98661 7
⎢ ⎥⎢ ⎥
⎣ -0.86185 -8.35227 1.61753 -2.92003 0.23093 2.23798 -0.98661 9.03431 ⎦ ⎣ 8 ⎦

5 6
⎡ ka ka ⎤
ka_aa ≔ ⎢
4,4 4,5
⎥ = ⎡ 72099.15012 32730.44140 ⎤ ⎡ 5 ⎤
⎢ ka ka ⎥ ⎢⎣ 32730.44140 82593.11573 ⎥⎦ ⎢⎣ 6 ⎥⎦
⎣ 5,4 5,5⎦

2
 Zadatak 3

Tačka numeričke integracije b

1 1
ξ ≔ ―― η ≔ -――
‾‾
3 ‾‾
3
1
J11 ≔ ―⋅ ((x1 ⋅ ((-1 + η)) + x2 ⋅ ((1 - η)) + x3 ⋅ ((1 + η)) + x4 ⋅ ((-1 - η)))) = 3.05000
4

1
J12 ≔ ―⋅ ((y1 ⋅ ((-1 + η)) + y2 ⋅ ((1 - η)) + y3 ⋅ ((1 + η)) + y4 ⋅ ((-1 - η)))) = -0.10566
4

1
J21 ≔ ―⋅ ((x1 ⋅ ((-1 + ξ)) + x2 ⋅ ((-1 - ξ)) + x3 ⋅ ((1 + ξ)) + x4 ⋅ ((1 - ξ)))) = 0.00000
4

1
J22 ≔ ―⋅ ((y1 ⋅ ((-1 + ξ)) + y2 ⋅ ((-1 - ξ)) + y3 ⋅ ((1 + ξ)) + y4 ⋅ ((1 - ξ)))) = 2.20566
4
detJ ≔ J11 ⋅ J22 - J12 ⋅ J21 = 6.72727
1
dN1dx ≔ ――― ⋅ ((y2 ⋅ ((1 - η)) + y3 ⋅ ((-ξ + η)) + y4 ⋅ ((-1 + ξ)))) = -0.13095
8 ⋅ detJ

1
dN1dy ≔ ――― ⋅ ((x2 ⋅ ((-1 + η)) + x3 ⋅ ((ξ - η)) + x4 ⋅ ((1 - ξ)))) = -0.04791
8 ⋅ detJ

1
dN2dx ≔ ――― ⋅ ((y1 ⋅ ((-1 + η)) + y3 ⋅ ((1 + ξ)) + y4 ⋅ ((-ξ - η)))) = 0.1231
8 ⋅ detJ

1
dN2dy ≔ ――― ⋅ ((x1 ⋅ ((1 - η)) + x3 ⋅ ((-1 - ξ)) + x4 ⋅ ((ξ + η)))) = -0.17878
8 ⋅ detJ

1
dN3dx ≔ ――― ⋅ ((y1 ⋅ ((ξ - η)) + y2 ⋅ ((-1 - ξ)) + y4 ⋅ ((1 + η)))) = 0.04084
8 ⋅ detJ

1
dN3dy ≔ ――― ⋅ ((x1 ⋅ ((-ξ + η)) + x2 ⋅ ((1 + ξ)) + x4 ⋅ ((-1 - η)))) = 0.17878
8 ⋅ detJ

1
dN4dx ≔ ――― ⋅ ((y1 ⋅ ((1 - ξ)) + y2 ⋅ ((ξ + η)) + y3 ⋅ ((-1 - η)))) = -0.03298
8 ⋅ detJ

1
dN4dy ≔ ――― ⋅ ((x1 ⋅ ((-1 + ξ)) + x2 ⋅ ((-ξ - η)) + x3 ⋅ ((1 + η)))) = 0.04791
8 ⋅ detJ

⎡ dN1dx 0 dN2dx 0 dN3dx 0 dN4dx 0 ⎤

Bb ≔ ⎢ 0 dN1dy 0 dN2dy 0 dN3dy 0 dN4dy ⎥
⎢ ⎥
⎣ dN1dy dN1dx dN2dy dN2dx dN3dy dN3dx dN4dy dN4dx ⎦

⎡ -0.13095 0.00000 0.12310 0.00000 0.04084 0.00000 -0.03298 0.00000 ⎤

Bb = ⎢ 0.00000 -0.04791 0.00000 -0.17878 0.00000 0.17878 0.00000 0.04791 ⎥
⎢ ⎥
⎣ -0.04791 -0.13095 -0.17878 0.12310 0.17878 0.04084 0.04791 -0.03298 ⎦

kb ≔ Bb T ⋅ D ⋅ Bb ⋅ h ⋅ detJ ⋅ w

3
 Zadatak 3

1 2 3 4 5 6 7 8
⎡ 5.69695 1.18692 -4.00287 0.73272 -2.76665 -1.72331 1.07257 -0.19633 ⎤ ⎡ 1 ⎤
⎢ 1.18692 2.88667 2.58118 0.66751 -3.07647 -3.37532 -0.69162 -0.17886 ⎥ ⎢ 2 ⎥
⎢ ⎥⎢ ⎥
⎢ -4.00287 2.58118 8.81014 -4.16398 -2.44660 0.46707 -2.36067 1.11574 ⎥ ⎢ 3 ⎥
kb ⎢ 0.73272 0.66751 -4.16398 11.99082 2.31553 -9.44541 1.11574 -3.21293 ⎥ ⎢ 4 ⎥
―― =
10 5 ⎢ -2.76665 -3.07647 -2.44660 2.31553 4.55768 1.38139 0.65556 -0.62044 ⎥ ⎢ 5 ⎥
⎢ -1.72331 -3.37532 0.46707 -9.44541 1.38139 10.28984 -0.12515 2.53089 ⎥ ⎢ 6 ⎥
⎢ ⎥⎢ ⎥
1.07257 -0.69162 -2.36067 1.11574 0.65556 -0.12515 0.63254 -0.29896 7
⎢ ⎥⎢ ⎥
⎣ -0.19633 -0.17886 1.11574 -3.21293 -0.62044 2.53089 -0.29896 0.86090 ⎦ ⎣ 8 ⎦

5 6
⎡ kb kb ⎤
kb_aa ≔ ⎢ 4,4 4,5
⎥ = ⎡ 455767.87581 138138.77205 ⎤ ⎡ 5 ⎤
⎢ kb kb ⎥ ⎢⎣ 138138.77205 1028983.96199 ⎥⎦ ⎢⎣ 6 ⎥⎦
⎣ 5,4 5,5⎦

Tačka numeričke integracije c

1 1
ξ ≔ ―― η ≔ ――
‾‾
3 ‾‾
3
1
J11 ≔ ―⋅ ((x1 ⋅ ((-1 + η)) + x2 ⋅ ((1 - η)) + x3 ⋅ ((1 + η)) + x4 ⋅ ((-1 - η)))) = 3.05000
4

1
J12 ≔ ―⋅ ((y1 ⋅ ((-1 + η)) + y2 ⋅ ((1 - η)) + y3 ⋅ ((1 + η)) + y4 ⋅ ((-1 - η)))) = -0.39434
4

1
J21 ≔ ―⋅ ((x1 ⋅ ((-1 + ξ)) + x2 ⋅ ((-1 - ξ)) + x3 ⋅ ((1 + ξ)) + x4 ⋅ ((1 - ξ)))) = 0.00000
4

1
J22 ≔ ―⋅ ((y1 ⋅ ((-1 + ξ)) + y2 ⋅ ((-1 - ξ)) + y3 ⋅ ((1 + ξ)) + y4 ⋅ ((1 - ξ)))) = 2.20566
4
detJ ≔ J11 ⋅ J22 - J12 ⋅ J21 = 6.72727
1
dN1dx ≔ ――― ⋅ ((y2 ⋅ ((1 - η)) + y3 ⋅ ((-ξ + η)) + y4 ⋅ ((-1 + ξ)))) = -0.04084
8 ⋅ detJ

1
dN1dy ≔ ――― ⋅ ((x2 ⋅ ((-1 + η)) + x3 ⋅ ((ξ - η)) + x4 ⋅ ((1 - ξ)))) = -0.04791
8 ⋅ detJ

1
dN2dx ≔ ――― ⋅ ((y1 ⋅ ((-1 + η)) + y3 ⋅ ((1 + ξ)) + y4 ⋅ ((-ξ - η)))) = 0.01153
8 ⋅ detJ

1
dN2dy ≔ ――― ⋅ ((x1 ⋅ ((1 - η)) + x3 ⋅ ((-1 - ξ)) + x4 ⋅ ((ξ + η)))) = -0.17878
8 ⋅ detJ

1
dN3dx ≔ ――― ⋅ ((y1 ⋅ ((ξ - η)) + y2 ⋅ ((-1 - ξ)) + y4 ⋅ ((1 + η)))) = 0.15241
8 ⋅ detJ

1
dN3dy ≔ ――― ⋅ ((x1 ⋅ ((-ξ + η)) + x2 ⋅ ((1 + ξ)) + x4 ⋅ ((-1 - η)))) = 0.17878
8 ⋅ detJ

4
 Zadatak 3

1
dN4dx ≔ ――― ⋅ ((y1 ⋅ ((1 - ξ)) + y2 ⋅ ((ξ + η)) + y3 ⋅ ((-1 - η)))) = -0.1231
8 ⋅ detJ

1
dN4dy ≔ ――― ⋅ ((x1 ⋅ ((-1 + ξ)) + x2 ⋅ ((-ξ - η)) + x3 ⋅ ((1 + η)))) = 0.04791
8 ⋅ detJ

⎡ dN1dx 0 dN2dx 0 dN3dx 0 dN4dx 0 ⎤

Bc ≔ ⎢ 0 dN1dy 0 dN2dy 0 dN3dy 0 dN4dy ⎥
⎢ ⎥
⎣ dN1dy dN1dx dN2dy dN2dx dN3dy dN3dx dN4dy dN4dx ⎦

⎡ -0.04084 0.00000 0.01153 0.00000 0.15241 0.00000 -0.12310 0.00000 ⎤

Bc = ⎢ 0.00000 -0.04791 0.00000 -0.17878 0.00000 0.17878 0.00000 0.04791 ⎥
⎢ ⎥
⎣ -0.04791 -0.04084 -0.17878 0.01153 0.17878 0.15241 0.04791 -0.12310 ⎦

kc ≔ Bc T ⋅ D ⋅ Bc ⋅ h ⋅ detJ ⋅ w

1 2 3 4 5 6 7 8
⎡ 0.81535 0.37014 0.93186 0.39080 -3.04294 -1.38139 1.29573 0.62044 ⎤ ⎡ 1 ⎤
⎢ 0.37014 0.93403 0.88610 2.64141 -1.38139 -3.48584 0.12515 -0.08959 ⎥ ⎢ 2 ⎥
⎢ ⎥⎢ ⎥
⎢ 0.93186 0.88610 4.07370 -0.38996 -3.47775 -3.30695 -1.52781 2.81082 ⎥ ⎢ 3 ⎥
kc ⎢ 0.39080 2.64141 -0.38996 10.09625 -1.45849 -9.85787 1.45765 -2.87979 ⎥ ⎢ 4 ⎥
―― =
10 5 ⎢ -3.04294 -1.38139 -3.47775 -1.45849 11.35642 5.15541 -4.83573 -2.31553 ⎥ ⎢ 5 ⎥
⎢ -1.38139 -3.48584 -3.30695 -9.85787 5.15541 13.00934 -0.46707 0.33437 ⎥ ⎢ 6 ⎥
⎢ ⎥⎢ ⎥
1.29573 0.12515 -1.52781 1.45765 -4.83573 -0.46707 5.06781 -1.11574 7
⎢ ⎥⎢ ⎥
⎣ 0.62044 -0.08959 2.81082 -2.87979 -2.31553 0.33437 -1.11574 2.63501 ⎦ ⎣ 8 ⎦

5 6
⎡ kc kc ⎤
kc_aa ≔ ⎢
4,4 4,5
⎥ = ⎡ 1135641.93740 515540.91577 ⎤ ⎡5⎤
⎢ kc kc ⎥ ⎢⎣ 515540.91577 1300933.58663 ⎥⎦ ⎢⎣ 6 ⎥⎦
⎣ 5,4 5,5⎦

Tačka numeričke integracije d

1 1
ξ ≔ -―― η ≔ ――
‾‾
3 ‾‾
3
1
J11 ≔ ―⋅ ((x1 ⋅ ((-1 + η)) + x2 ⋅ ((1 - η)) + x3 ⋅ ((1 + η)) + x4 ⋅ ((-1 - η)))) = 3.05000
4

1
J12 ≔ ―⋅ ((y1 ⋅ ((-1 + η)) + y2 ⋅ ((1 - η)) + y3 ⋅ ((1 + η)) + y4 ⋅ ((-1 - η)))) = -0.39434
4

1
J21 ≔ ―⋅ ((x1 ⋅ ((-1 + ξ)) + x2 ⋅ ((-1 - ξ)) + x3 ⋅ ((1 + ξ)) + x4 ⋅ ((1 - ξ)))) = 0.00000
4

1
J22 ≔ ―⋅ ((y1 ⋅ ((-1 + ξ)) + y2 ⋅ ((-1 - ξ)) + y3 ⋅ ((1 + ξ)) + y4 ⋅ ((1 - ξ)))) = 2.49434
4
detJ ≔ J11 ⋅ J22 - J12 ⋅ J21 = 7.60773
1
dN1dx ≔ ――― ⋅ ((y2 ⋅ ((1 - η)) + y3 ⋅ ((-ξ + η)) + y4 ⋅ ((-1 + ξ)))) = -0.05508
8 ⋅ detJ

5
 Zadatak 3

1
dN1dy ≔ ――― ⋅ ((x2 ⋅ ((-1 + η)) + x3 ⋅ ((ξ - η)) + x4 ⋅ ((1 - ξ)))) = -0.15809
8 ⋅ detJ

1
dN2dx ≔ ――― ⋅ ((y1 ⋅ ((-1 + η)) + y3 ⋅ ((1 + ξ)) + y4 ⋅ ((-ξ - η)))) = 0.02917
8 ⋅ detJ

1
dN2dy ≔ ――― ⋅ ((x1 ⋅ ((1 - η)) + x3 ⋅ ((-1 - ξ)) + x4 ⋅ ((ξ + η)))) = -0.04236
8 ⋅ detJ

1
dN3dx ≔ ――― ⋅ ((y1 ⋅ ((ξ - η)) + y2 ⋅ ((-1 - ξ)) + y4 ⋅ ((1 + η)))) = 0.13477
8 ⋅ detJ

1
dN3dy ≔ ――― ⋅ ((x1 ⋅ ((-ξ + η)) + x2 ⋅ ((1 + ξ)) + x4 ⋅ ((-1 - η)))) = 0.04236
8 ⋅ detJ

1
dN4dx ≔ ――― ⋅ ((y1 ⋅ ((1 - ξ)) + y2 ⋅ ((ξ + η)) + y3 ⋅ ((-1 - η)))) = -0.10885
8 ⋅ detJ

1
dN4dy ≔ ――― ⋅ ((x1 ⋅ ((-1 + ξ)) + x2 ⋅ ((-ξ - η)) + x3 ⋅ ((1 + η)))) = 0.15809
8 ⋅ detJ

⎡ dN1dx 0 dN2dx 0 dN3dx 0 dN4dx 0 ⎤

Bd ≔ ⎢ 0 dN1dy 0 dN2dy 0 dN3dy 0 dN4dy ⎥
⎢ ⎥
⎣ dN1dy dN1dx dN2dy dN2dx dN3dy dN3dx dN4dy dN4dx ⎦

⎡ -0.05508 0.00000 0.02917 0.00000 0.13477 0.00000 -0.10885 0.00000 ⎤

Bd = ⎢ 0.00000 -0.15809 0.00000 -0.04236 0.00000 0.04236 0.00000 0.15809 ⎥
⎢ ⎥
⎣ -0.15809 -0.05508 -0.04236 0.02917 0.04236 0.13477 0.15809 -0.10885 ⎦

kd ≔ Bd T ⋅ D ⋅ Bd ⋅ h ⋅ detJ ⋅ w

1 2 3 4 5 6 7 8
⎡ 4.64721 1.86329 0.38236 -0.49132 -3.60259 -3.20560 -1.42698 1.83362 ⎤ ⎡ 1 ⎤
⎢ 1.86329 9.34578 0.00398 2.15905 -1.85243 -3.44714 -0.01484 -8.05768 ⎥ ⎢ 2 ⎥
⎢ ⎥⎢ ⎥
⎢ 0.38236 0.00398 0.55933 -0.26436 1.14577 -0.72622 -2.08746 0.98661 ⎥ ⎢ 3 ⎥
kd ⎢ -0.49132 2.15905 -0.26436 0.76127 -0.23093 -0.07923 0.98661 -2.84109 ⎥ ⎢ 4 ⎥
―― =
10 5 ⎢ -3.60259 -1.85243 1.14577 -0.23093 6.73290 1.22152 -4.27608 0.86185 ⎥ ⎢ 5 ⎥
⎢ -3.20560 -3.44714 -0.72622 -0.07923 1.22152 3.23069 2.71031 0.29568 ⎥ ⎢ 6 ⎥
⎢ ⎥⎢ ⎥
-1.42698 -0.01484 -2.08746 0.98661 -4.27608 2.71031 7.79052 -3.68208 7
⎢ ⎥⎢ ⎥
⎣ 1.83362 -8.05768 0.98661 -2.84109 0.86185 0.29568 -3.68208 10.60310 ⎦ ⎣ 8 ⎦

5 6
⎡ kd kd ⎤
kd_aa ≔ ⎢
4,4 4,5
⎥ = ⎡ 673289.90164 122151.67024 ⎤ ⎡ 5 ⎤
⎢ kd kd ⎥ ⎢⎣ 122151.67024 323069.41634 ⎥⎦ ⎢⎣ 6 ⎥⎦
⎣ 5 , 4 5,5⎦

6
 Zadatak 3

Matrica krutosti KE

k ≔ ka + kb + kc + kd

1 2 3 4 5 6 7 8
⎡ 21.202 7.979 -7.877 -1.843 -12.103 -7.532 -1.222 1.396 ⎤ ⎡ 1 ⎤
⎢ 7.979 24.670 2.844 5.399 -7.532 -13.391 -3.292 -16.678 ⎥ ⎢ 2 ⎥
⎢ ⎥ ⎢ ⎥
⎢ -7.877 2.844 19.526 -5.977 -3.388 -3.398 -8.261 6.531 ⎥ ⎢ 3 ⎥
k -1.843 5.399 -5.977 25.819 1.289 -19.364 6.531 -11.854 ⎥ ⎢ 4 ⎥
―― =⎢
10 5 ⎢ -12.103 -7.532 -3.388 1.289 23.368 8.086 -7.877 -1.843 ⎥ ⎢ 5 ⎥
⎢ -7.532 -13.391 -3.398 -19.364 8.086 27.356 2.844 5.399 ⎥ ⎢ 6 ⎥
⎢ ⎥ ⎢ ⎥
-1.222 -3.292 -8.261 6.531 -7.877 2.844 17.359 -6.083 7
⎢ ⎥ ⎢ ⎥
⎣ 1.396 -16.678 6.531 -11.854 -1.843 5.399 -6.083 23.133 ⎦ ⎣ 8 ⎦
Sistem KE

Kaa ≔ ka_aa + kb_aa + kc_aa + kd_aa

5 6
Kaa ⎡ 23.36799 8.08562 ⎤ ⎡5⎤ ⎡ 1000 ⎤ ⎡ 5 ⎤
―― =⎢ ⎥ ⎢⎣ 6 ⎥⎦ Sa ≔ ⎢
10 5 ⎣ 8.08562 27.35580 ⎦ ⎣ 1000 ⎥⎦ ⎢⎣ 6 ⎥⎦

⎡ 0.0003357919 ⎤ ⎡ 5 ⎤
da ≔ Kaa -1 ⋅ Sa = ⎢
⎣ 0.0002663024 ⎥⎦ ⎢⎣ 6 ⎥⎦

Komponente pomeranja tačke m

⎡ 0 ⎤ ξm ≔ 0 ηm ≔ 0
⎢ 0 ⎥ ⎡ 0.0000000000 ⎤ ⎡1⎤
⎢ ⎥ ⎢ 0.0000000000 ⎥ ⎢2⎥
⎢ 0 ⎥ ⎢ 0.0000000000 ⎥ ⎢ ⎥ ξ1 ≔ -1 η1 ≔ -1
⎢ 0 ⎥ ⎢ ⎥ ⎢3⎥
⎢ 0.0000000000 ⎥ ⎢4⎥
dke ≔ ⎢ da ⎥ = ξ2 ≔ 1 η2 ≔ -1
⎢ 0 ⎥ ⎢ 0.0003357919 ⎥ ⎢5⎥
⎢ da ⎥ ⎢ 0.0002663024 ⎥ ⎢6⎥
⎢ 1 ⎥ ⎢ 0.0000000000 ⎥ ⎢ ⎥ ξ3 ≔ 1 η3 ≔ 1
7
⎢ 0 ⎥ ⎢ ⎥ ⎢ ⎥
⎣ 0.0000000000 ⎦ ⎣8⎦
⎢⎣ 0 ⎥⎦ ξ4 ≔ -1 η4 ≔ 1

1 1
N1 ((ξ , η)) ≔ ―⋅ ((1 + ξ1 ⋅ ξ)) ⋅ ((1 + η1 ⋅ η)) N2 ((ξ , η)) ≔ ―⋅ ((1 + ξ2 ⋅ ξ)) ⋅ ((1 + η2 ⋅ η))
4 4

1 1
N3 ((ξ , η)) ≔ ―⋅ ((1 + ξ3 ⋅ ξ)) ⋅ ((1 + η3 ⋅ η)) N4 ((ξ , η)) ≔ ―⋅ ((1 + ξ4 ⋅ ξ)) ⋅ ((1 + η4 ⋅ η))
4 4

u_m ≔ N3 ((ξm , ηm)) ⋅ dke = 0.0000839480 v_m ≔ N3 ((ξm , ηm)) ⋅ dke = 0.0000665756
4 5

Dekartove koordinate tačke m

xm ≔ N1 ((ξm , ηm)) ⋅ x1 + N2 ((ξm , ηm)) ⋅ x2 + N3 ((ξm , ηm)) ⋅ x3 + N4 ((ξm , ηm)) ⋅ x4 = 3.05

ym ≔ N1 ((ξm , ηm)) ⋅ y1 + N2 ((ξm , ηm)) ⋅ y2 + N3 ((ξm , ηm)) ⋅ y3 + N4 ((ξm , ηm)) ⋅ y4 = 2.35

7
 Zadatak 3

Komponente deformacije i napona u tački numeričke integracije a

⎡ 0.0000121258 ⎤ ⎡ 449.435 ⎤
εa ≔ Ba ⋅ dke = ⎢ 0.0000112808 ⎥ σa ≔ D ⋅ εa = ⎢ 428.311 ⎥
⎢ ⎥ ⎢ ⎥
⎣ 0.0000238409 ⎦ ⎣ 298.011 ⎦
-ν -ν
εaz ≔ ――⋅ ⎛εa + εa ⎞ = -0.0000058516 ili εaz ≔ ―― ⋅ ⎛σa + σa ⎞ = -0.0000058516
1-ν ⎝ 0 1⎠ E ⎝ 0 1⎠

Komponente deformacije i napona u tački numeričke integracije b

⎡ 0.0000137128 ⎤ ⎡ 726.091 ⎤
εb ≔ Bb ⋅ dke = ⎢ 0.0000476107 ⎥ σb ≔ D ⋅ εb = ⎢ 1573.538 ⎥
⎢ ⎥ ⎢ ⎥
⎣ 0.0000709093 ⎦ ⎣ 886.366 ⎦
-ν -ν
εbz ≔ ――⋅ ⎛εb + εb ⎞ = -0.0000153309 ili εbz ≔ ―― ⋅ ⎛σb + σb ⎞ = -0.0000153309
1-ν ⎝ 0 1⎠ E ⎝ 0 1⎠

Komponente deformacije i napona u tački numeričke integracije c

⎡ 0.0000511768 ⎤ ⎡ 1896.840 ⎤
εc ≔ Bc ⋅ dke = ⎢ 0.0000476107 ⎥ σc ≔ D ⋅ εc = ⎢ 1807.688 ⎥
⎢ ⎥ ⎢ ⎥
⎣ 0.0001006204 ⎦ ⎣ 1257.755 ⎦
-ν -ν
εcz ≔ ――⋅ ⎛εc + εc ⎞ = -0.0000246969 ili εcz ≔ ―― ⋅ ⎛σc + σc ⎞ = -0.0000246969
1-ν ⎝ 0 1⎠ E ⎝ 0 1⎠

Komponente deformacije i napona u tački numeričke integracije d

⎡ 0.0000452540 ⎤ ⎡ 1484.691 ⎤
εd ≔ Bd ⋅ dke = ⎢ 0.0000112808 ⎥ σd ≔ D ⋅ εd = ⎢ 635.363 ⎥
⎢ ⎥ ⎢ ⎥
⎣ 0.0000501135 ⎦ ⎣ 626.418 ⎦
-ν -ν
εdz ≔ ――⋅ ⎛εd + εd ⎞ = -0.0000141337 ili εdz ≔ ―― ⋅ ⎛σd + σd ⎞ = -0.0000141337
1-ν ⎝ 0 1⎠ E ⎝ 0 1⎠

8
 Zadatak 4

Odrediti komponente pomeranja, deformacije i napona u tački "c". Za numeričku simulaciju

formirati diskretni model sa jednim pravougaonim konačnim elementom koji ima četiri čvora
(ravansko stanje napona). Podaci: E = 30 GPa; ν = 0,2; h = 0,15 m; L1 = 6,0 m; L2 = 4,0 m; P1 = 500 kN i
P2 = 300 kN.

E ≔ 30 ⋅ 10 6

ν ≔ 0.2

h ≔ 0.15

L1 ≔ 6

L2 ≔ 4

L1 L2
čvorovi - mreža: 1, 2, 3, 4 a ≔ ―― =3 b ≔ ―― =2
čvorovi - KE: 1, 2, 3, 4 2 2
SS GKS: 1, 2, 3, 4, 5, 6, 7, 8 Nepoznata pomeranja: 3, 4, 5 i 6
SS KE: 1, 2, 3, 4, 5, 6, 7, 8

⎡ 1 ν 0 ⎤ ⎡ 31250000.000 6250000.000 0.000 ⎤

E ⎢ν 1 ⎥ = ⎢ 6250000.000 31250000.000
D ≔ ――― ⋅
2 ⎢
0 0.000 ⎥
1-ν ⎥ ⎢ ⎥
1-ν ⎥ ⎣
⎢ 0 0 ―― 0.000 0.000 12500000.000 ⎦
⎢⎣ 2 ⎥⎦

Napomena:
određeni su samo elementi matrice krutosti KE uz nepoznata pomeranja

h
k33 ≔ ――― ⋅ ⎛4 ⋅ ⎛b 2 ⋅ D + a 2 ⋅ D ⎞⎞ = 1979166.7
12 ⋅ a ⋅ b ⎜⎝ ⎜⎝ 0,0 2 , 2⎟
⎠⎟⎠

h
k34 ≔ ――― ⋅ ⎛-3 ⋅ a ⋅ b ⋅ ⎛D + D ⎞⎞ = -703125.0 k43 ≔ k34
12 ⋅ a ⋅ b ⎝ ⎝ 0,1 2 , 2⎠⎠

h
k35 ≔ ――― ⋅ ⎛2 ⋅ b 2 ⋅ D - 4 ⋅ a 2 ⋅ D ⎞ = -416666.7 k53 ≔ k35
12 ⋅ a ⋅ b ⎜⎝ 0,0 2 , 2⎟
⎠

h
k36 ≔ ――― ⋅ ⎛3 ⋅ a ⋅ b ⋅ ⎛D - D ⎞⎞ = -234375.0 k63 ≔ k36
12 ⋅ a ⋅ b ⎝ ⎝ 0,1 2 , 2⎠⎠

h
k44 ≔ ――― ⋅ ⎛4 ⋅ ⎛a 2 ⋅ D + b 2 ⋅ D ⎞⎞ = 2760416.7
12 ⋅ a ⋅ b ⎜⎝ ⎜⎝ 1,1 2 , 2⎟
⎠⎟⎠

h
k45 ≔ ――― ⋅ ⎛3 ⋅ a ⋅ b ⋅ ⎛-D + D ⎞⎞ = 234375.0 k54 ≔ k45
12 ⋅ a ⋅ b ⎝ ⎝ 1,0 2 , 2⎠⎠

h
k46 ≔ ――― ⋅ ⎛-4 ⋅ a 2 ⋅ D + 2 ⋅ b 2 ⋅ D ⎞ = -2135416.7 k64 ≔ k46
12 ⋅ a ⋅ b ⎜⎝ 1,1 2 , 2⎟
⎠
k55 ≔ k33 = 1979166.7 k56 ≔ -k34 = 703125.0 k65 ≔ k56 = 703125.0

k66 ≔ k44 = 2760416.7

3 4 5 6
1
 Zadatak 4

3 4 5 6
⎡ k33 k34 k35 k36 ⎤ ⎡ 1979166.7 -703125.0 -416666.7 -234375.0 ⎤ ⎡ 3 ⎤
⎢ k43 k44 k45 k46 ⎥ ⎢ -703125.0 2760416.7 234375.0 -2135416.7 ⎥ ⎢ 4 ⎥
k_aa ≔ ⎢ ⎥=⎢ ⎥ ⎢ ⎥
⎢ k53 k54 k55 k56 ⎥ ⎢ -416666.7 234375.0 1979166.7 703125.0 ⎥ ⎢ 5 ⎥
⎣ k63 k64 k65 k66 ⎦ ⎣ -234375.0 -2135416.7 703125.0 2760416.7 ⎦ ⎣ 6 ⎦
3 4 5 6
⎡ 1979166.7 -703125.0 -416666.7 -234375.0 ⎤ ⎡ 3 ⎤
⎢ -703125.0 2760416.7 234375.0 -2135416.7 ⎥ ⎢ 4 ⎥
Kaa ≔ k_aa = ⎢ ⎥⎢ ⎥
⎢ -416666.7 234375.0 1979166.7 703125.0 ⎥ ⎢ 5 ⎥
⎣ -234375.0 -2135416.7 703125.0 2760416.7 ⎦ ⎣ 6 ⎦

⎡ 0.0 ⎤ ⎡3⎤ ⎡ -0.0004474595 ⎤ ⎡3⎤

⎢ -500.0 ⎥ ⎢4⎥ ⎢ -0.0012389001 ⎥ ⎢4⎥
Sa ≔ ⎢ ⎥ ⎢ ⎥ da ≔ Kaa -1 ⋅ Sa = ⎢ ⎥ ⎢ ⎥
⎢ 300.0 ⎥ ⎢5⎥ ⎢ 0.0006135931 ⎥ ⎢5⎥
⎣ 0 ⎦ ⎣6⎦ ⎣ -0.0011526788 ⎦ ⎣6⎦

⎡ 0 ⎤
⎢ 0 ⎥ ⎡ ⎤ ⎡1⎤
⎢ ⎥ ⎢ 0.0000000000 ⎥ ⎢2⎥
⎢ da0 ⎥ ⎢ 0.0000000000 ⎥ 1 ⎛ x⎞ ⎛
N2 ((x , y)) ≔ ―⋅ ⎜1 + ―
y⎞
⎢ ⎥ ⎢ -0.0004474595 ⎥
⎢ ⎥ ⎟ ⋅ ⎜1 - ―⎟
⎢3⎥ 4 ⎝ a⎠ ⎝ b⎠
⎢ da1 ⎥ ⎢ -0.0012389001 ⎥ ⎢4⎥
dke ≔ ⎢ ⎥=⎢ ⎥ ⎢5⎥ 1 ⎛ x⎞ ⎛ y⎞
⎢ da2 ⎥ ⎢ 0.0006135931 ⎥ ⎢6⎥ N3 ((x , y)) ≔ ―⋅ ⎜1 + ―⎟ ⋅ ⎜1 + ―⎟
⎢ ⎥ ⎢ -0.0011526788 4 ⎝ a⎠ ⎝ b⎠
da ⎥ ⎢ ⎥
⎢ 3 ⎥ ⎢ 0.0000000000 ⎥ 7
⎢ ⎥
⎢ 0 ⎥ ⎣ 0.0000000000 ⎦ ⎣8⎦
⎢ ⎥
⎣ 0 ⎦

xc ≔ 0 yc ≔ 0

u ≔ N2 ((xc , yc)) ⋅ dke + N3 ((xc , yc)) ⋅ dke = 0.0000415334

2 4

v ≔ N2 ((xc , yc)) ⋅ dke + N3 ((xc , yc)) ⋅ dke = -0.0005978947

3 5

⎡ -b + yc 0 b - yc 0 b + yc 0 -b - yc 0 ⎤
1 ⎢
B ≔ ――― ⋅ 0 -a + xc 0 -a - xc 0 a + xc 0 a - xc ⎥
4 ⋅ a ⋅ b ⎢ -a - xc -b + yc -a - xc b - yc a + xc b + yc a - xc -b - yc ⎥
⎣ ⎦

⎡ -0.08333 0.00000 0.08333 0.00000 0.08333 0.00000 -0.08333 0.00000 ⎤

B = ⎢ 0.00000 -0.12500 0.00000 -0.12500 0.00000 0.12500 0.00000 0.12500 ⎥
⎢ ⎥
⎣ -0.12500 -0.08333 -0.12500 0.08333 0.12500 0.08333 0.12500 -0.08333 ⎦

⎡ 0.0000138445 ⎤ ⎡ 500.000 ⎤ ⎡ 75.000 ⎤

ε ≔ B ⋅ dke = ⎢ 0.0000107777 ⎥ σ ≔ D ⋅ ε = ⎢ 423.330 ⎥ sile ≔ σ ⋅ h = ⎢ 63.499 ⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎣ -0.0000666667 ⎦ ⎣ -833.333 ⎦ ⎣ -125.000 ⎦

-ν -ν
εz ≔ ――⋅ ⎛ε + ε ⎞ = -0.0000061555 ili εz ≔ ―― ⋅ ⎛σ + σ ⎞ = -0.0000061555
1-ν ⎝ 0 1⎠ E ⎝ 0 1⎠

2
 Zadatak 5

Odrediti raspodelu kinematičkih i statičkih veličina u konačnim elementima 1 i 2. Zid (KE 1)

diskretizovati sa jednim trougaonim konačnim elementom sa tri čvora (CST konačni element;
ravansko stanje napona).
Podaci:
Ÿ L1 = 6,0 m; L2 = 6,0 m; q = 10,0 MPa.
Ÿ KE 1: E = 30 GPa; ν = 0,2; h = 10,0 cm.
Ÿ KE 2: E = 200 GPa; A = 30,0 cm2.

x1 ≔ 0 y1 ≔ 0 x2 ≔ 6 y2 ≔ 0 x3 ≔ 0 y3 ≔ 6 x4 ≔ 6 y4 ≔ 6

KE 1 (CST KE): KE 2 (štapni KE):

E_CST ≔ 30 ⋅ 10 6 E_S ≔ 200 ⋅ 10 6 Nepoznata pomeranja: 5 i 6

30
ν ≔ 0.2 h ≔ 0.1 A ≔ ―― = 0.0030
100 2

KE 1 (CST KE)
čvorovi - mreža: 1, 2, 3
čvorovi - KE: 1, 2, 3
SS GKS: 1, 2, 3, 4, 5, 6
SS KE: 1, 2, 3, 4, 5, 6

⎡ 1 ν 0 ⎤ ⎡ 31250000.0 6250000.0 0.0 ⎤

E_CST ⎢ ⎥ ⎢
D ≔ ――― ⋅ ν 1 0 = 6250000.0 31250000.0 0.0 ⎥
1-ν 2 ⎢ ⎥ ⎢ ⎥
1-ν ⎥ ⎣
⎢ 0 0 ―― 0.0 0.0 12500000.0 ⎦
⎢⎣ 2 ⎥⎦

a1 ≔ x2 ⋅ y3 - x3 ⋅ y2 = 36 a2 ≔ x3 ⋅ y1 - x1 ⋅ y3 = 0.0 a3 ≔ x1 ⋅ y2 - x2 ⋅ y1 = 0.0
1
Ae ≔ ―⋅ ((a1 + a2 + a3)) = 18
2
b1 ≔ y2 - y3 = -6.0 b2 ≔ y3 - y1 = 6.0 b3 ≔ y1 - y2 = 0.0

c1 ≔ x3 - x2 = -6.0 c2 ≔ x1 - x3 = 0.0 c3 ≔ x2 - x1 = 6.0

1
 Zadatak 5

Matrica krutosti

Napomena:
određeni su samo elementi matrice krutosti KE uz nepoznata pomeranja

h
k55 ≔ ――⋅ ⎛b3 2 ⋅ D + c3 2 ⋅ D ⎞ = 625000.0
4 ⋅ Ae ⎜⎝ 0,0 2 , 2⎟
⎠

h
k56 ≔ ――⋅ ⎛b3 ⋅ c3 ⋅ ⎛D + D ⎞⎞ = 0.0 k65 ≔ k56 = 0.0
4 ⋅ Ae ⎝ ⎝ 0,1 2 , 2⎠⎠

h
k66 ≔ ――⋅ ⎛c3 2 ⋅ D + b3 2 ⋅ D ⎞ = 1562500.0
4 ⋅ Ae ⎜⎝ 1,1 2 , 2⎟
⎠
5 6
⎡ k55 k56 ⎤ ⎡ 625000.0 0.0 ⎤ ⎡5⎤
k_CST_aa ≔ ⎢ ⎥ =⎢
⎣ k65 k66 ⎦ ⎣ 0.0 1562500.0 ⎥⎦ ⎢⎣ 6 ⎥⎦

ili određivanje cele matrice krutosti

⎡ b1 0 b2 0 b3 0 ⎤ ⎡ -0.167 0.000 0.167 0.000 0.000 0.000 ⎤

1
B_CST ≔ ――⋅ ⎢ 0 c1 0 c2 0 c3 ⎥ = ⎢ 0.000 -0.167 0.000 0.000 0.000 0.167 ⎥
2 ⋅ Ae ⎢ c1 b1 c2 b2 c3 b3 ⎥ ⎢ -0.167 -0.167 0.000 0.167 0.167 0.000 ⎥
⎣ ⎦ ⎣ ⎦

k_CST ≔ B_CST T ⋅ D ⋅ B_CST ⋅ h ⋅ Ae

1 2 3 4 5 6
⎡ 21.87500 9.37500 -15.62500 -6.25000 -6.25000 -3.12500 ⎤ ⎡1⎤
⎢ 9.37500 21.87500 -3.12500 -6.25000 -6.25000 -15.62500 ⎥ ⎢2⎥
⎢ ⎥ ⎢ ⎥
k_CST
――― = ⎢ -15.62500 -3.12500 15.62500 0.00000 0.00000 3.12500 ⎥ ⎢3⎥
10 5 ⎢ -6.25000 -6.25000 0.00000 6.25000 6.25000 0.00000 ⎥ ⎢4⎥
⎢ -6.25000 -6.25000 0.00000 6.25000 6.25000 0.00000 ⎥ ⎢5⎥
⎢⎣ -3.12500 -15.62500 3.12500 0.00000 0.00000 15.62500 ⎥⎦ ⎢⎣ 6 ⎥⎦

Ekvivalentno opterećenje
L2
q ≔ 10000.0 L2 ≔ 6.0 Q3x ≔ q ⋅ h ⋅ ―― = 3000.0
2
ili
s
qx ((s)) ≔ q ⋅ h Ns3 ((s)) ≔ ――
L2
L2

Q3x ≔ ⌠
⌡ Ns3 ((s)) ⋅ qx ((s)) d s = 3000.0
0

⎡ Q3x ⎤ ⎡ 3000.0 ⎤ ⎡ 5 ⎤
Q_CST_a ≔ ⎢ =
⎣ 0 ⎥⎦ ⎢⎣ 0.0 ⎥⎦ ⎢⎣ 6 ⎥⎦

2
 Zadatak 5

KE 2 (štapni KE)
čvorovi - mreža: 3, 4
čvorovi - KE: 1, 2
SS GKS: 5, 6, 7, 8
SS KE: 1, 2, 3, 4

Matrica krutosti

Napomena:
određeni su samo elementi matrice krutosti KE uz nepoznata pomeranja

x4 - x3 y4 - y3
L ≔ ‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾‾
((x4 - x3)) 2 + ((y4 - y3)) 2 = 6.0 cos_α ≔ ――― = 1.0 sin_α ≔ ――― = 0.0
L L

E_S ⋅ A ⎛ E_S ⋅ A
k_S_11 ≔ ――― ⋅ ⎝cos_α 2 ⎞⎠ = 100000.0 k_S_12 ≔ ――― ⋅ ((cos_α ⋅ sin_α)) = 0.0
L L

E_S ⋅ A ⎛
k_S_21 ≔ k_S_12 = 0.0 k_S_22 ≔ ――― ⋅ ⎝sin_α 2 ⎞⎠ = 0.0
L

5 6
⎡ k_S_11 k_S_12 ⎤ ⎡ 100000.0 0.0 ⎤ ⎡5⎤
k_S_aa ≔ ⎢ =
⎣ k_S_21 k_S_22 ⎥⎦ ⎢⎣ 0.0 0.0 ⎥⎦ ⎢⎣ 6 ⎥⎦

Sistem KE
5 6
⎡ 725000.0 0.0 ⎤ ⎡ 5 ⎤ ⎡ 3000.0 ⎤ ⎡5⎤
Kaa ≔ k_CST_aa + k_S_aa = ⎢ Sa ≔ Q_CST_a = ⎢
⎣ 0.0 1562500.0 ⎥⎦ ⎢⎣ 6 ⎥⎦ ⎣ 0.0 ⎥⎦ ⎢⎣ 6 ⎥⎦

⎡ 0.0041379310 ⎤ ⎡ 5 ⎤
da ≔ Kaa -1 ⋅ Sa = ⎢
⎣ 0.0000000000 ⎥⎦ ⎢⎣ 6 ⎥⎦

Raspodela kinematičkih i statičkih veličina - KE 1

⎡ 0 ⎤
⎢ 0 ⎥ ⎡ 0.0000000000 ⎤ ⎡1⎤
⎢ ⎥ ⎢ 0.0000000000 ⎥ ⎢2⎥
⎢ 0 ⎥ ⎢ 0.0000000000 ⎥ ⎢ ⎥
d1 ≔ ⎢ 0 ⎥ = ⎢ ⎥ ⎢3⎥
⎢ da ⎥ ⎢ 0.0000000000 ⎥ ⎢4⎥
⎢ 0
⎥ ⎢ 0.0041379310 ⎥ ⎢5⎥
⎢ da ⎥ ⎢⎣ 0.0000000000 ⎥⎦ ⎢⎣ 6 ⎥⎦
⎣ 1⎦
float , 5
1 expand
N3 ((x , y)) ≔ ――⋅ ((a1 + b1 ⋅ x + c1 ⋅ y)) ――― → -0.16667 ⋅ y - 0.16667 ⋅ x + 1.0
2 ⋅ Ae
float , 5
expand
u ((x , y)) ≔ N3 ((x , y)) ⋅ d1 ――― → -0.00068967 ⋅ y - 0.00068967 ⋅ x + 0.0041379
4

float , 5
v ((x , y)) ≔ N3 ((x , y)) ⋅ d1 ――― → 0.0
5

3
 float , 5
v ((x , y)) ≔ N3 ((x , y)) ⋅ d1 ――― → 0.0 Zadatak 5
5

⎡ 0.0000000000 ⎤
Komponente deformacije: ε1 ≔ B_CST ⋅ d1 = ⎢ 0.0000000000 ⎥
⎢ ⎥
⎣ 0.0006896552 ⎦

⎡ b1 ⋅ D c1 ⋅ D b2 ⋅ D c2 ⋅ D b3 ⋅ D c3 ⋅ D ⎤
0,0 0,1 0,0 0,1 0,0 0,1
1 ⎢ ⎥
⎢
S ≔ ――⋅ b1 ⋅ D c1 ⋅ D b2 ⋅ D c2 ⋅ D b3 ⋅ D c3 ⋅ D ⎥
2 ⋅ Ae ⎢ 1,0 1,1 1,0 1,1 1,0 1,1
⎥
⎢ c1 ⋅ D b1 ⋅ D c2 ⋅ D b2 ⋅ D c3 ⋅ D b3 ⋅ D ⎥
⎣ 2,2 2,2 2,2 2,2 2,2 2,2⎦

⎡ -5208333.3 -1041666.7 5208333.3 0.0 0.0 1041666.7 ⎤

S = ⎢ -1041666.7 -5208333.3 1041666.7 0.0 0.0 5208333.3 ⎥
⎢ ⎥
⎣ -2083333.3 -2083333.3 0.0 2083333.3 2083333.3 0.0 ⎦

⎡ 0.000 ⎤ ⎡ 0.000 ⎤
Komponente napona: σ1 ≔ S ⋅ d1 = ⎢ 0.000 ⎥ ili σ1 ≔ D ⋅ ε1 = ⎢ 0.000 ⎥
⎢ ⎥ ⎢ ⎥
⎣ 8620.690 ⎦ ⎣ 8620.690 ⎦

⎡ 0.000 ⎤
sile ≔ σ1 ⋅ h = ⎢ 0.000 ⎥
⎢ ⎥
⎣ 862.069 ⎦
-ν -ν
εz ≔ ――⋅ ⎛ε1 + ε1 ⎞ = 0.0000000000 ili εz ≔ ――― ⋅ ⎛σ1 + σ1 ⎞ = 0.0000000000
1-ν ⎝ 0 1⎠ E_CST ⎝ 0 1⎠

Raspodela kinematičkih i statičkih veličina - KE 2

⎡ da ⎤
⎢ 0 ⎥ ⎡ 0.0041379310 ⎤ ⎡5⎤
⎢ da ⎥ ⎢ 0.0000000000 ⎥ ⎢6⎥
d_S_GKS ≔ ⎢ 1 ⎥ = ⎢ ⎥ ⎢ ⎥
⎢ 0.0000000000 ⎥ ⎢7⎥
⎢ 0.0 ⎥
⎢⎣ 0.0 ⎥⎦ ⎣ 0.0000000000 ⎦ ⎣8⎦

⎡ cos_α sin_α 0 0 ⎤ ⎡ 1.0 0.0 0.0 0.0 ⎤

T_S ≔ ⎢ =
⎣ 0 0 cos_α sin_α ⎥⎦ ⎢⎣ 0.0 0.0 1.0 0.0 ⎥⎦

⎡ 0.0041379310 ⎤
d_S_LKS ≔ T_S ⋅ d_S_GKS = ⎢ D_S ≔ E_S ⋅ A = 600000.0
⎣ 0.0000000000 ⎥⎦

⎡ 1 1⎤
B_S ≔ ⎢ -― ―⎥ = [[ -0.16667 0.16667 ]] N_S ≔ D_S ⋅ B_S ⋅ d_S_LKS = -413.793
⎣ L L⎦

E_S ⋅ A ⎡ 1 -1 ⎤ ⎡ 100000.0 -100000.0 ⎤

k_S_LKS ≔ ――― ⋅⎢ =
L ⎣ -1 1 ⎥⎦ ⎢⎣ -100000.0 100000.0 ⎥⎦

⎡ 413.793 ⎤
R_S ≔ k_S_LKS ⋅ d_S_LKS = ⎢
⎣ -413.793 ⎥⎦

4
 Zadatak 6

Odrediti komponente nepoznatih pomeranja za mrežu konačnih elemenata prikazanu na slici. Za

numeričku simulaciju formirati diskretni model sa jednim pravougaonim konačnim elementom koji
ima četiri čvora (KE 1) i jednim CST konačnim elementom (KE 2). Podaci: ravansko stanje napona;
E = 30 GPa; ν = 0,2; h = 0,2 m; L1 = 5,0 m; L2 = 5,0 m; L3 = 5,0 m; q = 20 MPa.

E ≔ 30 ⋅ 10 6 L1 ≔ 5
ν ≔ 0.2 L2 ≔ 5
h ≔ 0.2 L3 ≔ 5
q ≔ 20000

Nepoznata pomeranja: 3, 4, 5 i 6

⎡ 1 ν 0 ⎤ ⎡ 31250000.000 6250000.000 0.000 ⎤

E ⎢ ⎥ ⎢
D ≔ ――― ⋅ ν 1
2 ⎢
0 = 6250000.000 31250000.000 0.000 ⎥
1-ν ⎥ ⎢ ⎥
1-ν ⎥ ⎣
⎢ 0 0 ―― 0.000 0.000 12500000.000 ⎦
⎢⎣ 2 ⎥⎦

KE 1 - PRAVOUGAONI KE SA 4. ČVORA
čvorovi - mreža: 1, 2, 3, 4
čvorovi - KE: 1, 2, 3, 4
SS GKS: 1, 2, 3, 4, 5, 6, 7, 8
SS KE: 1, 2, 3, 4, 5, 6, 7, 8

L1
a ≔ ―― = 2.5
2

L2
b ≔ ―― = 2.5
2

Matrica krutosti

Napomena: određeni su samo elementi uz nepoznata pomeranja

h
k33 ≔ ――― ⋅ ⎛4 ⋅ ⎛b 2 ⋅ D + a 2 ⋅ D ⎞⎞ = 2916666.7
12 ⋅ a ⋅ b ⎜⎝ ⎜⎝ 0,0 2 , 2⎟
⎠⎟⎠

1
 Zadatak 6

h
k34 ≔ ――― ⋅ ⎛-3 ⋅ a ⋅ b ⋅ ⎛D + D ⎞⎞ = -937500.0 k43 ≔ k34 = -937500.0
12 ⋅ a ⋅ b ⎝ ⎝ 0,1 2 , 2⎠⎠

h
k35 ≔ ――― ⋅ ⎛2 ⋅ b 2 ⋅ D - 4 ⋅ a 2 ⋅ D ⎞ = 208333.3 k53 ≔ k35 = 208333.3
12 ⋅ a ⋅ b ⎜⎝ 0,0 2 , 2⎟
⎠

h
k36 ≔ ――― ⋅ ⎛3 ⋅ a ⋅ b ⋅ ⎛D - D ⎞⎞ = -312500.0 k63 ≔ k36 = -312500.0
12 ⋅ a ⋅ b ⎝ ⎝ 0,1 2 , 2⎠⎠

h
k44 ≔ ――― ⋅ ⎛4 ⋅ ⎛a 2 ⋅ D + b 2 ⋅ D ⎞⎞ = 2916666.7
12 ⋅ a ⋅ b ⎜⎝ ⎜⎝ 1,1 2 , 2⎟
⎠⎟⎠

h
k45 ≔ ――― ⋅ ⎛3 ⋅ a ⋅ b ⋅ ⎛-D + D ⎞⎞ = 312500.0 k54 ≔ k45 = 312500.0
12 ⋅ a ⋅ b ⎝ ⎝ 1,0 2 , 2⎠⎠

h
k46 ≔ ――― ⋅ ⎛-4 ⋅ a 2 ⋅ D + 2 ⋅ b 2 ⋅ D ⎞ = -1666666.7 k64 ≔ k46 = -1666666.7
12 ⋅ a ⋅ b ⎜⎝ 1,1 2 , 2⎟
⎠

k55 ≔ k33 = 2916666.7

k56 ≔ -k34 = 937500.0 k65 ≔ k56 = 937500.0

k66 ≔ k44 = 2916666.7

3 4 5 6
⎡ k33 k34 k35 k36 ⎤ ⎡ 2916666.7 -937500.0 208333.3 -312500.0 ⎤ ⎡ 3 ⎤
⎢ k43 k44 k45 k46 ⎥ ⎢ -937500.0 2916666.7 312500.0 -1666666.7 ⎥ ⎢ 4 ⎥
k_aa_KE1 ≔ ⎢ ⎥=⎢ ⎥ ⎢ ⎥
⎢ k53 k54 k55 k56 ⎥ ⎢ 208333.3 312500.0 2916666.7 937500.0 ⎥ ⎢ 5 ⎥
⎣ k63 k64 k65 k66 ⎦ ⎣ -312500.0 -1666666.7 937500.0 2916666.7 ⎦ ⎣ 6 ⎦

Ekvivalentno opterećenje
L1
Q3y_KE1 ≔ -q ⋅ h ⋅ ―― = -10000.0
2
ili
s
qy_KE1 ((s)) ≔ -q ⋅ h → -4000.0 Ns3_KE1 ((s)) ≔ ――
L1
L1

Q3y_KE1 ≔ ⌠
⌡ Ns3_KE1 ((s)) ⋅ qy_KE1 ((s)) d s = -10000.0
0

ili

y_opt_ivica ≔ b

1 ⎛ x⎞ ⎛ y_opt_ivica ⎞ explicit , ALL 1 ⎛ x ⎞ ⎛ 2.5 ⎞

N3_KE1 ((x)) ≔ ―⋅ ⎜1 + ―⎟ ⋅ ⎜1 + ――――⎟ ――――→ ―⋅ ⎜1 + ―― ⎟ ⋅ ⎜1 + ―― ⎟ → 0.2 ⋅ x + 0.5
4 ⎝ a⎠ ⎝ b ⎠ 4 ⎝ 2.5 ⎠ ⎝ 2.5 ⎠

1 ⎛ x⎞ ⎛ y_opt_ivica ⎞ explicit , ALL 1 ⎛ x ⎞ ⎛ 2.5 ⎞

N4_KE1 ((x)) ≔ ―⋅ ⎜1 - ―⎟ ⋅ ⎜1 + ――――⎟ ――――→ ―⋅ ⎜1 - ―― ⎟ ⋅ ⎜1 + ―― ⎟ → -0.2 ⋅ x + 0.5
4 ⎝ a⎠ ⎝ b ⎠ 4 ⎝ 2.5 ⎠ ⎝ 2.5 ⎠

2
 Zadatak 6

⎡ 0 0 0 0 N3_KE1 ((x)) 0 N4_KE1 ((x)) 0 ⎤

Ns_KE1 ((x)) ≔ ⎢
⎣ 0 0 0 0 0 N3_KE1 (x)
( ) 0 N4_KE1 ( ) ⎥⎦
(x)

⎡ 0 ⎤ ⎡ 0.0 ⎤
qb_KE1 ((x)) ≔ ⎢ =⎢
⎣ -q ⋅ h ⎦ ⎣ -4000.0 ⎥⎦
⎥

⎡ 0.0 ⎤ ⎡1⎤
⎢ 0.0 ⎥ ⎢2⎥
⎢ ⎥⎢ ⎥
⎢ 0.0 ⎥ ⎢3⎥
0.0 ⎥ ⎢4⎥
podInt ((x)) ≔ ((Ns_KE1 ((x)))) T ⋅ qb_KE1 ((x)) → ⎢
⎢ 0.0 ⎥ ⎢5⎥
⎢ -800.0 ⋅ x - 2000.0 ⎥ ⎢ 6 ⎥
⎢ ⎥⎢ ⎥
0.0 7
⎢ ⎥⎢ ⎥
⎣ 800.0 ⋅ x - 2000.0 ⎦ ⎣ 8 ⎦
a

Q3y_KE1 ≔ ⌠ podInt ((x)) d x = -10000.0

⎮ 5
⌡
-a

KE 2 - TROUGAONI KE SA 3. ČVORA - CST KE

čvorovi - mreža: 5, 3, 2
čvorovi - KE: 1, 2, 3
SS GKS: 9, 10, 5, 6, 3, 4
SS KE: 1, 2, 3, 4, 5, 6

Prvi čvor KE (čvor 5)

x1 ≔ L1 + L3 = 10.0 y1 ≔ L2 = 5.0

Drugi čvor KE (čvor 3)

x2 ≔ L1 = 5.0 y2 ≔ L2 = 5.0

Treći čvor KE (čvor 2)

x3 ≔ L1 = 5.0 y3 ≔ 0

a1 ≔ x2 ⋅ y3 - x3 ⋅ y2 = -25.0 a2 ≔ x3 ⋅ y1 - x1 ⋅ y3 = 25.0 a3 ≔ x1 ⋅ y2 - x2 ⋅ y1 = 25.0

1
Ae ≔ ―⋅ ((a1 + a2 + a3)) = 12.5
2
b1 ≔ y2 - y3 = 5.0 b2 ≔ y3 - y1 = -5.0 b3 ≔ y1 - y2 = 0.0

c1 ≔ x3 - x2 = 0.0 c2 ≔ x1 - x3 = 5.0 c3 ≔ x2 - x1 = -5.0

Matrica krutosti

Napomena: određeni su samo elementi matrice krutosti KE uz nepoznata pomeranja

h
k33 ≔ ――⋅ ⎛b2 2 ⋅ D + c2 2 ⋅ D ⎞ = 4375000.0
4 ⋅ Ae ⎜⎝ 0,0 2 , 2⎟
⎠

h
k34 ≔ ――⋅ ⎛b2 ⋅ c2 ⋅ ⎛D + D ⎞⎞ = -1875000.0 k43 ≔ k34 = -1875000.0
4 ⋅ Ae ⎝ ⎝ 0,1 2 , 2⎠⎠

3
 Zadatak 6

h
k35 ≔ ――⋅ ⎛b2 ⋅ b3 ⋅ D + c2 ⋅ c3 ⋅ D ⎞ = -1250000.0 k53 ≔ k35 = -1250000.0
4 ⋅ Ae ⎝ 0,0 2 , 2⎠

h
k36 ≔ ――⋅ ⎛b2 ⋅ c3 ⋅ D + b3 ⋅ c2 ⋅ D ⎞ = 625000.0 k63 ≔ k36 = 625000.0
4 ⋅ Ae ⎝ 0,1 2 , 2⎠

h
k44 ≔ ――⋅ ⎛c2 2 ⋅ D + b2 2 ⋅ D ⎞ = 4375000.0
4 ⋅ Ae ⎜⎝ 1,1 2 , 2⎟
⎠

h
k45 ≔ ――⋅ ⎛b3 ⋅ c2 ⋅ D + b2 ⋅ c3 ⋅ D ⎞ = 1250000.0 k54 ≔ k45 = 1250000.0
4 ⋅ Ae ⎝ 1,0 2 , 2⎠

h
k46 ≔ ――⋅ ⎛c2 ⋅ c3 ⋅ D + b2 ⋅ b3 ⋅ D ⎞ = -3125000.0 k64 ≔ k46 = -3125000.0
4 ⋅ Ae ⎝ 1,1 2 , 2⎠

h
k55 ≔ ――⋅ ⎛b3 2 ⋅ D + c3 2 ⋅ D ⎞ = 1250000.0
4 ⋅ Ae ⎜⎝ 0,0 2 , 2⎟
⎠

h
k56 ≔ ――⋅ ⎛b3 ⋅ c3 ⋅ ⎛D + D ⎞⎞ = 0.0 k65 ≔ k56 = 0.0
4 ⋅ Ae ⎝ ⎝ 0,1 2 , 2⎠⎠

h
k66 ≔ ――⋅ ⎛c3 2 ⋅ D + b3 2 ⋅ D ⎞ = 3125000.0
4 ⋅ Ae ⎜⎝ 1,1 2 , 2⎟
⎠

SS GKS: 9, 10, 5, 6, 3, 4
SS KE: 1, 2, 3, 4, 5, 6

5 6 3 4 3 4 5 6
⎡ k33 k34 k35 k36 ⎤ ⎡ 5 ⎤ ⎡ k55 k56 k53 k54 ⎤ ⎡ 3 ⎤
⎢ k43 k44 k45 k46 ⎥ ⎢ 6 ⎥ ⎢ k65 k66 k63 k64 ⎥ ⎢ 4 ⎥
k_aa_KE2 ≔ ⎢ ⎥⎢ ⎥ k_aa_KE2 ≔ ⎢ ⎥⎢ ⎥
⎢ k53 k54 k55 k56 ⎥ ⎢ 3 ⎥ ⎢ k35 k36 k33 k34 ⎥ ⎢ 5 ⎥
⎣ k63 k64 k65 k66 ⎦ ⎣ 4 ⎦ ⎣ k45 k46 k43 k44 ⎦ ⎣ 6 ⎦

3 4 5 6
⎡ 1250000.0 0.0 -1250000.0 1250000.0 ⎤ ⎡ 3 ⎤
⎢ 0.0 3125000.0 625000.0 -3125000.0 ⎥ ⎢ 4 ⎥
k_aa_KE2 = ⎢ ⎥ ⎢ ⎥
⎢ -1250000.0 625000.0 4375000.0 -1875000.0 ⎥ ⎢ 5 ⎥
⎣ 1250000.0 -3125000.0 -1875000.0 4375000.0 ⎦ ⎣ 6 ⎦
ili određivanje cele matrice krutosti

⎡ b1 0 b2 0 b3 0 ⎤ ⎡ 0.200 0.000 -0.200 0.000 0.000 0.000 ⎤

1
B_KE2 ≔ ――⋅ ⎢ 0 c1 0 c2 0 c3 ⎥ = ⎢ 0.000 0.000 0.000 0.200 0.000 -0.200 ⎥
2 ⋅ Ae ⎢ c1 b1 c2 b2 c3 b3 ⎥ ⎢ 0.000 0.200 0.200 -0.200 -0.200 0.000 ⎥
⎣ ⎦ ⎣ ⎦

k_KE2 ≔ B_KE2 T ⋅ D ⋅ B_KE2 ⋅ h ⋅ Ae

4
 Zadatak 6

1 2 5 6 3 4
⎡ 31.25000 0.00000 -31.25000 6.25000 0.00000 -6.25000 ⎤ ⎡ 1 ⎤
⎢ 0.00000 12.50000 12.50000 -12.50000 -12.50000 0.00000 ⎥ ⎢ 2 ⎥
⎢ ⎥⎢ ⎥
k_KE2 -31.25000 12.50000 43.75000 -18.75000 -12.50000 6.25000 ⎥ ⎢ 5 ⎥
――― = ⎢
10 5 ⎢ 6.25000 -12.50000 -18.75000 43.75000 12.50000 -31.25000 ⎥ ⎢ 6 ⎥
⎢ 0.00000 -12.50000 -12.50000 12.50000 12.50000 0.00000 ⎥ ⎢ 3 ⎥
⎢⎣ -6.25000 0.00000 6.25000 -31.25000 0.00000 31.25000 ⎥⎦ ⎢⎣ 4 ⎥⎦

Ekvivalentno opterećenje

L3 L3
Q3y_KE2 ≔ -q ⋅ h ⋅ ―― = -10000.0 Q5y_KE2 ≔ -q ⋅ h ⋅ ―― = -10000.0
2 2
ili
s
qy_KE2 ((s)) ≔ -q ⋅ h → -4000.0 Ns3_KE2 ((s)) ≔ 1 - ――
L3
L1

Q3y_KE2 ≔ ⌠
⌡ Ns3_KE2 ((s)) ⋅ qy_KE2 ((s)) d s = -10000.0
0

Sistem KE

3 4 5 6
⎡ 4166666.7 -937500.0 -1041666.7 937500.0 ⎤ ⎡3⎤
⎢ -937500.0 6041666.7 937500.0 -4791666.7 ⎥ ⎢4⎥
Kaa ≔ k_aa_KE1 + k_aa_KE2 = ⎢ ⎥ ⎢ ⎥
⎢ -1041666.7 937500.0 7291666.7 -937500.0 ⎥ ⎢5⎥
⎣ 937500.0 -4791666.7 -937500.0 7291666.7 ⎦ ⎣6⎦

⎡ 0 ⎤ ⎡ 0.0 ⎤ ⎡ 3 ⎤ ⎡ 0.0 ⎤ ⎡3⎤

⎢ 0 ⎥ ⎢ 0.0 ⎥ ⎢ 4 ⎥ ⎢ 0.0 ⎥ ⎢4⎥
Q_a ≔ ⎢ ⎥=⎢ ⎥⎢ ⎥ Sa ≔ Q_a = ⎢ ⎥ ⎢ ⎥
⎢ 0 ⎥ ⎢ 0.0 ⎥ ⎢ 5 ⎥ ⎢ 0.0 ⎥ ⎢5⎥
⎣ Q3y_KE1 + Q3y_KE2 ⎦ ⎣ -20000.0 ⎦ ⎣ 6 ⎦ ⎣ -20000.0 ⎦ ⎣6⎦

⎡ 0.000249513 ⎤ ⎡ 3 ⎤
⎢ -0.004501624 ⎥ ⎢ 4 ⎥
da ≔ Kaa -1 ⋅ Sa = ⎢ ⎥⎢ ⎥
⎢ -0.000124756 ⎥ ⎢ 5 ⎥
⎣ -0.005749188 ⎦ ⎣ 6 ⎦