Loading [MathJax]/jax/output/SVG/jax.js
 

The role of structural viscoelasticity in deformable porous media with incompressible constituents: Applications in biomechanics

  • Received: 24 July 2017 Accepted: 27 December 2017 Published: 01 August 2018
  • MSC : Primary: 35Q74, 35B44; Secondary: 35C10

  • The main goal of this work is to clarify and quantify, by means of mathematical analysis, the role of structural viscoelasticity in the biomechanical response of deformable porous media with incompressible constituents to sudden changes in external applied loads. Models of deformable porous media with incompressible constituents are often utilized to describe the behavior of biological tissues, such as cartilages, bones and engineered tissue scaffolds, where viscoelastic properties may change with age, disease or by design. Here, for the first time, we show that the fluid velocity within the medium could increase tremendously, even up to infinity, should the external applied load experience sudden changes in time and the structural viscoelasticity be too small. In particular, we consider a one-dimensional poro-visco-elastic model for which we derive explicit solutions in the cases where the external applied load is characterized by a step pulse or a trapezoidal pulse in time. By means of dimensional analysis, we identify some dimensionless parameters that can aid the design of structural properties and/or experimental conditions as to ensure that the fluid velocity within the medium remains bounded below a certain given threshold, thereby preventing potential tissue damage. The application to confined compression tests for biological tissues is discussed in detail. Interestingly, the loss of viscoelastic tissue properties has been associated with various disease conditions, such as atherosclerosis, Alzheimer’s disease and glaucoma. Thus, the findings of this work may be relevant to many applications in biology and medicine.

    Citation: Maurizio Verri, Giovanna Guidoboni, Lorena Bociu, Riccardo Sacco. The role of structural viscoelasticity in deformable porous media with incompressibleconstituents: Applications in biomechanics[J]. Mathematical Biosciences and Engineering, 2018, 15(4): 933-959. doi: 10.3934/mbe.2018042

    Related Papers:

    [1] Claas Hüter, Mingxuan Lin, Diego Schicchi, Martin Hunkel, Ulrich Prahl, Robert Spatschek . A multiscale perspective on the kinetics of solid state transformations with application to bainite formation. AIMS Materials Science, 2015, 2(4): 319-345. doi: 10.3934/matersci.2015.4.319
    [2] Aarne Pohjonen . Numerical simulations of gradient cooling technique for controlled production of differential microstructure in steel strip or plate. AIMS Materials Science, 2021, 8(6): 852-866. doi: 10.3934/matersci.2021052
    [3] Pavan Hiremath, Gowrishankar M. C., Manjunath Shettar, Sathyashankara Sharma, Jayashree P. K., Suhas Kowshik . Influence of normalizing post carburizing treatment on microstructure, mechanical properties and fracture behavior of low alloy gear steels. AIMS Materials Science, 2021, 8(5): 836-851. doi: 10.3934/matersci.2021051
    [4] Utpal K. Dhar, Md. Farabi Rahman, Mustafa Oguzhan Ayanoglu, Ahammad Abdullah . Welding on C67 steel grade sheet: Influence of the parameters and post welding heat treatment. AIMS Materials Science, 2023, 10(3): 453-464. doi: 10.3934/matersci.2023025
    [5] Zoubir Chaieb, Ould Mohamed Ouarda, Azzeddine Abderrahmane Raho, Mouhyddine Kadi-Hanifi . Effect of Fe and Si impurities on the precipitation kinetics of the GPB zones in the Al-3wt%Cu-1wt%Mg alloy. AIMS Materials Science, 2016, 3(4): 1443-1455. doi: 10.3934/matersci.2016.4.1443
    [6] Toshio Ogawa, Hiroyuki Dannoshita, Yoshitaka Adachi . Evaluation of tensile properties of ferrite single-phase low-carbon steel with different initial microstructures. AIMS Materials Science, 2019, 6(5): 798-805. doi: 10.3934/matersci.2019.5.798
    [7] M. P. Lavin-Lopez, L. Sanchez-Silva, J. L. Valverde, A. Romero . CVD-graphene growth on different polycrystalline transition metals. AIMS Materials Science, 2017, 4(1): 194-208. doi: 10.3934/matersci.2017.1.194
    [8] Na Ta, Lijun Zhang, Qin Li . Research on the oxidation sequence of Ni-Al-Pt alloy by combining experiments and thermodynamic calculations. AIMS Materials Science, 2024, 11(6): 1083-1095. doi: 10.3934/matersci.2024052
    [9] Sabah Senouci, Mouhyddine Kadi-Hanifi, Azzeddine Abderrahmane Raho . Effect of a plastic deformation on the Guinier-Preston zones growth in Al-5at%Ag alloy. AIMS Materials Science, 2017, 4(1): 16-26. doi: 10.3934/matersci.2017.1.16
    [10] Jun Wang, Hans Conrad . On the dynamic compared to static grain growth rate in 3 mole% yttria-stabilized tetragonal zirconia polycrystals (3 Y-TZP). AIMS Materials Science, 2016, 3(3): 1208-1221. doi: 10.3934/matersci.2016.3.1208
  • The main goal of this work is to clarify and quantify, by means of mathematical analysis, the role of structural viscoelasticity in the biomechanical response of deformable porous media with incompressible constituents to sudden changes in external applied loads. Models of deformable porous media with incompressible constituents are often utilized to describe the behavior of biological tissues, such as cartilages, bones and engineered tissue scaffolds, where viscoelastic properties may change with age, disease or by design. Here, for the first time, we show that the fluid velocity within the medium could increase tremendously, even up to infinity, should the external applied load experience sudden changes in time and the structural viscoelasticity be too small. In particular, we consider a one-dimensional poro-visco-elastic model for which we derive explicit solutions in the cases where the external applied load is characterized by a step pulse or a trapezoidal pulse in time. By means of dimensional analysis, we identify some dimensionless parameters that can aid the design of structural properties and/or experimental conditions as to ensure that the fluid velocity within the medium remains bounded below a certain given threshold, thereby preventing potential tissue damage. The application to confined compression tests for biological tissues is discussed in detail. Interestingly, the loss of viscoelastic tissue properties has been associated with various disease conditions, such as atherosclerosis, Alzheimer’s disease and glaucoma. Thus, the findings of this work may be relevant to many applications in biology and medicine.


    Bainitic press hardening is a key production technology for advanced high strength steels, specifically in the automotive industry [1]. Structural parts like A- and B-pillars, roof rails or bumpers are produced based on this method. Currently, it is of specific interest to the industry to take this important technology to the next level in terms of process optimisation. This technique, using a combination of interrupted cooling and quasi-isothermal holding in the bainitic parameter regime, offers highly desirable strength-failure-strain combinations [2]. The possibility to control the process tool temperature offers a vast amount of possible advances in process technology, and work pieces with locally varying thermoelastic properties for a wide range of applications are one of the central goals in this context. A challenge of central importance to the bainitic press hardening process is the understanding of the bainitic transformation.

    The mechanism of the bainite transformation is the subject of ongoing debate for decades, see e.g., [3,4,5,6]. The main controversy focuses on the following two arguments: (ⅰ) the nature of the bainitic ferrite/retained austenite interface and its capacity to support the displacive mechanism of bainitic ferrite growth and (ⅱ) the role of carbon transport during the bainite transformation, i.e., whether or not bainitic ferrite is initially supersaturated with carbon and then carbon partitioning takes place from the bainitic ferrite. Within this publication, we will focus on the transport-governed regime, i.e., the transformation is limited by carbon redistribution. The main branches of the scientific discussion about the details of the transport-governed regime can be summarised as follows.

    Considerations about the carbon partitioning published in [7,8] suggest that the ferritic component of bainitic microstructures can form with maximal carbon supersaturation, while the escape of carbon from ferrite to austenite after transformation is the essential aspect of the carbon transport. Alternatively, the transport of carbon during diffusional transformation of austenite to ferrite under metastable equilibrium conditions at the interface is also considered (e.g., [5]). An intermediate scenario, where ferrite grows with a partial supersaturation of carbon, with the remaining carbon partitioning into austenite or forming carbides, is described e.g., in [9,10,11]. All of these hypothesises effectively consider the competition between carbon escape and interface migration at high supersaturations, thus the transformation might exhibit a higher velocity than is expected by equilibrium partitioning of carbon. For our investigations, we focus on a regime where the diffusional transport of carbon from the supersaturated austenitic phase at local equilibrium is the rate-limiting mechanism.

    The following parts of this publication are organised in four sections; a basic theoretical treatment of solid-solid transitions in the displacive and the diffuse-displacive case is given in Section 2. We also introduce the basics of the boundary integral method, which is used for the investigations on the terminal subunit growth regime.

    In Section 3.1, we discuss the influence of coherency stresses and capillarity stresses on the growth of banitic subunits of small radii in the initial growth regime of the diffusional-displacive transition.

    In Section 3.2, we present calculations for the specific scenario of bainitic subunit growth frustration, which results form the decreasing temperature during the continuous cooling of the bainitic transformation.

    In Section 4, we discuss the context of other related growth regimes also relevant to elastically influenced solid-solid transformations.

    The work in this article is based on theoretical work and aims at elementary aspects of growth behaviour in bainitic subunits. As representative geometry, we present Figure 1.

    Figure 1.  The geometry of the growing bainitic subunit, approximated as lenticular ferritic inclusion in the retained austenite, with radius R and height h.

    The bainitic transformation is an example for diffusional or mixed diplacive-diffusional solid-solid transitions strongly influenced by elastic effects. An overview of elastical-diffusional transformations can be found e.g., in [12,13,14], in our studies reported in [15,16,17,18,19,20,21,22], we focused on transformations exhibiting coherent interfaces. Though the coupling of displacive and diffusional transitions has been investigated recently also emphasising engineering aspects, see [23], our focus in the course of this publication is on the possible dominance of elastic effects in the initial and terminal regime of subunit growth. When we refer to the terminal growth regime of a bainitic subunit, we consider the freezing-out carbon diffusion upon the continuous decrease of temperature during the progressing bainitic transformation.

    Within this manuscript, we focus on two distinct scenarios. First, the initial growth regime of the subunit is considered. The corresponding model is a sharp interface model, which aims at the prediction of the geometric properties of the growing pattern. Second, we investigate the terminal regime of subunit growth, where we apply the boundary integral method. As the main part of the work is associated to the terminal growth regime, we present here the theoretical basics of the boundary integral method in the context of the physical model we assume for the terminal growth regime. As the boundary integral method is especially useful in the study of steady-state moving fronts, we focus on these steady-state scenarios, here, as systematically integrated in the non-isothermal context of the bainitic transformation. The growth predictions obtained from such steady-state scenarios can be interpreted as valid during the persistence time scale which is basically limited by the variation of the governing parameters of the process.

    To introduce the boundary integral method, we focus on the diffusional aspect of the transformation. We represent the diffusion equation for a field u in a frame of reference which is co-moving along the y-direction symbolically via a linear differential operator L:

    L[u]=2u+2lDyu=0, (1)

    and use the Green's function G of this operator [15] to express the symbolic solution for the diffusion field via Lagrange's identity:

    GL[u]uL[G]=ˉJ(u,G) (2)

    where L is the adjoint operator of L and ˉJ(u,G) is a function of the diffusion field and the Green's function with derivations of order not higher than one in the case of the diffusion equation. Apparently, when we integrate the Lagrange identity and have no volume inhomogeneities, the integral contributions reduces to boundary terms. In our case, the inhomogeneity is the mass conservation of carbon at the moving interface, as we discuss in the subsequent description of the physical model. Finally, the whole description reduces to a fully interfacial problem as the obtained symbolic solution expresses the value of the diffusion field on each point at the interface when the local equilibrium condition is described, see also the following description of the physical model.

    For the terminal growth regime, the basis of our considerations is a coupled model of diffusional carbon transport and static elastic misfit stresses within the two adjoining phases, i.e., the parent austenite and the growing ferrite phase, together with the physically necessary boundary conditions at the progressing phase boundary. The diffusional transport is governed by a diffusion equation in the bulk, complemented by the mass conservation condition and the Gibbs-Thompson condition of local phase equilibrium. In total, we obtain the following three equations for the dimensionless carbon concentration field u=(cc0)/Δc, with c being the non-equilibrium carbon concentration field, c0 is a reference concentration and Δc=ceqγceqα denoting the difference in the equilibrium carbon concentrations of the γ-austenite and the α-ferrite phase at the respectively imposed constant temperature

    D2u=u/t,υn=DnΔu|int,u|int=Δdκ+TeqδFel/LmΔc. (3)

    By D we denote the carbon diffusion constant, Teq is the transition temperature, L is the effective latent heat [24,25]. The capillary length d is defined here as d=γTeq/LmΔc, where γ is the surface energy, κ is the interface curvature and the dimensionless undercooling is given by Δ=(ceqc0)/(cγcα). The slope of the coexistence line is denoted by m=dT/dc. In accordance with our previous studies on displacive transformations, we choose the class of invariant plain strain transitions here, as it allows for extended coherent interfaces. When we indicate the mother phase via α and the onsetting phase via β, the elastic contribution to the local equilibrium equation is defined as

    δFel=σ(0)ikϵ(β)ikE[(ϵ0ττ)2+(ϵ0ss)2+2ν(ϵ0ττ)(ϵ0ss)+2(1ν)(ϵ0sτ)2]2(1ν2). (4)

    Here the indices ik refer to the components of the strain and eigenstrain tensors, with ss, sτ and ττ denoting the two tangential components in natural coordinates. This is obtained by introduction of the appropriate free energy and assuming phase and mechanical equilibrium at the coherent interface, see [15] for details. We require that the eigenstrain tensor ϵ(0)ik allows to satisfy coherent interfaces, the eigenstrain stress tensor σ(0)ik is the stress tensor corresponding to Hooke's law for the eigenstrains ϵ(0)ik. Therefore, the first contribution on the r.h.s of Eq 4 involves the eigenstresses and the strains in the austenitic phase. Following this brief introduction of our theoretical basis, we continue with the discussion of the initial growth regime.

    The physical picture emerging from these equations is clear. The propagation of the interface releases carbon proportional to the normal component of the growth velocity υn, while the local transition point itself depends on the interfacial stresses arising from the coherent lattice transformation and the curvature. The transport of the carbon from the interface limits the propagation velocity. We note that this model represents the coupling of the elastostatic and the diffusional problem via the local equilibrium condition for the carbon concentration at the transformation interface alone.

    The transformation of austenite to bainite involves displacive deformations of the matrix, diffusion controlled growth of bainitic sheaves and the precipitation of carbides inside the ferrite (lower bainite) or inside the austenite (upper bainite). A successful prediction of the transformation kinetics requires consideration of the complex interaction of these processes: elastic/plastic deformation, carbon partitioning and carbide precipitation. The local driving force for the transformation depends sensitively on the carbon concentration, the local misfit-stresses and the temperature. During the transformation, the carbon in the supersaturated bainitic-ferrite can either be released in the surrounding austenite phase or be captured by the embedded carbides, which are competitive for the solute atoms.

    When we focus on the very initial regime of the solid solid transformation, involving a nucleus of nanometer size, the size dependence of the elastic parameters in this stage can lead to substantial elastic softening of the growing phase. Furthermore, on this scale, the bulk crystal symmetry structure can be suppressed due to the dominance of interfacial effects on the overall energy of the nucleus. Therefore, we approximate the growing inclusion by concerning only the hydrostatic, isotropic elastic response to stresses from the surrounding bulk phase. The assumption of isotropic elasticity exhibits the lowest expected error on the real elastic behaviour of the nucleus prior any information on the early nm growth regime. Thus, we consider this initial stage in analogy to the argumentation given in [12], which provides an approximative and analytically accessible description of the growth problem.

    We assume that the new phase, which contains N particles and occupies a volume W, is homogeneous with chemical potential μ and a pressure P (we assume that the external pressure is zero). Because of the conservation of mass, we have the following expression for the volume of the melt:

    W=Nvα=Nvγ+unds (5)

    where vα and vγ are the atomic volumes of the ferritic and austenitic phases, and the last term describes the change in the volume of the nucleus upon deformation (un is the normal component of the displacement vector at the interface). The subunit is assumed to be a very oblate lentil of radius R and height h R. In a first approximation, we can ignore the height of the lentil in solving the elastic problem. In this case the pressure P, exerted on the nucleus by the retained austenite, is given on a circular plane cut within the radius R. This problem is equivalent to a crack problem and its solution provides the normal component of the displacement vector at the interface, see [26]. In this scenario, the equilibrium Gibbs-Thomson condition at the interface reads as [27,28]

    μ(P)=vα[f0γ+P+ˉσ2rr+ˉσ2ϕϕ2νˉσrrˉσϕϕ2E+γκ] (6)

    Here f0γ is the free energy density of the undeformed (initial) crystal and we have introduced cylindrical coordinates, ˉσik is a normalised stress tensor, ν and E are the Poisson and Young coefficients, γ is the interfacial energy and κ is the curvature of the interface (considered positive for a convex solid). To focus on the role of the elastic effects of the transformation, Eq 6 approximates the iron chemical potential by reducing contributions depending on the carbon content to an effective lattice expansion which is taken into account via the atomic volumes vα, vβ.

    Following the argumentation in [29] and [12], the radius Rc of the critical nucleus is related to the pressure P by the Griffith formula, leading to a relation Δ(R)=Δ(11Rc/R)=uInt for the interfacial equilibrium condition. This relation together with mass conservation provides the characteristic length scale ρ and the height h of the inclusion

    h=ρR (7)
    ρ=128γ(1ν2)v2γ9πE(vαvγ)2 (8)

    Apparently, these expressions remain valid only when the elastic softening is pronounced and the assumption that neglecting the nonydrostatic aspect of the elastic behaviour in the inclusions is acceptable. However, within the approximative validity of these assumptions, we can consider the onset to be elastically dominated. The crucial point is the absence of any thermal quantities in ρ. This indicates that the pattern is dominated by elastic effects in this early growth regime. After the shape undergoes a Mullins-Sekerka instability, which takes place as derived in [30], it grows via the coupled diffusional-displacive mechanism, qualitatively comparable to the dendritic growth regime as we described it in [19]. In the next section, we introduce the regime of frustrated growth of bainitic subunits, which is represented by a limiting case of the elasto-chemical diffusion-limited growth, i.e., diffusion-limited growth mediated by an elastic selection.

    Within this section, we focus on the elasto-chemical diffusion-limited growth of a subunit in the terminal regime of carbon transport frustration, i.e., where an effective reduction of the thermodynamic driving force due to the size of the grown subunit can be assumed. The frustration of the carbon transport within the retained austenite is related to the continuous decrease of the overall temperature during the progression of the bainitic transformation. Upon a linear temperature decrease the carbon diffusion coefficient undergoes an exponential decay, as it obeys an Arrhenius-type temperature dependence. The combination of a reduced driving force and an exponentially decreasing diffusion coefficient suggest to study the zero velocity or zero diffusion limit.

    To take full advantage of our assumptions, the boundary integral method, as introduced in the methods section, see Eq 3 and in [19], is most suitable. The resulting formulation of the problem combines the governing equations of the elastostatic and transport problem into a closed representation for the unknown boundary, y(x), and the tip radius of curvature R:

    ΔdκR+TeqcpδFelL2=pπ+dxexp[p(y(x)y(x))]K0(pη(x,x)). (9)

    To simplify the representation, we introduce a splitting of the elastic contribution to the local equilibrium into a constant, hysteresis contribution Δel and the remaining, spatially varying contribution, for details we refer to [19]. Important here is to note that the Peclet number of the problem depends on the net driving force which takes the elastic hysteresis into account, i.e., ˜Δ=ΔΔel=pπexp(p)erfc(p). This includes that a strong elastic hysteresis leads to a substantially reduced driving force for the transformation, and we use approximations for the dependence of the effective driving force on the Peclet number and the Green's function of the problem as described in [19] to obtain the following representation:

    σκ+ΔelpΦ[σϵiku(α)ikEϵ2B(x)2(1υ2)]=1π+dxlog[(xx)2+(y(x)y(x))2(xx)2+(x22+x22)2]1/2. (10)

    Here we define Φ=2(1ν2)/(E(ϵ2yy+ϵ2zz+2νϵyyϵzz)), function B(x) includes the local contributions to the elastic free energy difference. The stability parameter of classical dendritic growth is σ=d0/pR. We note that the velocity and the radius of the asymptotically matching Ivantsov parabola are defined as

    υdD=2π2σ˜Δ4, (11)
    Rd=πσ˜Δ2. (12)

    As mentioned, the specific scenario we apply the boundary integral formulation for is a scenario of very slow diffusion, so that the interface propagates slowly. Then, the diffusion is sufficiently low to have the elastic effects dominate the transition. In the case of very small Δel/p, one shall exploit the possibilities of analytical methods, e.g., apply the complex matching method which yielded the solution of the classical dendritic problem, but for the limit of p/Δel 1, which represents the scenario of transport and frustrated transformation in the subunits, analytical approaches in the aforementioned spirit are not at hand, suggesting specific numerical consideration here.

    Therefore, we approximate this frustrated transformation regime via the limit of zero velocity and diffusion. Formally, we can reach the regime of a possibly negligible release of latent heat at the interface when υ0. The Peclet number is given for a fixed ˜Δ1 as

    υR2D=p˜Δ2, (13)

    which shows that for a given radius of the subunit, i.e., if R(˜Δ)1, the limiting case of a decreasing diffusion coefficient D0 is equivalent to υ0. As the exponential dependence of the carbon diffusion coefficient on the temperature provides a pragmatically reasonable approximation to the strong decay of carbon diffusion, we will consider if a scenario with R(˜Δ)1 can be predicted.

    We introduce a representation of the problem which suits this approach, namely

    dΔelRκ+Φ[σϵiku(α)ikEϵ2B(x,x)2(1ν2)]=1πpΔel+dxlog[(xx)2+(y(x)y(x))2(xx)2+(x22+x22)2]1/2. (14)
    Φ=2(1ν2)Eϵ2(ˉϵ2yy+ˉϵ2zzˉϵyyˉϵ2zz)

    and the corresponding solvability parameter reads as

    dΔelR=σel=1R2γ(1ν2)Eϵ2(ˉϵ2yy+ˉϵ2zzˉϵyyˉϵ2zz) (15)

    We note that we have introduced a solvability parameter which describes the growth regime without involving any diffusional transport parameters. This reflects the possible elastic selection on the formal level. We implemented a symmetric diffusion model instead of the one-sided model which is typically used for highly asymmetric carbon transport. Therefore, we note that the dependence of the stability parameter on the control parameter in the symmetric and one-sided model is qualitatively identical. The symmetric model implementation is numerically easier to handle in the regime of p/Δel0. The resulting calculations for pure shear can cover the regime up to p/Δel=0.05, but even at this already small value, an extrapolation to p/Δel0 remains inconclusive. As the plot in Figure 2 shows, we can at least assume that the hidden behaviour of the selection does not obey a linear law. Exemplarily, we discuss the result of a potentially present finite σel and the result for the case σel(p/Δel)1/2.

    Figure 2.  The close up of the regime up to p/Δel=0.5. The solid line was obtained in the asymptotic regime and matches the results obtained for Δel/p vs σ before, represented by the dashed line. The results for the asymptotic regime Δel/p remain inconclusive concerning the presence of a selection for zero velocity.

    A finite value of σel at zero p/Δel would lead to an elastic selection as

    R=dΔelσel1 (16)
    υ=2ΔelσelπDd˜Δ2˜Δ2, (17)

    where we note that the net driving force relative to the elastic hysteresis in our model is assumed to be reduced by the frustration of the transformation. A possible square-root dependence of the solvability parameter on p/Δel would lead to scaling laws for the transformation velocity and the radius as

    R1˜Δ, (18)
    υ~Δ3. (19)

    Despite of the inconclusive numerical results for the zero velocity limit concerning a finite σel, the results show that an elastically dominated terminal growth regime of the subunits is possible.

    We presented our recent work on the possibilities of elastically dominated regimes in diffusion-limited solid state transformations. Referencing the bainitic transformation, specifically the scenario of sheaf growth via subunit growth into the retained austenite, the very initial and the terminal growth regime are considered. For the early growth regime, we find that under strong assumptions for the elastic softening and effective isotropy of the onsetting phase, the problem can be considered analogous to a crack formation problem. This is in line with previously published investigations on melt inclusion growth published in [29]. When the Mullins-Sekerka instability sets the transition from this early, elastically governed growth regime to the coupled diffusive-displacive growth regime, described in [15,16,18,19], the transport of carbon in the retained austenite is assumed to be the limiting factor of subunit growth. Taking into account the decrease of the carbon diffusion coefficient with decreasing temperature suggests to consider the effect of a reduced diffusion coefficient. In combination with elastically frustrated growth, i.e., reduction of the effective driving force, the reduced diffusion coefficient suggests to study a zero velocity limit. The presented numerical results show the possibility of an elastically dominated terminal growth regime of subunits.

    N. T. would like to acknowledge financial support from Innovation Foundation For Postgraduate and Fundamental Research Funds of Central South University (Grant No. 1053320171074), Changsha, China. K. W. would like to acknowledge financial support from Chinese Scholarship Council (Grant No. 201706370221), China. X. Y., M. F. and C. H. acknowledge the support of the Deutsche Forschungsgemeinschaft under Project No. SPP 1713.

    All authors declare no conflict of interest in this paper.

    [1] [ P. Augat,S. Schorlemmer, The role of cortical bone and its microstructure in bone strength, Age and Ageing, 35 (2006): ii27-ii31.
    [2] [ H. T. Banks,K. Bekele-Maxwell,L. Bociu,M. Noorman,G. Guidoboni, Sensitivity analysis in poro-elastic and poro-visco-elastic models with respect to boundary data, Quart. Appl. Math. In press, 75 (2017): 697-735.
    [3] [ H. Barucq,M. Madaune-Tort,P. Saint-Macary, On nonlinear Biot's consolidation models, Nonlinear Anal Theory Methods Appl., 63 (2005): e985-e995, Invited Talks from the Fourth World Congress of Nonlinear Analysts (WCNA 2004).
    [4] [ H. Barucq,M. Madaune-Tort,P. Saint-Macary, Some existence-uniqueness results for a class of one-dimensional nonlinear Biot models, Nonlinear Anal Theory Methods Appl., 61 (2005): 591-612.
    [5] [ M. A. Biot, General theory of three-dimensional consolidation, J. Appl. Phys., 12 (1941): 155-164.
    [6] [ L. Bociu,G. Guidoboni,R. Sacco,J. Webster, Analysis of nonlinear poro-elastic and poro-viscoelastic models, Arch. Rational Mech. Anal., 222 (2016): 1445-1519.
    [7] [ S. Canic,J. Tambaca,G. Guidoboni,A. Mikelic,C. J. Hartley,D. Rosenstrauch, Modeling viscoelastic behavior of arterial walls and their interaction with pulsatile blood flow, SIAM J. Appl. Math., 67 (2006): 164-193.
    [8] [ Y. Cao,S. Chen,A. J. Meir, Analysis and numerical approximations of equations of nonlinear poroelasticity, Discrete Continuous Dyn Syst Ser B, 18 (2013): 1253-1273.
    [9] [ P. Causin,G. Guidoboni,A. Harris,D. Prada,R. Sacco,S. Terragni, A poroelastic model for the perfusion of the lamina cribrosa in the optic nerve head, Math. Biosci., 257 (2014): 33-41.
    [10] [ P. Causin,R. Sacco,M. Verri, A multiscale approach in the computational modeling of the biophysical environment in artificial cartilage tissue regeneration, Biomech. Model. Mechanobiol, 12 (2013): 763-780.
    [11] [ D. Chapelle,J. Sainte-Marie,J.-F. Gerbeau,I. Vignon-Clementel, A poroelastic model valid in large strains with applications to perfusion in cardiac modeling, Comput. Mech., 46 (2010): 91-101.
    [12] [ O. Coussy, Poromechanics, John Wiley & Sons Ltd, 2004.
    [13] [ S. C. Cowin, Bone poroelasticity, J. Biomech., 32 (1999): 217-238.
    [14] [ R. de Boer, Theory of Porous Media. Highlights in the Historical Development and Current State, Springer, Berlin/New York, 2000.
    [15] [ E. Detournay,A. H.-D. Cheng, Poroelastic response of a borehole in a non-hydrostatic stress field, Int J Rock Mech Min Sci Geomech Abstr, 25 (1988): 171-182.
    [16] [ E. Detournay,A. H.-D. Cheng, Fundamentals of poroelasticity, Comprehensive rock engineering, 2 (1993): 113-171.
    [17] [ J. C. Downs,J. K. Suh,K. A. Thomas,A. J. Bellezza,R. T. Hart,C. F. Burgoyne, Viscoelastic material properties of the peripapillary sclera in normal and early-glaucoma monkey eyes, Invest. Ophthalmol. Vis. Sci., 46 (2005): 540-546.
    [18] [ J. W. Freeman,M. D. Woods,D. A. Cromer,L. D. Wright,C. T. Laurencin, Tissue engineering of the anterior cruciate ligament: The viscoelastic behavior and cell viability of a novel braid-twist scaffold, J Biomater Sci Polym Ed, 20 (2009): 1709-1728.
    [19] [ I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, Elsevier/Academic Press, Amsterdam, 2015.
    [20] [ A. P. Guérin,B. Pannier,S. J. Marchais,G. M. London, Arterial structure and function in end-stage renal disease, Curr. Hypertens. Rep., 10 (2008): 107-111.
    [21] [ W. M. Lai,J. S. Hou,V. C. Mow, A triphasic theory for the swelling and deformation behaviors of articular cartilage, ASME J. Biomech. Eng., 113 (1991): 245-258.
    [22] [ V. C. Mow,S. C. Kuei,W. M. Lai,C. G. Armstrong, Biphasic creep and stress relaxation of articular cartilage in compression: Theory and experiments, ASME J. Biomech. Eng., 102 (1980): 73-84.
    [23] [ H. T. Nia,L. Han,Y. Li,C. Ortiz,A. Grodzinsky, Poroelasticity of cartilage at the nanoscale, Biophys. J., 101 (2011): 2304-2313.
    [24] [ M. S. Osidak,E. O. Osidak,M. A. Akhmanova,S. P. Domogatsky,A. S. Domogatskaya, Fibrillar, fibril-associated and basement membrane collagens of the arterial wall: Architecture, elasticity and remodeling under stress, Curr. Pharm. Des., 21 (2015): 1124-1133.
    [25] [ S. Owczarek, A Galerkin method for Biot consolidation model, Math. Mech. Solids, 15 (2010): 42-56.
    [26] [ N. Özkaya, M. Nordin, D. Goldsheyder and D. Leger, Fundamentals of Biomechanics. Equilibrium, Motion, and Deformation, Springer, New York, 1999.
    [27] [ P. J. Phillips,M. F. Wheeler, A coupling of mixed and continuous Galerkin finite-element methods for poroelasticity Ⅰ: The continuous in time case, Comput Geosci, 11 (2007): 131-144.
    [28] [ P. J. Phillips,M. F. Wheeler, A coupling of mixed and continuous Galerkin finite-element methods for poroelasticity Ⅱ: The continuous in time case, Comput Geosci, 11 (2007): 145-158.
    [29] [ P. J. Phillips,M. F. Wheeler, A coupling of mixed and continuous Galerkin finite-element methods for poroelasticity, Comput Geosci, 12 (2008): 417-435.
    [30] [ L. Recha-Sancho, F. T. Moutos, J. Abell, F. Guilak and C. E. Semino, Dedifferentiated human articular chondrocytes redifferentiate to a cartilage-like tissue phenotype in a poly(ε-caprolactone)/self-assembling peptide composite scaffold, Materials, 9 (2016), p472.
    [31] [ R. Sacco,P. Causin,P. Zunino,M. T. Raimondi, A multiphysics/multiscale 2D numerical simulation of scaffold-based cartilage regeneration under interstitial perfusion in a bioreactor, Biomech. Model. Mechanobiol., 10 (2011): 577-589.
    [32] [ I. Sack,B. Beierbach,J. Wuerfel,D. Klatt,U. Hamhaber,S. Papazoglou,P. Martus,J. Braun, The impact of aging and gender on brain viscoelasticity, NeuroImage, 46 (2009): 652-657.
    [33] [ A. P. S. Selvadurai, On the mechanics of damage-susceptible poroelastic media, Key Engineering Materials, 251/252 (2003): 363-374.
    [34] [ A. Settari and D. A. Walters, Advances in Coupled Geomechanical and Reservoir Modeling with Applications to Reservoir Compaction, Technical report, SPE Reservoir Simulation Symposium, 1999.
    [35] [ R. E. Showalter, Diffusion in poro-elastic media, J. Math. Anal. Appl., 251 (2000): 310-340.
    [36] [ R. E. Showalter,N. Su, Partially saturated flow in a poroelastic medium, Discrete Continuous Dyn Syst Ser B, 1 (2001): 403-420.
    [37] [ M. A. Soltz,G. A. Ateshian, Experimental verification and theoretical prediction of cartilage interstitial fluid pressurization at an impermeable contact interface in confined compression, J. Biomech., 31 (1998): 927-934.
    [38] [ K. Terzaghi, Erdbaumechanik auf Bodenphysikalischer Grundlage, Deuticke, Wien, 1925.
    [39] [ A. A. Tofangchi Mahyari, Computational Modelling of Fracture and Damage in Poroelastic Media, PhD thesis, Department of Civil Engineering and Applied Mechanics, McGill University, 1997.
    [40] [ A. Zenisek, The existence and uniqueness theorem in Biot's consolidation theory, Apl Mat, 29 (1984): 194-211.
  • This article has been cited by:

    1. Mingxuan Lin, Carina Zimmermann, Kai Wang, Martin Hunkel, Ulrich Prahl, Robert Spatschek, Modeling Bainitic Transformations during Press Hardening, 2021, 14, 1996-1944, 654, 10.3390/ma14030654
  • Reader Comments
  • © 2018 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4620) PDF downloads(632) Cited by(18)

Figures and Tables

Figures(15)  /  Tables(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog