Parameter Identification Methods for Hyperelastic and Hyper-Viscoelastic Models

In this paper, the Ogden model is employed to characterize the hyperelastic properties of rubber, and on the basis of a pattern search algorithm and the Levenberg-Marquardt algorithm, a professional method that can realize the comprehensive fitting of the uniaxial tension, biaxial tension, planar tension, and simple shear experimental data of hyperelastic materials was developed. The experiment data from Treloar (1944) was fitted very well, and the determined parameters by using this method were proven correct and practical in the numerical verification in ANSYS. Then, the constitutive model of the hyper-viscoelastic materials combining the Ogden model with the generalized Maxwell model was explained in detail, and the parameter identification method was also proposed by using the pattern search method. Then, three groups of relaxation tests of uniaxial tension and four groups of simple shear tests with different loading velocities were conducted to obtain the corresponding virtual experiment data. After discussing the constraints and initial setting values for the undetermined parameters, these virtual data of different loading histories were respectively employed to identify the parameters in the hyper-elastic model, and the accuracy and the reliability of the estimated parameters were also verified in ANSYS.


Introduction
With the rapid development of modern industrialization, rubbers are one of the most remarkable materials having a wide range of applications in civil engineering, aerospace engineering, mechanical engineering, automotive engineering, etc. In order to meet various requirements of the industry, special fillers, like carbon black or silica, with different proportions are usually added during vulcanization for improving the strength and toughness properties, which in turn makes it difficult to accurately characterize the mechanical a properties of rubbers [1]. Rubbers usually present a number of interesting features like hyperelasticity and hyper-viscoelasticity, many studies and experiments have been performed to investigate these features. For the former one, the stress-strain relationship of the hyperelasticity can be illuminated by a strain energy function W, and some mathematical models have been developed for the function, among which the strain-invariant-based models mainly include the Mooney-Rivlin model [2], the Yeoh model [3], and the Gent model [4], the principal-stretch-based models mainly include the Ogden model [5]. All the models have been verified by experiment results, for example Treloar's experiment [6], and the Ogden model breaks through the limitation of the strain energy function being the even power of the stretches, which is capable of more accurately fitting the experimental data when rubbers undergo extremely large deformations. Comparisons of the performances of 20 hyperelastic models are presented in the article from Marckmann and Verron [7]. More detailed reviews on the hyperelasticity of rubber can be accessed in several literatures [8][9][10]. The identification of the parameters in theses hyperelastic models has also gained much attention. For example, Twizell and Ogden [11] employed the Levenberg-Marquardt optimization algorithm to determine the material constants in Ogden models, Saleeb et al. [12,13] developed a professional scheme named COMPARE to get the parameters by using uniaxial tension, biaxial tension, and planar tension data.
The constitutive relations for the viscoelasticity can be characterized by combining the elastic components with viscous components. In the past decades, numerous studies focused on investigating the viscoelasticity properties of rubber-like materials and many models explaining such behaviors were proposed. In the range of small strain, the basis of linear viscoelasticity theory was presented by Coleman et al. [14]. A relaxation-modulus constitutive function for the incompressible isotropic materials based on the results of uniaxial and equal biaxial stress relaxation tests was developed by McGuirt [15]. The phenomenological theory of linear viscoelasticity was proved powerful for modelling the static and dynamic responses of rubber when subjected to small strain [16,17]. For the finite deformation of viscoelastic materials, two important hypotheses were assumed. One is a split of the free energy of viscoelastic solid into an equilibrium and a non-equilibrium part, and the other one is a multiplicative decomposition of the deformation gradient tensor into elastic and inelastic parts [18]. In such a framework, Huber and Tsakmak is generalized the application range of the so-called three-parameter solids to the finite deformation [19]. Haupt et al. proposed a fairly simple identification method for the constants in a finite viscoelasticity model [20]. Yoshida et al. presented a constitutive model comprised by an elastoplastic body with a strain-dependent isotropic hardening law and a hyperelastic body with a damage model [21]. Amin et al. proposed an improved hyperelastic constitutive equation for rubbers in compression and shear regimes [1,22,23], the parameters in the equation were identified by experiment data and a nonlinear viscous coefficient was introduced to represent the rate-dependent properties. In recent years, a novel testing technique, named nanoidentation, to identify the material parameters for thin polymer layers was developed, a comprehensive review and the realization of the method for hyperelastic and viscoelastic models are presented by Chen et al. [24][25][26] In this paper, the Ogden model is taken to characterize the hyperelastic properties of rubber-like materials. A professional method based on the pattern search algorithm (PSA) [27,28] and the Levenberg-Marquardt algorithm (LMA) [29] was developed to fit the unixial tension, biaxial tension, planar tension, and the simple shear experimental data of hyperelastic materials. Then, combining the generalized Maxwell model and the Ogden model as the constitutive model, the method that can identify the parameters in the hyper-viscoelastic models by using experimental data with different loading histories was also developed [24,25] and three groups of virtual relaxation tests and four groups of simple shear tests were employed to verify the accuracy and reliability of the proposed method.

