Placement 2025 Scholarship: Your Future Starts Here | 6 Guaranteed Job Interviews | 1̶0̶0̶ 29 seats only. Apply Now

00D 00H 00M 00S

Menu

Executive Programs

Workshops

Projects

Blogs

Careers

Placements

Student Reviews


For Business


More

Academic Training

Informative Articles

Find Jobs

We are Hiring!


All Courses

Choose a category

Mechanical

Electrical

Civil

Computer Science

Electronics

Offline Program

All Courses

All Courses

logo

CHOOSE A CATEGORY

Mechanical

Electrical

Civil

Computer Science

Electronics

Offline Program

Top Job Leading Courses

Automotive

CFD

FEA

Design

MBD

Med Tech

Courses by Software

Design

Solver

Automation

Vehicle Dynamics

CFD Solver

Preprocessor

Courses by Semester

First Year

Second Year

Third Year

Fourth Year

Courses by Domain

Automotive

CFD

Design

FEA

Tool-focused Courses

Design

Solver

Automation

Preprocessor

CFD Solver

Vehicle Dynamics

Machine learning

Machine Learning and AI

POPULAR COURSES

coursePost Graduate Program in Hybrid Electric Vehicle Design and Analysis
coursePost Graduate Program in Computational Fluid Dynamics
coursePost Graduate Program in CAD
coursePost Graduate Program in CAE
coursePost Graduate Program in Manufacturing Design
coursePost Graduate Program in Computational Design and Pre-processing
coursePost Graduate Program in Complete Passenger Car Design & Product Development
Executive Programs
Workshops
For Business

Success Stories

Placements

Student Reviews

More

Projects

Blogs

Academic Training

Find Jobs

Informative Articles

We're Hiring!

phone+91 9342691281Log in
  1. Home/
  2. Anish Augustine/
  3. Week - 10 Hyperelastic Material Models

Week - 10 Hyperelastic Material Models

