Scribd
Upload
45 views5 pages

Fem PDF

1) The deflections at each joint are calculated by solving the global stiffness matrix. 2) The reactions at nodes A, D, and H are 2840 N, 13940 N, and 115100 N respectively. 3) The average normal stresses in elements 2, 5, and 8 are -1.4354 MPa, -1.4354 MPa, and 1.15 MPa respectively.

Uploaded by

v ahmad
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
45 views5 pages

Fem PDF

1) The deflections at each joint are calculated by solving the global stiffness matrix. 2) The reactions at nodes A, D, and H are 2840 N, 13940 N, and 115100 N respectively. 3) The average normal stresses in elements 2, 5, and 8 are -1.4354 MPa, -1.4354 MPa, and 1.15 MPa respectively.

Uploaded by

v ahmad
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

FEM

(H.W2)

Ahmad nedal saleh


0160328
First: We discretize the problem into nodes and elements.

Second: We develop equations for elements and solve for each element.
Third: Assemble elements into a Global Matrix and Apply Boundary Conditions. I did this using MATLAB
code.
clc
clear all
close all
k1=2.5*10^6;
k2=1.9*10^6;
k3=1.79*10^6;
k4=2.41*10^6;
k5=1.52*10^6;
l1=0.5;
m1=0.5;
x1=0.5;
l4=0.899317755;
m4=0.100682245;
x4=0.300907512;
l7=0.360001715;
m7=0.639998285;
x7=0.4800005;
m9=m7;
l9=l7;
x9=-x7;
l10=0.899317755;
m10=0.100682245;
x10=-0.300907512;
l12=0.5;
m12=0.5;
x12=-0.5;
k=[k3*1+k1 x1 -k3*l1 -x1 -k1 0 0 0 0 0 0 0 0 0 0 0;
k3*x1 k3*m1 -k3*x1 -k3*m1 0 0 0 0 0 0 0 0 0 0 0 0;
-k3*l1 k3*x1 k3*l1+k4*l4 k4*x4+k3*x1 0 0 0 0 -k4*l4 -k4*x4 0 0 0 0
0 0;
-k3*x1 -k3*m1 k4*m4+k3*x1 k4*m4+k3*m1-k1 0 k1 0 0 -k4*x4 -k4*m4 0
0 0 0 0 0;
-k1 0 0 0 2*k1+k5*l7 k5*x7 -k1 0 -k5*l7 -k5*x7 0 0 0 0 0 0;
0 0 0 k1 k5*x7 k5*m7-k1 0 0 -k5*x7 -k5*m7 0 0 0 0 0 0;
0 0 0 0 -k1 0 2*k1 0 0 0 -k1 0 0 0 0 0;
0 0 0 0 0 0 0 k2 0 -k2 0 0 0 0 0 0;
0 0 -k4*l4 -k4*x4 -k5*l7 -k5*x7 0 0 k5*l7+k4*l4-k4*l10-k5*l9
k5*x7+k4*x4-k4*x10-k5*x9 k5*l9 k5*x9 k4*l10 k4*x10 0 0;
0 0 -k4*x4 -k4*m4 -k5*x7 -k5*m7 0 -k2 k4*x4-k4*x10+k5*x7-k5*x9
k4*m4+k2-k4*m10+k5*m7-k5*m9 k5*x9 k5*m9 k4*x10 k4*m10 0 0;
0 0 0 0 0 0 -k1 0 k5*l9 k5*x9 -k5*l9+2*k1 -k5*x9 0 0 -k1 0;
0 0 0 0 0 0 0 0 k5*x9 k5*m9 -k5*x9 -k5*m9+k1 0 -k1 0 0;
0 0 0 0 0 0 0 0 k4*l10 k4*x10 0 0 -k4*l10-k3*l12 -k4*x10-k3*x12
k3*l12 k3*x12;
0 0 0 0 0 0 0 0 k4*x10 k4*m10 0 -k1 -k4*x10-k3*x12 -k3*m12+k1-
k4*m10 k3*x12 k3*m12;
0 0 0 0 0 0 0 0 0 0 -k1 0 k3*l12 k3*x12 -k3*l12 -k3*x12;
0 0 0 0 0 0 0 0 0 0 0 0 k3*x12 k3*m12 -k3*x12 -k3*m12];
k(2,:)=[];
k(8:9,:)=[];
k(12:13,:)=[];

