Defining Isotropic Viscoelasticity
For isotropic frequency domain viscoelasticity, Abaqus assumes that the shear (deviatoric) and volumetric behaviors are independent in multiaxial stress states.
Defining the Shear Behavior
Consider a shear test at small strain, in which a harmonically varying shear strain γ is applied:
where γ0 is the amplitude, i=√-1 , ω is the circular frequency, and t is time. We assume that the specimen has been oscillating for a very long time so that a steady-state solution is obtained. The solution for the shear stress then has the form
where Gs and Gℓ are the shear storage and loss moduli. These moduli can be expressed in terms of the (complex) Fourier transform g*(ω) of the nondimensional shear relaxation function g(t)=GR(t)G∞-1 :
where GR(t) is the time-dependent shear relaxation modulus, ℜ(g*) and ℑ(g*) are the real and imaginary parts of g*(ω) , and G∞ is the long-term shear modulus. See Frequency domain viscoelasticity for details.
The above equation states that the material responds to steady-state harmonic strain with a stress of magnitude Gsγ0 that is in phase with the strain and a stress of magnitude Gℓγ0 that lags the excitation by 90∘ . Therefore, we can regard the factor
as the complex, frequency-dependent shear modulus of the steadily vibrating material. The absolute magnitude of the stress response is
and the phase lag of the stress response is
Measurements of |τ| and ϕ as functions of frequency in an experiment can, therefore, be used to define Gs and Gℓ and, thus, ℜ(g*) and ℑ(g*) as functions of frequency.
Unless stated otherwise explicitly, all modulus measurements are assumed to be “true” quantities.
Defining the Volumetric Behavior
For isotropic viscoelasticity, the frequency dependence of the shear (deviatoric) and volumetric behaviors are independent. The volumetric behavior is defined by the bulk storage and loss moduli Ks(ω) and Kℓ(ω) . Similar to the shear moduli, these moduli can also be expressed in terms of the (complex) Fourier transform k*(ω) of the nondimensional bulk relaxation function k(t) :
where K∞ is the long-term elastic bulk modulus.
Large-Strain Viscoelasticity
The linearized vibrations can also be associated with an elastomeric material whose long-term (elastic) response is nonlinear and involves finite strains (a hyperelastic material). We can retain the simplicity of the steady-state small-amplitude vibration response analysis in this case by assuming that the linear expression for the shear stress still governs the system, except that now the long-term shear modulus G∞ can vary with the amount of static prestrain ˉγ :
The essential simplification implied by this assumption is that the frequency-dependent part of the material's response, defined by the Fourier transform g*(ω) of the relaxation function, is not affected by the magnitude of the prestrain. Therefore, strain and frequency effects are separated, which is a reasonable approximation for many materials.
Another implication of the above assumption is that the anisotropy of the viscoelastic moduli has the same strain dependence as the anisotropy of the long-term elastic moduli. Therefore, the viscoelastic behavior in all deformed states can be characterized by measuring the (isotropic) viscoelastic moduli in the undeformed state.
In situations where the above assumptions are not reasonable, the data can be specified based on measurements at the prestrain level about which the steady-state dynamic response is desired. In this case you must measure Gs , Gℓ , and G∞ (likewise Ks , Kℓ , and K∞ ) at the prestrain level of interest. Alternatively, the viscoelastic data can be given directly in terms of uniaxial and volumetric storage and loss moduli that might be specified as functions of frequency and prestrain (see Direct Specification of Storage and Loss Moduli for Large-Strain Isotropic Viscoelasticity below.)
The generalization of these concepts to arbitrary three-dimensional deformations is provided in Abaqus/Standard by assuming that the frequency-dependent material behavior has two independent components: one associated with shear (deviatoric) straining and the other associated with volumetric straining. In the general case of a compressible material, the model is, therefore, defined for kinematically small perturbations about a predeformed state as
and
where
- S
-
is the deviatoric stress, S=σ+pI ;
- p
-
is the equivalent pressure stress, p=-13trace(σ) ;
- Δ▽(JS)
-
is the part of the stress increment caused by incremental straining (as distinct from the part of the stress increment caused by incremental rotation of the preexisting stress with respect to the coordinate system);
- J
-
is the ratio of current to original volume;
- Δe
-
is the (small) incremental deviatoric strain, Δe=Δε-13ΔεvolI ;
- ˙e
-
is the deviatoric strain rate, ˙e=˙ε-13˙εvolI ;
- Δεvol
-
is the (small) incremental volumetric strain, Δεvol=trace(Δε) ;
- ˙εvol
-
is the rate of volumetric strain, ˙εvol=trace(˙ε) ;
- CS|0
-
is the deviatoric tangent elasticity matrix of the material in its predeformed state (for example, C1212 is the tangent shear modulus of the prestrained material);
- Q|0
-
is the volumetric strain-rate/deviatoric stress-rate tangent elasticity matrix of the material in its predeformed state; and
- K|0
-
is the tangent bulk modulus of the predeformed material.
For a fully incompressible material, only the deviatoric terms in the first constitutive equation above remain and the viscoelastic behavior is completely defined by g*(ω) .