Parameter Identification Method for the Hyperelastic Model
The constitutive models representing the hyperelastic properties of rubbers mainly include the statistical models, the strain-invariant-based models, and the principal-stretch-based models, among which the Ogden model is able to fit the experiment data when rubbers undergo large deformations (λ < 7). In the following sections, the parameter identification method for the hyperelastic model and the numerical verification of these identified parameters are presented in detail.

Theories of the Ogden Model
According to the theories of the continuum mechanics, there exists a strain energy function W for the hyperelasticity properties of rubbers. The stresses can be obtained by the partial derivative of the variable W with respect to the strain.
where C is the right Cauchy-Green deformation gradient tensor, S is the second Piola-Kirchhoff stress. After transforming the second Piola-Kirchhoff stress to the Cauchy stress [30], we will get where F is the deformation gradient tensor, B is the left Cauchy-Green deformation gradient tensor, I 1 , I 2 , I 3 are the three strain invariants of B, and p is the undetermined hydrostatic pressure which can be decided by the underlying equilibrium and boundary conditions of the particular problem.
In view of the nearly incompressible properties of the volume of rubbers, the strain energy of the Ogden model can be defined with Equation (3) without regard for the contribution of the volume strain.
In Equation (3), n is the number of terms considered for the Ogden model, µ n , α n are the undetermined parameters of the model, λ i is the principal stretch.
In the uniaxial tension experiment, the deformation gradient tensor F and the left Cauchy-Green deformation gradient tensor B are described as Taking the boundary condition T 22 = T 33 = 0 into consideration, the corresponding nominal stress for uniaxial tension is For the biaxial tension experiment, the boundary conditions are T 22 = T 33 , T 11 = 0, λ 2 = λ 3 , and the nominal stress is The boundary condition of the planar tension experiment is T 33 = 0, λ 2 = 1, and the nominal stress is As for the simple shear experiment, the direction of the principal stretch is not consistent with the direction of the applied deformation, rather it involves a rotation of axes [1]. The left Cauchy-Green deformation tensor B is expressed as Appl. Sci. 2016, 6, 386 4 of 13 and the nominal shear stress is shown in Equation (9).
where the principal stretch λ 1 can be obtained by calculating the eigenvalues of the tensor B.

The Object Function for Parameter Identification
Previous studies have shown that the constitutive parameters of a hyperelastic model identified from a particular deformation mode may be invalid for other modes. For example, Charlton et al. claimed that the parameters identified from the uniaxial test data fail in predicting biaxial or planar tension responses [31]. So, it is of much significance to take the four types of experiment data listed above into consideration. This goal can be achieved by the least square procedure to make the deviations of the experiment data and the fitted data least. The object function is defined as where ω q is the weight of different experiment type, if there exists no experiment data for a specific type, we just need to set its ω q to be zero, if we want to emphasize the experiment data of a specific type, we can set its ω q larger than others'. n q represents the number of the experiment data for each type of experiment, ( f q ) i is the theoretical value which is calculated by Equations (5)- (7) and (9), (σ q ) i is the experiment data.

