Research Article | Open Access

# A Nonmonotonous Damage Model to Characterize Mullins and Residual Strain Effects of Rubber Strings Subjected to Transverse Vibrations

**Academic Editor:**Jun Zhang

#### Abstract

This work focuses on the formulation of a constitutive equation to predict Mullins and residual strain effects of buna-N, silicone, and neoprene rubber strings subjected to small transverse vibrations. The nonmonotone behavior exhibited by experimental data is captured by the proposed material model through the inclusion of a phenomenological non-monotonous softening function that depends on the strain intensity between loading and unloading cycles. It is shown that theoretical predictions compare well with uniaxial experimental data collected from transverse vibration tests.

#### 1. Introduction

It is well known that the response of rubber-like materials when subjected to cyclic loading, the load needed to stretch the material sample to a given stretch value during the second loading cycle in tension, compression, or shear deformation states, to say a few, is smaller than that required to reach the same stretch during the primary loading cycle. This reduction in the stress magnitude is known as the Mullins effect [1, 2]. To characterize this effect, there exist in the literature many different material models that have been developed to explain material damage mechanisms. See, for instance, the papers of Govindjee and Simo [3], Ogden and Roxburgh [4], Johnson and Beatty [5], Beatty and Krishnaswamy [6], Elías-Zúñiga and Beatty [7], Holzapfel et al. [8], de Souza Neto et al. [9], Elías-Zúñiga [10], Dorfmann and Ogden [11], Marckmann et al. [12], Diani et al. [13], Kazakevičiūtė-Makovska and Kačianauskas [14], de Tommasi et al. [15], Rickaby and Scott [16], Córdova et al. [17], and references cited therein for an overview of the main features of these models. Another effect that has been observed in rubber-like materials during the unloading process is related to the so-called healing effect in which a fraction of the broken cross-links can be reformed [18]. It is also known that the fraction percentage of reformed cross-links depends on the maximum and minimum strains of the loading cycles. Due to this healing effect, the material behavior upon unloading-reloading cycles is not elastic because of the hysteresis effects [19]. Other phenomena such as stress-induced anisotropy as well as crystallization in cyclic loading of polymeric-like materials are discussed in [13] and [20], respectively.

In the present work, we neglect healing and anisotropic material effects and use a constitutive material model to characterize softening and residual strain effects of rubber-like cords subjected to small transverse vibrations test. In particular, we characterize the nonmonotonous behavior shown by experimental data by using a nonmonotonous stress-softening function that depends on the amount of strain intensity. To consider residual strain effects, the procedure developed by Holzapfel et al. in [8] is used to develop a strain energy density function which depends on a damage parameter variable that becomes active only during unloading conditions. We will show that the nonmonotonous damage function describes well experimental data when the normalized transverse frequency ratio is plotted against the normalized strain intensity ratio [10, 21].

The paper is organized as follows. In Section 2, we begin with a brief review of the required equations for finite deformations of an incompressible elastic material. In Section 3, we investigate, by introducing a nonmonotonous damage softening function, the nonmonotonous behavior of experimental data when normalized uniaxial stresses and transverse vibration frequency ratios are plotted against the normalized strain intensity ratio [10, 21]. In Section 4, we have used the Ogden and Roxburgh pseudoelasticity theory to account for residual strain effects by introducing slight modifications to the model developed by Holzapfel and coworkers in [8]. In Section 5, we use the non-Gaussian average-stretch, full-network model of arbitrarily oriented molecular chains [22] to derive the constitutive equations that describe the small transverse vibrations of rubber cords. A comparison of the results corresponding to simulated and experimental data is done in Section 6. Finally, in Section 7, we make some closing remarks related to our derived stress-softening model that accounts for residual strains.

#### 2. Theoretical Aspects

We consider the deformation of an incompressible elastic body which in its natural configuration occupies the region . A material particle is assumed to be in its undeformed reference configuration of a body at the place . After a prescribed deformation, the body occupies the region , the current configuration, and the particle at moves to the place in a common rectangular Cartesian frame with origin and orthonormal basis . Thus, the Cauchy-Green deformation tensor has the form where , are the associated orthonormal principal directions, is the usual deformation gradient, and denote the principal stretches in . Note that the magnitude of the strain intensity at a material point denoted by , is defined by , where tr denotes the trace operation. In the undeformed state , the identity tensor and ; otherwise, for all isochoric deformations [6]. We note that for all , equality holding when and only when , the undeformed state. Since the principal invariants of are defined by the magnitude of the strain intensity is given as

