[go: up one dir, main page]
Scribd
Upload
76 views14 pages

This Supplementary Lab Sheet Must Be Read With The Lab Manual - Students Are Required To Write Their Scripts, Plotted Graph and

This document contains instructions and problems for a supplementary lab sheet for CES 513/512. Problem 1 provides the equation and parameters for a linearly increasing load on a beam. It asks the student to (a) use MATLAB to graph the bending moment and estimate the distance(s) where it equals 2e6 kNcm, and (b) use the bisection method in MATLAB to iteratively estimate this distance(s) within 0.5% accuracy. It then asks the student to (c) determine the theoretical solution to 6 decimal places and compute the absolute true error compared to the estimated solution.

Uploaded by

Hazim Azri
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
76 views14 pages

This Supplementary Lab Sheet Must Be Read With The Lab Manual - Students Are Required To Write Their Scripts, Plotted Graph and

This document contains instructions and problems for a supplementary lab sheet for CES 513/512. Problem 1 provides the equation and parameters for a linearly increasing load on a beam. It asks the student to (a) use MATLAB to graph the bending moment and estimate the distance(s) where it equals 2e6 kNcm, and (b) use the bisection method in MATLAB to iteratively estimate this distance(s) within 0.5% accuracy. It then asks the student to (c) determine the theoretical solution to 6 decimal places and compute the absolute true error compared to the estimated solution.

Uploaded by

Hazim Azri
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/ 14

NAME OF STUDENT: MOHAMMAD ZAKWAN BIN ZULKIFLI

GROUP: PEC2216B1
DATE: 3-10-2017

CES 513/512 SUPPLEMENTARY LAB SHEET DOCUMENT (LAB 3)

INSTRUCTIONS:

• THIS SUPPLEMENTARY LAB SHEET MUST BE READ WITH THE LAB MANUAL
BEFORE THE LAB ACTIVITIES
• STUDENTS ARE REQUIRED TO WRITE THEIR SCRIPTS, PLOTTED GRAPH AND
SOLUTIONS IN THIS LAB SHEET. CONVERT YOUR DOCUMENTS TO PDF
FORMAT BEFORE EMAIL IT TO YOUR CORRESPONDING LECTURER)

PROBLEM 1

Figure Q1 shows a linearly increasing load of 1.75 kN/cm acting on a beam with
span of 450 cm. The equation of the resulting beam’s deflection, y, is given in
Equation (1) as follows;

W=1.75 kN/cm

L=450 cm

Figure Q1

Equation (1)

It is noted that EI is the beam flexural rigidity and it is given as E = 50000 kN/cm 2 and
I = 30000 cm4 respectively.
a) Using graph plotted by MATLAB, estimate the distance(s) on the beam, x, so
that the beam bending moment, M, is equals to 2 × 106𝑘𝑁𝑐𝑚. Given that 𝑀 =

W = 1.75; L = 450;
fplot(@(x) 2e6-(W/120)*(12*L^2*x - 20*x^3),[0,450],'k')
hold on
fplot(@(x) 0,[0,450],'r')
b) Using the Bisection Method algorithm learned in the lecture session, write a
MATLAB script (complete with suitable comment lines) on estimating the
distance(s) of x giving bending moment, M, equals to 2 × 106𝑘𝑁𝑐𝑚. . Report
the solution according the following format and the absolute relative error must
be less than 0.5%. Select your own lower and upper limit of the guessed root.

Iteration xL xU xM f(xm) epsilon

eps_s =0.001 eps_a

=200

L=450 ; W=1.75

XL=305 XU=350 m=@(x)2e6-

(W/120)*(12*L^2*x-20*x^3) mxu= m(XU)

mxl= m(XL) sign_changes = m(XU)*m(XL)

if sign_changes <0

disp('A root lies inside [xl,xu]') xr=(XL+XU)/2;

else

disp('No root in [xl.xu]') disp('Then

script will ended')

return

end

i=1; while

(eps_a>eps_s)

xrold=xr;

xr=(XL+XU)/2;

mxu=m(XU); mxr=m(xr);

sign_changes2=mxu*mxr;

if (sign_changes2<0)
XL=xr;xrnew=(XL+XU)/2;

else

XU=xr;

xrnew=(XL+XU)/2;

end eps_a=abs((xrnew-xrold)/xrnew)*100;

i=1+i; disp([i XL XU xr m(xr)

eps_a]) end

solution:

eps_s =

1.0000e-03

eps_a =

200

W=

1.7500
XL =

305

XU =

350

m=

@(x)2e6-
(W/120)*(12*L^2*x-
20*x^3)

mxu =

2.1021e+06
mxl =

-5.3309e+05

sign_changes =

-1.1206e+12

A root lies inside [xl,xu]

1.0e+05 *

Columns 1 through 4

0.0000 0.0031
0.0033 0.0033

Columns 5 through 6
6.3943 0.0000

1.0e+04 *

Columns 1 through 4

0.0003 0.0305
0.0316 0.0316

Columns 5 through 6

1.8146 0.0005

1.0e+05 *

Columns 1 through 4

0.0000 0.0031
0.0032 0.0031
Columns 5 through 6

-2.6607 0.0000

1.0e+05 *

Columns 1 through 4

0.0001 0.0031
0.0032 0.0031

Columns 5 through 6

-1.2613 0.0000

1.0e+04 *

Columns 1 through 4

0.0006 0.0315
0.0316 0.0315
Columns 5 through 6

-5.4538 0.0001

1.0e+04 *

Columns 1 through 4

0.0007 0.0316
0.0316 0.0316

Columns 5 through 6

-1.8332 0.0000

Columns 1 through 4

8.0000 315.8984
316.2500 315.8984
Columns 5 through 6

-127.2823 0.1668

1.0e+03 *

Columns 1 through 4

0.0090 0.3159
0.3161 0.3161

Columns 5 through 6

9.0009 0.0000

1.0e+03 *

Columns 1 through 4

0.0100 0.3159
0.3160 0.3160
Columns 5 through 6

4.4346 0.0000

1.0e+03 *

Columns 1 through 4

0.0110 0.3159
0.3159 0.3159

Columns 5 through 6

2.1531 0.0000

1.0e+03 *

Columns 1 through 4

0.0120 0.3159
0.3159 0.3159

Columns 5 through 6

1.0128 0.0000

Columns 1 through 4

13.0000 315.8984
315.9094 315.9094

Columns 5 through 6

442.7255 0.0052

Columns 1 through 4

14.0000 315.8984
315.9039 315.9039

Columns 5 through 6
157.7132 0.0026

Columns 1 through 4

15.0000 315.8984
315.9012 315.9012

Columns 5 through 6

15.2134 0.0013

Columns 1 through 4

16.0000 315.8998
315.9012 315.8998

Columns 5 through 6

-56.0350 0.0002
c) Determine analytically the theoretical solution of x correct to 6 decimal points.
Using MATLAB, compute the absolute true error of the estimated x from b).

m = @(x) 2e6-(1.75/120)*(12*450^2*x - 20*x^3);

format long g solution1 = fzero(m,60)

solution2 = fzero(m,315)

Solution:

solution1 = 58.0471711584301
solution2 = 315.900890850594

You might also like