A Stabilized Finite Element Formulation of Non-Smooth Contact

The computational modeling of many engineering problems using the Finite Element method involves the modeling of two or more bodies that meet through an interface. The interface can be physical, as in multi-physics and contact problems, or purely numerical, as in the coupling of non-conforming meshes. The most critical part of the modeling process is to ensure geometric compatibility and a complete transfer of surface tractions between the different components at the connecting interfaces. Contact problems are a special family of interaction problems where the bodies on either side of the interface may separate freely or connect with each other, depending on the direction of motion. This type of behavior can be observed in complex civil, mechanical, bio-mechanical or aerospace structural components, and, on a smaller scale, in the interaction of different constituents in heterogeneous and composite materials and in the opening and closing of cracks in fracture mechanics. Popular contact modeling techniques rely on geometric projections to detect and resolve overlapping or mass interpenetration between two or more contacting bodies. Such approaches have been shown to have two major drawbacks: they are not suitable for contact at highly nonlinear surfaces and sharp corners where smooth normal projections are not feasible, and they fail to guarantee a complete and accurate transfer of pressure across the interface. This dissertation presents a novel formulation for the modeling of contact problems that possesses the ability to resolve complicated contact scenarios effectively, while being simpler to implement and more widely applicable than currently available methods. We show that the formulation boils down to a node-to-surface gap function that works effectively for non-smooth contact. The numerical implementation using the midpoint rule shows the need to guarantee the conservation of the total energy during impact, for which a Lagrange multiplier method is used. We propose a local enrichment of the interface and a simple stabilization procedure based on the discontinuous Galerkin method to guarantee an accurate transfer of the pressure field. The result is a robust interface formulation for contact problems and the coupling of non-conforming meshes.

Name thesis:degree_name

Ph.D.

Level thesis:degree_level

Dissertation

Discipline thesis:degree_discipline

Civil Engineering

Grantor

University of Illinois at Urbana-Champaign

Year dc:date

2014

Author dc:creator

• Haikal, Ghadir

Contributors dc:contributor

• Hjelmstad, Keith D.

• Engineering, Civil

Identifier

(UMI)AAI3362916

OAI identifier oai:identifier

oai:www.ideals.illinois.edu:2142/72183