The Identified Parameters and the Numerical Verification
In this paper, the uniaxial tension, biaxial tension, and planar tension data from Treloar are taken as the experiment data [6], and the Ogden model is chosen as the hyperelastic model owing to its favorable performances under large deformations. In order to solve Equation (11), the PSA is employed as the first step for its remarkable ability in searching for the global optimum results. Then, the L-M algorithm is chosen to further optimize the results from the first step due to the probable poor efficiency of the PSA when approaching the optimum results. All the procedures are realized in MATLAB [32].
The number of terms considered for the Ogden model are respectively set as three and four, and seven groups of different initial values of undetermined parameters for each one are iterated in parallel. For the limited space and the similarity between the results of different groups, only the ones with the least iterations and its corresponding deviation S 1 are chosen and listed in Table 1. The estimated parameters presented by Treloar are presented as well [6]. The unit of µ i in Table 1 is Mpa, and the item 'PSA 169000+' means more than 169 thousand iterations in PSA, and the 'LMA 48' means 48 iterations in LMA after the previous PSA procedure. Moreover, the computation time for the whole procedure with an i7 CPU is about 6 h.
As the Table 1 shows, these three groups of parameters are all capable of fitting the experiment data very well, and the last group is the best owing to its least deviation 0.012., However, despite thousands of iterations of the PSA and LMA, the two sets of the identified parameters in Table 1 may not be the optimal results of Equation (12), because there may be local minimum points extremely close to the global minimum one, and these local optimal points have not been excluded by the PSA for their relatively excellent, not merely most excellent, competitions and they are identified by the LMA. Even so, these identified parameters, being optimal or near-optimal, are still very useful due to enough accuracy and controlled computation cost. Additionally, It should be mentioned that the PSA takes more than 200 thousand iterations to reduce the S 1 from 0.1 to 0.0598, while the LMA only takes 65 iterations to reduce the deviation from 0.06 to 0.012, which strongly reflects the inefficiency of the PSA when approaching the global optimal point, and taking the LMA as a succeeding procedure to further optimize results is practical and efficient. In the last iteration of the LMA, the sensitivity values for the items ∂S/∂u i and ∂S/∂α i are all more than 10 9 , which means accurate estimations of all the parameters. For the purpose of verifying the estimated parameters, the last group of the identified parameters is taken as the example to predict the response of rubbers in different types of experiments in ANSYS. The results of the numerical simulation, the fitted data obtained by the use of the proceeding algorithm, along with the experiment data from Treloar are all plotted in Figure 1. PSA takes more than 200 thousand iterations to reduce the S1 from 0.1 to 0.0598, while the LMA only takes 65 iterations to reduce the deviation from 0.06 to 0.012, which strongly reflects the inefficiency of the PSA when approaching the global optimal point, and taking the LMA as a succeeding procedure to further optimize results is practical and efficient. In the last iteration of the LMA, the sensitivity values for the items / u i S ∂ ∂ and / ∂ ∂α i S are all more than 10 9 , which means accurate estimations of all the parameters. For the purpose of verifying the estimated parameters, the last group of the identified parameters is taken as the example to predict the response of rubbers in different types of experiments in ANSYS. The results of the numerical simulation, the fitted data obtained by the use of the proceeding algorithm, along with the experiment data from Treloar are all plotted in Figure 1. As illustrated in the Figure 1, the results identified by the PSA and LMA fit the experiment data very well. Besides, in the numerical simulation tests, the results of uniaxial tension and planar tension coincide well with the experiment data (when the principal stretch is larger than 5, the FE analysis of uniaxial tension tests cannot converge due to some uncertain reasons). However, the results of biaxial tension tests differ a bit from the experiment data, the reason for the errors lies in that the boundary conditions of the biaxial tension experiment cannot be accurately satisfied when the principal stretch ranges from 1.5 to 1.9 in the simulation test.
The FE initial models of uniaxial tension, planar tension, and biaxial tension experiment and their deformed shapes with the stretch being 5, 6.14, and 5.04 are presented respectively in  As illustrated in the Figure 1, the results identified by the PSA and LMA fit the experiment data very well. Besides, in the numerical simulation tests, the results of uniaxial tension and planar tension coincide well with the experiment data (when the principal stretch is larger than 5, the FE analysis of uniaxial tension tests cannot converge due to some uncertain reasons). However, the results of biaxial tension tests differ a bit from the experiment data, the reason for the errors lies in that the boundary conditions of the biaxial tension experiment cannot be accurately satisfied when the principal stretch ranges from 1.5 to 1.9 in the simulation test.
The FE initial models of uniaxial tension, planar tension, and biaxial tension experiment and their deformed shapes with the stretch being 5, 6.14, and 5.04 are presented respectively in Figure 2.

