Advanced search
Start date
Betweenand
Related content
(Reference retrieved automatically from Web of Science through information on FAPESP grant and its corresponding number as mentioned in the publication by the authors.)

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

Full text
Author(s):
Aguiar, Adair R. [1] ; Fosdick, Roger L. [2] ; Sanchez, Jesus A. G. [1]
Total Authors: 3
Affiliation:
[1] Univ Sao Paulo, Sao Carlos Sch Engn, Dept Struct Engn, BR-13566590 Sao Carlos, SP - Brazil
[2] Univ Minnesota, Minneapolis, MN 55455 - USA
Total Affiliations: 2
Document type: Journal article
Source: JOURNAL OF MECHANICS OF MATERIALS AND STRUCTURES; v. 3, n. 8, p. 1403-1427, OCT 2008.
Web of Science Citations: 3
Abstract

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)

FAPESP's process: 07/03753-9 - Mechanical behavior of solids in the vicinity of singular points
Grantee:Adair Roberto Aguiar
Support Opportunities: Regular Research Grants