Scribd
Upload
15 views

Assignment4 Ce23mtech11029

This document contains code for analyzing the force-displacement behavior of a three member frame structure under varying axial loads. It calculates the stiffness matrices of each member and assembles them into a global stiffness matrix to solve for displacements under different load cases.

Uploaded by

2006 A.Saikumar
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
15 views

Assignment4 Ce23mtech11029

This document contains code for analyzing the force-displacement behavior of a three member frame structure under varying axial loads. It calculates the stiffness matrices of each member and assembles them into a global stiffness matrix to solve for displacements under different load cases.

Uploaded by

2006 A.Saikumar
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/ 18

10/2/23 7:50 PM C:\Users\DELL\Desktop\ASM-ASS...\ASS4A.

m 1 of 3

% Force (P) vs lateral displacement behavior if V is equal to 1 Newton


clear all;
format short e;

for n=1:151
%ISMB500 Section (COLUMNS & BEAM)
A1 = 11074;A2 = 11074;A3 = 11074;
I1 = 452183000;I2 = 452183000;I3 = 452183000;
L1 = 20000;L2 =10000;L3 =20000;
E1= 200000;E2= 200000;E3= 200000;

P=(n-1)*10000; % increase axial force in increments of 10kN.

% Linear stiffness matrix of column member


SM01 = [A1*E1/L1 0 0 -A1*E1/L1 0,
0 12*E1*I1/L1^3 6*E1*I1/L1^2 0 -12*E1*I1/L1^3,
0 6*E1*I1/L1^2 4*E1*I1/L1 0 -6*E1*I1/L1^2,
-A1*E1/L1 0 0 A1*E1/L1 0,
0 -12*E1*I1/L1^3 -6*E1*I1/L1^2 0 12*E1*I1/L1^3];

% Geometric stiffness matrix of a column


SG1 = P/L1*[0 0 0 0 0,
0 6/5 L1/10 0 -6/5,
0 L1/10 2*L1^2/15 0 -L1/10,
0 0 0 0 0 ,
0 -6/5 -L1/10 0 6/5];

% Combines stiffness matrix of a column


SM1 = SM01-SG1;
SM2 = [A2*E2/L2 0 -A2*E2/L2 0,
0 12*E2*I2/L2^3 0 -12*E2*I2/L2^3,
-A2*E2/L2 0 A2*E2/L2 0,
0 -12*E2*I2/L2^3 0 12*E2*I2/L2^3];

SM03 = [A3*E3/L3 0 -A3*E3/L3 0 0,


0 12*E3*I3/L3^3 0 -12*E3*I3/L3^3 6*E3*I3/L3^2,
-A3*E3/L3 0 A3*E3/L3 0 0,
0 -12*E3*I3/L3^3 0 12*E3*I3/L3^3 -6*E3*I3/L3^2,
0 6*E3*I3/L3^2 0 -6*E3*I3/L3^2 4*E3*I3/L3];

SG3 = P/L3*[0 0 0 0 0,
0 6/5 0 -6/5 L3/10,
0 0 0 0 0,
0 -6/5 0 6/5 -L3/10,
0 L3/10 0 -L3/10 2*L3^2/15];
10/2/23 7:50 PM C:\Users\DELL\Desktop\ASM-ASS...\ASS4A.m 2 of 3

SM3 = SM03-2*SG3;

% Transformation matrix for a vertical upward column (member 1)


RT1=[0 1 0 0 0,
-1 0 0 0 0,
0 0 1 0 0,
0 0 0 0 1,
0 0 0 -1 0];

% Transformation matrix for a vertical downward column (member 3)


RT3=[0 -1 0 0 0,
1 0 0 0 0,
0 0 0 -1 0,
0 0 1 0 0,
0 0 0 0 1];

% Transformed stiffness matrix for first column (member 1)


SMS_C1 = RT1'*SM1*RT1;

% Transformed stiffness matrix for second column (member 3)


SMS_C3 = RT3'*SM3*RT3;

% no geomtric stiffness matrix required for the beam as there is negligible


% axial force on the beam

S_Total (1:10,1:10) = 0; % initialize the total stiffness matrix

% Assemble the global stiffness matrix


S_Total(1:5,1:5)=S_Total(1:5,1:5)+ SMS_C1;
S_Total(4:7,4:7)=S_Total(4:7,4:7)+ SM2;
S_Total(6:10,6:10)=S_Total(6:10,6:10)+SMS_C3;

%SWAY (FDOF = 4, 5, 6, 7, 8, 9)
SFF(1:4,1:4)= 0;
SFF(1:4,1:4)= S_Total(4:7,4:7);

AF = [1 -2*P 0 -P]'; % force vector


DF=inv(SFF)*AF; % calculate displacements at FDOFs

Force(n) = P;
Deflection(n) = 1*DF(1);

