CES Computational Analysis Methods for Engineers Assignment 1:

CES 512 / Assignment 1 2015 CES 512 - Computational Analysis Methods for Engineers Assignment 1: This assignment consist a civil engineering proble...
Author: Rhoda Patterson
0 downloads 2 Views 307KB Size
CES 512 / Assignment 1

2015

CES 512 - Computational Analysis Methods for Engineers Assignment 1:

This assignment consist a civil engineering problems that need to be solved by using knowledge of roots of equation, Numerical solution of simultaneous equations and Optimization. Assignment outcome: At the end of this assignment, student should be able to: i. usebracketing and open methods to solve civil engineering problemsmanually and using software. ii. use methods learn in Numerical solution of simultaneous equations to solve civil engineering problems manually and using software. iii. Use optimization to solve the constrained problem. Tasks: Question 1

Figure Q1(a)

Figure Q1(b)

CES512-Assignment 1-MAM-270815

CES 512 / Assignment 1

2015

Figure 1(a) shows the overhanging beam loaded by uniform distribution load, W and Point load, P. It is supported by pinned and roller at A and B respectively. Figure 1(b) represents the bending moment diagram of the beam which developed based on the load cases in Table Q1. Table Q1 Load Case 1 2

W(kN/m) 1 2.5

P (kN) 1 4

a) By considering distance, x from node A, derive the equation of reaction (V A & V B ) and moment, M AB in term of W, P, and xfor beam span AB.(CO1-PO2, C5) (2 marks) b) Based on derived equation in (a), derive the function, f(x)=0 to correlate the load cases in Table Q1 for the intercept point of bending moment as depicted in Figure Q1(b). From the derived function f(x), is it possible to determine the exact solution for the function? (CO1-PO2, C5) (3 marks) c) By using the function, f(x) established in (b), Use the bisection method to determine for the position, x at where the bending moment in both cases is intercept. Carry out the iteration up to approximate percent relative error πœ€πœ€π‘Žπ‘Ž is less than πœ€πœ€π‘ π‘  = 0.1%. Calculate the true percentage error, πœ€πœ€π‘‘π‘‘ as well. Tabulate your answer. [Hint: the initial guesses of lower (x l ) and upper (x u ) of x need to be determined based on rough estimation in bending moment diagram as plotted in Fig 1(b).](CO1-PO2, C4) (5 marks) d) For the same condition as in (c), apply the Newton-Raphson method to solve the problem as mentioned in (c). Initial guesses, x o must be assumed based on rough estimation the plotted graph in bending moment diagram as plotted in Fig 1(b).]Tabulate your answer.(CO1-PO2, C4) (3 marks) e) Compare your answers obtained in (c) and (d) by using the built-in-function fzero and writing programming (script M-file) in MAT-LAB.(CO2-PO4, C2) (2 marks) f) Comment and discuss for the following point: (CO2-PO4, C2) i. ii.

What is your opinion to improve the accuracy of the final answer and shorten the iteration by using both methods in (c) & (d) towards the exact solution in (b). What is the importance of graphical method in finding solution using methods in (c) and (d) (5 marks)

CES512-Assignment 1-MAM-270815

CES 512 / Assignment 1

2015

Question 2

Figure Q2 Member 2 of the truss shown in Figure Q2 is subjected to an increase in temperature of 83.3oC. This temperature consequent the axial displacements in entire members in the truss and the system of linear equations to determine the elongation of the members are given by: 0.135 D 1 + 0.035 D 2 = - 0.000689325 0.035 D 1 + 0.135 D 2 - 0.1 D 4 = - 0.000689325 0.135 D 3 – 0.035 D 4 + 0.035 D 5 = 0 -0.1 D 2 – 0.035 D 3 + 0.135 D 4 – 0.035 D 5 = 0 0.035 D 3 – 0.035 D 4 + 0.135 D 5 = 0 a) Manually solve for the member displacement for the truss using Gauss Elimination (Consider up to 5 decimal places) (CO1:PO2) – C3 : 10m (10 marks) b) By using the built in function in the MAT-LAB, compare your results with the manualas calculated in (a). (CO2:PO4) – C2 : 10m (5 marks)

CES512-Assignment 1-MAM-270815

CES 512 / Assignment 1

2015

Question 3 P kN w kN/m

E, I L/2

L/2

Figure Q3 Figure Q3 shows the simply supported beam imposed by point load, P and UDL, w. Thecorresponding total displacement at the mid-span of beam can be determined by following equation:

𝛿𝛿𝑑𝑑 = 𝛿𝛿𝑀𝑀 + 𝛿𝛿𝑝𝑝

5𝑀𝑀𝐿𝐿4 𝑃𝑃𝐿𝐿3 𝛿𝛿𝑑𝑑 = + 384𝐸𝐸𝐸𝐸 48𝐸𝐸𝐸𝐸

Where Ξ΄ w = displacement due to UDL, Ξ΄ p = displacement due to point load

If the total beams displacement, Ξ΄ t is constrained byconditions as shown in Table Q3: Table Q3 Constrained condition

Capacity

Point Load, P (kN)

40

Uniform distribution load UDL, w (kNm)

10

Shear force at support, V (kN)

30

Bending moment at mid-span, M (kNm)

60

Determine the value of point load, P and UDL, w at maximum deflection, Ξ΄ t of the beam to satisfy the condition in Table Q3 by: a) Formulate the linear programming problem.(CO1:PO2) – C5 : 5m CES512-Assignment 1-MAM-270815

b) c) d) e)

CES 512 / Assignment 1

2015

Solve the problem byusing Graphical method.(CO1:PO2) – C3 : 6m Solve the problem with an excel solver. (CO1:PO2) – C3 : 4m Discuss on the overall finding from the result analysis obtained.(CO2:PO4) – C2 : 5m Determine the percentage contribution of the Ξ΄ w and Ξ΄ p towards the total displacement of the beam,Ξ΄ t and conclude your finding. (CO2:PO4) – C6 : 5m f) Compare the total displacement, Ξ΄ t obtained in the part (b & c) with the allowable displacement of the beam,Ξ΄ a =3 x 10-3 m. If the displacement, Ξ΄ t is exceeded the allowable displacement, what is your suggestion to satisfy the requirement without violent the loading P and w.(CO2:PO4) - C5 : 5m (given L= 5m, E= 200GPa,andI = 1.07 x 10-3 m3)

CES512-Assignment 1-MAM-270815