All Courses
All Courses
Courses by Software
Courses by Semester
Courses by Domain
Tool-focused Courses
Machine learning
POPULAR COURSES
Success Stories
Aim: 1) To calculate the Mooney Rivlin and Ogden material constants and compare both using stress-strain data from a Dogbone specimen tensile test with 100 percent strain. 2) The given material data is the engineering stress-strain in MPa/(mm/mm). 3) The comparison should be shown from the d3hsp file and using simulation.…
Kumarasamy N
updated on 01 Jul 2021
Aim:
1) To calculate the Mooney Rivlin and Ogden material constants and compare both using stress-strain data from a Dogbone specimen tensile test with 100 percent strain.
2) The given material data is the engineering stress-strain in MPa/(mm/mm).
3) The comparison should be shown from the d3hsp file and using simulation.
Theory:
Hyper-elastic Material:
Hyperelastic material models are regularly used to represent the large deformation behaviour of materials with FEA. They are commonly used to model mechanical behaviours of unfilled/filled elastomers. In addition to elastomers, hyperelastic material models are also used to approximate the material behaviour of biological tissues, polymeric foams, etc. Linearly elastic materials are described through two material constants (like Young’s modulus and Poisson ratio). In contrast, hyperelastic materials are described through a strain-energy density function. The strain-energy density can be used to derive a nonlinear constitutive model (i.e., stresses as a function of large strain deformation measures like deformation) which are derived by researches and some acceptable models we are using popularly such as Neo-Hookean, Mooney-Rivlin and Signorini models, Ogden, Yeoh, Arruda-Boyce.
Mooney Rivlin:
The strain energy density function for an incompressible Mooney–Rivlin material is defined as
Where Cij">Cij are empirically determined material constants, I1 and I2 are the first and the second invariant. By varying the Index values of i and j we can change the constitutive equations and use them appropriately for our component. Similarly, we will be having multiple coefficients C. This changing of the index will be basically changing the order of the polynomial of the energy density equation. We need to find out the coefficients of these material models for a different order of the polynomial and extract the best fit curve data from the D3HSp file and compare it with the material curve.
Procedure:
1. Experimental Data:
1) The given data of Engg stress and Engg strain of hyperelastic material is shown in the image given below are plotted in excel.
2. Part Definition:
1) To verify and validate the given material data with material models of Mooney-Rivlin and Ogden a dogbone specimen is taken for the analysis.
2) To find the Mooney-Rivlin and Ogden constants the value of N is set to 1-3. The value of density and Poisson’s ratio is taken as the generic value for the rubber material.
For N>0, data from a uniaxial test are used.
SGL Specimen gauge length
SW Specimen width
ST Specimen thickness
3) The value of SGL, SW, and ST is set to unity (1.0), then the curve LCID1 is defined using *DEFINE_CURVE engineering stress versus engineering strain.
4) The section property of the dogbone specimen is assigned as a shell element with 2 mm thickness and ELFORM=2.
5) The Material and Section are assigned to the dogbone specimen part.
3. Boundary Condition:
1) The nodes at the fixed end are constrained in the X and Z direction only and the Y direction is not constrained because of lateral expansion during the tensile test.
2) The nodes at the middle is constrained in the Y direction since the neutral axis passes through the middle of the specimen in the X-direction.
3) The nodes of the pulling end are assigned with a boundary prescribed motion in X direction using displacement load curve LCID.
4. Control Cards:
1) *CONTROL_IMPLICIT_GENERAL
Activate implicit analysis and define associated control parameters. This keyword is required for all implicit analyses
2) *CONTROL_IMPLICIT_SOLUTION
The linear equation solver performs the CPU-intensive stiffness matrix inversion
3) *CONTROL_IMPLICIT_AUTO card has been defined to adjusts the time step size
4) *CONTROL_TERMINATION:
This card is used to mention the termination time of the simulation.
5. Analysis setup for Material model validation:
1) The d3hsp file for each material model is opened in notepad++ and the material constants of Mooney-Rivlin and Ogden as well as final fit data of stretch and Engg. stress is obtained for each case.
2) From the output files of material verification, the d3hsp file is opened in notepad++. The Mooney-Rivlin constants obtained are
C10 = c1 = 0.4419E-01
C01 = c2 = 0.3686E-01
For material model validation, the value of N is changed to 0 in the material card. These values are inputted to the material card MAT_77_H. The data obtained from the d3hsp file is an extension (stretch) and True stress which is converted to Engg. strain and Engg, stress to verify with experimental data.
3) From the d3hsp output file of the Ogden material model the Ogden constants obtained are
`mu1`=2.2769006165419E-01
`alpha1`=1.3427238801312E+00
4) These values are inputted to the material card MAT_77_O. The data obtained from the d3hsp file is Engg Stress and stretch ratio which is converted to Engg strain to verify with experimental data.
5) The keywords used for material model validation is the same as that of material model verification except for the changes to the material card(N=0) and addition of database extent binary to compute elastic strain. The keyword file is saved and made to run in the LS-DYNA program manager to get the requested output file.
Results and Discussion:
Animation Contour for Von mises Stress:
The maximum value of v-m stress developed for the Mooney-Rivlin material model(N=0) is 0.4341 MPa, whereas for the Ogden material model(N=0) maximum v-m stress is equal to 0.4308 MPa. Hence, the maximum v-m stress value developed is equal for both the material model.
X-Stress Plot:
The maximum X-stress developed for Mooney-Rivlin material model(N=0) is 0.4346 MPa.
The maximum X-stress developed for Mooney-Rivlin material model(N=1,2,3) is 0.4342,0.4295,0.4297 MPa.
maximum X-stress developed for Ogden material model(N=0) is 0.4312 MPa.
maximum X-stress developed for Ogden material model(N=1,2,3) is 0.4301,0.4302,0.4303 MPa.
X-Strain Plot:
The maximum Lower Ipt X-strain developed for Mooney-Rivlin material model(N=0) is 7.347e-01 and maximum Lower Ipt X-strain developed for Ogden material model(N=0) is 7.357e-01.
From the data obtained from the d3hsp file for different polynomial values i.e, (N=1,2,3) curves are plotted with Engg. stress vs Engg. strain in excel is as shown in the below images. It is observed from the graph, that the curves for different polynomial are superimposing and are better adjusted with the experimental data. The Ogden material model gives a close fit compared to the Mooney-Rivlin material model.
Comparison of Mooney-Rivlin and Ogden material model: From this data, we can easily found that Young's Modulus value for hyperelastic and Ogden material. The Ogden material model gives a close fit compared to the Mooney-Rivlin material model.
CONCLUSION: - The Mooney-Rivlin and Ogden material constants are obtained from the d3hsp file by inputting the value of polynomial N=1 in the material card MAT_77_H and MAT_77_O
- The stretch stress curve data obtained from d3hsp for both hyperelastic constitutive models were verified with experimental material data and found to be a close fit.
- The simulation results of the tensile test for both material models with experimental data showed a slight deviation due to the high non-linear behavior of the given hyperelastic material and the generic values of density and Poisson’s ratio.
- From the simulation results, both hyperelastic constitutive material models were found to be valid to describe the behavior of the given material data.
- The Ogden material model has more polynomials ranging from 1 to 8 and is easy to use and gives better adjustment to experimental data at the cost of computational time.
- Learned to verify and validate the hyperelastic material model for the given material data.
Leave a comment
Thanks for choosing to leave a comment. Please keep in mind that all the comments are moderated as per our comment policy, and your email will not be published for privacy reasons. Please leave a personal & meaningful conversation.
Other comments...
Bird Strike Simulation
AIM: To simulate a non-linear transient dynamic problem of a Bird Strike with a Jet Engine Blade and evaluate the results. OBJECTIVES: 1. The given file contains all the 4 parts which as to separate and save part wise. 2. The suitable material is to be created and saved in a subfolder. 3. The…
12 Jul 2021 07:17 AM IST
week-11 Joint creation and Demonstration
AIM: To demonstrate spherical, revolute, cylindrical, and translational joints between two rigid bodies and two deformable bodies. OBJECTIVE: 1. Demonstrate the different types of joints like Spherical, Revolute, Cylindrical & Translational joints. 2. For all the joints simulate with rigid & elastic…
06 Jul 2021 04:17 PM IST
Week - 10 Hyperelastic Material Models
Aim: 1) To calculate the Mooney Rivlin and Ogden material constants and compare both using stress-strain data from a Dogbone specimen tensile test with 100 percent strain. 2) The given material data is the engineering stress-strain in MPa/(mm/mm). 3) The comparison should be shown from the d3hsp file and using simulation.…
01 Jul 2021 12:24 PM IST
Week - 9 Material Modeling from Raw Data
AIM : 1) Using the diagram of the true stress-strain curve given for graphite iron casting, creating a material model for the Dogbone specimen. 2) Validating the material model by performing a tensile test on a Dogbone specimen. OBJECTIVE : By Using the diagram of the true stress-strain curve given for graphite…
28 Jun 2021 01:13 PM IST
Related Courses
Skill-Lync offers industry relevant advanced engineering courses for engineering students by partnering with industry experts.
© 2025 Skill-Lync Inc. All Rights Reserved.