An Anisotropic Hyperelastic Constitutive Model with Bending Stiffness Interaction forCord-RubberComposites: Comparison of Simulation Results with Experimental Data

Based on the invariant theory of continuummechanics by Spencer, the strain energy depends on deformation, fiber direction, and the gradients of the fiber direction in the deformed configuration. -e 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 reinforcedmaterial has tensile deformation and compression deformation.-e bending stiffness should be taken into consideration, and it is further verified by the bending simulation.


Introduction
e theory of finite deformation of fiber-reinforced 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 fiber-reinforced composites. is 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 fiber-reinforced materials [4] and intervertebral disc tissue in biology, and it achieved success. e 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. e 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 two-dimensional elastic surface theory [7]. e 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 fiber-reinforced composites with fiber bending stiffness under azimuthal shear [9]. Soldatos researched the invention of methods appropriate for characterization of fiber-reinforced 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].
ere is a premise that the fiber is infinitely soft under uniaxial compressive condition. is 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. e comparison between the simulation results and the bending experimental data of the rubber-cord composite proves the correctness of the model.

Constitutive Model of Conventional Anisotropic
Hyperelastic Fiber-Reinforced Composites. According to Spencer's theory of anisotropic fiber-reinforced composites [3], the Helmholtz strain energy function (per unit reference volume) W depends on the deformation and fiber direction a 0 : where the tensor F is deformation gradient t, the tensor C � F T · F is the right Cauchy-Green deformation, and the tensor A (a 0 ⊗ a 0 ) is of order two. e Helmholtz strain energy function per unit volume W can also be written in the form of invariants for C and A: W � W I 1 (C), I 2 (C), I 3 (C), I 4 (C, A), I 5 (C, A) , I 1 (C) � tr(C), Since the deformation gradient F can be decomposed into volume-preserving deformation and volumetric deformation, the strain energy function W 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. J − �� I 3 � det(F) is the deformation volume ratio, and I 1 , I 2 , I 4 , and I 5 are the corresponding partial invariants: (4) In order to facilitate parameter fitting, the strain energy can also be given by where According to the above analysis, for anisotropic fiber-reinforced composites that consider compressibility, a common form of strain energy per unit volume has the form [8-10] By means of chain rule, the second-order 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 form where where p � (dW vol /dJ), e Cauchy stress tensor is defined to be 2 Mathematical Problems in Engineering

Constitutive Model of Fiber-Reinforced Composites with
Bending Stiffness. e 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. e 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 G � zb/zX. By analogy with relation (12), we may postulate the free energy per unit reference energy: e free energy must be unchanged if the fiber-reinforced composite undergoes a rotation described by the proper orthogonal tensor Q. We obtain the following expressions: where Q is arbitrary orthogonal tensor. e tensors F and G can be written in the form of a scalar product, respectively: Hence, we can get Γ � Λ T C − 1 Λ using the Cayley-Hamilton formulation [11]: Let the tensor Γ be represented by the invariants of Λ. Finally, the function W that can be represented as Λ, C, A is expressed as According to the formula in Zheng [12], the free energy W consists of 33 invariants (Table 1) Λ, C, A 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 A(zb/zX) � A · G � κ 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 F, κ can be written in the form of product, respectively: is W can be expressed as a function of C, A, K, κ 2 : 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).

Mathematical Problems in Engineering
In summary, the total strain energy with the bending stiffness is in the form of where the free energy W G is the energy contributed to the invariants I 6 ∼ I 11 which is associated with the bending stiffness.

Determination of the Strain Energy Function and Corresponding Parameters of V-Belt
e V-belt shown in Figure 1 is a V-belt 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. e belt has anisotropic characteristics. Ishikawa et al. carried out a series of experiments and simulation studies [13]. e following is mainly based on the experimental data for the parameter fitting and simulation of the constitutive model.

Strain Energy and Parameters of Matrix Rubber.
Matrix rubber is generally regarded as a hyperelastic material. In order to obtain accurate material parameters, the V-belt is subjected to transverse uniaxial tensile test. e stress-strain 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. e typical Mooney-Rivlin model is used to compare the mechanical properties of the matrix material. e fitting curve is shown in Figure 2. e Mooney-Rivlin strain energy form of the rubber with compressibility is as follows: e 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,

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 force-displacement curve and finally express the tensile Mathematical Problems in Engineering strain energy as a function I 4 − 1. e contribution of energy I 5 is related to the fiber-matrix interaction. It is written according to formula (6) as follows: Data fitting was performed using MATLAB programming. e fitted curve is shown in Figure 3. e 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].

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.

Numerical Simulation
Combining the material parameters obtained by the previous parameter, the corresponding finite element program is compiled to simulate the bending of the rubber-cord composite, i.e., V-belt as shown in Figure 4. e cuboid sample of size 60 mm × 25 mm × 7 mm is cut from the V-belt, 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. e construction model is only fit for simple stretching and compression and is not fit for bending conditions. When C 6 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. is proves the correctness of the hyperelastic constitutive model with bending stiffness.

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. e parameters of matrix materials were fitted with the data by transverse uniaxial tensile test, and the parameters of the matrix were determined. en, 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 fiber-reinforced belt. e results of experiments are in good agreement with the numerical simulations.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.