k(:,2)=[];
k(:,8:9)=[];
k(:,12:13)=[];

f=[0 -50000 0 0 -30000 0 -30000 0 -90000 0 0]';


kk=inv(k);
u=kk*f;

u2=[0.0079;
0;
-0.0126;
-0.0114;
-0.0091;
-0.0033;
-0.0131;
-0.0158;
0;
0;
-0.0172;
0.0601;
0.0358;
0.0677;
0;0]

K2=[k3*1+k1 x1 -k3*l1 -x1 -k1 0 0 0 0 0 0 0 0 0 0 0;


k3*x1 k3*m1 -k3*x1 -k3*m1 0 0 0 0 0 0 0 0 0 0 0 0;
-k3*l1 k3*x1 k3*l1+k4*l4 k4*x4+k3*x1 0 0 0 0 -k4*l4 -k4*x4 0 0 0 0
0 0;
-k3*x1 -k3*m1 k4*m4+k3*x1 k4*m4+k3*m1-k1 0 k1 0 0 -k4*x4 -k4*m4 0
0 0 0 0 0;
-k1 0 0 0 2*k1+k5*l7 k5*x7 -k1 0 -k5*l7 -k5*x7 0 0 0 0 0 0;
0 0 0 k1 k5*x7 k5*m7-k1 0 0 -k5*x7 -k5*m7 0 0 0 0 0 0;
0 0 0 0 -k1 0 2*k1 0 0 0 -k1 0 0 0 0 0;
0 0 0 0 0 0 0 k2 0 -k2 0 0 0 0 0 0;
0 0 -k4*l4 -k4*x4 -k5*l7 -k5*x7 0 0 k5*l7+k4*l4-k4*l10-k5*l9
k5*x7+k4*x4-k4*x10-k5*x9 k5*l9 k5*x9 k4*l10 k4*x10 0 0;
0 0 -k4*x4 -k4*m4 -k5*x7 -k5*m7 0 -k2 k4*x4-k4*x10+k5*x7-k5*x9
k4*m4+k2-k4*m10+k5*m7-k5*m9 k5*x9 k5*m9 k4*x10 k4*m10 0 0;
0 0 0 0 0 0 -k1 0 k5*l9 k5*x9 -k5*l9+2*k1 -k5*x9 0 0 -k1 0;
0 0 0 0 0 0 0 0 k5*x9 k5*m9 -k5*x9 -k5*m9+k1 0 -k1 0 0;
0 0 0 0 0 0 0 0 k4*l10 k4*x10 0 0 -k4*l10-k3*l12 -k4*x10-k3*x12
k3*l12 k3*x12;
0 0 0 0 0 0 0 0 k4*x10 k4*m10 0 -k1 -k4*x10-k3*x12 -k3*m12+k1-
k4*m10 k3*x12 k3*m12;
0 0 0 0 0 0 0 0 0 0 -k1 0 k3*l12 k3*x12 -k3*l12 -k3*x12;
0 0 0 0 0 0 0 0 0 0 0 0 k3*x12 k3*m12 -k3*x12 -k3*m12]';

f2=[0 0 -50000 0 0 -30000 0 -30000 0 0 0 -90000 0 0 0 0]';


R=K2*u2-f2
(1) Determine the deflection of each joint:

Ux,y of element 1 Ux,y of element 5


0.0079 0 0 0
Ux,y of element 2 Ux,y of element 6
-0.0126 -0.0114 -0.0172 0.0601
Ux,y of element 3 Ux,y of element 7
-0.0091 -0.0033 0.0358 0.0677
Ux,y of element 4 Ux,y of element 8
-0.0131 -0.0158 0 0

(2) The Reactions at A, D and H respectively

Node 1 Node 5 Node 8


0 2840 13940 32060 36060 115100

(3) The Average Normal Stress in any three members (MPa):

Element 2 Element 5 Element 8


-1.4354 -1.4354 1.15

You might also like