Finite element modeling of human skin using an isotropic, nonlinear elastic constitutive model

The collagen network in skin is largely responsible for the nonlinear mechanical stress-strain response of skin. We hypothesize that the force-stretch response of collagen is governed by the entropics of long-chain molecules. We show that a constitutive model derived from the statistical mechanics of long-chain molecules, corresponding to the fibrous collagen network in skin, captures the mechanical response of skin. A connection between the physiologically meaningful parameters of network molecular chain density and free length of collagen fibers and the constitutively significant parameters of initial modulus and limiting stretch is thus established. The relevant constitutive law is shown to have predictive capabilities related to skin histology by replicating in vivo and in vitro experimental results. From finite element simulations, this modeling approach predicts that the collagen network in hypertrophic scars is more dense and the constituent collagen fibers have shorter free lengths than in healthy skin. Additionally, the model is shown to predict that as rat skin ages, collagen network density increases and fiber free length decreases. The importance of knowledge of the in situ stress state for analyzing skin response and validating constitutive laws is also demonstrated.     

Pipette aspiration technique for the measurement of nonlinear and anisotropic mechanical properties of blood vessel walls under biaxial stretch

A pipette aspiration technique was proposed for the measurement of nonlinear mechanical properties of arteries under biaxial stretching. A cross-shaped specimen of porcine thoracic aorta whose principal axes corresponded with the axial and circumferential directions of the aortic walls was excised. The intraluminal surface of the specimen was aspirated with a circular cross-sectioned glass pipette while the specimen was stretching in the axial and circumferential directions in 10% increments. The elastic modulus agreed with the incremental elastic modulus obtained through a conventional pressure-diameter test of the same specimen to within an error of 30% at a circumferential stretch ratio below 1.3 and an axial stretch ratio of 1.0, 1.1 or 1.2, which represent lower range of physiological stretch ratios for the porcine aorta. A rectangular cross-sectioned pipette was utilized to measure anisotropic properties of the specimen under biaxial stretching. When aspirated with such a pipette, the specimens' elastic properties along the length of the rectangular pipette cross section can be neglected. The elastic modulus was found to increase rapidly when the specimen was stretched in the direction of the pipette's width. Thus, pipette aspiration should have many advantages such as well measurement of the local nonlinear and anisotropic mechanical properties of blood vessel walls.     

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.     

A nonlinear viscoelastic finite element model of polyethylene

A nonlinear viscoelastic finite element model of ultra-high molecular weight polyethylene (UHMWPE) was developed in this study. Eight cylindrical specimens were machined from ram extruded UHMWPE bar stock (GUR 1020) and tested under constant compression at 7% strain for 100 sec. The stress strain data during the initial ramp up to 7% strain was utilized to model the instantaneous stress-strain response using a Mooney-Rivlin material model. The viscoelastic behavior was modeled using the time-dependent relaxation in stress seen after the initial maximum stress was achieved using a stored energy formulation. A cylindrical model of similar dimensions was created using a finite element analysis software program. The cylinder was made up of hexahedral elements, which were given the material properties utilizing the instantaneous stress-strain curve and the energy-relaxation curve obtained from the experimental data. The cylinder was compressed between two flat rigid bodies that simulated the fixtures of the testing machine. Experimental stress-relaxation, creep and dynamic testing data were then used to validate the model. The mean error for predicted versus experimental data for stress relaxation at different strain levels was 4.2%. The mean error for the creep test was 7% and for dynamic test was 5.4%. Finally, dynamic loading in a hip arthroplasty was modeled and validated experimentally with an error of 8%. This study establishes a working finite element material model of UHMWPE that can be utilized to simulate a variety of postoperative arthroplasty conditions.     

Transient finite element modeling of functional electrical stimulation

Transcutaneous functional electrical stimulation is commonly used for strengthening muscle. However, transient effects during stimulation are not yet well explored. The effect of an amplitude change of the stimulation can be described by static model, but there is no differency for different pulse duration. The aim of this study is to present the finite element (FE) model of a transient electrical stimulation on the forearm. Discrete FE equations were derived by using a standard Galerkin procedure. Different tissue conductive and dielectric properties are fitted using least square method and trial and error analysis from experimental measurement. This study showed that FE modeling of electrical stimulation can give the spatial-temporal distribution of applied current in the forearm. Three different cases were modeled with the same geometry but with different input of the current pulse, in order to fit the tissue properties by using transient FE analysis. All three cases were compared with experimental measurements of intramuscular voltage on one volunteer.     

