WJM  Vol.11 No.7 , July 2021
Implementation of Particle Swarm Optimization Algorithm in Matlab Code for Hyperelastic Characterization
Abstract: The purpose of this paper is to demonstrate the applicability of Particle Swarm Optimization algorithm to determine material parameters in incompressible isotropic elastic strain-energy functions using combined tension and torsion loading. Simulation of rubber behavior was conducted from the governing equations of the deformation of a cylinder composed of isotropic hyperelastic incompressible materials. Four different forms of strain-energy function were considered based respectively on polynomial, exponential and logarithmic terms to reproduce load force (N) and torque (M) trends using natural rubber experimental data. After highlighting the minimization of the objective function generated in the fitting process, the study revealed that a particle swarm optimization algorithm could be successfully used to identify the best material parameters and characterize the behavior of rubber-like hyperelastic materials.

1. Introduction

The study of the mechanical behavior of materials is of great importance to various industrial applications such as automation, aeronautics and biomechanics. The wide use of hyperelastic materials in industrial applications accounts for the high research interest expressed both in their characterization through experimental techniques and in the study of accurate constitutive models capable of describing their particular behaviour. Three criteria are necessary to determine the real behavior of a material: the choice of an adapted mechanical test, the choice of a constitutive law and the strategy of optimization of the material’s parameters. The optimization of hyperelastic parameters is an essential step in the designing of industrial parts. Several numerical methods suitable for the collection of experimental data are available in the literature notably Levenberg-Marquardt algorithm [1] [2] [3] [4], Beda-Chevalier method [5] [6] and the most recent methods such as Genetic Algorithm [7] [8] [9] (GA) and Particle Swarm Optimization method (PSO) [10] [11]. It has been shown that the PSO has better performance in finding the correct optimal solution than GA from a stochastic point of view [12]. The authors compared the performance of PSO and GA in terms of their accuracy, their robustness vis-à-vis population size and algorithm coefficients. Their study also revealed that not only does PSO have a higher probability to find the correct optimal solution than GA but is preferable to GA when time is a limiting factor. Without making a comparative study of optimization methods, this work seeks to demonstrate the applicability of the powerful PSO method to optimize the parameters of hyperelastic constitutive models.

The choice of an appropriate test is also essential to predict the mechanical behavior of hyperelastic materials. Materials’ parameters were thus optimized using a combined tension-torsion test. Simultaneous tension-torsion experiments are a relevant alternative to consider multiple loading conditions [13] and guarantee the model’s reliability to reproduce complex loading conditions. Several studies show that this type of test is well adapted to complex problems [14] [15] [16]. Lectez et al. [17] for example, use the tension-torsion test to characterize the behavior of elastometric engine mounts. Humphrey et al. [18] studied the finite extension and torsion of a nonlinear homogenous transversely isotropic incompressible solid cylinder. They proposed an analytical solution of the problem specifically for experimental correlation to identify a constitutive model for papillary muscles. These characterizations of rubbery materials on combined tension-torsional tests are reported in the literature. Emphasis is on the problem of torsion superimposed on the axial extension of a circular cylinder composed of an incompressible isotropic hyperelastic material [19] - [26].

This paper aims to implement PSO algorithm in Matlab code in order to determine the material parameters of the hyperelastic incompressible isotropic models using combined tension-torsion test data. The adopted method consists of reproducing trends of load force (N) and torque (M) of selected models using a solid circular cylinder’s large deformations. To investigate PSO algorithm’s appropriateness, the identification step of the material parameters is presented based on the correlation between the torsional couple (rather than the conventional methods of stress correlation or strain energy density correlation) rebuilt from test data and the theory. The identified parameters have been validated from the load force curves. Based on the combined tension-torsion test, the predictions of the selected models concur with experimental data.

2. Fundamentals

Homogeneous large strain elastic materials such as elastomers and biological tissues are conveniently modeled by hyperelastic constitutive equations. Practically, there exists a strain energy density function W, such that the stress tensor its derivative with respect to the strain [27]. Let us denote B the left Cauchy-Green strain tensor defined as F F T , F being the deformation gradient tensor and T the transposition. Based on the assumption that elastomers are both isotropic and incompressible, this strain-energy function only depends on the two first invariants of the left B, i.e.