#### 3. Characterization of the Nonmonotonous Material Behavior

For an incompressible and isotropic elastic material, the corresponding time-independent Cauchy stress constitutive equation has the form in which is the Cauchy stress, is an undetermined pressure, and , , denote the virgin material response functions related to the strain energy function , per unit reference volume, in accordance with wherein [22]. By using (4), Elías-Zúñiga and Beatty in [7] proposed a damage type model to describe the stress-softened material behavior of the form in which denotes the Cauchy stress in the stress-softened material, denotes the maximum previous strain at which the material is unloaded from the primary path, and is an isotropic softening function at the damage level on the interval . They assumed that this softening function is a monotone increasing function of the strain intensity that satisfies the conditions

Based on this assumption, Elas-Zúñiga and Beatty proposed the following softening function: where is a dimensionless positive material softening parameter. After substituting (8) into (6), Elas-Zúñiga and Beatty obtained the following stress-softening phenomenological material model:

Theoretical predictions provided by (9) were computed in [7, 10] and compared to experimental data for uniaxial, pure shear, and equibiaxial deformation states. There, the theoretical predictions showed reasonably good agreement with experimental data not only for the virgin loading path, but also for the reloading paths. However, in a recent publication Kazakevičiūtė-Makovska [23] observed that the experimental data when plotted as the normalized stress components on versus the stretch ratio showed nonmonotonous behavior with the characteristics S-shaped form and she concluded that, because of the variations in shapes of the curves for different deformation cycles, different values of the softening parameters were needed to fit experimental data for a particular choice of the softening function. Moreover, she showed that the softening function given by (8) fails in predicting the nonmonotonous behavior exhibited by experimental data collected by Mullins and Tobin [2], Cheng and Chen [24], and Mars and Fatemi [25] at higher stretch values. In this sense, de Tommasi et al. showed the importance of microscopic inhomogeneity to describe known experimental effects observed in amorphous materials such as the transition from diffuse to localized damage as the distribution properties are varied [26]. In fact, they showed that the monotone stress-stretch loading curve behavior is mainly due to a diffuse damage mechanism. They also considered that amorphous materials can be characterized by unstable strain domains, which gives the possibility of having homogeneous or localized damage with nonmonotone primary loading curve.

To confirm these observations, we use uniaxial tension data from a strain-controlled experiment (buna6) by Johnson on a buna-N elastomer in which, for each loading data point, there is a corresponding data point for the same amount of stretch on the unloading curve [7]. Figure 1 illustrates the normalized stress components versus the normalized strain intensity ratio curves in which residual strains are considered. We can see from Figure 1 that the experimental data exhibits monotone behavior. Nevertheless, there is a wide variation between theoretical predictions provided by the corresponding uniaxial stress component given by (9) and experimental data, especially for small values of near to . Next, we use Elías-Zúñiga and Beatty strain-controlled experimental data on silicone elastomer cords subjected to the small amplitude transverse vibrations and plot the transverse vibrational frequency data ratio, , versus the normalized strain intensity ratio, where represents the transverse vibrational frequency of the virgin string, and the corresponding frequency of the stress-softening cord [21]. As we can see from Figure 2, the collected experimental data exhibit nonmonotone behavior. Therefore, it is clear that theoretical predictions computed from (9) must deviate from experimental data since this equation considers the softening function (8) which has monotonous behavior. In Figures 1 and 2, the color dots represent experimental data and the dashed black lines are theoretical predictions. These observations and the conclusions drawn by Kazakevičiūtė-Makovska in [23] and de Tommasi et al. in [26] motivate us to consider as a nonmonotone softening function that must satisfy the conditions given by (7). Therefore, we use here the nonmonotonous softening function [27] where is a positive softening material parameter and and are positive scaling constants chosen to best fit experimental data for a given rubber-like material. Then, the substitution of (10) into (6) yields the following material model:

Notice that (11) does not take into account permanent set effects. Therefore, in the next section, we will use pseudoelasticity theory to obtain an expression that will include residual strain effects.

#### 4. Inclusion of Residual Strain Effects

