ENS 5147 - Advanced Structural Analysis and Design
Part 1 - Civil Engineering Numerical Analysis
Submission: before 4.45pm, Friday 21st April 2017. No late submission accepted.
H(m) = 0.3(m).
1. [40 %] Treating the bar illustrated above as being a one dimensional axial bar where plane sections remain plane and parallel, analyse the bar using the Galerkin finite element method. The coding solution for linear h-refinement and its manual solution for a limited number of linear elements are on Blackboard for your reference.
Consider h‐refinement using three‐noded quadratic elements, produce the manual solution for when the domain includes of 2 equal spacing quadratic elements. Generate the outcomes to produce:
a. Nodal displacements
b. Displacement at the left quarter of each element.
c. Normal stress at the left quarter of each element.
2. [10 %] The bar illustrated above is solved as a two dimensional plane stress problem. This task allows you to start with FEM programming at very beginning level. You will be given full coding solution and asked to extract a number of results given by running the problem FEM program.
Modify the code to yield 2 equal spacing 4‐noded quadrilateral elements along x-axis. Generate the outcomes to produce:
a. Each element coding loop: indexArray, K, f.
b. Global K, f, u.
c. Nodal stresses.
3. [50 %] Extend the two‐noded linear Matlab program provided to produce the corresponding solution for three‐noded quadratic elements. You need to amend all sub function m-files, and the main m-file. Debug your code to confirm if it works properly before put all files in a zip folder for USB submission.