W = W ( I 1 , I 2 ) , (1)


I 1 = tr B and I 2 = I 1 2 tr ( B 2 ) . (2)

The corresponding Cauchy stress tensor σ is expressed as follows:

σ = p I + 2 W I 1 B 2 W I 2 B 1 , (3)

where p is a hydrostatic pressure term associated with the incompressibility constraint and I is the identity tensor. The equations describing the tension-torsion deformation of an isochoric cylinder were derived by Rivlin and Saunders [19]. We consider a deformed circular cylinder of radius R, length l, subject to a combination of a tension load N and a torsion torque M. The deformation consists in an axial displacement u in the direction ez (Figure 1). Let us assume that ( r , θ , z ) is the undeformed cylindrical coordinate and ( R , Θ , Z ) , the coordinate in the deformed configuration. This can be written thus:

Figure 1. Circular cylinder subject to combined tension and torsion.

R = r λ , Θ = θ + Ψ λ z and Z = λ z , (4)

where the constant of proportionality Ψ is termed the amount of torsion, θ the angle through which one end of the cylinder is rotated with respect to the other and λ the global extension is given by

θ = Ψ l = Ψ λ L . (5)

Considering a cylinder subject to both extension and torsion, the deformation gradient tensor is:

F = [ λ 1 / 2 0 0 0 λ 1 / 2 λ 1 / 2 Ψ r 0 0 λ ] , (6)

and the left Cauchy-Green tensor B is:

B = [ λ 1 0 0 0 λ 1 + λ Ψ 2 r 2 λ 3 / 2 Ψ r 0 λ 3 / 2 Ψ r λ 2 ] . (7)

Thus, the main strain invariants are reduced to

I 1 = λ Ψ 2 r 2 + 2 λ + λ 2 and I 2 = Ψ 2 r 2 + 1 λ 2 + 2 λ . (8)

The resultant applied load force and torque necessary to maintain the deformation are given by:

M = 4 π Ψ 0 a ( W I 1 + 1 λ W I 2 ) r 3 d r , (9)


N = 4 π ( λ λ 2 ) 0 a ( W I 1 + 1 λ W I 2 ) r d r 2 π Ψ 2 0 a ( W I 1 + 2 λ W I 2 ) r 3 d r . (10)

3. Hyperelastic Constitutive Models

A broad collection of models are available in the literature. In this study four phenomenological models were used models:

Mooney-Rivlin model [28]

The mooney-rivlin model is rather simple and uses a reduced amount of parameters. It deals with moderate nonlinear and simple material behavior. Strain-energy density function based on two invariants of the right Cauchy-Green deformation tensor is defined as:

W = C 10 ( I 1 3 ) + C 01 ( I 2 3 ) , (11)

where C 10 and C 01 are the material parameters.

Yeoh model [29]

The Yeoh model is a reduced third-order polynomial based only on the first invariant of the right Cauchy-Green tensor. This model deals with high non-linear material behavior using only three parameters.

W = C 10 ( I 1 3 ) + C 20 ( I 1 3 ) 2 + C 30 ( I 1 3 ) 3 , (12)

where C 10 , C 20 , and C 30 are the material parameters.

Gent and Thomas model [30]

Considering the general form proposed by Rivlin and Saunders [17], Gent and Thomas proposed the following empirical strain energy function which involves only two material parameters and a logarithmic function of the second invariant:

W = K 1 ( I 1 3 ) + K 2 log I 2 3 , (13)

where K 1 and K 2 are the material parameters.

Fung model [31]

Yuan-Cheng Fung, the Founder of Modern Biomechanics, formulated the famous exponential strain constitutive equation for preconditioned soft tissues. The expression of this strain energy density function is based only on the first invariant of the right Cauchy-Green tensor

W = C ( exp β ( I 1 3 ) 1 ) , (14)

