To define the gasket behavior by a gasket behavior model, you must assign a gasket section definition to a region of your model and assign the name of a gasket behavior definition to the gasket section definition. The gasket behavior for this region is defined entirely by the properties specified by the gasket behavior definition referring to the same name.

*GASKET SECTION, ELSET=name, BEHAVIOR=name
*GASKET BEHAVIOR, NAME=name

The thickness-direction, transverse shear, and membrane behaviors are defined to be uncoupled. Each behavior is specified independently.

You must specify the thickness-direction behavior. You can specify multiple thickness-direction behaviors to define the loading and unloading characteristics. You can obtain an average contact pressure output when the thickness-direction behavior is defined as force or force per unit length versus closure.

The transverse shear and membrane behaviors are optional for gasket elements that have all displacement degrees of freedom active at their nodes. You can define one or both of these behaviors.

When thermal and rate-dependent effects are important, you can define thermal expansion and creep behavior for gaskets; user subroutines UEXPAN and CREEP can be used to define these behaviors.

You cannot specify density for gasket elements since they have no mass matrix.

*GASKET BEHAVIOR, NAME=name
*GASKET THICKNESS BEHAVIOR
*GASKET ELASTICITY
*GASKET CONTACT AREA
*EXPANSION
*CREEP
*DEPVAR
*USER OUTPUT VARIABLES

To define the thickness-direction behavior of gaskets, ABAQUS/Standard offers a nonlinear elastic model with damage and a nonlinear elastic-plastic model with the possibility of considering creep effects. Thermal effects in the thickness direction can also be accounted for.

ABAQUS/Standard measures the thickness-direction deformation as the closure between the bottom and top faces of the gasket element; therefore, the thickness-direction behavior must always be defined in terms of closure. The closure is the sum of the elastic closure, plastic closure, creep closure, thermal closure, plus any initial gap in the thickness direction. As explained below, the behavior can be defined as pressure versus closure, force versus closure, or force per unit length versus closure. In all cases the thickness-direction behavior can be defined as a function of temperature and/or field variables.

*GASKET THICKNESS BEHAVIOR, DEPENDENCIES

*GASKET THICKNESS BEHAVIOR, VARIABLE=STRESS

*GASKET THICKNESS BEHAVIOR, VARIABLE=FORCE

Defining a nonlinear elastic model with damage

The nonlinear elastic model with damage is illustrated in Figure 26.6.6–2.

Figure 26.6.6–2 Elastic model with damage.

As the gasket is compressed, the pressure (or force, or force per unit length) follows the path given by the loading curve. If the gasket is unloaded, for example at point B, the pressure follows the unloading curve . Reloading after unloading follows the unloading curve until the loading is such that the closure becomes greater than , after which the loading path follows the loading curve . The arrows shown in the figure illustrate the loading/unloading paths of this model.

Defining the loading curve

To define the loading curve in piecewise linear form, you provide data points of pressure versus elastic closure, starting with point A. For negative elastic closures, the model gives zero pressure (or force). For closures larger than the last user-specified closure, the pressure-closure relationship is extrapolated based on the last slope computed from the user-specified data.

Input File Usage:

*GASKET THICKNESS BEHAVIOR, TYPE=DAMAGE, DIRECTION=LOADING

ABAQUS/CAE Usage:

Property module: material editor: OtherGasketGasket Thickness Behavior: Type: Damage, Loading

Defining the unloading curve

To define the unloading curves (, , and so on), you provide data points of pressure (or force) versus elastic closure up to a given maximum closure (, or , and so on). You can specify as many unloading curves as are necessary. Each unloading curve always starts at point A, the point of zero pressure for zero elastic closure, since the damaged elasticity model does not allow any permanent deformation. If unloading occurs from a maximum closure for which an unloading curve is not specified, the unloading is interpolated from neighboring unloading curves. The unloading curves are stored in normalized form so that they intersect the loading curve at a unit stress (or unit force) for a unit elastic closure, and the interpolation occurs between these normalized curves. If unloading curves are not specified, the loading/unloading will follow the loading curve.

Input File Usage:

*GASKET THICKNESS BEHAVIOR, TYPE=DAMAGE, DIRECTION=UNLOADING

ABAQUS/CAE Usage:

