Finite element comparison for a geologically motivated benchmark

Geologic deformation in three dimensions can be modeled using finite element analysis. In choosing the elements used to solve a model it is important to consider the accuracy of the solution and the computational intensity. The results for models using six element types and six element side lengths are compared for the accuracy of the displacements calculated by the solution and the number of nodes required, as a proxy for computational intensity. Elements that allow higher order solutions are much more accurate than elements that only allow linear interpolation of the stresses and displacements between nodes; however the number of nodes required is five times greater. Free-form meshes do not significantly improve the performance of tetrahedra for the models tested, but could be accurate enough to model curved problem geometries. Comparisons for other models, such as a thrust fault, can be made using a twodimensional simplification of the three-dimensional problem. If three-dimensional comparisons are required it is important to choose a model that has an analytical solution.

Department dc:contributor.department

Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences.

Grantor dc:publisher

Massachusetts Institute of Technology

Year dc:date.issued

2003

Author dc:creator

• Peters, Teresa Baker, 1981-

Advisor dc:contributor.advisor

• Bradford H. Hager.

• Earth, Atmospheric, and Planetary Sciences.

Statement dc:rights

• MIT theses are protected by copyright. They may be viewed, downloaded, or printed from this source but further reproduction or distribution in any format is prohibited without written permission.

Licence dc:rights.uri

• http://dspace.mit.edu/handle/1721.1/7582

Language dc:language.iso

eng

Handle dc:identifier.uri

http://hdl.handle.net/1721.1/114317

OAI identifier oai:identifier

oai:dspace.mit.edu:1721.1/114317