where C and β are the material constants.

4. PSO Algorithm for Hyperelastic Characterization

The ability of the different material models to reproduce any deformation curve depends on the definition of the values of the parameters associated to each model. Furthermore, the choice of these parameters depends on how a particular experiment is conducted and how the experimental data are processed [32] [33] [34] [35]. For an elastic material subject to large strains, the usual approach is to directly approximate the constants that appear in the mathematical model by employing numerical optimization techniques to minimize the residual between the stress-strain relation (for case of classical tests), the torsional couple or load force (for case of tension-torsion test, see for example [36] [37] ) and the experimental data. Then, identification consists in relating the theoretical solution M m o d e l resulting from a model with the experimental curves represented by a couple of the point ( M e x p , Ψ e x p ) , with i = 1 , 2 , 3 , , n the number of the experimental points. To quantify the gap between the two responses, the objective function f o b j is generally defined as:

f o b j = i = 1 n ( ( Ψ i ) e x p ( Ψ i ) m o d e l ) 2 . (15)

A PSO algorithm will be used to minimize the objective function and identify the material parameters to obtain a proper fitness of the curves. The load-force curves simulation response is fitted to validate the specified parameters by the PSO method employed.

The PSO algorithm was formulated by Edward and Kennedy in 1995 [38]. This algorithm is a stochastic optimization technique that is well adapted to the optimization of nonlinear functions in multidimensional space. It has been applied to several real world problems [39]. It is a class of swarm intelligence algorithms that attempts to reproduce the social behavior of a generic population. PSO population is called a swarm and each individual in the population of PSO is called a particle. All of the particles have fitness values evaluated by the fitness function to be optimized. They have velocities that direct the ?ying of the particles. The system is initialized with a population of random solutions and searches for optima by updating generations. To seek the optimal solution, each particle moves in the direction to its previously best ( P b e s t ) position and the global best ( G b e s t ) position in the swarm. The speed and position of each particle change according to the following equality (see [40] ):

V i , j k + 1 = w V i , j k + c 1 r 1 ( P b e s t i , j k X i , j k ) + c 2 r 2 ( G b e s t i , j k X i , j k ) , (16)

X i , j k + 1 = X i , j k + V i , j k + 1 . (17)

The index i varies from 1 to n, whereas the index j varies from 1 to D. n is population size and D is component of each individual of population. P b e s t i , j k represents personal best jth component of ith individual, whereas G b e s t j k represents jth component of the best individual of population up to iteration k. w, the inertia weight used to balance the global exploration and local exploitation, r 1 and r 2 are uniformly distributed random variables within range [ 0,1 ] , and c 1 and c 2 are positive constant parameters called acceleration coefficients.

Application of the PSO to the Mooney-Rivlin model.

Set Parameters

F: objective fonction.

n = 100 : population size.

w min = 0.4 : minimal inertia weigh.

w max = 0.9 : maximal inertia weigh.

c 1 = 0.5 : acceleration coefficient of particule.

c 2 = 1 : acceleration coefficient of essaim.

V i t max = 1.5 : maximal velocity.

V i t min = 1.5 : minimal velocity.

I max = 100 : maximum number of iteration.

C 01 min = 0 : lower bound of variable C 01 .

C 01 max = 5 : upper bound of variable C 01 .

C 10 min = 0 : lower bound of variable C 10 .

C 10 max = 5 : upper bound of variable C 10 .

Initialize population of particles having positions and velocities

C 01 = C 01 min + ( C 01 max C 01 min ) r a n d ( 1 , n )

V C 01 = V i t min C 01 r a n d ( 1 , n )

C 10 = C 10 min + ( C 10 max C 10 min ) r a n d ( 1 , n )

V C 10 = V i t min C 01 r a n d ( 1 , n )


Set iteration k = 1 .

Calculate fitness of particles F ( C p ) < F ( C g ) ; and find the index of the best particle p = g .

Update velocity and position of particles

