We present a mathematical model for the reorientation of fibers in a soft, fiber-reinforced, fluid-saturated porous medium describing a hypothetical biological tissue. We consider two types of remodeling that, at different scales, concur in determining the structural reorganization of the solid phase of the medium: one pertains to the development of plastic-like distortions, which are introduced as the macroscopic manifestation of processes occurring at lower scales, but not resolved explicitly; the other one, originating at the mesoscale, concerns the capability of the fibers of reorienting in the extracellular matrix. This latter form of remodeling is studied as a Langevin-like process in which two main agencies are recognized: a drift term, which is given by the effects that the deformation of the extracellular matrix exerts on the fibers, and a noise term, which accounts for the interactions among nearby fibers. We employ the framework of the Principle of Virtual Power to present our model, in which we “free” the kinematics of the system from the constraints of isochoricity of remodeling and incompressibility of the mixture by appending the Chetaev forms of these constraints to the Principle of Virtual Power. Then, we specialize our model to articular cartilage. We do this for comparing our results with some experimental curves describing the distribution of collagen fibers in a sample of articular cartilage, and for studying, through the simulation of a uniaxial compression test, the interplay between fluid flow and the two aforementioned forms of remodeling. Our main result is the establishment of a framework that captures effectively the mechanical coupling between the fiber distribution specific to a given medium and the mechanical stimuli exerted on such distribution. Quite differently from other works on this subject, our framework is capable of accounting for the stochastic effects of the fibers on the overall evolution of the considered tissue.
GundaSGiammariniARamírez-TorresA, et al. Fractionalization of Forchheimer’s correction to Darcy’s law in porous media in large deformations. Math Mech Solids. Epub ahead of print 2August2024. DOI: 10.1177/10812865241252577.
2.
AteshianG.On the theory of reactive mixtures for modeling biological growth. Biomech Model Mechan2007; 6(6): 423–445.
3.
HolmesMMowV.The nonlinear characteristics of soft gels and hydrated connective tissues in ultrafiltration. J Biomech1990; 23: 1145–1156.
4.
LoretBSimõesF. A framework for deformation, generalized diffusion, mass transfer and growth in multi-species multi-phase biological tissues. Eur J Mech A Solids2005; 24: 757–781.
5.
AmbrosiDPreziosiLVitaleG.The insight of mixtures theory for growth and remodeling. Z Angew Math Phys2010; 61: 177–191.
6.
GrilloAFedericoSWittumG.Growth, mass transfer, and remodeling in fiber-reinforced, multi-constituent materials. Int J Non Linear Mech2012; 47: 388–401.
HolzapfelGGasserTOgdenR.A new constitutive framework for arterial wall mechanics and a comparative study of material models. J Elasticity2000; 61(1–3): 1–48.
9.
PandolfiAManganielloF.A model for the human cornea: constitutive formulation and numerical analysis. Biomech Model Mechan2006; 5(4): 237–246.
10.
PandolfiAHolzapfelGA.Three-dimensional modeling and computational analysis of the human cornea considering distributed collagen fibril orientations. J Biomech Eng2008; 130(6): 061006.
11.
GentlemanELayANDickersonDA, et al. Mechanical characterization of collagen fibers and scaffolds for tissue engineering. Biomaterials2003; 24(21): 3805–3813.
12.
MerodioJOgdenR.Mechanical response of fiber-reinforced incompressible non-linearly elastic solids. Int J Non Linear Mech2005; 40(2–3): 213–227.
13.
GasserTOgdenRHolzapfelG.Hyperelastic modelling of arterial layers with distributed collagen fibre orientations. J R Soc Interface2006; 3: 15–35.
14.
FedericoSGrilloALa RosaG, et al. A transversely isotropic, transversely homogeneous microstructural-statistical model of articular cartilage. J Biomech2005; 38: 2008–2018.
SpencerASoldatosK.Finite deformations of fibre-reinforced elastic solids with fibre bending stiffness. Int J Non Linear Mech2007; 42(2): 355–368.
17.
HaritonIdeBottonGGasserT, et al. Stress-driven collagen fiber remodeling in arterial walls. Biomech Model Mechan2007; 6(3): 163–175.
18.
HolzapfelGAOgdenRW.Constitutive modelling of arteries. Proc R Soc A Math Phys Eng Sci2010; 466(2118): 1551–1597.
19.
HolzapfelGANiestrawskaJAOgdenRW, et al. Modelling non-symmetric collagen fibre dispersion in arterial walls. J Roy Soc Interface2015; 12(106): 20150188.
20.
HolzapfelGOgdenR.On fiber dispersion models: exclusion of compressed fibers and spurious model comparisons. J Elasticity2017; 129(1–2): 49–68.
21.
MelnikADa RochaHBGorielyA.On the modeling of fiber dispersion in fiber-reinforced elastic materials. Int J Non Linear Mech2015; 75: 92–106.
22.
ItskovMAkselN.A class of orthotropic and transversely isotropic hyperelastic constitutive models based on a polyconvex strain energy function. Int J Solids Struct2004; 41(14): 3833–3848.
23.
NolanDGowerADestradeM, et al. A robust anisotropic hyperelastic formulation for the modelling of soft tissue. J Mech Behav Biomed Mater2014; 39: 48–60.
24.
LanirY.Constitutive equations for fibrous connective tissues. J Biomech1983; 16: 1–12.
25.
FedericoSHerzogW.On the anisotropy and inhomogeneity of permeability in articular cartilage. Biomech Model Mechan2008; 7: 367–378.
26.
FedericoSHerzogW.Towards an analytical model of soft biological tissues. J Biomech2008; 41: 3309–3313.
27.
FedericoSGasserT.Non-linear elasticity of biological tissues with statistical fibre orientation. J R Soc Interface2010; 7: 955–966.
28.
FedericoSGrilloA.Elasticity and permeability of porous fibre-reinforced materials under large deformations. Mech Mater2012; 44: 58–71.
29.
TomicAGrilloAFedericoS.Poroelastic materials reinforced by statistically oriented fibres — numerical implementation and application to articular cartilage. IMA J Appl Math2014; 79: 1027–1059.
30.
GizziAPandolfiAVastaM.A generalized statistical approach for modeling fiber-reinforced materials. J Eng Math2017; 109(1): 211–226.
31.
GrilloACarfagnaMFedericoS.An Allen-Cahn approach to the remodelling of fibre-reinforced anisotropic materials. J Eng Math2018; 109(1): 139–172.
32.
CrevacoreEDi StefanoSGrilloA.Coupling among deformation, fluid flow, structural reorganisation and fibre reorientation in fibre-reinforced, transversely isotropic biological tissues. Int J Non Linear Mech2019; 111: 1–13.
33.
MalekiMHashlamounKHerzogW, et al. Effect of structural distortions on articular cartilage permeability under large deformations. Biomech Model Mechan2019; 19(1): 317–334.
34.
VastaMGizziAPandolfiA.On three- and two-dimensional fiber distributed models of biological tissues. Probabilist Eng Mech2014; 37: 170–179.
35.
HolzapfelGOgdenRSherifovaS.On fibre dispersion modelling of soft biological tissues: a review. Proc R Soc A Math Phys Eng Sci2019; 475(2224): 1–22.
36.
QuinnTDierickxPGrodzinskyA.Glycosaminoglycan network geometry may contribute to anisotropic hydraulic permeability in cartilage under compression. J Biomech2001; 34(11): 1483–1490.
37.
PierceDMRickenTHolzapfelGA.A hyperelastic biphasic fibre-reinforced model of articular cartilage considering distributed collagen fibre orientations: continuum basis, computational aspects and applications. Comput Method Biomec2013; 16(12): 1344–1361.
38.
HashlamounKGrilloAFedericoS.Efficient evaluation of the material response of tissues reinforced by statistically oriented fibres. Z Angew Math Phys2016; 67: 113–145.
39.
MaroudasABulloughP.Permeability of articular cartilage. Nature1968; 219(5160): 1260–1261.
40.
EpsteinMMauginG.Thermomechanics of volumetric growth in uniform bodies. Int J Plasticity2000; 16(7–8): 951–978.
41.
GarikipatiKArrudaEGroshK, et al. A continuum treatment of growth in biological tissue: the coupling of mass transport and mechanics. J Mech Phys Solids2004; 52: 1595–1625.
42.
WilsonWDriessenNvan DonkelaarCC, et al. Prediction of collagen orientation in articular cartilage by a collagen remodeling algorithm. Osteoarthr Cartilage2006; 14: 1196–1202.
43.
AmbrosiDAteshianGArrudaE, et al. Perspectives on biological growth and remodeling. J Mech Phys Solids2011; 59(4): 863–883.
44.
LinWJIafratiMDPeattieRA, et al. Growth and remodeling with application to abdominal aortic aneurysms. J Eng Math2017; 109(1): 113–137.
45.
SoldatosKP.Foundation of polar linear elasticity for fibre-reinforced materials. J Elasticity2013; 114(2): 155–178.
46.
SoldatosKPShariffMHBMMerodioJ.On the constitution of polar fiber-reinforced materials. Mech Adv Mater Struc2020; 28(21): 2255–2266.
47.
DriessenNWilsonWBoutenC, et al. A computational model for collagen fibre remodelling in the arterial wall. J Theor Biol2004; 226: 53–64.
48.
MenzelA.Modelling of anisotropic growth in biological tissues — a new approach and computational aspects. Biomech Model Mechan2005; 3: 147–171.
49.
MenzelA.A fibre reorientation model for orthotropic multiplicative growth: configurational driving stresses, kinematics-based reorientation and algorithmic aspects. Biomech Model Mechan2007; 6(5): 303–320.
50.
GizziACyronCJFalcinelliC, et al. Evolution of fiber distributions in homogenized constrained mixture models of soft tissue growth and remodeling: uniaxial loading. J Mech Phys Solids2024; 183: 105491.
51.
GrilloADi StefanoS.A formulation of volumetric growth as a mechanical problem subjected to non-holonomic and rheonomic constraint. Math Mech Solids2023; 28: 2215–2241.
52.
GrilloADi StefanoS. Addendum to “a formulation of volumetric growth as a mechanical problem subjected to non-holonomic and rheonomic constraint”. Math Mech Solids (2024); 29(1): 62–70.
53.
MartinsJPiresESalvadoR, et al. A numerical model of passive and active behavior of skeletal muscles. Comput Method Appl M1998; 151(3–4): 419–433.
54.
BlemkerSSPinskyPMDelpSL.A 3D model of muscle reveals the causes of nonuniform strains in the biceps brachii. J Biomech2005; 38(4): 657–665.
55.
GiantesioGMusestiARiccobelliD.A comparison between active strain and active stress in transversely isotropic hyperelastic materials. J Elasticity2018; 137(1): 63–82.
56.
OlssonTKlarbringA.Residual stresses in soft tissue as a consequence of growth and remodeling: application to an arterial geometry. Eur J Mech A/Solids2008; 27(6): 959–974.
57.
WilsonWHuygheJvan DonkelaarCC.A composition-based cartilage model for the assessment of compositional changes during cartilage damage and adaptation. Osteoarthr Cartilage2006; 14: 554–560.
58.
PierceDMRickenTHolzapfelGA.A hyperelastic biphasic fibre-reinforced model of articular cartilage considering distributed collagen fibre orientations: continuum basis, computational aspects and applications. Comput Method Biomec2013; 16(12): 1344–1361.
59.
LiKOgdenRHolzapfelGA.Computational method for excluding fibers under compression in modeling soft fibrous solids. Eur J Mech A/Solids2016; 57: 178–193.
60.
GarciaDZyssetPCharleboisM, et al. A three-dimensional elastic plastic damage constitutive law for bone tissue. Biomech Model Mechan2009; 8(2): 149–165.
61.
Di StefanoSCarfagnaMKnodelMM, et al. Anelastic reorganisation of fibre-reinforced biological tissues. Comput Vis Sci2019; 20(3–6): 95–109.
62.
MićunovićM.Thermomechanics of viscoplasticity. New York: Springer, 2009.
63.
DiCarloAQuiligottiS.Growth and balance. Mech Res Commun2002; 29(6): 449–456.
64.
Di StefanoSRamírez-TorresAPentaR, et al. Self-influenced growth through evolving material inhomogeneities. Int J Non Linear Mech2018; 106: 174–187.
65.
CiancioVDolfinMFrancavigliaM, et al. Uniform materials and the multiplicative decomposition of the deformation gradient in finite elasto-plasticity. J Non-Equilib Thermodyn2008; 33(3): 199–234.
66.
MarsdenJHughesT.Mathematical foundations of elasticity. Mineola, NY: Dover Publications, Inc., 1983.
67.
QuiligottiSMauginGdell’IsolaF.An Eshelbian approach to the nonlinear mechanics of constrained solid-fluid mixtures. Acta Mech2003; 160: 45–60.
68.
CermelliPFriedESellersS.Configurational stress, yield and flow in rate-independent plasticity. Proc R Soc Math Phys Eng Sci2001; 457: 1447–1467.
69.
GiorgioIdell’IsolaFAndreausU, et al. An orthotropic continuum model with substructure evolution for describing bone remodeling: an interpretation of the primary mechanism behind Wolff’s law. Biomech Model Mech2023; 22(6): 2135–2152.
70.
GantmacherF.Lectures in analytical mechanics. Moscow, Russia: MIR Publishers, 1975.
71.
LlibreJRamírezRSadovskaiaN.Nonlinear dynamics. a new approach to the vakonomic mechanics. Nonlinear Dynam2014; 78: 2219–2247.
72.
ParsL.A treatise on analytical dynamics. London: Heinemann, 1965.
73.
HassanizadehSM.Derivation of basic equations of mass transport in porous media, part 1: macroscopic balance laws. Adv Water Resour1986; 9(4): 196–206.
74.
HassanizadehS.Derivation of basic equations of mass transport in porous media, part 2: generalized Darcy’s and Fick’s laws. Adv Water Resour1986; 9: 207–222.
75.
BearJBachmatY.Introduction to modeling of transport phenomena in porous media. Dordrecht: Kluwer, 1990.
76.
BennethumLMuradMCushmanJ.Macroscale thermodynamics and the chemical potential for swelling porous media. Transport Porous Med2000; 39(2): 187–225.
77.
AteshianGHumphreyJ.Continuum mixture models of biological growth and remodeling: past successes and future opportunities. Annu Rev Biomed Eng2012; 14(1): 97–111.
78.
HolzapfelGAOgdenRW.An arterial constitutive model accounting for collagen content and cross-linking. J Mech Phys Solids2020; 136: 103682.
79.
MooEKTanskaPFedericoS, et al. Collagen fibres determine the crack morphology in articular cartilage. Acta Biomater2021; 126: 301–314.
80.
MauginGAMetrikineAV.Mechanics of generalized continua: one hundred years after the Cosserats. New York: Springer, 2010.
GrilloAWittumGTomicA, et al. Remodelling in statistically oriented fibre-reinforced materials and biological tissues. Math Mech Solids2015; 20(9): 1107–1129.
83.
LoyNPreziosiL.A statistical mechanics approach to describe cell reorientation under stretch. B Math Biol2023; 85(7): 60.
84.
FedericoSHerzogW.On the permeability of fibre-reinforced porous materials. Int J Solids Struct2008; 45: 2160–2172.
GrilloACarfagnaMFedericoS.Non-Darcian flow in fibre-reinforced biological tissues. Meccanica2017; 52: 3299–3320.
87.
CarfagnaMGrilloA.The spherical design algorithm in the numerical simulation of biological tissues with statistical fibre-reinforcement. Comput Vision Sci2017; 18: 157–184.
88.
GrilloADi StefanoSRamírez-TorresA, et al. A study of growth and remodeling in isotropic tissues, based on the Anand-Aslan-Chester theory of strain-gradient plasticity. GAMM-Mitt2019; 42(4): e201900015.
89.
DestradeMDonaldBMMurphyJG, et al. At least three invariants are necessary to model the mechanical response of incompressible, transversely isotropic materials. Comput Mech2013; 52(4): 959–969.
90.
HamdaouiMEMerodioJOgdenRW, et al. Finite elastic deformations of transversely isotropic circular cylindrical tubes. Int J Solids Struct2014; 51(5): 1188–1196.
91.
SimoJ. A framework for finite strain elastoplasticity based on maximum plastic dissipation and the multiplicative decomposition: part i. Continuum formulation. Comput Method Appl M1988; 66(2): 199–219.
92.
GreenAENaghdiPM.A general theory of an elastic-plastic continuum. Arch Ration Mech An1965; 18(4): 251–281.
93.
WuJZHerzogW.Finite element simulation of location- and time-dependent mechanical behavior of chondrocytes in unconfined compression tests. Ann Biomed Eng2000; 28(3): 318–330.
GrilloADi StefanoS.An a posteriori approach to the mechanics of volumetric growth. Math Mech Complex Syst2023; 11(1): 57–86.
96.
Pericak-SpectorKASpectorSJ.On the representation theorem for linear, isotropic tensor functions. J Elasticity1995; 39(2): 181–185.
97.
HolzapfelGAOgdenRW.Constitutive modelling of passive myocardium: a structurally based framework for material characterization. Philos Trans R Soc Math Phys Eng Sci2009; 367(1902): 3445–3475.
98.
RiskenH.The Fokker-Planck equation: methods of solution and applications. Berlin: Springer, 1996.
99.
JordanRKinderlehrerDOttoF.The variational formulation of the Fokker–Planck equation. SIAM J Math Anal1998; 29(1): 1–17.
100.
AspdenRMHukinsDW.Collagen organization in articular cartilage, determined by x-ray diffraction, and its relationship to tissue function. Proc R Soc Lond B Biol Sci1981; 212(1188): 299–304.
101.
ClarkABarclayLMatyasJ, et al. In situ chondrocyte deformation with physiological compression of the feline patellofemoral joint. J Biomech2003; 36(4): 553–568.
102.
MollenhauerJAurichMMuehlemanC, et al. X-ray diffraction of the molecular substructure of human articular cartilage. Connect Tissue Res2003; 44(5): 201–207.
103.
MansourJ. Biomechanics of cartilage. In: OatisCA (ed.) Kinesiology: the mechanics and pathomechanics of human movement. Baltimore, MD: Wolters Kluwer Health, 2013, pp. 69–83.
104.
HedlundHMengarelli-WidholmSReinholtFP, et al. Stereologic studies on collagen in bovine articular cartilage. APMIS1993; 101(1–6): 133–140.
105.
ChenSFalcovitzYSchneidermanR, et al. Depth-dependent compressive properties of normal aged human femoral head articular cartilage: relationship to fixed charge density. Osteoarthr Cartilage2001; 9(6): 561–569.
HardinRSloaneN.Mclaren’s improved snub cube and other new spherical designs in three dimensions. Discrete Comput Geom1996; 15: 429–441.
108.
MarkowichPAVillaniC.On the trend to equilibrium for the Fokker-Planck equation: an interplay between physics and functional analysis. Mat Contemp2000; 19: 1–31.
109.
MansfieldJCMandaliaVTomsA, et al. Collagen reorganization in cartilage under strain probed by polarization sensitive second harmonic generation microscopy. J Roy Soc Interface2019; 16(150): 20180611.
110.
GiversoCDi StefanoSGrilloA, et al. A three dimensional model of multicellular aggregate compression. Soft Matter2019; 15(48): 10005–10019.
111.
GrilloACarfagnaMFedericoS.The Darcy-Forchheimer law for modelling fluid flow in biological tissues. Theor Appli Mech2014; 41(4): 283–322.
112.
PinsGDHuangEKChristiansenDL, et al. Effects of static axial strain on the tensile properties and failure mechanisms of self-assembled collagen fibers. J Appl Polym Sci1997; 63(11): 1429–1440.
113.
HolmesMH.Finite deformation of soft tissue: analysis of a mixture model in uni-axial compression. J Biomech Eng1986; 108(4): 372–381.
HolzapfelGAOgdenRW.A damage model for collagen fibres with an application to collagenous soft tissues. Proc R Soc A Math Phys Eng Sci2020; 476(2236): 20190821.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
