A Robust Solution Procedure for Hyperelastic Solids with Large Boundary Deformation

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Shontz, Suzanne M., Vavasis, Stephen A.
Formato: Preprint
Publicado: 2006
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866917022356471808
author Shontz, Suzanne M.
Vavasis, Stephen A.
author_facet Shontz, Suzanne M.
Vavasis, Stephen A.
contents Compressible Mooney-Rivlin theory has been used to model hyperelastic solids, such as rubber and porous polymers, and more recently for the modeling of soft tissues for biomedical tissues, undergoing large elastic deformations. We propose a solution procedure for Lagrangian finite element discretization of a static nonlinear compressible Mooney-Rivlin hyperelastic solid. We consider the case in which the boundary condition is a large prescribed deformation, so that mesh tangling becomes an obstacle for straightforward algorithms. Our solution procedure involves a largely geometric procedure to untangle the mesh: solution of a sequence of linear systems to obtain initial guesses for interior nodal positions for which no element is inverted. After the mesh is untangled, we take Newton iterations to converge to a mechanical equilibrium. The Newton iterations are safeguarded by a line search similar to one used in optimization. Our computational results indicate that the algorithm is up to 70 times faster than a straightforward Newton continuation procedure and is also more robust (i.e., able to tolerate much larger deformations). For a few extremely large deformations, the deformed mesh could only be computed through the use of an expensive Newton continuation method while using a tight convergence tolerance and taking very small steps.
format Preprint
id arxiv_https___arxiv_org_abs_cs_0609001
institution arXiv
publishDate 2006
record_format arxiv
spellingShingle A Robust Solution Procedure for Hyperelastic Solids with Large Boundary Deformation
Shontz, Suzanne M.
Vavasis, Stephen A.
Numerical Analysis
Computational Engineering, Finance, and Science
G.1.0; G.1.5; G.1.8; J.2
Compressible Mooney-Rivlin theory has been used to model hyperelastic solids, such as rubber and porous polymers, and more recently for the modeling of soft tissues for biomedical tissues, undergoing large elastic deformations. We propose a solution procedure for Lagrangian finite element discretization of a static nonlinear compressible Mooney-Rivlin hyperelastic solid. We consider the case in which the boundary condition is a large prescribed deformation, so that mesh tangling becomes an obstacle for straightforward algorithms. Our solution procedure involves a largely geometric procedure to untangle the mesh: solution of a sequence of linear systems to obtain initial guesses for interior nodal positions for which no element is inverted. After the mesh is untangled, we take Newton iterations to converge to a mechanical equilibrium. The Newton iterations are safeguarded by a line search similar to one used in optimization. Our computational results indicate that the algorithm is up to 70 times faster than a straightforward Newton continuation procedure and is also more robust (i.e., able to tolerate much larger deformations). For a few extremely large deformations, the deformed mesh could only be computed through the use of an expensive Newton continuation method while using a tight convergence tolerance and taking very small steps.
title A Robust Solution Procedure for Hyperelastic Solids with Large Boundary Deformation
topic Numerical Analysis
Computational Engineering, Finance, and Science
G.1.0; G.1.5; G.1.8; J.2
url https://arxiv.org/abs/cs/0609001