V C 01 i , j k + 1 = w V C 01 i , j k + c 1 r a n d ( 1 ) ( C p 01 b o n i , j k C 01 i , j k ) + c 2 r a n d ( 1 ) ( C g 01 b o n i , j k C 01 i , j k )

C 01 i , j k + 1 = C 01 i , j k + V C 01 i , j k + 1

V C 10 i , j k + 1 = w V C 10 i , j k + c 1 r a n d ( 1 ) ( C p 10 b o n i , j k C 10 i , j k ) + c 2 r a n d ( 1 ) ( C g 10 b o n i , j k C 10 i , j k )

C 10 i , j k + 1 = C 10 i , j k + V C 10 i , j k + 1

with w = w min + ( w max w min ) i t e r / I max .

Evaluate fitness

F ( C p b o n e n d ) = F ( C g b o n e n d )

If k < I max then k = k + 1 and goto step update velocity and position of particles else goto.

Print optimal solutions

C p b o n e n d = C g b o n e n d


The methodology presented in this section is summarized in the flowchart presented in Figure 2.

5. Results and Discussion

5.1. Pure Torsion Application

The application of combined finite axial stretch and finite torsion to a solid right cylinder of length l and radius is described in Rivlin’s works (1949). The cylinder is first subjected to a simple extension in which its length increases and its radius decreases. Then, it is subjected to a simple torsion.

In the pure torsion case, the cylinder is maintained unextended, so that λ = 1 , Equations (9) and (10) become:

M = 4 π Ψ 0 a ( W I 1 + W I 2 ) r 3 d r . (18)


N = 2 π Ψ 2 0 a ( W I 1 + 2 W I 2 ) r 3 d r (19)

The analytical expressions of the torsional couple (M) and load force (N) calculated for each model in this case are summarized in Table 1.

Figure 2. PSO Flowchart.

Table 1. Analytical solutions of load force N and torque M calculated on pure-torsion test.

5.1.1. Optimization of Material Parameters

The optimization of the material parameters already described above is carried out in the case of pure torsion. The experimental data are those of Rivlin [41] and the analytical expressions of the torque are those given in Table 1. An optimization particle swarm algorithm was used to fit material constants’ best values by minimizing the objective function.

Fitness evaluation

Finally, the algorithm must evaluate the fitness level between the experimental and theoretical curves. In this case, it is clear that the fitness function must be related to the error between the curve created with a set of optimized parameters and the data obtained with the real test. The process of fitness evaluation is aimed at determining the margin of error. So, when the error of the curve is near to zero, fitness is near to 100.

Relative error is defined by

| Δ M | M = | M e x p M m o d e l | M e x p (20)

where M e x p is the experimental torque and M m o d e l , the theoretical predicted response of the model. The calculation of the relative errors generated in the fitting process presented in Figure 3 reveals that the relative errors are 0.012, 0.015, 0.016 and 0.016 for Mooney-Rivlin, Fung, Gent-Thomas and Yeoh models, respectively.

Both generated and approximated curves are presented in Figure 4 to show case the results obtained. From Figure 4, it can be observed that all the curves almost coincide along all the strain domain. It can thus be concluded that if

Figure 3. Plot of the relative error.

Figure 4. Best fit material parameters optimization: torsional torque curves on pure torsion test.

Table 2. Results of the identified parameters.

enough amount of data is used to fit the curve, the results will converge and correctly match the input data. The quest in this work for a proper match between the experimental and theoretical torsional torque curve of each model, led to the determination of the values of the material parameters which are presented in Table 2.

5.1.2. Verification of the Identified Material Parameters

To verify the accuracy of the optimized parameters, the identified constants in Table 1 were employed to predict the corresponding responses of the load force curve (N). The resulted calculated were compared with the data and plotted as presented in Figure 5. Based on the identified parameters, Figure 5 shows agreement between the experimental and calculated results curves. This suggests that the proposed identification method by using PSO is reliable and efficient.

5.2. Characterization of Hyperelastic Models by Combined Tension-Torsion Test

Prediction of Tension-Torsion Behavior

