Description

Coding Assignment 3
Consider the 3D elastostatics problem. Find u such that
PDE σij,j + fi = 0 in Ω
Constitutive relation
Kinematic relation
Neumann b.c. σijnj = hi on ∂Ωhi
Dirichlet b.c. ui = ugi on ∂Ωui
Consider a three-dimensional domain defined by x1 = [0,1] m; x2 = [0,1] m; x3 = [0,1] m (i.e. the unit cube). Use E = 2.0e11 Pa and ν = 0.3. Assume traction hi = 0 N.m−2 on all surfaces where no other conditions are specified. Use linear basis functions and a 10 x 10 x 10 element mesh for submission.
Apply the following boundary conditions:
h1 = h2 = 0,h3 = 1.0e9 ∗ x1 Pa on the face x3 = 1 m and u1 = u2 = u3 = 0 m on the face x3 = 0 m.
Coding Instructions: You will receive the following files:
• main3.cc (the source file)
• FEM3.h (the template header file)
• writeSolutions.h
• CMakeLists.txt
You will also receive the .vtk solution file for a 5 x 5 x 5 mesh. As before, most of your coding will be done in FEM3.h. Do not modify any function names or the names of any class data structure. The only parts of main3.cc that should be modified are the function inputs that define the number of elements in the mesh. Your FEM3.h file must run with the given main3.cc file. Nothing should be changed in writeSolutions.h or CMakeLists.txt.
Submission Instructions: You should submit (through the Coursera website) a .zip file (name it CA3.zip) containing your FEM3.h file and the .h5 solution file. You will have one .h5 solution file. The .h5 file will contain the solution vector. Do not submit any files other than those listed here.
1