Aim-
Given is a Ogden material card which is hyperelastic (rubber). Create a 10x10x10mm solid mesh and run a pure tensile run on it with stretch ratio of 5. Compare different element formulations of solids 1,-1,2,-2 and plot engineering stress vs stretch values.
Explanation-
Since this simulation is happening at a really slow pace (quasi static) the inertial effects can be neglected we can consider implicit for the analysis.
Since this a hyperelastic material even though we are running it in implicit , we need steps and iterations to achieve the equilibruim.
Element formulations used are 1,-1,2,-2 for solids.
ELFORM 1 - underintegrated constant stress ; needs hourglass stabilisation; efficient and accurate
ELFORM 2 - Reduced integrated ; no hour glass stabilisation needed; slower than ELFORM 1; shear locking can occur
ELFORM(-1 )- Similiar to ELFORM 2 , but shearlocking is accounted for; efficient formulation (CPU time)
ELFORM (-2) - Similiar to -1 , more accurate, more computational cost.
Procedure-
- We have to open the given LS-Dyna keyword (.k file) file in LS-PrePost, using option File>Open>LS-Dyna Keyword File as shown in below snap.
- Now we will open Ogden.k file and create the cube of size 10mm X 10mm X 10mm.
- First we will start with section card, to create section card from keyword>>all>>section>>solid>>here we will input id, elform type as shown below.
- As per given cases we have to change ELFORM as 1, -1 & -2, accordingly as shown in above snap or we can do it form .k file after saving it.
- Now we will check the MAT card given to us.
- Finally we have to assign the section and material to the part.
- Here we will fix the Bottom of the Cube and give the motion to the set of nodes on the Upper face.
- To fix the Bottom face, we will use different set of nodes and using single point constraint boundary, we can allow or constraint the nodes with different degrees of freedom.
- Boundary conditions using single point constrained for different set of nodes will be as given below:
- As shown in above image, constraints are given for the different set of nodes in a way that cube can shrink in X and Y direction upon tension in Z-Direction and also stick to its position.
- Thus , this boundary condition can be relatable to the real-time boundary condition.
- To apply the tensile load in Z-Direction to the cube, we have to give the motion to the set of node on the upper face of the cube.
- As shown in above image, we have given DOF = 3 , to allow translation in Z-Direction translation
- VAD = 2 is given for Displacement.
- We will define curve for displacement is given below, displacement is linearly applied at constat rate.
We know that Stretch = 1 + Engg strain
=> 5 = 1 + (change in length / 10)
=> Displacement(change in length) = 40mm
LS-Dyna employs both implicit and explicit methods to solve Static, Quasi-Static and Dynamic problems.
- These problems can differ in Time Increment approach.
- In implicit each Time increment has to converge, but we can set quite long time increments.
- In explicit, convergence for each time step is not needed.
- In implicit method, at each time increment, global equilibrium is established, due to iterative method.
- Implicit method is used for accuracy and can be applied for long duration problems.
- Also implicit methods are used when the deformation is not so large or the deformation is happening in a gradual manner.
- Here, we are stretch the cube up to 40 mm, is large deformation, but for the hyper elastic materials, this deformation can be small.
- So, implicit scheme is best for this type of analysis.
- Since this is the nonlinear implicit analysis, we need to create some implicit cards to activate implicit method like implicit auto, implicit general etc as given below:
- In *CONTOL_IMPLICIT_GENERAL, IMFLAG =1 is given to activate the implicit analysis method and initial Time step value is given as 0.1.
- In *CONTROL_IMPLICIT_AUTO, values of DTMIN and DTMAX are not defined, so that the solver can change the time step size as needed automatically.
- In CONTROL_TERINATION, ENDTIM is set to 10. Here termination does not represent the time but gives the value at which step the analysis will end.
- Now we will provide database output parameters to get the required results as d3plot and ascii plot files.
- Now we have to check the model from model checker.
- Since there is no error we can proceed further to save and run the keyword file.
- All simulations ended with Normal Termination. Therefore, all the models are simulated successfully.
- Now we can open this files one by one using LS-post processor.
Results & Plots-
Stresses-Z-stress
- We can see the stress results for different element formulations.
- Maximum stress we obtain is 4 Mpa along the Z-direction. As the tension is applied in Z-direction.
- We can observe the same value of maximum stress for all cases, but the distribution of stress is different in all cases.
Strains-Z-strain-
- Here also the maximum strain in Z-Direction is 1.6028, which is same for all cases.
Calculation of Nominal Stress & Stretch Ratio:
Now we will see the calculation of the Nominal Stress and Stretch Ratio for all cases.
- After getting the values of True strain and True stress from the plots, we can use the XY data and export it as “CSV” file and open it in EXCEL.
- Here all the data of True Stress and rue Strain are collected for all cases using an element results from the middle of the cube, which is Element 440.
- So all the results will be regarding Element ID: 440
- Here Stress results for Element 440 will be from Z-Stress
- Strain results for Element 440 will be from Lower Ipt Z-Strain.
Stretch Ratio:
- Here we can get the Stretch ratio from the True strain.
- To get the Stretch ratio, we can use formula as given below:
Nominal Stress:
- Here Nominal stress is the Engineering stress in Z-Direction.
- To get the Nominal Stress, we can use the formula as given below:
Calculations and the result data for all cases are same as given below:
Time |
True Strain |
True Stress |
Stretch Ratio |
Engg. Stress |
0.00E+00 |
0.00E+00 |
-2.37E-16 |
1.00E+00 |
1.00 |
-2.36978E-16 |
1.00E-01 |
7.69E-02 |
9.43E-02 |
1.08E+00 |
1.08 |
0.087286104 |
2.58E-01 |
1.88E-01 |
2.36E-01 |
1.21E+00 |
1.21 |
0.195760021 |
5.10E-01 |
3.42E-01 |
4.48E-01 |
1.41E+00 |
1.41 |
0.318247186 |
9.08E-01 |
5.45E-01 |
7.62E-01 |
1.72E+00 |
1.72 |
0.442103686 |
1.54E+00 |
8.00E-01 |
1.24E+00 |
2.23E+00 |
2.23 |
0.554992224 |
2.54E+00 |
1.10E+00 |
1.98E+00 |
3.02E+00 |
3.02 |
0.65511009 |
3.54E+00 |
1.34E+00 |
2.75E+00 |
3.81E+00 |
3.81 |
0.72297703 |
4.54E+00 |
1.53E+00 |
3.58E+00 |
4.60E+00 |
4.60 |
0.778281171 |
5.54E+00 |
1.60E+00 |
3.97E+00 |
4.97E+00 |
4.97 |
0.799731704 |
6.54E+00 |
1.60E+00 |
3.97E+00 |
4.97E+00 |
4.97 |
0.799720067 |
7.54E+00 |
1.60E+00 |
3.97E+00 |
4.97E+00 |
4.97 |
0.799722299 |
8.54E+00 |
1.60E+00 |
3.97E+00 |
4.97E+00 |
4.97 |
0.799723088 |
9.54E+00 |
1.60E+00 |
3.97E+00 |
4.97E+00 |
4.97 |
0.799722204 |
1.00E+01 |
1.60E+00 |
3.97E+00 |
4.97E+00 |
4.97 |
0.799722008 |
- From these values of stress and stretch ratios, we can create the graph of Nominal Stress Vs Stretch ratio for all element formulations.
Plots:
Plots of Nominal Stress Vs Stretch Ratio for Element formulation 1, 2,-1, and -2 as per case will be as given below:
CASE_1_Element Formulation = (1):
CASE_2_Element Formulation = (2):
CASE_3_Element Formulation = (-1):
CASE_4_Element Formulation = (-2):
Comparison:
- Plots obtained from the simulation results, we can observe the Stress and Strain results are same for different solid element formulations.
- We can observe same stretch ratio for all different element formulations.
- But we can see the distribution of stress over the length of the cube is in different manner while stretching. This is due to different behavior of element formulations.
- From above Snaps, we can observe different cycles and computational time for different element formulation.
- We can see that in case of element formulation -2, it is taking highest time for computation.
Conclusion:
- Since this is the uniaxial tensile test, we are getting similar results for different element formulations.
- Also the deformation we have given to the cube was not so large for this type of hyper elastic material, so the difference cannot be observed.
- Since the element formulation (1) is constant stress solid element, which default, cosserat point element is invoked, which is good to propose the strain energy function.
- Element formulation (2) is selective reduced integrated solid element, which give constant pressure throughout the element to avoid pressure locking, but type (1) is more accurate the type (2) for implicit problems.
- Element formulations type (-1) and (-2) are fully integrated elements intended for elements with poor aspect ratios. Among which type (-1) is more efficient and type (-2) is more accurate formulation.
- Since Element formulation (-1) has some side effects related to some particular deformation modes , we can conclude that Element formulation (-2) is more accurate and mostly preferred for the solid elements in Implicit problems.
Animations-
Z-Stress-
ELFORM_1
ELFORM_2