In the case of combined tension and torsion, the uniaxial extension is kept constant during torsion. The analytical expressions of the torque M of the deformed cylinder are presented with each model in Table 3 and the simulation results of

Figure 5. Verification of identified material parameters: load force curves on pure torsion test.

Table 3. Analytical solutions of torque M calculated on combined tension and torsion loading.

the torque in Figures 6-9.

To investigate the prediction power of the mechanical behaviour of hyperelastic models using combined tension and torsion test, the material parameters were determined for the first time by the PSO algorithm. This mechanical study shows that exponential and logarithmic strain energy functions do not give a good prediction of the tension-torsion test. The Fung (Figure 6) and Gent-Thomas (Figure 7)

Figure 6. Prediction of tension-torsion curve by Fung model.

Figure 7. Prediction of tension-torsion curve by Gent-Thomas model.

Figure 8. Prediction of tension-torsion curve by Yeoh model.

Figure 9. Prediction of tension-torsion curve by Mooney-Rivlin model.

models failed at small torsion angles. It is at intermediate angles that these two models reproduce the tensile-torsional test curve. The polynomial models are good models for the prediction of the tensile-torsional stress. Yeoh’s model produces the tensile-torsional test up to a value of Ψ = 0.7 rad / cm , then changes concavity (see Figure 8) due to the absence of I 2 in this model. Among the selected models, the Mooney-Rivlin model, the most accurate model to reproduce the shape of the experimental curve, as shown in Figure 9. The PSO method used in this work thus allowed for a qualitative mechanical characterization of the hyperelastic models.

6. Conclusions

The aim of this work was to identify the parameters of hyperelastic incompressible isotropic materials models using PSO algorithm through the combined tension and torsion loading test. Four classical hyperelastic models were used based on the simulation of their torque and force curves to determine that is the best theoretical model capable of reproducing the mechanical behavior for rubber-like materials under tension-torsion loading. An optimization particle swarm algorithm was used to fit material constants’ best values on the pure torsion test. These parameters were then validated and used to simulate the tension-torsion behavior. All in all, the results of the study revealed that:

1) The combined tension and torsion test’s mechanical behavior can be accurately adjusted by using the polynomial hyperelastic models with a good agreement with the Rivlin’s data for the torsionnal couple curve;

2) The particle swarm optimization algorithm can be successfully used for hyperelastic characterization and for the generation of best material parameters;

3) The models capable of better reproducing the mechanical behavior of rubber materials from the combined tension and torsion test depend on the two invariants of Cauchy-Green I 1 and I 2 .

Appendix A. Matlab Code: Identification of Material Parameters by PSO Algorithm

Cite this paper: Dya, T. , Blaise, B. , Betchewe, G. and Alidou, M. (2021) Implementation of Particle Swarm Optimization Algorithm in Matlab Code for Hyperelastic Characterization. World Journal of Mechanics, 11, 146-163. doi: 10.4236/wjm.2021.117011.

[1]   Martins, P.A.L.S., Natal Jorge, R.M. and Ferreira, A.J.M. (2006) A Comparative Study of Several Material Models for Prediction of Hyperelastic Properties: Application to Silicone-Rubber and Soft Tissues. Strain, 42, 135-147.

[2]   Wu, Y., Wang, H. and Li, A. (2016) Parameter Identification Methods for Hyperelastic and Hyper-Viscoelastic Models. Applied Sciences, 6, 386.

[3]   Nunes, L.C.S. (2011) Mechanical Characterization of Hyperelastic Polydimethylsiloxane by Simple Shear Test. Materials Science and Engineering, 528, 1799-1804.

[4]   Rachik, M., Schmidtt, F., Reuge, N., Le Maoult, Y. and Abbeé, F. (2001) Elastomer Biaxial Characterization Using Bubble Inflation Technique. II: Numerical Investigation of Some Constitutive Models. Polymer Engineering and Science, 41, 532-541.

[5]   Beda, T. (2005) Optimizing the Ogden Strain Energy Expression of Rubber Materials. Journal of Engineering Materials and Technology, 127, 351-353.

