ProductsAbaqus/StandardAbaqus/Explicit Material behaviorThe parallel rheological framework allows definition of a nonlinear viscoelastic-elastoplastic model consisting of multiple networks connected in parallel, as shown in Figure 1. Figure 1. Nonlinear viscoelastic-elastoplastic model with multiple parallel
networks.
The number of viscoelastic networks, N, can be arbitrary; however, at most one equilibrium network (network in Figure 1) is allowed in the model. The equilibrium network response might be purely elastic or elastoplastic. In addition, it might include Mullins effect to predict material softening. The definition of the equilibrium network is optional. If it is not defined, the stress in the material will relax completely over time. The model can be used to predict complex behavior of materials subjected to finite strains, which cannot be modeled accurately using other models available in Abaqus. An example of such complex behavior is depicted in Figure 2, which shows normalized stress relaxation curves for three different strain levels. This behavior can be modeled accurately using the nonlinear viscoelastic model depicted in Figure 3, which can be defined within the framework; but it cannot be captured with the linear viscoelastic model (see Time domain viscoelasticity). In the latter case, the three curves would coincide. Figure 2. Normalized stress relaxation curves for three different strain
levels.
Figure 3. Nonlinear viscoelastic model with multiple parallel networks.
Elastic behaviorThe elastic part of the response for all the networks is specified using the hyperelastic material model. Any of the hyperelastic models available in Abaqus can be used (see Hyperelastic behavior of rubberlike materials). The same hyperelastic material definition is used for all the networks, scaled by a stiffness ratio specific to each network. Consequently, only one hyperelastic material definition is required by the model along with the stiffness ratio for each network. The elastic response can be specified by defining either the instantaneous response or the long-term response. Equilibrium network behaviorIn addition to the elastic response described above, the response of the equilibrium network can include plasticity and Mullins effect to predict material softening. If the plastic response is defined using isotropic hardening, the response in the equilibrium network is equivalent to that of the permanent set model available in Abaqus (see Permanent set in rubberlike materials for a detailed description of the model). In Abaqus/Standard the nonlinear kinematic hardening model with multiple backstresses can be specified in addition to isotropic plastic hardening. The nonlinear kinematic hardening model is a generalization of the model used for metal plasticity. See Models for metals subjected to cyclic loading for a detailed description of the model, with the difference that the Cauchy stress is replaced with the Kirchhoff stress in the current formulation. Viscous behaviorViscous behavior must be defined for each viscoelastic network. It is modeled by assuming the multiplicative split of the deformation gradient and the existence of the creep potential, , from which the flow rule is derived. In the multiplicative split the deformation gradient is expressed as where is the elastic part of the deformation gradient (representing the hyperelastic behavior) and is the creep part of the deformation gradient (representing the stress-free intermediate configuration). The creep potential is assumed to have the general form where is the Cauchy stress. If the potential is specified, the flow rule can be obtained from where is the symmetric part of the velocity gradient, , expressed in the current configuration and is the proportionality factor. In this model the creep potential is given by and the proportionality factor is taken as , where is the equivalent deviatoric Cauchy stress and is the equivalent creep stain rate. In this case the flow rule has the form or, equivalently where is the Kirchhoff stress, is the determinant of , is the deviatoric Cauchy stress, is the deviatoric Kirchhoff stress, and . To complete the derivation, the evolution law for must be provided. In this model can be determined from either a power-law strain hardening model or a hyperbolic-sine model. Power law modelThe power law model is available in the form where
Power-law strain hardening modelThe power-law strain hardening model is available in the form where
Hyperbolic-sine law modelThe hyperbolic-sine law is available in the form where
Bergstrom-Boyce modelThe Bergstrom-Boyce model is available in the form where and
The response of the network defined by the Bergstrom-Boyce model is very similar to the response of the time-dependent network in the hysteresis model (see Hysteresis in elastomers). However, there are also important differences between the models. In the Bergstrom-Boyce model the equivalent Kirchhoff stress is used instead of the equivalent Cauchy stress, which is used in the hysteresis model. (The two stress measures become equivalent for the case of incompressible materials.) In addition, the material parameters, A, in the models differ by a factor of . The parameter in the hysteresis model must be multiplied by to make the parameters equivalent. User-defined model in Abaqus/StandardA user-defined creep model is available of the following general form: where and
The tensor is defined in UCREEPNETWORK. Thermal expansionOnly isotropic thermal expansion is permitted with nonlinear viscoelastic materials (Thermal expansion). Defining viscoelastic responseThe nonlinear viscoelastic response is defined by specifying the identifier, stiffness ratio, and creep law for each viscoelastic network. Specifying network identifierEach viscoelastic network in the material model must be assigned a unique network identifier or network id. The network identifiers must be consecutive integers starting with 1. The order in which they are specified is not important. Input File Usage Use the following option to specify the network identifier: VISCOELASTIC, NONLINEAR, NETWORKID=networkId Defining the stiffness ratioThe contribution of each network to the overall response of the material is determined by the value of the stiffness ratio, , which is used to scale the elastic response of the network material. The sum of the stiffness ratios of the viscoelastic networks must be smaller than or equal to 1. If the sum of the ratios is equal to 1, the purely elastic equilibrium network is not created. If the sum of the ratios is smaller than 1, the equilibrium network is created with a stiffness ratio, , equal to where denotes the number of viscoelastic networks and is the stiffness ratio of network . You can specify the stiffness ratio to remain constant during the analysis or to vary as a function of temperature and predefined field variables. Defining a constant stiffness ratioYou can specify that the stiffness ratio remains constant during the analysis: Input File Usage VISCOELASTIC, NONLINEAR, SRATIO=ratio Defining a temperature- and field-variable dependent stiffness ratioAlternatively, you can define the stiffness ratio as a function of temperature and predefined field variables. Input File Usage NETWORK STIFFNESS RATIO, N=nNetworks, DEPENDENCIES=n Specifying the creep lawThe definition of creep behavior in Abaqus/Standard is completed by specifying the creep law. Power law creep modelThe power law model is defined by specifying five material parameters: , n, m, a, and . The parameter must be positive. It is introduced for dimensional consistency, and its default value is 1.0. For physically reasonable behavior and n must be positive, a must be nonnegative (the default is 0.0), and −1 < m ≤ 0. Input File Usage VISCOELASTIC, NONLINEAR, LAW=POWER LAW Strain hardening power law creep modelThe strain hardening law is defined by specifying three material parameters: A, n, and m. For physically reasonable behavior A and n must be positive and −1 < m ≤ 0. Input File Usage VISCOELASTIC, NONLINEAR, LAW=STRAIN Hyperbolic sine creep modelThe hyperbolic sine creep law is specified by providing three nonnegative parameters: A, B, and n. Input File Usage VISCOELASTIC, NONLINEAR, LAW=HYPERB Bergstrom-Boyce creep modelThe Bergstrom-Boyce creep law is specified by providing four parameters: A, m, C, and E. The parameters A and E must be nonnegative, the parameter m must be positive, and the parameter C must lie in . Input File Usage VISCOELASTIC, NONLINEAR, LAW=BERGSTROM-BOYCE User-defined creep modelAn alternative method for defining the creep law involves using user subroutine UCREEPNETWORK in Abaqus/Standard or VUCREEPNETWORK in Abaqus/Explicit. Optionally, you can specify the number of property values needed as data in the user subroutine. Input File Usage VISCOELASTIC, NONLINEAR, LAW=USER, PROPERTIES=n Numerical difficultiesDepending on the choice of units, the value of A in the strain power-law, hyperbolic-sine, and Bergstrom-Boyce models may be very small for typical creep strain rates. If A is less than 10−27, numerical difficulties can cause errors in the material calculations; therefore, a different system of units should be used to avoid such difficulties in the calculation of creep strain increments. In such cases it is recommended to replace the strain power-law model with the power law model, which does not have the limitation described above. The strain power-law model is a special case of the power law model obtained by setting , , and . Thermo-rheologically simple temperature effectsThermo-rheologically simple temperature effects can be included for each viscoelastic network. In this case the creep law is modified and takes the following form: where and denote the reduced time and the shift function, respectively. The reduced time is related to the actual time through the integral differential equation Abaqus supports the following forms of the shift function: the Williams-Landel-Ferry (WLF) form and the Arrhenius form (see Thermo-rheologically simple temperature effects). In addition, user-defined forms can be specified in Abaqus/Standard. User-defined form in Abaqus/StandardAn alternative method for specifying the shift function involves using user subroutine UTRSNETWORK. Optionally, you can specify the number of property values needed as data in the user subroutine. Input File Usage TRS, DEFINITION=USER, PROPERTIES=n Material response in different analysis stepsIn Abaqus/Standard the material is active during all stress/displacement procedure types. However, the creep effects are taken into account only in quasi-static (Quasi-static analysis), coupled temperature-displacement (Fully coupled thermal-stress analysis), direct-integration implicit dynamic (Implicit dynamic analysis using direct integration), and steady-state transport (Steady-state transport analysis) analyses. If the material is used in a steady-state transport analysis, it cannot include plasticity. In other stress/displacement procedures the evolution of the state variables is suppressed and the creep strain remains unchanged. In Abaqus/Explicit the creep effects are always active. ElementsThe parallel rheological framework is available with continuum elements that include mechanical behavior (elements that have displacement degrees of freedom), except for one-dimensional elements. The parallel rheological framework is also supported with elements that use the plane stress formulation such as solid plane stress elements, membranes, and shells. However, those elements are not supported with compressible materials. If a compressible material is specified with plane stress elements, Abaqus will modify the material to make it incompressible and issue an informational message. OutputIn addition to the standard output identifiers available in Abaqus (Abaqus/Standard output variable identifiers and Abaqus/Explicit output variable identifiers), the following variables have special meaning for the nonlinear viscoelastic material model:
In the above definitions denotes the stiffness ratio for network , denotes the number of viscoelastic networks, the subscript or superscript is used to denote network quantities, and the network is assumed to be the purely elastic network. If plasticity is specified in the equilibrium network, the standard output identifiers available in Abaqus corresponding to other isotropic and kinematic hardening plasticity models can be obtained for this model as well. In addition, if the Mullins effect is used in the model, the output variables available for the Mullins effect model (see Mullins effect) can be requested. References
|