A meta-model analysis of a finite element simulation for defining poroelastic properties of intervertebral discs

Finite element analysis is an effective tool to evaluate the material properties of living tissue. For an interactive optimization procedure, the finite element analysis usually needs many simulations to reach a reasonable solution. The meta-model analysis of finite element simulation can be used to reduce the computation of a structure with complex geometry or a material with composite constitutive equations. The intervertebral disc is a complex, heterogeneous, and hydrated porous structure. A poroelastic finite element model can be used to observe the fluid transferring, pressure deviation, and other properties within the disc. Defining reasonable poroelastic material properties of the anulus fibrosus and nucleus pulposus is critical for the quality of the simulation. We developed a material property updating protocol, which is basically a fitting algorithm consisted of finite element simulations and a quadratic response surface regression. This protocol was used to find the material properties, such as the hydraulic permeability, elastic modulus, and Poisson’s ratio, of intact and degenerated porcine discs. The results showed that the in vitro disc experimental deformations were well fitted with limited finite element simulations and a quadratic response surface regression. The comparison of material properties of intact and degenerated discs showed that the hydraulic permeability significantly decreased but Poisson’s ratio significantly increased for the degenerated discs. This study shows that the developed protocol is efficient and effective in defining material properties of a complex structure such as the intervertebral disc.

[1]  J. Laible,et al.  Refinement of Elastic, Poroelastic, and Osmotic Tissue Properties of Intervertebral Disks to Analyze Behavior in Compression , 2010, Annals of Biomedical Engineering.

[2]  G Josse,et al.  Dual-parameter optimisation of the elastic properties of skin , 2012, Computer methods in biomechanics and biomedical engineering.

[3]  Jaw-Lin Wang,et al.  Rheology of Intervertebral Disc: An Ex Vivo Study on the Effect of Loading History, Loading Magnitude, Fatigue Loading, and Disc Degeneration , 2010, Spine.

[4]  J. Lotz,et al.  Effect of Frozen Storage on the Creep Behavior of Human Intervertebral Discs , 2001, Spine.

[5]  A Ratcliffe,et al.  Compressive mechanical properties of the human anulus fibrosus and their relationship to biochemical composition. , 1994, Spine.

[6]  Lutz Claes,et al.  Creep associated changes in intervertebral disc bulging obtained with a laser scanning device. , 2007, Clinical biomechanics.

[7]  JD Jan Janssen,et al.  A validation of the quadriphasic mixture theory for intervertebral disc tissue , 1997 .

[8]  S. Delp,et al.  New resource for the computation of cartilage biphasic material properties with the interpolant response surface method , 2009, Computer methods in biomechanics and biomedical engineering.

[9]  Kaifeng Liu,et al.  Mechanical characterization of soft viscoelastic gels via indentation and optimization-based inverse finite element analysis. , 2009, Journal of the mechanical behavior of biomedical materials.

[10]  V C Mow,et al.  The anisotropic hydraulic permeability of human lumbar anulus fibrosus. Influence of age, degeneration, direction, and water content. , 1999, Spine.

[11]  Zhen Zhao,et al.  Optimization of drawbead design in sheet forming using one step finite element method coupled with response surface methodology , 2006 .

[12]  D. Hillel.  Introduction to environmental soil physics , 1982 .

[13]  Lutz Claes,et al.  Application of a new calibration method for a three-dimensional finite element model of a human lumbar annulus fibrosus. , 2006, Clinical biomechanics.

[14]  J. D. Janssen,et al.  Nonhomogeneous Permeability of Canine Anulus Fibrosus , 1997, Spine.

[15]  M. Biot.  General Theory of Three‐Dimensional Consolidation , 1941 .

[16]  V C Mow,et al.  Degeneration affects the anisotropic and nonlinear behaviors of human anulus fibrosus in compression. , 1998, Journal of biomechanics.

[17]  Fulin Lei,et al.  Inverse analysis of constitutive models: biological soft tissues. , 2007, Journal of biomechanics.