Property module: material editor: OtherGasketGasket Thickness Behavior: Type: Damage, Unloading, toggle on Include user-specified unloading curves

Defining the behavior for elements with an initial gap

For cases when the load in the gasket does not increase as soon as the gasket is compressed (see Figure 26.6.6–3), you can specify an initial gap as part of the gasket section property definition (see “Defining the gasket element's initial geometry,” Section 26.6.4) and define the loading/unloading curves as if the initial gap were not present (the case of Figure 26.6.6–2). This method is convenient when many gasket elements refer to the same gasket behavior and the only difference is the initial gap.

Figure 26.6.6–3 Elastic model with damage and initial gap.

Defining a nonlinear elastic-plastic model

The nonlinear elastic-plastic model is illustrated in Figure 26.6.6–4.

Figure 26.6.6–4 Elastic-plastic model.

As the gasket is compressed, the pressure (or force) follows the path given by the loading curve . The loading curve is a nonlinear elastic curve until point B is reached. At point B the slope of the loading curves decreases by more than 10%, which is assumed to correspond with the onset of plastic deformation. The value of 10% was chosen as a reasonable minimum value that can be expected at the onset of yield. If yield starts at a point at which no decrease in the slope occurs, numerical difficulties may occur. If the elastic part of the loading curve has a changing slope, the curve should be defined such that the slope does not decrease by more than 10% at any given point. After point B plastic deformation starts taking place. If unloading occurs before point B is reached, unloading will take place along the initial loading curve. Once loading has gone beyond point B, unloading will take place along an unloading curve such as curve . The unloading is assumed to be entirely elastic. The amount of closure at point D represents the plastic closure for the unloading curve . Reloading after unloading follows the same curve until the gasket yields, after which the loading curve is followed. Plastic deformation takes place until the last point M on the loading curve is reached. Beyond point M, the curve is followed for both loading and unloading; this behavior represents the behavior of a crushed gasket, which is assumed to be entirely elastic and can be specified in a piecewise-linear fashion, even beyond point M. The arrows shown in the figure illustrate the loading/unloading paths for the elastic-plastic model.

ABAQUS/Standard will automatically convert the curves so that the unloading curves become curves of pressure (or force) versus elastic closure for a given plastic closure. The loading curve will be transformed into an elastic loading/unloading curve defined at zero plastic closure (the portion of the curve) and a yield curve (the portion of the curve). By default, the onset of yield (point B) will be obtained as soon as the slope of the loading curve decreases by 10% from the maximum slope recorded up to that point while traveling along the loading curve from point A to point M. ABAQUS/Standard offers two alternatives to allow you to override this default method of determining the onset of yield as described below. If only a loading curve is provided, the unloading will be based on the curve , independent of the level of plasticity.

Defining the loading curve

To define the loading curve in piecewise linear form, you provide data points of pressure (or force, or force per unit length) versus closure (where closure represents the elastic plus the plastic closure), starting with point A. The last closure value given represents the closure at which the gasket is assumed crushed (point M in Figure 26.6.6–4); at this point, the maximum permanent deformation is reached. For negative closures the model gives zero pressure (or force).

To override the default method of determining the onset of yield, you can specify either a value for the decrease in slope other than the default value of 10% or the closure value at which onset of yield occurs. The specified value must correspond to a point on the loading curve at which the slope decreases.

Input File Usage:

Use the following option to define the loading curve and use the default method for determining the onset of yield:

*GASKET THICKNESS BEHAVIOR, TYPE=ELASTIC-PLASTIC, DIRECTION=LOADING Use the following option to define the loading curve and specify a nondefault value for the decrease in slope that defines the onset of yield: *GASKET THICKNESS BEHAVIOR, TYPE=ELASTIC-PLASTIC, DIRECTION=LOADING, SLOPE DROP=drop Use the following option to define the loading curve and specify the closure value that defines the onset of yield: *GASKET THICKNESS BEHAVIOR, TYPE=ELASTIC-PLASTIC, DIRECTION=LOADING, YIELD ONSET=closure_value

ABAQUS/CAE Usage:

Property module: material editor: OtherGasketGasket Thickness Behavior: Type: Elastic-Plastic, Loading, Yield onset method: Relative slope drop drop or Yield onset method: Closure value closure_value

Defining the unloading curve

To define the unloading curves (, , and so on), you provide data points of pressure (or force, or force per unit length) versus closure (elastic plus plastic) for each given plastic closure (closure at points D, F, and so on) in ascending values of closure. You can specify as many unloading curves as are necessary. If unloading occurs at a plastic closure for which an unloading curve is not specified, the unloading curve is interpolated from neighboring unloading curves. If no unloading curves are specified, unloading is assumed to follow a curve similar to the initial nonlinear elastic segment of the loading curve. The unloading curves are stored in normalized form so that they intersect the yield curve at a unit stress (or unit force) for a unit elastic closure, and the interpolation occurs between these normalized curves.

If the loading curve includes highly nonlinear behavior after the onset of yield, the interpolated unloading may give unreasonable behavior (such as the interpolated unloading path crossing over the user-defined loading curve). You should specify as many user-defined unloading curves as are needed to create regions for which interpolated unloading response is appropriate. For example, Figure 26.6.6–5 illustrates a loading curve that includes a sharp decrease in the hardening slope well after the onset of yield.

Figure 26.6.6–5 Elastic-plastic behavior with complex loading curve.

In this case it is insufficient to specify only one unloading curve at the gasket crush point (the end of the loading data). If unloading were to take place from point C, the unloading path would cross over the loading path. At least one additional unloading curve is required, after the sharp decrease in hardening slope, to prevent the interpolated unloading path crossing the loading curve.

Input File Usage:

*GASKET THICKNESS BEHAVIOR, TYPE=ELASTIC-PLASTIC, DIRECTION=UNLOADING

ABAQUS/CAE Usage:

Property module: material editor: OtherGasketGasket Thickness Behavior: Type: Elastic-Plastic, Unloading, toggle on Include user-specified unloading curves

Defining the behavior for elements with an initial gap

For cases when the load in the gasket does not increase as soon as the gasket is compressed (see Figure 26.6.6–6), you can specify an initial gap as part of the gasket section property definition (see “Defining the gasket element's initial geometry,” Section 26.6.4) and define the loading/unloading curves as if the initial gap were not present (the case of Figure 26.6.6–4). This method is convenient when many gasket elements refer to the same gasket behavior and the only difference is the initial gap.

Figure 26.6.6–6 Elastic-plastic model with initial gap.

*GASKET THICKNESS BEHAVIOR, DIRECTION=LOADING, 
TENSILE STIFFNESS FACTOR=factor

You can define the elastic transverse shear stiffness of the gasket. ABAQUS/Standard measures the relative displacement between the bottom and top of the gasket element along the local 2- or 3-directions to define the transverse shear in the gasket. Therefore, you should always define the elastic transverse stiffness as stress (or force, or force per unit length) per unit displacement. You can specify the stiffness as a function of temperature and field variables. The same stiffness is used for the shear in the 1–2 plane and the shear in the 1–3 plane. For each set of temperature and/or field variables, the first slope of the initial loading curve for the gasket's thickness-direction behavior will be used to compute the transverse shear stiffness if the transverse shear behavior is not defined explicitly.

*GASKET ELASTICITY, COMPONENT=TRANSVERSE SHEAR, DEPENDENCIES

*GASKET ELASTICITY, COMPONENT=TRANSVERSE SHEAR, VARIABLE=STRESS

*GASKET ELASTICITY, COMPONENT=TRANSVERSE SHEAR, VARIABLE=FORCE

You can define the linear elastic behavior of the gasket by giving Young's modulus and Poisson's ratio. These data can be provided as a function of temperature and/or field variables. If you do not specify the linear elastic behavior of the gasket, the gasket has no membrane stiffness. In this case you must ensure that the nodes of the elements are restrained adequately in the directions orthogonal to the thickness direction of the gasket.

*GASKET ELASTICITY, COMPONENT=MEMBRANE, DEPENDENCIES

You can define isotropic thermal expansion to specify the same coefficient of thermal expansion for the membrane and thickness-direction behaviors.

Alternatively, you can define orthotropic thermal expansion to specify three different coefficients of thermal expansion. The first coefficient will apply to the thermal expansion of the gasket in the thickness direction; the other two coefficients will apply to the expansion of the gasket in the local 2- and 3-directions, respectively.

The membrane thermal strains, , are obtained as explained in “Thermal expansion,” Section 20.1.2. ABAQUS/Standard computes the thermal closure for the thickness direction as

If user subroutine UEXPAN is used to define the thermal expansion of the gasket, the incremental thermal strains must be provided in the subroutine. The thermal closure will be obtained from the thermal strain in the thickness direction, as described above.

*EXPANSION, TYPE=ISO
*EXPANSION, TYPE=ORTHO

*EXPANSION, TYPE=ISO, USER
*EXPANSION, TYPE=ORTHO, USER

You can define creep behavior in the thickness direction of the gasket only when the elastic-plastic model (see “Defining a nonlinear elastic-plastic model” above) is used. The creep closure rate will be obtained as

If user subroutine CREEP is used to define the rate-dependent thickness-direction response of the gasket, the compressive creep strain increment must be provided in the subroutine. The creep closure will be obtained from the creep strain, as described above.

*CREEP

*CREEP, LAW=USER

You can define viscoelastic behavior in the thickness direction of the gasket only when the elastic-damage model (see “Defining a nonlinear elastic model with damage” above) is used. Only frequency domain viscoelastic behavior is supported. This behavior is useful for modeling the steady-state dynamic response of automotive components with gaskets about some pre-loaded base state, such as would be obtained at the end of a nonlinear sealing analysis, to determine the noise-vibration-harshness (NVH) characteristics of the system.

During the nonlinear sealing analysis step the frequency-domain viscoelastic behavior is ignored, and the material response is determined by the long-term elastic properties of the material. It is generally accepted (Zubeck and Marlow, 2002) that the dynamic stiffness and damping characteristics of automotive components such as gaskets and grommets vary with the frequency of excitation as well as the level of preload. These structural properties also depend on the geometry and the level of confinement of the gasket. This capability allows the direct specification of such dynamic properties as quantified by the effective storage and loss moduli in the thickness-direction, as tabular functions of the frequency of excitation and the level of preload. The preload is quantified by the amount of closure in the base state about which the steady-state dynamic response is desired.

In determining the dynamic response of the gasket, the long-term elastic response is assumed to be defined by the nonlinear elastic model with damage. The steady-state dynamic response is assumed to be a perturbation about a base state defined by this elastic damage behavior at a certain value of closure. The viscoelastic response can be specified using two approaches, as discussed below.

*VISCOELASTIC, TYPE=TRACTION, PRELOAD=UNIAXIAL

*VISCOELASTIC, TYPE=TRACTION

When the thickness-direction behavior of the gasket is defined in terms of force or force per unit length versus closure, ABAQUS/Standard will provide the thickness-direction force or force per unit length as output variable S11. In this case you can define either a contact width or contact area versus closure curve that will be used to obtain the average “contact” pressure at each integration point as output variable CS11. This average pressure considers the changing contact area that occurs as a result of the deformation of a gasket, as shown in Figure 26.6.6–1. The closure used for input of this curve corresponds to the total mechanical closure, defined as the sum of the elastic, plastic, and creep closures.

When two- and three-dimensional link gasket elements are used, you should specify the contact area versus mechanical closure in tabular form. When axisymmetric link and three-dimensional line elements are used, you should specify the contact width versus mechanical closure in tabular form. A typical curve is shown in Figure 26.6.6–7.

*GASKET CONTACT AREA, DEPENDENCIES

Output variable E is usually used in ABAQUS/Standard to output strain. For gasket elements with behavior defined by a gasket behavior model this output variable has thickness-direction and transverse shear components with units of displacement and membrane strains. Output variable NE is used to output an effective strain. The effective strain components are computed as follows:

The output variables THE, PE, or CE can also be used for gasket elements to output generalized thermal strains, plastic strains, or creep strains, respectively.

For all stress/strain output variables the 11-component refers to the through-thickness direction; the 22-, 33- and 23-components refer to two direct and one shear membrane component, respectively; the remaining 12- and 13-components refer to the transverse shear components. For details about these definitions, see “Gasket elements: overview,” Section 26.6.1.

The output of the elastic strain energy (output variable ALLSE) also contains the energy due to damage or change in elasticity as a function of plasticity. Therefore, this energy is usually not fully recoverable.