Finte Elerment Analysis

NamevSau
ASsignmet-2
Roll No,- ME 20BlO4
f ='lm

k-sm

Ruleigh - RAz Neltod:

EA du -P
Bourdery ondons; uto) -0, d l>l

TOlay plental

Tolal polertal

Coob-o-l)
Trail eoution:
otonS,
Appy

du q+ 2ct
 o-6-o-l)tCGt 20,1)°dr

(o6-ol1) CG2c,
(2s

25c4Y -i0G61)
2
3

- 2"5C4

Co-6-o1)Cqi16,)dt

olGr]
6

2cL - 2c,
t [L- 3.33L>JGc
+[20L- 2sL]G
-I·25c
Tp = [ISL
 L-284+20ug-ç)
-2l3l2]t
+ [u0
2 [&o-2o] - 24.92

t -h2.8s+
[u0- 213:27 G8-4O82-32
Glo]Ct
2[1280- 124-5

1280

= Obeeell

e-eeetNY Cin m)
Dis placemen
t
U)=O3234

I3.4|

pa)
q.uat 13.411 Ch
 OOxO3
E2A b
(2300X
O
342xlo
IO3 YOx L
200000
xIO 8"4 x
0enoting,
have aslo we
Bo
m
DistonCe
2
Npa12300 Es
000 200o Ei
MPC
n(s) A
embes node
fosmuha. Easily problem
ie Tho
me aefining by
=1230ONPa E
 chose penaltq consant

Modfed eAffnesS qiven by

b
-a

-b

btc
 - I681S60

\68" 4Sb -Qo0S

\684564:305
-Q"305
- \684560

-0-sSmm

200000AoCoouaxo
ECa-a2)

Tabe i 12300x0 C-ous}o

FoY Slement 0

200o
EAoT

EA)45T

L
foon elemeni
Ga
-6

Globcl Nals

OO
-244

-0ol| 232 P
tO-ll232-O|sUSG
O-\SUG - 0-30|qOY
o.30|l9Oy

using limiation ApDoreh,

-O-422

l-40 38
 6.03276
-00%66a

-6)
-uO 438^1Q, = -O-O% 6604

- 21239 mm
-0-00213Am

= EB

* lo-6

200KO

26 2-03o165x O pa

00

Eleme :
Yess in

6co63

- 221·24
6COxIg3

600
Sess in elere -

83No

-6
83xo

= 22.9S3SI0

poat ReoionSi

-21gicqN

-5S,08g"SlN

24-S0Orn SOomn

)
6-2000OOHpa 3
 NOdeA NOde LCOSo)
1

1 3

FO e r e ;

A 73

o
FOS Slcnen
2 S
-o-64 to87
-O6
-o6y +0u3
|tO- -0-36 J6

FoS Hefnert 3

-3

-o-64 40-% o
 S
+0uB -O36
3

S
- ob4
-036

=0,

USing tltoinatfon appcach:

|.64
-)00Oo

- 2 0Y

-1667xoY

-l-667

-O"75q25 m
 by S3:.
nodos I 3
sleronH connecAod

[-a
le

Sx0xo3

[1o433)
+0"4r0-43xot

pa 333-MP
 A:30Omm 60omm

20k

LCoeo) mcsingJ
Hement Nde! Node 2
-0-6

Foi lene1: 3
-O-/8-O"Bb
2x0:n4 +B6
3

-o36 O"64 4
-0"644
O48

FOS elerent 2
-0-64
-Oe

-064 6
C36
fo elernont

-0R -0-64
10x o3

-0U -o36

Gobal Nair 6
3 -O-6u
o"64
6"6s -08S3
o-64 -0853
2.061

-6
3
o64 -b853

-o64
o-64

-o64
-08s3
-0"64

Appich,
ey using etinircden

6xl0
2·067
 -o:0rmm

I=-020. o-OS2 ¥I0 m

eSS in lerNert ;

le

200x0 -o- -o-6 o-&

2C0x(o? cot)Co-os) +(o8) Co- |14))

o-264x\o [o-0312 0-34]

te -
 lo3 co6)(o-cs»)

O-2y10 -o-032t o-1392]

2"6Pg

200o9

o-261o,:0

28.836 Hpa
clc
clear all
E 200e9;
f= 77000;
a =0.1;
b = 0.025;
n

%n = 8
%n = 16;
1 = 1/n;
A = zeros (n, 1);
for i = 1:2:2*n
A((i+1)/2,1) = (((2*n-i)*a + i*b)/double (2*n) )^2;
end
K= zeros (n+1, n+1) ;
Q(1,1) = 0;
S = zeros (n+1, 1);
S(1,1) = f"A(1,1)*1/2; S(n+1,1) = f*A(n,1)*1/2;
for 1 2:n
S(1,1) - f*A(i-1,1)#l/2 + f*A(i,1)*1/2;
end
K(1,1) = E"A(1)/l;
for i = 1:n
if 1 n
K(1,i+1) = -E"A(1)/1;
K(i+1, i) -E*A(1)/1; K(i+1, i+1) = E"A(1)/1 + E*A(i+1)/1;
else
K(i,i+1) -E*A(1)/1;
K(1+1,1) = -*A(1)/1; K(1+1, i+1) = EA(1)/1;
end
end
X - transpose (S(2: n+1,1))/K(2:n+1,2: n+1);
R = (K(1, 2: n+1) *transpose(X)) - s(1, 1);
Q = zeros (n+1, 1);
for i = 2:n+1
Q(1) - x(i-1);
end
A % AREA OF EACH ELEMENT

