Biological soft tissues and cells may be subjected to mechanical as well as chemical (osmotic) loading under their natural physiological environment or various experimental conditions. The interaction of mechanical and chemical effects may be very significant under some of these conditions, yet the highly nonlinear nature of the set of governing equations describing these mechanisms poses a challenge for the modeling of such phenomena. This study formulated and implemented a finite element algorithm for analyzing mechanochemical events in neutral deformable porous media under finite deformation. The algorithm employed the framework of mixture theory to model the porous permeable solid matrix and interstitial fluid, where the fluid consists of a mixture of solvent and solute. A special emphasis was placed on solute-solid matrix interactions, such as solute exclusion from a fraction of the matrix pore space (solubility) and frictional momentum exchange that produces solute hindrance and pumping under certain dynamic loading conditions. The finite element formulation implemented full coupling of mechanical and chemical effects, providing a framework where material properties and response functions may depend on solid matrix strain as well as solute concentration. The implementation was validated using selected canonical problems for which analytical or alternative numerical solutions exist. This finite element code includes a number of unique features that enhance the modeling of mechanochemical phenomena in biological tissues. The code is available in the public domain, open source finite element program FEBio (http://mrl.sci.utah.edu/software).
Issue Section:
Research Papers
References
1.
Crank
, J.
, 1979, The Mathematics of Diffusion
, 2nd ed., Clarendon Press
, Oxford
.2.
Deen
, W. M.
, 1987, “Hindered Transport of Large Molecules in Liquid-Filled Pores
,” AIChE J.
, 33
(9
), pp. 1409
–1425
.3.
Ogston
, A. G.
, and Phelps
, C. F.
, 1961, “The Partition of Solutes Between Buffer Solutions and Solutions Containing Hyaluronic Acid
,” Biochem. J.
, 78
, pp. 827
–833
.4.
Laurent
, T. C.
, and Killander
, J.
, 1963, “A Theory of Gel Filtration and Its Experimental Verification
,” J. Chromatogr.
, 14
, pp. 317
–330
.5.
Albro
, M. B.
, Chahine
, N. O.
, Caligaris
, M.
, Wei
, V. I.
, Likhitpanichkul
, M.
, Ng
, K. W.
, Hung
, C. T.
, and Ateshian
, G. A.
, 2007, “Osmotic Loading of Spherical Gels: A Biomimetic Study of Hindered Transport in the Cell Protoplasm
,” ASME J. Biomech. Eng.
, 129
(4
), pp. 503
–510
.6.
Quinn
, T. M.
, Morel
, V.
, and Meister
, J. J.
, 2001, “Static Compression of Articular Cartilage Can Reduce Solute Diffusivity and Partitioning: Implications for the Chondrocyte Biological Response
,” J. Biomech.
, 34
(11
), pp. 1463
–1469
.7.
Gu
, W. Y.
, Yao
, H.
, Vega
, A. L.
, and Flagler
, D.
, 2004, “Diffusivity of Ions in Agarose Gels and Intervertebral Disc: Effect of Porosity
,” Ann. Biomed. Eng.
, 32
(12
), pp. 1710
–1717
.8.
Lai
, W. M.
, Hou
, J. S.
, and Mow
, V. C.
, 1991, “A Triphasic Theory for the Swelling and Deformation Behaviors of Articular Cartilage
,” ASME J. Biomech. Eng.
, 113
(3
), pp. 245
–258
.9.
Huyghe
, J. M.
, and Janssen
, J. D.
, 1997, “Quadriphasic Mechanics of Swelling Incompressible Porous Media
,” Int. Eng. Sci.
, 35
(8
), pp. 793
–802
.10.
Gu
, W. Y.
, Lai
, W. M.
, and Mow
, V. C.
, 1998, “A Mixture Theory for Charged-Hydrated Soft Tissues Containing Multi-electrolytes: Passive Transport and Swelling Behaviors
,” ASME J. Biomech. Eng.
, 120
(2
), pp. 169
–180
.11.
Gu
, W. Y.
, Lai
, W. M.
, and Mow
, V. C.
, 1993, “Transport of Fluid and Ions Through a Porous-Permeable Charged-Hydrated Tissue, and Streaming Potential Data on Normal Bovine Articular Cartilage
,” J. Biomech.
, 26
(6
), pp. 709
–723
.12.
Lai
, W. M.
, Mow
, V. C.
, Sun
, D. D.
, and Ateshian
, G. A.
, 2000, “On the Electric Potentials Inside a Charged Soft Hydrated Biological Tissue: Streaming Potential Versus Diffusion Potential
,” ASME J. Biomech. Eng.
, 122
(4
), pp. 336
–346
.13.
Huyghe
, J. M.
, Houben
, G. B.
, Drost
, M. R.
, and van Donkelaar
, C. C.
, 2004, “An Ionised/Nonionised Dual Porosity Model of Intervertebral Disc Tissue: Experimental Quantification of Parameters
,” Biomech. Model Mechanobiol.
, 2
(4
), pp. 3
–19
.14.
Gu
, W. Y.
, Lai
, W. M.
, and Mow
, V. C.
, 1997, “A Triphasic Analysis of Negative Osmotic Flows Through Charged Hydrated Soft Tissues
,” J. Biomech.
, 30
(1
), pp. 71
–78
.15.
Simon
, B. R.
, Liable
, J. P.
, Pflaster
, D.
, Yuan
, Y.
, and Krag
, M. H.
, 1996, “A Poroelastic Finite Element Formulation Including Transport and Swelling in Soft Tissue Structures
,” ASME J. Biomech. Eng.
, 118
(1
), pp. 1
–9
.16.
Frijns
, A. J. H.
, Huyghe
, J. M.
, and Janssen
, J. D.
, 1997, “Validation of the Quadriphasic Mixture Theory for Intervertebral Disc Tissue
,” Int. J. Eng. Sci.
, 35
(15
), pp. 1419
–1429
.17.
Sun
, D. N.
, Gu
, W. Y.
, Guo
, X. E.
, Lai
, W. M.
, and Mow
, V. C.
, 1999, “A Mixed Finite Element Formulation of Triphasic Mechano-electrochemical Theory for Charged, Hydrated Biological Soft Tissues
,” Int. J. Numer. Methods Eng.
, 45
(10
), pp. 1375
–1402
.18.
Kaasschieter
, E. F.
, Frijns
, A. J. H.
, and Huyghe
, J. M.
, 2003, “Mixed Finite Element Modelling of Cartilaginous Tissues
,” Math. Comput. Simul.
, 61
(3–6
), pp. 549
–560
.19.
Yao
, H.
, and Gu
, W. Y.
, 2007, “Three-Dimensional Inhomogeneous Triphasic Finite-Element Analysis of Physical Signals and Solute Transport in Human Intervertebral Disc Under Axial Compression
,” J. Biomech.
, 40
(9
), pp. 2071
–2077
.20.
Magnier
, C.
, Boiron
, O.
, Wendling-Mansuy
, S.
, Chabrand
, P.
, and Deplano
, V.
, 2009, “Nutrient Distribution and Metabolism in the Intervertebral Disc in the Unloaded State: A Parametric Study
,” J. Biomech.
, 42
(2
), pp. 100
–108
.21.
van Loon
, R.
, Huyghe
, J. M.
, Wijlaars
, M. W.
, and Baaijens
, F. P. T.
, 2003, “3D FE Implementation of an Incompressible Quadriphasic Mixture Model
,” Int. J. Numer. Methods Eng.
, 57
(9
), pp. 1243
–1258
.22.
Wu
, J. Z.
, and Herzog
, W.
, 2002, “Simulating the Swelling and Deformation Behaviour in Soft Tissues Using a Convective Thermal Analogy
,” Biomed. Eng. Online
, 1
, p. 8
.23.
Azeloglu
, E. U.
, Albro
, M. B.
, Thimmappa
, V. A.
, Ateshian
, G. A.
, and Costa
, K. D.
, 2008, “Heterogeneous Transmural Proteoglycan Distribution Provides a Mechanism for Regulating Residual Stresses in the Aorta
,” Am. J. Physiol. Heart Circ. Physiol.
, 294
(3
), pp. H1197
–H1205
.24.
Ateshian
, G. A.
, Rajan
, V.
, Chahine
, N. O.
, Canal
, C. E.
, and Hung
, C. T.
, 2009, “Modeling the Matrix of Articular Cartilage Using a Continuous Fiber Angular Distribution Predicts Many Observed Phenomena
,” ASME J. Biomech. Eng.
, 131
(6
), p. 061003
.25.
Sengers
, B. G.
, Oomens
, C. W.
, and Baaijens
, F. P.
, 2004, “An Integrated Finite-Element Approach to Mechanics, Transport and Biosynthesis in Tissue Engineering
,” ASME J. Biomech. Eng.
, 126
(1
), pp. 82
–91
.26.
Steck
, R.
, Niederer
, P.
, and Knothe Tate
, M. L.
, 2003, “A Finite Element Analysis for the Prediction of Load-Induced Fluid Flow and Mechanochemical Transduction in Bone
,” J. Theor. Biol.
, 220
(2
), pp. 249
–259
.27.
Zhang
, L.
, and Szeri
, A.
, 2005, “Transport of Neutral Solute in Articular Cartilage: Effects of Loading and Particle Size
,” Proc. R. Soc. London
, 461
(2059), pp. 2021
–2042
.28.
Basser
, P. J.
, 1995, “Inferring Microstructural Features and the Physiological State of Tissues From Diffusion-Weighted Images
,” NMR Biomed.
, 8
(7–8
), pp. 333
–344
.29.
Basser
, P. J.
, and Pierpaoli
, C.
, 1996, “Microstructural and Physiological Features of Tissues Elucidated by Quantitative-Diffusion-Tensor MRI
,” J. Magn. Reson., Ser. B
111
(3
), pp. 209
–219
.30.
Leddy
, H. A.
, Haider
, M. A.
, and Guilak
, F.
, 2006, “Diffusional Anisotropy in Collagenous Tissues: Fluorescence Imaging of Continuous Point Photobleaching
,” Biophys. J.
, 91
(1
), pp. 311
–316
.31.
Ellegood
, J.
, McKay
, R. T.
, Hanstock
, C. C.
, and Beaulieu
, C.
, 2007, “Anisotropic Diffusion of Metabolites in Peripheral Nerve Using Diffusion Weighted Magnetic Resonance Spectroscopy at Ultra-high Field
,” J. Magn. Reson.
, 184
(1
), pp. 20
–28
.32.
Mansour
, J. M.
, and Mow
, V. C.
, 1976, “The Permeability of Articular Cartilage Under Compressive Strain and at High Pressures
,” J. Bone Joint Surg. AM.
, 58
(4
), pp. 509
–516
.33.
Quinn
, T. M.
, Kocian
, P.
, and Meister
, J. J.
, 2000, “Static Compression Is Associated With Decreased Diffusivity of Dextrans in Cartilage Explants
,” Arch. Biochem. Biophys.
, 384
(2
), pp. 327
–34
.34.
Gu
, W. Y.
, Yao
, H.
, Huang
, C. Y.
, and Cheung
, H. S.
, 2003, “New Insight Into Deformation-Dependent Hydraulic Permeability of Gels and Cartilage, and Dynamic Behavior of Agarose Gels in Confined Compression
,” J. Appl. Biomech.
, 36
(4
), pp. 593
–598
.35.
Ateshian
, G. A.
, and Weiss
, J. A.
, 2010, “Anisotropic Hydraulic Permeability Under Finite Deformation
,” ASME J. Biomech. Eng.
, 132
(11
), p. 111004
.36.
Mackie
, J. S.
, and Meares
, P.
, 1955, “The Diffusion of Electrolytes in a Cation-Exchange Resin Membrane. I. Theoretical
,” Proc. R. Soc. London
, 232
(1191), pp. 498
–509
.37.
Evans
, R. C.
, and Quinn
, T. M.
, 2006, “Solute Convection in Dynamically Compressed Cartilage
,” J. Appl. Biomech.
, 39
(6
), pp. 1048
–1055
.38.
Mauck
, R. L.
, Hung
, C. T.
, and Ateshian
, G. A.
, 2003, “Modeling of Neutral Solute Transport in a Dynamically Loaded Porous Permeable Gel: Implications for Articular Cartilage Biosynthesis and Tissue Engineering
,” ASME J. Biomech. Eng.
, 125
(5
), pp. 602
–614
.39.
Albro
, M. B.
, Chahine
, N. O.
, Li
, R.
, Yeager
, K.
, Hung
, C. T.
, and Ateshian
, G. A.
, 2008, “Dynamic Loading of Deformable Porous Media Can Induce Active Solute Transport
,” J. Biomech.
, 41
(15
), pp. 3152
–3157
.40.
Albro
, M. B.
, Li
, R.
, Banerjee
, R. E.
, Hung
, C. T.
, and Ateshian
, G. A.
, 2010, “Validation of Theoretical Framework Explaining Active Solute Uptake in Dynamically Loaded Porous Media
,” J. Biomech.
, 43
(12
), pp. 2267
–2273
.41.
Ateshian
, G. A.
, Likhitpanichkul
, M.
, and Hung
, C. T.
, 2006, “A Mixture Theory Analysis for Passive Transport in Osmotic Loading of Cells
,” J. Biomech.
, 39
(3
), pp. 464
–475
.42.
Weiss
, T. F.
, 1996, Cellular Biophysics
, MIT Press
, Cambridge, MA
.43.
Albro
, M. B.
, Petersen
, L. E.
, Li
, R.
, Hung
, C. T.
, and Ateshian
, G. A.
, 2009, “Influence of the Partitioning of Osmolytes by the Cytoplasm on the Passive Response of Cells to Osmotic Loading
,” Biophys. J.
, 97
(11
), pp. 2886
–2893
.44.
Lazzara
, M. J.
, and Deen
, W. M.
, 2004, “Effects of Concentration on the Partitioning of Macromolecule Mixtures in Agarose Gels
,” J. Colloid Interface Sci.
, 272
(2
), pp. 288
–297
.45.
Albro
, M. B.
, Rajan
, V.
, Li
, R.
, Hung
, C. T.
, and Ateshian
, G. A.
, 2009, “Characterization of the Concentration-Dependence of Solute Diffusivity and Partitioning in a Model Dextran-Agarose Transport System
,” Cell Mol. Bioeng.
, 2
(3
), pp. 295
–305
.46.
Tinoco
, I.
, Sauer
, K.
, and Wang
, J. C.
, 1995, Physical Chemistry: Principles and Applications in Bbiological Sciences
, 3rd ed., Prentice Hall
, Englewood Cliffs, NJ
.47.
McNaught
, A. D.
, and Wilkinson
, A.
, 1997, Compendium of Chemical Terminology: IUPAC Recommendations
, 2nd ed., Blackwell Science
, Oxford
.48.
Truesdell
, C.
, and Toupin
, R.
, 1960, The Classical Field Theories, Vol. III of Handbuch der Physik
, Springer
, Heidelberg
.49.
Bowen
, R.
, 1976, Theory of Mixtures Vol. 3 of Continuum Physics
, Academic
, New York
.50.
Yao
, H.
, and Gu
, W. Y.
, 2006, “Physical Signals and Solute Transport in Human Intervertebral Disc During Compressive Stress Relaxation: 3D Finite Element Analysis
,” Biorheology
, 43
(3–4
), pp. 323
–35
.51.
Zhang
, L.
, Gardiner
, B. S.
, Smith
, D. W.
, Pivonka
, P.
, and Grodzinsky
, A.
, 2007, “The Effect of Cyclic Deformation and Solute Binding on Solute Transport in Cartilage
,” Arch. Biochem. Biophys.
, 457
(1
), pp. 47
–56
.52.
Ateshian
, G. A.
, 2007, “On the Theory of Reactive Mixtures for Modeling Biological Growth
,” Biomech. Model. Mechanobiol.
, 6
(6
), pp. 423
–445
.53.
Lai
, W. M.
, Mow
, V. C.
, and Roth
, V.
, 1981, “Effects of Nonlinear Strain-Dependent Permeability and Rate of Compression on the Stress Behavior of Articular Cartilage
,” ASME J. Biomech. Eng.
, 103
(2
), pp. 61
–66
.54.
Holmes
, M. H.
, and Mow
, V. C.
, 1990, “The Nonlinear Characteristics of Soft Gels and Hydrated Connective Tissues in Ultrafiltration
,” J. Biomech.
, 23
(11
), pp. 1145
–1156
.55.
Chen
, A. C.
, Bae
, W. C.
, Schinagl
, R. M.
, and Sah
, R. L.
, 2001, “Depth- and Strain-Dependent Mechanical and Electromechanical Properties of Full-Thickness Bovine Articular Cartilage in Confined Compression
,” J. Biomech.
, 34
(1
), pp. 1
–12
.56.
Eringen
, A.
, and Ingram
, J.
, 1965, “Continuum Theory of Chemically Reacting Media - 1[prime]
,” Int. J. Comput. Eng. Sci.
, 3
, pp. 197
–212
.57.
Katzir-Katchalsky
, A.
, and Curran
, P. F.
, 1965, Nonequilibrium Thermodynamics in Biophysics
, Harvard University Press
, Cambridge
.58.
Bonet
, J.
, and Wood
, R. D.
, 1997, Nonlinear Continuum Mechanics for Finite Element Analysis
, Cambridge University Press
, New York
.59.
Marsden
, J. E.
, and Hughes
, T. J. R.
, 1994, Mathematical Foundations of Elasticity
, Dover
, New York
.60.
Albro
, M. B.
, Banerjee
, R. E.
, Li
, R.
, Oungoulian
, S. R.
, Chen
, B.
, Del Palomar
, A. P.
, Hung
, C. T.
, and Ateshian
, G. A.
, 2011, “Dynamic Loading of Immature Epiphyseal Cartilage Pumps Nutrients Out of Vascular Canals
,” J. Biomech
, 44
(9
), pp. 1654
–1659
.61.
Armstrong
, C. G.
, Lai
, W. M.
, and Mow
, V. C.
, 1984, “An Analysis of the Unconfined Compression of Articular Cartilage
,” ASME J. Biomech. Eng.
, 106
(2
), pp. 165
–173
.62.
Lucio
, A. D.
, Santos
, R. A.
, and Mesquita
, O. N.
, 2003, “Measurements and Modeling of Water Transport and Osmoregulation in a Single Kidney Cell Using Optical Tweezers and Videomicroscopy
,” Phys. Rev. E
, 68
, p. 041906
.63.
Mow
, V. C.
, Kuei
, S. C.
, Lai
, W. M.
, and Armstrong
, C. G.
, 1980, “Biphasic Creep and Stress Relaxation of Articular Cartilage in Compression: Theory and Experiments
,” ASME J. Biomech. Eng.
, 102
(1
), pp. 73
–84
.64.
Curnier
, A.
, He
, Q.-C.
, and Zysset
, P.
, 1995, “Conewise Linear Elastic Materials
,” J. Elasticity
, 37
(1
), pp. 1
–38
Copyright © 2011
by American Society of Mechanical Engineers
You do not currently have access to this content.
