Republic of the Philippines

CENTRAL LUZON STATE UNIVERSITY

College of Engineering

DEPARTMENT OF CIVIL ENGINEERING

CENGR 3140 LAB

(Numerical Solutions to CE Problem)

Full Name/s: JABONILLA, FLONIE R. ID Number: 22-1260

Laboratory Work No. 2

SOLVING SYSTEM OF LINEAR EQUATIONS

OBJECTIVE:

a. Recall basic concepts and formulas to create system of linear equations in solving
civil engineering problems.
b. Illustrate how to use MATLAB in solving civil engineering problems involving
systems of linear equations.

EQUIPMENT:

a. Reference books
b. Personal computer

PROCEDURES:

a. Choose your partner for this laboratory work. Your instructor will assign a specific
field/topic in civil engineering. Select one civil engineering problem within the
assigned topic that is solved using system of linear equations with minimum of four
equations (no repetition of problems).
b. Explain the concepts, theories, and/or laws involved in the problem and how did you
come up with the equations.
c. Solve the unknown values using one of the four numerical methods in solving system
of linear equations.
d. To verify your answers, solve the unknown values using MATLAB.
 STRUCTURAL THEORY NO.2
Problem:

Determine the moments at each joint of the frame shown in the figure.

A SLOPE DEFLECTION METHOD FOOR

BUILDING FRAMES SAMPLE PROBLEM
FROM THE BOOK STRUCTURAL
TH
ANALYSIS 8 EDITION BY R.C.
HIBBELER

R. C. HIBBELER
Independently Published, Sep 8, 2020 - Science - 704
pages
Structural Analysis, 8th, provides readers with a clear and
thorough presentation of the theory and application of
structural analysis as it applies to trusses, beams, and
frames. Emphasis is placed on teaching readers to both
model and analyze a structure. Procedures for Analysis,
Hibbeler's problem-solving methodologies, provides
readers with a logical, orderly method to follow when
applying theory.

Discussion:
A frame will sidesway, or be displaced to the side, when the structure itself or the
loading acting on it is nonsymmetric.

S O L U T I O N:
The frame will sway/drift to the left because the moment due to the applied load is
larger on the right side than the left. Therefore, the yielded direction of the unknown
horizontal displacements ∆ of the joints should be negative since the column will rotate in a
counterclockwise direction.

Numerical Methods Calculations:

Method Used: Slope Deflection Method

Slope-Deflection Equations (per Member)

Joint A3
 1 1 1 3
MA3B3= EI θ A 3+ EI θ B 3−32 MA3A2= EI θ A 3+ EI θ A 2− EI ∆2
2 4 2 8

Joint B3
1 1 1 3
MB3A3= EI θ A 3+ EI θ B 3 + 48 MB3B2= EI θ B 3 + EI θB 2− EI ∆ 2
4 2 2 8

JointA2
1 3 1 1
MA2A3= EI θ A 3+ EI θ A 2− EI ∆2 MA2B2= EI θ A 2+ EI θ B 2−64
2 8 2 4

4 5
MA2A1= EI θ A 2− EI ∆ 1
5 25

JointB2
1 1 1 3
MB2A2= EI θ A 2+ EI θ B 2 +96 MB2B3= EI θ B 3+ EI θB 2− EI ∆ 2
4 2 2 8

2 1
MB2B1= EI θ B 2− EI ∆1
3 6

Joint A1
2 6
MA1A2= EI θ A 2− EI ∆1
5 25

Joint B1
1 1
MB1B2= EI θ B 2− EI ∆1
3 6

Moment Equilibrium Conditions (per Joint/Node)

ΣMA3 = 0; MA3B3 + MA3A2 = 0

1 3 1 3
EI θ A 3+ EI θB 3 + EI θ A 2− EI ∆2=32 eq.1
4 2 2 8

ΣMB3 = 0; MB3A3 + MB3B2 = 0

3 1 1 3
EI θ A 3+ EI θB 3 + EI θ B 2− EI ∆2=−48 eq.2
2 4 2 8

ΣMA2 = 0; MA2A3 + MA2B2 + MA2A1 = 0

1 23 1 3 6
EI θ A 3+ EI θ A 2 + EI θB 2 − EI ∆ 2− EI ∆ 1=64 eq.3
2 10 4 8 25

ΣMB2 = 0; MB2A2 + MB2B3 + MB2B1 = 0

