Scribd
Upload
54 views4 pages

Ervin

Uploaded by

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

Ervin

Uploaded by

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

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: ERVIN JOHNSON I. ID Number: 21-1457


LAGAZO

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.
RESULTS: (use the format below)

(DISCIPLINE/BRANCH OF CE)
Problem:
A simple truss bridge is supported by two pin supports at points A and B. A vertical
load of 20 kN is applied at joint C. Determine the internal forces in members CD, DE,
and EF.
Image:

simple truss bridge with points A, B, C, D, E, and F

(Include a screenshot or picture of the source of your civil engineering problem.


Properly cite your source/reference in APA format.)

Discussion:
Solution:
1. Determine Reactions:
o Calculate the vertical reactions at supports A and B.
2. Cut Section:
o Pass a cutting plane through members CD, DE, and EF.
3. Analyze Section:
o Draw the free body diagram of the section to the left or right of the cut.
o Apply the equations of equilibrium (sum of forces in x and y directions,
and sum of moments) to solve for the unknown internal forces in
members CD, DE, and EF.
Note:
 The direction of the internal forces (tension or compression) can be assumed.
If the calculated value is positive, the member is in tension; if negative, it is in
compression.
 The choice of which side of the cut to analyze depends on which side
provides the simplest equations of equilibrium.
MATLAB Approach:
1. We define the truss properties like member lengths, applied load, and
reactions at supports A and B (calculated using equilibrium for the whole
truss).
2. We choose a section to analyze (here, to the left of CD-DE-EF).
3. We define vectors for the internal forces in CD and DE.
4. We create a free body diagram of the chosen section, showing the reactions,
applied load, and unknown internal forces.
5. We apply the equations of equilibrium:
o Sum of forces in y-direction to solve for F_DE(2).
o Moment balance about point D to solve for F_DE(1).
6. The code displays the calculated internal forces for CD and DE.

% Define truss properties


L_CD = 4; % Length of member CD (m)
L_DE = 3; % Length of member DE (m)
L_EF = 4; % Length of member EF (m)
P = 20e3; % Applied load at C (N)

% Since supports A and B are pins, reaction forces will be vertical


only
Ry_A = P/2; % Vertical reaction at A (N)
Ry_B = P/2; % Vertical reaction at B (N)

% Choose section to analyze (let's analyze to the left of cut CD-DE-


EF)
F_CD = zeros(2, 1); % Internal force in CD (tension positive)
F_DE = zeros(2, 1); % Internal force in DE (tension positive)

% Free body diagram of the section


figure;
hold on;
arrow([0 0],[0 Ry_A],'g','MaxHeadSize',10); % Reaction A

Unrecognized function or variable 'arrow'.

arrow([L_CD 0],[0 -P],'r','MaxHeadSize',10); % Applied load


arrow(F_CD, 'b', 'MaxHeadSize', 10); % Internal force CD (unknown)
arrow(F_DE, 'b', 'MaxHeadSize', 10); % Internal force DE (unknown)
text(L_CD/2, -P/2, '20 kN', 'HorizontalAlignment', 'center');
text(L_CD/3, Ry_A/2, sprintf('%.2f kN', Ry_A/1000),
'HorizontalAlignment', 'center');
axis([-1 L_CD+1 -P/2 Ry_A+P/5]);
xlabel('x (m)');
ylabel('y (N)');
title('Free Body Diagram (Section Left of Cut)');
hold off;

% Solve for internal forces using equilibrium equations


% Sum of forces in y-direction
F_DE(2) = Ry_A - P;

% Moment balance about point D


moment_DE = (F_CD(1) * L_CD) - (P * (L_CD - L_DE));
F_DE(1) = moment_DE / L_DE;

% Display results
disp('Internal forces (N):');
disp(['CD: ', num2str(F_CD)]);
disp(['DE: ', num2str(F_DE)]);

You might also like