Research Article  Open Access
Shulei Sun, Wenguo Chen, "An Anisotropic Hyperelastic Constitutive Model with Bending Stiffness Interaction for CordRubber Composites: Comparison of Simulation Results with Experimental Data", Mathematical Problems in Engineering, vol. 2020, Article ID 6750369, 7 pages, 2020. https://doi.org/10.1155/2020/6750369
An Anisotropic Hyperelastic Constitutive Model with Bending Stiffness Interaction for CordRubber Composites: Comparison of Simulation Results with Experimental Data
Abstract
Based on the invariant theory of continuum mechanics by Spencer, the strain energy depends on deformation, fiber direction, and the gradients of the fiber direction in the deformed configuration. The resulting extended theory is very complicated and brings a nonsymmetric stress and couple stress. By introducing the gradient of fiber vector in the current configuration, the strain energy function can be decomposed into volumetric, isochoric, anisotropic, and bending deformation energy. Due to the particularity of bending deformation, the reinforced material has tensile deformation and compression deformation. The bending stiffness should be taken into consideration, and it is further verified by the bending simulation.
1. Introduction
The theory of finite deformation of fiberreinforced composites based on the continuum mechanics theory was first proposed and established by Adkins and Rivlin [1]. Spencer used the invariant method [2, 3] to represent the fiber direction as a unit vector in the reference configuration. Spencer further refined the theory of finite deformation of fiberreinforced composites. This method is widely used in the finite deformation of composite materials, especially in the industrial and biomechanical fields, i.e., air compression springs composed of fiberreinforced materials [4] and intervertebral disc tissue in biology, and it achieved success. The results demonstrate the effectiveness of this model. Peng et al. proposed a new constitutive equation [5], which divides the energy function (per unit reference volume) of human intervertebral disc annulus into three parts: matrix, fiber, and the interaction of the two, and determines the material parameters, experiments, and simulation through step by step experiments. The data are in good agreement. Another important application in biomechanics is the establishment of arterial wall constitutive [6]. Steigmann developed a model for the mechanics of woven fabrics in the framework of twodimensional elastic surface theory [7]. The implementation of an enhanced modelling approach for fiberreinforced composites is presented which may, in addition to the directional dependency induced by the fibers, allow the capturing of the fiber bending stiffness [8]. Comparison of simulation results with analytical solutions was made with respect to the fiberreinforced composites with fiber bending stiffness under azimuthal shear [9]. Soldatos researched the invention of methods appropriate for characterization of fiberreinforced materials that exhibit polar material behavior due to fiber bending resistance [10]. A linear model, framed in the setting of the second strain gradient theory, is presented for the mechanics of an elastic solid reinforced with fibers resistant to flexure [11].
However, when this homogenous constitutive model based on continuum mechanics is applied to the bending model, it is found that the constitutive model cannot accurately describe the mechanical behavior of bending [3]. There is a premise that the fiber is infinitely soft under uniaxial compressive condition. This is a valid assumption in most cases, such as stretching and compression. In the conventional theory, there are no parameters in terms of dimensions, and thus it cannot explain the size effect in any aspect. For this reason, Spencer established an invariantbased model that considers the fiber size effect [3], but the model is too complex for practical application.
Based on the invariant theory established by Spencer [12], this paper attempts to establish an anisotropic hyperplastic constitutive model with bending stiffness, i.e., increasing the parameter dependent on the fiber direction vector gradient in current deformation and establishing the relationship with curvature. The comparison between the simulation results and the bending experimental data of the rubbercord composite proves the correctness of the model.
2. RubberCord Constitutive Model
2.1. Constitutive Model of Conventional Anisotropic Hyperelastic FiberReinforced Composites
According to Spencer’s theory of anisotropic fiberreinforced composites [3], the Helmholtz strain energy function (per unit reference volume) depends on the deformation and fiber direction :where the tensor F is deformation gradient t, the tensor is the right Cauchy–Green deformation, and the tensor A () is of order two.
The Helmholtz strain energy function per unit volume can also be written in the form of invariants for C and A:
Since the deformation gradient F can be decomposed into volumepreserving deformation and volumetric deformation, the strain energy function may be expressed as a function of the following three parts [6]:where W_{vol}, W_{iso}, and W_{ani} denote the elastic deformation energy caused by volume deformation, isotropic deformation, and anisotropic deformation, respectively.
is the deformation volume ratio, and , , , and are the corresponding partial invariants:
In order to facilitate parameter fitting, the strain energy can also be given bywhere equivalent to formula (3) W_{vol} + W_{iso}, W_{F} + W_{FM} equivalent to W_{ani} formula (3) W_{ani}. According to the above analysis, for anisotropic fiberreinforced composites that consider compressibility, a common form of strain energy per unit volume has the form [8–10]
By means of chain rule, the secondorder Piola–Kirchhoff stress tensor can be obtained by deriving the strain energy function with respect to the right Cauchy–Green variable C. We write
In particular, we consider the incompressible materials and obtain the formwhere .where ,
The Cauchy stress tensor is defined to be
2.2. Constitutive Model of FiberReinforced Composites with Bending Stiffness
The form of strain energy in a conventional constitutive model depends only on both the deformation gradient F and the fiber direction tensor A. For a material reinforced by a single family of fibers, the fiber direction is defined by a unit vector field A in the reference configuration and a unit vector field a in the deformed configuration. A single family of fibers here mainly refers to the direction of fiber arrangement and single means a direction for fibers. The fibers are convected with the material, and thus
When considering the contribution with the bending stiffness, we need to increase the fiber direction vector b and the tensor G, which is given by . By analogy with relation (12), we may postulate the free energy per unit reference energy:
The free energy must be unchanged if the fiberreinforced composite undergoes a rotation described by the proper orthogonal tensor Q. We obtain the following expressions:where Q is arbitrary orthogonal tensor. The tensors F and G can be written in the form of a scalar product, respectively:
Hence, we can get using the Cayley–Hamilton formulation [11]:
Let the tensor be represented by the invariants of . Finally, the function W that can be represented as is expressed as
According to the formula in Zheng [12], the free energy W consists of 33 invariants (Table 1) about isotropic, which is too complicated. In order to simplify the constitutive equation, it is now assumed that the W curvature of the fiber is more strictly dependent, i.e., the directional derivative deformation gradient of the fiber direction vector A after deformation and the unit vector before deformation F, i.e.,

