Modelling and convergence in arterial wall simulations using a parallel FETI solution strategy

Arterial walls are characterised by nearly incompressible, anisotropic, hyperelastic material behaviour. Several polyconvex material functions representing such materials are considered and adjusted to experimental data. For all of these functions and for different parameter sets numerical simulations using a three-dimensional model of a diseased artery are performed. A finite element tearing and interconnecting-dual primal domain decomposition algorithm is used to solve the linearised systems of equations. The numerical performance of the different models is discussed with respect to convergence of the linear and nonlinear solvers.     

Incorporating tissue anisotropy and heterogeneity in finite element models of trabecular bone altered predicted local stress distributions

Trabecular bone is composed of organized mineralized collagen fibrils, which results in heterogeneous and anisotropic mechanical properties at the tissue level. Recently, biomechanical models computing stresses and strains in trabecular bone have indicated a significant effect of tissue heterogeneity on predicted stresses and strains. However, the effect of the tissue-level mechanical anisotropy on the trabecular bone biomechanical response is unknown. Here, a computational method was established to automatically impose physiologically relevant orientation inherent in trabecular bone tissue on a trabecular bone microscale finite element model. Spatially varying tissue-level anisotropic elastic properties were then applied according to the bone mineral density and the local tissue orientation. The model was used to test the hypothesis that anisotropy in both homogeneous and heterogeneous models alters the predicted distribution of stress invariants. Linear elastic finite element computations were performed on a 3 mm cube model isolated from a microcomputed tomography scan of human trabecular bone from the distal femur. Hydrostatic stress and von Mises equivalent stress were recorded at every element, and the distributions of these values were analyzed. Anisotropy reduced the range of hydrostatic stress in both tension and compression more strongly than the associated increase in von Mises equivalent stress. The effect of anisotropy was independent of the spatial redistribution high compressive stresses due to tissue elastic heterogeneity. Tissue anisotropy and heterogeneity are likely important mechanisms to protect bone from failure and should be included for stress analyses in trabecular bone.     

Mechanical stresses in abdominal aortic aneurysms: influence of diameter, asymmetry, and material anisotropy

Biomechanical studies suggest that one determinant of abdominal aortic aneurysm (AAA) rupture is related to the stress in the wall. In this regard, a reliable and accurate stress analysis of an in vivo AAA requires a suitable 3D constitutive model. To date, stress analysis conducted on AAA is mainly driven by isotropic tissue models. However, recent biaxial tensile tests performed on AAA tissue samples demonstrate the anisotropic nature of this tissue. The purpose of this work is to study the influence of geometry and material anisotropy on the magnitude and distribution of the peak wall stress in AAAs. Three-dimensional computer models of symmetric and asymmetric AAAs were generated in which the maximum diameter and length of the aneurysm were individually controlled. A five parameter exponential type structural strain-energy function was used to model the anisotropic behavior of the AAA tissue. The anisotropy is determined by the orientation of the collagen fibers (one parameter of the model). The results suggest that shorter aneurysms are more critical when asymmetries are present. They show a strong influence of the material anisotropy on the magnitude and distribution of the peak stress. Results confirm that the relative aneurysm length and the degree of aneurysmal asymmetry should be considered in a rupture risk decision criterion for AAAs.     

A Gauss-Kronrod-Trapezoidal integration scheme for modeling biological tissues with continuous fiber distributions

Fibrous biological tissues may be modeled using a continuous fiber distribution (CFD) to capture tension–compression nonlinearity, anisotropic fiber distributions, and load-induced anisotropy. The CFD framework requires spherical integration of weighted individual fiber responses, with fibers contributing to the stress response only when they are in tension. The common method for performing this integration employs the discretization of the unit sphere into a polyhedron with nearly uniform triangular faces (finite element integration or FEI scheme). Although FEI has proven to be more accurate and efficient than integration using spherical coordinates, it presents three major drawbacks: First, the number of elements on the unit sphere needed to achieve satisfactory accuracy becomes a significant computational cost in a finite element (FE) analysis. Second, fibers may not be in tension in some regions on the unit sphere, where the integration becomes a waste. Third, if tensed fiber bundles span a small region compared to the area of the elements on the sphere, a significant discretization error arises. This study presents an integration scheme specialized to the CFD framework, which significantly mitigates the first drawback of the FEI scheme, while eliminating the second and third completely. Here, integration is performed only over the regions of the unit sphere where fibers are in tension. Gauss–Kronrod quadrature is used across latitudes and the trapezoidal scheme across longitudes. Over a wide range of strain states, fiber material properties, and fiber angular distributions, results demonstrate that this new scheme always outperforms FEI, sometimes by orders of magnitude in the number of computational steps and relative accuracy of the stress calculation.     