1 1 13 3 1
EI θ B 3+ EI θ A 2 + EI θB 2− EI ∆ 2− EI ∆1=−96 eq.4
2 4 6 8 6
Force Equilibrium Conditions (per Floor/Level)

ΣV2nd Floor = 0; VA2A3 + VB2B3 = 0

3 3 3 3 3
EI θ A 3+ EI θB 3+ EI θ A 2 + EI θ B 2− EI ∆2=0 eq.5
2 2 2 2 2

ΣV1st Floor = 0; VA1A2 + VB1B2 = 0

36 3 41
EI θ A 2+ 5 EI θB 2− EI ∆1=0 eq.6
5 75

Numerical Methods Calculations:

Method Used: CRAMER’S RULE

1 3 1 3
EI θ A 3+ EI θB 3 + EI θ A 2− EI ∆2=32 eq.1
4 2 2 8

3 1 1 3
EI θ A 3+ EI θB 3 + EI θ B 2− EI ∆2=−48 eq.2
2 4 2 8

1 23 1 3 6
EI θ A 3+ EI θ A 2 + EI θ B 2− EI ∆2− EI ∆1=64 eq.3
2 10 4 8 25

1 1 13 3 1
EI θ B 3+ EI θ A 2 + EI θ B 2− EI ∆2− EI ∆1=−96 eq.4
2 4 6 8 6

3 3 3 3 3
EI θ A 3+ EI θB 3+ EI θ A 2 + EI θ B 2− EI ∆2=0 eq.5
2 2 2 2 2

36 341
EI θ A 2 +5 EI θB 2 − EI ∆1=0 eq.6
5 75

ѲA3 ѲB3 ѲA2 ѲB2 ∆2 ∆1 C

3/2 ¼ ½ 0 -3/8 0 32
¼ 3/2 0 ½ -3/8 0 -48
½ 0 23/10 ¼ -3/8 -6/25 64
0 ½ ¼ 13/6 -3/8 -1/6 -96
3/2 3/2 3/2 3/2 -3/2 0 0
0 0 36/5 5 0 -341/75 0
D= 16.9280

ѲA3 --- C; D= 99.2435

ѲB3 --- C; D= -496.9038
ѲA2 --- C; D= 355.4358
ѲB2 --- C; D= -862.6825
 ∆2---C; D=226.2267
∆1---C; D= 12.8612

THEREFORE:

EI θA3= 5. 862668698141

EI θB3= −29. 3538876568219

EI θA2= 20. 9968666389349

EI θB2= −50. 9617424736569

EI ∆2= −53. 4560947934025

EI ∆1= −22. 7928018844472

MATLAB Calculations:
Method Used: LEFT DIVISION

a=[3/2 1/4 1/2 0 -3/8 0; 1/4 3/2 0 1/2 -3/8 0]

a=

1.5000 0.2500 0.5000 0 -0.3750 0

0.2500 1.5000 0 0.5000 -0.3750 0

a=[3/2 1/4 1/2 0 -3/8 0; 1/4 3/2 0 1/2 -3/8 0]

a=

1.5000 0.2500 0.5000 0 -0.3750 0

0.2500 1.5000 0 0.5000 -0.3750 0

a=[3/2 1/4 1/2 0 -3/8 0; 1/4 3/2 0 1/2 -3/8 0; 1/2 0 23/10 1/4 -3/8 -6/25; 0 1/2 1/4 13/6 -3/8 -1/6; 3/2 3/2
3/2 3/2 -3/2 0; 0 0 36/5 5 0 -341/75]

a=

1.5000 0.2500 0.5000 0 -0.3750 0

0.2500 1.5000 0 0.5000 -0.3750 0
0.5000 0 2.3000 0.2500 -0.3750 -0.2400
0 0.5000 0.2500 2.1667 -0.3750 -0.1667
1.5000 1.5000 1.5000 1.5000 -1.5000 0
0 0 7.2000 5.0000 0 -4.5467
 b=[32; -48; 64; -96; 0; 0]

b=

32
-48
64
-96
0
0

x=a\b

x=

5.8627
-29.3539
20.9969
-50.9617
-53.4561
-22.7928

CONCLUSION:
As we can see, MATLAB also can solve things like this. Like we used to think
that those problems are only can be used by manual computation only, but in this
laboratory, it shows that we can also solved those problems by numerical solutions
using softwares like MATLAB, GTU OCTAVE and many more softwares.