Merging of intersecting triangulations for finite element modeling

Surface mesh generation over intersecting triangulations is a problem common to many branches of biomechanics. A new strategy for merging intersecting triangulations is described. The basis of the method is that object surfaces are represented as the zero-level iso-surface of the distance-to-surface function defined on a background grid. Thus, the triangulation of intersecting objects reduces to the extraction of an iso-surface from an unstructured grid. In a first step, a regular background mesh is constructed. For each point of the background grid, the closest distance to the surface of each object is computed. Background points are then classified as external or internal by checking the direction of the surface normal at the closest location and assigned a positive or negative distance, respectively. Finally, the zero-level iso-surface is constructed. This is the final triangulation of the intersecting objects. The overall accuracy is enhanced by adaptive refinement of the background grid elements. The resulting surface models are used as support surfaces to generate three-dimensional grids for finite element analysis. The algorithms are demonstrated by merging arterial branches independently reconstructed from contrast-enhanced magnetic resonance images and by adding extra features such as vascular stents. Although the methodology is presented in the context of finite element analysis of blood flow, the algorithms are general and can be applied in other areas as well.     

Porohyperelastic finite element modeling of abdominal aortic aneurysms

Abdominal aortic aneurysm (AAA) is the gradual weakening and dilation of the infrarenal aorta. This disease is progressive, asymptomatic, and can eventually lead to rupture--a catastrophic event leading to massive internal bleeding and possibly death. The mechanical environment present in AAA is currently thought to be important in disease initiation, progression, and diagnosis. In this study, we utilize porohyperelastic (PHE) finite element models (FEMs) to investigate how such modeling can be used to better understand the local biomechanical environment in AAA. A 3D hypothetical AAA was constructed with a preferential anterior bulge assuming both the intraluminal thrombus (ILT) and the AAA wall act as porous materials. A parametric study was performed to investigate how physiologically meaningful variations in AAA wall and ILT hydraulic permeabilities affect luminal interstitial fluid velocities and wall stresses within an AAA. A corresponding hyperelastic (HE) simulation was also run in order to be able to compare stress values between PHE and HE simulations. The effect of AAA size on local interstitial fluid velocity was also investigated by simulating maximum diameters (5.5 cm, 4.5 cm, and 3.5 cm) at the baseline values of ILT and AAA wall permeability. Finally, a cyclic PHE simulation was utilized to study the variation in local fluid velocities as a result of a physiologic pulsatile blood pressure. While the ILT hydraulic permeability was found to have minimal affect on interstitial velocities, our simulations demonstrated a 28% increase and a 20% decrease in luminal interstitial fluid velocity as a result of a 1 standard deviation increase and decrease in AAA wall hydraulic permeability, respectively. Peak interstitial velocities in all simulations occurred on the luminal surface adjacent to the region of maximum diameter. These values increased with increasing AAA size. PHE simulations resulted in 19.4%, 40.1%, and 81.0% increases in peak maximum principal wall stresses in comparison to HE simulations for maximum diameters of 35 mm, 45 mm, and 55 mm, respectively. The pulsatile AAA PHE FEM demonstrated a complex interstitial fluid velocity field the direction of which alternated in to and out of the luminal layer of the ILT. The biomechanical environment within both the aneurysmal wall and the ILT is involved in AAA pathogenesis and rupture. Assuming these tissues to be porohyperelastic materials may provide additional insight into the complex solid and fluid forces acting on the cells responsible for aneurysmal remodeling and weakening.     

Dynamic finite element modeling of poroviscoelastic soft tissue