Prediction of fibre architecture and adaptation in diseased carotid bifurcations

Many studies have used patient-specific finite element models to estimate the stress environment in atherosclerotic plaques, attempting to correlate the magnitude of stress to plaque vulnerability. In complex geometries, few studies have incorporated the anisotropic material response of arterial tissue. This paper presents a fibre remodelling algorithm to predict the fibre architecture, and thus anisotropic material response in four patient-specific models of the carotid bifurcation. The change in fibre architecture during disease progression and its affect on the stress environment in the plaque were predicted. The mean fibre directions were assumed to lie at an angle between the two positive principal strain directions. The angle and the degree of dispersion were assumed to depend on the ratio of principal strain values. Results were compared with experimental observations and other numerical studies. In non-branching regions of each model, the typical double helix arterial fibre pattern was predicted while at the bifurcation and in regions of plaque burden, more complex fibre architectures were found. The predicted change in fibre architecture in the arterial tissue during plaque progression was found to alter the stress environment in the plaque. This suggests that the specimen-specific anisotropic response of the tissue should be taken into account to accurately predict stresses in the plaque. Since determination of the fibre architecture in vivo is a difficult task, the system presented here provides a useful method of estimating the fibre architecture in complex arterial geometries.     

Determination of stress levels and profiles in the periodontal ligament by means of an improved three-dimensional finite element model for various types of orthodontic and natural force systems

This study was undertaken to investigate the stress-strain levels and distribution within the periodontal ligament for various types of physiological and orthodontic force systems, assuming that the bone resorption process, leading to tooth movements, is partly controlled by those conditions. Two finite element models were developed, simulating a full and partial mandibular morphology, respectively. Both models were based on morphology and physical parameters of human autopsy material. The effect of changing material parameters and structure, type of boundary conditions, calculation method and fineness of the model on the stress levels and profiles in the periodontal ligament was evaluated by a series of tests. A structure optimization technique was used to investigate the load bearing characteristics of the mandible and the influence of the anisotropic material properties of both the mandible and the segment. A ‘multiple modelling’ technique based on both the mandible and the segment was developed to test various types of boundary conditions in the analysis of the segment. Results presented as ‘stress profiles’ showing the correlation between the applied force system and the stress distribution in the periodontal ligament, based on the improved finite element models, were established.     

A microstructurally based continuum model of cartilage viscoelasticity and permeability incorporating measured statistical fiber orientations

The remarkable mechanical properties of cartilage derive from an interplay of isotropically distributed, densely packed and negatively charged proteoglycans; a highly anisotropic and inhomogeneously oriented fiber network of collagens; and an interstitial electrolytic fluid. We propose a new 3D finite strain constitutive model capable of simultaneously addressing both solid (reinforcement) and fluid (permeability) dependence of the tissue’s mechanical response on the patient-specific collagen fiber network. To represent fiber reinforcement, we integrate the strain energies of single collagen fibers—weighted by an orientation distribution function (ODF) defined over a unit sphere—over the distributed fiber orientations in 3D. We define the anisotropic intrinsic permeability of the tissue with a structure tensor based again on the integration of the local ODF over all spatial fiber orientations. By design, our modeling formulation accepts structural data on patient-specific collagen fiber networks as determined via diffusion tensor MRI. We implement our new model in 3D large strain finite elements and study the distributions of interstitial fluid pressure, fluid pressure load support and shear stress within a cartilage sample under indentation. Results show that the fiber network dramatically increases interstitial fluid pressure and focuses it near the surface. Inhomogeneity in the tissue’s composition also increases fluid pressure and reduces shear stress in the solid. Finally, a biphasic neo-Hookean material model, as is available in commercial finite element codes, does not capture important features of the intra-tissue response, e.g., distributions of interstitial fluid pressure and principal shear stress.     

Mechanical characterisation of in vivo human skin using a 3D force-sensitive micro-robot and finite element analysis