[6]   Beda, T. (2006) Combining Approach in Stages with Least Squares for Fits of Data in Hyperelasticity. Comptes Rendus Mecanique, 334, 628-633.

[7]   Fernández, J.R., López-Campos, J.A., Segade, A. and Vilán, J.A. (2018) A Genetic Algorithm for the Characterization of Hyperelastic Materials. Applied Mathematics and Computation, 329, 239-250.

[8]   Blaise, B.B., Betchewe, G. and Beda, T. (2019) Optimization of the Model of Ogden Energy by the Genetic Algorithm Method. Applied Rheology, 29, 21-29.

[9]   López-Campos, J.A., Segade, A., Casarejos, E., Fernández, J.R. and Días, G.R. (2019) Hyperelastic Characterization Oriented to Finite Element Applications Using Genetic Algorithms. Advances in Engineering Software, 133, 52-59.

[10]   Alam, M.N. (2016) Particle Swarm Optimization: Algorithm and Its Codes in MATLAB.

[11]   Ramzanpour, M., Hosseini-Farid, M., Ziejewski, M. and Karami, G. (2019) Particle Swarm Optimization Method for Hyperelastic Characterization of Soft Tissues. International Mechanical Engineering Congress and Exposition, Salt Lake City, 11-14 November 2019.

[12]   Duan, Y., Harley, R.G. and Habetler, T.G. (2009) Comparison of Particle Swarm Optimization and Genetic Algorithm in the Design of Permanent Magnet Motors. IEEE 6th International Power Electronics and Motion Control Conference, Wuhan, 17-20 May 2009, 822-825.

[13]   Rivlin, R.S. and Saunders, D.W. (1951) Large Elastic Deformations of Isotropic Materials VII. Experiments on the Deformation of Rubber. Philosophical Transactions of the Royal Society A, 243, 251-288.

[14]   Haupt, P. and Sedlan, K. (2001) Viscoplasticity of Elastomeric Materials. Experimental Facts and Constitutive Modeling. Archive of Applied Mechanics, 71, 89-109.

[15]   Mars, W.V. and Fatemi, A. (2004) A Novel Specimen for Investigating the Mechanical Behavior of Elastomers under Multiaxial Loading Conditions. Experimental Mechanics, 44, 136-146.

[16]   Suphadon, N., Thomas, A.G. and Busfield, J.J.C. (2009) Viscoelastic Behavior of Rubber under a Complex Loading. Polymer, 113, 693-699.

[17]   Lectez, A.S., Verron, E., Huneau, B., Béranger, A.S. and Le Brazidec, F. (2013) Characterization of Elastomers under Simultaneous Tension and Torsion for Application to Engine Mounts. In: Constitutive Models for Rubber VIII, CRC Press, Boca Raton, 585-590.

[18]   Humphrey, J.D. (2002) Cardiovascular Solid Mechanics. Springer-Verlag, New York, 758-766.

[19]   Rivlin, R.S. (1949) Large Elastic Deformations of Isotropic Materials VI. Further Results in the Theory of Torsion, Shear and Flexure. Philosophical Transactions of the Royal Society A, 242, 173-195.

[20]   Rivlin, R.S. and Saunders, D.W. (1951) Large Elastic Deformations of Isotropic Materials VII. Experiments on the Deformation of Rubber. Philosophical Transactions of the Royal Society A, 243, 251-288.

[21]   Sigaeva, T. and Czekanski, A. (2016) Coupling of Surface Effect and Hyperelasticity in Combined Tension and Torsion Deformations of a Circular Cylinder. International Journal of Solids and Structures, 85-86, 172-179.

[22]   Suphadon, N. and Busfield, J.J.C. (2009) Elastic Behaviour of Rubber Cylinders under Combined Torsion and Tension Loading. Plastics Rubber and Composites, 38, 337.

[23]   Kashaev, R.M. (2008) On Tension-Torsion Testing of Solid Cylindrical Specimens. Letters on Materials, 8, 346-352.