Parameter Identification Method for the Hyper-Viscoelastic Model
In the previous section, the parameter identification method for the hyperelastic model was explained in detail, this model is suitable for the numerical simulation of natural rubbers with low damping. As for the high damping rubbers, they can be simulated by combining the Ogden model with the generalized Maxwell model. In the following, theories of the hyper-viscoelastic model and the parameter identification method for the model were expounded.

Theories of the Hyper-Viscoelastic Model
The integration of the relaxation modulus for viscoelasticity is expressed in Equation (12)

Parameter Identification Method for the Hyper-Viscoelastic Model
In the previous section, the parameter identification method for the hyperelastic model was explained in detail, this model is suitable for the numerical simulation of natural rubbers with low damping. As for the high damping rubbers, they can be simulated by combining the Ogden model with the generalized Maxwell model. In the following, theories of the hyper-viscoelastic model and the parameter identification method for the model were expounded.

Theories of the Hyper-Viscoelastic Model
The integration of the relaxation modulus for viscoelasticity is expressed in Equation (12)  where Y(t − ξ) is the relaxation modulus function for characterizing the viscoelastic properties, and dε(ξ) is the micro-strain of materials. Generally, when ξ < 0, the outer action is assumed to be 0, so Equation (12) will be changed to Equation (13) can be calculated by dividing the time t into n equal portions, then we will get When ξ ∈ (t i , t i+1 ), the Young modulus of materials is assumed to be a constant E i , and the Prony series is introduced to characterize the relaxation properties. Besides, the strain loading speed is also deemed as a constant λ i , then Equation (14) will become (15) where µ k is the weight of modulus corresponding to the relaxation time t k , m is the number of terms considered for the Prony series, µ 0 is the weight of the equilibrium modulus. It is obvious for the equation m ∑ k=0 µ k = 1 to be true, and the physical model of Equation (15) is presented in Figure 3. where ( ) − ξ Y t is the relaxation modulus function for characterizing the viscoelastic properties, and ( ) ε ξ d is the micro-strain of materials. Generally, when 0 ξ < , the outer action is assumed to be 0, so Equation (12) will be changed to Equation (13) can be calculated by dividing the time t into n equal portions, then we will get When 1 ( , ) + ξ∈ i i t t , the Young modulus of materials is assumed to be a constant i E , and the Prony series is introduced to characterize the relaxation properties. Besides, the strain loading speed is also deemed as a constant λ i , then Equation (14) will become where μ k is the weight of modulus corresponding to the relaxation time k t , m is the number of terms considered for the Prony series, 0 μ is the weight of the equilibrium modulus. It is obvious for the equation to be true, and the physical model of Equation (15) is presented in Figure 3. This physical model includes two main features, the first one is that it substitutes the stretchdependent modulus spring for the constant modulus spring of the generalized Maxwell model, which is based on an assumption that the modulus corresponded to the viscosity coefficient η k is obtained by multiplying the instantaneous modulus i E and the weight k u . The second one is that the instantaneous modulus i E can be decided by the hyperelasticity properties of materials, and the weight k u is the undetermined parameters which are time independent.
If we stipulate that In Figure 3, E i , µ 0 , µ k shares the same meaning with that in Equation (15), η k represents the viscosity coefficient and η k = µ k E i t k , t k is the relaxation time included in Equation (15).
This physical model includes two main features, the first one is that it substitutes the stretch-dependent modulus spring for the constant modulus spring of the generalized Maxwell model, which is based on an assumption that the modulus corresponded to the viscosity coefficient η k is obtained by multiplying the instantaneous modulus E i and the weight u k . The second one is that the instantaneous modulus E i can be decided by the hyperelasticity properties of materials, and the weight u k is the undetermined parameters which are time independent.
If we stipulate that t k = 10 k−1 t 1 , which is adopted by the FE software LS-DYNA [33], then, σ(t)can be expressed as µ k e −(t−ξ)/(10 k−1 t 1 ) )dξ (16) As mentioned before, E i represents the instantaneous modulus of materials within the time range from t i to t i+1 . It can be solved by the following procedures. Firstly, the fast stretch-stress experiment data is obtained by conducting experiment at sufficiently high speed to exclude the viscous effect. Then they will be fitted by the Ogden model to decide the undetermined parameters, which has been explained in the previous section. After acquiring the stretch-stress curve based on the Ogden model, E i will be solved by the derivation of the stress σ with respect to the principal stretch λ, namely, ∂σ ∂λ t=t i , or taking the secant modulus as a substitution for the tangent modulus, namely,