The complex mechanical properties of skin have been the subject of much study in recent years. Several experimental methods developed to measure the mechanical properties of skin in vivo, such as suction or torsion, are unable to measure skin’s anisotropic characteristics. An experiment characterising the mechanical properties of in vivo human skin using a novel force-sensitive micro-robot is presented. The micro-robot applied in-plane deformations to the anterior forearm and the posterior upper arm. The behaviour of the skin in each area is highly nonlinear, anisotropic, and viscoelastic. The response of the upper arm skin is very dependent on the orientation of the arm. A finite element model consisting of an Ogden strain energy function and quasi-linear viscoelasticity was developed to simulate the experiments. An orthogonal initial stress field, representing the in vivo skin tension, was used as an additional model parameter. The model simulated the experiments accurately with an error-of-fit of 17.5% for the anterior lower forearm area, 6.5% for the anterior upper forearm and 9.3% for the posterior upper arm. The maximum in vivo tension in each area determined by the model was 6.2 Nm⁻¹ in the anterior lower forearm, 11.4 Nm⁻¹ in anterior upper forearm and 5.6 Nm⁻¹ in the posterior upper arm. The results also show that a finite element model with a neo-Hookean strain energy function cannot simulate the experiments with the same accuracy.     

Dynamic finite element implementation of nonlinear,anisotropic hyperelastic biological membranes

We present a novel method for the implementation of hyperelastic finite strain, non-linear strain-energy functions for biological membranes in an explicit finite element environment. The technique is implemented in LS-DYNA but may also be implemented in any suitable non-linear explicit code. The constitutive equations are implemented on the foundation of a co-rotational uniformly reduced Hughes-Liu shell. This shell is based on an updated-Lagrangian formulation suitable for relating Cauchy stress to the rate-of-deformation, i.e. hypo-elasticity. To accommodate finite deformation hyper-elastic formulations, a co-rotational deformation gradient is assembled over time, resulting in a formulation suitable for pseudo-hyperelastic constitutive equations that are standard assumptions in biomechanics. Our method was validated by comparison with (1) an analytic solution to a spherically-symmetric dynamic membrane inflation problem, incorporating a Mooney-Rivlin hyperelastic equation and (2) with previously published finite element solutions to a non-linear transversely isotropic inflation problem. Finally, we implemented a transversely isotropic strain-energy function for mitral valve tissue. The method is simple and accurate and is believed to be generally useful for anyone who wishes to model biologic membranes with an experimentally driven strain-energy function.     

A nonlinear anisotropic model for porcine aortic heart valves

The anisotropic property of porcine aortic valve leaflet has potentially significant effects on its mechanical behaviour and the failure mechanisms. However, due to its complex nature, testing and modelling the anisotropic porcine aortic valves remains a continuing challenge to date. This study has developed a nonlinear anisotropic finite element model for porcine heart valves. The model is based on the uniaxial experimental data of porcine aortic heart valve leaflet and the properties of nonlinear composite material. A finite element code is developed to solve this problem using the 8-node super-parameter nonlinear shells and the update Lagrangian method. The stress distribution and deformation of the porcine aortic valves with either uniform and non-uniform thicknesses in closed phase and loaded condition are calculated. The results showed significant changes in the stress distributions due to the anisotropic property of the leaflets. Compared with the isotropic valve at the same loading condition, it is found that the site of the peak stress of the anisotropic leaflet is different; the maximum longitudinal normal stress is increased, but the maximum transversal normal stress and in-plane shear stress are reduced. We conclude that it is very important to consider the anisotropic property of the porcine heart valves in order to understand the failure mechanism of such valves in vivo.     

Bone morphology allows estimation of loading history in a murine model of bone adaptation

Bone adapts its morphology (density/micro- architecture) in response to the local loading conditions in such a way that a uniform tissue loading is achieved (‘Wolff’s law’). This paradigm has been used as a basis for bone remodeling simulations to predict the formation and adaptation of trabecular bone. However, in order to predict bone architectural changes in patients, the physiological external loading conditions, to which the bone was adapted, need to be determined. In the present study, we developed a novel bone loading estimation method to predict such external loading conditions by calculating the loading history that produces the most uniform bone tissue loading. We applied this method to murine caudal vertebrae of two groups that were in vivo loaded by either 0 or 8 N, respectively. Plausible load cases were sequentially applied to micro-finite element models of the mice vertebrae, and scaling factors were calculated for each load case to derive the most uniform tissue strain-energy density when all scaled load cases are applied simultaneously. The bone loading estimation method was able to predict the difference in loading history of the two groups and the correct load magnitude for the loaded group. This result suggests that the bone loading history can be estimated from its morphology and that such a method could be useful for predicting the loading history for bone remodeling studies or at sites where measurements are difficult, as in bone in vivo or fossil bones.     

