Physical signals and solute transport in human intervertebral disc during compressive stress relaxation: 3D finite element analysis.

A 3D finite element model for charged hydrated soft tissues containing charged/uncharged solutes was developed based on the multi-phasic mechano-electrochemical mixture theory (Lai et al., J. Biomech. Eng. 113 (1991), 245-258; Gu et al., J. Biomech. Eng. 120 (1998), 169-180). This model was applied to analyze the mechanical, chemical and electrical signals within the human intervertebral disc during an unconfined compressive stress relaxation test. The effects of tissue composition [e.g., water content and fixed charge density (FCD)] on the physical signals and the transport rate of fluid, ions and nutrients were investigated. The numerical simulation showed that, during disc compression, the fluid pressurization was more pronounced at the center (nucleus) region of the disc while the effective (von Mises) stress was higher at the outer (annulus) region. Parametric analyses revealed that the decrease in initial tissue water content (0.7-0.8) increased the peak stress and relaxation time due to the reduction of permeability, causing greater fluid pressurization effect. The electrical signals within the disc were more sensitive to FCD than tissue porosity, and mechanical loading affected the large solute (e.g., growth factor) transport significantly, but not for small solute (e.g., glucose). Moreover, this study confirmed that the interstitial fluid pressurization plays an important role in the load support mechanism of IVD by sharing more than 40% of the total load during disc compression. This study is important for understanding disc biomechanics, disc nutrition and disc mechanobiology.

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

[2]  D Pflaster,et al.  A poroelastic finite element formulation including transport and swelling in soft tissue structures. , 1996, Journal of biomechanical engineering.

[3]  S. Eisenstein,et al.  Biochemical changes in intervertebral disc degeneration. , 1981, Biochimica et biophysica acta.

[4]  Martin H Krag,et al.  Influence of fixed charge density magnitude and distribution on the intervertebral disc: applications of a poroelastic and chemical electric (PEACE) model. , 2003, Journal of biomechanical engineering.

[5]  Hai Yao,et al.  Physical Signals and Solute Transport in Cartilage Under Dynamic Unconfined Compression: Finite Element Analysis , 2004, Annals of Biomedical Engineering.

[6]  W. Deen.  Hindered transport of large molecules in liquid‐filled pores , 1987 .

[7]  V. Mow,et al.  A MIXED FINITE ELEMENT FORMULATION OF TRIPHASIC MECHANO-ELECTROCHEMICAL THEORY FOR CHARGED, HYDRATED BIOLOGICAL SOFT TISSUES , 1999 .

[8]  A R Hargens,et al.  Intervertebral disc nutrition. Diffusion versus convection. , 1986, Clinical orthopaedics and related research.

[9]  Jd Jan Janssen,et al.  Quadriphasic mechanics of swelling incompressible porous media , 1997 .

[10]  A B Schultz,et al.  Nonlinear behavior of the human intervertebral disc under axial load. , 1976, Journal of Biomechanics.

[11]  W M Lai,et al.  A mixture theory for charged-hydrated soft tissues containing multi-electrolytes: passive transport and swelling behaviors. , 1998, Journal of biomechanical engineering.

[12]  A. Grodzinsky,et al.  Contributions of fluid convection and electrical migration to transport in cartilage: relevance to loading. , 1996, Archives of biochemistry and biophysics.

[13]  W M Lai,et al.  Transport of fluid and ions through a porous-permeable charged-hydrated tissue, and streaming potential data on normal bovine articular cartilage. , 1993, Journal of biomechanics.

[14]  J P Laible,et al.  A Poroelastic-Swelling Finite Element Model With Application to the Intervertebral Disc , 1993, Spine.

[15]  A. Maroudas,et al.  Diffusion of small solutes into the intervertebral disc: as in vivo study. , 1978, Biorheology.

[16]  G Ray,et al.  Systems identification for material properties of the intervertebral joint. , 1978, Journal of biomechanics.

[17]  J. Meister,et al.  Static compression of articular cartilage can reduce solute diffusivity and partitioning: implications for the chondrocyte biological response. , 2001, Journal of biomechanics.

[18]  Farshid Guilak,et al.  Site-Specific Molecular Diffusion in Articular Cartilage Measured using Fluorescence Recovery after Photobleaching , 2003, Annals of Biomedical Engineering.

[19]  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.

[20]  Keita Ito,et al.  Fluid flow and convective transport of solutes within the intervertebral disc. , 2004, Journal of biomechanics.

[21]  H. Yao,et al.  Effects of Hydration and Fixed Charge Density on Fluid Transport in Charged Hydrated Soft Tissues , 2003, Annals of Biomedical Engineering.

[22]  A B Schultz,et al.  Finite element stress analysis of an intervertebral disc. , 1974, Journal of biomechanics.

[23]  A B Schultz,et al.  Mechanical response of a simple finite element model of the intervertebral disc under complex loading. , 1984, Journal of biomechanics.

[24]  H. Yao,et al.  Diffusivity of Ions in Agarose Gels and Intervertebral Disc: Effect of Porosity , 2004, Annals of Biomedical Engineering.

[25]  C. Whyne,et al.  Parametric finite element analysis of vertebral bodies affected by tumors. , 2001, Journal of biomechanics.

[26]  R L Spilker,et al.  Mechanical behavior of a simple model of an intervertebral disk under compressive loading. , 1980, Journal of biomechanics.

[27]  K. Kaneda,et al.  A cineradiographic study on the lumbar disc deformation during flexion and extension of the trunk. , 1995, Clinical biomechanics.

[28]  P. Meares,et al.  The diffusion of electrolytes in a cation-exchange resin membrane I. Theoretical , 1955, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[29]  A. M. Ahmed,et al.  Stress analysis of the lumbar disc-body unit in compression. A three-dimensional nonlinear finite element study. , 1984, Spine.

[30]  V. Mow,et al.  Fixed negative charges modulate mechanical behaviours and electrical signals in articular cartilage under unconfined compression — a triphasic paradigm , 2002 .

[31]  Y K Liu,et al.  A three-dimensional nonlinear finite element model of lumbar intervertebral joint in torsion. , 1987, Journal of biomechanical engineering.

[32]  W C Hutton,et al.  The effect of fluid loss on the viscoelastic behavior of the lumbar intervertebral disc in compression. , 1998, Journal of biomechanical engineering.

[33]  H. Cheung,et al.  New insight into deformation-dependent hydraulic permeability of gels and cartilage, and dynamic behavior of agarose gels in confined compression. , 2003, Journal of biomechanics.

[34]  A Shirazi-Adl,et al.  A finite element study of a lumbar motion segment subjected to pure sagittal plane moments. , 1986, Journal of biomechanics.

[35]  Hutton Wc,et al.  Do bending, twisting, and diurnal fluid changes in the disc affect the propensity to prolapse? A viscoelastic finite element model , 1996 .

[36]  B R Simon,et al.  Structural models for human spinal motion segments based on a poroelastic view of the intervertebral disk. , 1985, Journal of biomechanical engineering.

[37]  Gerard A. Ateshian,et al.  Interstitial Fluid Pressurization During Confined Compression Cyclical Loading of Articular Cartilage , 2000, Annals of Biomedical Engineering.

[38]  A. Shirazi-Adl,et al.  Dynamics of Human Lumbar Intervertebral Joints: Experimental and Finite‐Element Investigations , 1992, Spine.

[39]  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.

[40]  W M Lai,et al.  A triphasic theory for the swelling and deformation behaviors of articular cartilage. , 1991, Journal of biomechanical engineering.