A 4x1
0.0082
0.0052
8.0028
0.0012

S % NODAL FORCE VECTOR

S= 5x1
79.0491
128.7720
76.8872
38.5376
11.3733

K % STIFFNESS MATRIX
K = 5x5

1010 x
0.6570 -0.6570
-0.6570 1.0703 -0.4133
-0.4133 0.6391 -0.2258
-0.2258 0.3203 -0.0945
-0.0945 0.0945

Q % NODAL DISPLACEMENT MATRIX

Q= 5x1
10-6 x

0.0389
0.0696
0.0917
0.1037

R % REACTION FORCE

R= -334.6191
clc
clear all
E = 200e9;
f= 77000;
a = 0.1;
b = 0.025;
%n = 4 ;
n =8
%n = 16;
1 = 1/0;
A = zeros(n,1);
for i = 1:2:2*n
A((i+1)/2, 1) = (((2*n-i) *a + 1*b)/double(2*n))^2;
end
K = zeros (n+1, n+1) ;
Q(1,1) 0 ;
S = zeros (n+1,1);
S(1,1) = f*A(1,1)*1/2; S(n+1,1) = f*A(n,1)*1/2;
for 1 = 2:n
S(1,1) - f*A(i-1,1)#l/2 + f*A(i,1)*1/2;
end
K(1,1) = E*A(1)/1;
for i = 1:n

K(i, i+1) = -E*A(i)/1;

K(i+1, i) -E*A(1)/1; K(i+1, i+1) = E"A(1)/1 + E*A(i+1)/1;
else
K(i,i+1) -E*A(i)/l;
K(1+1,1) = -*A(1)/1; K(1+1, i+1) = E"A(1)/1;
end
end
X = transpose (S(2: n+1,1))/K(2:n+1,2: n+1);
R= (K(1, 2:n+1)*transpose(X)) - s(1,1);
Q = zeros (n+1, 1);
for i = 2:n+1
Q(1) = Xx(i-1);
end
A % AREA OF EACH ELEMENT

A 8x1
0.0091
0.0074
0.0059
0.0045
0.0033
0.0023
0.0015
0.0009

S % NODAL FORCE VECTOR

S = 9x1
 43.7190
79.2606
63.7515
49.9344
37.8091
27.3758
18.6343
11.5848
4.2415

K % STIFFNESS MATRIX

K= 9x9
1010
1.4535 -1.4535
-1.4535 2.6352 -1.1816
-1.1816 2.1195 -0.9379
-0.9379 1.6602 -0.7223
-0.7223 1.2570 -0.5348
-0.5348 0.9102 -0.3754
-0.3754 0.6195 -0.2441
-0.2441 0.3852 -0.1410
-0.1410 8.1410

0 % NODAL DISPLACEMENT MATRIX

Q = 9x1
10- x

0.2013
0.3818
0.5413
0.6793
0.7949
0.8867
0.9515
0.9816

R % REACTION FORCE

R = -336.3110
clc
clear all
E = 200e9;
f= 77000;
a = 0.1;
b = 0.025;
%n 4;
%n = 8 ;
n 16;
1 = 1/n;
A = zeros (n, 1);
for i = 1:2:2*n
A((i+1)/2,1) = (((2*n-1) *a + 1*b)/double(2*n) )^2;
end
K = zeros (n+1, n+1);
Q(1,1) - 8;
S = zeros (n+1, 1);
S(1,1) = f*A(1,1) *1/2; S(n+1,1) = ftA(n,1) *1/2;
for i = 2:n
S(1,1) - f*A(i-1,1) *1/2 + f*A(1,1) *1/2;
end
K(1,1) = E"A(1)/1;
for i = 1:n
if 1 n
K(1,i+1) = -E*A(1)/1;
K(i+1, i) -E*A(1)/1; K(i+1, i+1) = E"A(1)/1 + E*A(i+1)/1;
else
K(i,i+1) = -E*A(i)/1;
K(1+1,1) = -*A(1)/1; K(1+1, i+1) = EA(1)/1;
end
end
X = transpose (S(2: n+1, 1) )/K(2:n+1,2: n+1);
R= (K(1, 2:n+1)*transpose(X)) - s(1,1);
Q = zeros (n+1, 1);
for i = 2:n+1
Q(1) - x(i-1);
end
A % AREA OF EACH ELEMENT

A 16x1
0.0095
Ø.0086
0.0078
0.007
0.0062
e.0055
B.0048
0.0042
e.0036
e.e031

S % NODAL FORCE VECTOR

S = 17x1
 22.9478
43.7455
39.551
35.5680
31.7965
28.2364
24.8879
21.7508
18.8253
16.1112

K % STIFFNESS MATRIX

K = 17x17

101 x
3.0518 -3.0518

11144
-3.0518 5.8176 -2.7658
-2.7658 5.2598 -2.4939
-2.4939 4.7301 -2.2361
-2.2361 4.2285 -1.9924
-1.9924 3.7551 -1.7627
-1,7627 3.3098 -1.5471
-1.5471 2.8926 -1.3455
-1.3455 2.5035 -1

-1.158O 2

0 % NODAL DISPLACEMENT MATRIX

Q = 17x1
10- x

0.1028
0.2005
0.2929
0.3800
0.4619
0.5384
0.6095
0.6751
0.7351

R % REACTION FORCE

R = -336. 7340