[24]   Kirkinis, E. and Ogden, R.W. (2002) On Extension and Torsion of a Compressible Elastic Circular Cylinder. Mathematics and Mechanics of Solids, 7, 373-392.

[25]   Michelle, S., Fatt, H. and Ouyang, X. (2007) Integral-Based Constitutive Equation for Rubber at High Strain Rates. International Journal of Solids and Structures, 44, 6491-6506.

[26]   Lectez, A.S., Verron, E. and Huneau, B. (2014) How to Identify a Hyperelastic Constitutive Equation for Rubber-Like Materials with Multiaxial Tension-Torsion Experiments. International Journal of Non-Linear Mechanics, 65, 260-270.

[27]   Treloar, L.R.G., Hopkins, H.G., Rivlin, R.S. and Ball, J.M. (1976) The Mechanics of Rubber Elasticity [and Discussions]. Proceedings of the Royal Society A: Mathematical Physical and Engineering Sciences, 351, 301-330.

[28]   Rivlin, R.S. (1948) Large Elastic Deformations of Isotropic Materials, I, II, III, Fundamental Concepts. Philosophical Transactions of the Royal Society A, 240, 459.

[29]   Yeoh, O.H. (1990) Characterization of Elastic Properties of Carbon-Black-Filled Rubber Vulcanizates. Rubber Chemistry and Technology, 63, 792-805.

[30]   Gent, A.N. and Thomas, A.G. (1958) Forms for the Stored (Strain) Energy Function for Vulcanized Rubber. Journal of Polymer Science, 28, 625-628.

[31]   Fung, Y.C.B. (1967) Elasticity of Soft Tissues in Sample Elongation. American Journal of Physiology, 213, 1532-1544.

[32]   Chatelin, S., Constantinesco, A. and Willinger, R. (2010) Fifty Years of Brain Tissue Mechanical Testing: From in Vitro to in Vivo Investigations. Biorheology, 47, 255-276.

[33]   Evans, S.L. (2017) How Can We Measure the Mechanical Properties of Soft Tissues? In: Avril, S. and Evans, S.L., Eds., Material Parameter Identification and Inverse Problems in Soft Tissue Biomechanics, Springer, London, 67-83.

[34]   Staber, B. and Guilleminot, J. (2017) Stochastic Hyperelastic Constitutive Laws and Identification Procedure for Soft Biological Tissues with Intrinsic Variability. Journal of the Mechanical Behavior of Biomedical Materials, 65, 743-752.

[35]   Wex, C., Arndt, S., Stoll, A., Bruns, C. and Kupriyanova, Y. (2015) Isotropic Incompressible Hyperelastic Models for Modelling the Mechanical Behaviour of Biological Tissues: A Review. Biomedical Engineering/Biomedizinische Technik, 60, 577-592.

[36]   Hartmann, S. (2001) Numerical Studies on the Identification of the Material Parameters of Rivlin’s Hyperelasticity Using Tension-Torsion Tests. Acta Mechanica, 148, 129-155.

[37]   Lectez, A.S. (2014) Comportement multiaxial de pièces élastomères précontraintes: Application aux suspensions moteur. Mécanique des matériaux [physics.class-ph] Ecole Centrale de Nantes, 151-170.

[38]   Kennedy, J. and Eberhart, R.C. (1995) Particle Swarm Optimization. Proceedings IEEE Conference on Neural Networks, Perth, 27 November-1 December 1995, 1942-1948.

[39]   Boehner, K., DePaula, R., Dourish, P. and Sengers, P. (2007) How Emotion Is Made and Measured. International Journal of Human-Computer Studies, 65, 275-291.

[40]   Shi, Y. and Eberhart, R. (1998) A Modified Particle Swarm Optimizer. IEEE International Conference on Evolutionary Computation, Anchorage, 4-9 May 1998, 69-73.

[41]   Rivlin, R.S. (1948) Large Elastic Deformations of Isotropic Materials III. Some Simple Problems in Cylindrical Polar Coordinates. Philosophical Transactions of the Royal Society A, 240, 509-525.