A STUDY OF PENALTY FORMULATIONS USED IN THE NUMERICAL APPROXIMATION OF A RADIALLY SYMMETRIC ELASTICITY PROBLEM

Autor(es)

Aguiar, Adair R. ^[1] ; Fosdick, Roger L. ^[2] ; Sanchez, Jesus A. G. ^[1]

Afiliação do(s) autor(es)

^[1] Univ Sao Paulo, Sao Carlos Sch Engn, Dept Struct Engn, BR-13566590 Sao Carlos, SP - Brazil

^[2] Univ Minnesota, Minneapolis, MN 55455 - USA

• Tipo de documento: Artigo Científico
• Fonte: JOURNAL OF MECHANICS OF MATERIALS AND STRUCTURES; v. 3, n. 8, p. 1403-1427, OCT 2008.
• Citações Web of Science: 3

Resumo

We consider a class of two-dimensional problems in classical linear elasticity for which material overlapping occurs in the absence of singularities. Of course, material overlapping is not physically realistic, and one possible way to prevent it uses a constrained minimization theory. In this theory, a minimization problem consists of minimizing the total potential energy of a linear elastic body subject to the constraint that the deformation field must be locally invertible. Here, we use an interior and an exterior penalty formulation of the minimization problem together with both a standard finite element method and classical nonlinear programming techniques to compute the minimizers. We compare both formulations by solving a plane problem numerically in the context of the constrained minimization theory. The problem has a closed-form solution, which is used to validate the numerical results. This solution is regular everywhere, including the boundary. In particular, we show numerical results which indicate that, for a fixed finite element mesh, the sequences of numerical solutions obtained with both the interior and the exterior penalty formulations converge to the same limit function as the penalization is enforced. This limit function yields an approximate deformation field to the plane problem that is locally invertible at all points in the domain. As the mesh is refined, this field converges to the exact solution of the plane problem. (AU)

• Processo FAPESP: 07/03753-9 - Comportamento mecânico de sólidos na vizinhança de pontos singulares
• Beneficiário: Adair Roberto Aguiar
• Modalidade de apoio: Auxílio à Pesquisa - Regular