By analogy with the F and G, the tensor can be written in the form of product, respectively:
This W can be expressed as a function of C, :
According to the representation theory of the tensor function [12], a total of 11 invariants can be expressed, where the first five invariants are the same as in the formula (2).where can be obtained by equation (19). From equation (19),
In summary, the total strain energy with the bending stiffness is in the form ofwhere the free energy W_{G} is the energy contributed to the invariants which is associated with the bending stiffness.
3. Determination of the Strain Energy Function and Corresponding Parameters of VBelt
The Vbelt shown in Figure 1 is a Vbelt widely used in transmissions of automobiles, machines, and machine tools. It has a complicated structure and contains a fabric and a core wire in addition to rubber. The belt has anisotropic characteristics. Ishikawa et al. carried out a series of experiments and simulation studies [13]. The following is mainly based on the experimental data for the parameter fitting and simulation of the constitutive model.
3.1. Strain Energy and Parameters of Matrix Rubber
Matrix rubber is generally regarded as a hyperelastic material. In order to obtain accurate material parameters, the Vbelt is subjected to transverse uniaxial tensile test. The stressstrain relationship can be approximated as the performance of the hyperelastic matrix rubber. Using the data fitting module of the ABAQUS hyperelastic part, considering the compressibility, Poisson’s ratio of the rubber is set to 0.47. The typical Mooney–Rivlin model is used to compare the mechanical properties of the matrix material. The fitting curve is shown in Figure 2. The Mooney–Rivlin strain energy form of the rubber with compressibility is as follows:
The parameters obtained by fitting the transverse uniaxial stretching of the above curve are D = 0.01763 MPa^{−1}, C_{10} = −0.875 MPa, and C_{01} = 4.3478 MPa,
3.2. Determination of the Strain Energy Form and Parameters of the Fiber Part
Tensile strain and fiber elongation have the following approximate relationship:where is the elongation of the fiber; we use equation (23) to obtain the relationship between stress and stress through the forcedisplacement curve and finally express the tensile strain energy as a function . The contribution of energy is related to the fibermatrix interaction. It is written according to formula (6) as follows:
Data fitting was performed using MATLAB programming. The fitted curve is shown in Figure 3.
The final material parameters are C_{2} = 7.9718, C_{3} = −6.7227 MPa, C_{4} = 0.086, and C_{5} = −0.027 MPa. From the results of the above fitting, the free energy contributed by the invariant I_{4} is much larger than the energy contributed by I_{5}. For a more detailed discussion, refer to Peng et al. [5].
3.3. Determination of Strain Energy and Parameters Contributed by Fiber Bending Stiffness
It can be seen from equation (13) that the invariants are all related to the curvature of the fiber. To simplify the formula, it is assumed that the deformation of the fiber is linear elastic, ignoring the term of the quadratic term and above, so that the strain energy contributing to the bending stiffness is only considered. I_{9} which is in equation 20,
Due to the influence of the bending stiffness, the generation of the couple stress and the stress matrix are no longer symmetrical, and MATLAB programming is used to solve the problem. Due to the lack of corresponding experimental data to determine the specific values C_{6}, separate energy function W_{G} is used to study the effect on bending deformation.
4. Numerical Simulation
Combining the material parameters obtained by the previous parameter, the corresponding finite element program is compiled to simulate the bending of the rubbercord composite, i.e., Vbelt as shown in Figure 4. The cuboid sample of size is cut from the Vbelt, and the fiber direction is a group of longitudinal reinforcement.
Figure 5 shows the finite element model with a loading speed of 10 mm/min. Because of the symmetry of the model, we take half of the model to reduce the amount of calculation. As can be seen from Figure 5, when the lower portion of the sample is in a stretched state, the portion in contact with the indenter is in a compressed state. Figure 6 shows a graph of the mass force of the rigid indenter and the corresponding displacement. It can be seen that when the bending stiffness is zero, the curve obtained by the simulation is greatly different from the experimental curve, which also indicates the conventional continuous medium. The construction model is only fit for simple stretching and compression and is not fit for bending conditions. When is 5 MPa and 22.3 MPa, respectively, it can be seen that the curve obtained by the simulation and the curve of the experiment are gradually consistent. This proves the correctness of the hyperelastic constitutive model with bending stiffness.
5. Conclusions
Both the theoretical and simulation results show that the constitutive model based on the conventional continuum mechanics has achieved great success, but it is only fit for the tensile or compression deformation and is not suitable for the bending deformation.
A set of parameter identification schemes was proposed. The parameters of matrix materials were fitted with the data by transverse uniaxial tensile test, and the parameters of the matrix were determined. Then, the uniaxial test along the cord direction is used to fit the other parameters of the constitutive equation, and finally all the material parameters of the constitutive equation except the bending stiffness are determined.
Based on the invariant theory of Spencer, we introduce the fiber vector and propose a new hyperelastic constitutive model with the interaction of bending stiffness for anisotropic fiberreinforced belt. The results of experiments are in good agreement with the numerical simulations.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was supported by the National Natural Science Foundation of China (11172171), the Science and Technology Project of Guizhou Province (no. 20161063), the Joint Fund of Guizhou Province (no. 20167105), and the Scientific Research Startup Project of HighLevel Talent for GIT (XJGC20190928).
References
 J. Adkins and R. Rivlin, “Large elastic deformations of isotropic materials X. Reinforcement by inextensible cords,” Philosophical Transactions of the Royal Society of London Series A, Mathematical and Physical Sciences, vol. 248, no. 944, pp. 201–223, 1955. View at: Publisher Site  Google Scholar
 A. J. M. Spencer, Deformations of FiberReinforced Materials, Clarendon Press, Oxford, UK, 1972.
 A. J. M. Spencer and K. P. Soldatos, “Finite deformations of fibrereinforced elastic solids with fibre bending stiffness,” International Journal of Nonlinear Mechanics, vol. 42, no. 2, pp. 355–368, 2007. View at: Publisher Site  Google Scholar
 T. H. Nam and B. Marvalová, “Deformation analysis of inflated cylindrical membrane of composite with rubber matrix reinforced by cords,” in Proceedings of the XXI International Congress of Theoretical and Applied Mechanics ICTAM 2004, Warsaw, Poland, August 2004. View at: Google Scholar
 X. Q. Peng, Z. Y. Guo, and B. Moran, “An anisotropic hyperelastic constitutive model with fibermatrix shear interaction for the human annulus fibrosus,” Journal of Applied Mechanics, vol. 73, no. 5, pp. 815–824, 2006. View at: Publisher Site  Google Scholar
 G. A. Holzapfel, T. C. Gasser, and R. W. Ogden, “A new constitutive framework for arterial wall mechanics and a comparative study of material models,” Journal of Elasticity and the Physical Science of Solids, vol. 61, pp. 1–48, 2000. View at: Google Scholar
 D. J. Steigmann, “Continuum theory for elastic sheets formed by inextensible crossed elasticae,” International Journal of NonLinear Mechanics, vol. 106, pp. 324–329, 2018. View at: Publisher Site  Google Scholar
 T. Asmanoglo and A. Menzel, “A multifield finite element approach for the modelling of fiberreinforced composites with fiberbending stiffness,” Computer Methods in Applied Mechanics and Engineering, vol. 317, pp. 1037–1067, 2017. View at: Publisher Site  Google Scholar
 T. Asmanoglo and A. Menzel, “Fiberreinforced composites with fiberbending stiffness under azimuthal shear–comparison of simulation results with analytical solutions,” International Journal of NonLinear Mechanics, vol. 91, pp. 128–139, 2017. View at: Publisher Site  Google Scholar
 K. P. Soldatos, “On the characterisation of polar fibrous composites when fibers resist bending,” International Journal of Solids and Structures, vol. 148149, pp. 35–43, 2018. View at: Publisher Site  Google Scholar
 S. E. S. Bolouri, C. I. Kim, and S. Yang, “Linear theory for the mechanics of thirdgradient continua reinforced with fibers resistance to flexure,” Mathematics and Mechanics of Solids, vol. 25, no. 4, pp. 937–960, 2019. View at: Publisher Site  Google Scholar
 Q.S. Zheng, “Theory of representations for tensor functions—a unified invariant approach to constitutive equations,” Applied Mechanics Reviews, vol. 47, no. 11, pp. 545–587, 1994. View at: Publisher Site  Google Scholar
 S. Ishikawa, A. Tokuda, and H. Kotera, “Numerical simulation for fiber reinforced rubber,” Journal of Computational Science and Technology, vol. 2, no. 4, pp. 587–596, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Shulei Sun and Wenguo Chen. 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.