Clinical evidences relative to biomechanical factors have demonstrated their important contribution to the behaviour of soft tissues. Finite element (FE) analysis is used to study the mechanical behaviour of soft tissue because it can provide numerical solutions to problems that are intractable to analytic solutions. This study focuses on the development of a FE model of a poroelastic biological tissue, which incorporates the viscoelastic material behaviour, finite deformation and inertial effect. The FE formulation is based on the weak form derived from the governing equation, and Newmark-beta method as well as Newton's method is incorporated into the implicit non-linear solutions. One-dimensional analytical solutions were used to verify the theoretical formulation and the numerical implementation of the proposed model. This study was further extended to analyze two-dimensional biomechanical models and the results clearly demonstrate the importance of including finite deformation, viscoelasticity and inertial effects.     

Computer-integrated finite element modeling of human middle ear

 The objective of this study was to produce an improved finite element (FE) model of the human middle ear and to compare the model with human data. We began with a systematic and accurate geometric modeling technique for reconstructing the middle ear from serial sections of a freshly frozen temporal bone. A geometric model of a human middle ear was constructed in a computer-aided design (CAD) environment with particular attention to geometry and microanatomy. Using the geometric model, a working FE model of the human middle ear was created using previously published material properties of middle ear components. This working FE model was finalized by a cross-calibration technique, comparing its predicted stapes footplate displacements with laser Doppler interferometry measurements from fresh temporal bones. The final FE model was shown to be reasonable in predicting the ossicular mechanics of the human middle ear. Received: 18 February 2002 / Accepted: 6 June 2002 The preparation of temporal bone histological sections of Robert K. Dyer, Jr., MD is gratefully recognized. The Whitaker Foundation supported this work (Research Grant RG-98-0305).     

Evolving biaxial mechanical properties of mouse carotid arteries in hypertension

Quantifying the time course of load-induced changes in arterial wall geometry, microstructure, and properties is fundamental to developing mathematical models of growth and remodeling. Arteries adapt to altered pressure and flow by modifying wall thickness, inner diameter, and axial length via marked cell and matrix turnover. To estimate particular biomaterial implications of such adaptations, we used a 4-fiber family constitutive relation to quantify passive biaxial mechanical behaviors of mouse carotid arteries 0 (control), 7-10, 10-14, or 35-56 days after an aortic arch banding surgery that increased pulse pressure and pulsatile flow in the right carotid artery. In vivo circumferential and axial stretches at mean arterial pressure were, for example, 11% and 26% lower, respectively, in hypertensive carotids 35-56 days after banding than in normotensive controls; this finding is consistent with observations that hypertension decreases distensibility. Interestingly, the strain energy W stored in the carotids at individual in vivo conditions was also less in hypertensive compared with normotensive carotids. For example, at 35-56 days after banding, W was 24%, 39%, and 47% of normal values at diastolic, mean, and systolic pressures, respectively. The energy stored during the cardiac cycle, W(sys)-W(dias), also tended to be less, but this reduction did not reach significance. When computed at normal in vivo values of biaxial stretch, however, W was well above normal for the hypertensive carotids. This net increase resulted from an overall increase in the collagen-related anisotropic contribution to W despite a decrease in the elastin-related isotropic contribution. The latter was consistent with observed decreases in the mass fraction of elastin.     

A novel two-layer, coupled finite element approach for modeling the nonlinear elastic and viscoelastic behavior of human erythrocytes

A novel finite element approach is presented to simulate the mechanical behavior of human red blood cells (RBC, erythrocytes). As the RBC membrane comprises a phospholipid bilayer with an intervening protein network, we propose to model the membrane with two distinct layers. The fairly complex characteristics of the very thin lipid bilayer are represented by special incompressible solid shell elements and an anisotropic viscoelastic constitutive model. Properties of the protein network are modeled with an isotropic hyperelastic third-order material. The elastic behavior of the model is validated with existing optical tweezers studies with quasi-static deformations. Employing material parameters consistent with literature, simulation results are in excellent agreement with experimental data. Available models in literature neglect either the surface area conservation of the RBC membrane or realistic loading conditions of the optical tweezers experiments. The importance of these modeling assumptions, that are both included in this study, are discussed and their influence quantified. For the simulation of the dynamic motion of RBC, the model is extended to incorporate the cytoplasm. This is realized with a monolithic fully coupled fluid-structure interaction simulation, where the fluid is described by the incompressible Navier–Stokes equations in an arbitrary Lagrangian Eulerian framework. It is shown that both membrane viscosity and cytoplasm viscosity have significant influence on simulation results. Characteristic recovery times and energy dissipation for varying strain rates in dynamic laser trap experiments are calculated for the first time and are found to be comparable with experimental data.     

3D finite element modeling of pelvic organ prolapse

Objectives: The purpose of this study is to develop a validated 3D finite element model of the pelvic floor system which can offer insights into the mechanics of anterior vaginal wall prolapse and have the ability to assess biomedical device treatment methods. The finite element results should accurately mimic the clinical findings of prolapse due to intra-abdominal pressure (IAP) and soft tissues impairment conditions. Methods: A 3D model of pelvic system was created in Creo Parametric 2.0 based on MRI Images, which included uterus, cervix, vagina, cardinal ligaments, uterosacral ligaments, and a simplified levator plate and rectum. The geometrical model was imported into ANSYS Workbench 14.5. Mechanical properties of soft tissues were based on experimental data of tensile test results from current literature. Studies were conducted for IAP loadings on the vaginal wall and uterus, increasing from lowest to extreme values. Results: Anterior vaginal wall collapse occurred at an IAP value corresponding to maximal valsalva and showed similar collapsed shape as clinical findings. Prolapse conditions exhibited high sensitivity to vaginal wall stiffness, whereas healthy tissues was found to support the vagina against prolapse. Ligament impairment was found to have only a secondary effect on prolapse.     

Explicit finite element modeling of total knee replacement mechanics

Joint kinematics and contact mechanics dictate the success of current total knee replacement (TKR) devices. Efficient computer models present an effective way of evaluating these characteristics. Predicted contact stress and area due to articulations at the tibio-femoral and patello-femoral interfaces indicate potential clinical performance. Previous finite element (FE) knee models have generally been used to predict contact stresses and/or areas during static or quasi-static loading conditions. Explicit dynamic FE analyses have recently been used to efficiently predict TKR kinematics and contact mechanics during dynamic loading conditions. The objective of this study was to develop and experimentally validate an explicit FE TKR model that incorporates tibio-femoral and patello-femoral articulations. For computational efficiency, we developed rigid body analyses that can reasonably reproduce the kinematics, contact pressure distribution, and contact area of a fully deformable system. Results from the deformable model showed that the patello-femoral and tibio-femoral kinematics were in good agreement with experimental knee simulator measurements. Kinematic results from the rigid body analyses were nearly identical to those from the fully deformable model, and the contact pressure and contact area correlation was acceptable given the great reduction in analysis time. Component mesh density studied had little effect on the predicted kinematics, particularly for the patellar component, and small effects on the predicted contact pressure and area. These analyses have shown that, at low computational cost, a force-control dynamic simulation of a gait cycle can yield useful and predictable results.     

Extramuscular myofascial force transmission: experiments and finite element modeling

The specific purpose of the present study was to show that extramuscular myofascial force transmission exclusively has substantial effects on muscular mechanics. Muscle forces exerted at proximal and distal tendons of the rat extensor digitorium longus (EDL) were measured simultaneously, in two conditions (1) with intact extramuscular connections (2) after dissecting the muscles' extramuscular connections to a maximum extent without endangering circulation and innervation (as in most in situ muscle experiments). A finite element model of EDL including the muscles' extramuscular connections was used to assess the effects of extramuscular myofascial force transmission on muscular mechanics, primarily to test if such effects lead to distribution of length of sarcomeres within muscle fibers. In condition (1), EDL isometric forces measured at the distal and proximal tendons were significantly different (F(dist) > F(prox), DeltaF approximates maximally 40% of the proximal force). The model results show that extramuscular myofascial force transmission causes distributions of strain in the fiber direction (shortening in the proximal, lengthening in the distal ends of fibers) at higher lengths. This indicates significant length distributions of sarcomeres arranged in series within muscle fibers. Stress distributions found are in agreement with the higher distal force measured, meaning that the muscle fiber is no longer the unit exerting equal forces at both ends. Experimental results obtained in condition (2) showed no significant changes in the length-force characteristics (i.e., proximo-distal force differences were maintained). This shows that a muscle in situ has to be distinguished from a muscle that is truly isolated in which case the force difference has to be zero. We conclude that extramuscular myofascial force transmission has major effects on muscle functioning.     

