Hyper-elastic formulations are usually employed to represent the constitutive response of biological tissues. Particularly these formulations are often encountered in models assessing the arterial wall behavior. A great number of distinct hyper-elastic laws can be found in the specialized literature devised to represent, at different extents, the mechanical aspects of the arterial tissues. However, no one can be, a priori, considered more suitable than the others. The choice depends on the specific applications and on the type of mechanical response the corresponding models are intended to account for. Consequently, it is convenient to take at hand the possibility of implementing and rapid prototyping different constitutive laws in an easy and reliable manner. Therefore, in this work we describe how to implement a generic Finite Element framework capable to accommodate practically any hyper-elastic material law. This is carried out using a spatial variational formulation for the momentum equation which is linearized by means of a Newton-Raphson scheme. The iterative algorithm is such that for a given load, the equilibrium is reached in the deformed spatial configuration. The main feature of our approach is based on the evaluation of the second order stress tensor and of the fourth order constitutive tangent tensor using finite differences. That is, given a strain energy potential we compute, by means of a second order finite difference centered scheme, the first (stress) and second (tangent matrix) derivatives. In this way, a generic computational implementation in the context of Finite Elements is achieved, making possible to change the material behavior just changing the procedure that evaluates the elastic function and reusing the entire numerical element structure. The developments are carried out for cuasi-incompressible materials, and some implementation issues are presented and discussed. The method is validated for two common constitutive laws of the arterial tissues including the mechanical response of the arterial wall considered as a fiber-reinforced multilayer material.
[1]
Pablo J. Blanco,et al.
Assessing the influence of heart rate in local hemodynamics through coupled 3D‐1D‐0D models
,
2010
.
[2]
R. Ogden,et al.
A New Constitutive Framework for Arterial Wall Mechanics and a Comparative Study of Material Models
,
2000
.
[3]
N. Stergiopulos,et al.
Residual strain effects on the stress field in a thick wall finite element model of the human carotid bifurcation.
,
1996,
Journal of biomechanics.
[4]
M. Gurtin,et al.
An introduction to continuum mechanics
,
1981
.
[5]
Jay D. Humphrey,et al.
Structure, Mechanical Properties, and Mechanics of Intracranial Saccular Aneurysms
,
2000
.
[6]
Kozaburo Hayashi,et al.
A strain energy function for arteries accounting for wall composition and structure.
,
2004,
Journal of biomechanics.
[7]
Pablo J. Blanco,et al.
Multidimensional modelling for the carotid artery blood flow
,
2006
.
[8]
P. Blanco,et al.
A unified variational approach for coupling 3D-1D models and its blood flow applications
,
2007
.
[9]
P. Blanco,et al.
On the potentialities of 3D-1D coupled models in hemodynamics simulations.
,
2009,
Journal of biomechanics.