[18]  Mohammad Nikkhoo,et al.  DYNAMIC RESPONSES OF INTERVERTEBRAL DISC DURING STATIC CREEP AND DYNAMIC CYCLIC LOADING: A PARAMETRIC POROELASTIC FINITE ELEMENT ANALYSIS , 2013 .

[19]  D. Elliott,et al.  Effects of Degeneration on the Biphasic Material Properties of Human Nucleus Pulposus in Confined Compression , 2005, Spine.

[20]  M Parnianpour,et al.  Predictive equations to estimate spinal loads in symmetric lifting tasks. , 2011, Journal of biomechanics.

[21]  K. Miller,et al.  Total Lagrangian explicit dynamics finite element algorithm for computing soft tissue deformation , 2006 .

[22]  Gunnar B J Andersson,et al.  Inclusion of regional poroelastic material properties better predicts biomechanical behavior of lumbar discs subjected to dynamic loading. , 2007, Journal of biomechanics.

[23]  D. Lacroix,et al.  Comparison of four methods to simulate swelling in poroelastic finite element models of intervertebral discs. , 2011, Journal of the mechanical behavior of biomedical materials.

[24]  J M Huyghe,et al.  Experimental and model determination of human intervertebral disc osmoviscoelasticity , 2008, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[25]  M. Parnianpour,et al.  Search for critical loading condition of the spine–A meta analysis of a nonlinear viscoelastic finite element model , 2005, Computer methods in biomechanics and biomedical engineering.

[26]  Delphine Périé,et al.  Confined compression experiments on bovine nucleus pulposus and annulus fibrosus: sensitivity of the experiment in the determination of compressive modulus and hydraulic permeability. , 2005, Journal of biomechanics.

[27]  Wei Shyy,et al.  Response surface techniques for diffuser shape optimization , 1997 .

[28]  B R Simon,et al.  1985 Volvo Award in Biomechanics: Poroelastic Dynamic Structural Models of Rhesus Spinal Motion Segments , 1985, Spine.

[29]  L. Setton,et al.  Anisotropic and inhomogeneous tensile behavior of the human anulus fibrosus: experimental measurement and material model predictions. , 2001, Journal of biomechanical engineering.

[30]  N Dhillon,et al.  The internal mechanics of the intervertebral disc under cyclic loading. , 2002, Journal of biomechanics.

[31]  J M Huyghe,et al.  Confined compression of canine annulus fibrosus under chemical and mechanical loading. , 1995, Journal of biomechanical engineering.

[32]  S. Laporte,et al.  The non-linear response of a muscle in transverse compression: assessment of geometry influence using a finite element model , 2012, Computer methods in biomechanics and biomedical engineering.

[33]  P. Riches,et al.  Determination of the strain-dependent hydraulic permeability of the compressed bovine nucleus pulposus. , 2008, Journal of Biomechanics.

[34]  G A Ateshian,et al.  Experimental verification and theoretical prediction of cartilage interstitial fluid pressurization at an impermeable contact interface in confined compression. , 1998, Journal of biomechanics.

[35]  Suvranu De,et al.  An efficient soft tissue characterization algorithm from in vivo indentation experiments for medical simulation , 2008, The international journal of medical robotics + computer assisted surgery : MRCAS.

[36]  A Shirazi-Adl,et al.  Poroelastic creep response analysis of a lumbar motion segment in compression. , 1996, Journal of biomechanics.

[37]  H. Yao,et al.  Effects of Swelling Pressure and Hydraulic Permeability on Dynamic Compressive Behavior of Lumbar Annulus Fibrosus , 2002, Annals of Biomedical Engineering.

[38]  Aboulfazl Shirazi-Adl,et al.  Response analysis of the lumbar spine during regular daily activities--a finite element analysis. , 2010, Journal of biomechanics.

[39]  N. Filipovic,et al.  A finite element solution procedure for porous medium with fluid flow and electromechanical coupling , 1998 .