Abstract
Effective properties of non-aging linear viscoelastic and hierarchical composites are investigated via a three-scale asymptotic homogenization method. In this approach, we consider the assumption of a generalized periodicity in the different structural levels and their characterization through the so-called stratified functions. The expressions for the associated local and homogenized problems, and the effective coefficients are derived at each level of organization by using the correspondence principle and the Laplace-Carson transform. Considering isotropic components and a perfect contact at the interfaces between the constituents, analytical solutions, in the Laplace-Carson space, are found for the local problems and the effective coefficients are computed. An interconversion procedure between the effective relaxation modulus and the effective creep compliance is carried out for obtaining information about both viscoelastic properties. The numerical inversion to the original temporal space is also performed. Finally, we exploit the potential of the approach and study the overall properties of a hierarchical viscoelastic composite structure representing the dermis.
1. Introduction
In the scientific literature there exist several works focusing on the development of micromechanical techniques to predict the macroscopic properties of composite materials. An excellent review on these methods can be found in Sevostianov and Giraud [Citation1]; Davit et al. [Citation2]. In particular, multiscale asymptotic homogenization methods take advantage of the information available at the smaller scales of a given heterogeneous medium to predict the effective properties at its larger scales (see e.g. Lukkassen and Milton [Citation3]; Penta and Gerisch [Citation4]; Ramírez-Torres et al. [Citation5]; Ramírez-Torres et al. [Citation6]; Yang et al. [Citation7]; Dong et al. [Citation8]). This homogenization procedure requires the solution of cell problems with input data corresponding to the homogenized material properties resulting in previous steps. In addition, the use of more general periodic or stratification functions (see e.g. Tsalis, Chatzigeorgiou, and Charalambakis [Citation9]; Tsalis et al. [Citation10]) permits to describe the different length scales of the composite materials. They can be applied to homogenization problems related to hierarchical real-world composites as human tissue (see e.g. Penta et al. [Citation11] for bone and tendons) that generally form geometries which cannot be produced by the repetition of one elementary volume.
In the present work, the three-scale Asymptotic Homogenization Method (AHM) introduced in Ramírez-Torres et al. [Citation5]; Ramírez-Torres et al. [Citation6]; Ramírez-Torres et al. [Citation12] is used to model a non-aging linear viscoelastic composite material with generalized periodicity and two hierarchical levels. This three-scale asymptotic approach has been used for modeling hierarchical laminated and fiber-reinforced elastic composites in Ramírez-Torres et al. [Citation5] and Ramírez-Torres et al. [Citation6]; Ramírez-Torres et al. [Citation12], respectively. Also, there are several other works dealing with the modeling of multiscale hierarchical heterogeneous media. For example, from the theoretical point of view we find the works Allaire and Briane [Citation13]; Bensoussan, Papanicolau, and Lions [Citation14]; Trucu, Chaplain, and Marciniak-Czochra [Citation15]. Regarding these results, a rigorous study about the multiscale convergence is developed in Bensoussan, Papanicolau, and Lions [Citation14]; here the authors consider that and are the local variables. Subsequently, and applied to the heat equation for composites, the method of reiterate homogenization is well founded in Allaire and Briane [Citation13]. An additional overview of reiterated homogenization is presented in Trucu, Chaplain, and Marciniak-Czochra [Citation15] via a three-scale convergence approach where the asymptotic parameters independently approach zero. More recently, in Dong et al. [Citation8]; Yang et al. [Citation7, Citation16], the authors study the properties of thermo-mechanical, non-aging and aging viscoelastic composites with multiple spatial scales by using a three-scale asymptotic expansion and a periodic layout of the heterogeneities in the structures. They also provide a finite element algorithm based on inverse Laplace transform and the three-scale asymptotic homogenization to obtain the numerical results.
The article presents some novelties with respect to previous works of ours [Citation5, Citation17] and, to the best of our knowledge, with others in the literature [Citation7, Citation16, Citation18–20]). Specifically, we generalize the results obtained in [Citation5] for linear elasticity by extending them to a non-aging, linear viscoelasticity framework; we deal with the stratified functions in the homogenization procedure; an analytical solution for the local problems associated with each hierarchical level is obtained; the expressions for the effective coefficients for hierarchical laminated composites with generalized periodicity, isotropic components and perfect contact at the interfaces are provided; and the methodology is used to model the overall behavior of the dermis in the skin. A better understanding of this tissue has a real impact on biomedical applications and in the development of modern technology such as flexible instance electronics, soft robotics and prosthetics (see Ref. [Citation21]). Many researches are focused on the study of the viscoelastic, non-linear hyperelastic and anisotropic behavior of the structural components of the skin (see Ref. [Citation22, Citation23]) and the relation with its hierarchical structure (see Ref. [Citation21, Citation24]). An inherent characteristic of the mechanical properties of the skin lies in the fact that its properties vary according to the skin type, age, gender, location, and other environmental factors Panchal et al. [Citation25]. In particular, skin has three main layers: the epidermis, the dermis and the hypodermis, in addition to the stratum corneum that acts as the outermost barrier of human skin Muha et al. [Citation26]. Our research is focused on the dermis which represents the most important mechanical and thermal unit of the skin, as it occupies around 90% of the skin’s thickness (see Ref. [Citation27]). In the present work, the dermis is assumed to behave as a non-aging linear viscoelastic and hierarchical composite material, and thus, the investigation of its effective properties is based on the correspondence principle and the Laplace transform (see Ref. [Citation28]). The procedure, see for instance To et al. [Citation29] and Yang et al. [Citation16], consists in the change of the convolution constitutive law by a fictitious elastic one in the Laplace domain. Then, the inversion of the Laplace transform is considered to derive the effective behavior in the time domain.
The article is organized as follows. The mathematical framework of the problem and the homogenization steps are illustrated in Sections 2 and 3, respectively. In Section 4, we present the analytical solution of the local problems and the calculation of the effective coefficients. In Section 5, we highlight the potential of the current approach in the investigation of the effective relaxation and creep compliance properties of the dermis.
2. The linear viscoelastic problem for hierarchical structures
2.1. Separation of scales
Let us consider the dimensionless, scaling parameters and defined as follows, (see Ref. [Citation12]) (1) (1)
In addition, the well-separated characteristic length scales and L introduce three dimensionless spatial variables in the reference configuration, (2) (2) where x is said to be the physical spatial variable, whereas and represent the macroscopic, mesoscopic and the microscopic dimensionless spatial variables, respectively.
By using Eq. (Equation1(1) (1) ), and can be related through the expressions, (3) (3) (4) (4)
Now, by defining a field over the region of interest, the separation of scales allows to rephrase the space dependence of as and taking into account Eqs. Equation(1)–(4), the spatial derivative of can be computed as follows, (5) (5)
In the present work, we consider the generalized periodicity (see Ref. [Citation9, Citation10]) as an additional condition in our model, and thus, we rewrite the three dimensionless spatial variables in Eq. (Equation2(2) (2) ) as follows (6) (6) where with and the property represent the so-called stratified functions. It is worthy to note that the parametric equations describe the structural organization of the constituents at the and hierarchical levels, respectively.
Therefore, we have that (7) (7) (8) (8)
In addition, the spatial derivative of is given by the formula (9) (9)
By following this approach, all equations should be written in dimensionless form. In the literature, the switch to auxiliary variables and is often omitted. However, as shown for example in Auriault, Boutin, and Geindreau [Citation30], both paths are equivalent. Therefore, in the article, the analysis is carried out directly in a system of physical variables x, y and z, where (10) (10)
Similar ideas are developed in Ramírez-Torres et al. [Citation31] to the case of two scales.
2.2. The physical model
Let us denote by a linear multiscale viscoelastic composite material (see ) characterized by a generalized periodicity at the different levels of organization.
First stage: -structural level:
The domain Ω is considered to be a two-phase quasi-periodic composite such that and (see ). It is assumed that where N1 denotes the number of inclusions, and we denote with the interface between and In this description, Y represents the unitary periodic cell and the its constituents fulfill the constraint and The interface between Y1 and Y2 is denoted with ΓY.
Second stage: -structural level
At this structural level we consider that, at the same time, each phase () is a two-phase quasi-periodic composite material, where and (see ). Besides, and represents the interface between and Besides, Z is the corresponding unitary periodic cell, in which we enforce that the constituents satisfy the following relation and The interface between Z1 and Z2 is indicated with ΓZ.
At this point, we extend the concept of generalized or quasi-periodicity, given in Tsalis, Chatzigeorgiou, and Charalambakis [Citation9]; Tsalis et al. [Citation10], to the case of hierarchical composites with three scales. The definition is presented as follows,
Definition 2.1.
(Generalized periodicity) A composite material with two structural levels of organization is said to have generalized periodicity (or quasi-periodicity) if there exists a global spatial coordinate system x and functions with (see Eq. (Equation10(10) (10) )) such that the operator which relates the second-order stress tensor and the second-order strain tensor by the relation is regular in x, and periodic in y and z.
This concept will permit us to study different structural effects in the composite.
Remark 1.
From now on, we consider the following notation where is assumed to be periodic in y and z.
2.3. Formulation of the problem
Considering that the constitutive response of all the constituents of the composite body is linear viscoelastic, and ignoring inertial terms, the problem in Ω reads, (11a) (11a) (11b) (11b) (11c) (11c) (11d) (11d) where, by taking into account Definition 1, is given by the following scale-dependent constitutive law (see Ref. [Citation32]), (12) (12)
In Eq. (Equation12(12) (12) ), denotes the second-order strain tensor determined by the formula (13) (13) and u is the displacement field. Moreover, f denotes the action of external volume forces that satisfies and are the prescribed displacement and traction on the boundary with and n is the outward unit vector normal to the surface The relaxation modulus is denoted by with the following symmetry properties Furthermore, we assume that and that it is positively defined, i.e. the relation is verified for all symmetric real valued second-order tensors and some positive constant β.
We further impose continuity conditions for displacements and tractions on both and i.e. (14) (14) where the outward unit vectors to the surfaces and are represented by and respectively. The operator denotes the “jump” of across the interface between the constituents.
2.4. Re-formulation of the problem in Laplace-Carson space
The scale-dependent constitutive law Eq. (Equation12(12) (12) ) corresponds to the special form of a non-aging linear viscoelastic materials where we emphasize the specific functional dependency on (see Ref. [Citation33]). Therefore, the viscoelastic problem given by Eqs. Equation(11a–14) can be transformed into an elastic one using the Laplace-Carson transform. It is defined for all real numbers as (see Appendix B in Christensen [Citation32]) where the variable p denotes the Laplace-Carson space.
The aforementioned transformation is known as the correspondence principle and permits us to rewrite the system Eqs. Equation(11a–11d) in the Laplace-Carson space, (15a) (15a) (15b) (15b) (15c) (15c) (15d) (15d)
Furthermore, the interface contact conditions become, (16) (16)
Remark 2.
In what follows the functions depending on the variable p (e.g. ) are defined in the Laplace-Carson space and the homogenization process is developed in Laplace-Carson space.
3. The three-scale asymptotic homogenization approach
In this section we recall the three-scale asymptotic homogenization approach introduced in Ramírez-Torres et al. [Citation5] and we adapt it to the investigation of the overall behavior of hierarchical, linear viscoelastic composites with generalized periodicity of its internal structure.
We first notice that, since each field and material property are assumed to be regular in x and periodic in y and z, the separation of scales together with Eq. (Equation10(10) (10) ) imply that, (17) (17) and using a similar idea, the EquationEq. (13)(13) (13) becomes, (18) (18) where (19a) (19a) (19b) (19b)
3.1. Three-scale homogenization procedure
Following Ramírez-Torres et al. [Citation5], the solution to the problem Eqs. Equation(15a–16) is given as follows, (20) (20) where is defined as (21) (21)
Here and subsequently (unless necessary), the variable dependence is dropped out for convenience.
Replacing expansion Eq. (Equation20(20) (20) ) into Eq. (Equation15a(15a) (15a) ) and using the relations Eqs. (Equation17(17) (17) ) and (Equation18(18) (18) ), it is obtained (22) (22) and analogously, the interface conditions Eq. (Equation16(16) (16) ) read, (23a) (23a) (23b) (23b) (23c) (23c)
Before proceeding, the following cell average operators are introduced, (24) (24) where and represents the volume of the periodic cell Z and Y, respectively.
3.1.1. First level of the hierarchical structure
Grouping in powers of ε2 Eqs. Equation(22–23b), and after some algebraic manipulations, the sequence of problems depicted below can be solved in recursive form. A summary for each problem is proposed.
Problem for (25a) (25a) (25b) (25b) (25c) (25c)
Since the right hand side of Eq. (Equation25a(25a) (25a) ) is zero, is a solution of Eq. (Equation25a(25a) (25a) ) if and only if it is constant in relation to the variable z, i.e. (26) (26)
The validity of Eq. Equation(26)(26) (26) is analyzed in Tsalis, Chatzigeorgiou, and Charalambakis [Citation9]. Thus, taking into account Eq. (Equation26(26) (26) ) and Eq. (Equation21(21) (21) ) is obtained (27a) (27a) (27b) (27b)
Problem for
Using the fact that the problem reads (28a) (28a) (28b) (28b) (28c) (28c)
The existence and uniqueness of a solution for problem Eqs. Equation(28a–28c) is proved in Tsalis, Chatzigeorgiou, and Charalambakis [Citation9] (Proposition 1, page 18), which is based on Lax-Milgram’s theorem and Korn’s inequality (see Ref. [Citation34]). The divergence theorem and the periodicity of in the variable z ensure the fulfillment of the solvability conditions in Eq. Equation(28a)(28a) (28a) .
We can use the method of separation of variables to obtain a general solution for Eq. Equation(28a–28c) as follows (29) (29) where the third order tensor is periodic in y and z and (30) (30)
Then, substituting Eq. (Equation29(29) (29) ) into Eq. Equation(28a–28c), and after some simplifications, the following auxiliary problem, which we call the -cell problem, is obtained (31a) (31a) (31b) (31b) (31c) (31c)
The expression in Eq. Equation(31a)(31a) (31a) is defined taking into account Eq. Equation(19b)(19b) (19b) as follows (32) (32)
Problem for
From Eq. (Equation22(22) (22) ), the expressions are grouped in powers of Using Eq. (Equation29(29) (29) ), the problem is given by (33) (33)
Averaging Eq. Equation(33)(33) (33) and taking into account the divergence theorem and the periodicity of the involved functions in the variable z (see Ref. [Citation9, Citation34]), we obtain that (34) (34) where (35) (35)
is the effective coefficient at the -structural level.
3.1.2. Second level of the hierarchical structure
In the previous Section, we have carried out an analysis for finding the homogenized properties at the -structural level. Now, in order to find the overall behavior of the hierarchical composite, we follow a similar procedure where the results found so far become the new input values for the new problems.
Using relations Eqs. (Equation21(21) (21) ) and (Equation30(30) (30) ) into Eq. Equation(34)(34) (34) , it can be verified that (36) (36)
Considering the representation Eq. (Equation21(21) (21) ) and relation Eq. Equation(23a)(23a) (23a) , the following equation is obtained (37) (37)
Substituting Eq. (Equation21(21) (21) ) and Eqs. Equation(29)–(30) into Eq. (Equation23c(23c) (23c) ) and taking into account that is z-constant, (38) (38)
Working with Eqs. Equation(36–38) and grouping in powers of the following sequence of problems are analyzed.
Problem for (39) (39) where the following result is reached (see Ref. [Citation9]) (40) (40)
From Eqs. Equation(37–38) and applying the cell average operator over (41) (41)
Problem for
At this point, using Eqs. Equation(36–38), the fact that and applying the cell average operator over the unit cell Z, we obtain that (42a) (42a) (42b) (42b) (42c) (42c)
Similarly to what was said in the previous Section, the existence and uniqueness of a solution for the problem given by Eqs. Equation(42a–42c) is guaranteed (see Proposition 1 in Tsalis, Chatzigeorgiou, and Charalambakis [Citation9] for more details). Then, a general solution can be given by (43) (43)
By substituting Eq. (Equation43(43) (43) ) in Eqs. Equation(42a–42c), the third order tensor χklm, periodic in y, is the solution of the following problem referred to as the -cell problem, (44a) (44a) (44b) (44b) (44c) (44c)
Analogously to Eq. (Equation32(32) (32) ), the expression in Eq. Equation(44a)(44a) (44a) is defined, (45) (45)
Problem for
Equating in powers of (46) (46)
Averaging EquationEq. (46)(46) (46) over the unit cell Y, the homogenized problem is obtained and it can be written in the form (47) (47) where (48) (48) is the expression for the effective relaxation modulus of the hierarchical composite media. Furthermore, from Eqs. Equation(15b–15d), the boundary conditions associated to the problem Eq. (Equation47(47) (47) ) are (49) (49) (50) (50) and the initial condition is given by (51) (51)
Remark 3.
In the previous sections and as a consequence of the application of the three-scale Asymptotic Homogenization Method, we were able to obtain the homogenized problem Eqs. Equation(47–51) that describes the overall behavior of a hierarchical, non-aging, linear viscoelastic composite material with generalized periodicity; the respective cell problems Eqs. Equation(31a–31c) and Eqs. Equation(44a–44c) and the effective coefficients Eqs. (Equation35(35) (35) ) and (Equation48(48) (48) ) for each level of organization. The procedure can be used in order to solve a linear elastic problem (t = 0) in composites with the same structural characteristics. Also, we can recover the classical results of two-scale as particular cases of the formulation (see e.g. Cruz-González et al. [Citation17]).
4. Effective coefficients for hierarchical laminated composites with generalized periodicity
Laminated composites materials can be described by using a stratified function such that with n = 2, 3 (see Ref. [Citation9]). Here, we consider the general case when the stratified functions are given by i.e. with Taking into account Voigt’s notation, the symmetry properties of the relaxation modulus and expressions Eqs. (Equation32(32) (32) ) and (Equation45(45) (45) ), the effective coefficients Eqs. (Equation35(35) (35) ) and (Equation48(48) (48) ) can be written, respectively, as follows (52a) (52a) (52b) (52b) for
Additionally, the local problem Eq. Equation(31a)(31a) (31a) can be transformed into the following partial differential equations system by using the Voigt’s notation (53) (53) where and the expressions for the components are (54a) (54a) (54b) (54b) (54c) (54c) (54d) (54d)
Analogously, the local problem Eq. Equation(44a)(44a) (44a) is expressed as follows, (55) (55) where (56a) (56a) (56b) (56b) (56c) (56c) (56d) (56d)
Systems of Eqs. Equation(53–54d) and Equation(55–56d) can be solved analytically by integrating each equation with respect to the local variables and by determining the constants of integration using the corresponding cell average operator (see Ref. [Citation35] for more details). The resolution scheme to obtain the effective viscoelastic properties of a hierarchical composite can be summarized as following:
To solve the system Eqs. Equation(53–54d), and then the expressions for with i = 1, 2, 3, can be substituted into Eq. Equation(52a)(52a) (52a) .
At this point, it is possible to solve the system given in Eq. Equation(55–56d) and to obtain the expressions for with j = 1, 2, 3.
Finally, the effective viscoelastic properties can be obtained by using formula Eq. Equation(52b)(52b) (52b) .
It is worthy to remark that the expressions to the local problems Eqs. Equation(53–54d) and Eqs. Equation(55–56d) are developed for anisotropic composites.
5. An approach for modeling the mechanical properties of the dermis
The dermis is considered as a multilayer collagen-reinforced structure (see Ref. [Citation27]). Two of its principal components are collagen and ground substance (also known as supporting matrix). On the one hand, the collagen contributes to 75% of the fat-free dry mass, 18%-30% of the volume of the dermis and can be considered to behave as a linear elastic material (see Ref. [Citation36, Citation37]). On the other hand, the ground substance is responsible for the viscoelastic behavior of the dermis and comprises about 20% of the dry weight of skin and makes up to between 70% and 90% of the skin’s volume. It is a gel-like substance containing a class of chemicals including glycosaminoglycans (GAG), proteoglycans and glycoproteins (see Ref. [Citation37]). In this section, an approach for modeling the dermis as a hierarchical viscoelastic composite material is proposed. We follow the structural ideas depicted in Figure 3 of Sherman et al. [Citation21]. It is stated that the stiffness of collagenous materials is a consequence of the properties, arrangement, and geometric distribution of collagen fibrils, which are embedded in the ground substance (see Ref. [Citation36]) and forming the dermis. These structures are long strands with wavy effects and are arranged in wavy parallel fibers, which are flattened in the plane of the dermis and determine its properties.
In this work, we present the problematic from the point of view of laminated composite materials (as represented in ). We notice that, in doing this, we take inspiration in the work of Sherman et al. [Citation21]. In fact, (ii) corresponds to Figure 3(b) of Sherman et al. [Citation21] and corresponds to Figure 3(d) of Sherman et al. [Citation21], respectively.
Moreover, we consider the dermis as a two-phase composite made of collagen (elastic constituent) and ground substance (viscoelastic constituent) organized in a hierarchical form. The mechanical properties of the ground substance can be taken approximately from McBride et al. [Citation38]; Panchal et al. [Citation25]; Xu and Lu [Citation37] and the collagen’s mechanical properties are obtained from Benítez and Montáns [Citation27]; Xu and Lu [Citation37]. They are listed in .
The viscoelastic constituent (ground substance) is modeled using the normalized Prony series (see Ref. [Citation36]), (57) (57) where G(t) is the time dependent relaxation shear modulus, τn denote the relaxation times, Gn represent the modulus coefficients and is the long-term modulus (see ). For the sake of simplicity in the model, only two term in the Prony series are considered (see Ref. [Citation36]).
The notation Vi with i = 1, 2 represents the volumetric fractions of each constituent (see ). The values of the volume fraction in the -structural level for the fiber and ground substance (see (i)) are found in pages 10 and 14 of Xu Citation2011, respectively). However, in the case of the values of the volume fraction at the -structural level for the fibril and ground substance see (ii)), to the best of our knowledge and understanding, there is no experimental data available and therefore we take approximate values having into account the transmission electron micrograph of collagen fibrils shown in Figure 3(b) of Sherman et al. [Citation21].
The stratification functions (58) (58) describe the wavy effect in the meso- and nano-structures, respectively (see ).
5.1. Results and discussion
The previously developed methodology and the corresponding data for the dermis constituents allow to estimate the effective relaxation module and creep compliance of this biological tissue. In the formulation, the dermis was considered as a non-aging linear viscoelastic and hierarchical composite with isotropic constituents. The presence of a generalized periodicity at the different length scales, determined by the stratified functions with and the consideration of x2, x3 as principal axes (see (ii)) have an influence in the anisotropic character of the structure. In this sense, the homogenized viscoelastic tensor given in Eq. (52b) has a monoclinic symmetry (13 independent viscoelastic coefficients). This result is not true when the stratification function coincides with the identity function where it is obtained instead a material with transversal isotropic symmetry.
For the completeness of the analysis, we present the mathematical relationship between the effective relaxation modulus and the effective creep compliance, given in the Laplace-Carson space, (59) (59) where I is the fourth order identity tensor. The expression Eq. (Equation59(59) (59) ) is reported in Hashin [Citation39] and it is used to compute the effective creep compliance once is known.
The final results in the temporal space are displayed in and . The MATLAB’s functions INVLAP and GAVSTEH developed by Hollenbeck [Citation40] and Srigutomo [Citation41], respectively are used in the inversion of Laplace-Carson transform. The algorithms can transform functions of complex variable where α is a real exponent. They can also transform functions which contain rational, irrational and transcendent expressions. As a negative aspect, they present problems close to zero.
In , we focus on the behavior of some chosen coefficients: As observed, the curves present an asymptotic behavior, which have a good agreement with the physical relaxation and creep response in viscoelastic materials for long times. Moreover, the monotony changes and the convexities for and in the same components, for instance, and are opposite. This behavior is a consequence of the Eq. Equation(59)(59) (59) and the fact that and are inverse tensors. Similar patterns have been remarked in other works (see Ref. Citation42, Citation43]).
6. Conclusions
A general methodology in the use of the three-scale Asymptotic Homogenization Method (AHM) is considered for modeling non-aging, quasi-periodic, hierarchical and linear viscoelastic composite materials. The present work takes into account the stratified functions in the homogenization procedure allowing the study with more general periodic structures. The analytical solution for the local problems and the expressions of the effective coefficients for hierarchical laminated composites with anisotropic components and perfect contact at the interfaces are derived. Also, the interconnection between the effective relaxation modulus and the effective creep compliance is performed. The approach is applied to study the overall viscoelastic behavior of the dermis. The results shown in and could have a special interest in clinical applications, for instance, in the study of the epidural anesthesia procedure and in the development of haptic devices used in surgical operations. Its study is important in order to provide a location as accurate as possible for the insertion of the needle and the mechanical response of the skin (see Ref. [Citation44]). The results can be used in the automation of the needle insertion procedure (see Ref. [Citation45]).
The theoretical framework that we have developed here can be applied to large variety of hierarchical tissues which exhibit a viscoelastic mechanical response and are characterized by a complex, tortuous geometry, such as solid tumors, see, e.g. Netti et al. [Citation46]; Penta and Ambrosi [Citation47]; Ramírez-Torres et al. [Citation48]. A further development also includes the generalization to a viscoplastic framework by considering remodeling aspects Ramírez-Torres et al. [Citation31]. Although in the most general case numerical simulations cannot be avoided, application of our generalized framework would enable us to encode most of the geometrical complexity analytically, rather than enforce it directly in the setup of the computational grid, thus leading to a substantial reduction of the computational cost.
Also, the article is open to several improvements; for instance, with the framework established to layered materials, it is possible to work with more realistic geometries (e.g. fiber-reinforced composites) by using the stratified functions in the general form with (see Ref. [Citation10]). It represents a novel point of view in order to make comparisons with framework entirely developed to fiber elastic [Citation12] and viscoelastic [Citation16, Citation42] reinforced composites. Besides, the results can be extended to the case of imperfect contact conditions by following the methodology in Guinovart-Sanjuán et al. [Citation35] and considering viscoelastic interfaces (see e.g. Ref. [Citation49]) for each structural level. Another natural step is to generalize the present work to aging viscoelastic solids, see e.g. Sanahuja [Citation50]. This way, the model will be applicable to much larger variety of biophysical (diseased and healthy) scenarios of interest.
Acknowledgements
OLCG kindly thanks to Ecole Doctorale no. 353 de L’Université de Aix Marseille and L’équipe Matériaux & Structures du Laboratoire de Mécanique et d’Acoustique LMA - UMR 7031 AMU - CNRS - Centrale Marseille 4 impasse Nikola Tesla CS 40006 13453 Marseille Cedex 13, France. RRR thanks to Departamento de Matemáticas y Mecánica, IIMAS and PREI-DGAPA at UNAM, for its support to his research project and the Aix-Marseille Université. The authors are thankful to Ana Pérez Arteaga and Ramiro Chávez Tovar for technical assistance. ART kindly acknowledges the Dipartimento di Scienze Matematiche (DISMA) “G.L. Lagrange” of the Politecnico di Torino, “Dipartimento di Eccellenza 2018–2022”, Project code: E11G18000350001. RP is partially supported by EPSRC grant (EP/S030875/1). Also, FJS thanks funding from DGAPA, UNAM.
References
- I. Sevostianov and A. Giraud, Generalization of Maxwell homogenization scheme for elastic material containing inhomogeneities of diverse shape, Int. J. Eng. Sci., vol. 64, pp. 23–36, 2013. DOI: https://doi.org/10.1016/j.ijengsci.2012.12.004.
- Y. Davit, et al., Homogenization via formal multiscale asymptotics and volume averaging: How do the two techniques compare? Adv. Water Resource., vol. 62, pp. 178–206, 2013. DOI: https://doi.org/10.1016/j.advwatres.2013.09.006.
- D. Lukkassen and G. W. Milton, On hierarchical structures and reiterated homogenization. Proceedings of the International Conference in Honour of Jaak Peetre on his 65th birthday: Function Spaces, Interpolation Theory and Related Topics [Series: De Gruyter Proceedings in Mathematics], pp. 355–368, August 17–22, Lund, Sweden, 2000.
- R. Penta and A. Gerisch, The asymptotic homogenization elasticity tensor properties for composites with material discontinuities, Continuum Mech. Thermodyn., vol. 29, no. 1, pp. 187–206, 2017. DOI: https://doi.org/10.1007/s00161-016-x.
- A. Ramírez-Torres, S. D. Stefano, A. Grillo, R. Rodríguez-Ramos, J. Merodio, and J. Penta, Three scales asymptotic homogenization and its application to layered hierarchical hard tissues, Int. J. Solids Struct., vol. 130–131, pp. 190–198, 2018a. DOI: https://doi.org/10.1016/j.ijsolstr.2017.09.035.
- A. Ramírez-Torres, P. Penta, R. Rodríguez-Ramos, and A. Grillo, Homogenized out-of-plane shear response of three-scale fiberreinforced composites, Comput. Visual. Sci., vol. 20, no. 3–6, pp. 85–93, 2019a. DOI: https://doi.org/10.1007/s00791-018-0301-6.
- Z. Yang, Y. Sun, Y. Liu, T. Guan, and H. Dong, A three-scale asymptotic expansion for predicting viscoelastic properties of composites with multiple configuration, Eur. J. Mech. A. Solids., vol. 76, pp. 235–246, 2019b. DOI: https://doi.org/10.1016/j.euromechsol.2019.04.016.
- H. Dong, X. Zheng, J. Cui, Y. Nie, Z. Yang, and Z. Yang, High-order three-scale computational method for dynamic thermo-mechanical problems of composite structures with multiple spatial scales, Int. J. Solids Struct., vol. 169, pp. 95–121, 2019. DOI: https://doi.org/10.1016/j.ijsolstr.2019.04.017.
- D. Tsalis, G. Chatzigeorgiou, and N. Charalambakis, Homogenization of structures with generalized periodicity, Composites Part B. Eng., vol. 43, no. 6, pp. 2495–2512, 2012. DOI: https://doi.org/10.1016/j.compositesb.2012.01.054.
- D. Tsalis, T. Baxevanis, G. Chatzigeorgiou, and N. Charalambakis, Homogenization of elastoplastic composites with generalized periodicity in the microstructure, Int. J. Plasti., vol. 51, pp. 161–187, 2013. DOI: https://doi.org/10.1016/j.ijplas.2013.05.006.
- R. Penta, K. Raum, Q. Grimal, S. Schrof, and A. Gerisch, Can a continuous mineral foam explain the stiffening of aged bone tissue? A micromechanical approach to mineral fusion in musculoskeletal tissues, Bioinspir. Biomim., vol. 11, no. 3, pp. 035004, 2016. DOI: https://doi.org/10.1088/1748-3190/11/3/035004.
- A. Ramírez-Torres, et al., Effective properties of hierarchical fiber-reinforced composites via a three-scale asymptotic homogenization approach, Math. Mech. Solid., vol. 24, no. 11, pp. 3554–3574, 2019b. DOI: https://doi.org/10.1177/1081286519847687.
- G. Allaire and M. Briane, Multiscale convergence and reiterated homogenisation, Proc. R. Soc. Edin. Sec. A. Math., vol. 126, no. 2, pp. 297–342, 1996. DOI: https://doi.org/10.1017/S0308210500022757.
- A. Bensoussan, G. Papanicolau, and J.-L. Lions, Asymptotic Analysis for Periodic Structures, 1st ed., Vol. 5. Elsevier, North-Holland, 1978.
- D. Trucu, M. Chaplain, and A. Marciniak-Czochra, Three-scale convergence for processes in heterogeneous media, Applicable Anal., vol. 91, no. 7, pp. 1351–1373, 2012. DOI: https://doi.org/10.1080/00036811.2011.569498.
- Z. Yang, Y. Sun, J. Cui, and J. Ge, A three-scale asymptotic analysis for ageing linear viscoelastic problems of composites with multiple configurations, Appl. Mathe. Model., vol. 71, pp. 223–242, 2019a. DOI: https://doi.org/10.1016/j.apm.2019.02.021.
- O. L. Cruz-González, et al., An approach for modeling non-ageing linear viscoelastic composites with general periodicity, Composite Structures., vol. 223, pp. 110927, 2019. DOI: https://doi.org/10.1016/j.compstruct.2019.110927.
- Y.-M. Yi, S.-H. Park, and S.-K. Youn, Asymptotic homogenization of viscoelastic composites with periodic microstructures, Int. J. Solids Struct., vol. 35, no. 17, pp. 2039–2055, 1998. DOI: https://doi.org/10.1016/S0020-7683(97)00166-2.
- S. Liu, K.-Z. Chen, and X.-A. Feng, Prediction of viscoelastic property of layered materials, Int. J. Solids Struct., vol. 41, no. 13, pp. 3675–3688, 2004. DOI: https://doi.org/10.1016/j.ijsolstr.2004.01.015.
- A. B. Tran, J. Yvonnet, Q.-C. He, C. Toulemonde, and J. Sanahuja, A simple computational homogenization method for structures made of linear heterogeneous viscoelastic materials, Comput. Method. Appl. Mech. Eng., vol. 200, no. 45–46, pp. 2956–2970, 2011. DOI: https://doi.org/10.1016/j.cma.2011.06.012.
- V. R. Sherman, Y. Tang, S. Zhao, W. Yang, and M. A. Meyers, Structural characterization and viscoelastic constitutive modeling of skin, Acta Biomater., vol. 53, pp. 460–469, 2017. DOI: https://doi.org/10.1016/j.actbio.2017.02.011.
- A. N. Annaidh, K. Bruyère, M. Destrade, M. D. Gilchrist, and M. Otténio, Characterization of the anisotropic mechanical properties of excised human skin, J. Mech. Beh. Biomed. Materi., vol. 5, no. 1, pp. 139–148, 2012. DOI: https://doi.org/10.1016/j.jmbbm.2011.08.016.
- D. Malhotra, et al., Linear viscoelastic and microstructural properties of native male human skin and in vitro 3D reconstructed skin models, J. Mech. Behav. Biomed. Mater., vol. 90, pp. 644–654, 2019. DOI: https://doi.org/10.1016/j.jmbbm.2018.11.013.
- W. Yang, et al., On the tear resistance of skin, Nat. Commun., vol. 6, no. 6649, pp. 6, 2015. DOI: https://doi.org/10.1038/ncomms7649.
- R. Panchal, L. Horton, P. Poozesh, J. Baqersad, and M. Nasiriavanaki, Vibration analysis of healthy skin: toward a noninvasive skin diagnosis methodology, J. Biomed. Opt., vol. 24, no. 1, 2019. DOI: https://doi.org/10.1117/1.JBO.24.1.015001.
- I. Muha, A. Naegel, S. Sabine, A. Grillo, M. Heisig, and G. Wittum, Effective diffusivity in membranes with tetrakaidekahedral cells and implications for the permeability of human stratum corneum, J. Membr. Sci., vol. 368, no. 1-2, pp. 18–25, 2011. DOI: https://doi.org/10.1016/j.memsci.2010.10.020.
- J. M. Benítez and F. J. Montáns, The mechanical behavior of skin: Structures and models for the finite element analysis, Comput. Struct., vol. 190, pp. 75–107, 2017. DOI: https://doi.org/10.1016/j.compstruc.2017.05.003.
- Z. Hashin, Viscoelastic Behavior of Heterogeneous Media, J. Appl. Mech., vol. 32, no. 3, pp. 630–636, 1965. DOI: https://doi.org/10.1115/1.3627270.
- Q.-D. To, S.-T. Nguyen, G. Bonnet and M.-N. Vu, Overall viscoelastic properties of 2D and two-phase periodic composites constituted of elliptical and rectangular heterogeneities, Eur. J. Mech. A. Solids., vol. 64, pp. 186–201, 2017. DOI: https://doi.org/10.1016/j.euromechsol.2017.03.004.
- J. L. Auriault, C. Boutin, and C. Geindreau, Homogenization of Coupled Phenomena in Heterogenous Media, Wiley, New York, 2009.
- A. Ramírez-Torres, et al., An asymptotic homogenization approach to the microstructural evolution of heterogeneous media, Int. J. Non- Linear Mech., vol. 106, pp. 245–257, 2018b. DOI: https://doi.org/10.1016/j.ijnonlinmec.2018.06.012.
- R. M. Christensen, Theory of Viscoelasticity - 2nd Edition an Introduction, Academic Press, Cambridge, MA, 1982.
- S. Maghous and G. J. Creus, Periodic homogenization in thermoviscoelasticity: case of multilayered media with ageing, Int. J. Solids Struct., vol. 40, no. 4, pp. 851–870, 2003. DOI: https://doi.org/10.1016/S0020-7683(02)00549-8.
- E. Sanchez-Palencia, Non-Homogeneous Media and Vibration Theory. Lecture Notes in Physics. Berlin: Springer-Verlag, 1980.
- D. Guinovart-Sanjuán, et al., Analysis of effective elastic properties for shell with complex geometrical shapes, Composite Struct., vol. 203, pp. 278–285, 2018. DOI: https://doi.org/10.1016/j.compstruct.2018.07.036.
- H. Joodaki and M. B. Panzer, Skin mechanical properties and modeling: A review, Proc. Inst. Mech. Eng. H., vol. 232, no. 4, pp. 323–343, 2018. DOI: https://doi.org/10.1177/0954411918759801.
- F. Xu and T. Lu, Introduction to Skin Bio-Thermo Mechanics and Thermal Pain, Vol. 7. New York: Springer, 2011.
- A. McBride, S. Bargmann, D. Pond, and G. Limbert, Thermoelastic modelling of the skin at finite deformations, J. Therm. Biol., vol. 62, pp. 201–209, 2016. DOI: https://doi.org/10.1016/j.jtherbio.2016.06.017.
- Z. Hashin, Theory of fiber reinforced materials. NASA Contractor Report, Washington DC, USA, Rep. NASA CR-1974, 1972.
- K. J. Hollenbeck, Invlap. M: A Matlab function for numerical inversion of Laplace transforms by the Hoog algorithm. 1998. http://www.isva.dtu.dk/staff/karl/invlap.html.
- W. Srigutomo, Gaver-Stehfest algorithm for inverse Laplace transform, 2006. www.mathworks.com/matlabcentral/fileexchange/9987.
- M. A. Araújo-Cavalcante and S. P. Cavalcanti-Marques, Homogenization of periodic materials with viscoelastic phases using the generalized FVDAM theory, Comput. Mater. Sci., vol. 87, pp. 43–53, 2014. DOI: https://doi.org/10.1016/j.commatsci.2014.01.053.
- M. A. Araújo-Cavalcante and S. P. Cavalcanti-Marques, Microstructure effects in wavy-multilayers with viscoelastic phases, Eur. J. Mech. A. Solids., vol. 64, pp. 178–185, 2017. DOI: https://doi.org/10.1016/j.euromechsol.2017.03.003.
- D. Gao, Y. Lei, and B. Yao, Analysis of dynamic tissue deformation during needle insertion into soft tissue, IFAC Proc., vol. 46, no. 5, pp. 684–691, 2013.
- N. Abolhassani, R. Patel, and M. Moallem, Experimental study of robotic needle insertion in soft tissue, Int. Congr. Ser., vol. 1268, pp. 797–802, 2004. DOI: https://doi.org/10.1016/j.ics.2004.03.110.
- P. A. Netti, D. A. Berk, M. A. Swartz, A. J. Grodzinsky, and R. K. Jain, Role of extracellular matrix assembly in interstitial transport in solid tumors, Cancer Res., vol. 60, no. 9, pp. 2497–2503, 2000.
- R. Penta and D. Ambrosi, The role of the microvascular tortuosity in tumor transport phenomena, J. Theoret. Biol., vol. 364, pp. 80–97, 2015. DOI: https://doi.org/10.1016/j.jtbi.2014.08.007.
- A. Ramírez-Torres, et al., The role of malignant tissue on the thermal distribution of cancerous breast, J. Theoret. Biol., vol. 426, pp. 152–161, 2017. DOI: https://doi.org/10.1016/j.jtbi.2017.05.031.
- L. Daridon, C. Licht, S. Orankitjaroen, and S. Pagano, Periodic homogenization for Kelvin-Voigt viscoelastic media with a Kelvin-Voigt viscoelastic interphase, Eur. J. Mech. A. Solids., vol. 58, pp. 163–171, 2016. DOI: https://doi.org/10.1016/j.euromechsol.2015.12.007.
- J. Sanahuja, Effective behaviour of ageing linear viscoelastic composites: Homogenization approach, Int. J. Solids Struct., vol. 50, no. 19, pp. 2846–2856, 2013. DOI: https://doi.org/10.1016/j.ijsolstr.2013.04.023.