The Parameter Identification Algorithm for the Hyper-Viscoelastic Model
According to Equation (16), the number of the undetermined parameters is m + 1, including u k and t 1 . Similarly to Equation (11), the object function for the parameter identification of hyper-viscoelastic model can be defined as follows, where p represents the number of the tests with different loading histories, n j is the number of the experiment data for the jth group of experiment, T i is the experiment data, f i is the theoretical value which is calculated by Equation (16). In order to get the global optimum results for the undetermined parameters in Equation (17), the global optimal searching algorithm, namely pattern search algorithm, is also employed in this section.

Different Loading Histories and Virtual Experiment Data
Similarly, with the identification in the Section 2.2, it is necessary to incorporate different loading histories to guarantee the accuracy of the identified multi-parameters in Equation (17). Hence, the relaxation tests of uniaxial tension and the simple shear tests with different loading velocities are taken as examples in this study. Due to the lack of appropriate data for these experiments, we followed the method applied in Chen et al. [24][25][26] and got the corresponding virtual experiment data by using ANSYS, in which we assume the hyperelastic parameters are the same with the last group in Table 1, and the viscoelastic part is characterized by the following Prony series, P(ξ) = 0.4e −(t−ξ)/10 + 0.1e −(t−ξ)/100 (18) In the virtual experiment, the finite element (FE) model for the relaxation tests is identical with that in Figure 2a, and the model for the simple shear tests is illustrated in Figure 4. All the material constants are determined by the above discussions. experiment data is obtained by conducting experiment at sufficiently high speed to exclude the viscous effect. Then they will be fitted by the Ogden model to decide the undetermined parameters, which has been explained in the previous section. After acquiring the stretch-stress curve based on the Ogden model, i E will be solved by the derivation of the stress σ with respect to the principal stretch λ, namely, or taking the secant modulus as a substitution for the tangent modulus, namely,

The Parameter Identification Algorithm for the Hyper-Viscoelastic Model
According to Equation (16), the number of the undetermined parameters is m + 1, including k u and 1 t . Similarly to Equation (11), the object function for the parameter identification of hyperviscoelastic model can be defined as follows,  (17) where p represents the number of the tests with different loading histories, j n is the number of the experiment data for the jth group of experiment, i T is the experiment data, i f is the theoretical value which is calculated by Equation (16). In order to get the global optimum results for the undetermined parameters in Equation (17), the global optimal searching algorithm, namely pattern search algorithm, is also employed in this section.

Different Loading Histories and Virtual Experiment Data
Similarly, with the identification in the Section 2.2, it is necessary to incorporate different loading histories to guarantee the accuracy of the identified multi-parameters in Equation (17). Hence, the relaxation tests of uniaxial tension and the simple shear tests with different loading velocities are taken as examples in this study. Due to the lack of appropriate data for these experiments, we followed the method applied in Chen et al. [24][25][26] and got the corresponding virtual experiment data by using ANSYS, in which we assume the hyperelastic parameters are the same with the last group in Table 1, and the viscoelastic part is characterized by the following Prony series, In the virtual experiment, the finite element (FE) model for the relaxation tests is identical with that in Figure 2a, and the model for the simple shear tests is illustrated in Figure 4. All the material constants are determined by the above discussions.  Then, three groups of data for relaxation tests and four groups of simple shear experiment data with different loading velocities are obtained. As shown in Figure 5, the maximum stretch is 1.3 in the relaxation tests with the loading time being respectively 3 s, 6 s, and 30 s, and the hold time being Then, three groups of data for relaxation tests and four groups of simple shear experiment data with different loading velocities are obtained. As shown in Figure 5, the maximum stretch is 1.3 in the relaxation tests with the loading time being respectively 3 s, 6 s, and 30 s, and the hold time being 150 s. For the simple shear tests, the maximum shear strain is 1.0 with the loading velocities being 0.5/s, 0.1/s, 0.01/s, and 0.001/s.

