2.2.3 Mullins effect

Products: ABAQUS/Standard  ABAQUS/Explicit  

I. 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 Energy Potential FormPrimary Hyperelastic CoefficientsMullins Effect Parameters
Compressible Arruda-Boyce = 200.0, = 5.0,  = 0.001r = 1.1, m = 100.0, = 0.1
Compressible Ogden = 160.0, = 2.0, = 40.0, = –2.0, = 0.001r = 5.0, m = 220.0, = 0.1
Incompressible Ogden = 160.0, = 2.0, = 40.0, = –2.0r = 5.0, m = 220.0
User-defined hyperelastic materialSame as the compressible Yeoh model.
Compressible Van der Waals = 200.0, = 10.0, = 0.1, = 0.0, D = 0.001r = 3.0, m = 100.0, = 0.1
Compressible Yeoh = 1.326, = –0.326,  = 0.1319, = 0.000725r = 1.1, m = 100.0, = 0.1
Incompressible Yeoh = 1.326, = –0.326,  = 0.1319r = 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. In certain situations the value of the parameter m is found to be significantly different from the expected value (based on which the test data were originally generated). This occurs in situations where the deformation level leads to a relatively large value of the maximum deviatoric strain energy density, , such that the value of dominates over the value of m. The overall fit is found to be very good in these cases.

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 ).

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_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.

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.

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.

II. Energy dissipation in elastomeric foams

Elements tested

CPS4R    C3D8R    T3D2   

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 effect of energy dissipation in elastomeric foams. 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. 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. The results are then imported to ABAQUS/Explicit, where the element is unloaded. These results are again imported back into ABAQUS/Standard, where the element is loaded to deformation levels higher than the prior loading. Finally, the last set of results are imported from ABAQUS/Standard to ABAQUS/Standard, and then the element is 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:

Coefficients for Primary Elastomeric Foam BehaviorMullins Effect Parameters
= –1048.43, = 0.3025, = 532.20, = 0.3958, =517.027, =0.2135, = 0.2, = 0.2, = 0.2r = 1.75, m = 0.3, = 0.6

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. In certain situations the value of the parameter m is found to be significantly different from the expected value (based on which the test data were originally generated). This occurs in situations where the deformation level leads to a relatively large value of the maximum strain energy density, , such that the value of dominates over the value of m. The overall fit is found to be very good in these cases.

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

mmecdo1cut_hfoam.inp

T3D2 element, cyclic uniaxial tension.

mmecdo2cut_hfoam.inp

CPS4R element, cyclic uniaxial tension.

mmecdo3cut_hfoam.inp

C3D8R element, cyclic uniaxial tension.

mmecdo3cbt_hfoam.inp

C3D8R element, cyclic biaxial tension.

mmecdo3cpt_hfoam.inp

C3D8R element, cyclic planar loading.

mmetdo3cut_hfoam.inp

Calibration test with uniaxial unloading test data, C3D8R element, cyclic uniaxial tension.

mmetdo3cbt_hfoam.inp

Calibration test with biaxial unloading test data, C3D8R element, cyclic biaxial tension.

mmetdo3cpt_hfoam.inp

Calibration test with planar unloading test data, C3D8R element, cyclic planar tension.

mmetdo3cpt_m_hfoam.inp

Calibration test with unloading test data from uniaxial, biaxial, and planar tests and with the value of the parameter m fixed; C3D8R element; cyclic planar tension.

x_mmecdo1cut_hfoam.inp

Explicit dynamic test, T3D2 element, cyclic uniaxial tension.

x_mmecdo2cut_hfoam.inp

Explicit dynamic test, CPS4R element, cyclic uniaxial tension.

x_mmecdo3cut_hfoam.inp

Explicit dynamic test, C3D8R element, cyclic uniaxial tension.

x_mmecdo3cbt_hfoam.inp

Explicit dynamic test, C3D8R element, cyclic biaxial tension.

x_mmecdo3cpt_hfoam.inp

Explicit dynamic test, C3D8R element, cyclic planar loading.

x_mmetdo3cut_hfoam.inp

Calibration test with uniaxial unloading test data, explicit dynamic, C3D8R element, cyclic uniaxial tension.

x_mmetdo3cbt_hfoam.inp

Calibration test with biaxial unloading test data, explicit dynamic, C3D8R element, cyclic biaxial tension.

x_mmetdo3cpt_hfoam.inp

Calibration test with planar unloading test data, explicit dynamic, C3D8R element, cyclic planar tension.

x_mmetdo3cpt_m_hfoam.inp

Calibration test with unloading test data from uniaxial, biaxial, and planar tests and with the value of the parameter m fixed; explicit dynamic; C3D8R element; cyclic planar tension.

sx_s_mullins_hfoam.inp

Base problem for carrying out import from ABAQUS/Standard to ABAQUS/Explicit, C3D8R element, cyclic uniaxial tension.

sx_x_mullins_hfoam_n_y.inp

Explicit dynamic continuation of sx_s_mullins_hfoam.inp with only the state imported, C3D8R element, cyclic uniaxial tension.

xs_s_mullins_hfoam_n_y.inp

Import into ABAQUS/Standard from sx_x_mullins_hfoam_n_y.inp with only the state imported, C3D8R element, cyclic uniaxial tension.

ss_mullins_hfoam_n_y.inp

ABAQUS/Standard to ABAQUS/Standard import from xs_s_mullins_hfoam_n_y.inp with only the state imported, C3D8R element, cyclic uniaxial tension.