end
10/2/23 7:50 PM C:\Users\DELL\Desktop\ASM-ASS...\ASS4A.m 3 of 3

plot(Deflection,Force)
10/2/23 7:58 PM C:\Users\DELL\Desktop\ASM-ASS...\ASS4B.m 1 of 3

% 4-B: Force (V) vs lateral displacement behavior if P=50 kN


clear all;
format short e;

for n=1:100
%ISMB500 Section (COLUMNS & BEAM)
A1 = 11074;A2 = 11074;A3 = 11074;
I1 = 452183000;I2 = 452183000;I3 = 452183000;
L1 =20000;L2 =10000;L3 =20000;
E1= 200000;E2= 200000;E3= 200000;
P = 50000;
V=(n-1)*1; % increase axial force in increments of 10kN.

% Linear stiffness matrix of column member


SM01 = [A1*E1/L1 0 0 -A1*E1/L1 0,
0 12*E1*I1/L1^3 6*E1*I1/L1^2 0 -12*E1*I1/L1^3,
0 6*E1*I1/L1^2 4*E1*I1/L1 0 -6*E1*I1/L1^2,
-A1*E1/L1 0 0 A1*E1/L1 0,
0 -12*E1*I1/L1^3 -6*E1*I1/L1^2 0 12*E1*I1/L1^3];

% Geometric stiffness matrix of a column


SG1 = P/L1*[0 0 0 0 0,
0 6/5 L1/10 0 -6/5,
0 L1/10 2*L1^2/15 0 -L1/10,
0 0 0 0 0 ,
0 -6/5 -L1/10 0 6/5];
SM1 = SM01-SG1;

SM02 = [A2*E2/L2 0 -A2*E2/L2 0,


0 12*E2*I2/L2^3 0 -12*E2*I2/L2^3,
-A2*E2/L2 0 A2*E2/L2 0,
0 -12*E2*I2/L2^3 0 12*E2*I2/L2^3];

SG2 = V/L2*[0 0 0 0,
0 6/5 0 -6/5,
0 0 0 0,
0 -6/5 0 6/5];
SM_B2 = SM02-SG2;

SM03 = [A3*E3/L3 0 -A3*E3/L3 0 0,


0 12*E3*I3/L3^3 0 -12*E3*I3/L3^3 6*E3*I3/L3^2,
-A3*E3/L3 0 A3*E3/L3 0 0,
0 -12*E3*I3/L3^3 0 12*E3*I3/L3^3 -6*E3*I3/L3^2,
0 6*E3*I3/L3^2 0 -6*E3*I3/L3^2 4*E3*I3/L3];

SG3 = 2*P/L3*[0 0 0 0 0,
0 6/5 0 -6/5 L3/10,
10/2/23 7:58 PM C:\Users\DELL\Desktop\ASM-ASS...\ASS4B.m 2 of 3

0 0 0 0 0,
0 -6/5 0 6/5 -L3/10,
0 L3/10 0 -L3/10 2*L3^2/15];

SM3 = SM03-SG3;

% Transformation matrix for a vertical upward column (member 1)


RT1=[0 1 0 0 0,
-1 0 0 0 0,
0 0 1 0 0,
0 0 0 0 1,
0 0 0 -1 0];

% Transformation matrix for a vertical downward column (member 3)


RT3=[0 -1 0 0 0,
1 0 0 0 0,
0 0 0 -1 0,
0 0 1 0 0,
0 0 0 0 1];

% Transformed stiffness matrix for first column (member 1)


SMS_C1 = RT1'*SM1*RT1;
% Transformed stiffness matrix for Beam (member 2)
SM_B2 = SM02-SG2;
% Transformed stiffness matrix for second column (member 3)
SMS_C3 = RT3'*SM3*RT3;

% no geomtric stiffness matrix required for the beam as there is negligible


% axial force on the beam

S_Total (1:10,1:10) = 0; % initialize the total stiffness matrix

% Assemble the global stiffness matrix


S_Total(1:5,1:5)=S_Total(1:5,1:5)+ SMS_C1;
S_Total(4:7,4:7)=S_Total(4:7,4:7)+ SM_B2;
S_Total(6:10,6:10)=S_Total(6:10,6:10)+SMS_C3;

%SWAY (FDOF = 4, 5, 6, 7, 8, 9)
SFF(1:4,1:4)= 0;
SFF(1:4,1:4)= S_Total(4:7,4:7);

AF = [V -2*P 0 -P]'; % force vector


DF=inv(SFF)*AF; % calculate displacements at FDOFs
10/2/23 7:58 PM C:\Users\DELL\Desktop\ASM-ASS...\ASS4B.m 3 of 3

Force(n) = V;
Deflection(n) = 1*DF(1);

end

plot(Deflection,Force)

You might also like