Constraints and Initial Values of Parameters
Before presenting results, it is necessary to clarify the constraints as well as the initial setting values for the undetermined parameters in the identified procedure. As for the weight k u , as well as the summation of all k u except 0 u , they should be always within the range from 0 to 1. While for in which the 0 σ is the instantaneous stress and the ∞ σ is the equilibrium stress.
The creep experiment of the material characterized by Equation (18) and the last group of parameters in Table 1 is conducted in ANSYS, and the main relaxation time is 12.54 s.
As for the number of the terms considered for the Prony series, theoretically, the experiment data in Figure 5 will be more accurately fitted with a higher number, however, this may lead to uncontrolled computation cost and sometimes unreasonable estimated parameters. In this study, the number to be decided is suggested to make the largest relaxation time in Equation (16) the same order of magnitude with the largest experiment time. For example, the largest experiment time in Figure 5 is 1000 s, so, the number of the terms considered for the Prony series is 3 since the corresponding largest relaxation time now is 1254 s.
In view of all the restrictions presented above, as well as the consideration of testing the reliability of the pattern search algorithm, four groups of initial values for the undetermined parameters are presented in Table 2. Besides, taking the third group as an example, the linear constraint equations for the parameters during iterations are shown in Equation (20).

Constraints and Initial Values of Parameters
Before presenting results, it is necessary to clarify the constraints as well as the initial setting values for the undetermined parameters in the identified procedure. As for the weight u k , as well as the summation of all u k except u 0 , they should be always within the range from 0 to 1. While for the relaxation time t 1 , it should be positive and a suggested initial value is the main relaxation time of materials, which can be determined by the point corresponding to the relaxed stress σ 0.368 in creep tests, here the σ 0.368 is defined as in which the σ 0 is the instantaneous stress and the σ ∞ is the equilibrium stress. The creep experiment of the material characterized by Equation (18) and the last group of parameters in Table 1 is conducted in ANSYS, and the main relaxation time is 12.54 s.
As for the number of the terms considered for the Prony series, theoretically, the experiment data in Figure 5 will be more accurately fitted with a higher number, however, this may lead to uncontrolled computation cost and sometimes unreasonable estimated parameters. In this study, the number to be decided is suggested to make the largest relaxation time in Equation (16) the same order of magnitude with the largest experiment time. For example, the largest experiment time in Figure 5 is 1000 s, so, the number of the terms considered for the Prony series is 3 since the corresponding largest relaxation time now is 1254 s.
In view of all the restrictions presented above, as well as the consideration of testing the reliability of the pattern search algorithm, four groups of initial values for the undetermined parameters are presented in Table 2. Besides, taking the third group as an example, the linear constraint equations for the parameters during iterations are shown in Equation (20).

Uniaxial Relaxation Tests
The identified parameters by each group of relaxation test are presented in Table 3. The original values for the undetermined parameters t 1 , u 1 , u 2 , and u 3 are respectively 10, 0.4, 0.1, and 0.

T load (s) = 3 s T load (s) = 6 s
Group As shown in Table 3, all the groups with different loading histories and initial values are able to identify the parameters in the hyper-viscoelastic models. The identified weight in each group is nearly identical with its original one, while the weights u 2 , and u 3 in some groups, are a bit poorly determined. For the identified relaxation time t 1 , it becomes more accurate with the increasing of the loading time T load, and in the group of the T load being 30 s, the identified relaxation time is about 3.6% lower than its true value. Subsquently, the deviation S 2 of each group decreases as the T load increases with the least one being 2.87 × 10 −5 . Besides, for the results with the same loading time, the identified parameters are little affected by the initial setting value and the number of terms for the Ogden model, which indicates a strong reliability of the PSA.

