Mullins effect in elastomers
Elements tested
SAX1
CPS4R
CPE4R
CPE4RH
C3D8R
C3D8RH
T2D2
Problem description
The problems in this set can be broadly classified under three categories. The first category of problems consists of simple displacement- or load-controlled cyclic tests to verify the Mullins effect, with the primary response defined by different strain energy potential functions. The tests consist of a single element that is cyclically loaded to a maximum strain (stress) level, then unloaded to zero strain (stress). This is followed by further reloading to levels of strain (stress) that are higher than those reached during the loading segment of the first cycle, followed again by unloading to zero strain (stress). The tests in this section use parts and assemblies.
The second category of problems is intended for testing the calibration capabilities for determining the Mullins effect coefficients. The problems use unloading test data that were generated by running a model with specified values of the Mullins effect coefficients. The calibration capability is meant to recover the specified values of the Mullins effect coefficients. These tests use different loading states, such as uniaxial tension, biaxial tension, and planar tension.
The third category of problems tests the import capability with the Mullins effect. All tests in this section are set up with a uniaxial stress state. The tests consist of first loading a single element in Abaqus/Standard and unloading it. The results are then imported into Abaqus/Explicit, where the element is loaded to deformation levels higher than the original loading and then unloaded. These results are again imported back into Abaqus/Standard, where the element is loaded to deformation levels higher than the prior loading and then unloaded. Finally, the last set of results are imported from Abaqus/Standard to Abaqus/Standard, and the element is further deformed and unloaded. The above series of tests includes problems that import both the state and the reference configuration, problems that import only the state, and problems that import neither the state nor the reference configuration.
Material:
The following material data are used for the first category of tests:
Strain energypotential form | Primary hyperelastic coefficients | Mullins effect parameters |
---|---|---|
Compressible Arruda-Boyce | = 200.0, = 5.0, = 0.001 | r = 1.1, m = 100.0, = 0.1 |
Compressible Ogden | = 160.0, = 2.0, = 40.0, = –2.0, = 0.001 | r = 5.0, m = 220.0, = 0.1 |
Incompressible Ogden | = 160.0, = 2.0, = 40.0, = –2.0 | r = 5.0, m = 220.0 |
User-defined hyperelastic material | Same as the compressible Yeoh model | |
Compressible Van der Waals | = 200.0, = 10.0, = 0.1, = 0.0, D = 0.001 | r = 3.0, m = 100.0, = 0.1 |
Compressible Yeoh | = 1.326, = –0.326, = 0.1319, = 0.000725 | r = 1.1, m = 100.0, = 0.1 |
Incompressible Yeoh | = 1.326, = –0.326, = 0.1319 | r = 1.1, m = 100.0 |
For the second and third category of tests the primary material response is defined using the incompressible Yeoh potential with the deviatoric coefficients as given above. For the second category of tests the unloading test data are generated for uniaxial, biaxial, and planar stress states using the following values for the Mullins effect parameters: r = 1.25, m = 0.01, and = 0.9. These parameters are also used to define the Mullins effect in the third category of tests.
Loading:
The first category of problems includes both displacement- and force-controlled loading. The second and third categories of problems are carried out under only displacement-controlled loading.
Results and discussion
For the first category of problems the results of the Abaqus/Standard and Abaqus/Explicit numerical simulations are in good agreement with the analytical results.
For the second category of problems, which tests the calibration of the Mullins effect parameters, it is observed that the parameters r and are always captured accurately. A good fit for m is obtained in situations where the deformation level leads to a relatively large value of maximum deviatoric strain energy density, , such that the value of dominates over the value of m.
For the final category of problems, which tests the import capability, the response after each import of results is as expected. When the state is imported, further deformation upon import shows the appropriate level of stress softening. On the other hand, when the state is not imported, no stress softening is observed.
Input files
- mmecdo2cut_arruda.inp
-
Compressible Arruda-Boyce model, CPE4RH element, cyclic uniaxial tension.
- mmecdo2cut_vdw.inp
-
Compressible Van der Waals model, CPE4RH element, cyclic uniaxial tension.
- mmecdo2cut_yeoh.inp
-
Compressible Yeoh model, SAX1 element, cyclic uniaxial tension.
- mmecdo2cut.inp
-
Compressible Yeoh model, CPS4R element, cyclic uniaxial tension, tests temperature- and field-variable-dependent Mullins effect material properties.
- mmecdo2cut_po.inp
-
Tests the POST OUTPUT capability for the damage-related output variables. This job needs the restart file from the job mmecdo2cut.inp.
- mmecoo2cut_yeoh.inp
-
Incompressible Yeoh model, SAX1 element, cyclic uniaxial tension.
- mmecoo2cut_user.inp
-
Incompressible Yeoh model, CPS4R element, cyclic uniaxial tension. The Mullins effect is implemented with user subroutine UMULLINS; the use of solution- dependent state variables in UMULLINS is also tested (the solution-dependent state variables are used to provide a nonzero initial value of ).
- mmecoo2cut_user.f
-
User subroutine UMULLINS to define damage variable for the Mullins effect material model.
- mmecdo3cut_ogden.inp
-
Compressible Ogden model, C3D8RH element, cyclic uniaxial tension.
- mmecoo3cut_ogden.inp
-
Incompressible Ogden model, C3D8RH element, cyclic uniaxial tension.
- mmecdo3cut_user.inp
-
Compressible user-defined hyperelastic material, C3D8RH element, cyclic uniaxial tension, user subroutine UHYPER provided in the file mmecdo3cut_user.f.
- mmecdo3cut_user.f
-
User subroutine UHYPER to define a compressible hyperelastic material model.
- mmecdo3cut_yeoh.inp
-
Compressible Yeoh model, C3D8RH element, cyclic uniaxial tension.
- mmecdo3cut_yeoh_load.inp
-
Compressible Yeoh model, C3D8RH element, cyclic uniaxial tension with load control.
- mmecdo3ctu.inp
-
Compressible Yeoh model, C3D8RH element, triaxial tension followed by unloading, further loading in uniaxial tension, and unloading. The purpose of this test is to demonstrate that a purely volumetric deformation does not result in any damage. This problem also tests a static linear perturbation analysis about a damaged base state.
- mmecoo3cut_yeoh.inp
-
Incompressible Yeoh model, C3D8RH element, cyclic uniaxial tension.
- mmecdo2cut_arruda_visco.inp
-
Compressible Arruda-Boyce model with viscoelasticity, CPE4RH element, cyclic uniaxial tension.
- mmecdo3cut_ogden_visco.inp
-
Compressible Ogden model with viscoelasticity, C3D8RH element, cyclic uniaxial tension.
- mmecoo2cut_yeoh_visco.inp
-
Incompressible Yeoh model with viscoelasticity, SAX1 element, cyclic uniaxial tension.
- mmecoo3cut_ogden_visco.inp
-
Incompressible Ogden model with viscoelasticity, C3D8RH element, cyclic uniaxial tension.
- neoh_mullins_ve.inp
-
Neo-Hookean model with viscoelasticity and Mullins effect; C3D8R, CPE4R, and CPS4R elements; uniaxial loading-unloading at different strain levels.
- x_mmecdo2cut_arruda.inp
-
Explicit dynamic test with compressible Arruda-Boyce model, CPE4RH element, cyclic uniaxial tension.
- x_mmecdo2cut_vdw.inp
-
Explicit dynamic test with compressible Van der Waals model, CPE4RH element, cyclic uniaxial tension.
- x_mmecdo3cut_ogden.inp
-
Explicit dynamic test with compressible Ogden model, C3D8RH element, cyclic uniaxial tension.
- x_mmecdo2cut_yeoh.inp
-
Explicit dynamic test with compressible Yeoh model, SAX1 element, cyclic uniaxial tension.
- x_mmecdo3cut_yeoh.inp
-
Explicit dynamic test with compressible Yeoh model, C3D8RH element, cyclic uniaxial tension.
- x_mmecdo2cut_visarruda.inp
-
Explicit dynamic test with compressible Arruda-Boyce model and viscoelasticity, CPE4RH element, cyclic uniaxial tension.
- x_mmecdo2cut_visvdw.inp
-
Explicit dynamic test with compressible Van der Waals model and viscoelasticity, CPE4RH element, cyclic uniaxial tension.
- x_mmecdo3cut_visogden.inp
-
Explicit dynamic test with compressible Ogden model and viscoelasticity, C3D8RH element, cyclic uniaxial tension.
- x_mmecdo3cut_visyeoh.inp
-
Explicit dynamic test with compressible Yeoh model and viscoelasticity, C3D8RH element, cyclic uniaxial tension.
- x_mmecdo3cut_user.inp
-
Explicit dynamic test with compressible Yeoh model and viscoelasticity, C3D8RH element, cyclic uniaxial tension. The Mullins effect is implemented with user subroutine VUMULLINS.
- x_mmecdo3cut_user.f
-
User subroutine VUMULLINS to define damage variable for the Mullins effect material model.
- mmetdo3cut.inp
-
Calibration test with uniaxial unloading test data, C3D8RH element, cyclic uniaxial tension.
- mmetdo3cut_marlow.inp
-
Calibration test with uniaxial unloading test data, C3D8RH element, cyclic uniaxial tension, Marlow model.
- mmetdo3cbt.inp
-
Calibration test with biaxial unloading test data, C3D8RH element, cyclic biaxial tension.
- mmetdo3cpt.inp
-
Calibration test with planar unloading test data, C3D8RH element, cyclic planar tension.
- mmetdo3cpt_mult.inp
-
Calibration test with unloading test data from uniaxial, biaxial, and planar tests; C3D8RH element; cyclic planar tension.
- mmetdo3cpt_r.inp
-
Calibration test with unloading test data from uniaxial, biaxial, and planar tests and with the value of the parameter r fixed; C3D8RH element; cyclic planar tension.
- mmetdo3cpt_m.inp
-
Calibration test with unloading test data from uniaxial, biaxial, and planar tests and with the value of the parameter m fixed; C3D8RH element; cyclic planar tension.
- mmetdo3cpt_beta.inp
-
Calibration test with unloading test data from uniaxial, biaxial, and planar tests and with the value of the parameter fixed; C3D8RH element; cyclic planar tension.
- sx_s_mullins.inp
-
Base problem for carrying out import from Abaqus/Standard to Abaqus/Explicit, C3D8RH element, cyclic uniaxial tension.
- sx_x_mullins_y_y.inp
-
Explicit dynamic continuation of sx_s_mullins.inp with both the reference configuration and the state imported, C3D8RH element, cyclic uniaxial tension.
- sx_x_mullins_n_y.inp
-
Explicit dynamic continuation of sx_s_mullins.inp with only the state imported, C3D8RH element, cyclic uniaxial tension.
- sx_x_mullins_n_n.inp
-
Explicit dynamic continuation of sx_s_mullins.inp without importing the state or the reference configuration, C3D8RH element, cyclic uniaxial tension.
- xs_s_mullins_y_y.inp
-
Import into Abaqus/Standard from sx_x_mullins_y_y.inp with both the state and the reference configuration imported, C3D8RH element, cyclic uniaxial tension.
- xs_s_mullins_n_y.inp
-
Import into Abaqus/Standard from sx_x_mullins_n_y.inp with only the state imported, C3D8RH element, cyclic uniaxial tension.
- xs_s_mullins_n_n.inp
-
Import into Abaqus/Standard from sx_x_mullins_n_n.inp without importing the state or the reference configuration, C3D8RH element, cyclic uniaxial tension.
- ss_mullins_y_y.inp
-
Abaqus/Standard to Abaqus/Standard import from xs_s_mullins_y_y.inp with both the state and the reference configuration imported, C3D8RH element, cyclic uniaxial tension.
- ss_mullins_n_y.inp
-
Abaqus/Standard to Abaqus/Standard import from xs_s_mullins_n_y.inp with only the state imported, C3D8RH element, cyclic uniaxial tension.
- ss_mullins_n_n.inp
-
Abaqus/Standard to Abaqus/Standard import from xs_s_mullins_n_n.inp without importing the state or the reference configuration, C3D8RH element, cyclic uniaxial tension.
- mmecdo1cut_marlow.inp
-
Compressible Marlow model, T2D2 element, cyclic uniaxial tension, tests temperature- and field-variable-dependent Mullins effect material properties.
- mmecdo2cut_marlow.inp
-
Compressible Marlow model, CPS4R element, cyclic uniaxial tension, tests temperature- and field-variable-dependent Mullins effect material properties.
- mmecdo3cut_marlow.inp
-
Compressible Marlow model, C3D8RH element, cyclic uniaxial tension.