Exercise 2: Tensile Test of Elastomer

In this exercise, you will perform a uni-axial tensile test on a rubber strip (2mmx5mmx50mm).

Before you begin, in the menu bar click File > New > Session to clear your session from Exercise 1: Model Cantilever Beam Bending.

Hyper-elastic material constants have been sourced from reference[3].

Hyper-elastic materials are large strain materials when compared to metals. In case of hyper-elastic materials the non-linear relation between stress and strain is derived from a strain energy density function. Currently MotionSolve supports three hyper-elastic material models: Neo-Hookean, Mooney–Rivlin, and Yeoh.

Add a New Material Property

In this step, you will add a new material property.

  1. From the menu bar, click Model > Materials.
  2. In the Material Properties dialog, for Elasticity type choose Hyper Elastic.
  3. Click the Add button.


    Figure 1. Adding new HyperElastic material
  4. In the Add a Material Property dialog, specify a Label of Yeoh Model and a Variable name as propmat_yeoh.
  5. For Source for values, select Rubber (Yeoh).


    Figure 2. Selecting source values for new material
  6. Click OK to add the material property.
  7. In the Material Properties dialog, specify the material constant values.
    1. Element shear modulus (c10): 0.545235
    2. Element shear modulus (c20): 0.0610498
    3. Element shear modulus (c30): −0.000802537
    4. Poisson ratio (): 0.48
    5. Density ( ): 1.1e-6
    6. Elastic strain limit (εL): 2.0


    Figure 3. Specifying material constant values
  8. Click Close.

Model the Rubber Strip

In this step, you will model the rubber strip for the tensile test.

  1. Create two point entities for the rubber strip length profile using the details in Table 1:
    Table 1.
    S.No Label Variable name x y z
    1 Rubber End 1 p_rub_end1 0.0 0.0 0.0
    2 Rubber End 2 p_rub_end2 50.0 0.0 0.0
  2. Use the (Create Points Along A Vector) macro to create 9 intermediate points between Rubber End 1 and Rubber end 2.


    Figure 4. Creating intermediate points using “Create Points along a Vector” macro
  3. On the Model-Reference toolbar, click the (Body) icon.
  4. In the dialog, add a new NLFE body with a Label Rubber Strip and a Variable name nlfeb_rubber_yeoh.
  5. Configure the Properties tab.
    1. Specify the Type as Beam.
    2. Specify the Cross-section as Bar.
    3. For dim1, enter 2.
    4. For dim2, enter 25.
    5. For the Material Property, select the Yeoh Model that you created previously.


    Figure 5. Specifying beam properties
  6. Configure the Connectivity tab.
    1. Append 8 more points to the displayed table.
    2. Activate the first collector, and select the point Rubber End 1.
    3. Select the next available intermediate points sequentially for the other collectors, with the last collector being linked to Rubber End 2.


      Figure 6. Rubber strip model

Add Constraints

In this step, you will create constraints for the rubber strip model.

  1. Create a Fixed joint and one end of the rubber strip and a Translational joint at the other using the specifics in Table 2:
    Table 2.
    S Label Variable name Type Body 1 Body 2 Origin(s) Orientation Method Reference 1 Reference 2
    1 Fix Joint j_fix Fixed Joint Rubber Strip Ground Body Rubber End 1      
    2 Translation Joint j_trans Translational Joint Rubber Strip Ground Body Rubber End 2 Alignment axis (Vector) Global X  
  2. Create a motion on the Translational joint to apply axial pull. Use the specifics in Table 3:
    Table 3.
    S.No Label Variable Name Define Motion Joint Property
    1 Axial Motion mot_axial On Joint Translational Joint Displacement
  3. In the Properties tab for the motion, for Define by: choose Expression. In the Expression field, enter `15*time`.


    Figure 7. Motion expression
  4. In the Project Browser, under the Datasets folder, click Solver Gravity and change the Gravity option to Off.
  5. On the Standard toolbar, click (Save Model) and save the model as rubber_strip.mdl to your <working directory>.

Add Outputs

Now you will create outputs to measure engineering strain and engineering stress values.

Engineering Strain =

Engineering Stress =

  1. Right-click the (Outputs) icon.
  2. In the dialog, specify a Label of Eng strain (F2), Eng Stress (F3), and click OK.
  3. In the panel, for F2 type the expression `(DM({j_trans.i.idstring},{j_fix.i.idstring})-50)/50`. For F3, type the expression `MOTION({mot_axial.idstring},{0},{2},{0})/50`.


    Figure 8. Output requests
    Note: In the expression F2, the solver function DM() measures the distance magnitude between two markers; the I marker of the Translation Joint and the I marker of the Fix Joint. Expression F3 uses the solver function MOTION() which measures the reaction force due to the imposed motion Axial Motion.

Solve and Post-Process the Model

  1. Click the (Run) panel icon.
  2. Specify the MotionSolve file name as rubber_yeoh.xml.
  3. Specify the Simulation type as Transient, the End time as 4 seconds, and the Print interval as 0.01.
  4. Click Run.
  5. As the solver is executing, a warning message similar to the one shown below may be displayed:
    WARNING: Maximum vonMises strain exceeded maximum strain (YS) specified 
             for NLFE element BEAM12 (id=20000009) on Body_Flexible (id=30102)
             at time=3.922E+00
             Maximum strain Computed : 2.002E+00
             Maximum strain Specified: 2.000E+00
             Future warning for yield strain violation suppressed

    This message states that the maximum vonMises strain in your NLFE component exceeded what was specified (2.0) at time 3.92s. This message lets you know if your component is deforming more than what you would expect it to, which allows you to inspect your results and make corrections in modeling your system if required.

  6. After the simulation is completed, click Animate to view the animation in HyperView.
    You can use the (Start/Pause Animation) button to play the animation.
  7. Go to the Contour panel, select NLFE Stress (t), XX and click Apply.


    Figure 9. Stress contour
  8. Return to the MotionView Run panel and click the Plot button to load the .abf file in HyperGraph.
  9. Plot Engineering Stress vs Engineering Strain by selecting the data shown in Table 4 and Table 5 in HyperGraph.
    Table 4.
    X-Axis Data
    X Type Expression
    X Request REQ/70000000 Eng strain(F2), Eng Stress(F3)
    X Component F2
    Table 5.
    Y-Axis Data
    Y Type Expression
    Y Request REQ/70000000 Eng strain(F2), Eng Stress(F3)
    Y Component F3


    Figure 10. Stress versus strain curve
    Note: The animation shows the true stress due to cross-section deformation.
  10. Click to save the model.
  11. Click and save the session as hyperelastic.mvw.

References

JUSSI T, SOPANEN and AKI M. MIKKOLA:

Description of Elastic Forces in Absolute Nodal Coordinate Formulation. Journal of Nonlinear Dynamics 34: 53– 74, 2003.

Oleg Dmitrochenko:

Finite elements using absolute nodal coordinates for large deformation flexible multibody dynamics. Proceedings of the Third International Conference on Advanced Computational Methods in Engineering (ACOMEN 2005).

Sung Pil Jung, TaeWon Park, Won Sun Chung:

Dynamic analysis of rubber like material using absolute nodal coordinate formulation based on the non-linear constitutive law. Journal of Nonlinear Dyn (2011) 63: 149–157.