USER’S MANUAL FOR FRAME ANALYSIS MATLAB PROGRAM
INPUT DATA:
40kN
60kN/m
400kN
2 2 3
3m
1 3
1
4
1m 3m
INPUT DATA: Units:
Length=(cm), load=(kN), Elasticity modulus=(kN/cm2), area=(cm2)
nm=3; number of members.
nj=4; number of joints.
e=20000; modulus of elasticity
m1=[1 2 3]; first end number of members;
m2=[2 3 4]; second end number of members
mt=[1 1 1]; member type
ar=[50 30 50]; areas of members;
at=[10000 6000 10000];moment of inertias of members
x=[0 0 400 400]; x-coordinates of joints;
y=[0 300 300 0]; y-coordinates of joints
jr=[0 0 0 1 1 1 1 1 1 0 0 0]; joint restraints
p=[0 400 0 0 0 0 0 0]; external joint load vector
mloads=1; if there are no loads on members mloads=0 else mloads =1
ndl=1; total # of distributed loads acting on various members in the frame
npl=1; total # of point loads acting of various members in the frame
nal=0; total # of axial loads acting of various members in the frame
dld=[2 -0.6 -0.6 0 400]; input data for distributed loads
pld =[2 -40 100];input data for point loads
ald = [0 0 0]; input data for axial loads
Distributed loads:
dld(i,1): member number
dld(i,2): qa
dld(i,3): qb
dld(i,4): bℓ
dld(i,5): c P
Point loads:
pld(i,1): member number
pld(i,2): P a
pld(i,3): a
P
Axial loads:
ald(i,1): member number a
ald(i,2): P
ald(i,3): a
OUTPUT:
No hinge connection at the ends of member 2 : mt=[1 1 1]
joint displacements
joint no x-disp y-disp rotation
1 0.000 0.0000e+00 -7.7636e-02
2 19.330 4.5000e-02 -3.8030e-02
3 19.181 -1.2900e-01 -3.0340e-02
4 0.000 0.0000e+00 -8.0734e-02
member forces
member end axial force shear force moment
1 1 -150.000 1.7602e+02 -2.9104e-11
1 2 150.000 -1.7602e+02 5.2807e+04
2 1 223.976 -1.5000e+02 -5.2807e+04
2 2 -223.976 4.3000e+02 -6.7193e+04
3 1 430.000 2.2398e+02 6.7193e+04
3 2 -430.000 -2.2398e+02 -2.9104e-11
60kN/m
400kN
Hinge at the first end of member 2: 2 2 3
3m
1 3
mt=[1 2 1]
1
4
1m 3m
joint displacements
joint no x-disp y-disp rotation
1 0.000 0.0000e+00 0.0000e+00
2 5.803 -2.3081e-02 -2.9013e-02
3 5.622 -6.0919e-02 -1.7150e-02
4 0.000 0.0000e+00 0.0000e+00
member forces
member end axial force shear force moment
1 1 76.938 1.2895e+02 3.8684e+04
1 2 -76.938 -1.2895e+02 1.4552e-11
2 1 271.054 7.6938e+01 0.0000e+00
2 2 -271.054 2.0306e+02 -2.9225e+04
3 1 203.062 2.7105e+02 2.9225e+04
3 2 -203.062 -2.7105e+02 5.2091e+04
60kN/m
400kN
Hinge at the second end of member 2:
2 2 3
mt=[1 3 1]
3m
1 3
1
4
1m 3m
joint displacements
joint no x-disp y-disp rotation
1 0.000 0.0000e+00 0.0000e+00
2 7.463 -3.4389e-02 -3.2008e-02
3 7.354 -4.9611e-02 -3.6769e-02
4 0.000 0.0000e+00 0.0000e+00
member forces
member end axial force shear force moment
1 1 114.629 2.3658e+02 5.6826e+04
1 2 -114.629 -2.3658e+02 1.4148e+04
2 1 163.419 1.1463e+02 -1.4148e+04
2 2 -163.419 1.6537e+02 0.0000e+00
3 1 165.371 1.6342e+02 0.0000e+00
3 2 -165.371 -1.6342e+02 4.9026e+04
Hinge at both ends of member 2:
40kN
60kN/m
mt = [1 4 1] 400kN
2 3
2
3m
1 3
1
4
1m 3m
joint displacements
joint no x-disp y-disp rotation
1 0.000 0.0000e+00 0.0000e+00
2 9.066 -7.2000e-02 -4.5331e-02
3 8.934 -7.2000e-02 -4.4669e-02
4 0.000 0.0000e+00 0.0000e+00
member forces
member end axial force shear force moment
1 1 240.000 2.0147e+02 6.0441e+04
1 2 -240.000 -2.0147e+02 0.0000e+00
2 1 198.529 2.4000e+02 0.0000e+00
2 2 -198.529 2.4000e+02 0.0000e+00
3 1 240.000 1.9853e+02 0.0000e+00
3 2 -240.000 -1.9853e+02 5.9559e+04
