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.
BestBAGuilakFSettonLA. Compressive mechanical properties of the human anulus fibrosus and their relationship to biochemical composition. Spine (Phila Pa 1976)1994; 19(2): 212–221.
2.
GuWYMaoXGFosterRJ. The anisotropic hydraulic permeability of human lumbar anulus fibrosus: influence of age, degeneration, direction, and water content. Spine (Phila Pa 1976)1999; 24(23): 2449–2455.
3.
IatridisJCSettonLAFosterRJ. Degeneration affects the anisotropic and nonlinear behaviors of human anulus fibrosus in compression. J Biomech1998; 31(6): 535–544.
4.
JohannessenWElliottDM. Effects of degeneration on the biphasic material properties of human nucleus pulposus in confined compression. Spine (Phila Pa 1976)2005; 30(24): E724–E729.
5.
PerieDKordaDIatridisJC. Confined compression experiments on bovine nucleus pulposus and annulus fibrosus: sensitivity of the experiment in the determination of compressive modulus and hydraulic permeability. J Biomech2005; 38(11): 2164–2171.
6.
DhillonNBassECLotzJC. Effect of frozen storage on the creep behavior of human intervertebral discs. Spine (Phila Pa 1976)2001; 26(8): 883–888.
7.
KuoYWWangJL. Rheology of intervertebral disc: an ex vivo study on the effect of loading history, loading magnitude, fatigue loading, and disc degeneration. Spine (Phila Pa 1976)2010; 35(16): E743–E752.
8.
StokesIALaibleJPGardner-MorseMG. Refinement of elastic, poroelastic, and osmotic tissue properties of intervertebral disks to analyze behavior in compression. Ann Biomed Eng2011; 39(1): 122–131.
9.
HeuerFSchmittHSchmidtH. Creep associated changes in intervertebral disc bulging obtained with a laser scanning device. Clin Biomech (Bristol, Avon)2007; 22(7): 737–744.
10.
MillerKJoldesGLanceD. Total Lagrangian explicit dynamics finite element algorithm for computing soft tissue deformation. Commun Numer Meth En2006; 23(2): 121–134.
11.
KimJAhnBDeS. An efficient soft tissue characterization algorithm from in vivo indentation experiments for medical simulation. Int J Med Robot2008; 4(3): 277–285.
SchroederYElliottDMWilsonW. Experimental and model determination of human intervertebral disc osmoviscoelasticity. J Orthop Res2008; 26(8): 1141–1146.
14.
GrasLLMittonDCrevier-DenoixN. The non-linear response of a muscle in transverse compression: assessment of geometry influence using a finite element model. Comput Methods Biomech Biomed Engin2012; 15(1): 13–21.
15.
DelalleauAJosseGLagardeJM. Dual-parameter optimisation of the elastic properties of skin. Comput Methods Biomech Biomed Engin2012; 15(1): 83–92.
16.
LiuKVanLandinghamMROvaertTC. Mechanical characterization of soft viscoelastic gels via indentation and optimization-based inverse finite element analysis. J Mech Behav Biomed Mater2009; 2(4), 355–362; discussion 362–363.
WangJLShirazi-AdlAParnianpourM. Search for critical loading condition of the spine—a meta analysis of a nonlinear viscoelastic finite element model. Comput Methods Biomech Biomed Engin2005; 8(5): 323–330.
19.
ArjmandNPlamondonAShirazi-AdlA. Predictive equations to estimate spinal loads in symmetric lifting tasks. J Biomech2011; 44(1): 84–91.
20.
TangBSunJZhaoZ. Optimization of drawbead design in sheet forming using one step finite element method coupled with response surface methodology. Int J Adv Manuf Tech2006; 31(3–4): 225–234.
21.
KeenanKEKourtisLCBesierTF. New resource for the computation of cartilage biphasic material properties with the interpolant response surface method. Comput Methods Biomech Biomed Engin2009; 12(4): 415–422.
22.
HsuYCChangYCLinMH. Dynamic loading is more effective than forced circulation in transporting small weight molecule in an in-vitro whole disc culturing system. In: ORS 2011 annual meeting, Orthopaedic Research Society, Long Beach, CA, 13-16 January, 2011, Poster No. 609.
23.
SimonBRWuJSCarltonMW. Poroelastic dynamic structural models of rhesus spinal motion segments. Spine (Phila Pa 1976)1985; 10(6): 494–507.
24.
BiotMA. General theory of three-dimensional consolidation. J Appl Phys1941; 12(2): 155–164.
25.
KojicMFilipovicNVulovicS. A finite element solution procedure for porous medium with fluid flow and electromechanical coupling. Commun Numer Meth En1998; 14(4): 381–392.
26.
NikkhooMHaghpanahiMParnianpourM. Dynamic responses of intervertebral disc during static creep and dynamic cyclic loading: a parametric poroelastic finite element analysis. Biomed Eng App Bas C, DOI: 10.1142/S1016237213500130.
27.
LiS. Response of human intervertebral discs to prolonged axial loading and low-frequency vibration. Chicago, IL: University of Illinois at Chicago, 1994.
28.
SchmidtHHeuerFSimonU. Application of a new calibration method for a three-dimensional finite element model of a human lumbar annulus fibrosus. Clin Biomech2006; 21(4): 337–344.
WilliamsJRNatarajanRNAnderssonGB. Inclusion of regional poroelastic material properties better predicts biomechanical behavior of lumbar discs subjected to dynamic loading. J Biomech2007; 40(9): 1981–1987.
31.
GalbuseraFSchmidtHNoaillyJ. Comparison of four methods to simulate swelling in poroelastic finite element models of intervertebral discs. J Mech Behav Biomed Mater2011; 4(7): 1234–1241.
32.
ArgoubiMShirazi-AdlA. Poroelastic creep response analysis of a lumbar motion segment in compression. J Biomech1996; 29(10): 1331–1339.
33.
SchmidtHShirazi-AdlAGalbuseraF. Response analysis of the lumbar spine during regular daily activities—a finite element analysis. J Biomech2010; 43(10): 1849–1856.
34.
DrostMRWillemsPSnijdersH. Confined compression of canine annulus fibrosus under chemical and mechanical loading. J Biomech Eng1995; 117(4): 390–396.
35.
FrijnsAJHHuygheJMJanssenJD. A validation of the quadriphasic mixture theory for intervertebral disc tissue. Int J Eng Sci1997; 35(15): 1419–1429.
36.
HoubenGBDrostMRHuygheJM. Nonhomogeneous permeability of canine anulus fibrosus. Spine (Phila Pa 1976)1997; 22(1): 7–16.
37.
SoltzMAAteshianGA. Experimental verification and theoretical prediction of cartilage interstitial fluid pressurization at an impermeable contact interface in confined compression. J Biomech1998; 31(10): 927–934.
38.
ElliottDMSettonLA. Anisotropic and inhomogeneous tensile behavior of the human anulus fibrosus: experimental measurement and material model predictions. J Biomech Eng2001; 123(3): 256–263.
39.
RichesPEDhillonNLotzJ. The internal mechanics of the intervertebral disc under cyclic loading. J Biomech2002; 35(9): 1263–1271.
40.
YaoHJustizMAFlaglerD. Effects of swelling pressure and hydraulic permeability on dynamic compressive behavior of lumbar annulus fibrosus. Ann Biomed Eng2002; 30(10): 1234–1241.
41.
HeneghanPRichesPE. Determination of the strain-dependent hydraulic permeability of the compressed bovine nucleus pulposus. J Biomech2008; 41(4): 903–906.
