Files
김경종 b7f84e1c0f
Tests / Hermetic test suite (push) Has been cancelled
Tests / Skill frontmatter validation (push) Has been cancelled
add documents
2026-05-29 15:59:56 +09:00

363 lines
21 KiB
Markdown
Raw Permalink Blame History

This file contains ambiguous Unicode characters
This file contains Unicode characters that might be confused with other characters. If you think that this is intentional, you can safely ignore this warning. Use the Escape button to reveal them.
<!-- source-page: 131 -->
# 22.5.2 HYPERELASTIC BEHAVIOR IN ELASTOMERIC FOAMS
Products: Abaqus/Standard Abaqus/Explicit Abaqus/CAE
# References
• “Material library: overview,” Section 21.1.1
• “Elastic behavior: overview,” Section 22.1.1
• “Energy dissipation in elastomeric foams,” Section 22.6.2
• \*HYPERFOAM
• \*UNIAXIAL TEST DATA
• \*BIAXIAL TEST DATA
• \*PLANAR TEST DATA
• \*VOLUMETRIC TEST DATA
• \*SIMPLE SHEAR TEST DATA
• \*MULLINS EFFECT
• “Creating a hyperfoam material model” in “Defining elasticity,” Section 12.9.1 of the Abaqus/CAE Users Guide, in the HTML version of this guide
# Overview
The elastomeric foam material model:
• is isotropic and nonlinear;
• is valid for cellular solids whose porosity permits very large volumetric changes;
• optionally allows the specification of energy dissipation and stress softening effects (see “Energy dissipation in elastomeric foams,” Section 22.6.2);
• can deform elastically to large strains, up to 90% strain in compression; and
• requires that geometric nonlinearity be accounted for during the analysis step (see “Defining an analysis,” Section 6.1.2, and “General and linear perturbation procedures,” Section 6.1.3), since it is intended for finite-strain applications.
Abaqus/Explicit also provides a separate foam material model intended to capture the strain-rate sensitive behavior of low-density elastomeric foams such as used in crash and impact applications (see “Low-density foams,” Section 22.9.1).
# Mechanical behavior of elastomeric foams
Cellular solids are made up of interconnected networks of solid struts or plates that form the edges and faces of cells. Foams are made up of polyhedral cells that pack in three dimensions. The foam cells can be either open (e.g., sponge) or closed (e.g., flotation foam). Common examples of elastomeric
<!-- source-page: 132 -->
foam materials are cellular polymers such as cushions, padding, and packaging materials that utilize the excellent energy absorption properties of foams: the energy absorbed by foams is substantially greater than that absorbed by ordinary stiff elastic materials for a certain stress level.
Another class of foam materials is crushable foams, which undergo permanent (plastic) deformation. Crushable foams are discussed in “Crushable foam plasticity models,” Section 23.3.5.
Foams are commonly loaded in compression. Figure 22.5.21 shows a typical compressive stressstrain curve.
![](images/page-132_ee3f890578500594c0e769fc6e8107794fcbb227b52f006d0d8b009cccef1d85.jpg)
<details>
<summary>text_image</summary>
STRESS
Densification
Plateau: Elastic buckling
of cell walls
Cell wall bending
STRAIN
</details>
Figure 22.5.21 Typical compressive stress-strain curve.
Three stages can be distinguished during compression:
1. At small strains ( 5%) the foam deforms in a linear elastic manner due to cell wall bending.
2. The next stage is a plateau of deformation at almost constant stress, caused by the elastic buckling of the columns or plates that make up the cell edges or walls. In closed cells the enclosed gas pressure and membrane stretching increase the level and slope of the plateau.
3. Finally, a region of densification occurs, where the cell walls crush together, resulting in a rapid increase of compressive stress. Ultimate compressive nominal strains of 0.7 to 0.9 are typical.
The tensile deformation mechanisms for small strains are similar to the compression mechanisms, but they differ for large strains. Figure 22.5.22 shows a typical tensile stress-strain curve. There are two stages during tension:
1. At small strains the foam deforms in a linear, elastic manner as a result of cell wall bending, similar to that in compression.
2. The cell walls rotate and align, resulting in rising stiffness. The walls are substantially aligned at a tensile strain of about . Further stretching results in increased axial strains in the walls.
<!-- source-page: 133 -->
![](images/page-133_49eaa8e8e69577c2669bf339f83b5cc0aa589268a14627b079187c5e5093577a.jpg)
<details>
<summary>text_image</summary>
STRESS
Cell wall alignment
Cell wall bending
STRAIN
</details>
Figure 22.5.22 Typical tensile stress-strain curve.
At small strains for both compression and tension, the average experimentally observed Poissons ratio, $\nu ,$ of foams is 1/3. At larger strains it is commonly observed that Poissons ratio is effectively zero during compression: the buckling of the cell walls does not result in any significant lateral deformation. However, is nonzero during tension, which is a result of the alignment and stretching of the cell walls.
The manufacture of foams often results in cells with different principal dimensions. This shape anisotropy results in different loading responses in different directions. However, the hyperfoam model does not take this kind of initial anisotropy into account.
# Strain energy potential
In the elastomeric foam material model the elastic behavior of the foams is based on the strain energy function
$$
U = \sum_ {i = 1} ^ {N} \frac {2 \mu_ {i}}{\alpha_ {i} ^ {2}} \left[ \hat {\lambda} _ {1} ^ {\alpha_ {i}} + \hat {\lambda} _ {2} ^ {\alpha_ {i}} + \hat {\lambda} _ {3} ^ {\alpha_ {i}} - 3 + \frac {1}{\beta_ {i}} \left((J ^ {e \ell}) ^ {- \alpha_ {i} \beta_ {i}} - 1\right) \right],
$$
where N is a material parameter; $\mu _ { i } , \alpha _ { i }$ , and $\beta _ { i }$ are temperature-dependent material parameters;
$$
\hat {\lambda} _ {i} = (J ^ {t h}) ^ {- \frac {1}{3}} \lambda_ {i} \quad \rightarrow \quad \hat {\lambda} _ {1} \hat {\lambda} _ {2} \hat {\lambda} _ {3} = J ^ {e \ell};
$$
and $\lambda _ { i }$ are the principal stretches. The elastic and thermal volume ratios, $J ^ { e \ell }$ and $J ^ { t h }$ , are defined below. The coefficients $\mu _ { i }$ are related to the initial shear modulus, $\mu _ { 0 }$ , by
<!-- source-page: 134 -->
$$
\mu_ {0} = \sum_ {i = 1} ^ {N} \mu_ {i},
$$
while the initial bulk modulus, $K _ { 0 }$ , follows from
$$
K _ {0} = \sum_ {i = 1} ^ {N} 2 \mu_ {i} (\frac {1}{3} + \beta_ {i}).
$$
For each term in the energy function, the coefficient $\beta _ { i }$ determines the degree of compressibility. $\beta _ { i }$ is related to the Poissons ratio, $\nu _ { i } ,$ by the expressions
$$
\beta_ {i} = \frac {\nu_ {i}}{1 - 2 \nu_ {i}}, \qquad \nu_ {i} = \frac {\beta_ {i}}{1 + 2 \beta_ {i}}.
$$
Thus, if $\beta _ { i }$ is the same for all terms, we have a single effective Poissons ratio, . This effective Poissons ratio is valid for finite values of the logarithmic principal strains $\varepsilon _ { 1 } , \varepsilon _ { 2 } , \varepsilon _ { 3 } ,$ ; in uniaxial tension $\varepsilon _ { 2 } = \varepsilon _ { 3 } =$ $- \nu \varepsilon _ { 1 }$ .
# Thermal expansion
Only isotropic thermal expansion is permitted with the hyperfoam material model.
The elastic volume ratio, $J ^ { e \ell }$ , relates the total volume ratio (current volume/reference volume), J, and the thermal volume ratio, $J ^ { t h }$ :
$$
J ^ {e \ell} = \frac {J}{J ^ {t h}}.
$$
$J ^ { t h }$ is given by
$$
J ^ {t h} = (1 + \varepsilon^ {t h}) ^ {3},
$$
where $\varepsilon ^ { t h }$ is the linear thermal expansion strain that is obtained from the temperature and the isotropic thermal expansion coefficient (“Thermal expansion,” Section 26.1.2).
# Determining the hyperfoam material parameters
The response of the material is defined by the parameters in the strain energy function, $U _ { \star }$ these parameters must be determined to use the hyperfoam model. Two methods are provided for defining the material parameters: you can specify the hyperfoam material parameters directly or specify test data and allow Abaqus to calculate the material parameters.
The elastic response of a viscoelastic material (“Time domain viscoelasticity,” Section 22.7.1) can be specified by defining either the instantaneous response or the long-term response of such a material. To define the instantaneous response, the experiments outlined in the “Experimental tests” section that
<!-- source-page: 135 -->
follows have to be performed within time spans much shorter than the characteristic relaxation time of the material.
Input File Usage: \*HYPERFOAM, MODULI=INSTANTANEOUS
Abaqus/CAE Usage: Property module: material editor: Mechanical→Elasticity→Hyperfoam: Moduli time scale (for viscoelasticity): Instantaneous
If, on the other hand, the long-term elastic response is used, data from experiments have to be collected after time spans much longer than the characteristic relaxation time of the viscoelastic material. Long-term elastic response is the default elastic material behavior.
Input File Usage: \*HYPERFOAM, MODULI=LONG TERM
Abaqus/CAE Usage: Property module: material editor: Mechanical→Elasticity→Hyperfoam: Moduli time scale (for viscoelasticity): Long-term
# Direct specification
When the parameters $N , \mu _ { i } , \alpha _ { i }$ , and $\nu _ { i }$ are specified directly, they can be functions of temperature.
The default value of $\nu _ { i }$ is zero, which corresponds to an effective Poissons ratio of zero. The incompressible limit corresponds to all $\nu _ { i } ~ \to ~ 0 . 5$ . However, this material model should not be used for approximately incompressible materials: use of the hyperelastic model (“Hyperelastic behavior of rubberlike materials,” Section 22.5.1) is recommended if the effective Poissons ratio $\nu _ { e f f } > 0 . 4 5$ .
Input File Usage: \*HYPERFOAM, N=n ( )
Abaqus/CAE Usage: Property module: material editor: Mechanical→Elasticity→Hyperfoam: Strain energy potential order: $\pmb { \mathrm { ~ \mathscr ~ { ~ n ~ } ~ } } ( n \leq \mathrm { ~ \bf ~ 6 ~ } )$ ; optionally, toggle on Use temperature-dependent data
# Test data specification
The value of N and the experimental stress-strain data can be specified for up to five simple tests: uniaxial, equibiaxial, simple shear, planar, and volumetric. Abaqus contains a capability for obtaining the $\mu _ { i } , \alpha _ { i }$ , and $\beta _ { i }$ for the hyperfoam model with up to six terms (N=6) directly from test data. Poisson effects can be included either by means of a constant Poissons ratio or through specification of volumetric test data and/or lateral strains in the other test data.
It is important to recognize that the properties of foam materials can vary significantly from one batch to another. Therefore, all of the experiments should be performed on specimens taken from the same batch of material.
This method does not allow the properties to be temperature dependent.
As many data points as required can be entered from each test. Abaqus will then compute $\mu _ { i } , \alpha _ { i } ,$ , and, if necessary, $\nu _ { i } .$ . The technique uses a least squares fit to the experimental data so that the relative error in the nominal stress is minimized.
It is recommended that data from the uniaxial, biaxial, and simple shear tests (on samples taken from the same piece of material) be included and that the data points cover the range of nominal strains expected to arise in the actual loading. The planar and volumetric tests are optional.
<!-- source-page: 136 -->
For all tests the strain data, including the lateral strain data, should be given as nominal strain values (change in length per unit of original length). For the uniaxial, equibiaxial, simple shear, and planar tests, stress data are given as nominal stress values (force per unit of original cross-sectional area). The tests allow for both compression and tension data; compressive stresses and strains should be entered as negative values. For the volumetric tests the stress data are given as pressure values.
<table><tr><td>Input File Usage:</td><td>Use the first option to define an effective Poissons ratio ( $\nu_i = \nu$ for all $i$ ), or use the second option to define the lateral strains as part of the test data input:*HYPERFOAM, N=n, POISSON=ν, TEST DATA INPUT ( $n \leq 6$ )*HYPERFOAM, N=n, TEST DATA INPUT ( $n \leq 6$ ).In addition, use at least one and up to five of these additional options to give the experimental stress-strain data (see “Experimental tests” below):*UNIAXIAL TEST DATA*BIAXIAL TEST DATA*PLANAR TEST DATA*SIMPLE SHEAR TEST DATA*VOLUMETRIC TEST DATA</td></tr><tr><td>Abaqus/CAE Usage:</td><td>Property module: material editor: Mechanical→Elasticity→Hyperfoam: toggle on Use test data; Strain energy potential order: $n$ ( $n \leq 6$ ); optionally, toggle on Use constant Poisson&#x27;s ratio: and enter a value for the effective Poisson&#x27;s ratio ( $\nu_i = \nu$ for all $i$ )In addition, use at least one and up to five of the suboptions to give the experimental stress-strain data (see “Experimental tests” below):Suboptions→Uniaxial Test DataSuboptions→Biaxial Test DataSuboptions→Planar Test DataSuboptions→Simple Shear Test DataSuboptions→Volumetric Test Data</td></tr></table>
# Experimental tests
For a homogeneous material, homogeneous deformation modes suffice to characterize the material parameters. Abaqus accepts test data from the following deformation modes:
• Uniaxial tension and compression
• Equibiaxial tension and compression
• Planar tension and compression (pure shear)
• Simple shear
• Volumetric tension and compression
<!-- source-page: 137 -->
The stress-strain relations are defined in terms of the nominal stress (the force divided by the original, undeformed area) and the nominal, or engineering, strains, $\epsilon _ { i }$ . The principal stretches, $\lambda _ { i } ,$ , are related to the principal nominal strains, $\epsilon _ { i } ,$ , by
$$
\lambda_ {i} = 1 + \epsilon_ {i}.
$$
# Uniaxial, equibiaxial, and planar tests
The deformation gradient, expressed in the principal directions of stretch, is
$$
\mathbf {F} = \left[ \begin{array}{c c c} \lambda_ {1} & 0 & 0 \\ 0 & \lambda_ {2} & 0 \\ 0 & 0 & \lambda_ {3} \end{array} \right],
$$
where $\lambda _ { 1 } , \lambda _ { 2 }$ , and $\lambda _ { 3 }$ are the principal stretches: the ratios of current length to length in the original configuration in the principal directions of a material fiber. The deformation modes are characterized in terms of the principal stretches, $\lambda _ { i }$ , and the volume ratio, $J = \operatorname* { d e t } ( \mathbf { F } )$ . The elastomeric foams are not incompressible, so that $J = \lambda _ { 1 } \lambda _ { 2 } \lambda _ { 3 } \neq 1$ . The transverse stretches, $\lambda _ { 2 }$ and/or $\lambda _ { 3 }$ , are independently specified in the test data either as individual values from the measured lateral deformations or through the definition of an effective Poissons ratio.
The three deformation modes use a single form of the nominal stress-stretch relation,
$$
T _ {L} = \frac {\partial U}{\partial \lambda_ {L}} = \frac {2}{\lambda_ {L}} \sum_ {i = 1} ^ {N} \frac {\mu_ {i}}{\alpha_ {i}} (\lambda_ {L} ^ {\alpha_ {i}} - J ^ {- \alpha_ {i} \beta_ {i}}),
$$
where $T _ { L }$ is the nominal stress and $\lambda _ { L }$ is the stretch in the loading direction. Because of the compressible behavior, the planar mode does not result in a state of pure shear. In fact, if the effective Poissons ratio is zero, planar deformation is identical to uniaxial deformation.
# Uniaxial mode
In uniaxial mode $\lambda _ { 1 } = \lambda _ { U } , \lambda _ { 2 } = \lambda _ { 3 }$ , and $J = \lambda _ { U } \lambda _ { 2 } ^ { 2 }$ .
Input File Usage: \*UNIAXIAL TEST DATA
Abaqus/CAE Usage: Property module: material editor: Mechanical→Elasticity→Hyperfoam: toggle on Use test data, Suboptions→Uniaxial Test Data
# Equibiaxial mode
In equibiaxial mode $\lambda _ { 1 } = \lambda _ { 2 } = \lambda _ { B }$ and $J = \lambda _ { B } ^ { 2 } \lambda _ { 3 }$
Input File Usage: \*BIAXIAL TEST DATA
Abaqus/CAE Usage: Property module: material editor: Mechanical→Elasticity→Hyperfoam: toggle on Use test data, Suboptions→Biaxial Test Data
<!-- source-page: 138 -->
# Planar mode
In planar mode $\lambda _ { 1 } = \lambda _ { P } , \lambda _ { 2 } = 1$ , and $J = \lambda _ { P } \lambda _ { 3 }$ . Planar test data must be augmented by either uniaxial or biaxial test data.
Input File Usage: \*PLANAR TEST DATA
Abaqus/CAE Usage: Property module: material editor: Mechanical→Elasticity→Hyperfoam: toggle on Use test data, Suboptions→Planar Test Data
# Simple shear tests
Simple shear is described by the deformation gradient
$$
\mathbf {F} = \left[ \begin{array}{c c c} 1 & \gamma & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{array} \right],
$$
where $\gamma$ is the shear strain. For this deformation $J = \operatorname* { d e t } ( \mathbf { F } ) = 1$ . A schematic illustration of simple shear deformation is shown in Figure 22.5.23.
![](images/page-138_9ce950973eca3aa5b75f6049b544fddc4fed88996340f6d3cd3253859d862409.jpg)
<details>
<summary>text_image</summary>
shear strain
γ = Δx/h
Δx
σ₂₂=Tₜ
τ=Tₛ
σ₁₁
l
F
fixed distance h
</details>
Figure 22.5.23 Simple shear test.
The nominal shear stress, $T _ { S }$ , is
$$
T _ {S} = \frac {\partial U}{\partial \gamma} = \sum_ {j = 1} ^ {2} \biggl \{\frac {2 \gamma}{2 (\lambda_ {j} ^ {2} - 1) - \gamma^ {2}} \sum_ {i = 1} ^ {N} \frac {\mu_ {i}}{\alpha_ {i}} (\lambda_ {j} ^ {\alpha_ {i}} - 1) \biggr \},
$$
where $\lambda _ { j }$ are the principal stretches in the plane of shearing, related to the shear strain by
<!-- source-page: 139 -->
$$
\lambda_ {1, 2} = \sqrt {1 + \frac {\gamma^ {2}}{2} \pm \gamma \sqrt {1 + \frac {\gamma^ {2}}{4}}}.
$$
The stretch in the direction perpendicular to the shear plane is $\lambda _ { 3 } = 1$ The transverse (tensile) stress, $T _ { T }$ , developed during simple shear deformation due to the Poynting effect, is
$$
T _ {T} = \frac {\partial U}{\partial \varepsilon} = \sum_ {j = 1} ^ {2} \left\{\frac {2 (\lambda_ {j} ^ {2} - 1)}{2 \lambda_ {j} ^ {4} - \lambda_ {j} ^ {2} (\gamma^ {2} + 2)} \sum_ {i = 1} ^ {N} \frac {\mu_ {i}}{\alpha_ {i}} (\lambda_ {j} ^ {\alpha_ {i}} - 1) \right\}.
$$
Input File Usage: \*SIMPLE SHEAR TEST DATA
Abaqus/CAE Usage: Property module: material editor: Mechanical→Elasticity→Hyperfoam: toggle on Use test data, Suboptions→Simple Shear Test Data
# Volumetric tests
The deformation gradient, , is the same for volumetric tests as for uniaxial tests. The volumetric deformation mode consists of all principal stretches being equal;
$$
\lambda_ {1} = \lambda_ {2} = \lambda_ {3} = \lambda_ {V}, J = \lambda_ {V} ^ {3}.
$$
The pressure-volumetric ratio relation is
$$
- p = \frac {\partial U}{\partial J} = \frac {2}{J} \sum_ {i = 1} ^ {N} \frac {\mu_ {i}}{\alpha_ {i}} (J ^ {\frac {1}{3} \alpha_ {i}} - J ^ {- \alpha_ {i} \beta_ {i}}).
$$
A volumetric compression test is illustrated in Figure 22.5.24. The pressure exerted on the foam specimen is the hydrostatic pressure of the fluid, and the decrease in the specimen volume is equal to the additional fluid entering the pressure chamber. The specimen is sealed against fluid penetration.
Input File Usage: \*VOLUMETRIC TEST DATA
Abaqus/CAE Usage: Property module: material editor: Mechanical→Elasticity→Hyperfoam: toggle on Use test data, Suboptions→Volumetric Test Data
# Difference between compression and tension deformation
For small strains ( 5%) foams behave similarly for both compression and tension. However, at large strains the deformation mechanisms differ for compression (buckling and crushing) and tension (alignment and stretching). Therefore, accurate hyperfoam modeling requires that the experimental data used to define the material parameters correspond to the dominant deformation modes of the problem being analyzed. If compression dominates, the pertinent tests are:
• Uniaxial compression
• Simple shear
<!-- source-page: 140 -->
![](images/page-140_3518ddc6c9fd57f9038e33b0ce821439f001ec31b0429b849782b44c8fb16442.jpg)
<details>
<summary>text_image</summary>
volumetric gauge
pressure gauge
valve
pump
fluid
foam
rigid pressure chamber
</details>
Figure 22.5.24 Volumetric compression test.
• Planar compression (if Poissons ratio $\nu \neq 0 )$
• Volumetric compression (if Poissons ratio $\nu \neq 0 )$
If tension dominates, the pertinent tests are:
• Uniaxial tension
• Simple shear
• Biaxial tension (if Poissons ratio $\nu \neq 0 )$
• Planar tension (if Poissons ratio $\nu \neq 0 )$
Lateral strain data can also be used to define the compressibility of the foam. Measurement of the lateral strains may make other tests redundant; for example, providing lateral strains for a uniaxial test eliminates the need for a volumetric test. However, if volumetric test data are provided in addition to the lateral strain data for other tests, both the volumetric test data and the lateral strain data will be used in determining the compressibility of the foam. The hyperfoam model may not accurately fit Poissons ratio if it varies significantly between compression and tension.