Mechanical APDL Material Reference - Ansys
Mechanical APDL Material Reference - Ansys
ANSYS, ANSYS Workbench, Ansoft, AUTODYN, EKM, Engineering Knowledge Manager, CFX, FLUENT, HFSS and any
and all ANSYS, Inc. brand, product, service and feature names, logos and slogans are registered trademarks or
trademarks of ANSYS, Inc. or its subsidiaries in the United States or other countries. ICEM CFD is a trademark used
by ANSYS, Inc. under license. CFX is a trademark of Sony Corporation in Japan. All other brand, product, service
and feature names or trademarks are the property of their respective owners.
Disclaimer Notice
THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION INCLUDE TRADE SECRETS AND ARE CONFID-
ENTIAL AND PROPRIETARY PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS. The software products
and documentation are furnished by ANSYS, Inc., its subsidiaries, or affiliates under a software license agreement
that contains provisions concerning non-disclosure, copying, length and nature of use, compliance with exporting
laws, warranties, disclaimers, limitations of liability, and remedies, and other provisions. The software products
and documentation may be used, disclosed, transferred, or copied only in accordance with the terms and conditions
of that software license agreement.
ANSYS, Inc. is certified to ISO 9001:2008.
Third-Party Software
See the legal information in the product help files for the complete Legal Notice for ANSYS proprietary software
and third-party software. If you are unable to access the Legal Notice, please contact ANSYS, Inc.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. iii
ANSYS Mechanical APDL Material Reference
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
iv of ANSYS, Inc. and its subsidiaries and affiliates.
ANSYS Mechanical APDL Material Reference
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. v
ANSYS Mechanical APDL Material Reference
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
vi of ANSYS, Inc. and its subsidiaries and affiliates.
ANSYS Mechanical APDL Material Reference
List of Figures
3.1. Generalized Maxwell Solid in One Dimension ........................................................................................ 47
3.2. Sphere Discretization by 42 Microplanes ............................................................................................... 58
3.3. Damage Parameter d Depending on the Equivalent Strain Energy ......................................................... 60
3.4. Stress-strain Behavior at Uniaxial Tension .............................................................................................. 60
3.5. Pseudoelasticity (PE) and Shape Memory Effect (SME) ............................................................................ 89
3.6. Typical Superelasticity Behavior ............................................................................................................. 90
3.7. Idealized Stress-Strain Diagram of Superelastic Behavior ........................................................................ 91
3.8. Admissible Paths for Elastic Behavior and Phase Transformations ........................................................... 96
List of Tables
3.1. Linear Material Property Descriptions .................................................................................................... 15
3.2. Implicit Creep Equations ....................................................................................................................... 77
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. vii
ANSYS Mechanical APDL Material Reference
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
viii of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 1: Introduction to Material Models
Material models (also called constitutive models), are the mathematical representation of a material's
response to an applied load.
Typical model classes include the relationships between stress-strain, heat flux-temperature gradient,
voltage-strain, and current-voltage, but also include more general behaviors such as friction and
bonding, and response due to changes in the physical environment such as thermal expansion and
swelling.
This reference provides information about material model behavior and application, including details
about the load-response relationship and the necessary information required to use the material models
in an analysis. The models are grouped based on the degrees of freedom that, directly or indirectly,
give the loading function that serves as the input for the material model.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 1
Chapter 1: Introduction to Material Models
Type Description
Magnetic Gives the magnetomotive force in response to the magnetic
flux.
Conductivity For electric and magnetic current conductivities that model
the relationship between the respective field and its flux.
Permittivity and Per- Gives the energy storage in a material in response to an
meability electromagnetic field.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
2 of ANSYS, Inc. and its subsidiaries and affiliates.
Material Parameters
Type Description
Loss Represents the energy lost in electromagnetic and dielectric
materials in response to changes in electromagnetic fields.
Although the models are coupled, they remain distinct and give the same load-response behavior.
However, the piezoelectric and piezoresistive materials are electromechanical coupled models that give
a strain in response to a voltage and also produce a voltage in response to straining.
In some cases, the parameters can be a function of physical field quantities such as temperature, fre-
quency or time or interaction quantities such as normal pressure, relative distance, or relative velocity.
Matching the model to the actual behavior can be challenging; therefore, some built-in curve-fitting
methods are available that use minimization to select a set of parameters that give a close fit to measured
material behavior. The curve-fitting methods help you to select material parameters for creep, hypere-
lastic, viscoelastic, and some plastic models.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 3
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
4 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 2: Element Support for Material Models
Following is a list of available material models and the elements that support each material.
Material models are specified via the TB,Lab command, where Lab represents the material model label
(shortcut name).
Hill plasticity:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 5
Chapter 2: Element Support for Material Models
Hill plasticity:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
6 of ANSYS, Inc. and its subsidiaries and affiliates.
Label Material Model Elements
(Lab)
CREEP Creep Implicit creep with von Mises or Hill potential:
SOLID62, SOLID65
CZM Cohesive zone CONTA171, CONTA172, CONTA173, CONTA174, CONTA175,
CONTA176,CONTA177,INTER202,INTER203,INTER204,INTER205
DISCRETE Explicit spring- COMBI165
damper (discrete)
DMGE Damage evolution Progressive damage evolution (MPDG option):
law
LINK180, SHELL181, PLANE182, PLANE183, SOL-
ID185, SOLID186, SOLID187, BEAM188, BEAM189,
SOLSH190, SHELL208, SHELL209, SOLID272,
SOLID273, SHELL281, SOLID285, PIPE288,
PIPE289, ELBOW290
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 7
Chapter 2: Element Support for Material Models
Hill plasticity:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
8 of ANSYS, Inc. and its subsidiaries and affiliates.
Label Material Model Elements
(Lab)
LINK180, SHELL181, PLANE182, PLANE183, SOL-
ID185, SOLID186, SOLID187, BEAM188, BEAM189,
SOLSH190, SHELL208, SHELL209, PLANE223,
SOLID226, SOLID227, REINF264, REINF265, SOL-
ID272, SOLID273, SHELL281, SOLID285, EL-
BOW290
LSEM Anisotropic electric HF119, HF120
and magnetic loss
tangents
MELAS Multilinear elasticity SOLID62, SOLID65
MISO Multilinear isotropic von Mises plasticity:
hardening
SOLID62, SOLID65, LINK180, SHELL181,PLANE182,
PLANE183, SOLID185, SOLID186, SOLID187,
BEAM188, BEAM189, SOLSH190, SHELL208,
SHELL209, PLANE223, SOLID226, SOLID227, RE-
INF264, REINF265, SOLID272, SOLID273,
SHELL281, SOLID285, PIPE288, PIPE289, EL-
BOW290
Hill plasticity:
Hill plasticity:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 9
Chapter 2: Element Support for Material Models
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
10 of ANSYS, Inc. and its subsidiaries and affiliates.
Label Material Model Elements
(Lab)
PLANE223, SOLID226, SOLID227, PIPE288, PIPE289, ELBOW290,
REINF264, REINF265, SOLID272, SOLID273, SHELL281, SOLID285
SMA Shape memory alloy PLANE182, PLANE183, PLANE223 (with plane strain or axisym-
metric stress states),SOLID185,SOLID186,SOLID187,SOLSH190,
SOLID226, SOLID227, SOLID272, SOLID273, SOLID285
STATE State variables (user- LINK180,SHELL181,PLANE182,PLANE183,SOLID185,SOLID186,
defined) SOLID187,BEAM188,BEAM189,SOLSH190,SHELL208,SHELL209,
REINF264,REINF265,SOLID272,SOLID273,SHELL281,SOLID285,
PIPE288, PIPE289, ELBOW290
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 11
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
12 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 3: Material Models
This document describes all material models available for implicit analysis, including information about
material data table input (TB). The following material model topics are available:
3.1. Understanding Material Data Tables
3.2. Experimental Data
3.3. Linear Material Properties
3.4. Plasticity Models
3.5. Hyperelastic Material Models
3.6. Viscoelastic Material Model
3.7. Microplane Material Model
3.8. Porous Media
3.9. Magnetic Material Model
3.10. High-Frequency Electromagnetic Material Models
3.11. Anisotropic Elastic Material Model
3.12. Piezoelectric Material Model
3.13. Piezoresistive Material Model
3.14. Anisotropic Electric Permittivity Material Model
3.15. Rate-Dependent Plastic (Viscoplastic) Material Models
3.16. Gasket Material Model
3.17. Creep Equations
3.18. Swelling Model
3.19. Shape Memory Alloy (SMA) Material Model
3.20. MPC184 Joint Material Models
3.21. Contact Friction
3.22. Cohesive Zone Material Model
3.23. Fluid Material Models
3.24. User-Defined Material Model
3.25. Material Strength Limits
3.26. Damage Initiation Criteria
3.27. Damage Evolution Law
For a list of the elements that support each material model, see "Element Support for Material Models".
For related information, see "Nonlinear Structural Analysis" in the Structural Analysis Guide.
For information about explicit dynamics material models, including detailed data table input, see Ma-
terial Models in the ANSYS LS-DYNA User's Guide.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 13
Chapter 3: Material Models
For some element types, the data table is used for special element input data other than material
properties. The form of the data table (referred to as the TB table) depends upon the data being defined:
• Where the form is peculiar to only one element type, the table is described with the element in "Element
Library".
• If the form applies to more than one element, it is described here and referenced in the element de-
scription.
Initiate the experimental data table, then specify the appropriate experimental data type (TBOPT), as
follows:
Enter the data (TBPT) for each data point. Each data point entered consists of the independent variable
followed by one or more dependent variables. The specific definition of the input points depends on
the requirements of the material model associated with the experimental data.
You can also define experimental data as a function of field variables. Field-dependent data are entered
by preceding a set of experimental data (TBFIELD) to define the value of the field.
Conversely, if properties needed for a thermal analysis (such as KXX) are temperature-dependent, the
problem is nonlinear. Properties such as stress-strain data are described as nonlinear properties because
an analysis with these properties requires an iterative solution.
Linear material properties that are required for an element, but which are not defined, use default values.
(The exception is that EX and KXX must be input with a nonzero value where applicable.) Any additional
material properties are ignored.
The X, Y, and Z portions of the material property labels refer to the element coordinate system. In
general, if a material is isotropic, only the “X” and possibly the “XY” term is input.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
14 of ANSYS, Inc. and its subsidiaries and affiliates.
Linear Material Properties
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 15
Chapter 3: Material Models
MP,
Lab Units Description
Value
REFT Temp Reference temperature (as a property) (see also TREF)
Coefficient of friction (or, for FLUID29 and FLUID30 ele-
MU None
ments, boundary admittance)
ALPD None Mass matrix multiplier for damping (also see ALPHAD)
BETD None Stiffness matrix multiplier for damping (also see BETAD)
DMPR None Constant material damping coefficient
DENS Mass/Vol Mass density
KXX Thermal conductivity, element x direction
Heat*Length/
KYY Thermal conductivity, element y direction
(Time*Area*Temp)
KZZ Thermal conductivity, element z direction
C Heat/Mass*Temp Specific heat
ENTH Heat/Vol Enthalpy ( DENS*C d(Temp))
Heat /
HF Convection (or film) coefficient
(Time*Area*Temp)
EMIS None Emissivity
Heat/Time Heat generation rate for thermal mass element MASS71
QRATE Fraction of plastic work converted to heat (Taylor-Quinney
None coefficient) for coupled-field elements PLANE223, SOLID226,
and SOLID227
VISC Force*Time/ Length2 Viscosity
Sonic velocity (FLUID29, FLUID30, FLUID129, and FLUID130
SONC Length/Time
elements only)
MURX Magnetic relative permeability, element x direction
MURY None Magnetic relative permeability, element y direction
MURZ Magnetic relative permeability, element z direction
MGXX Magnetic coercive force, element x direction
MGYY Current/Length Magnetic coercive force, element y direction
MGZZ Magnetic coercive force, element z direction
RSVX Electrical resistivity, element x direction
RSVY Resistance*Area/Length Electrical resistivity, element y direction
RSVZ Electrical resistivity, element z direction
PERX Electric relative permittivity, element x direction
PERY None Electric relative permittivity, element y direction
PERZ Electric relative permittivity, element z direction
LSST None Dielectric loss tangent
SBKX Seebeck coefficient, element x direction
SBKY Voltage/Temp Seebeck coefficient, element y direction
SBKZ Seebeck coefficient, element z direction
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
16 of ANSYS, Inc. and its subsidiaries and affiliates.
Linear Material Properties
MP,
Lab Units Description
Value
DXX Diffusion coefficient, element x direction
2
DYY Length /Time Diffusion coefficient, element y direction
DZZ Diffusion coefficient, element z direction
3
CREF Mass/Length Saturated concentration
3
CSAT Mass/Length Reference concentration
BETX Coefficient of diffusion expansion, element x direction
3
BETY Length /Mass Coefficient of diffusion expansion, element y direction
BETZ Coefficient of diffusion expansion, element z direction
Poisson's ratio (PRXY or NUXY) defaults to 0.3. If a zero value is desired, input PRXY or
NUXY with a zero or blank value. Poisson's ratio should not be ≥ 0.5 nor ≤ -1.0.
The shear modulus (GXY) defaults to EX/(2(1+NUXY)). If GXY is input, it must match EX/(2
(1+NUXY)). The sole purpose for inputting GXY is to ensure consistency with the other
two properties.
EX, EY, EZ, (PRXY, PRYZ, PRXZ, or NUXY, NUYZ, NUXZ), GXY, GYZ, and GXZ must all be
input if the element type uses the material property. There are no defaults. For example,
if only EX and EY are input (with different values) to a plane stress element, The program
generates an error message indicating that the material is orthotropic and that GXY and
NUXY are also needed.
Poisson's ratio may be input in either major (PRXY, PRYZ, PRXZ) or minor (NUXY, NUYZ,
NUXZ) form, but not both for a particular material. The major form is converted to the
minor form during the solve operation (SOLVE). Solution output is in terms of the minor
form, regardless of how the data was input. If zero values are desired, input the labels
with a zero (or blank) value.
For axisymmetric analyses, the X, Y, and Z labels refer to the radial (R), axial (Z), and
hoop (θ) directions, respectively. Orthotropic properties given in the R, Z, θ system should
be input as follows: EX = ER, EY = EZ, and EZ = E θ. An additional transformation is re-
quired for Poisson's ratios. If the given R, Z, θ properties are column-normalized (see the
Mechanical APDL Theory Reference), NUXY = NURZ, NUYZ = NUZ θ = (ET/EZ) *NU θZ, and
NUXZ = NUR θ. If the given R, Z, θ properties are row-normalized, NUXY = (EZ/ER)*NURZ,
NUYZ = (E θ/EZ)*NUZ θ = NU θZ, and NUXZ = (E θ/ER)*NUR θ.
For all other orthotropic material properties (including ALPX, ALPY, and ALPZ), the X, Y,
and Z part of the label (as in KXX, KYY, and KZZ) refers to the direction (in the element
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 17
Chapter 3: Material Models
coordinate system) in which that particular property acts. The Y and Z directions of the
properties default to the X direction (for example, KYY and KZZ default to KXX) to reduce
the amount of input required.
3.3.3. Damping
Material dependent mass and stiffness damping (MP,ALPD and MP,BETD) is an additional method of
including damping for dynamic analyses and is useful when different parts of the model have different
damping values.
ALPD and BETD are not assumed to be temperature dependent and are always evaluated at T = 0.0.
Special-purpose elements, such as MATRIX27 and FLUID29, generally do not require damping. However,
if material property ALPD and BETD are specified for these elements, the value will be used to create
the damping matrix at solution time.
See Damping Matrices in the Mechanical APDL Theory Reference for more details about the damping
formulation.
See Damping in the Structural Analysis Guide for more information about DMPR.
The effects of thermal expansion can be accounted for in three different (and mutually exclusive) ways:
When you use ALPX to enter values for the secant coefficient of thermal expansion (αse), the program
interprets those values as secant or mean values, taken with respect to some common datum or defin-
ition temperature. For instance, suppose you measured thermal strains in a test laboratory, starting at
23°C, and took readings at 200°, 400°, 600°, 800°, and 1000°. When you plot this strain-temperature
data, you could input this directly using THSX. The slopes of the secants to the strain-temperature curve
would be the mean (or secant) values of the coefficient of thermal expansion, defined with respect to
the common temperature of 23° (To). You can also input the instantaneous coefficient of thermal expan-
sion (αin, using CTEX). The slopes of the tangents to this curve represent the instantaneous values.
Hence, the figure below shows how the alternate ways of inputting coefficients of thermal expansion
relate to each other.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
18 of ANSYS, Inc. and its subsidiaries and affiliates.
Linear Material Properties
εth
αin
αse
To Tn T
where:
If the material property data is in terms of instantaneous values of α, then the program will convert
those instantaneous values into secant values as follows:
∫α
α =
−
where:
If the material property data is in terms of thermal strain, the program will convert those strains into
secant values of coefficients of thermal expansion as follows:
ε
α =
− rf
where:
If necessary, the data is shifted so that the thermal strain is zero when Tn = Tref. If a data point at Tref
exists, a discontinuity in αse may be generated at Tn = Tref. This can be avoided by ensuring that the
slopes of εith on both sides of Tref match.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 19
Chapter 3: Material Models
If the αse values are based upon a definition temperature other than TREF, then you need to convert
those values to TREF. This can be done using the MPAMOD command. Also see the Mechanical APDL
Theory Reference.
3.3.5. Emissivity
EMIS defaults to 1.0 if not defined. However, if EMIS is set to zero or blank, EMIS is taken to be 0.0.
• Heat-transfer elements: All properties are evaluated at the centroid (except for the specific heat or
enthalpy, which is evaluated at the integration points).
• Structural elements: All properties are evaluated at the integration points.
• All other elements: All properties are evaluated at the centroid.
If the temperature of the centroid or integration point falls below or rises above the defined temperature
range of tabular data, ANSYS assumes the defined extreme minimum or maximum value, respectively,
for the material property outside the defined range.
For a list of the elements that support this material model, see "Element Support for Material Models".
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
20 of ANSYS, Inc. and its subsidiaries and affiliates.
Plasticity Models
Select the material behavior option via menu path Main Menu> Preprocessor> Material Props> Ma-
terial Models, or via a TB,Lab command, as follows:
Most options require a uniaxial stress-strain curve to be input. The exceptions are CHABOCHE, NLISO,
HILL, DP, and USER.
Most options must have elastically isotropic (EX = EY = EZ) materials. The exceptions are HILL, ANISO,
and USER.
Required values not included in the data table are assumed to be zero. If the data table is not defined
(or contains all zero values), the material is assumed to be linear.
The material behavior options are described below. See the Mechanical APDL Theory Reference for more
detail.
3.4.1. Bilinear Kinematic Hardening
3.4.2. Multilinear Kinematic Hardening
3.4.3. Nonlinear Kinematic Hardening
3.4.4. Bilinear Isotropic Hardening
3.4.5. Multilinear Isotropic Hardening
3.4.6. Nonlinear Isotropic Hardening
3.4.7. Anisotropic
3.4.8. Hill's Anisotropy
3.4.9. Drucker-Prager
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 21
Chapter 3: Material Models
Initialize the stress-strain table with TB,BKIN. For each stress-strain curve, define the temperature (TB-
TEMP), then define C1 and C2 (TBDATA). You can define up to six temperature-dependent stress-strain
curves (NTEMP = 6 maximum on the TB command) in this manner. The constants C1 and C2 are:
BKIN can be used with the TBOPT option. In this case, TBOPT takes two arguments. For TB,BKIN,,,,0,
there is no stress relaxation with an increase in temperature. This option is not recommended for non-
isothermal problems. For TB,BKIN,,,,1, Rice's hardening rule is applied (which does take stress relaxation
with temperature increase into account). This is the default.
See Plastic Material Options in the Structural Analysis Guide for more information on this material option.
You can combine this option with other material options to simulate more complex material behaviors.
See Material Model Combinations (p. 151) for further information.
The two options use the Besseling model (see the Mechanical APDL Theory Reference), also called the
sublayer or overlay model. The material response is represented by multiple layers of perfectly plastic
material; the total response is obtained by weighted average behavior of all the layers. Individual weights
are derived from the uniaxial stress-strain curve. The uniaxial behavior is described by a piece-wise linear
"total stress-total strain curve", starting at the origin, with positive stress and strain values. The slope
of the first segment of the curve must correspond to the elastic modulus of the material and no segment
slope should be larger. The slope of the stress-strain curve is assumed to be zero beyond the last user-
defined stress-strain data point. In the following, the option TB,KINH is described first, followed by that
of TB,MKIN.
The KINH option is recommended because layers are scaled (Rice's model), providing better represent-
ations. The KINH option allows you to define up to 40 temperature-dependent stress-strain curves. If
you define more than one stress-strain curve for temperature-dependent properties, then each curve
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
22 of ANSYS, Inc. and its subsidiaries and affiliates.
Plasticity Models
should contain the same number of points (up to a maximum of 20 points in each curve). The assumption
is that the corresponding points on the different stress-strain curves represent the temperature dependent
yield behavior of a particular sublayer.
For stress vs. total strain input, initialize the stress-strain table with TB,KINH. For stress vs. plastic strain
input, initialize the stress-strain table with either TB,KINH,,,,PLASTIC or TB,PLASTIC,,,,KINH. Input the
temperature of the first curve with the TBTEMP, then input stress and strain values using the TBPT.
Input the remaining temperatures and stress-strain values using the same sequence (TBTEMP followed
by TBPT).
See Plastic Material Options in the Structural Analysis Guide for more information on this material option.
You can combine this option with other material options to simulate more complex material behaviors.
See Material Model Combinations (p. 151) for further information.
The curve defined with the MKIN option is continuous from the origin with a maximum of five total
stress-total strain points. The slope of the first segment of the curve must correspond to the elastic
modulus of the material and no segment slope should be larger.
Initialize the stress-strain table with TB,MKIN, followed by a special form of the TBTEMP command
(TBTEMP,,STRAIN) to indicate that strains are defined next. The constants (C1-C5), entered on the next
TBDATA command, are the five corresponding strain values (the origin strain is not input). The temper-
ature of the first curve is then input with TBTEMP, followed by the TBDATA command with the constants
C1-C5 representing the five stresses corresponding to the strains at that temperature. You can define
up to five temperature-dependent stress-strain curves (NTEMP = 5 max on the TB command) with the
TBTEMP command.
MKIN can also be used in conjunction with the TBOPT option (TB,MKIN,,,,TBOPT). TBOPT has the fol-
lowing three valid arguments:
0 - No stress relaxation with temperature increase (this is not recommended for nonisothermal
problems); also produces thermal ratcheting.
1 - Recalculate total plastic strain using new weight factors of the subvolume.
2 - Scale layer plastic strains to keep total plastic strain constant; agrees with Rice's model (TB, BKIN
with TBOPT = 1). Produces stable stress-strain cycles.
Note
The mechanical behavior of the TB,KINH option is the same as TB,MKIN with TBOPT =
2.
See Plastic Material Options in the Structural Analysis Guide for more information on this material option.
You can combine this option with other material options to simulate more complex material behaviors.
See Material Model Combinations (p. 151) for further information.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 23
Chapter 3: Material Models
n n pl pl
αɺ = ∑ αɺ i = ∑ iεɺ − γiαiεɺ + i θα
ɺ
i
i i i θ
where:
σ ε = σ − =
where:
Initialize the data table with TB,CHABOCHE. For each set of data, define the temperature (TBTEMP),
then define C1 through Cm (TBDATA), where m = 1 + 2NPTS. The maximum number of constants, m is
11, which corresponds to 5 kinematic models (NPTS = 5 on the TB command). The default value for m
is 3, which corresponds to one kinematic model (NPTS = 1). You can define up to 1000 temperature-
dependent constants ([NTEMP x m ≤ 1000] maximum on the TB command) in this manner. The constants
C1 through C(1 + 2NPTS) are:
Constant Meaning
C1 k = Yield stress
C2 C1 = Material constant for first kinematic model
C3 γ1 = Material constant for first kinematic model
C4 C2 = Material constant for second kinematic
model
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
24 of ANSYS, Inc. and its subsidiaries and affiliates.
Plasticity Models
Constant Meaning
C5 γ2 = Material constant for second kinematic
model
... ...
C(2NPTS) CNPTS = Material constant for last kinematic
model
C(1 + γNPTS = Material constant for last kinematic
2NPTS) model
k, and all C and γ values in the right column are material constants in the Chaboche model (see the
Mechanical APDL Theory Reference for details).
See Plastic Material Options in the Structural Analysis Guide for more information on this material option.
You can combine this option with other material options to simulate more complex material behaviors.
See Material Model Combinations (p. 151) for further information.
Initialize the stress-strain table with TB,BISO. For each stress-strain curve, define the temperature (TB-
TEMP), then define C1 and C2 (TBDATA). Define up to six temperature-dependent stress-strain curves
(NTEMP = 6 max on the TB command) in this manner. The constants C1 and C2 are:
See Plastic Material Options in the Structural Analysis Guide for more information on this material option.
You can combine this option with other material options to simulate more complex material behaviors.
See Material Model Combinations (p. 151) for further information.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 25
Chapter 3: Material Models
You can specify up to 20 temperature-dependent stress-strain curves. For stress vs. total strain input,
initialize the curves with TB,MISO. For stress vs. plastic strain input, initialize the curves with
TB,PLASTIC,,,,MISO. Input the temperature for the first curve (TBTEMP), followed by up to 100 stress-
strain points (the origin stress-strain point is not input) (TBPT). Define up to 20 temperature-dependent
stress-strain curves (NTEMP = 20, maximum on the TB command) in this manner. The constants (X, Y)
entered on the TBPT command (two per command) are:
See Plastic Material Options in the Structural Analysis Guide for more information on this material option.
You can combine this option with other material options to simulate more complex material behaviors.
See Material Model Combinations (p. 151) for further information.
Once you initialize the data table with TB,NLISO, you define the temperature via the TBTEMP command,
and then define either C1 through C4 (TBOPT = VOCE) or C1 and C2 (TBOPT = POWER) via the TBDATA
command. You can define up to twenty temperature-dependent stress-strain curves (NTEMP = 20
maximum on the TB command) in this manner.
See Plastic Material Options in the Structural Analysis Guide for more information on this material option.
3.4.7. Anisotropic
This option (TB,ANISO) allows for different stress-strain behavior in the material x, y, and z directions
as well as different behavior in tension and compression (see Anisotropic Elastic Material Model (p. 69)).
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
26 of ANSYS, Inc. and its subsidiaries and affiliates.
Plasticity Models
A modified von Mises yield criterion is used to determine yielding. The theory is an extension of Hill's
formulation as noted in the Mechanical APDL Theory Reference. This option is not recommended for
cyclic or highly nonproportional load histories since work hardening is assumed. The principal axes of
anisotropy coincide with the material (or element) coordinate system and are assumed not to change
over the load history.
The material behavior is described by the uniaxial tensile and compressive stress-strain curves in three
orthogonal directions and the shear stress-engineering shear strain curves in the corresponding directions.
A bilinear response in each direction is assumed. The initial slope of the curve is taken as the elastic
moduli of the material. At the specified yield stress, the curve continues along the second slope defined
by the tangent modulus (having the same units as the elastic modulus). The tangent modulus cannot
be less than zero or greater than the elastic modulus. Temperature dependent curves cannot be input.
All values must be input as no defaults are defined. Input the magnitude of the yield stresses (without
signs). No yield stress can have a zero value. The tensile x-direction is used as the reference curve for
output quantities SEPL and EPEQ.
Initialize the stress-strain table with TB,ANISO. You can define up to 18 constants with TBDATA com-
mands. The constants (C1-C18) entered on TBDATA commands (6 per command) are:
See Plastic Material Options in the Structural Analysis Guide for more information on this material option.
• Rate-independent anisotropic plasticity with isotropic hardening -- TB,HILL combined with TB,BISO
or TB,MISO or TB,NLISO.
• Rate-independent anisotropic plasticity with kinematic hardening -- TB,HILL combined with TB,BKIN
or TB,MKIN or TB,KINH or TB,CHAB.
• Rate-independent anisotropic plasticity with combined hardening -- TB,HILL combined with TB,CHAB
and TB,BISO or TB,MISO or TB,NLISO.
• Rate-dependent anisotropic plasticity (anisotropic viscoplasticity) with isotropic hardening -- TB,HILL
combined with TB,BISO or TB,MISO or TB,NLISO and TB,RATE.
• Anisotropic creep -- TB,HILL combined with TB,CREEP (implicit).
• Anisotropic creep and anisotropic plasticity with isotropic hardening -- TB,HILL combined with
TB,CREEP and TB,BISO or TB,MISO or TB,NLISO (implicit).
• Anisotropic creep and anisotropic plasticity with kinematic hardening -- TB,HILL combined with
TB,CREEP and TB,BKIN (implicit)
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 27
Chapter 3: Material Models
See Material Model Combinations (p. 151) for more information on combining the HILL option with the
plasticity and creep options.
The HILL option's material behavior is described by six constants that define the stress ratios in different
directions (see the Mechanical APDL Theory Reference for details).
Initialize the data table with TB,HILL. For each set of data, you then define the temperatures using the
TBTEMP command, then define C1 through C6 using the TBDATA command. The input must then be
followed by the TB command again, but with one of the plasticity and / or creep options.
For each set of data, you then define the temperature using the TBTEMP command, and then define
the constants using the TBDATA command.
For plasticity, rij is the ratio of the yield stress in the ij direction, to the yield stress specified for the
plasticity input as part of the TB command.
For creep, rij is the ratio of the creep strain in the ij direction to the reference value calculated by the
implicit creep equation.
3.4.9. Drucker-Prager
This option (TB,DP) is applicable to granular (frictional) material such as soils, rock, and concrete and
uses the outer cone approximation to the Mohr-Coulomb law (see the Mechanical APDL Theory Reference).
The input consists of only three constants:
The amount of dilatancy (the increase in material volume due to yielding) can be controlled with the
dilatancy angle. If the dilatancy angle is equal to the friction angle, the flow rule is associative. If the
dilatancy angle is zero (or less than the friction angle), there is no (or less of an) increase in material
volume when yielding and the flow rule is nonassociated. Temperature-dependent curves are not allowed.
Initialize the constant table with TB,DP. You can define up to three constants with TBDATA commands.
The constants (C1-C3) entered on TBDATA are:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
28 of ANSYS, Inc. and its subsidiaries and affiliates.
Plasticity Models
See Plastic Material Options in the Structural Analysis Guide for more information on this material option.
where:
where:
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 29
Chapter 3: Material Models
You can use any combination of the yield and potential functions listed (above), but one of each must
be specified. The cap model, however, is a special EDP case that requires exclusive use of combined
cap potential and yield functions. More information on the EDP cap model is provided in the following
section.
When plasticity is defined by TB,MISO, TB,BISO, TB,NLISO, or TB,PLAS, that definition overrides the yield
stress definition you define using TB,EDP and TBDATA.
See the EDP argument and associated specifications in the TB command, and also The Extended
Drucker-Prager Model in the Mechanical APDL Theory Reference for more information.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
30 of ANSYS, Inc. and its subsidiaries and affiliates.
Plasticity Models
C10
D
C11
2
See Cap Model in the Mechanical APDL Theory Reference for more information on the constant conven-
tions.
You can input the cohesion-related shear hardening constant σi by using TB,MISO, TB,BISO, TB,NLISO,
or TB,PLAS. This input regulates the relationship between the cohesion and the effective deviatoric
plastic strain.
Note
Calibrating the cap constants σi, βY, A, αY, βF, αF and the hardening input for σi differs signi-
ficantly from the other EDP options. The cap parameters are all defined in relation to I1 and
J2, while the other EDP coefficients are defined according to p and q.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 31
Chapter 3: Material Models
You can combine certain TBOPT values to address specific model behaviors. The following combinations
are valid:
2
Although all three constants (q1, q2, q3) are independent input parameters, setting q3 = q1 is recom-
mended, unless specific data is available.
Isotropic hardening can be modeled using TB, BISO, TB, MISO, TB, NLISO or TB, PLASTIC.
The material behavior is described by a piece-wise linear stress-strain curve, starting at the origin, with
positive stress and strain values. The curve is continuous from the origin through 100 (max) stress-strain
points. Successive slopes can be greater than the preceding slope; however, no slope can be greater
than the elastic modulus of the material. The slope of the first curve segment usually corresponds to
the elastic modulus of the material, although the elastic modulus can be input as greater than the first
slope to ensure that all slopes are less than or equal to the elastic modulus.
Specify up to 20 temperature-dependent stress-strain curves. Initialize the curves with TB,MELAS. The
temperature for the first curve is input with TBTEMP, followed by TBPT commands for up to 100 stress-
strain points (the origin stress-strain point is not input). You can define up to 20 temperature- dependent
stress-strain curves (NTEMP = 20 max on the TB command) in this manner. The constants (X, Y) entered
on TBPT (two per command) are:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
32 of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelastic Material Models
See Multilinear Elasticity in the Structural Analysis Guide for more information on this material option.
Initiate the cast iron material model with TB,CAST. Activate the stress-strain table in tension using
TB,UNIAXIAL with the TENSION option, then enter the stress-strain relation using the TBPT command.
Activate the stress-strain table in compression using theTB, UNIAXIAL with the COMPRESSION option,
then enter the stress-strain relation using the TBPT command. The slope of the stress-strain curve is
assumed to be zero beyond the last user-defined stress-strain data point.
The NROPT,UNSYM command should be used at the solution level as the flow rule is not associated
and the material Jacobian matrix is unsymmetric.
Initialize the database with TB,CAST. For each set of data, define the temperature using TBTEMP, then
define the constant C1.
See Plastic Material Options in the Structural Analysis Guide for more information on this material option.
Hyperelasticity options are available via the TBOPT argument on the TB,HYPER command.
Several forms of strain energy potentials describe the hyperelasticity of materials. These are based on
either strain invariants or principal stretches. The behavior of materials is assumed to be incompressible
or nearly incompressible.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 33
Chapter 3: Material Models
For information about other hyperelastic material models, see Special Hyperelastic Material Models (p. 43).
2 3
= µ 1− + 1 − + 1 −
λL2 λL4
2−
+ 4 5 +
1 − + 1 − −
λL6 λL8
where:
The constants µ, λL and d are defined by C1, C2, and C3 using the TBDATA command.
For a list of elements that can be used with this material option, see "Element Support for Material
Models".
See Arruda-Boyce Hyperelastic Option in the Structural Analysis Guide for more information on this ma-
terial option.
µ
= + −
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
34 of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelastic Material Models
= µ
The model has only one constant µ and is defined by C1 using the TBDATA command.
See Blatz-Ko Foam Hyperelastic Option in the Structural Analysis Guide for more information on this
material option.
Following the material data table command (TB), specify the material constant values via the TBDATA
command , as shown in this example:
TB,HYPER,1,,5,ETUBE ! Hyperelastic material, 1 temperature,
! 5 material constants, and the extended tube option
TBDATA,1,0.25, 0.8,1.0,0.5,1.0e-5 ! Five material constant values (C1 through C5)
For more information, see the documentation for the TB,HYPER command, and Extended Tube Model
in the Mechanical APDL Theory Reference.
2−
=−
µ m − 1− + −
m
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 35
Chapter 3: Material Models
The constants µ, Jm, and d are defined by C1, C2, and C3 using the TBDATA command.
For a list of elements that can be used with this material option, see "Element Support for Material
Models".
See Gent Hyperelastic Option in the Structural Analysis Guide for more information on this material option.
The TB,HYPER,,,,MOONEY option allows you to define 2, 3, 5, or 9 parameter Mooney-Rivlin models using
NPTS = 2, 3, 5, or 9, respectively.
For NPTS = 2 (2 parameter Mooney-Rivlin option, which is also the default), the form of the strain energy
potential is:
= 0 − + 0 2− + − 2
where:
µ= +
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
36 of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelastic Material Models
The constants c10, c01, and d are defined by C1, C2, and C3 using the TBDATA command.
For NPTS = 3 (3 parameter Mooney-Rivlin option, which is also the default), the form of the strain energy
potential is:
= 10 1 − + 01 2 − + 11 1 − 2− + − 2
The constants c10, c01, c11; and d are defined by C1, C2, C3, and C4 using the TBDATA command.
For NPTS = 5 (5 parameter Mooney-Rivlin option), the form of the strain energy potential is:
= − +
− + −
+ − + −
− + −
The constants c10, c01, c20, c11, c02, and d are material constants defined by C1, C2, C3, C4, C5, and C6
using the TBDATA command.
For NPTS = 9 (9 parameter Mooney-Rivlin option), the form of the strain energy potential is:
= − + − + −
+ − + 3
− + − 3 −
+ − + 3+ −
− + − − 3 −
The constants c10, c01, c20, c11, c02, c30, c21, c12, c03, and d are material constants defined by C1, C2, C3,
C4, C5, C6, C7, C8, C9, and C10 using the TBDATA command.
See Mooney-Rivlin Hyperelastic Option (TB,HYPER) in the Structural Analysis Guide for more information
on this material option.
µ
= − + −
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 37
Chapter 3: Material Models
See Neo-Hookean Hyperelastic Option in the Structural Analysis Guide for more information on this
material option.
N
µ α α α α
N
µ −α β
=∑ i / 3
λ +λ +λ − +∑ i
−
i =
1
α i
1 2 3
i =
1
α β
i i
where:
For this material option, the volumetric and deviatoric terms are tightly coupled. Hence, this model is
meant to simulate highly compressible elastomers.
In general there is no limitation on the value of N. (See the TB command.) A higher value of N can
provide a better fit to the exact solution. It may however cause numerical difficulties in fitting the ma-
terial constants. For this reason, very high values of N are not recommended.
∑µα
µ= =
= ∑µα
+β
=
For N = 1, α1 = –2, µ1 = -µ, and β1 = 0.5, the Ogden foam option is equivalent to the Blatz-Ko option.
The constants µi, αi and βi are defined using the TBDATA command in the following order:
For N (NPTS) = 1:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
38 of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelastic Material Models
µ1, α1, β1
For N (NPTS) = 2:
For N (NPTS) = 3:
For N (NPTS) = k:
µ1, α1, µ2, α2, ..., µk, αk, β1, β2, ..., βk
See Ogden Compressible Foam Hyperelastic Option in the Structural Analysis Guide for more information
on this material option.
µ α
λ + λα + λα −
N N
=∑ i
+ ∑ − 2k
i = α 1 i
1 2 3
k =
1 k
where:
λ p
( = ,,) = devatorc rncal stretches, defned as λ =
p
λ p
In general there is no limitation on the value of N. (See the TB command.) A higher value of N can
provide a better fit to the exact solution. It may however cause numerical difficulties in fitting the ma-
terial constants. For this reason, very high values of N are not recommended.
µ= ∑α µ
=
For N = 1 and α1 = 2, the Ogden option is equivalent to the Neo-Hookean option. For N = 2, α1 = 2,
and α2 = -2, the Ogden option is equivalent to the 2 parameter Mooney-Rivlin option.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 39
Chapter 3: Material Models
The constants µp, αp and dp are defined using the TBDATA command in the following order:
For N (NPTS) = 1:
µ1, α1, d1
For N (NPTS) = 2:
For N (NPTS) = 3:
For N (NPTS) = k:
µ1, α1, µ2, α2, ..., µk, αk, d1, d2, ..., dk
See Ogden Hyperelastic Option in the Structural Analysis Guide for more information on this material
option.
N N
= ∑ i j+ − 2k
ij 1 − 2− ∑
i + j =1 k =1 k
where:
=
J = determinant of the elastic deformation gradient F
N, cij, and d = material constants
In general there is no limitation on the value of N. (See the TB command.) A higher value of N can
provide a better fit to the exact solution. It may however cause a numerical difficulty in fitting the ma-
terial constants, and it also requests enough data to cover the whole range of deformation for which
you may be interested. For these reasons, a very high value of N is not recommended.
µ= 0 + 0
=
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
40 of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelastic Material Models
For N = 1 and c01 = 0, the polynomial form option is equivalent to the Neo-Hookean option. For N =
1, it is equivalent to the 2 parameter Mooney-Rivlin option. For N = 2, it is equivalent to the 5 parameter
Mooney-Rivlin option, and for N = 3, it is equivalent to the 9 parameter Mooney-Rivlin option.
The constants cij and d are defined using the TBDATA command in the following order:
For N (NPTS) = 1:
c10, c01, d1
For N (NPTS) = 2:
For N (NPTS) = 3:
c10, c01, c20, c11, c02, c30, c21, c12, c03, d1, d2, d3
For N (NPTS) = k:
c10, c01, c20, c11, c02, c30, c21, c12, c03, ..., ck0, c(k-1)1, ..., c0k, d1, d2, ..., dk
See Polynomial Form Hyperelastic Option in the Structural Analysis Guide for more information on this
material option.
The response functions (first derivatives of the hyperelastic potential) are used to determine hyperelastic
constitutive behavior of the material. In general, the stiffness matrix requires derivatives of the response
functions (second derivatives of the hyperelastic potential).
The method for determining the derivatives is ill-conditioned near the zero stress-strain point; therefore,
a deformation limit is used, below which the stiffness matrix is calculated with only the response func-
tions. The deformation measure is δ = I1 - 3, where I1 is the first invariant of the Cauchy-Green deform-
ation tensor.
The stiffness matrix is then calculated with only the response functions if δ < C1, where C1 is the mater-
ial constant deformation limit (default 1 x 10-5).
The remaining material parameters are for the volumetric strain energy potential, given by
N
= ∑ ( − )2k
k =1 k
where N is the NPTS value (TB,HYPER,,,,RESPONSE) and dk represents the material constants incompress-
ibility parameters (default 0.0) and J is the volume ratio. Use of experimental volumetric data requires
NPTS = 0. Incompressible behavior results if all dk = 0 or NPTS = 0 with no experimental volumetric
data.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 41
Chapter 3: Material Models
N N
= ∑ i0 i+ − 2k
1− ∑
i =1 k =1 k
where:
In general there is no limitation on the value of N. (See the TB command.) A higher value of N can
provide a better fit to the exact solution. It may however cause a numerical difficulty in fitting the ma-
terial constants, and it also requests enough data to cover the whole range of deformation for which
you may be interested. For these reasons, a very high value of N is not recommended.
µ=
=
The constants ci0 and dk are defined using the TBDATA command in the following order:
For N (NPTS) = 1:
c10, d1
For N (NPTS) = 2:
For N (NPTS) = 3:
For N (NPTS) = k:
See Yeoh Hyperelastic Option in the Structural Analysis Guide for more information on this material option.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
42 of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelastic Material Models
Two strain energy potentials, as forms of polynomial or exponential function, are available for charac-
terizing the isochoric part of strain energy potential.
You can use anisotropic hyperelasticity to model elastomers with reinforcements, and for biomedical
materials such as muscles or arteries.
= v + d ⊗ ⊗
= ⋅ − 2
3 3 6
i j
⊗ ⊗ = ∑ i 1− + ∑ j − + ∑ k 4− k
i =1 j =1 k =
6 6 6 6
+∑ l 5 − l+ ∑ m 6 − m+ ∑ n 7− n+ ∑ o 8 −ς o
l= m= n= o =
( ⊗ ⊗ ) = ∑= ( − ) + ∑= ( − ) +
(
− )
− +
(
− ) −
Use TB,AHYPER,,TBOPT to define the isochoric part, material directions and the volumetric part. Only
one TB table can be defined for each option. You can either define polynomial or exponential strain
energy potential.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 43
Chapter 3: Material Models
You can enter temperature-dependent data for anisotropic hyperelastic material via the TBTEMP com-
mand. For the first temperature curve, issue TB, AHYPER,,,TBOPT, then input the first temperature (TB-
TEMP). The subsequent TBDATA command inputs the data.
The program interpolates the temperature data to the material points automatically using linear inter-
polation. When the temperature is out of the specified range, the closest temperature point is used.
For more information, see the TB command, and Anisotropic Hyperelasticity in the Mechanical APDL
Theory Reference.
This model requires seven material constants input for the isochoric (TBOPT = ISO) option and one
material constant for the volumetric potential (TBOPT = PVOL) option. Issue the TBDATA data table
command to input the constant values in the order shown:
Isochoric
TB,BB,,,,ISO
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
44 of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelastic Material Models
The default optional material constant is ε = 1 x 10-5. However, if TBNPT > 7 or TBNPT is unspecified,
the table value is used instead. If the table value is zero or exceeds 1 x 10-3, the default constant value
is used.
Volumetric Potential
TB,BB,,,,PVOL
TB,CDM,MAT,NTEMPS,NPTS,TBOPT
TBOPT = PSE2
Constant Meaning Property
C1 r Damage variable para-
meter
C2 m Damage variable para-
meter
C3 β Damage variable para-
meter
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 45
Chapter 3: Material Models
See User-Defined Hyperelastic Option (TB,HYPER,,,,USER) in the Structural Analysis Guide for more inform-
ation on this material option.
The elastic behavior is rate-independent and represents the recoverable deformation due to mechanical
loading. The viscous behavior is rate-dependent and represents dissipative mechanisms within the
material.
A wide range of materials (such as polymers, glassy materials, soils, biologic tissue, and textiles) exhibit
viscoelastic behavior.
Following are descriptions of the viscoelastic constitutive models, which include both small- and large-
deformation formulations. Also presented is time-temperature superposition for thermorheologically
simple materials and a harmonic domain viscoelastic model.
3.6.1. Viscoelastic Formulation
3.6.2.Time-Temperature Superposition
3.6.3. Harmonic Viscoelasticity
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
46 of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelastic Material Model
The spring stiffnesses are µi, the dashpot viscosities are ηi , and the relaxation time is defined as the
ratio of viscosity to stiffness, τi = ηi / µi.
In three dimensions, the constitutive model for a generalized Maxwell model is given by:
t t ∆
σ= ∫ −τ τ+ ∫ −τ τ (3–1)
0 τ 0 τ
where:
σ = Cauchy stress
e = deviatoric strain
∆ = volumetric strain
τ = past time
I = identity tensor
and G(t) and K(t) are the Prony series shear and bulk-relaxation moduli, respectively:
n G
= αG∞ + ∑ αi −
i =1 τG (3–2)
i
K
= αK∞ + ∑ α −
= τK (3–3)
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 47
Chapter 3: Material Models
For use in the incremental finite element procedure, the solution for Equation 3–1 (p. 47) at t1 = t0 +
∆t is:
∆ t −τ
i 1 = i( 0) − + αG − 1 τ
τG t∫ 0 i τG τ (3–4)
i i
∆ −τ ∆
= ( ) − + ∫ αK
− τ (3–5)
τK τK τ
where si and pi are the deviatoric and pressure components, respectively, of the Cauchy stress for each
Maxwell element.
∆
∆
= ( ) − +
τ α − (
τ
− ) (3–6)
∆ ∆
= ( ) − + α
− (∆ − ∆ ) (3–7)
τ
τ
The model requires input of the parameters in Equation 3–2 (p. 47) and Equation 3–3 (p. 47). The relax-
ation moduli at t = 0 are obtained from the elasticity parameters input using the MP command or via
an elastic data table (TB,ELASTIC). The Prony series relative moduli and relaxation times are input via a
Prony data table (TB,PRONY), and separate data tables are necessary for specifying the bulk and shear
Prony parameters.
For the shear Prony data table, TBOPT = SHEAR, NPTS = nG, and the constants in the data table follow
this pattern:
For the bulk Prony data table, TBOPT = BULK, NPTS = nK, and the constants in the data table follow
this pattern:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
48 of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelastic Material Model
To account for large displacement, the model is formulated in the co-rotated configuration using the
co-rotated deviatoric stress Σ = RTsR, where R is the rotation obtained from the polar decomposition
of the deformation gradient. The pressure component of the Cauchy stress does not need to account
for the material rotation and uses the same formulation as the small-deformation model.
∆ ∆
i = ∆ Σi ( 0 ) ∆
T − +
τG
G
0 αi −
τG ( 1 −∆ 0 ∆
T
) (3–8)
i i
Parameter input for this model resembles the input requirements for the small-deformation viscoelastic
model.
The linear structure of the formulation is provided by the generalized Maxwell model. Extension to
large-deformation requires only a hyperelastic model for the springs in the Maxwell elements. Hypere-
lasticity is defined by a strain energy potential where, for isotropic materials:
= + (3–9)
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 49
Chapter 3: Material Models
∂
= (3–10)
∂
i
and the large-deformation stress update for the Maxwell element stresses is given by:
′ ′
∆ ∆
( )=1
( ) 0
−
τ
+α
G
−
τ
− (3–11)
G
G
t
t
∆ ∆
( )= ( )
−
τ
+α
K
−
τ
− (3–12)
K
K
where:
′
=
deviatoric component of Si
=
pressure component of Si
An anisotropic hyperelastic model can also be used for Equation 3–9 (p. 49) , in which case the form of
the Maxwell element stress updates are unchanged.
This model requires the Prony series parameters to be input via the Prony data table (as described in
Small Deformation (p. 46)). The hyperelastic parameters for this model are input via a hyperelastic data
table (TB,HYPER). For more information, see Hyperelastic Material Models (p. 33).
= r
ξ (3–13)
where:
T = current temperature
Tr = constant reference temperature
ξ = shifted time given by ξ = t / A(T), where A(T) = shift function.
The constitutive equations are solved in the shifted time scale. This method has the potential to reduce
the experimental effort required to determine the material parameters but requires the determination
of the shift function. The following forms of the shift function are available:
3.6.2.1. Williams-Landel-Ferry Shift Function
3.6.2.2.Tool-Narayanaswamy Shift Function
3.6.2.3. User-Defined Shift Function
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
50 of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelastic Material Model
( − )
( )= 1 r
(3–14)
+( − )
10
2 r
where C1 and C2 are material parameters. (The shift function is often given in the literature with the
opposite sign.)
The parameters are input via a shift function data table (TB,SHIFT).
For the Williams-Landel-Ferry shift function, TBOPT = WLF, and the required input constants are:
( )= − (3–15)
For the Tool-Narayanaswamy shift function, TBOPT = TN, and the required input constants are:
The second form of the Tool-Narayanaswamy shift function includes an evolving fictive temperature.
The fictive temperature is used to model material processes that contain an intrinsic equilibrium tem-
perature that is different from the ambient temperature of the material. The shift function is given by:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 51
Chapter 3: Material Models
−
( )= − − (3–16)
r F
where:
X = weight parameter
TF = fictive temperature.
n
= ∑ fi fi (3–17)
i =1
where:
τ 0 + ∆ 0
= 0 (3–18)
τ + ∆
where:
τ =
fictive temperature relaxation time
0 (superscript) = values from the previous time step
The fictive temperature model modifies the volumetric thermal strain model and gives an incremental
thermal strain as:
where αg and αl are the glass and liquid coefficients, respectively, of thermal expansion given by:
α = α + α
+ α 2 2 + α 3 3 + α4 4 (3–20)
α = α + α + α + α + α (3–21)
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
52 of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelastic Material Model
For the Tool-Narayanaswamy with fictive temperature shift function, TBOPT = FICT, NPTS = nf, and the
required input constants are:
Given the input parameters, the routine must evolve the internal state variables, then return the current
and half-step shifted time.
Assuming that the strain varies harmonically and that all transient effects have subsided, Equa-
tion 3–1 (p. 47) has the form:
′ ′′ ′ ′′
σ= + Ω +δ + + ∆ Ω +δ (3–22)
where:
Comparing Equation 3–22 (p. 53) to the harmonic equation of motion, the material stiffness is due to
the storage moduli and the material damping matrix is due to the loss moduli divided by the frequency.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 53
Chapter 3: Material Models
n αG τGΩ 2 n αG τGΩ
′
= 0 − ∑ αiG + i i ′′
= 0∑ i i
+ τiGΩ 2 i =1 + τiGΩ 2
(3–23)
i =1
αK K
τ Ω
αK τK Ω
′
= − ∑ αK + ′′
= ∑
+ τK = + τK Ω
(3–24)
= Ω
Input of the Prony series parameters for a viscoelastic material in harmonic analyses follows the input
method for viscoelasticity in the time domain detailed above.
The points for the experimental data table (input via the TBPT command) have frequency as the inde-
pendent variable, and the dependent variables are the real component, imaginary component, and
tan(δ). If the imaginary component is empty or zero for the data point, the tan(δ) value is used to de-
termine it; otherwise tan(δ) is not used.
Input complex shear modulus on an experimental data table (TB,EXPE) with TBOPT = GMODULUS. The
data points are defined by:
Position Value
1 Ω
2 ′
3 ′′
4 ′′
δ =
′
Input complex bulk modulus on an experimental data table (TB,EXPE) with TBOPT = KMODULUS. The
data points are defined by:
Position Value
1 Ω
2 ′
3 ′′
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
54 of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelastic Material Model
Position Value
4 ′′
δ =
′
Input complex tensile modulus on an experimental data table (TB,EXPE) with TBOPT = EMODULUS.
The data points are defined by:
Position Value
1 Ω
2 ′
3 ′′
4 ′′
δ =
′
Input complex Poisson's ratio on an experimental data table (TB,EXPE) with TBOPT = NUXY. The data
points are defined by:
Position Value
1 Ω
2
ν′
3
ν′′
4
ν′′
δ =
ν′
Using experimental data to define the complex constitutive model requires elastic constants (defined
via MP or by an elastic data table [TB,ELASTIC]). The elastic constants are unused if two sets of complex
modulus experimental data are defined. This model also requires an empty Prony data table (TB,PRONY)
with TBOPT = EXPERIMENTAL.
Two elastic constants are required to define the complex constitutive model. If only one set of experi-
mental data for a complex modulus is defined, the Poisson's ratio (defined via MP or by elastic data
table) is used as the second elastic constant.
Frequency-temperature superposition can be used with either the Prony series complex modulus or
any of the experimental data for complex moduli or Poisson's ratio.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 55
Chapter 3: Material Models
3.6.3.4. Stress
The magnitude of the real and imaginary stress components are obtained from expanding Equa-
tion 3–22 (p. 53) and using the storage and loss moduli from either the Prony series parameters or the
experimental data:
σ = ( ′ ′
) (
− ′′ ′′ + ′ ′
∆ − ′′ ∆′′ ) (3–25)
σ = ( ′ ′′
) (
+ ′′ ′ + ′ ′′
∆ + ′′ ∆′ ) (3–26)
where:
The model is well suited for simulating engineering materials consisting of various aggregate composi-
tions with differing properties (for example, concrete modeling, in which rock and sand are embedded
in a weak matrix of cements).
The following topics concerning the microplane material model are available:
3.7.1. Microplane Modeling
3.7.2. Material Models with Degradation and Damage
3.7.3. Material Parameters Definition and Example Input
3.7.4. Element Support for the Microplane Material Model
3.7.5. Learning More About Microplane Material Modeling
1. Apply a kinematic constraint to relate the macroscopic strain tensors to their microplane counterparts.
2. Define the constitutive laws on the microplane levels, where unidirectional constitutive equations
(such as stress and strain components) are applied on each microplane.
3. Relate the homogenization process on the material point level to derive the overall material response.
(Homogenization is based on the principle of energy equivalence.)
The microplane material model formulation is based on the assumption that microscopic free energy
Ψmic on the microplane level exists and that the integral of Ψmic over all microplanes is equivalent to a
macroscopic free Helmholtz energy Ψmac [3 (p. 62)], expressed as:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
56 of ANSYS, Inc. and its subsidiaries and affiliates.
Microplane Material Model
ψmac = ∫ψ
mic
Ω
πΩ
The factor π results from the integration of the sphere of unit radius with respect to the area Ω.
The strains and stresses at microplanes are additively decomposed into volumetric and deviatoric parts,
respectively, based on the volumetric-deviatoric (V-D) split.
= D + εv
ε = =
where V is the second-order volumetric projection tensor and 1 the second-order identity tensor.
ε = = ⋅∏− ⋅ ⊗ = ⋅ ∏ de
where Π is the fourth-order identity tensor and the vector n describes the normal on the microsphere
(microplane).
ε= ∫ ( ε + T ⋅ ε ) Ω
πΩ
∂ψ
σ= ∫ Ω= ∫ ( σ + ⋅ σ
) Ω
π Ω ∂ε πΩ
where σv and σD are the scalar volumetric stress and the deviatoric stress tensor on the microsphere,
σ = σ σ = σ
and .
∂ψ
= =
∂ε
and
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 57
Chapter 3: Material Models
∂È mic mic
D =
∂εD
= D
where Kmic and Gmic are microplane elasticity parameters and can be interpreted as a sort of microplane
bulk and shear modulus.
The integrals of the macroscopic strain equation and the derived stresses equation are solved via nu-
merical integration:
Np
= ∫ Ω= ∑( )⋅
πΩ π =1
3.7.1.1. Discretization
Discretization is the transfer from the microsphere to microplanes which describe the approximate form
of the sphere. Forty-two microplanes are used for the numerical integration. Due to the symmetry of
the microplanes (where every other plane has the same normal direction), 21 microplanes are considered
and summarized.[3 (p. 62)]
≤
The damage parameter d mic
is the normalized damage variable
{ ≤ }.
The stresses are derived by:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
58 of ANSYS, Inc. and its subsidiaries and affiliates.
Microplane Material Model
∂ mic
σ= ∫
πΩ ∂
Ω= ∫
πΩ
( v+
T⋅
D Ω)
where
( ) (
= −
)
and
( ) (
= −
)
.
ηa = 0 1 + 22
11 + 2 2
where I1 is the first invariant of the strain tensor ε, J2 is the second invariant of the deviatoric part of
the strain tensor ε, and k0, k1, and k2 are material parameters that characterize the form of damage
function. The equivalent strain function implies the Mises-Hencky-Huber criterion for k0 = k1 = 0, and
k2 = 1, and the Drucker-Prager-criterion for k0 > 0, k1 = 0, and k2 = 1.
= − γ o − α
+ α
⋅ (β
( γ
o − η
) )
η
where αmic defines the maximal degradation, βmic determines the rate of damage evolution, and
γ
characterizes the equivalent strain energy on which the material damaging starts (damage starting
boundary).
The following figure shows the evolution of the damage variable d as a function of equivalent strain
energy ηmic for the implemented exponential damage model:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 59
Chapter 3: Material Models
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
60 of ANSYS, Inc. and its subsidiaries and affiliates.
Microplane Material Model
γmic
The material parameters in the model are: E, ν, k0, k1, k2, 0 , and βmic.
E is Young’s modulus and ν is Poisson’s ratio. Both are microplane elastic properties and are defined
via the MP command.
γ
The parameters k0, k1, k2, , and βmic are defined via the TB command (TB,MPLANE).
TB,MPLAN,MAT,NTEMP,NPTS,TBOPT
TBDATA,1,C1,C2,C3,C4,C5,C6
γ
C4 Critical equivalent-strain-energy
density
C5 αmic Maximum damage parameter
mic
C6 β Scale for rate of damage
• PLANE182 and PLANE183 with plane strain and axisymmetric stress states
• SOLID185, SOLID186, and SOLID187
The model can be used with reinforcing elements REINF264 and REINF265 to model reinforced concrete.
1. Bazant, Z. P., P.G. Gambarova . “Crack Shear in Concrete: Crack Band Microplane Model.” Journal of
Structural Engineering . 110 (1984): 2015-2036.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 61
Chapter 3: Material Models
2. Bazant, Z. P., B. H. Oh. “Microplane Model for Progressive Fracture of Concrete and Rock.” Journal for
Engineering Mechanics . 111 (1985): 559-582.
3. Leukart, M., E. Ramm. “A Comparison of Damage Models Formulated on Different Material Scales.”
Computational Materials Science. 28.3-4 (2003): 749-762.
Initialize the curves with the TB,BH command. Use TBPT commands to define up to 500 points (H, B).
The constants (X, Y) entered on TBPT (two per command) are:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
62 of ANSYS, Inc. and its subsidiaries and affiliates.
High-Frequency Electromagnetic Material Models
Specify the system of units (MKS or user defined) with EMUNIT, which also determines the value of the
permeability of free space. This value is used with the relative permeability property values (MP) to es-
tablish absolute permeability values. The defaults (also obtained for Lab = MKS) are MKS units and
free-space permeability of 4 πE-7 Henries/meter. You can specify Lab = MUZRO to define any system
of units, then input free-space permeability.
For more information about this material option, see Additional Guidelines for Defining Regional Mater-
ial Properties and Real Constants in the Low-Frequency Electromagnetic Analysis Guide
For a list of the elements that support 3-D high-frequency material properties, see "Element Support
for Material Models".
11 12 13
= 21
22 23
31 32 33
Constant Meaning
C1-C9 X11, X22, X33, X12, X23, X13, X21, X32, X31
If Xij is 0 where i and j are indexes, then Xji must also be zero.
For TB, DPER and TB,MUR the diagonal elements cannot be zero.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 63
Chapter 3: Material Models
The matrix is give by the following equations for a static magnetic field (Ho) in the x, y, and z directions,
respectively.
µɺ r = µr − κ (Fo H in x-di ection)
κ µr
µ κ
µɺ = y
− κ µ
µ −κ
µɺ = κ
µ z
where:
ωm ω0 + ωα
µ = +
2
ω0 + ωα − ω2
ω!ω
κ=
#
ω" + ωα − ω#
ω$ = γµ$ %
ω& = γµ' s
α = γµ* ∆ ω
γ = γ+
and
Con- Meaning
stant
C1 Saturation magnetization 4πMs (Gauss) (no de-
fault).
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
64 of ANSYS, Inc. and its subsidiaries and affiliates.
High-Frequency Electromagnetic Material Models
Con- Meaning
stant
C2 Lande g-factor (1.8 to 2.5, defaults to 2.0).
C3 Resonance line width ∆H (Oe) (defaults to 0).
C4 Internal dc magnetic field Ho (Oe) (no default).
C5 Direction of Ho.
0 - z direction (default)
1 - y direction (default)
2 - x direction (default)
C6 Sign of off-diagonal element of permeab-
ility matrix.
C6 ≥ 0,
µr −κ
µɺ r = κ
µr (default)
<
µ
κ
µɺ
= − κ
µ
ε′1 − δ1
[ε ] = ε ε′2 − δ2
ε′3 − δ3
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 65
Chapter 3: Material Models
µr′,1 − δm1
[µ ] = µ0 µr′,2 − δm2
µr′,3 − δm3
where:
Constant Meaning
C1 Electric loss tangent in the X direction, tanδ1
C2 Electric loss tangent in the Y direction, tanδ2
C3 Electric loss tangent in the Z direction, tanδ3
C4 Magnetic loss tangent in the X direction,
tanδm1
C5 Magnetic loss tangent in the Y direction,
tanδm2
C6 Magnetic loss tangent in the Z direction,
tanδm3
εs − ε∞
εcoplex = ε∞ +
+ ωτ
where:
τ = relaxation time
ε∞ = optical permittivity
εrs = static permittivity
ω = working angular frequency
The frequency-dependent lossy dielectric is characterized by a dielectric constant (εr) and a loss tangent
(tanδ) at two frequencies. In terms of Maxwell's equations, the real parts of the dielectric constant and
the conductivity are given by:
ε
− ε ∞
ε = ε ∞ +
+ ωτ
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
66 of ANSYS, Inc. and its subsidiaries and affiliates.
High-Frequency Electromagnetic Material Models
where:
σ0 = DC conductivity
ε0 = free space permittivity
εr and σ are determined by the four parameters: εrs, σ0, ε ∞ , and τ. Experimentally measured values of
the dielectric constant and loss tangent are usually available at two frequencies:
The lower frequency data is considered static or DC values. Accordingly, the static permittivity and DC
conductivity are given by:
εrs = εr1
σ0 = ω1ε0εr1tanδ1
If ε∞ is known from experimental measurements, the Debye's model can then be completely defined
by calculating the relaxation time τ by:
− −
τ=
ω
where:
= ωε ε − ε∞ σ and σ = ωε ε δ − σ
Knowing τ, the material characteristics can be determined over the entire frequency range.
If ε∞ is unknown, the following equations (based on the higher frequency data) can be solved simul-
taneously for ε∞ and τ.
ε ε − ε∞
τ=
σ
− σ
+ ω τ
ε
∞ = ε
− ε
− ε
ω τ
where:
σ2 = ω2ε0εr2tanδ2
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 67
Chapter 3: Material Models
Con- Meaning
stant
C1 Lower frequency (f1) at which measured data is con-
sidered static or DC values.
C2 Higher frequency (f2) at which measured data is avail-
able.
C3 Relative permittivity at lower frequency (εr1).
C4 Relative permittivity at lower frequency (εr2).
C5
Relative permittivity at optical frequency ( εr∞ ). Input
if known. Calculated if it is not known.
C6 DC conductivity (σ0). It does not have to be defined if
the loss tangent at lower frequency is defined.
C7 Loss tangent at lower frequency (tanδ1). It does not have
to be defined if the DC conductivity is defined.
C8 Loss tangent at higher frequency (tanδ2).
The element must lie in an X-Y plane and a 3 x 3 matrix relates terms ordered by x, y, and z.
xx
= yy
zz
Constant Meaning
C1-C9 Xxx, Xyy, Xzz, Xxy, 0, 0, Xyx, 0, 0
Note
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
68 of ANSYS, Inc. and its subsidiaries and affiliates.
Piezoelectric Material Model
The full 6 x 6 elastic coefficient matrix [D] relates terms ordered x, y, z, xy, yz, xz via 21 constants as
shown below.
11
21 22
31 32 33 Symmetric
41 42 43 44
51 52 53 54 55
61 62 63 64 65 66
For 2-D problems, a 4 x 4 matrix relates terms ordered x, y, z, xy via 10 constants (D11, D21, D22, D31,
D32, D33, D41, D42, D43, D44). Note, the order of the vector is expected as {x, y, z, xy, yz, xz}, whereas for
some published materials the order is given as {x, y, z, yz, xz, xy}. This difference requires the "D" matrix
terms to be converted to the expected format. The "D" matrix can be defined in either "stiffness" form
(with units of Force/Area operating on the strain vector) or in "compliance" form (with units of the inverse
of Force/Area operating on the stress vector), whichever is more convenient. Select a form using TBOPT
on the TB command. Both forms use the same data table input as described below.
Enter the constants of the elastic coefficient matrix in the data table via the TB family of commands.
Initialize the constant table with TB,ANEL. Define the temperature with TBTEMP, followed by up to 21
constants input with TBDATA commands. The matrix may be input in either stiffness or flexibility form,
based on the TBOPT value. For the coupled-field elements, temperature- dependent matrix terms are
not allowed. You can define up to six temperature-dependent sets of constants (NTEMP = 6 max on the
TB command) in this manner. Matrix terms are linearly interpolated between temperature points. The
constants (C1-C21) entered on TBDATA (6 per command) are:
Con- Meaning
stant
C1-C6 Terms D11, D21, D31, D41, D51, D61
C7-C12 Terms D22, D32, D42, D52, D62, D33
C13-C18 Terms D43, D53, D63, D44, D54, D64
C19-C21 Terms D55, D65, D66
For a list of the elements that support this material model, see "Element Support for Material Models".
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 69
Chapter 3: Material Models
Input the dielectric constants either by specifying orthotropic dielectric permittivity (PERX, PERY, PERZ)
on the MP command or by specifying the terms of the anisotropic permittivity matrix [ε] on the TB,DPER
command. The values input on the MP command will be interpreted as permittivity at constant strain
[εS]. Using TB,DPER, you can specify either permittivity at constant strain [εS] (TBOPT = 0), or permittivity
at constant stress [εT] (TBOPT = 1).
Input the elastic coefficient matrix [c] either by specifying the stiffness constants (EX, EY, etc.) with MP
commands, or by specifying the terms of the anisotropic elasticity matrix with TB commands as described
in Anisotropic (p. 26).
You can define the piezoelectric matrix in [e] form (piezoelectric stress matrix) or in [d] form (piezoelectric
strain matrix). The [e] matrix is typically associated with the input of the anisotropic elasticity in the
form of the stiffness matrix [c], and the permittivity at constant strain [εS]. The [d] matrix is associated
with the input of compliance matrix [s] and permittivity at constant stress [εT]. Select the appropriate
matrix form for your analysis using the TB,PIEZ command.
The full 6 x 3 piezoelectric matrix relates terms x, y, z, xy, yz, xz to x, y, z via 18 constants as shown:
11 12 13
21 22 23
31 32 33
41 42 43
51 52 53
61 62 63
For 2-D problems, a 4 x 2 matrix relates terms ordered x, y, z, xy via 8 constants (e11, e12, e21, e22, e31,
e32, e41, e42). The order of the vector is expected as {x, y, z, xy, yz, xz}, whereas for some published ma-
terials the order is given as {x, y, z, yz, xz, xy}. This difference requires the piezoelectric matrix terms to
be converted to the expected format.
Use the TB commands to enter the constants of the piezoelectric matrix in the data table. Initialize the
constant table with TB,PIEZ. You can define up to 18 constants (C1-C18) with TBDATA commands (6
per command):
Con- Meaning
stant
C1-C6 Terms e11, e12, e13, e21, e22, e23
C7-C12 Terms e31, e32, e33, e41, e42, e43
C13-C18 Terms e51, e52, e53, e61, e62, e63
See Piezoelectric Analysis in the Coupled-Field Analysis Guide for more information on this material
model.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
70 of ANSYS, Inc. and its subsidiaries and affiliates.
Piezoresistive Material Model
You can define the piezoresistive matrix either in the form of piezoresistive stress matrix [π] (TBOPT =
0) or piezoresistive strain matrix [m] (TBOPT = 1).
The piezoresistive stress matrix [π] uses stress to calculate the change in electric resistivity due to
piezoresistive effect, while the piezoresistive strain matrix [m] (TBOPT = 1) uses strain to calculate the
change in electric resistivity. See Piezoresistivity in the Mechanical APDL Theory Reference for more in-
formation.
The full 6x6 piezoresistive matrix relates the x, y, z, xy, yz, xz terms of stress to the x, y, z, xy, yz, xz terms
of electric resistivity via 36 constants:
Con- Meaning
stant
C1-C6 Terms π11, π12, π13, π14, π15, π16
C7-C12 Terms π21, π22, π23, π24, π25, π26
C13-C18 Terms π31, π32, π33, π34, π35, π36
C19-C24 Terms π41, π42, π43, π44, π45, π46
C25-C30 Terms π51, π52, π53, π54, π55, π56
C31-C36 Terms π61, π62, π63, π64, π65, π66
For 2-D problems, a 4x4 matrix relates terms ordered x, y, z, xy via 16 constants.
Con- Meaning
stant
C1-C4 Terms π11, π12, π13, π14
C7-C10 Terms π21, π22, π23, π24
C13-C16 Terms π31, π32, π33, π34
C19-C22 Terms π41, π42, π43, π44
The order of the vector is expected as {x, y, z, xy, yz, xz}, whereas for some published materials the order
is given as {x, y, z, yz, xz, xy}. This difference requires the piezoresistive matrix terms to be converted
to the expected format.
See Piezoresistive Analysis in the Coupled-Field Analysis Guide for more information on this material
model.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 71
Chapter 3: Material Models
The program converts matrix [εT] to [εS] using piezoelectric strain and stress matrices.
The full 3x3 electric permittivity matrix relates x, y, z components of electric field to the x, y, z components
of electric flux density via 6 constants:
ε
11 ε12 ε13
ε22 ε23
sym
ε33
Constant Meaning
C1-C6 ε11, ε22, ε33, ε12, ε23, ε13
For 2-D problems, a 2x2 matrix relates terms ordered x, y via 3 constants (ε11 ε22 ε12):
Constant Meaning
C1, C2, C4 ε11, ε22, ε12
For further information about rate-dependent plastic (viscoplastic) material options, see Rate-Dependent
Plasticity in the Mechanical APDL Theory Reference and Viscoplasticity in the Structural Analysis Guide.
σ
/
εɺ pl = γ −
σo
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
72 of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Dependent Plastic (Viscoplastic) Material Models
1/ m
ɺε = γ σ −
pl
σo
In both cases, σ is the material yield stress, εɺ is the equivalent plastic strain rate, m is the strain rate
hardening parameter, γ is the material viscosity parameter, and σo is the static yield stress of material.
σo is a function of some hardening parameter and can be defined by isotropic plasticity (for example,
TB,BISO). As γ approaches ∞ , or m approaches zero, or εɺ approaches zero, the solution approaches
the static (rate-independent) solution.
When m is very small, the Peirce model has less difficulty converging as compared to the Perzyna
model.
For details, see Rate-Dependent Plasticity in the Mechanical APDL Theory Reference.
The two material constants for the Perzyna and Peirce models (defined by the TBDATA) are:
Con- Meaning
stant
C1 m - Material strain rate hardening para-
meter
C2 γ - Material viscosity parameter
−b ε
σ
= +
ε + ∞ −
where .
This six material constants in the EVH option are input via the data table command (TBDATA) in the
order shown:
Con- Meaning
stant
C1 Ko = Material strain hardening paramet-
er
C2 Ro = Material strain hardening parameter
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 73
Chapter 3: Material Models
Con- Meaning
stant
C3 R∞ = Material strain hardening paramet-
er
C4 b = Material strain hardening parameter
C5 m = Material strain rate hardening
parameter
C6 K = Material viscosity parameter
Details for the Anand option appear in Anand Viscoplasticity Option in the Mechanical APDL Theory
Reference.
This option requires nine material constants input via the data table command (TBDATA) in the order
shown:
1. Initialize the data table via TB,RATE and specify the model option (TBOPT) as described above.
2. Define the temperature (TBTEMP) for each set of data.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
74 of ANSYS, Inc. and its subsidiaries and affiliates.
Gasket Material Model
You can define up to nine material constants via TBDATA commands, but only six constants per
command. For the Anand model, therefore, you must issue the TBDATA command more than
once.
The data table command (TBDATA) must be defined for the same material number to specify the
static hardening behavior of the materials (rate-independent and isotropic).
You can define some general parameters including the initial gap, stable stiffness for numerical stabil-
ization, and stress cap for a gasket in tension. You can also directly input data for the experimentally
measured complex pressure closure curves for the gaskets.
Sub-options are also available to define gasket unloading behavior including linear and nonlinear un-
loading. Linear unloading simplifies the input by defining the starting closure at the compression curves
and the slope. Nonlinear unloading option allows you to directly input unloading curves to more accur-
ately model the gasket unloading behavior. When no unloading curves are defined, the material beha-
vior follows the compression curve while it is unloaded.
Enter the general parameters and the pressure closure behavior data via the TBOPT option on the
TB,GASKET command. Input the material data (TBDATA or TBPT) as shown in the following table:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 75
Chapter 3: Material Models
1. Stable stiffness is used for numerical stabilization such as the case when the gasket is opened up
and thus no stiffness is contributed to the element nodes, which in turn may cause numerical
difficulty.
2. Multiple curves may be required to define the complex nonlinear unloading behavior of a gasket
material.
When there are several nonlinear unloading curves defined, the program requires that the
starting point of each unloading curve be on the compression curve to ensure the gasket
unloading behavior is correctly simulated. Though it is not a requirement that the temperature
dependency of unloading data be the same as the compression data, when there is a missing
temperature, the program uses linear interpolation to obtain the material data of the missing
temperature. This may result in a mismatch between the compression data and the unloading
data. Therefore, it is generally recommended that the number of temperatures and temper-
ature points be the same for each unloading curve and compression curve.
When using the material GUI to enter data for the nonlinear unloading curves, an indicator
at the top of the dialog box states the number of the unloading curve whose data is currently
displayed along with the total number of unloading curves defined for the particular material
(example: Curve number 2/5). To enter data for the multiple unloading curves, type the data
for the first unloading curve, then click on the Add Curve button and type the data for the
second curve. Repeat this procedure for entering data for the remaining curves. Click the Del
Curve button if you want to remove the curve whose data is currently displayed. Click the
> button to view the data for the next curve in the sequence, or click the < button to view
the data for the previous curve in the sequence. To insert a curve at a particular location in
the sequence, click on the > or < buttons to move to the curve before the insertion location
point and click on the Add Curve button. For example, if the data for Curve number 2/5 is
currently displayed and you click on the Add Curve button, the dialog box changes to allow
you to enter data for Curve number 3/6. You can define a total of 100 nonlinear unloading
curves per material.
You can enter temperature-dependent data (TBTEMP) for any of the gasket data types. For the first
temperature curve, issue TB,GASKET,,,,TBOPT, then input the first temperature using TBTEMP, followed
by the data using either TBDATA or TBPT depending on the value of TBOPT as shown in the table.
The program automatically interpolates the temperature data to the material points using linear inter-
polation. When the temperature is out of the specified range, the closest temperature point is used.
For more information, see Gasket Material in the Mechanical APDL Theory Reference.
For a detailed description of the gasket joint simulation capability, see "Gasket Joints Simulation" in the
Structural Analysis Guide.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
76 of ANSYS, Inc. and its subsidiaries and affiliates.
Creep Equations
The creep strain rate, εcr , can be a function of stress, strain, temperature, and neutron flux level. Libraries
ɺ
of creep strain rate equations are included under the Implicit Creep Equations (p. 77) and Explicit Creep
Equations (p. 79) sections. Enter the constants shown in these equations using TB,CREEP and TBDATA
as described below. These equations (expressed in incremental form) are characteristic of materials
being used in creep design applications (see the Mechanical APDL Theory Reference for details).
For a list of the elements that support creep behavior, see "Element Support for Material Models".
• Primary creep
• Secondary creep
• Irradiation induced creep
You can define the combined effects of more than one type of creep using the implicit equations specified
by TBOPT = 11 or 12, the explicit equations, or a user-defined creep equation.
The program analyzes creep using the implicit and the explicit time-integration method. The implicit
method is robust, fast, accurate, and recommended for general use, especially with problems involving
large creep strain and large deformation. It has provisions for including temperature-dependent constants.
The program can model pure creep, creep with isotropic hardening plasticity, and creep with kinematic
hardening plasticity, using both von Mises and Hill potentials. See Material Model Combinations (p. 151)
for further information. Since the creep and plasticity are modeled simultaneously (no superposition),
the implicit method is more accurate and efficient than the explicit method. Temperature dependency
can also be incorporated by the Arrhenius function (see the Mechanical APDL Theory Reference for details).
The explicit method is useful for cases involving very small time steps, such as in transient analyses.
There are no provisions for temperature-dependent constants, nor simultaneous modeling of creep
with any other material models such as plasticity. However, there is temperature dependency using the
Arrhenius function, and you can combine explicit creep with other plasticity options using non-simul-
taneous modeling (superposition). In these cases, the program first performs the plastic analysis, then
the creep calculation.
The terms implicit and explicit, as applied to creep, have no relationship to “explicit dynamics,” or any
elements referred to as “explicit elements.”
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 77
Chapter 3: Material Models
Creep Mod-
el
Name Equation Type
(TB-
OPT)
Time Harden-
2 εɺ cr = 1σC2 C3 −C4/ T C1>0 Primary
ing
Generalized Ex- C1>0,
3 εɺ = σ −t , = 5 σ − Primary
ponential C5>0
Generalized
4 εɺ
= σ + + 6 7 −8
C1>0 Primary
Graham
∂εg
εɺ gh = j εg = + εɺ m
∂ + Primary
Rational poly-
12 inσ io C2>0 + Sec-
nomial εɺ m = k σ
ondary
i i
= l εɺ mp σiq = js εɺ muuσiuv
x −y { z
εwx =
Generalized
13 Time Harden- = |σ + } σ} + ~ σ~ Primary
ing
= + σ
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
78 of ANSYS, Inc. and its subsidiaries and affiliates.
Creep Equations
T = temperature (absolute). The offset temperature (from TOFFST), is internally added to all
temperatures for convenience.
C1 through C12 = constants defined by the TBDATA command
t = time at end of substep
e = natural logarithm base
You can define the user creep option by setting TBOPT = 100, and using TB,STATE to specify the
number of state variables for the user creep subroutine. See the Guide to ANSYS User Programmable
Features for more information. The RATE command is necessary to activate implicit creep for specific
elements (see the RATE command description for details). The RATE command has no effect for explicit
creep.
For temperature-dependent constants, define the temperature using TBTEMP for each set of data. Then,
define constants C1 through Cm using TBDATA (where m is the number of constants, and depends on
the creep model you choose).
The following example shows how you would define the implicit creep model represented by TBOPT
= 1 at two temperature points.
TB,CREEP,1,,,1 !Activate creep data table, specify creep model 1
TBTEMP,100 !Define first temperature
TBDATA,1,c11,c12,c13,c14 !Creep constants c11, c12, c13, c14 at first temp.
TBTEMP,200 !Define second temperature
TBDATA,1,c21,c22,c23,c24 !Creep constants c21, c22, c23, c24 at second temp.
Coefficients are linearly interpolated for temperatures that fall between user defined TBTEMP values.
For some creep models, where the change in coefficients spans several orders of magnitude, this linear
interpolation might introduce inaccuracies in solution results. Use enough curves to accurately capture
the temperature dependency. Also, consider using the curve fitting subroutine to calculate a temperature
dependent coefficient that includes the Arrhenius term.
When a temperature is outside the range of defined temperature values, the program uses the coefficients
defined for the constant temperature.
For a list of elements that can be used with this material option, see "Element Support for Material
Models".
See Creep in the Structural Analysis Guide for more information on this material option.
Specify primary creep with constant C6. Primary Explicit Creep Equation for C6 = 0 (p. 80), through Primary
Explicit Creep Equation for C6 = 100 (p. 86), show the available equations. You select an equation with
the appropriate value of C6 (0 to 15). If C1 ≤ 0, or if T + Toffset ≤ 0, no primary creep is computed.
Specify secondary creep with constant C12. Secondary Explicit Creep Equation for C12 = 0 (p. 86) and
Secondary Explicit Creep Equation for C12 = 1 (p. 87) show the available equations. You select an equation
with the appropriate value of C12 (0 or 1). If C7 ≤ 0, or if T + Toffset ≤ 0, no secondary creep is computed.
Also, primary creep equations C6 = 9, 10, 11, 13, 14, and 15 bypass any secondary creep equations since
secondary effects are included in the primary part.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 79
Chapter 3: Material Models
Specify irradiation induced creep with constant C66. Irradiation Induced Explicit Creep Equation for C66 =
5 (p. 87) shows the single equation currently available; select it with C66 = 5. This equation can be used
in conjunction with equations C6 = 0 to 11. The constants should be entered into the data table as in-
dicated by their subscripts. If C55 ≤ 0 and C61 ≤ 0, or if T + Toffset ≤ 0, no irradiation induced creep is
computed.
A linear stepping function is used to calculate the change in the creep strain within a time step (∆ εcr
= ( εcr )(∆t)). The creep strain rate is evaluated at the condition corresponding to the beginning of the
ɺ
time interval and is assumed to remain constant over the time interval. If the time step is less than 1.0e-
6, then no creep strain increment is computed. Primary equivalent stresses and strains are used to
evaluate the creep strain rate. For highly nonlinear creep strain vs. time curves, use a small time step
if you are using the explicit creep algorithm. A creep time step optimization procedure is available for
automatically increasing the time step whenever possible. A nonlinear stepping function (based on an
exponential decay) is also available (C11 = 1) but should be used with caution since it can underestimate
the total creep strain where primary stresses dominate. This function is available only for creep equations
C6 = 0, 1 and 2. Temperatures used in the creep equations should be based on an absolute scale
(TOFFST).
Use the BF or BFE commands to enter temperature and fluence values. The input fluence (Φt) includes
the integrated effect of time and time explicitly input is not used in the fluence calculation. Also, for
the usual case of a constant flux (Φ), the fluence should be linearly ramp changed.
Temperature dependent creep constants are not permitted for explicit creep. You can incorporate
other creep options by setting C6 = 100. See the Guide to ANSYS User Programmable Features for more
information.
The following example shows how you would use the explicit creep equation defined by C6 = 1.
TB,CREEP,1 !Activate creep data table
TBDATA,1,c1,c2,c3,c4,,1 !Creep constants c1, c2, c3, c4 for equation C6=1
The explicit creep constants that you enter with the TBDATA are:
Constant Meaning
C1-CN Constants C1, C2, C3, etc. (as defined in Primary Explicit Creep
Equation for C6 = 0 to Irradiation Induced Explicit Creep Equation
for C66 = 5) These are obtained by curve fitting test results for
your material to the equation you choose. Exceptions are
defined below.
εɺ = 1σC2ε C3 −C4/ T
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
80 of ANSYS, Inc. and its subsidiaries and affiliates.
Creep Equations
εɺ cr = 1σC2 C3 −C4/ T
εɺ = σ −t
where:
= 5 σ −
∂ε
εɺ
=
∂
ε = ε x − − s + ε − − + εɺ m
where:
εx = 0 for σ ≤ C2
εx = G + H σ for C2 < σ ≤ C3
C2 = 6000 psi (default), C3 = 25000 psi (default)
This double exponential equation is valid for Annealed 304 Stainless Steel over a temperature range
from 800 to 1100°F. The equation, known as the Blackburn creep equation when C1 = 1, is described
completely in the High Alloy Steels. The first two terms describe the primary creep strain and the last
term describes the secondary creep strain.
To use this equation, input a nonzero value for C1, C6 = 9.0, and C7 = 0.0. Temperatures should be in
°R (or °F with Toffset = 460.0). Conversion to °K for the built-in property tables is done internally. If the
temperature is below the valid range, no creep is computed. Time should be in hours and stress in psi.
The valid stress range is 6,000 - 25,000 psi.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 81
Chapter 3: Material Models
εc + + εɺ m
+
where:
This standard rational polynomial creep equation is valid for Annealed 304 SS over a temperature range
from 427°C to 704°C. The equation is described completely in the High Alloy Steels. The first term de-
scribes the primary creep strain. The last term describes the secondary creep strain. The average "lot
constant" is used to calculate ε .
ɺ
To use this equation, input C1 = 1.0, C4 = 1.0, C6 = 9.0, and C7 = 0.0. Temperature must be in °C and
Toffset must be 273 (because of the built-in property tables). If the temperature is below the valid range,
no creep is computed. Also, time must be in hours and stress in Megapascals (MPa).
Various hardening rules governing the rate of change of creep strain during load reversal may be selected
with the C5 value: 0.0 - time hardening, 1.0 - total creep strain hardening, 2.0 - primary creep strain
hardening. These options are available only with the standard rational polynomial creep equation.
This standard rational polynomial equation is the same as described above except that temperature
must be in °F, Toffset must be 460, and stress must be in psi. The equivalent valid temperature range is
800 - 1300°F.
1 ∂
∂ε
εɺ r =
This equation, also described in High Alloy Steels, differs from the Annealed 304 SS equation in that
the built-in property tables are for Annealed 316 SS, the valid stress range is 4000 - 30,000 psi, C2 defaults
to 4000 psi, C3 defaults to 30,000 psi, and the equation is called with C6 = 10.0 instead of C6 = 9.0.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
82 of ANSYS, Inc. and its subsidiaries and affiliates.
Creep Equations
This standard rational polynomial equation, also described in High Alloy Steels, differs from the Annealed
304 SS equation in that the built-in property tables are for Annealed 316 SS, the valid temperature
range is 482 - 704°C, and the equation is called with C6 = 10.0 instead of C6 = 9.0. The hardening rules
for load reversal described for the C6 = 9.0 standard Rational Polynomial creep equation are also available.
The average "lot constant" from High Alloy Steels is used in the calculation of εm .
ɺ
This standard rational polynomial equation is the same as described above except that the temperatures
must be in °F, Toffset must be 460, and the stress must be in psi (with a valid range from 0.0 to 24220
psi). The equivalent valid temperature range is 900 - 1300°F.
∂εc
εɺ cr = 1
∂
ε = + εɺ
+
This modified rational polynomial equation is valid for Annealed 2 1/4 Cr -1 Mo Low Alloy steel over a
temperature range of 700 - 1100°F. The equation is described completely in the Low Alloy Steels. The
first term describes the primary creep strain and the last term describes the secondary creep strain. No
modification is made for plastic strains.
To use this equation, input C1 = 1.0, C6 = 11.0, and C7 = 0.0. Temperatures must be in °R (or °F with
Toffset = 460.0). Conversion to °K for the built-in property tables is done internally. If the temperature is
below the valid range, no creep is computed. Time should be in hours and stress in psi. Valid stress
range is 1000 - 65,000 psi.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 83
Chapter 3: Material Models
εc + + εɺ m
+
where:
This standard rational polynomial creep equation is valid for Annealed 2 1/4 Cr - 1 Mo Low Alloy Steel
over a temperature range from 371°C to 593°C. The equation is described completely in the Low Alloy
Steels. The first term describes the primary creep strain and the last term describes the secondary creep
strain. No tertiary creep strain is calculated. Only Type I (and not Type II) creep is supported. No modi-
fication is made for plastic strains.
To use this equation, input C1 = 1.0, C4 = 1.0, C6 = 11.0, and C7 = 0.0. Temperatures must be in °C and
Toffset must be 273 (because of the built-in property tables). If the temperature is below the valid range,
no creep is computed. Also, time must be in hours and stress in Megapascals (MPa). The hardening
rules for load reversal described for the C6 = 9.0 standard Rational Polynomial creep equation are also
available.
This standard rational polynomial equation is the same as described above except that temperatures
must be in °F, Toffset must be 460, and stress must be in psi. The equivalent valid temperature range is
700 - 1100°F.
εɺ r = N (M−1)
1σ
where:
C1 = Scaling constant
M, N, K = Function of temperature (determined by linear interpolation within table) as follows:
Constant Meaning
C5 Number of temperature values to describe
M, N, or K function (2 minimum, 6 maximum)
C49 First absolute temperature value
C50 Second absolute temperature value
...
C48 + C5 C5th absolute temperature value
C48 + C5 + 1 First M value
...
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
84 of ANSYS, Inc. and its subsidiaries and affiliates.
Creep Equations
Constant Meaning
C48 + 2C5 C5th M value
C48 + 2C5 C5th M value
...
C48 + 2C5 C5th M value
C48 + 2C5 + 1 First N value
...
C48 + 3C5 C5th N value
C48 + 3C5 + 1 First K value
...
This power function creep law having temperature dependent coefficients is similar to Equation C6 =
1.0 except with C1 = f1(T), C2 = f2(T), C3 = f3(T), and C4 = 0. Temperatures must not be input in decreasing
order.
εacc
εɺ cr =
εB
acc σ
A ( 3 A + 2B + C)
where:
εacc = creep strain accumulated to this time (calculated by the program). Internally set to 1 x
10-5 at the first substep with nonzero time to prevent division by zero.
A = C1/T
B = C2/T + C3
C = C4/T + C5
This equation is often referred to as the Sterling Power Function creep equation. Constant C7 should
be 0.0. Constant C1 should not be 0.0, unless no creep is to be calculated.
∂ε
εɺ = 1
∂
where:
εc = cpt/(1+pt) + εm
ɺ
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 85
Chapter 3: Material Models
This creep law is valid for Annealed 316 SS over a temperature range from 800°F to 1300°F. The equation
is similar to that given for C6 = 10.0 and is also described in High Alloy Steels.
To use equation, input C1 = 1.0 and C6 = 14.0. Temperatures should be in °R (or °F with Toffset = 460).
Time should be in hours. Constants are only valid for English units (pounds and inches). Valid temper-
ature range: 800° - 1300°F. Maximum stress allowed for ec calculation: 45,000 psi; minimum stress: 0.0
psi. If T + Toffset < 1160, no creep is computed.
∂εc
εɺ cr = 1
∂
where:
ε + + εɺ m
+
C3σ C4
εɺ = 2 σ ( ust not be negative)
= 7 εɺ 8 σ9
= 0 εɺ σ
This rational polynomial creep equation is a generalized form of the standard rational polynomial
equations given as C6 = 9.0, 10.0, and 11.0 (C4 = 1.0 and 2.0). This equation reduces to the standard
equations for isothermal cases. The hardening rules for load reversal described for the C6 = 9.0 standard
Rational Polynomial creep equation are also available.
εɺ = σ / − / T
where:
σ = equivalent stress
T = temperature (absolute). The offset temperature (from TOFFST), is internally added to all
temperatures for convenience.
t = time
e = natural logarithm base
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
86 of ANSYS, Inc. and its subsidiaries and affiliates.
Swelling Model
εɺ cr = 7 σC8 −C10 / T
where:
B = FG + C63
−
=
9 + −
= − −φ 2
σ = equivalent stress
T = temperature (absolute). The offset temperature (from TOFFST) is internally added to all
temperatures for convenience.
Φt0.5 = neutron fluence (input on BF or BFE command)
e = natural logarithm base
t = time
This irradiation induced creep equation is valid for 20% Cold Worked 316 SS over a temperature range
from 700° to 1300°F. Constants 56, 57, 58 and 62 must be positive if the B term is included.
See Creep in the Structural Analysis Guide for more information on this material option.
Irradiation-induced swelling and creep apply to metal alloys that are exposed to nuclear radiation.
However, the swelling equations and the fluence input may be completely unrelated to nuclear swelling.
You can also model other types of swelling behavior, such as moisture-induced volume expansion.
where ε is the total mechanical strain, εel is the elastic strain, εpl is the plastic strain, and εsw is the
swelling strain.
You can combine swelling strain with other material models such as plasticity and creep; however, you
can use swelling with any hyperelasticity or anisotropic hyperelasticity material model.
Irradiation-induced swelling is generally accompanied by irradiation creep for metals and composites,
such as silicon carbide (SiC). The irradiation-induced swelling strain rate may depend on temperature,
time, fluence (the flux x time), and stress, such as:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 87
Chapter 3: Material Models
εɺ sw εsw ( Φt σ)
where t is time, T is the temperature, Φt is the fluence, and σ is the stress. Temperatures used in the
swelling equations should be based on an absolute scale (TOFFST). Specify temperature and fluence
values via the BF or BFE command.
• Linear swelling defines swelling strain rate as a function of fluence rate, expressed as:
εɺ Φ
ɺ
ε 1+ 2Φ + 3 − 4 Φ −
• A user-defined swelling option is available if you wish to create your own swelling function. For more
information, see userswstrain in the Guide to ANSYS User Programmable Features.
For highly nonlinear swelling strain vs. fluence curves, it is good practice to use a small fluence step for
better accuracy and solution stability. If time is changing, a constant flux requires a linearly changing
fluence (because the swelling model uses fluence [Φt] rather than flux [Φ]).
Initialize the swelling table (TB,SWELL) with the desired data table option (TBOPT), as follows:
Issue the TBDATA command to enter the swelling table constants (up to six per command), as shown
in the table.
For a list of the elements that you can use with the swelling model, see "Element Support for Material
Models"
For more information about this material model, see Swelling in the Structural Analysis Guide.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
88 of ANSYS, Inc. and its subsidiaries and affiliates.
Shape Memory Alloy (SMA) Material Model
Such distinct material behavior is due to the material microstructure in which there exists two different
crystallographic structures, one characterized by austenite (A), and another one by martensite (M).
Austenite is the crystallographically more-ordered phase, and martensite is the crystallographically less-
ordered phase. The key characteristic of an SMA is the occurrence of a martensitic phase transformation.
Typically, the austenite is stable at high temperatures and low stress, while the martensite is stable at
low temperatures and high stress. The reversible martensitic phase transformation results in unique
effects: the pseudoelasticity (PE) and the shape memory effect (SME).
As shown by (a) in the following figure, whenever σL is positive, the specimen recovers its original shape
completely and returns to a stress-free configuration (PE).
As shown by (b) in the figure, when σL is negative, residual strains (E and E') can be observed after un-
loading into a stress-free configuration. If the material is heated, then eventually σL becomes positive;
however, the admissible configuration under a stress-free state points to A. The material therefore un-
dergoes an inverse transformation process (SME).
Nitinol
A typical shape memory alloy is Nitinol, a nickel titanium (Ni-Ti) alloy discovered in the 1960s at
the U.S. Naval Ordnance Laboratory (NOL). The acronym NiTi-NOL (or Nitinol) has since been
commonly used when referring to Ni-Ti-based shape memory alloys.
Two SMA material model options (accessed via TB,SMA) are available, one for simulating superelastic
behavior and the other for simulating the shape memory effect behavior of shape memory alloys.
The material option for superelasticity is based on Auricchio et al. [1 (p. 97)] in which the material un-
dergoes large-deformation without showing permanent deformation under isothermal conditions, as
shown by (a) in Figure 3.5 (p. 89). The material option for the shape memory effect is based on the 3-
D thermomechanical model for stress-induced solid phase transformations [2 (p. 97)] [3 (p. 97)] [4 (p. 97)].
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 89
Chapter 3: Material Models
Two of the phase transformations are considered here: A->S and S->A. The material is composed of
two phases, the austenite (A) and the martensite (S). Two internal variables, the martensite fraction (ξS)
and the austenite fraction (ξA), are introduced. One of them is a dependent variable, and they are assumed
to satisfy the relation expressed as:
ξS + ξ A =
The material behavior is assumed to be isotropic. The pressure dependency of the phase transformation
is modeled by introducing the Drucker-Prager loading function, as follows:
= + α
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
90 of ANSYS, Inc. and its subsidiaries and affiliates.
Shape Memory Alloy (SMA) Material Model
=σ
= σ
where α is the material parameter, σ is the stress, and 1 is the identity tensor.
AS ɺ
− −ξ → transormation
− AS
ξɺ S = f
SA ξ ɺ
→ trransormation
S SA
− f
where:
+ α
= σ
= σ + α
σ
d σ
σ
∫
σ
σ
σ
∫
ε
εL
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 91
Chapter 3: Material Models
AS < < AS
AS i ɺ f
= >
otherwise
< <
ɺ
= <
= σ +α
The material parameter α characterizes the material response in tension and compression. If tensile and
compressive behaviors are the same, then α = 0. For a uniaxial tension-compression test, α can be related
to the initial value of austenite to martensite phase transformation in tension and compression
σc σ
( , respectively) as:
σ!"# − σ $"#
α=
σ!"# + σ $"#
σ= ε − ε %&
ɺ '( = ξɺ ∂
L
∂σ
)* +
where D is the elastic stiffness tensor, is the transformation strain tensor, and is the material
parameter shown in Figure 3.7 (p. 91).
The superelastic SMA option is described by six constants that define the stress-strain behavior in
loading and unloading for the uniaxial stress-state.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
92 of ANSYS, Inc. and its subsidiaries and affiliates.
Shape Memory Alloy (SMA) Material Model
For each data set, define the temperature (TBTEMP), then define constants C1 through C6 (TBDATA).
You can define up to 99 sets of temperature-dependent constants in this manner.
σsAS
C1 Starting stress value for the forward
phase transformation
σf
C2 Final stress value for the forward phase
transformation
σ
C3 Starting stress value for the reverse
phase transformation
σ
C4 Final stress value for the reverse phase
transformation
C5 Maximum residual strain
L
C6 α Parameter measuring the difference
between material responses in tension
and compression
2
Ψε ε tr = ( ε − ε tr ) ( ε − ε tr ) + τM ( ) ε'tr + ε'tr + Ι ε
( ε'tr )
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 93
Chapter 3: Material Models
τM(T) = a positive and monotonically increasing function of the temperature as 〈β(T - T0)〉+ in which
〈·〉+ is the positive part of the argument (also known as Maxwell stress).
β = material parameter
T = temperature
T0 = temperature below which no twinned martensite is observed
h = material parameter related to the hardening of the material during the phase transformation
ε'tr ( ) = indicator function introduced to satisfy the constraint on the transformation norm
≤ ε ≤ εL
Ι ε ( ε ) =
+∞
∂ψ
σ=
∂ε
∂ψ
∈−
∂ε
Stresses, strains, and the transformation strains are then related as follows:
σ= (ε − ε )
=σ−
= σ
where S is the deviatoric stress and p is the volumetric stress (also called hydrostatic pressure)
ε
= − τM ( )+ ε + γ
ε
where γ is defined by
≤ ε < ε
³ =
≥ ε = ε
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
94 of ANSYS, Inc. and its subsidiaries and affiliates.
Shape Memory Alloy (SMA) Material Model
where L
is a maximum transformation strain.
Numerous experimental tests show an asymmetric behavior of SMA in tension and compression, and
suggest describing SMA as an isotropic material with a Prager-Lode-type limit surface. Accordingly, the
following yield function is assumed:
tr
= 2
+ 3
−
2
where Xtr is the transformation stress, J2 and J3 are the second and third invariants of transformation
stress, m is a material parameter related to Lode dependency, and R is the elastic domain radius.
= (
)
= (
)
The evolution of transformation strain is defined as:
∂
εɺ = ξε
ɺ
∂σ
where ξ is an internal variable and is called as transformation strain multiplier. ξ and F(Xtr) must satisfy
the classical Kuhn-Tucker conditions, as follows:
ξɺ ≥
ξɺ (
)=
The elastic properties of austenite and martensite phase differ. During the transformation phase, the
elastic stiffness tensor of material varies with the deformation. The elastic stiffness tensor L is therefore
ε'
= ( S
− A )+ A
ε
where DA is the elastic stiffness tensor of austenite phase, and DS is the elastic stiffness tensor of
martensite phase. The Poisson’s ratio of the austenite phase is assumed to be the same as the
martensite phase. When the material is in its austenite phase, D = DA, and when the material undergoes
full transformation (martensite phase), D = DS.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 95
Chapter 3: Material Models
Figure 3.8 Admissible Paths for Elastic Behavior and Phase Transformations
The austenite phase is associated with the horizontal region abcd. Mixtures of phases are related to the
surface cdef. The martensite phase is represented by the horizontal region efgh. Point c corresponds
to the nucleation of the martensite phase. Phase transformations take place only along line cf, where
3 =
tr + 2
tr . Saturated phase transformations are represented by paths on line fg. The horizontal
region efgh contains elastic processes except, of course, those on line fg.
A backward Euler integration scheme is used to solve the stress update and the consistent tangent
stiffness matrix required by the finite element solution for obtaining a robust nonlinear solution.
The shape memory effect option is described by seven constants that define the stress-strain behavior
of material in loading and unloading cycles for the uniaxial stress-state and thermal loading.
For each data set, define the temperature (TBTEMP), then define constants C1 through C7 (TBDATA).
You can define up to 99 sets of temperature-dependent constants in this manner.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
96 of ANSYS, Inc. and its subsidiaries and affiliates.
MPC184 Joint Material Models
C6 Em Martensite modulus
C7 m Lode dependency parameter
Example 3.3 Defining Shape Memory Effect Properties of the Austenite Phase
MP,EX,1,60000.0
MP,NUXY,1,0.36
Define SMA material properties
TB,SMA,1,,,MEFF
TBDATA,1,1000, 223, 50, 2.1, 0.04, 45000
TBDATA,7,0.05
The TB command may be repeated with the same material ID number to specify both the stiffness and
damping behavior.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 97
Chapter 3: Material Models
11
21 22
31 32 33
41 42 43 44
51 52 53 54 55
61 62 63 64 65 66
Enter the stiffness or damping coefficient of the matrix in the data table with TB set of commands.
Initialize the constant table with TB,JOIN,,,STIF (for stiffness behavior) or TB,JOIN,,,DAMP (for damping
behavior). Define the temperature with TBTEMP, followed by the relevant constants input with TBDATA
commands. Matrix terms are linearly interpolated between temperature points. Based on the joint type,
the relevant constant specification is as follows:
The following example shows how you would define the uncoupled linear elastic stiffness behavior for
a universal joint at the two available components of relative motion, with two temperature points:
TB,JOIN,1,2,,STIF ! Activate JOIN material model with linear elastic stiffness
TBTEMP,100.0 ! Define first temperature
TBDATA,16,D44 ! Define constant D44 in the local ROTX direction
TBDATA,21,D66 ! Define constant D66 in the local ROTZ direction
TBTEMP,200.0 ! Define second temperature
TBDATA,16,D44 ! Define constant D44 in the local ROTX direction.
TBDATA,21,D66 ! Define constant D66 in the local ROTZ direction.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
98 of ANSYS, Inc. and its subsidiaries and affiliates.
MPC184 Joint Material Models
Use the TBPT command to specify the data points or specify the name of a function that defines the
curve on the TB command. (Use the Function Tool to generate the specified function.) The values may
be temperature-dependent.
You can specify nonlinear damping behavior in a similar manner by supplying velocity versus damping
force (or moment).
The appropriate TBOPT labels for each joint element type are shown in the following tables. For a de-
scription of each TBOPT label, see JOIN Specifications in the TB command documentation.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 99
Chapter 3: Material Models
The following example illustrates the specification of nonlinear stiffness behavior for a revolute joint
that has only one available component of relative motion (the rotation around the axis of revolution).
Two temperature points are specified.
TB,JOIN,1,2,2,JNS4
TBTEMP,100.
TBPT,,rotation_value_1,moment_value_1
TBPT,,rotation_value_2,moment_value_2
TBTEMP,200.0
TBPT,,rotation_value_1,moment_value_1
TBPT,,rotation_value_2,moment_value_2
Example Consider a function where the damping force varies with temperature and relative velocity:
Define the function using the Function Editor, then retrieve and load it using the Function Loader. (The
editor and the loader are both components of the Function Tool.)
Assuming a function name of dampfunc, you can then use the TB command to define the joint mater-
ial:
TB, JOIN, 1, , , JND4, , %dampfunc%
For more information about the Function Tool utility, see "Using the Function Tool" in the Basic Analysis
Guide.
There are three options for defining the Coulomb friction coefficient.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
100 of ANSYS, Inc. and its subsidiaries and affiliates.
MPC184 Joint Material Models
• Define a single value of the Coulomb friction coefficient by specifying TBOPT = MUSx, where the value
of x depends on the joint under consideration. Use the TBDATA command to specify the value of the
friction coefficient.
• Define the Coulomb friction coefficient as a function of the sliding velocity. Use TBOPT = MUSx (as
stated above) and use the TBPT command to specify the data values.
• Use the exponential law for friction behavior. Specify TBOPT = EXPx, where the value of x depends on
the joint under consideration, and use the TBDATA command to specify the values required for the
exponential law. In this case, the TBDATA command format is:
where µs is the coefficient of friction in the static regime, µd is the coefficient of friction in the dy-
namic regime, and c is the decay coefficient.
• The maximum allowable value of critical force/moment can be specified using TBOPT = TMXx, where
x depends on the joint under consideration.
Elastic Slip
• The elastic slip can be specified by setting TBOPT = SLx, where x depends on the joint under consid-
eration.
• If the stick-stiffness value is not specified, then this value along with the critical force/moment is used
to determine the stick-stiffness.
• If the elastic slip is not specified, then a default value is computed for stick-stiffness calculations if ne-
cessary.
Stick-Stiffness
• A stick-stiffness value can be specified for controlling the behavior in the stick regime when friction
behavior is specified. Use TBOPT = SKx, where x depends on the joint under consideration.
• If the stick-stiffness value is not specified, then the following procedure is adopted:
– If both maximum force/moment and elastic slip are specified, then the stick-stiffness is calculated
from these values.
– If only maximum force/moment is specified, then a default elastic slip is computed and then the
stick-stiffness is calculated.
– If only the elastic slip is specified, then the stick-stiffness value is computed based on the current
normal force/moment (Friction Coefficient * Normal Force or Moment/elastic-slip).
• If the forces that are generated during a joint assembly have to be modeled, the interference fit
force/moment can be specified using TBOPT = FIx, where x depends on the joint under consideration.
This force/moment will contribute to the normal force/moment in friction calculations.
The appropriate TBOPT labels (TB command) for each joint element type are shown in the table below:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 101
Chapter 3: Material Models
The following examples illustrate how to specify Coulomb friction parameters for various scenarios.
Example 1 Specifying a single value of coefficient of friction and other friction parameters for an x-
axis revolute joint.
TB, JOIN, 1, , , MUS4 ! Label for friction coefficient
TBDATA, 1, 0.1 ! Value of coefficient of friction
TB, JOIN, 1, , , SK4 ! Label for stick-stiffness
TBDATA, 1, 3.0E4 ! Value for stick-stiffness
TB, JOIN, 1, , , FI4 ! Label for interference fit force
TBDATA, 1, 10000.00 ! Value for interference fit force
Example 2 Specifying temperature dependent friction coefficient and other friction parameters for
a z-axis revolution joint.
TB, JOIN, 1,2 , 1, MUS6 ! 2 temp points, 2 data points and label for friction coefficient
TBTEMP, 10 ! 1st temperature
TBDATA, 1, 0.15 ! Value of coefficient of friction
TBTEMP, 20 ! 2nd temperature
TBDATA, 1, 0.1 ! Value of coefficient of friction
!
TB, JOIN, 1, , , SK4 ! Label for stick-stiffness
TBDATA, 1, 3.0E4 ! Value for stick-stiffness
TB, JOIN, 1, , , FI4 ! Label for interference fit force
TBDATA, 1, 10000.00 ! Value for interference fit force
Example 3 Specifying the exponential law for friction and other friction parameters for a z-axis revolute
joint.
TB, JOIN, 1, , , EXP6 ! Label for friction coefficient
TBDATA, 1, 0.4, 0.2, 0.5 ! Static friction coeff, dynamic friction coeff, decay constant
!
TB, JOIN, 1, , , SK6 ! Label for stick-stiffness
TBDATA, 1, 3.0E4 ! Value for stick-stiffness
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
102 of ANSYS, Inc. and its subsidiaries and affiliates.
Contact Friction
To define a coefficient of friction that is dependent on temperature, time, normal pressure, sliding dis-
tance, or sliding relative velocity, use the TBFIELD command. Suitable combinations of up to two fields
can be used to define dependency, for example, temperature and sliding distance as shown below:
TB,FRIC,1,,,ISO ! Activate isotropic friction model
TBFIELD,TEMP,100.0 ! Define first value of temperature
TBFIELD,SLDI,0.1 ! Define first value of sliding distance
TBDATA,1,MU ! Define coefficient of friction
TBFIELD,SLDI,0.5 ! Define second value of sliding distance
TBDATA,1,MU ! Define coefficient of friction
TBFIELD,TEMP,200.0 ! Define second value of temperature
TBFIELD,SLDI,0.2 ! Define first value of sliding distance
TBDATA,1,MU ! Define coefficient of friction
TBFIELD,SLDI,0.7 ! Define second value of sliding distance
TBDATA,1,MU ! Define coefficient of friction
See Understanding Field Variable Interpolation (p. 155) for more information on the interpolation scheme
used for field-dependent material properties defined using TBFIELD.
To define a coefficient of friction that is dependent on temperature only, use the TBTEMP command
as shown below:
TB,FRIC,1,2,,ISO ! Activate isotropic friction model
TBTEMP,100.0 ! Define first temperature
TBDATA,1,MU ! Define coefficient of friction at temp 100.0
TBTEMP,200.0 ! Define second temperature
TBDATA,1,MU ! Define coefficient of friction at temp 200.0
Alternatively, you can use MU on the MP command to specify the isotropic friction. Use the MPTEMP
command to define MU as a function of temperature. See Linear Material Properties (p. 14) for details.
Issue the TB,FRIC command with TBOPT = ORTHO to define orthotropic friction, and specify the coeffi-
cients of friction, MU1 and MU2, on the TBDATA command.
To define a coefficient of friction that is dependent on temperature, time, normal pressure, sliding dis-
tance, or sliding relative velocity, use the TBFIELD command. Suitable combinations of up to two fields
can be used to define dependency, for example, sliding relative velocity and normal pressure as shown
below:
TB,FRIC,1,,,ORTHO ! Activate orthotropic friction model
TBFIELD,SLRV,10.0 ! Define first value of sliding relative velocity
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 103
Chapter 3: Material Models
See Understanding Field Variable Interpolation (p. 155) for more information on the interpolation scheme
used for field-dependent material properties defined using TBFIELD.
To define a coefficient of friction that is dependent on temperature only, use the TBTEMP command
as shown below:
TB,FRIC,1,2,,ORTHO ! Activate orthotropic friction model
TBTEMP,100.0 ! Define first temperature
TBDATA,1,MU1,MU2 ! Define coefficients of friction at temp 100.0
TBTEMP,200.0 ! Define second temperature
TBDATA,1,MU1,MU2 ! Define coefficients of friction at temp 200.0
/SOLUTION
!* LOAD STEP 1
...
TIME,1
SOLVE
!* LOAD STEP 2
...
TIME,2
SOLVE
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
104 of ANSYS, Inc. and its subsidiaries and affiliates.
Cohesive Zone Material Model
The USERFRIC subroutine can only be used with penalty-based tangential contact (KEYOPT(2) = 0, 1,
or 3 on the contact element).
To specify user-defined friction, use the TB,FRIC command with TBOPT = USER and specify the friction
properties on the TBDATA command, as shown below. Also, use the USERFRIC subroutine to program
the friction model.
TB,FRIC,1,,2,USER ! Activate user defined friction model
TBDATA,1,PROP1,PROP2 ! Define friction properties
Field variables specified with the TBFIELD command are not available for TB,FRIC,,,,USER.
For more detailed information about cohesive zone materials, see Cohesive Zone Material (CZM) Model
in the Mechanical APDL Theory Reference.
To define a temperature dependent material, use the TBTEMP command as shown below:
TB,CZM,1,2,,EXPO ! Activate exponential material model
TBTEMP,100.0 ! Define first temperature
TBDATA,1, max, n, t ! Define material constants at temp 100.0
TBTEMP,200.0 ! Define second temperature
TBDATA,1, max, n, t ! Define material constants at temp 200.0
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 105
Chapter 3: Material Models
To define a temperature-dependent material, issue the TBTEMP command as shown in the following
example input fragment:
TB,CZM,1,2,,BILI ! Activate bilinear CZM material model
!
! Define first temperature
!
TBTEMP,100.0
!
! Define Mode I dominated material constants at temp 100.0:
!
δ δ
!TBDATA,1, max, ,- max,
,
!
! Define second temperature
!
TBTEMP,200.0
TBTEMP,200.0
!
! Define Mode I dominated material constants at temp 200.0:
δ
δ
TBDATA,1, max, ,- max, ,
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
106 of ANSYS, Inc. and its subsidiaries and affiliates.
Cohesive Zone Material Model
To define bilinear material behavior with tractions and separation distances, issue the TB,CZM,,,,CBDD
comman, then specify the following material constants via the TBDATA command:
To define a temperature dependent material, use the TBTEMP command as shown below:
TB,CZM,1,2,,CBDD ! Activate bilinear material model with tractions
! and separation distances
TBTEMP,100.0 ! Define first temperature
δ δ
TBDATA,1, max, , max, , , ! Define material constants at temp 100.0
TBTEMP,200.0 ! Define second temperature
δ δ
TBDATA,1, max, , max, , , ! Define material constants at temp 200.0
Use the TB,CZM command with TBOPT = CBDE to define bilinear material behavior with tractions and
critical fracture energies, and specify the following material constants using the TBDATA command.
To define a temperature dependent material, use the TBTEMP command as shown below:
TB,CZM,1,2,,CBDE ! Activate bilinear material model with
! tractions and facture energies
TBTEMP,100.0 ! Define first temperature
TBDATA,1, max,Gcn, max,Gct, , ! Define material constants at temp 100.0
TBTEMP,200.0 ! Define second temperature
TBDATA,1, max,Gcn, max,Gct, , ! Define material constants at temp 200.0
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 107
Chapter 3: Material Models
There are three ways to define material data for compressible fluids: liquid, gas, or pressure-volume
data.
Liquid
Use the TB,FLUID command with TBOPT = LIQUID to define material behavior for a liquid, and specify
the following material constants using the TBDATA command:
You can define a temperature dependent liquid material with up to 20 temperatures (NTEMP = 20 max
on the TB command) by using the TBTEMP command, as shown in the example below:
TB,FLUID,1,2,,LIQUID ! Activate liquid material model
TBTEMP,100.0 ! Define first temperature
TBDATA,1,K, , 0f ! Define material constants at temp 100.0
TBTEMP,200.0 ! Define second temperature
TBDATA,1,K, , 0f ! Define material constants at temp 200.0
When specifying temperature dependent density values for a liquid, keep in mind that the current
density (ρf) for hydrostatic fluid elements is computed at each iteration as a function of pressure change
(∆P), bulk modulus (K), coefficient of thermal expansion (α), and temperature change (∆T). A reference
temperature may be input using the TREF or MP,REFT command. For details on how the current
density is calculated, refer to Liquid in the Mechanical APDL Theory Reference.
Gas
Use the TB,FLUID command with TBOPT = GAS to define material behavior for a gas, and specify the
following material constant using the TBDATA command:
You can define a temperature dependent gas material with up to 20 temperatures (NTEMP = 20 max
on the TB command) by using the TBTEMP command, as shown in the example below:
TB,FLUID,1,2,,GAS ! Activate gas material model
TBTEMP,100.0 ! Define first temperature
TBDATA,1, 0f ! Define material constants at temp 100.0
TBTEMP,200.0 ! Define second temperature
TBDATA,1, 0f ! Define material constants at temp 200.0
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
108 of ANSYS, Inc. and its subsidiaries and affiliates.
User-Defined Material Model
When specifying temperature dependent density values for a gas, keep in mind that the current density
(ρf ) for hydrostatic fluid elements is computed at each iteration based on the Ideal Gas Law. For details
on how the current density is calculated, refer to Gas in the Mechanical APDL Theory Reference.
To use the Ideal Gas Law, you also need to define a reference pressure (input as real constant PREF)
and a reference temperature (input with the TREF or MP,REFT command) with temperature offset (input
with the TOFFST command).
Pressure-Volume Data
Use the TB,FLUID command with TBOPT = PVDATA to define compressible fluid behavior in terms of
a pressure-volume curve. You can specify up to 20 temperature-dependent pressure-volume curves
(NTEMP = 20 max on the TB command). The temperature for the first curve is input with TBTEMP, fol-
lowed by TBPT commands for up to 100 pressure-volume data points. The data points (X, Y) entered
on TBPT are:
Con- Meaning
stant
X Pressure value
Y Corresponding volume
value
The pressure-volume data point must be defined in terms of total pressure and total volume of the
fluid in the containing vessel.
Input for the user-defined option is determined by constants which you define. The TB,USER command
initializes the constant table. The constants are defined via TBDATA commands (six per command). The
number of constants can be any combination of the number of temperatures (NTEMP) and the number
of data points per temperature (NPTS), to a maximum of NTEMP x NPTS = 1000.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 109
Chapter 3: Material Models
Example 3.5 Initializing the Values of State Variables for a User-Defined Material
TB,STATE,1,,4, ! Define material 1, which
! has four state variables.
TBDATA,1,C1,C2,C3,C4, ! Initialize the four state variables.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
110 of ANSYS, Inc. and its subsidiaries and affiliates.
Material Strength Limits
The following table summarizes the applicable strength-limit constants for each failure criterion:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 111
Chapter 3: Material Models
Specify the TBOPT value on the TB,DMGI command to correspond to failure criteria (TBOPT = 1 or
FCRT).
The following table shows the coefficient values addressed for the available TBOPT values:
1 -- Maximum strain
2 -- Maximum stress
3 -- Puck
4 -- Hashin
5 -- LaRc03
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
112 of ANSYS, Inc. and its subsidiaries and affiliates.
Damage Evolution Law
6 -- LaRc04
11 -- User-defined #1
...
19 -- User-defined #9
To complete the material damage definition, it is also necessary to specify a compatible material damage
evolution law (TB,DMGE). Without a damage evolution law, the damage-initiation criteria have no effect
on the material. The following table summarizes the compatible damage evolution laws with specific
damage-initiation criteria:
Compatible TB,DMGE,,,,TBOPT
Option
1 or MPDG
TB,DGMI,TBOPT 1 or Y
option FCRT
Specify the TBOPT value on the TB,DMGE command to correspond to the instant stiffness reduction
(TBOPT = 1 or MPDG).
The following table shows the coefficient values addressed for the available TBOPT values:
TBOPT = 1 or
Con- MPDG
stants
NPTS = 4
C1 Tensile fiber stiff-
ness reduction
[1]
C2 Compressive
fiber stiffness re-
duction [1]
C3 Tensile matrix
stiffness reduc-
tion [1]
C4 Compressive
matrix stiffness
reduction [1]
To complete the material damage definition, it is also necessary to specify a compatible material damage
initiation criterion (TB,DMGI). Without a damage initiation criterion, the damage evolution law has no
effect on the material. The following table summarizes the compatible damage-initiation criteria with
specific damage evolution laws:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 113
Chapter 3: Material Models
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
114 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 4: Explicit Dynamics Materials
Material properties used in explicit dynamic analyses (ANSYS LS-DYNA User's Guide program) differ
somewhat from those used in implicit analyses (described in Linear Material Properties (p. 14) and Ma-
terial Models (p. 13).)
Most explicit dynamics material models require data table input. A data table is a series of constants
that are interpreted when they are used. Data tables are always associated with a material number and
are most often used to define nonlinear material data (that is, stress-strain curves). The form of the data
table (referred to as the TB table) depends on the material model being defined.
For a complete description of all explicit dynamics material models, including detailed data table input,
see Material Models in the ANSYS LS-DYNA User's Guide.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 115
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
116 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 5: Material Curve Fitting
Material curve fitting allows you to derive coefficients from experimental data that you provide for your
material. Curve fitting involves comparing your experimental data to certain preexisting nonlinear ma-
terial models. With this capability, you compare experimental data versus program-calculated data for
different nonlinear models; based on those comparisons, you can determine the best material model
to use during solution.
Curve fitting is based on the data table configurations outlined in the TB command. The data manipu-
lations and constructions are performed via the TBFT command.
You perform curve fitting either interactively (GUI) or via batch commands. You input your experimental
data, choose a model from one of nine supplied hyperelastic models, perform a regression analysis,
graphically view the curve-fitting results, compare the fits to the experimental data, and write the fitted
coefficients to the database as nonlinear data table commands for the subsequent finite element analyses.
Hyperelastic models can define three types of behavior: purely incompressible, nearly incompressible,
and compressible. Hyperelastic curve fitting is based on the HYPER option of the TB command.
1 Step 1. Prepare Experimental Data The experimental data must be a plain text file delimited by
a space or a comma.
2 Step 2. Input the Experimental Data The experimental data can be read in by browsing to the file
location (GUI) or by specifying the filename and path (batch)
on the command line.
3 Step 3. Select a Material Model Op- The material options for the applicable curve-fitting regimen
tion are defined in the TB command. Nine hyperelastic models
are supported. After you pick a model, you can still change
to another model if an ideal fit is not realized.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 117
Chapter 5: Material Curve Fitting
4 Step 4. Initialize the Coefficients Depending on the model, hyperelastic curve fitting can be a
linear or nonlinear regression process.The hyperelastic ma-
terial models, along with the associated process for each are
listed in Table 5.3: Hyperelastic Curve-Fitting Model Types.
5 Step 5. Specify Control Parameters You will specify the type of error norm to be used to generate
and Solve the curve fit.
6 Step 6. Plot Your Experimental Data You review and verify the results by comparing the experi-
and Analyze mental data and the regression errors. If the results you obtain
are not acceptable, repeat steps 3 to 5 to perform a new curve-
fitting solution.
7 Step 7.Write Data to the TB Com- Write your curve-fitting results to the database in the TB
mand command table format.
• Case 1 - Totally Incompressible Models (see Table 5.1: Experimental Details for Case 1 and 2 Models and
Blatz-Ko (p. 118))
• Case 2 - Nearly Incompressible Models (see Table 5.1: Experimental Details for Case 1 and 2 Models and
Blatz-Ko (p. 118))
• Case 3 - Compressible Models (see Table 5.2: Experimental Details for Case 3 Models (p. 118))
The types of data required for each of these cases is defined in the tables below.
Table 5.1 Experimental Details for Case 1 and 2 Models and Blatz-Ko
Experimental Type Column 1 Column 2 Column 3
Uniaxial Test Engineering Strain Engineering Stress
Biaxial Test Engineering Strain Engineering Stress
Planar/Shear Test Engineering Strain Engineering Stress
(in loading direc-
tion)
Simple Shear Test Engineering Shear Engineering Shear (Optional) Engineering Nor-
Strain Stress mal Stress (normal to the
edge of shear stress)
Volumetric Test Volume Ratio (J) Hydrostatic Pressure
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
118 of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelastic Material Curve Fitting
All stresses output in POST1/POST26 are true stresses and logarithmic strains.
For temperature-dependent curve fitting, specify your temperature values at the top of the experimental
data using the
/temp,value
line. This header format specifies the attribute (temp) and its value (100). An example of a typical data
input using these attributes follows:
/temp,100
0.9703 60.00
0.9412 118.2
0.9127 175.2
0.8847 231.1
Adding this header introduces a temperature attribute of 100 degrees. You can add additional data
sets at other temperatures, in additional files. One file can have data at only one temperature.
For compressible materials, the curve-fitting tool's default behavior is to solve only for stress as a function
of strain and lateral strain. To force the curve-fitting tool to also fit experimental lateral strain data to
generate the coefficients for the Ogden compressible foam model, add the line /USEL,1 near the top
of the experimental data file. This option is valid for uniaxial, biaxial and shear test data.
5.1.3.1. Batch
TBFT,EADD,ID,Option1,Option2,Option3,Option4
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 119
Chapter 5: Material Curve Fitting
5.1.3.2. GUI
The Material Properties GUI provides an input field where you can type in the filename of your data
file, and also include the appropriate path. You can also browse to a file in a specified location. Separate
input is performed for each Option1 value (UNIA, BIAX, SHEA, etc.).
1. The number of coefficients is usually the sum of the number of deviatoric coefficients and the number
of volumetric coefficients.
2. The number of coefficients for a polynomial will be dependent on the polynomial order N.
Number of Coefficients =∑ + +
=1
Blatz-Ko and Ogden hyper-foam are compressible models. For Ogden hyper-foam, the experimental
data you supply will require additional fields.
For more information about the hyperelastic models available for curve fitting, see Hyperelastic Mater-
ial Models in the Material Reference.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
120 of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelastic Material Curve Fitting
where:
You can also fix (hold constant) your coefficients (Option4 = FIX on the TBFT command). You specify
a value for a coefficient and keep it unchanged, while allowing the other coefficients to be operated
on. You can then release the coefficient you fixed and operate on it after the others are optimized. By
default, all of the coefficients are free to vary. The capability to fix coefficients applies only to nonlinear
curve fits (as listed in Table 5.3: Hyperelastic Curve-Fitting Model Types (p. 120)).
You can also specify tref = all and initiate multiple solves to evaluate coefficients at all available discrete
temperature values. In this case, for data at three temperatures (t1, t2, and t3), a single TBFT,SOLVE
entry will initiate three different solve operations at those three discrete temperatures.
With temperature dependency on and the reference temperature set to a particular value, a TBFT,SOLVE
command solves for coefficients only at that temperature. To solve for coefficients at other temperatures,
you set the reference temperature to each desired discrete temperature value and solve again.
You can initialize the coefficients before or after turning temperature dependency on. If the coefficients
are initialized before turning temperature dependency on, the specified coefficients will become the
initial coefficients for all future solves for that particular model. These coefficients are, however, overridden
when temperature dependency is turned on and another set of initial values are specified at discrete
temperature values. The curve-fitting tool looks for the initial coefficients at a particular temperature.
If no coefficients are specified at a discrete temperature, the initial coefficients set before temperature
dependency was turned on are used.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 121
Chapter 5: Material Curve Fitting
5.1.5.1. Batch
TBFT,SET,ID,HYPER,Option2,Option3,Option4,Option5
where:
where:
where:
5.1.5.2. GUI
The GUI automatically updates your coefficient tables depending on the model you pick. You can
modify individual coefficients to initialize them at values you believe are more appropriate.
1. Number of iterations
2. Residual tolerance
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
122 of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelastic Material Curve Fitting
The solution stops when both the residual tolerance and the coefficient change tolerance of your error
norm are met, or if the number of iterations criteria is met. When you use nonlinear regression, you
must initialize your coefficients.
5.1.6.1. Batch
TBFT,SOLVE,ID,HYPER,Option2,Option3,Option4, ..., Option7
where:
Other solution parameters are available. See the TBFT command for details.
5.1.6.2. GUI
The GUI lets you specify all of your control parameters (error norm, solution control parameters, and
the solver options) interactively. You select the appropriate options from the provided menus, and solve
to generate the coefficients. You can change the parameters and repeat the solution as necessary to
ensure an accurate result. The unused options are disabled whenever necessary.
You should try to reserve column one for the variable that you would like to see vary in the plot. For
example, if your data contains time, temperature, stress, and creep strain, and you wish to see the creep
strain vary as a function of time at different temperatures and stresses in the plot, you add your exper-
imental data using multiple TBFT,EADD commands (or the corresponding GUI operations) by splitting
the file into multiple experimental files as prescribed earlier, one for each combination of temperature
and stress.
5.1.7.1. GUI
Use the GRAPH button to plot the data. Your plots will show columns 2 and above as separate curves,
plotted as a function of column 1. The data in column 1 is always the X-axis. By default, all the experi-
ments are plotted in the GUI window. To view specific data and its corresponding fitting result, you
can click the right mouse button (RMB) on the specific data set, and pick a desired option to view the
results. Other RMB plotting utilities can be found for different data fields in the curve-fitting GUI window.
Use RMB functions to Zoom, Fit, Save Plot to File, View/Hide Legend and View/Hide Grid. Two or more
fitted functions can also be compared in the same plot. For example, you can view Mooney2 Uniaxial
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 123
Chapter 5: Material Curve Fitting
and Mooney9 Uniaxial plots directly on top of each other, provided you have already solved for both
of these functions. RMB also allows you to set the number of points used to generate the plot, and also
change the minimum X value and the maximum X value in a plot. You can also hide a particular curve
within a graph, turn the legend and/or axis displays on and off, and switch the scales between log scale
and regular scale. Use the middle-mouse button to eliminate a specific curve and clarify or refine the
remaining curve.
5.1.7.2. Review/Verify
The two factors you consider in determining results acceptability are visual fit and the error norm/residual
values. When you plot the curve, the error norm/residual values are printed in the curve-fitting GUI
window. Error norm values help you determine the quality of curve fitting and whether to accept the
results, but are not always the best indicator of a valid curve fit. Plotting the curves and visually assessing
the result is usually the best indication. If the results are unacceptable, you may want to go back to
step 3 and solve again, either by picking a different model, increasing the order, or redefining your
initial values of the coefficients or other control parameters. You can continue to use your original ex-
perimental data, repeating step 3 through step 7 until you get an acceptable solution.
5.1.8.1. Batch
TBFT,FSET,ID,HYPER,Option2,Option3
where:
5.1.8.2. GUI
Once you complete the process and update your material data properties with the representative curve
data, you are returned to the material properties dialog. The curve data can now be accessed for the
full range of material behavior.
Curve fitting is performed either interactively or via batch commands. You input the experimental data,
define the order of Prony series expansion, perform nonlinear regression, view the curve-fitting results
graphically, compare the experimental data, and write the fitted coefficients to the database as nonlinear
data table commands for the subsequent finite elements analyses. You can fit shear modulus and/or
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
124 of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelastic Material Curve Fitting
bulk modulus and/or shift functions, along with discrete temperature dependencies for multiple data
sets.
1 Step 1. Prepare Experimental Data The experimental data must be a plain text file delimited by
a space or a comma.
2 Step 2. Input the Data The experimental data can be read in from the GUI or via batch
commands, as a plain text file.
3 Step 3. Select a Material Model Op- This includes Prony series expansion of shear and/or bulk
tion moduli as well as shift function.The supported shift functions
include WLF and TN.
4 Step 4. Initialize the Coefficients Viscoelastic curve fitting is a nonlinear regression; the initial
value of your coefficients is very important for a successful
solution.
5 Step 5. Specify Control Parameters Specify the error norm to be used, the solution control para-
and Solve meters, and perform the nonlinear regression.
6 Step 6. Plot the Experimental Data Graphically view the curve-fitting results. Review and verify
and Analyze the results by comparing them with the experimental data
and the regression errors. If any factor is not acceptable, repeat
steps 3 to 7 to obtain a new curve-fitting solution.
7 Step 7.Write Data to the TB Com- Write your curve-fitting results in TB command format to the
mand database.
Your viscoelastic test data must be a plain text file with headers to define the test data. The data file
should be in table form, delimited by a space or a comma. The header defines the test data type and
the temperature for your test data. For viscoelastic curve fitting with multiple temperatures, you can
evaluate coefficients at each discrete temperature point and write it as a temperature-dependent Prony
data table, or you can use the Williams-Landau-Ferry (WLF) or Tool-Narayanaswamy (TN) shift functions
to account for the temperature dependency. (See Shift Functions in the Mechanical APDL Theory Reference.)
A separate data file must be provided for each discrete temperature. The viscoelastic test data can be
any of 4 data types, see Table 5.4: Viscoelastic Data Types and Abbreviations (p. 125).
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 125
Chapter 5: Material Curve Fitting
The headers are used to describe the data types that characterize the test data columns or attributes
of the data.
The following listing contains the appropriate headers, followed by the delimited data:
/temp,100 ! define temperature attribute
0.01 2992.53
1 2978.514207
2 2965.45541
4 2942.293214
6 2922.530649
8 2905.612202
10 2891.073456
20 2842.506984
40 2798.142793
60 2772.383729
80 2750.631843
100 2730.398114
200 2643.125432
400 2517.475394
600 2431.262053
800 2366.580897
1000 2313.955396
2000 2117.922594
4000 1833.734397
6000 1627.199197
8000 1470.6806
10000 1347.264527
20000 964.0141125
40000 586.1405449
60000 392.186777
80000 277.2706253
100000 202.0025278
200000 46.87056342
400000 2.669209118
600000 0.156653269
800000 0.0137224
1000000 0.005591539
You can include several data sets, including tests performed at different temperatures. Although different
data sets can have the same/or different temperature, each file can have only one temperature. Multiple
temperature data sets must be input with multiple files.
Two types of data may be required for viscoelastic curve fitting, either shear modulus vs. time and/or
bulk modulus vs. time. The data can also be a function of temperature, which can then be accounted
for by either WLF or TN shift function. You can use the GUI or batch command to input your experimental
data.
5.2.3.1. Batch
TBFT,EADD,ID,Option1,Option2,Option3,Option4
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
126 of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelastic Material Curve Fitting
The sdec coefficient refers to the shear modulus as a function of time. The bdec coefficient refers to
the bulk modulus as a function of time.
5.2.3.2. GUI
Click on the Add Dataset button and type the filename into the area provided. You can also browse
to a file in a specified location. Separate input is performed for each data type (Option1 = sdec, or
bdec)
First you define a case name to associate the set of coefficients for the Prony expansions with the shift
functions that characterize the material behavior. You can use the case name to define several different
options that characterize the same test data, and then to compare the curve-fitting results. To define
the material model, you must first define a case name, and then specify the order of shear and bulk
moduli and the type of the shift function(s), if required. You need to create an additional case to define
different shear order, bulk order or shift options. Once you create a case, the number of shear terms,
bulk terms, or shift options cannot be changed.
where:
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 127
Chapter 5: Material Curve Fitting
TBFT,FCASE,ID,FINI
The coefficients are ordered as shear terms first, then the bulk terms, and then the shift function. The
coefficients are ordered as A0G, A1G, τ1G, A2G, τ2G, … AnG, and τnG for shear modulus, and A0K, A1K, τ1K,
A2K, τ2K, … AnK, and τnK for bulk modulus.
A shift function must be used together with your shear and/or bulk modulus for temperature-dependent
experimental data.
The default coefficient is set to 1, but it is good practice to redefine the initial values before solving.
When initializing your coefficients, set AnKs to 1 and τnKs to time values that are equally distributed in
the log scale, spanning the data range from minimum to maximum time.
For example, consider the shear-decay versus time-data file. If the time values vary from 1 to 10000,
and if you use third-order Prony, logical guesses for τ1G , τ2G and τ3G that span this range could be τ1G
= 1, τ2G = 100, and τ3G = 10,000 (also (1), (10) and (10,000) or (1), (1,000) and (10,000), respectively).
The equations below describe the relationship between the Prony Coefficient (αn) and the corresponding
coefficient generated in curve fitting (An). K and G are the shear modulus and bulk modulus at t = 0.
This was done to keep all αnG and αnK values used in the TB tables positive.
A good guess for the WLF or TN parameter is the reference temperature you used during your partial
solve for shear and bulk. The index of the reference/base temperatures is the sum of NumShear +
NumBulk + 1.
You can also fix (hold constant) your coefficients. You specify a value for a coefficient and keep it un-
changed, while allowing the other coefficients to be operated on. You can then release the fixed coef-
ficient later if desired. By default, all of the coefficients are free to vary.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
128 of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelastic Material Curve Fitting
You estimate coefficients for temperature-dependent data either by using the shift function or by setting
the temperature-dependency flag and setting a reference temperature before solving for the coefficients.
You can set the reference temperature only to values specified using the /temp,value header line in
the experimental data.
You can also specify tref = all and initiate multiple solves to evaluate coefficients at all available discrete
temperature values. In this case, for data at three temperatures (t1, t2, and t3), a single TBFT,SOLVE
command initiates three separate solve operations at those three discrete temperature values, and
generate data at three corresponding discrete temperatures.
With temperature dependency specified and the reference temperature set to a particular value, a
TBFT,SOLVE command solves for coefficients only at that discrete temperature. To solve for coefficients
at other temperatures, set the reference temperature to each of the desired discrete temperature values
and solve again.
You can initialize the coefficients before or after activating temperature dependency. If the coefficients
initialize before setting temperature dependency, the specified coefficients become the initial coefficients
for all future solves for that particular model. These coefficients are, however, overridden when temper-
ature dependency is active and another set of values is specified at a discrete temperature value. The
curve-fitting tool looks for the initial coefficients at a particular temperature. If no coefficients are specified
at discrete temperature values, the initial coefficients set before temperature dependency was activated
are used.
where:
Use the TBFT,FIX command to fix a coefficient to a value set by the TBFT, SET command or to release
a previously fixed coefficient. By default, coefficients are not fixed.
TBFT,FIX,ID,CASE,Option2,Option3,Option4,Option5
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 129
Chapter 5: Material Curve Fitting
Temperature dependency uses Option4, and references data files you entered with the appropriate
temp header.
TBFT,SET,ID,VISCO,Option2,Option3,Option4,Option5
where:
The solution control parameters for nonlinear regression include number of iterations, residual tolerance,
and coefficient change tolerance. The solution stops when both residual tolerance of error norm and
coefficient change tolerance is met or if the number of iterations criteria is met. The coefficients are
updated when the solution is completed. In general it is very difficult to directly solve a complete case
including coefficients of the shear modulus, the bulk modulus and the shift function. Three solver options
including shear modulus only, bulk modulus, and shift function (or all) are provided to allow you to
solve only Prony coefficients of the shear modulus, Prony coefficients of the bulk modulus, and coeffi-
cients of the shift function. In many cases, however, the coefficients of shift function can't be solved
until the shear or bulk modulus are solved. It is normal for a solution to not converge at first, but to
stop when the maximum iterations criteria is reached. At that point, you should examine the curve-fitting
results and the solution history before proceeding any further. You can then adjust parameters and
resolve the problem whenever it is necessary.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
130 of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelastic Material Curve Fitting
Set the reference temperature at which your partial solution will be performed using TBFT,
SET,,,,,TREF, TX. Only data at temperature TX will be used to estimate shear coefficients
Solve
2. Solve the bulk coefficients (if there are any):
Set the reference temperature at which your partial solution will be performed using TBFT,
SET,,,,,TREF, TX. The reference temperature should be the same for both shear and bulk. Only
data at temperature TX will be used to estimate shear coefficients
TBFT,SOLVE
3. Solve the shift function (or all) coefficients:
TREF is not used when solving for all parameters. All temperature data is used to estimate
the coefficients.
TBFT,SOLVE
For GUI operations, when only the shear and bulk buttons are checked, only your shear coefficients are
solved. To solve for both shear and bulk, you must check all three buttons.
Set the reference temperature at which your partial solution will be performed using TBFT,
SET,,,,,TREF, TX. Only data at temperature TX will be used to estimate shear coefficients
TBFT,SOLVE
Set the reference temperature at which your partial solution will be performed using TBFT,
SET,,,,,TREF, TX. Only data at temperature TX will be used to estimate the shear coefficients.
TBFT,SOLVE
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 131
Chapter 5: Material Curve Fitting
Alternatively, you can solve for both shear and bulk data at the same time.
1. Set the temperature dependency flag using the command TBFT, SET,,,,,TDEP,1.
2. Solve for ALL coefficients.
Set the reference temperature at which your partial solution will be performed using TBFT,
SET,,,,,TREF, TX. Only data at temperature TX will be used to estimate shear and bulk coefficients
TBFT,SOLVE
where:
The SOLVE option allows you to specify procedure types, tolerances, and number of iterations.
TBFT,SOLVE,ID,CASE,Option2,Option3,Option4, ..., Option7
where:
Other solving parameters are available. See the TBFT command for details.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
132 of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelastic Material Curve Fitting
Right mouse button (RMB) functions allow you to Zoom, Fit, Save Plot to File, View/Hide Legend, and
View/Hide Grid. Two or more fitted functions can also be compared in the same plot. RMB also allows
you to see the number of points used to generate the plot, and also change the minimum X value and
the maximum X value in a plot. You can use the middle-mouse button (context sensitive) to hide a
particular curve within a graph.
Reviewing your curve-fitting result graphically is the only way to ensure a good fit. After plotting the
curve-fitting results, you can then review multiple plots and also verify the error norm/residual value
that is printed in the curve-fitting GUI window. This information helps you determine the quality of a
curve fit and decide whether to accept the results. If not, you may want to go back to step 3, solve
again by changing the order of the Prony series, redefining certain initial values of the coefficients, and
possibly other control parameters. You can continue to use your original data, repeating step 3 through
step 6 until you are satisfied with the solution.
The curve-fitting tool calculates coefficients that differ from those shown in the TB,PRONY table. The
following equations are used to convert the calculated coefficients to the TB,PRONY format:
Curve-Fitting Equation:
− t / τi
= 2
0
+ ∑ 2
i
i =1−N
Prony Equation:
= = α + ∑ α − τ
αi and τi (i = 1 to N) are the Prony coefficients (entered via the TB,PRONY) command. A0 is not entered
in the Prony table.
Conversion Procedure:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 133
Chapter 5: Material Curve Fitting
t =0 = ∑ 2
i
i = 0 −N
α =
=
where:
Option2 = Case
Option3 = Case name
The following topics concerning creep material curve fitting are available:
5.3.1. Understanding the Creep Material Curve-Fitting Process
5.3.2. Step 1. Prepare Experimental Data
5.3.3. Step 2. Input the Experimental Data
5.3.4. Step 3. Select a Material Model Option
5.3.5. Step 4. Initialize the Coefficients
5.3.6. Step 5. Specify Control Parameters and Solve
5.3.7. Step 6. Plot the Experimental Data and Analyze
5.3.8. Step 7. Write Data to the TB Command
5.3.9.Tips For Curve Fitting Creep Models
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
134 of ANSYS, Inc. and its subsidiaries and affiliates.
Creep Material Curve Fitting
The header is used to describe the data types that characterize the test data columns or attributes of
the data. Five different creep data types are available:
The header format to define each column's data type is /n, abbr, where n is the index of the data column
in the file, and abbr is the abbreviation for the type of data in the column, as described in Table 5.5: Creep
Data Types and Abbreviations (p. 135).
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 135
Chapter 5: Material Curve Fitting
When a particular column is unchanged over the loading history, you can define it as an attribute. For
instance, in the above example, the stress and temperature are constant throughout the range. You
define this data as an attribute.
The header format to define a data attribute is /attr, value, where attr is the data-type abbreviation, and
value is the value of the attribute. The constant stress and temperature values above can be written
into the file header, as follows:
/seqv,4000 ! indicate this creep has a constant stress of 4000
/temp,100 ! indicate this creep data is at a constant temperature of 100
/1,creq ! indicate first column is creep strain
/2,dcreq ! indicate second column is creep strain rate
0.00215869 0.000203055
0.00406109 0.000181314
0.00664691 0.000165303
0.0102068 0.000152217
0.0151416 0.000140946
0.0220102 0.000130945
Thirteen model types are available for creep curve fitting. The model you choose determines the exper-
imental data required for the curve-fitting process. The following table describes the creep data required
to perform curve fitting for each model type. Note that for strain hardening and modified strain
hardening, you need to input both creep strain and creep strain rate in the experimental data.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
136 of ANSYS, Inc. and its subsidiaries and affiliates.
Creep Material Curve Fitting
Each file is viewed as a data set, and can be a complete set of experimental test data or a part of a
series of files of experimental test data. You can include several data sets, such as tests performed at
different loading conditions, strain ranges, and temperatures.
Input your experimental data using either the batch method or the GUI method.
where:
See Table 5.7: Creep Models and Abbreviations (p. 138) to determine a starting point for the initial creep
model coefficients.
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 137
Chapter 5: Material Curve Fitting
The following table describes the creep models available and their abbreviated names for Option2
(above).
The experimental data must be consistent with the creep model you choose. See Table 5.6: Creep Model
and Data/Type Attribute (p. 136) for the data types required for each creep model.
Models with many parameters will sometimes converge more easily if you fix (hold constant) your
coefficients. You specify a value for a coefficient and keep it unchanged, while allowing the other
coefficients to be operated on. You can then release the fixed coefficients to obtain a solution. By default,
all of the coefficients are free to vary.
Although some creep models implicitly address temperature dependency, other models lack this cap-
ability. For both types of models, you can account for temperature dependency by generating coefficients
at discrete temperature values.
Solve for separate coefficients at each of your desired discrete temperatures. If the equation implicitly
supports temperature dependency already, that portion is eliminated by appropriately fixing the coef-
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
138 of ANSYS, Inc. and its subsidiaries and affiliates.
Creep Material Curve Fitting
ficient to a specific value. For example, if you have e(-C*T), then C is set to 0 to eliminate this term. To
do this, you activate temperature dependency by setting the tdep parameter to 1.
You can also specify tref = all and initiate multiple solves to evaluate coefficients at all available discrete
temperature values. In this case, for data at three temperatures (t1, t2, and t3), a single TBFT,SOLVE
entry will initiate three different solve operations at those three discrete temperatures.
With temperature dependency on and the reference temperature set to a particular value, a TBFT,SOLVE
command solves for coefficients only at that temperature. To solve for coefficients at other temperatures,
you set the reference temperature to each of the desired discrete temperature values and solve again.
You can initialize the coefficients before or after turning temperature dependency on. If the coefficients
are initialized before turning temperature dependency on, the specified coefficients will become the
initial coefficients for all future solves for that particular model. These coefficients are, however, overridden
when temperature dependency is turned on and another set of initial values is specified at discrete
temperatures. The curve-fitting tool looks for the initial coefficients at a particular temperature. If no
coefficients are specified at discrete temperatures, the initial coefficients set before temperature depend-
ency was turned on are used.
where:
where:
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 139
Chapter 5: Material Curve Fitting
Other available solve criteria are number of iterations, residual tolerance, and coefficient change tolerance.
The solution stops when both residual tolerance and coefficient change tolerance is met or if the
number of iterations criteria is met. The coefficients are updated after every iteration during the solve
process.
where:
Other solving parameters are available. See the TBFT command for details.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
140 of ANSYS, Inc. and its subsidiaries and affiliates.
Creep Material Curve Fitting
You should reserve column one for the variable that you would like to see vary in the plot. For example,
if your data contains time, temperature, stress and creep strain, you may wish to see the creep strain
vary as a function of time, at different temperatures and stresses in the plot. Add your experimental
data using multiple TBFT,EADD commands (or the corresponding GUI method). Split the file into multiple
experimental fields as prescribed earlier, one for each combination of temperature and stress.
Right mouse button (RMB) functions allow you to Zoom, Fit, Save Plot to File, View/Hide Legend and
View/Hide Grid. Two or more fitted functions can also be compared in the same plot. For example, you
can view Mooney2 Uniaxial and Mooney9 Uniaxial plots right on top of each other, provided both of
these function are already solved for. RMB also allows you to see the number of points used to generate
the plot, and also change the Minimum X Value and the Maximum X Value in a plot. You can use the
middle-mouse button (context sensitive) to hide a particular curve within a graph.
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 141
Chapter 5: Material Curve Fitting
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
142 of ANSYS, Inc. and its subsidiaries and affiliates.
Chaboche Material Curve Fitting
set C4 to zero. To keep the Sinh term within floating-point range, keep
c2σ close to one when you initialize the coefficients.
Exponential Form Exponential form has 3 coefficients, with C3 dedicated to temperature
dependency. If you do not have temperature-dependent data, set C3
to zero.To keep eσC2 within floating-point range, keep σ/C2 close to
one.
Norton Norton model has 3 coefficients, with C3 dedicated to temperature
dependency. If you do not have temperature-dependent data, set C3
to zero.
Prim+Sec Time Harden- Time hardening has 7 coefficients.This is a complex model. Here it is
ing advisable to solve for temperature independent data first and then
introduce temperature related data.
Prim+Sec Rational Poly- Rational polynomial is a very complex model for curve fitting, with
nomial 10 coefficients. If you find it hard to fit this data, it is advisable that
you split the experimental data into primary creep data and secondary
creep data. Primary creep data is the initial part of the curve that
covers the nonlinearity in the strain rate. Fit only the secondary data
by fixing C1 to 1 and then set all other coefficients except C2, C3 and
C4 to zero. Use a low value of C3 to keep 10C3σ within floating-point
range. Coefficients C5 to C10 in curve fitting refers to coefficients C7
to C12 in the implicit creep equation (rational polynomial).Then add
the primary creep data, release all coefficients, and solve.
Generalized Time Generalized time hardening has 6 coefficients. Set C6 to zero if you
Hardening have temperature independent data. When initializing coefficients
set C5σ close to 1 to avoid floating-point overflows.
Curve fitting is performed either interactively or via batch commands. You can fit uniaxial plastic strain
vs. stress data, along with discrete temperature dependencies for multiple data sets.
The following topics concerning Chaboche material curve fitting are available:
5.4.1. Understanding the Chaboche Material Curve-Fitting Process
5.4.2. Step 1. Prepare Experimental Data
5.4.3. Step 2. Input the Experimental Data
5.4.4. Step 3. Select a Material Model Option
5.4.5. Step 4. Initialize the Coefficients
5.4.6. Step 5. Specify Control Parameters and Solve
5.4.7. Step 6. Plot the Experimental Data and Analyze
5.4.8. Step 7. Write Data to the TB Command
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 143
Chapter 5: Material Curve Fitting
Your uniaxial test data must be a plain text file with headers to define the test data. The data file should
be in table format, delimited by spaces or commas. Headers can be used to describe the data types
that characterize the test data columns or additional attributes of the data.
For Chaboche curve fitting with multiple temperatures, you can evaluate coefficients at each discrete
temperature point and write it as a temperature-dependent Chaboche data table. A separate data file
is necessary for each discrete temperature.
Issue this command at the top of the experimental data file to specify the temperature for the experiment:
/temp,TempValue
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
144 of ANSYS, Inc. and its subsidiaries and affiliates.
Chaboche Material Curve Fitting
4.42E-03 3.22E+02
4.42E-03 7.78E+01
4.41E-03 -1.65E+02
3.54E-03 -2.31E+02
2.51E-03 -2.66E+02
1.40E-03 -2.84E+02
2.49E-04 -2.95E+02
-9.11E-04 -3.03E+02
-2.08E-03 -3.10E+02
-3.24E-03 -3.17E+02
-4.41E-03 -3.24E+02
-4.41E-03 -7.95E+01
-4.40E-03 1.63E+02
-3.53E-03 2.30E+02
-2.50E-03 2.65E+02
-1.39E-03 2.83E+02
-2.44E-04 2.94E+02
9.16E-04 3.02E+02
2.08E-03 3.09E+02
3.25E-03 3.16E+02
4.41E-03 3.23E+02
4.41E-03 7.87E+01
4.41E-03 -1.64E+02
3.53E-03 -2.31E+02
2.51E-03 -2.66E+02
1.40E-03 -2.84E+02
2.48E-04 -2.95E+02
Uniaxial test data can include loading, unloading, and cyclic loading.
For plasticity, experimental data is path-dependent and the stress-strain behavior depends on the history
of the loading and/or unloading.
Each file is viewed as a data set, and can be a complete set of experimental test data or a part of a
series of files of experimental test data. You can include several data sets, such as tests performed at
different stress levels and/or temperatures, when you perform creep curve fitting.
Input your experimental data using either the batch method or the GUI method.
where:
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 145
Chapter 5: Material Curve Fitting
To define the material model, you must first specify a case name, then specify the order of Chaboche
kinematic model.
After you create a case, the number of Chaboche terms cannot be changed without deleting the case.
Additional cases can be created to model Chaboche kinematic hardening with different orders.
The following syntax examples and argument descriptions illustrate a complete case definition:
TBFT,FCASE,ID,NEW,Option2,Option3 ! define case
where:
where:
The number of parameter depends on the order of the Chaboche model. For a Chaboche model of order
N, there are 2*N+1 coefficients. Coefficients 1 through 2*N are Chaboche parameters, and coefficient
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
146 of ANSYS, Inc. and its subsidiaries and affiliates.
Chaboche Material Curve Fitting
2*N+1 is the yield stress. The odd coefficients in the Chaboche model refer to the slope of the curve
and the even coefficients are decay function parameters. The number of terms depend on the complexity
of the curve. Evenly distributed estimates of the slope over the range of the curve can be used as initial
guesses for the curve-fitting process. Different terms can dominate at different parts of the curve. Initial
guesses for the decay parameters can typically be chosen one or two orders less than odd coefficients
(slope).
You can also fix (hold constant) your coefficients (TBFT,FIX). Specify a value for a coefficient and keep
it unchanged, while allowing the other coefficients to be operated on. You can then release the fixed
coefficient later if desired. By default, all of the coefficients are free to vary.
Estimate coefficients for temperature-dependent data by setting the temperature-dependency flag and
setting a reference temperature before solving for the coefficients. You can set the reference temperature
only to values specified via the /temp,TempValue header line in the experimental data.
You can also specify tref = all and initiate multiple solves to evaluate coefficients at all available
discrete temperature values. In this case, for data at three temperatures (t1, t2, and t3), a single
TBFT,SOLVE command initiates three separate solve operations at those three discrete temperature
values, and generates data at three corresponding discrete temperatures.
With temperature dependency specified and the reference temperature set to a specific value, a
TBFT,SOLVE command solves for coefficients only at that discrete temperature. To solve for coefficients
at other temperatures, set the reference temperature to each of the desired discrete temperature values
and solve again.
You can initialize the coefficients before or after activating temperature dependency. If the coefficients
initialize before setting temperature dependency, the specified coefficients become the initial coefficients
for all future solves for that particular model. These coefficients are, however, overridden when temper-
ature dependency is active and another set of values is specified at a discrete temperature value. The
curve-fitting tool looks for the initial coefficients at a particular temperature. If no coefficients are specified
at discrete temperature values, the initial coefficients set before temperature dependency was activated
are used.
where:
By default, coefficients are not fixed. To fix a coefficient to a value set via the TBFT,SET command, or
to release a previously fixed coefficient, issue the TBFT,FIX command.
TBFT,FIX,ID,CASE,Option2,Option3,Option4,Option5
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 147
Chapter 5: Material Curve Fitting
where:
Temperature dependency uses Option4 and references your specified data files with the appropriate
“temp” header:
TBFT,SET,ID,CASE,Option2,Option3,Option4,Option5
where:
The solution control parameters for nonlinear regression include number of iterations, residual tolerance,
and coefficient change tolerance. The solution stops when both residual tolerance of the error norm
and coefficient change tolerance is met, or if the number of iterations criteria is met. The coefficients
are updated when the solution is completed.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
148 of ANSYS, Inc. and its subsidiaries and affiliates.
Chaboche Material Curve Fitting
c. Solve (TBFT,SOLVE).
where:
where:
The SOLVE option allows you to specify procedure types, tolerances, and the number of iterations:
TBFT,SOLVE,ID,CASE,Option2,Option3,Option4, ... , Option7
where:
Other parameters for solving are available. See the TBFT command for more information.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 149
Chapter 5: Material Curve Fitting
Take advantage of the right-mouse-button (RMB) functions to zoom, fit, save your plot to a file, view
or hide objects, toggle between log scale and regular scale, and so on. With the middle mouse button,
you can eliminate specific curves from each window's display in order to view the remaining data more
clearly.
After plotting the curve-fitting results, you can then review multiple plots and also verify the error
norm/residual value displayed in the curve-fitting GUI window. This information helps you to determine
the quality of a curve fit and decide whether or not to accept the results.
If the curve-fitting results are unsatisfactory, you may want to go back to Step 3. Select a Material
Model Option and solve again by changing the order of the Chaboche model or other options, redefining
certain initial values of the coefficients, and possibly redefining other control parameters. You can
continue to use your original data, repeating step 3 through step 6 until you are satisfied with the
solution.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
150 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 6: Material Model Combinations
You can combine several of the material model options discussed in this chapter to simulate various
material behaviors. Table 6.1: Material Model Combination Possibilities (p. 151) presents the model options
you can combine along with the associated TB command labels, and links to sample input listings located
under Material Model Combinations in the Structural Analysis Guide.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 151
Chapter 6: Material Model Combinations
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
152 of ANSYS, Inc. and its subsidiaries and affiliates.
Model With ... Combina- Command, Label Link to Example
tion Type
Multilinear HILL and PLAS (MISO)
TB,HILL +
Anisotropic Plas- Combined Isotropic and CHAB Example
TB,PLAS,,,,MISO +
ticity Hardening and
TB,CHAB
Chaboche
Nonlinear HILL and NLISO and
Anisotropic Plas- Combined Isotropic TB,HILL + TB,NLISO + CHAB Example
ticity Hardening and TB,CHAB
Chaboche
Anisotropic Vis- Isotropic TB,HILL + TB,RATE + HILL and RATE and
Bilinear
coplasticity Hardening TB,BISO BISO Example
Anisotropic Vis- Isotropic TB,HILL + TB,RATE + HILL and RATE and
Multilinear
coplasticity Hardening TB,MISO MISO Example
Anisotropic Vis- Isotropic TB,HILL + TB,RATE + HILL and RATE and
Nonlinear
coplasticity Hardening TB,NLISO NLISO Example
Anisotropic HILL and CREEP Ex-
--- TB,HILL + TB,CREEP
Creep (Implicit) ample
Anisotropic HILL and CREEP and
Isotropic TB,HILL + TB,CREEP
Creep and Plasti- Bilinear BISO Example
Hardening + TB,BISO
city (Implicit)
Anisotropic HILL and CREEP and
Isotropic TB,HILL + TB,CREEP
Creep and Plasti- Multilinear MISO Example
Hardening + TB,MISO
city (Implicit)
Anisotropic HILL and CREEP and
Isotropic TB,HILL + TB,CREEP
Creep and Plasti- Multilinear PLAS (MISO) Example
Hardening + TB,MISO
city (Implicit)
Anisotropic HILL and CREEP and
Isotropic TB,HILL + TB,CREEP
Creep and Plasti- Nonlinear NLISO Example
Hardening + TB,NLISO
city (Implicit)
Anisotropic HILL and CREEP and
Kinematic
Creep and Plasti- Bilinear TB,HILL + BKIN Example
Hardening
city (Implicit)
Finite Strain AHYPER and PRONY
Anisotropic Hy- Nonlinear
Visco-Aniso- TB,AHYPER + (Anisotropic Hypere-
perelasticity and Anisotropic
tropic Hypere- TB,PRONY lasticity and Vis-
Viscoelasticity Elasticity
lasticity coelasticity) Example
HYPER and PRONY
Hyperelasticity Finite Strain Nonlinear
TB,HYPER + (Hyperelasticity and
and Viscoelasti- Visco-Hypere- Isotropic
TB,PRONY Viscoelasticity) Ex-
city lasticity Elasticity
ample
Extended Druck- EDP and CREEP and
Bilinear, Mul- TB,EDP + TB,CREEP +
er-Prager (EDP) Isotropic PLAS (MISO) Example
tilinear, or TB,BISO or TB,MISO
and Creep (Impli- Hardening
Nonlinear or TB,NLISO
cit)
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 153
Chapter 6: Material Model Combinations
Following are cross-reference links to other sections in the documentation that provide descriptions of
the individual material model options represented in the table above.
Additional information about most these material model options is available under Plastic Material
Models in the Structural Analysis Guide.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
154 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 7: Understanding Field Variable Interpolation
When you define field-dependent material properties via the TBFIELD and related TB commands, the
program uses linear interpolation between the data points you provide to determine specific material
property values. To do so, the program:
TEMPERATURE = 100.00
SLIDE DIST = 0.10000
TEMPERATURE = 100.00
SLIDE DIST = 0.50000
TEMPERATURE = 200.00
SLIDE DIST = 0.20000
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 155
Chapter 7: Understanding Field Variable Interpolation
TEMPERATURE = 200.00
SLIDE DIST = 0.70000
Sliding Distance
Temperature 0.1 0.2 0.5 0.7
100 0.3 0.5
200 0.2 0.1
When defining tabular data, the first specified field variable forms the rows of the table. The subsequent
variables form the columns. In this example, Temperature is the first defined field variable.
In this case, the user defined only four out of a possible eight grid locations.
To populate the interpolation search space, the program fills the missing grid points in each row from
left to right. If the first or subsequent grid locations of a row are not defined, the program uses the first
defined value within the row to backfill the grid. The program then fills any undefined locations within
the grid by linearly interpolating between defined points in each row. If the last value(s) along a row
are not defined, the program gives them the last previously defined value within that row.
Therefore, based on the defined field-dependent friction values, the program generates the following
grid automatically (where values in italics represent those provided by the program):
Sliding Distance
Temperature 0.1 0.2 0.5 0.7
100 0.3 0.35 0.5 0.5
200 0.2 0.2 0.14 0.1
= 0 +α 1− 0 Equation (1)
where
− Equation (2)
α=
−
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
156 of ANSYS, Inc. and its subsidiaries and affiliates.
Example: Two-Dimensional Interpolation
− Equation (3)
α= =
−
and solving for the interpolated values using Equation (1), we obtain
Equation (4)
= + − =
In this example, the program performs the first interpolation at a temperature of 100 and a sliding dis-
tance of 0.4, yielding the result of 0.45 (as shown in Equation (4)).
The program performs the second interpolation for a temperature of 200 and a sliding
distance of 0.4. In this case, we find that
− Equation (5)
α= =
−
and solving for the interpolated values using Equation (1), we obtain
Equation (6)
= + − =
Finally, the program performs a third interpolation between the temperature value of
100 and 200 at a sliding distance of 0.4.
− Equation (7)
α= =
−
and solving for the interpolated values using Equation (1), we obtain
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 157
Chapter 7: Understanding Field Variable Interpolation
Equation (8)
= + − =
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
158 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 8: GUI-Inaccessible Material Properties
The following material properties are not available via the material property menus of the interactive
GUI. You can specify them from the command line or input file. Subsequent graphic display and post-
processing is still displayed.
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 159
Release 14.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
160 of ANSYS, Inc. and its subsidiaries and affiliates.