A nonlinear finite element model of cartilage growth

The long range objective of this work is to develop a cartilage growth finite element model (CGFEM), based on the theories of growing mixtures that has the capability to depict the evolution of the anisotropic and inhomogeneous mechanical properties, residual stresses, and nonhomogeneities that are attained by native adult cartilage. The CGFEM developed here simulates isotropic in vitro growth of cartilage with and without mechanical stimulation. To accomplish this analysis a commercial finite element code (ABAQUS) is combined with an external program (MATLAB) to solve an incremental equilibrium boundary value problem representing one increment of growth. This procedure is repeated for as many increments as needed to simulate the desired growth protocol. A case study is presented utilizing a growth law dependent on the magnitude of the diffusive fluid velocity to simulate an in vitro dynamic confined compression loading protocol run for 2 weeks. The results include changes in tissue size and shape, nonhomogeneities that develop in the tissue, as well as the variation that occurs in the tissue constitutive behavior from growth.     

Finite element implementation of a new model of slight compressibility for transversely isotropic materials

Modelling transversely isotropic materials in finite strain problems is a complex task in biomechanics, and is usually addressed by using finite element (FE) simulations. The standard method developed to account for the quasi-incompressible nature of soft tissues is to decompose the strain energy function (SEF) into volumetric and deviatoric parts. However, this decomposition is only valid for fully incompressible materials, and its use for slightly compressible materials yields an unphysical response during the simulation of hydrostatic tension/compression of a transversely isotropic material. This paper presents the FE implementation as subroutines of a new volumetric model solving this deficiency in two FE codes: Abaqus and FEBio. This model also has the specificity of restoring the compatibility with small strain theory. The stress and elasticity tensors are first derived for a general SEF. This is followed by a successful convergence check using a particular SEF and a suite of single-element tests showing that this new model does not only correct the hydrostatic deficiency but may also affect stresses during shear tests (Poynting effect) and lateral stretches during uniaxial tests (Poisson's effect). These FE subroutines have numerous applications including the modelling of tendons, ligaments, heart tissue, etc. The biomechanics community should be aware of specificities of the standard model, and the new model should be used when accurate FE results are desired in the case of compressible materials.     

A calcium-driven mechanochemical model for prediction of force generation in smooth muscle

A new model for the mechanochemical response of smooth muscle is presented. The focus is on the res- ponse of the actin–myosin complex and on the related generation of force (or stress). The chemical (kinetic) model describes the cross-bridge interactions with the thin filament in which the calcium-dependent myosin phosphorylation is the only regulatory mechanism. The new mechanical model is based on Hill’s three-component model and it includes one internal state variable that describes the contraction/relaxation of the contractile units. It is characterized by a strain-energy function and an evolution law incorporating only a few material parameters with clear physical meaning. The proposed model satisfies the second law of thermodynamics. The results of the combined coupled model are broadly consistent with isometric and isotonic experiments on smooth muscle tissue. The simulations suggest that the matrix in which the actin–myosin complex is embedded does have a viscous property. It is straightforward for implementation into a finite element program in order to solve more complex boundary-value problems such as the control of short-term changes in lumen diameter of arteries due to mechanochemical signals.     

A new constitutive model for multi-layered collagenous tissues

Collagenous tissues such as the aneurysmal wall or the aorta are multi-layered structures with the mean fibre alignments distinguishing one layer from another. A constitutive representation of the multiple collagen layers is not yet developed, and hence the aim of the present study. The proposed model is based on the constitutive theory of finite elasticity and is characterized by an anisotropic strain-energy function which takes the material structure into account. The passive tissue behaviour is modelled and the related mechanical response is assumed to be dominated by elastin and collagen. While elastin is modelled by the neo-Hookean material the constitutive response of collagen is assumed to be transversely isotropic for each individual layer and based on an exponential function. The proposed constitutive function is polyconvex which ensures material stability. The model has five independent material parameters, each of which has a clear physical interpretation: the initial stiffnesses of the collagen fabric in the two principal directions, the shear modulus pertaining to the non-collagenous matrix material, a parameter describing the level of nonlinearity of the collagen fabric, and the angle between the principal directions of the collagen fabric and the reference coordinate system. An extension-inflation test of the adventitia of a human femoral artery is simulated by means of the finite element method and an error function is minimized by adjusting the material parameters yielding a good agreement between the model and the experimental data.     

The influence and biomechanical role of cartilage split line pattern on tibiofemoral cartilage stress distribution during the stance phase of gait

This study presents an evaluation of the role that cartilage fibre ‘split line’ orientation plays in informing femoral cartilage stress patterns. A two-stage model is presented consisting of a whole knee joint coupled to a tissue-level cartilage model for computational efficiency. The whole joint model may be easily customised to any MRI or CT geometry using free-form deformation. Three ‘split line’ patterns (medial–lateral, anterior–posterior and random) were implemented in a finite element model with constitutive properties referring to this ‘split line’ orientation as a finite element fibre field. The medial–lateral orientation was similar to anatomy and was derived from imaging studies. Model predictions showed that ‘split lines’ are formed along the line of maximum principal strains and may have a biomechanical role of protecting the cartilage by limiting the cartilage deformation to the area of higher cartilage thickness.     

A bimodular polyconvex anisotropic strain energy function for articular cartilage

A strain energy function for finite deformations is developed that has the capability to describe the nonlinear, anisotropic, and asymmetric mechanical response that is typical of articular cartilage. In particular, the bimodular feature is employed by including strain energy terms that are only mechanically active when the corresponding fiber directions are in tension. Furthermore, the strain energy function is a polyconvex function of the deformation gradient tensor so that it meets material stability criteria. A novel feature of the model is the use of bimodular and polyconvex "strong interaction terms" for the strain invariants of orthotropic materials. Several regression analyses are performed using a hypothetical experimental dataset that captures the anisotropic and asymmetric behavior of articular cartilage. The results suggest that the main advantage of a model employing the strong interaction terms is to provide the capability for modeling anisotropic and asymmetric Poisson's ratios, as well as axial stress-axial strain responses, in tension and compression for finite deformations.     

Stress profile of infant rib in the setting of child abuse: A finite element parametric study

The primary goal of this study is to advance our current understanding of infant rib injuries in the setting of child abuse. To this end, we employed finite element model simulations to determine the sensitivity of an infant rib's stress response to varying material properties and under varying degrees of anterior-posterior chest compression. Using high-resolution chest CT images obtained from a 6-day-old infant, we constructed a simplified geometric model consisting of bone and cartilage structures. To simulate the lateral gripping of an infant in child abuse, an anterior-posterior chest compression load was applied to cause increased stresses along the costovertebral articulation, a classic site for inflicted rib fractures. A sensitivity analysis was conducted to quantify the effects of varying Young's modulus and Poisson's ratio of the bones and cartilages. In addition, we varied the amount of anterior-posterior chest displacement to assess the sensitivity of this parameter to the rib's stress profile. We found that varying Young's modulus of the bone and cartilage not only changed the magnitude but also the shape profile of the rib's stress response. In contrast, varying the degree of chest compression only changed the magnitude of the stress response and not the shape profile. We also discovered that by varying Poisson's ratio of the bone and cartilage, no appreciable change was seen in the magnitude or the shape profile of the rib's stress response. Finite element modeling shows promise as a tool to elucidate the mechanisms of rib fractures in abused infants.     

Changes in the Fabric and Compliance Tensors of Cancellous Bone due to Trabecular Surface Remodeling,Predicted by a Digital Image-based Model

Fabric and compliance tensors of a cube of cancellous bone with a complicated three-dimensional trabecular structure were obtained for trabecular surface remodeling by using a digital image-based model combined with a large-scale finite element method. Using mean intercept length and a homogenization method, the fabric and compliance tensors were determined for the trabecular structure obtained in the computer remodeling simulation. The tensorial quantities obtained indicated that anisotropic structural changes occur in cancellous bone adapting to the compressive loading condition. There were good correlations between the fabric tensor, bone volume fraction, and compliance tensor in the remodeling process. The result demonstrates that changes in the structural and mechanical properties of cancellous bone are essentially anisotropic and should be expressed by tensorial quantities.