Hyperelastic materials are rheologically described by employing strain energy density functions which are related to the energy stored in the material; there are several hyperelastic material models proposed to predict elastomeric material behavior, but most of them do not take into account residual strains [4–7]. Holzapfel and coworkers proposed the following phenomenological model to include residual strain effects [8]:

Here, represents the strain energy function associated with the primary loading path that is supplemented by a damage mechanism energy function which depends on the principal stretches , and the discontinuous damage variables , , which are related to the maximum principal stretch values; thus, the modified strain energy function depends on the principal stretches. Based on this modified strain energy function , Holzapfel and coworkers derived a stress-softening material model by using the theory of pseudoelasticity developed in [4] and implemented it into a finite element program. Later, Dorfmann and Ogden proposed a model based on the theory of pseudoelasticity to capture the stress-softening and residual strain effects [28]. This model uses the monotonic dissipation energy functions and that are assumed to depend on the amount of stretch and evolve with the deformation material history. This model has the form where is a function that characterizes the residual strains, is a damage or softening variable, represents the residual strain variable, and represents the material strain energy density when and are inactive.

Here, we follow Holzapfel et al. phenomenological model given by (12) and assume that the residual strain energy that accounts for the damage material mechanism has the form where is the material shear modulus, represents the damage function, ,, are the maximum values of the principal stretches at which unloading begins on the primary loading path, and is a positive scaling constant that in general takes the value of one. By using the concepts of pseudo-elastic theory [4], we assume that [17]

Substitution of (14) into (15) yields from which

If is assumed to have the form then the integration of (18) provides where is a positive dimensionless material constant and is an integration constant. Notice from (17) and (18) that

On the primary loading path, is inactive, while on the unloading path, has the value given by (20). Substitution of (20) into (14) gives the following pseudo-elastic residual strain energy density per unit volume:

Therefore, the strain energy function that accounts for residual strains of an incompressible, isotropic, and hyperelastic material, in accordance with pseudo-elastic theory, has the form

In the next section, we will use (22) to derive a material model that will be used to characterize experimental data of elastomeric rubber cords subjected to small transverse vibrations.

#### 5. A Nonmonotonous Material Model

To illustrate the application of the softening function given by (10) as well as of the residual strain energy function (22), we first consider the Arruda-Boyce constitutive equation for an average-stretch, full-network of arbitrarily oriented molecular chains and modify this to include stress-softening and residual strain effects. Then, we will use Taylors’ equation to derive the material model that describes the dynamics response of rubber-like cords when subjected to small amplitude transverse vibrations.

First, let us begin with the modification of the average-stretch, full-network material model to include residual strain effects. For this material model, it is well known that the strain energy per unit volume for the loading path is given by where is the relative chain stretch defined by and represents the fully extended chain stretch, is the chain number of rigid links, each of length , and is the chain deformation that in the affine deformation is determined by is the inverse of the Langevin function which is defined as and is a constant that ensures that the strain energy vanishes in the undeformed state [22, 27]. Substitution of (23) into (22) provides the modified non-Gaussian pseudostrain energy per unit volume that accounts for residual strains on the unloading path; that is, where is an energy constant. This is the strain energy density of an elastic material for which is active.

Substituting of (23) into (4) and with the aid of (5), the Cauchy stress-stretch averaging network model components for the virgin material, for which , becomes

Eliminating the pressure from (28) gives where (no sum). The Cauchy stress-stretch constitutive equation for a stress-softening material can be obtained by substituting (27) into (4) and by using (5) and (11). This yields the following stress-softening components: where

Then, on elimination of from (30), we get that where, in general, (no sum).

Recalling that for an incompressible material, the engineering stress is related to the Cauchy stress by
then, the uniaxial *engineering stress-stretch relation for an average-stretch, full-network stress-softening material model that considers residual strain effects is* obtained by using (31), (32), and (33):

Here, and the relative chain stretch which can be obtained from (24) and (25) is given as

In equibiaxial extension, the engineering stress-softened relation is given by [27] where

In the case of a pure shear or plane strain compression deformation state, the engineering stress-softened constitutive equation is given by where

We next use Taylors’ equation [21] to derive the nonmonotonous material constitutive equation that describes the behavior of rubber cords subjected to transverse vibrations.

##### 5.1. Transverse Vibration of Elastomeric Cords