Simple Shear Tests
For the four groups of simple shear experiment data in Figure 5, they are not employed separately like the relaxation tests, instead all data with different shear velocities are combined together to determine the global optimal results. Similarly, the PSA is utilized and the identified parameters are shown in Table 4. It can be seen from Table 4 that the parameters in different groups are almost the same with each other, and the weights u 1 and u 2 are almost accurately identified with the relaxation time a bit deviated from its true value. However, an abnormal negative value −0.001, extremely close to 0, of the weight u 3 is observed even though we have set constraints for all the weights as shown in Equation (20), this unusual phenomenon may be an indication of the weight u 3 redundant for the problem.

Verification of the Identified Results
In order to verify the accuracy of the estimated parameters, one out of the four groups of the identified parameters in Tables 3 and 4 are employed in ANSYS to predict the corresponding responses of the hyper-viscoelastic materials. Comparisons of the results calculated by ANSYS and the data in Figure 5 are plotted in Figure 6.
In Figure 6, the solid line represents the original data and the scatter points represent the results from ANSYS by using the identified parameters. It is obvious that errors in every experiment are very small, which means the proposed analytical model for hyper-viscoelastic materials and the identification method by using PSA are reliable and efficient. It can be seen from Table 4 that the parameters in different groups are almost the same with each other, and the weights u 1 and u 2 are almost accurately identified with the relaxation time a bit deviated from its true value. However, an abnormal negative value −0.001, extremely close to 0, of the weight u3 is observed even though we have set constraints for all the weights as shown in Equation (20), this unusual phenomenon may be an indication of the weight u3 redundant for the problem.

Verification of the Identified Results
In order to verify the accuracy of the estimated parameters, one out of the four groups of the identified parameters in tables 3 and 4 are employed in ANSYS to predict the corresponding responses of the hyper-viscoelastic materials. Comparisons of the results calculated by ANSYS and the data in Figure 5 are plotted in Figure 6.
In Figure 6, the solid line represents the original data and the scatter points represent the results from ANSYS by using the identified parameters. It is obvious that errors in every experiment are very small, which means the proposed analytical model for hyper-viscoelastic materials and the identification method by using PSA are reliable and efficient.

Conclusions
In this study, we focused on the parameter identification of the hyperelastic and hyperviscoelastic materials, and the following conclusions have been made: • Based on the pattern search algorithm and the Levenberg Marquardt algorithm, the comprehensive fitting of the multi-type experiment data of hyperelastic models, including uniaxial tension, planar tension, biaxial tension, and simple shear experiment have been realized.

•
According to the fitting results of the experiment data from Treloar (1944) and the verification results of the numerical simulation, the method of combining the pattern search algorithm and Levenberg Marquardt algorithm proposed in this paper for identifying parameters in the Ogden model is practical and efficient.

Conclusions
In this study, we focused on the parameter identification of the hyperelastic and hyper-viscoelastic materials, and the following conclusions have been made:

•
Based on the pattern search algorithm and the Levenberg Marquardt algorithm, the comprehensive fitting of the multi-type experiment data of hyperelastic models, including uniaxial tension, planar tension, biaxial tension, and simple shear experiment have been realized. • According to the fitting results of the experiment data from Treloar (1944) and the verification results of the numerical simulation, the method of combining the pattern search algorithm and Levenberg Marquardt algorithm proposed in this paper for identifying parameters in the Ogden model is practical and efficient.

•
The parameter identification method for the hyper-viscoelastic materials which is characterized by combing the Ogden model with the generalized Maxwell model was developed, and the physical model of the hyper-viscoelastic model was also presented. • Three groups of relaxation tests and four groups of simple shear tests were conducted in ANSYS and the corresponding experiment data were obtained. After employing the data to determine the parameters, it can be concluded that the proposed analytical model of the hyper-viscoelastic model and the identification method by using a pattern search algorithm is reliable since all the experiment data are fitted very well with extremely low deviations, and the identified results are little affected by the initial values as well as the term considered for Maxwell models.