Defining frequency domain viscoelasticity

This type of material model describes frequency-dependent material behavior in small, steady-state harmonic oscillations. In these cases the dissipative losses caused by internal damping effects must be modeled in the frequency domain.

Context:

The dissipative part of the material behavior is defined by giving the real and imaginary parts of g* and k* (for compressible materials) as functions of frequency. The moduli can be defined as functions of the frequency in one of three ways: by a power law, by tabular input, or by a Prony series expression for the shear and bulk relaxation moduli.

For more information on frequency domain viscoelasticity, see Frequency Domain Viscoelasticity.

  1. From the menu bar in the Edit Material dialog box, select MechanicalElasticityViscoelastic.

    (For information on displaying the Edit Material dialog box, see Creating or editing a material.)

  2. Click the arrow to the right of the Domain field, and select Frequency.
  3. Click the arrow to the right of the Frequency field, and select the option of your choice for determining viscoelastic material parameters:

    • Select Formula to define the frequency dependence by the power law formulae.

    • Select Tabular to define the frequency response in tabular form. You must provide the real and imaginary parts of ωg* and ωk*—where ω is the circular frequency—as functions of frequency in cycles per time.

    • Select Prony if you want Abaqus to calculate the frequency dependence from a time domain Prony series description of the dimensionless shear and bulk relaxation moduli.

    • Select Creep test data if you want Abaqus to calculate the parameters in the Prony series from creep test data that you provide. If you select this option, you must enter shear test data and/or volumetric test data in the Test Data Editor.

    • Select Relaxation test data if you want Abaqus to calculate the Prony series parameters from relaxation test data. If you select this option, you must enter shear test data and/or volumetric test data in the Test Data Editor.

  4. If you selected Creep test data or Relaxation test data from the Frequency option list, you can specify two additional parameters related to the calibration of Prony series parameters:

    • Click the arrows to the right of the Maximum number of terms in the Prony series field to specify the maximum number of terms (N) in the Prony series. Abaqus will perform the least-squares fit from N=1 to N= NMAX until convergence is achieved for the lowest N with respect to the error tolerance.

    • In the Allowable average root-mean-square error field, enter the error tolerance of the data points in the least-squares fit.

  5. If you selected Tabular from the Frequency option list, you can specify two additional options:

    Type

    This parameter specifies whether you are defining continuum material properties or effective thickness-direction gasket properties.

    Select Isotropic if you are defining continuum material properties. This choice is appropriate when the viscoelasic material model is used for any continuum, structural, or special-purpose elements whose material response is modeled using continuum material properties.

    Select Traction if you are defining effective thickness-direction gasket properties. This option is supported only for gasket elements whose behavior is modeled directly using a gasket behavior model.

    Preload

    This parameter specifies the nature of preload used for defining frequency domain viscoelastic material properties or effective thickness-direction gasket properties.

    Select Uniaxial to specify that the material properties correspond to a uniaxial test.

    Select Volumetric to specify that the material properties correspond to a volumetric test. This setting cannot be used to define effective thickness-direction gasket properties.

    Select Uniaxial and Volumetric to specify that the material properties correspond to both types of tests. This setting cannot be used to define effective thickness-direction gasket properties.

    Select None if you choose not to specify a preload parameter.

  6. If you selected Formula from the list of Frequency options, enter the following data in the Data table:

    g1*real

    Real part of g1*.

    g1*imag

    Imaginary part of g1*.

    a

    Value of a.

    k1*real

    Real part of k1*.

    k1*imag

    Imaginary part of k1*. If the material is incompressible, this value is ignored.

    b

    Value of b. If the material is incompressible, this value is ignored.

    You may need to expand the dialog box to see all the columns in the Data table. For detailed information on how to enter data, see Entering tabular data.

  7. If you selected Tabular from the list of Frequency options, enter the data relevant to your Type and Preload selections (not all of the following parameters will apply):

    Omega g* real

    Real part of ωg*. (ω(g*)=G/G)

    Omega g* imag

    Imaginary part of ωg*. (ω(g*)=1-Gs/G)

    Omega k* real

    Real part of ωk*. (ω(k*)=K/K) If the material is incompressible, this value is ignored.

    Omega k* imag

    Imaginary part of ωk*. (ω(k*)=1-Ks/K) If the material is incompressible, this value is ignored.

    Frequency

    Frequency, f, in cycles per time.

    Loss Modulus

    Uniaxial or bulk loss modulus.

    Storage Modulus

    Uniaxial or bulk storage modulus.

    Uniaxial Strain

    Uniaxial nominal strain (defines the level of uniaxial preload).

    Volume Ratio

    Volume ratio, J (current volume/original volume; defines the level of volumetric preload).

    Normalized Loss Modulus

    Real part of ωk*. (ω(k*)=K/K) for thickness-direction gasket behavior with no preload.

    Normalized Shear Modulus

    Imaginary part of ωk*. (ω(k*)=1-Ks/K) for thickness-direction gasket behavior with no preload.

    You may need to expand the dialog box to see all the columns. For detailed information on how to enter data, see Entering tabular data.

  8. If you selected Prony from the list of Frequency options, enter the following data in the Data table:

    g_i Prony

    g¯iP, the modulus ratio in the first term in the Prony series expansion of the shear relaxation modulus.

    k_i Prony

    k¯iP, the modulus ratio in the first term in the Prony series expansion of the bulk relaxation modulus.

    tau_i Prony

    τi, the relaxation time for the first term in the Prony series expansion.

    You may need to expand the dialog box to see all the columns in the Data table. For detailed information on how to enter data, see Entering tabular data.

  9. If applicable, use the Test Data menu to specify test data from which to define viscoelastic behavior. See the following sections for details:

  10. Click OK to create the material and to close the Edit Material dialog box. Alternatively, you can select another material behavior to define from the menus in the Edit Material dialog box (see Browsing and modifying material behaviors, for more information).