Here, we briefly review some of the fundamental concepts that are used to describe the small amplitude transverse vibrations of rubber-like cords. Full details of the theoretical analysis can be found in [21].

###### 5.1.1. Basic Equations

The frequency of the fundamental mode of the small amplitude transverse vibrations about the stretched equilibrium position of an elastomeric cord with fixed ends is given by Taylor’s equation as where is the initial cord length, represents the uniform cross-sectional area, is the mass density per unit length, and is the Cauchy stress for the uniaxial loading of the virgin material given by (29). In the case of the stress-softening material, the transverse vibrations frequency can be determined by using the following equation: where represents the stress-softening Cauchy stress described by (32). Notice from previous remarks that . Therefore, the magnitude of the transverse vibrations frequency of the stress-softening material cord is always less than the corresponding frequency of the virgin material at each value of the strain intensity , equality holding when and only when . We now introduce the normalized small amplitude transverse vibration frequencies for the virgin and stress-softened materials: and . The corresponding constitutive equation of , for the softening material, can be obtained by substituting (34) into (42); this yields

Of course, the normalized virgin cord frequency is obtained from (43) with and .

#### 6. Theoretical Simulations

To assess the precision of the derived constitutive (43) which includes residual strains and has a nonmonotonous softening function that describes Mullins effect, we use experimental data collected during the small amplitude transverse vibrations of buna-N, silicone, and neoprene elastomers [21]. Notice that only four constitutive material constants and one fitting parameter are required to use (43), that is, the shear modulus , the chain number of links , the stress-softening parameter , the residual strain material constant , and the fitting parameter . In what follows, we will show that the value of for the normalized transverse vibration frequency model described in (43) provides good fit to experimental data.

We begin by using data collected from the transverse vibrations of a buna-N elastomer material cord [21]. Figure 3 illustrates the predicted curves obtained from (43) compare to experimental data. Notice from Figure 3(a) that the normalized transverse frequency predictions are in good agreement with experimental data for the several loading and unloading cycles that exhibit residual strains. The nonmonotonous characteristic behavior of the buna-N experimental data is shown in Figure 3(b) where the ratio of the transverse vibrations frequency of the softening and the virgin material is plotted versus the normalized strain intensity ratio . One can notice that there is some discrepancies between experimental data and theoretical predictions that could be due to viscoelastic material effects that were neglected during the derivation of the material model (43). In this case, the constitutive material constants used to best fit experimental data are Hz, , , and MPa.

**(a)**

**(b)**

Figures 4 and 5 illustrate data of stress softening collected from silicone and neoprene rubber strings, respectively. We can see that predicted curves, for the proposed material model defined by (43), stand in good agreement with experimental data for the several loading and unloading cycles. Moreover, the nonmonotonous behavior and residual strains depicted in Figures 4(b) and 5(b) are described well by our derived material model. In these figures, the dots represent the experimental data and the solid and dashed lines represent theoretical results obtained, respectively, from (29) and (43). The material constants for the material model are provided by a best fit analysis and these are listed in the figure captions.

**(a)**

**(b)**

**(a)**

**(b)**

#### 7. Conclusions

In this paper, we have used a nonmonotonous softening function that accounts for the microstructural material damage upon deformation from the natural, undistorted state of the virgin material. This phenomenological softening function is based on three parameters; one is a positive softening material parameter and the other two are basically positive scaling constants that in general have the values of and . We have shown that this softening function predicts well the nonmonotonous behavior observed in rubber-like materials when experimental data are plotted against the normalized strain intensity ratio. Nevertheless, there is a variation in the results for values of near to , which we believe is due to some viscoelastic and anisotropic effects.

The Arruda-Boyce constitutive equation for an average-stretch, full-network model of arbitrarily oriented molecular chains was modified to include the residual strain effect for the unloading paths as described by the simple constitutive relation (30). Then, we have used Taylors’ equation and derived the constitutive equation that describes the small transverse vibrations of elastomeric cords subjected to uniaxial extension. This model captures well the qualitative and quantitative behavior shown by experimental data collected from samples of buna-N, silicone, and neoprene elastomeric cords. In modeling the small transverse vibrations of the aforementioned rubber cords, we have only used four material constants: , the chain number of rigid links , the material softening parameter , and a positive material constant that is related to the pseudoelastic residual strain energy. Based on the accuracy of our proposed nonmonotonous model (43) to predict experimental data, we can conclude that the extent of damage of rubber-like materials, when subjected to small transverse vibrations, can be conveniently determined by considering the nonmonotonous behavior exhibited by experimental data. Furthermore, our proposed damage material mechanism model given by (14) can be applied to three-dimensional deformation states since the damage function is controlled by the strain energy associated with the primary deformation and not just by the specific deformation from which it is calculated.