Finite elasticity modeling of the biaxial and uniaxial properties of compliant vascular grafts

Compliant vascular grafts were modeled by finite elasticity theory. A linear, biaxial model satisfactorily described the stress-strain behavior in inflation tests, where the sample length was fixed longitudinally and inflated. The model was then used to predict the behavior in a longitudinal test (where the longitudinal stretch was varied while keeping the pressure zero), and in a uniaxial test of a circumferential strip. The model satisfactorily predicted the longitudinal Young's modulus measured in the longitudinal test. The model was less successful in predicting the circumferential Young's modulus measured in the uniaxial test, possibly because the state of stress was not purely uniaxial.     

In vivo characterization of mechanical tissue properties of internal organs using endoscopic microscopy and inverse finite element analysis

Knowledge about mechanical tissue properties is required for functional modelling and simulating of tissue and organ responses to external mechanical stress. To get the right properties especially for functional modelling of organs, tissue properties have to be determined in vivo. There are only few described methods for characterization of internal organ's tissue mechanics that can be applied in vivo. We introduce and evaluate a method to determine mechanical tissue properties, especially those of lung tissue, endoscopically. Inverse finite element analysis (utilizing a Neo-Hookean model for hyperelastic materials) and image processing algorithms are used to determine the shear modulus of a soft tissue. The resulting values for shear moduli were normally distributed. The shear modulus of the artificial tissue sample was determined with a relative error of 0.47% compared to the value obtained by uniaxial tensile test.     

Sensitivity of the tibio-femoral response to finite element modeling parameters

A generic finite element (FE) model of the lower limb was used to study the knee response in-vivo during a one-legged hop. The approach uses an explicit FE code and a combination of estimated muscle forces and measured three-dimensional tibio-femoral kinematics and ground reaction force as input to the FE model. The sensitivity of the simulated tibio-femoral response to variations of key geometric and material parameters was investigated by performing a total of 38 different simulations. The amplitudes of both kinematic and kinetic responses were affected by the change of these parameters. For the current approach, the results suggest that while cartilage mechanical and geometric properties are very important for the estimation of tibio-femoral cartilage pressure, they have limited effects on the overall kinematic response. The study may help to better define the relative importance of modeling parameters for the development of subject-specific models.     

A finite element model of burn injury in blood-perfused skin

The burn process resulting from the application of a hot, cylindrical source to the skin surface was modeled using the finite element technique. A rotationally symmetric 125-element mesh was defined within the tissue beneath and outlying to an applied heating disk. The disk temperature and duration of contact were varied, respectively, between 50 and 100 degrees C for up to 30 s. Natural convection with ambient air was assumed for areas of skin surface not in direct contact with the disk. The simulated thermal history was used in a damage integral model to calculate the extent and severity of injury in the radial and axial dimensions.     

Determination of the axial and circumferential mechanical properties of the skin tissue using experimental testing and constitutive modeling

The skin, being a multi-layered material, is responsible for protecting the human body from the mechanical, bacterial, and viral insults. The skin tissue may display different mechanical properties according to the anatomical locations of a body. However, these mechanical properties in different anatomical regions and at different loading directions (axial and circumferential) of the mice body to date have not been determined. In this study, the axial and circumferential loads were imposed on the mice skin samples. The elastic modulus and maximum stress of the skin tissues were measured before the failure occurred. The nonlinear mechanical behavior of the skin tissues was also computationally investigated through a suitable constitutive equation. Hyperelastic material model was calibrated using the experimental data. Regardless of the anatomic locations of the mice body, the results revealed significantly different mechanical properties in the axial and circumferential directions and, consequently, the mice skin tissue behaves like a pure anisotropic material. The highest elastic modulus was observed in the back skin under the circumferential direction (6.67 MPa), while the lowest one was seen in the abdomen skin under circumferential loading (0.80 MPa). The Ogden material model was narrowly captured the nonlinear mechanical response of the skin at different loading directions. The results help to understand the isotropic/anisotropic mechanical behavior of the skin tissue at different anatomical locations. They also have implications for a diversity of disciplines, i.e., dermatology, cosmetics industry, clinical decision making, and clinical intervention.