Coursework Plates and Shells: Analysis and Computation (4D9)
Coursework Plates and Shells: Analysis and Computation (4D9)
Problem description
The objective of this coursework is to implement a finite element code for static and dynamic analysis of plate structures. There is a MATLAB finite element library CEAKIT
(Computational Engineering Analysis Kit) to be used for this coursework, which can obtained from:
http://www-g.eng.cam.ac.uk/csml/teaching/4d9/CEAKIT.tar.gz
The dowloaded file CEAKIT.tar.gz can be unpacked with tar -xzf CEAKIT.tar.gz,
which will create the CEAKIT directory containing several *.m files. Running the driver
function Ceakit in MATLAB will plot a deflected plate.
Tasks to be completed
1. Study the Ceakit.m driver and describe with few sentences the purpose of each
function called.
1
2. Implement a function which computes the load vector for uniform pressure loading.
3. Figure 1 shows the geometry of a plate to be analyzed. The boundaries of the plate
are simply supported and the plate thickness is t = 0.2. The Youngs modulus
of the material is E = 35000 and the Poissons ratio is = 0.3. The plate is
loaded by uniform pressure loading of p = 0.003. Study the convergence of the
maximal displacements for fully and selectively reduced integrated finite elements.
The meshes to be used should have 4 4, 8 8, 16 16 and 32 32 elements.
15.0
10.0
10.0
Figure 1:
4. Implement a function for computing the element mass matrix of a plate finite element and extend Ceakit.m for assembling the global mass matrix.
5. Implement the implicit Newmark time integration scheme.
6. Next, the dynamics of the plate due to sudden uniform loading should be studied.
The mass density of the material is = 2000. Apply a sudden uniform loading of
p = 0.003 and plot the evolution of the maximal displacements over time.
Note, it is not sufficient just to submit a working MATLAB implementation. It is crucial that you submit a report, which addresses item by item each
point of the previous list. Do not forget to include the requested plots and
the implemented equations.