Finally, we have also found, by using our theoretical model, excellent agreement with stress-softened experimental data collected for other deformation states. The results of this new work will be reported in a subsequent paper.

#### Acknowledgment

This work was funded by Tecnológico de Monterrey—Campus Monterrey, through the Research Chair in Nanomaterials for Medical Devices and Research Chair in Intelligent Machines and from Consejo Nacional de Ciencia y Tecnología (Conacyt), México.

#### References

- L. J. Mullins, “Effect of stretching on the properties of rubber,”
*Journal of Rubber Research*, vol. 16, pp. 275–289, 1947. View at: Google Scholar - L. Mullins and N. R. Tobin, “Theoretical model for the elastic behavior of filledreinforced vulcanized rubbers,”
*Journal of Rubber Chemestry and Technology*, vol. 30, no. 2, pp. 555–571, 1957. View at: Publisher Site | Google Scholar - S. Govindjee and J. Simo, “A micro-mechanically based continuum damage model for carbon black-filled rubbers incorporating Mullins' effect,”
*Journal of the Mechanics and Physics of Solids*, vol. 39, no. 1, pp. 87–112, 1991. View at: Google Scholar - R. W. Ogden and D. G. Roxburgh, “A pseudo-elastic model for the Mullins effect in filled rubber,”
*Proceedings of the Royal Society A*, vol. 455, no. 1988, pp. 2861–2877, 1999. View at: Google Scholar - M. A. Johnson and M. F. Beatty, “The Mullins effect in uniaxial extension and its influence on the transverse vibration of a rubber string,”
*Continuum Mechanics and Thermodynamics*, vol. 5, no. 2, pp. 83–115, 1993. View at: Publisher Site | Google Scholar - M. F. Beatty and S. Krishnaswamy, “A theory of stress-softening in incompressible isotropic materials,”
*Journal of the Mechanics and Physics of Solids*, vol. 48, no. 9, pp. 1931–1965, 2000. View at: Publisher Site | Google Scholar - A. Elías-Zúñiga and M. F. Beatty, “A new phenomenological model for stress-softening in elastomers,”
*Zeitschrift fur Angewandte Mathematik und Physik*, vol. 53, no. 5, pp. 794–814, 2002. View at: Google Scholar - G. A. Holzapfel, M. Stadler, and R. W. Ogden, “Aspects of stress softening in filled rubbers incorporating residual strains,” in
*Proceedings of the First European Conference on Constitutive Models for Rubber*, A. Dorfmann and A. Muhr, Eds., pp. 189–193, Balkema, Rotterdam, The Netherlands, 1999. View at: Google Scholar - E. A. de Souza Neto, D. Perić, and D. R. J. Owen, “A phenomenological three-dimensional rate-idependent continuum damage model for highly filled polymers: formulation and computational aspects,”
*Journal of the Mechanics and Physics of Solids*, vol. 42, no. 10, pp. 1533–1550, 1994. View at: Google Scholar - A. Elías-Zúñiga, “A phenomenological energy-based model to characterize stress-softening effect in elastomers,”
*Polymer*, vol. 46, no. 10, pp. 3496–3506, 2005. View at: Publisher Site | Google Scholar - A. Dorfmann and R. W. Ogden, “A constitutive model for the Mullins effect with permanent set in particle-reinforced rubber,”
*International Journal of Solids and Structures*, vol. 41, no. 7, pp. 1855–1878, 2004. View at: Publisher Site | Google Scholar - G. Marckmann, E. Verron, L. Gornet, G. Chagnon, P. Charrier, and P. Fort, “A theory of network alteration for the Mullins effect,”
*Journal of the Mechanics and Physics of Solids*, vol. 50, no. 9, pp. 2011–2028, 2002. View at: Publisher Site | Google Scholar - J. Diani, M. Brieu, and J. M. Vacherand, “A damage directional constitutive model for Mullins effect with permanent set and induced anisotropy,”
*European Journal of Mechanics*, vol. 25, no. 3, pp. 483–496, 2006. View at: Publisher Site | Google Scholar - R. Kazakevičiūtė-Makovska and R. Kačianauskas, “Modelling of stress softening in elastomeric materials,”
*Mechanics Research Communications*, vol. 31, no. 4, pp. 395–403, 2004. View at: Publisher Site | Google Scholar - D. de Tommasi, G. Puglisi, and G. Saccomandi, “A micromechanics-based model for the Mullins effect,”
*Journal of Rheology*, vol. 50, no. 4, pp. 495–512, 2006. View at: Publisher Site | Google Scholar - S. R. Rickaby and N. H. Scott, “A cyclic stress softening model for the Mullins effect,”
*International Journal of Solids and Structures*, vol. 50, pp. 111–120, 2013. View at: Google Scholar - R. Córdova, A. Elías-Zúñiga, L. Elizalde et al., “Polymeric materials reinforced with multiwall carbon nanotubes: a constitutive material model,”
*Materials*, vol. 6, no. 7, pp. 2873–2891, 2013. View at: Google Scholar - P. D'Ambrosio, D. de Tommasi, D. Ferri, and G. Puglisi, “A phenomenological model for healing and hysteresis in rubber-like materials,”
*International Journal of Engineering Science*, vol. 46, no. 4, pp. 293–305, 2008. View at: Publisher Site | Google Scholar - D. de Tommasi, S. Marzano, G. Puglisi, and G. Zurlo, “Damage and healing effects in rubber-like balloons,”
*International Journal of Solids and Structures*, vol. 46, no. 22-23, pp. 3999–4005, 2009. View at: Publisher Site | Google Scholar - S. Toki, I. Sics, S. Ran et al., “New insights into structural development in natural rubber during uniaxial deformation by in situ synchrotron X-ray diffraction,”
*Macromolecules*, vol. 35, no. 17, pp. 6578–6584, 2002. View at: Publisher Site | Google Scholar - A. Elías-Zúñiga and M. F. Beatty, “Stress-softening effects in the transverse vibration of a non-Gaussian rubber string,”
*Meccanica*, vol. 38, no. 4, pp. 419–433, 2003. View at: Publisher Site | Google Scholar - M. F. Beatty, “An average-stretch full-network model for rubber elasticity,”
*Journal of Elasticity*, vol. 70, no. 1–3, pp. 65–86, 2003. View at: Publisher Site | Google Scholar - R. Kazakevičiūtė-Makovska, “Experimentally determined properties of softening functions in pseudo-elastic models of the Mullins effect,”
*International Journal of Solids and Structures*, vol. 44, no. 11-12, pp. 4145–4157, 2007. View at: Publisher Site | Google Scholar - M. Cheng and W. Chen, “Experimental investigation of the stress-stretch behavior of EPDM rubber with loading rate effects,”
*International Journal of Solids*, vol. 40, no. 18, pp. 4749–4768, 2003. View at: Google Scholar - W. V. Mars and A. Fatemi, “Observations of the constitutive response and characterization of filled natural rubber under monotonic and cyclic multiaxial stress states,”
*Journal of Engineering Materials and Technology*, vol. 126, no. 1, pp. 19–28, 2004. View at: Publisher Site | Google Scholar - D. de Tommasi, G. Puglisi, and G. Saccomandi, “Localized versus diffuse damage in amorphous materials,”
*Physical Review Letters*, vol. 100, no. 8, Article ID 085502, 4 pages, 2008. View at: Publisher Site | Google Scholar - A. Elías-Zúñiga and C. A. Rodríguez, “A non-monotonous damage function to characterize stress-softening effects with permanent set during inflation and deflation of rubber balloons,”
*International Journal of Engineering Science*, vol. 48, no. 12, pp. 1937–1943, 2010. View at: Publisher Site | Google Scholar - A. Dorfmann and R. W. Ogden, “A pseudo-elastic model for loading, partial unloading and reloading of particle-reinforced rubber,”
*International Journal of Solids and Structures*, vol. 40, no. 11, pp. 2699–2714, 2003. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2013 Alex Elías-Zúñiga et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.