Hyperelastic Material Modelling using LS-DYNA AIM: To calculate the Mooney Rivlin and Ogden material constants and compare the both using stress-strain data from a Dogbone specimen tensile test with 100 percent strain. The given material data is the engineering stress-strain in MPa/(mm/mm). The comparison should be shown…

    • Anish Augustine

      updated on 14 Jan 2021

    Hyperelastic Material Modelling using LS-DYNA

    AIM:

    1. To calculate the Mooney Rivlin and Ogden material constants and compare the 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.

    Note: The unit system used is g-mm-ms, N, MPa, N-mm.

    PROCEDURE:

    Analysis setup for Material model verification:

    1. Experimental Data:

    1

    Fig.1 Experimental data of Engg. stress and Engg. strain of hyperelastic material.

    2

    Fig.2 Experimental data plotted with Engg. stress vs Engg. strain.

    The given data of Engg. stress and Engg. strain of hyperelastic material is plotted in excel is as shown in fig.2.

    2. Hyperelastic Material Model:

    Hyperelastic constitutive models describe the elastic property of a material from the strain energy. The specific form of the strain energy controls the elastic material properties of the model. There are many different kinds of functions. They all try to follow the stress stretch curve for different loading cases. The simplest ones are built as a polynomial.

    The given material data refers to rubber which is generally considered to be fully incompressible since the bulk modulus greatly exceeds the shear modulus in magnitude. To model the rubber as an unconstrained material a hydrostatic work term, 𝑊H(𝐽), is included in the strain energy functional which is function of the relative volume, 𝐽. They are usually written in the following form,

    W(J1,J2,J)=∑np,q=0Cpq(J1−3)p(J2−3)q+WH(J)W(J1,J2,J)=∑p,q=0nCpq(J1-3)p(J2-3)q+WH(J)

    i) The Mooney-Rivlin rubber model is obtained by specifying 𝑛 = 1.

    W=C10(I1−3)+C01(I2−3)W=C10(I1-3)+C01(I2-3)

    where, C10 and C01 are Mooney-Rivlin constants. MAT_77_H material card is used to find these constants and obtain a curve matching with experimental data.

    ii) The Ogden rubber model is obtained from the equation,

    W(λ1,λ2,λ3)=∑np=1μpαp(λαp1+λαp2+λαp3−3)W(λ1,λ2,λ3)=∑p=1nμpαp(λ1αp+λ2αp+λ3αp-3)

    Where, μpμp and αpαp are Ogden constants. MAT_77_O material card is used to find these constants and obtain a curve matching with experimental data.

    3. Part definition:

    To verify and validate the given material data with material models of Mooney-Rivlin and Ogden a dogbone specimen is taken for the analysis.

    3

    Fig.3 Material card for hyperelastic rubber model with N=1

    4

    Fig.4 Material card for Ogden rubber model with N=1

    To find the Mooney-Rivlin and Ogden constants the value of N is set to 1. The value of density and Poisson’s ratio are taken as 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

    The value of SGL, SW, and ST are set to unity (1.0), then the curve LCID1 is engineering stress versus engineering strain as shown in fig. 5.

    5

    Fig.5 LCID1 curve for given data of Engg. stress vs Engg. strain.

    6

    Fig.6 Section property of dogbone specimen.

    The section property of dogbone specimen is assigned as shell element with 2 mm thickness and ELFORM=2.

    7

    Fig.7 Part definition of dogbone specimen.

    4. Boundary Condition:

    8

    Fig.8 The nodes at the fixed end constrained in X and Z direction.

    The nodes at the fixed end is constrained in X and Z direction only and Y direction is not constrained because of lateral expansion during tensile test.

    9

    Fig.9 The nodes at the middle constrained in Y direction.

    The nodes at the middle is constrained in Y direction since the neutral axis passes through the middle of the specimen in X direction.

    10.1

    Fig.10.1 Boundary prescribed motion in the pulling end of specimen.

    The nodes of the pulling end are assigned with a boundary prescribed motion in X direction using displacement load curve LCID as shown in fig. 10.2.

    10.2

    Fig.10.2 Displacement load curve.

    5. Control Functions:

    11

    Fig.11 Control functions for implicit analysis.

    The tensile test of given dogbone specimen is considered as quasi-static analysis, hence implicit analysis is carried out with necessary control implicit cards as shown in fig.11. The keyword file is saved and made to run in LS-DYNA program manager to get d3hsp file.

    Note: The keyword file for each material model is set up separately and made to run individually.

    II. Analysis setup for Material model validation:

    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.

    12

    Fig.12 d3hsp file of Mooney-Rivlin material model.

    From the output files of material verification, d3hsp file is opened in notepad++ as shown in fig. 12. The Mooney-Rivlin constants obtained are

    C10 = c1 = 0.1768E+00

    C01 = c2 = 0.1474E+00

    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 as shown in fig. 13. The data obtained from d3hsp file is extension (stretch) and True stress which is converted to Engg. strain and Engg, stress to verify with experimental data.

    13

    Fig.13 Material card for hyperelastic rubber model with N=0

    14

    Fig.14 d3hsp file of Ogden material model.

    From the d3hsp output file of Ogden material model the Ogden constants obtained are

    μ1=9.1076025585232E−01μ1=9.1076025585232E-01

    α1=1.3427238694706E+00α1=1.3427238694706E+00

    These values are inputted to the material card MAT_77_O as shown in fig. 15. The data obtained from d3hsp file is Engg. Stress and stretch ratio which is converted to Engg. strain to verify with experimental data.

    15

    Fig.15 Material card for Ogden rubber model with N=0

    16

    Fig.16 Keywords used for the simulation of material validation.

    The keywords used for material model validation is same as that of material model verification except the changes to material card(N=0) and addition of database extent binary to compute elastic strain is as shown in fig.16. The keyword file is saved and made to run in LS-DYNA program manager to get requested output file.

    Note: The simulation is carried out for different polynomial values of N=1,2,3 for each material model so as to get better adjustment to the experimental results.

    RESULTS AND DISCUSSIONS:

    1. The animation of Von-Mises stress contour.

    i) Mooney-Rivlin material model:

    17

    ii) Ogden material model:

    18

    The maximum value of v-m stress developed for Mooney-Rivlin material model(N=1) is 1.66801 MPa, whereas for Ogden material model(N=1) maximum v-m stress is equal to 1.65691 MPa. Hence, the maximum v-m stress value developed is equal for both the material model.

    2. The animation of Effective Plastic Strain contour.

    i) Mooney-Rivlin material model:

    19

    ii) Ogden material model:

    20

    The maximum Effective plastic strain developed for both material models is equal to 3.8.

    3. Material model verification:

    i) Mooney-Rivlin material model:

    17

    Fig.17 Plot of simulation results of constants verification with experimental data.

    ii) Ogden material model:

    18

    Fig.18 Plot of simulation results of constants verification with experimental data.

    From the data obtained from d3hsp file for different polynomial values i.e, (N=1,2,3) curves are plotted with Engg. stress vs Engg. strain in excel with comparison to Experimental data is as shown in fig.17 and Fig.18. It is observed from the graph, that the curves for different polynomial are superimposing and is better adjusted with the experimental data. The Ogden material model gives close fit compared to Mooney-Rivlin material model.

    4. Material model validation:

    i) Mooney-Rivlin material model:

    19

    Fig.19 Plot of simulation results of constants validation with experimental data.

    ii) Ogden material model:

    20

    Fig.20 Plot of simulation results of constants validation with experimental data.

    For validation, the plot of true stress vs true strain for a middle element of dogbone specimen is saved in .csv file format. The file is opened in excel and the values of true stress and strain are converted to Engg. stress and Engg. strain to plot the curves as shown in fig.19 and Fig. 20. Form the graph, there is a slight deviation of simulation results compared to experimental data for both the material models because the material is highly non-linear and generic values of density and Poisson’s ratio was considered for simulation.

    5. Comparison of Mooney-Rivlin and Ogden material model:

    21

    Fig.21 Plot of validation results of Mooney-Rivlin and Ogden constants with experimental data.

    From the graph, both material models are behaving similar for the given material data. The Ogden material model has polynomials ranging from 1 to 8 and is easy to use and gives better adjustment to experimental data.

    CONCLUSION:

    1. The Mooney-Rivlin and Ogden material constants are obtained from d3hsp file by inputting the value of polynomial N=1 in the material card MAT_77_H and MAT_77_O
    2. 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.
    3. The simulation results of tensile test for both material models with experimental data showed a slight deviation due to high non-linear behaviour of the given hyperelastic material and the generic values of density and Poisson’s ratio.
    4. From the simulation results, both hyperelastic constitutive material models were found to be valid to describe the behaviour of the given material data.
    5. 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 computional time.
    6. Learned to verify and validate hyperelastic material model for the given material data.

    Google Drive Link:Hyperelastic Material modelling using LS-DYNA

    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.

    Please  login to add a comment

    Other comments...

    No comments yet!
    Be the first to add a comment

    Read more Projects by Anish Augustine (38)

    Week 11 Car Crash simulation

    Objective:

    CAR CRASH SIMULATION USING ANSYS WORKBENCH OBJECTIVE 1. To simulate car crash for different thickness of car body, Case-1: Thickness=0.3 mm. Case-2: Thickness=0.7 mm. Case-3: Thickness=1.5 mm. 2. To find out Total deformation and Equivalent stress developed in car body for each case and compare the results. 1. THEORY 1.1…

    calendar

    14 Jul 2021 09:52 AM IST

    • CAE
    Read more

    Week 10 Bullet penetrating a Bucket Challenge

    Objective:

    SIMULATION OF BULLET PENETRATING INTO A BUCKET USING ANSYS WORKBENCH OBJECTIVE To simulate bullet penetrating into a bucket for different cases of bucket material, Case-1: Aluminium Alloy NL Case-2: Copper Alloy NL Case-3: Stainless Steel NL To find out Total deformation and Equivalent stress developed in bucket for…

    calendar

    19 Jun 2021 08:51 AM IST

      Read more

      Week 9 Tension and Torsion test challenge

      Objective:

      SIMULATION OF TENSION AND TORSION TEST ON A SPECIMEN USING ANSYS WORKBENCH OBJECTIVE To perform the tension and torsion test simulation on the specimen by following the necessary boundary conditions, For the tension test, one end of the specimen has to be displaced to 18mm while keeping the other end fixed. For the torsion…

      calendar

      11 Jun 2021 11:10 AM IST

        Read more

        Week 9 Machining with Planer Challenge

        Objective:

        EXPLICIT DYNAMIC ANALYSIS OF MACHINING WITH PLANER USING ANSYS WORKBENCH OBJECTIVE To perform explicit dynamic analysis of machining with planer for the following two different cases of cutting velocity, Case-1: Cutting velocity=20000 mm/s Case-2: Cutting velocity=15000 mm/s To find out Directional Deformation, Equivalent…

        calendar

        06 Jun 2021 03:39 AM IST

          Read more

          Schedule a counselling session

          Please enter your name
          Please enter a valid email
          Please enter a valid number

          Related Courses

          coursecard

          Design loads considered on bridges

          Recently launched

          10 Hours of Content

          coursecard

          Design of Steel Superstructure in Bridges

          Recently launched

          16 Hours of Content

          coursecard

          Design for Manufacturability (DFM)

          Recently launched

          11 Hours of Content

          coursecard

          CATIA for Medical Product Design

          Recently launched

          5 Hours of Content

          coursecardcoursetype

          Accelerated Career Program in Embedded Systems (On-Campus) Courseware Partner: IT-ITes SSC nasscom

          Recently launched

          0 Hours of Content

          Schedule a counselling session

          Please enter your name
          Please enter a valid email
          Please enter a valid number

          logo

          Skill-Lync offers industry relevant advanced engineering courses for engineering students by partnering with industry experts.

          https://d27yxarlh48w6q.cloudfront.net/web/v1/images/facebook.svghttps://d27yxarlh48w6q.cloudfront.net/web/v1/images/insta.svghttps://d27yxarlh48w6q.cloudfront.net/web/v1/images/twitter.svghttps://d27yxarlh48w6q.cloudfront.net/web/v1/images/youtube.svghttps://d27yxarlh48w6q.cloudfront.net/web/v1/images/linkedin.svg

          Our Company

          News & EventsBlogCareersGrievance RedressalSkill-Lync ReviewsTermsPrivacy PolicyBecome an Affiliate
          map
          EpowerX Learning Technologies Pvt Ltd.
          4th Floor, BLOCK-B, Velachery - Tambaram Main Rd, Ram Nagar South, Madipakkam, Chennai, Tamil Nadu 600042.
          mail
          info@skill-lync.com
          mail
          ITgrievance@skill-lync.com

          Top Individual Courses

          Computational Combustion Using Python and CanteraIntroduction to Physical Modeling using SimscapeIntroduction to Structural Analysis using ANSYS WorkbenchIntroduction to Structural Analysis using ANSYS Workbench

          Top PG Programs

          Post Graduate Program in Hybrid Electric Vehicle Design and AnalysisPost Graduate Program in Computational Fluid DynamicsPost Graduate Program in CADPost Graduate Program in Electric Vehicle Design & Development

          Skill-Lync Plus

          Executive Program in Electric Vehicle Embedded SoftwareExecutive Program in Electric Vehicle DesignExecutive Program in Cybersecurity

          Trending Blogs

          Heat Transfer Principles in Energy-Efficient Refrigerators and Air Conditioners Advanced Modeling and Result Visualization in Simscape Exploring Simulink and Library Browser in Simscape Advanced Simulink Tools and Libraries in SimscapeExploring Simulink Basics in Simscape

          © 2025 Skill-Lync Inc. All Rights Reserved.

                      Do You Want To Showcase Your Technical Skills?
                      Sign-Up for our projects.