ANSYS Mechanical APDL Material Reference
ANSYS Mechanical APDL Material Reference
ANSYS Mechanical APDL Material Reference
ANSYS, Inc.
Southpointe
2600 ANSYS Drive
Canonsburg, PA 15317
ansysinfo@ansys.com
http://www.ansys.com
(T) 724-746-3304
(F) 724-514-9494
Release 16.0
January 2015
ANSYS, Inc. is
certified to ISO
9001:2008.
Disclaimer Notice
THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION INCLUDE TRADE SECRETS AND ARE CONFIDENTIAL 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.
Published in the U.S.A.
Table of Contents
1. Introduction to Material Models ............................................................................................................. 1
1.1. Material Models for Displacement Applications ................................................................................. 1
1.2. Material Models for Temperature Applications ................................................................................... 2
1.3. Material Models for Electromagnetic Applications ............................................................................. 2
1.4. Material Models for Coupled Applications ......................................................................................... 3
1.5. Material Parameters .......................................................................................................................... 3
2. Material Model Element Support ........................................................................................................... 5
3. Material Models .................................................................................................................................... 13
3.1. Understanding Material Data Tables ................................................................................................ 13
3.2. Experimental Data .......................................................................................................................... 14
3.3. Linear Material Properties ............................................................................................................... 14
3.3.1. Defining Linear Material Properties ......................................................................................... 15
3.3.2. Stress-Strain Relationships ...................................................................................................... 17
3.3.3. Anisotropic Elasticity .............................................................................................................. 18
3.3.4. Damping ............................................................................................................................... 18
3.3.5. Thermal Expansion ................................................................................................................. 19
3.3.6. Emissivity ............................................................................................................................... 20
3.3.7. Specific Heat .......................................................................................................................... 20
3.3.8. Film Coefficients ..................................................................................................................... 21
3.3.9. Temperature Dependency ...................................................................................................... 21
3.3.10. How Material Properties Are Evaluated ................................................................................. 21
3.4. Rate-Independent Plasticity ............................................................................................................ 21
3.4.1. Understanding the Plasticity Models ....................................................................................... 22
3.4.1.1. Nomenclature ............................................................................................................... 23
3.4.1.2. Strain Decomposition .................................................................................................... 24
3.4.1.3.Yield Criterion ................................................................................................................ 24
3.4.1.4. Flow Rule ...................................................................................................................... 25
3.4.1.5. Hardening ..................................................................................................................... 26
3.4.1.6. Large Deformation ........................................................................................................ 27
3.4.1.7. Output .......................................................................................................................... 27
3.4.1.8. Resources ...................................................................................................................... 28
3.4.2. Isotropic Hardening ............................................................................................................... 29
3.4.2.1. Yield Criteria and Plastic Potentials ................................................................................. 30
3.4.2.1.1. Von Mises Yield Criterion ....................................................................................... 30
3.4.2.1.2. Hill Yield Criterion ................................................................................................. 31
3.4.2.1.2.1. Separated Hill Potentials for Plasticity and Creep ........................................... 33
3.4.2.2. General Isotropic Hardening Classes .............................................................................. 33
3.4.2.2.1. Bilinear Isotropic Hardening .................................................................................. 33
3.4.2.2.1.1. Defining the Bilinear Isotropic Hardening Model ........................................... 34
3.4.2.2.2. Multilinear Isotropic Hardening ............................................................................. 34
3.4.2.2.2.1. Defining the Multilinear Isotropic Hardening Model ...................................... 35
3.4.2.2.3. Nonlinear Isotropic Hardening .............................................................................. 36
3.4.2.2.3.1. Power Law Nonlinear Isotropic Hardening .................................................... 36
3.4.2.2.3.2. Voce Law Nonlinear Isotropic Hardening ....................................................... 36
3.4.2.2.4. Isotropic Hardening Static Recovery ...................................................................... 38
3.4.2.2.4.1. Defining the Isotropic Static Recovery ........................................................... 38
3.4.3. Kinematic Hardening ............................................................................................................. 39
3.4.3.1. Yield Criteria and Plastic Potentials ................................................................................. 40
3.4.3.2. General Kinematic Hardening Classes ............................................................................ 40
3.4.3.2.1. Bilinear Kinematic Hardening ................................................................................ 40
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
iii
Material Reference
3.4.3.2.1.1. Defining the Bilinear Kinematic Hardening Model ......................................... 41
3.4.3.2.2. Multilinear Kinematic Hardening ........................................................................... 41
3.4.3.2.2.1. Defining the Multilinear Kinematic Hardening Model .................................... 43
3.4.3.2.3. Nonlinear Kinematic Hardening ............................................................................ 44
3.4.3.2.3.1. Defining the Nonlinear Kinematic Hardening Model ..................................... 44
3.4.3.2.4. Kinematic Hardening Static Recovery .................................................................... 45
3.4.3.2.4.1. Defining the Kinematic Static Recovery ......................................................... 45
3.4.4. Generalized Hill ...................................................................................................................... 46
3.4.4.1. Defining the Generalized Hill Model ............................................................................... 48
3.4.5. Drucker-Prager ....................................................................................................................... 48
3.4.5.1. Classic Drucker-Prager ................................................................................................... 48
3.4.5.1.1. Defining the Classic Drucker-Prager Model ............................................................ 49
3.4.5.2. Extended Drucker-Prager (EDP) ...................................................................................... 49
3.4.5.2.1. EDP Yield Criteria Forms ........................................................................................ 49
3.4.5.2.1.1. Linear Form .................................................................................................. 50
3.4.5.2.1.2. Power Law Form ........................................................................................... 50
3.4.5.2.1.3. Hyperbolic Form .......................................................................................... 51
3.4.5.2.2. EDP Plastic Flow Potentials .................................................................................... 52
3.4.5.2.2.1. Linear Form .................................................................................................. 52
3.4.5.2.2.2. Power Law Form ........................................................................................... 52
3.4.5.2.2.3. Hyperbolic Form .......................................................................................... 52
3.4.5.2.3. Plastic Strain Increments for Flow Potentials .......................................................... 53
3.4.5.2.4. Example EDP Material Model Definitions ............................................................... 53
3.4.5.3. Extended Drucker-Prager Cap ........................................................................................ 54
3.4.5.3.1. Defining the EDP Cap Yield Criterion and Hardening .............................................. 56
3.4.5.3.2. Defining the EDP Cap Plastic Potential ................................................................... 57
3.4.5.3.3. Example EDP Cap Material Model Definition .......................................................... 57
3.4.6. Gurson ................................................................................................................................... 57
3.4.6.1. Void Volume Fraction ..................................................................................................... 58
3.4.6.2. Hardening ..................................................................................................................... 60
3.4.6.3. Defining the Gurson Material Model .............................................................................. 60
3.4.6.3.1. Defining the Gurson Base Model ........................................................................... 60
3.4.6.3.2. Defining Stress- or Strain-Controlled Nucleation .................................................... 60
3.4.6.3.3. Defining the Void Coalescence Behavior ............................................................... 61
3.4.6.3.4. Example Gurson Model Definition ......................................................................... 61
3.4.7. Cast Iron ................................................................................................................................ 62
3.4.7.1. Defining the Cast Iron Material Model ............................................................................ 63
3.5. Rate-Dependent Plasticity (Viscoplasticity) ...................................................................................... 64
3.5.1. Perzyna and Peirce Options .................................................................................................... 64
3.5.2. Exponential Visco-Hardening (EVH) Option ............................................................................. 65
3.5.3. Anand Option ........................................................................................................................ 66
3.5.4. Defining Rate-Dependent Plasticity (Viscoplasticity) ............................................................... 66
3.5.5. Creep ..................................................................................................................................... 67
3.5.5.1. Implicit Creep Equations ................................................................................................ 68
3.5.5.2. Explicit Creep Equations ................................................................................................ 69
3.5.5.2.1. Primary Explicit Creep Equation for C6 = 0 ............................................................. 70
3.5.5.2.2. Primary Explicit Creep Equation for C6 = 1 ............................................................. 71
3.5.5.2.3. Primary Explicit Creep Equation for C6 = 2 ............................................................. 71
3.5.5.2.4. Primary Explicit Creep Equation for C6 = 9 ............................................................. 71
3.5.5.2.4.1. Double Exponential Creep Equation (C4 = 0) ................................................. 71
3.5.5.2.4.2. Rational Polynomial Creep Equation with Metric Units (C4 = 1) ...................... 71
3.5.5.2.4.3. Rational Polynomial Creep Equation with English Units (C4 = 2) .................... 72
iv
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Reference
3.5.5.2.5. Primary Explicit Creep Equation for C6 = 10 ........................................................... 72
3.5.5.2.5.1. Double Exponential Creep Equation (C4 = 0) ................................................. 72
3.5.5.2.5.2. Rational Polynomial Creep Equation with Metric Units (C4 = 1) ...................... 72
3.5.5.2.5.3. Rational Polynomial Creep Equation with English Units (C4 = 2) .................... 73
3.5.5.2.6. Primary Explicit Creep Equation for C6 = 11 ........................................................... 73
3.5.5.2.6.1. Modified Rational Polynomial Creep Equation (C4 = 0) .................................. 73
3.5.5.2.6.2. Rational Polynomial Creep Equation with Metric Units (C4 = 1) ...................... 73
3.5.5.2.6.3. Rational Polynomial Creep Equation with English Units (C4 = 2) .................... 74
3.5.5.2.7. Primary Explicit Creep Equation for C6 = 12 ........................................................... 74
3.5.5.2.8. Primary Explicit Creep Equation for C6 Equals 13 ................................................... 75
3.5.5.2.9. Primary Explicit Creep Equation for C6 = 14 ........................................................... 75
3.5.5.2.10. Primary Explicit Creep Equation for C6 = 15 ......................................................... 76
3.5.5.2.11. Primary Explicit Creep Equation for C6 = 100 ....................................................... 76
3.5.5.2.12. Secondary Explicit Creep Equation for C12 = 0 ..................................................... 76
3.5.5.2.13. Secondary Explicit Creep Equation for C12 = 1 ..................................................... 76
3.5.5.2.14. Irradiation Induced Explicit Creep Equation for C66 = 5 ........................................ 76
3.6. Hyperelasticity ................................................................................................................................ 77
3.6.1. Arruda-Boyce Hyperelasticity .................................................................................................. 77
3.6.2. Blatz-Ko Foam Hyperelasticity ................................................................................................. 78
3.6.3. Extended Tube Hyperelasticity ............................................................................................... 78
3.6.4. Gent Hyperelasticity ............................................................................................................... 79
3.6.5. Mooney-Rivlin Hyperelasticity ................................................................................................ 79
3.6.6. Neo-Hookean Hyperelasticity ................................................................................................. 81
3.6.7. Ogden Hyperelasticity ............................................................................................................ 81
3.6.8. Ogden Compressible Foam Hyperelasticity ............................................................................. 82
3.6.9. Polynomial Form Hyperelasticity ............................................................................................. 83
3.6.10. Response Function Hyperelasticity ....................................................................................... 85
3.6.11. Yeoh Hyperelasticity ............................................................................................................. 85
3.6.12. Special Hyperelasticity .......................................................................................................... 86
3.6.12.1. Anisotropic Hyperelasticity .......................................................................................... 86
3.6.12.2. Bergstrom-Boyce Material ............................................................................................ 87
3.6.12.3. Mullins Effect ............................................................................................................... 88
3.6.12.4. User-Defined Hyperelastic Material .............................................................................. 89
3.7. Viscoelasticity ................................................................................................................................. 89
3.7.1. Viscoelastic Formulation ......................................................................................................... 90
3.7.1.1. Small Deformation ......................................................................................................... 90
3.7.1.2. Small Strain with Large Deformation .............................................................................. 92
3.7.1.3. Large Deformation ........................................................................................................ 92
3.7.1.4. Dissipation .................................................................................................................... 93
3.7.2. Time-Temperature Superposition ........................................................................................... 94
3.7.2.1. Williams-Landel-Ferry Shift Function .............................................................................. 94
3.7.2.2. Tool-Narayanaswamy Shift Function ............................................................................... 95
3.7.2.3. User-Defined Shift Function ........................................................................................... 96
3.7.3. Harmonic Viscoelasticity ......................................................................................................... 97
3.7.3.1. Prony Series Complex Modulus ...................................................................................... 97
3.7.3.2. Experimental Data Complex Modulus ............................................................................ 97
3.7.3.3. Frequency-Temperature Superposition .......................................................................... 99
3.7.3.4. Stress ............................................................................................................................ 99
3.8. Microplane ..................................................................................................................................... 99
3.8.1. Microplane Modeling ............................................................................................................. 99
3.8.1.1. Discretization .............................................................................................................. 101
3.8.2. Material Models with Degradation and Damage .................................................................... 101
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Reference
3.8.3. Material Parameters Definition and Example Input ................................................................ 103
3.8.4. Learning More About Microplane Material Modeling ............................................................. 104
3.9. Porous Media ................................................................................................................................ 104
3.9.1. Coupled Pore-Fluid Diffusion and Structural Model of Porous Media ...................................... 104
3.10. Electricity and Magnetism ........................................................................................................... 105
3.10.1. Piezoelectricity ................................................................................................................... 105
3.10.2. Piezoresistivity ................................................................................................................... 106
3.10.3. Magnetism ......................................................................................................................... 107
3.10.4. Anisotropic Electric Permittivity .......................................................................................... 107
3.11. Gasket ........................................................................................................................................ 108
3.12. Swelling ...................................................................................................................................... 110
3.13. Shape Memory Alloy (SMA) ......................................................................................................... 111
3.13.1. SMA Model for Superelasticity ............................................................................................ 112
3.13.1.1. Constitutive Model for Superelasticity ........................................................................ 113
3.13.1.2. Material Parameters for the Superelastic SMA Material Model ..................................... 115
3.13.2. SMA Material Model with Shape Memory Effect .................................................................. 115
3.13.2.1. The Constitutive Model for Shape Memory Effect ........................................................ 116
3.13.2.2. Material Parameters for the Shape Memory Effect Option ........................................... 118
3.13.3. Result Output of Solution Variables ..................................................................................... 119
3.13.4. Element Support for SMA ................................................................................................... 119
3.13.5. Learning More About Shape Memory Alloy ......................................................................... 120
3.14. MPC184 Joint .............................................................................................................................. 120
3.14.1. Linear Elastic Stiffness and Damping Behavior ..................................................................... 120
3.14.2. Nonlinear Elastic Stiffness and Damping Behavior ............................................................... 121
3.14.2.1. Specifying a Function Describing Nonlinear Stiffness Behavior .................................... 123
3.14.3. Frictional Behavior .............................................................................................................. 123
3.15. Contact Friction .......................................................................................................................... 125
3.15.1. Isotropic Friction ................................................................................................................ 126
3.15.2. Orthotropic Friction ............................................................................................................ 126
3.15.3. Redefining Friction Between Load Steps ............................................................................. 127
3.15.4. User-Defined Friction .......................................................................................................... 128
3.16. Contact Interaction ..................................................................................................................... 128
3.16.1. Interaction Options for General Contact Definitions ............................................................ 128
3.16.2. User-Defined Interaction .................................................................................................... 129
3.17. Cohesive Material Law ................................................................................................................. 130
3.17.1. Exponential Cohesive Zone Material for Interface Elements and Contact Elements ............... 130
3.17.2. Bilinear Cohesive Zone Material for Interface Elements and Contact Elements ...................... 130
3.17.3. Viscous Regularization of Cohesive Zone Material for Interface Elements and Contact Elements ........................................................................................................................................... 131
3.17.4. Cohesive Zone Material for Contact Elements ...................................................................... 132
3.17.5. Post-Debonding Behavior at the Contact Interface .............................................................. 133
3.18. Contact Surface Wear .................................................................................................................. 134
3.18.1. Archard Wear Model ........................................................................................................... 134
3.18.2. User-Defined Wear Model ................................................................................................... 135
3.19. Acoustics .................................................................................................................................... 135
3.19.1. Equivalent Fluid Model of Perforated Media ........................................................................ 135
3.19.1.1. Johnson-Champoux-Allard Equivalent Fluid Model of Perforated Media ...................... 135
3.19.1.2. Delany-Bazley Equivalent Fluid Model of Perforated Media ......................................... 136
3.19.1.3. Miki Equivalent Fluid Model of Perforated Media ........................................................ 137
3.19.1.4. Complex Impedance and Propagating-Constant Equivalent Fluid Model of Perforated
Media ..................................................................................................................................... 138
3.19.1.5. Complex Density and Velocity Equivalent Fluid Model of Perforated Media ................. 138
vi
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Reference
3.19.1.6. Transfer Admittance Matrix Model of Perforated Media ............................................... 139
3.19.1.7. Transfer Admittance Matrix Model of a Square or Hexagonal Grid Structure ................ 140
3.19.2. Acoustic Frequency-Dependent Materials ........................................................................... 141
3.19.3. Low Reduced Frequency (LRF) Model of Acoustic Viscous-Thermal Media ............................ 141
3.19.3.1. Thin Layer .................................................................................................................. 141
3.19.3.2.Tube with Rectangular Cross-Section .......................................................................... 141
3.19.3.3. Tube with Circular Cross-Section ................................................................................ 142
3.20. Fluids .......................................................................................................................................... 142
3.21. Custom Material Models .............................................................................................................. 143
3.21.1. User-Defined Material Model (UserMat) .............................................................................. 143
3.21.2. User-Defined Thermal Material Model (UserMatTh) ............................................................. 144
3.21.3. User-Defined Cohesive Material (UserCZM) ......................................................................... 144
3.21.4. Using State Variables with User-Defined Materials ............................................................... 145
3.21.4.1. Using State Variables with the UserMat Subroutine ..................................................... 145
3.21.4.2. Using State Variables with the UserMatTh Subroutine ................................................. 145
3.21.4.3. Using State Variables with the UserCZM Subroutine .................................................... 146
3.22. Material Strength Limits .............................................................................................................. 146
3.23. Material Damage ........................................................................................................................ 148
3.23.1. Damage Initiation Criteria ................................................................................................... 148
3.23.2. Damage Evolution Law ....................................................................................................... 149
3.23.2.1. Predicting Post-Damage Degradation of Brittle Anisotropic Materials ......................... 151
3.23.2.1.1. Damage Modes ................................................................................................. 153
4. Explicit Dynamics Materials ................................................................................................................ 155
5. Material Curve Fitting ......................................................................................................................... 157
5.1. Hyperelastic Material Curve Fitting ................................................................................................ 157
5.1.1. Understanding the Hyperelastic Material Curve-Fitting Process ............................................. 157
5.1.2. Step 1. Prepare Experimental Data ........................................................................................ 158
5.1.3. Step 2. Input the Experimental Data ...................................................................................... 159
5.1.3.1. Batch ........................................................................................................................... 159
5.1.3.2. GUI .............................................................................................................................. 160
5.1.4. Step 3. Select a Material Model Option .................................................................................. 160
5.1.4.1. Batch Method .............................................................................................................. 161
5.1.4.2. GUI Method ................................................................................................................. 161
5.1.5. Step 4. Initialize the Coefficients ............................................................................................ 161
5.1.5.1. Batch ........................................................................................................................... 162
5.1.5.2. GUI .............................................................................................................................. 162
5.1.6. Step 5. Specify Control Parameters and Solve ........................................................................ 162
5.1.6.1. Batch ........................................................................................................................... 163
5.1.6.2. GUI .............................................................................................................................. 163
5.1.7. Step 6. Plot Your Experimental Data and Analyze ................................................................... 163
5.1.7.1. GUI .............................................................................................................................. 163
5.1.7.2. Review/Verify .............................................................................................................. 164
5.1.8. Step 7. Write Data to the TB Command .................................................................................. 164
5.1.8.1. Batch ........................................................................................................................... 164
5.1.8.2. GUI .............................................................................................................................. 164
5.2. Viscoelastic Material Curve Fitting ................................................................................................. 164
5.2.1. Understanding the Viscoelastic Material Curve-Fitting Process .............................................. 165
5.2.2. Step 1. Prepare Experimental Data ........................................................................................ 165
5.2.3. Step 2. Input the Data ........................................................................................................... 166
5.2.3.1. Batch ........................................................................................................................... 167
5.2.3.2. GUI .............................................................................................................................. 167
5.2.4. Step 3. Select a Material Model Option .................................................................................. 167
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
vii
Material Reference
5.2.4.1. Batch Method .............................................................................................................. 167
5.2.4.2. GUI Method ................................................................................................................. 168
5.2.5. Step 4. Initialize the Coefficients ............................................................................................ 168
5.2.5.1. Batch Method ............................................................................................................. 169
5.2.5.2. GUI Method ................................................................................................................. 170
5.2.6. Step 5. Specify Control Parameters and Solve ........................................................................ 170
5.2.6.1. Temperature-Dependent Solutions Using the Shift Function ......................................... 171
5.2.6.2. Temperature-Dependent Solutions Without the Shift Function ..................................... 171
5.2.6.3. Batch Method .............................................................................................................. 172
5.2.6.4. GUI Method ................................................................................................................. 173
5.2.7. Step 6. Plot the Experimental Data and Analyze ..................................................................... 173
5.2.7.1. Analyze Your Curves for Proper Fit ................................................................................ 173
5.2.8. Step 7. Write Data to the TB Command .................................................................................. 173
5.2.8.1. Batch Method .............................................................................................................. 174
5.2.8.2. GUI Method ................................................................................................................. 175
5.3. Creep Material Curve Fitting .......................................................................................................... 175
5.3.1. Understanding the Creep Material Curve-Fitting Process ....................................................... 175
5.3.2. Step 1. Prepare Experimental Data ........................................................................................ 176
5.3.3. Step 2. Input the Experimental Data ...................................................................................... 177
5.3.3.1. Batch Method .............................................................................................................. 177
5.3.3.2. GUI Method ................................................................................................................. 178
5.3.4. Step 3. Select a Material Model Option .................................................................................. 178
5.3.4.1. Batch Method .............................................................................................................. 178
5.3.4.2. GUI Method ................................................................................................................. 179
5.3.5. Step 4. Initialize the Coefficients ............................................................................................ 179
5.3.5.1. Batch Method .............................................................................................................. 180
5.3.5.2. GUI Method ................................................................................................................. 180
5.3.6. Step 5. Specify Control Parameters and Solve ........................................................................ 180
5.3.6.1. Batch Method .............................................................................................................. 181
5.3.6.2. GUI Method ................................................................................................................. 181
5.3.7. Step 6. Plot the Experimental Data and Analyze ..................................................................... 181
5.3.7.1. GUI Method ................................................................................................................. 181
5.3.7.2. Analyze Your Curves for Proper Fit ................................................................................ 182
5.3.8. Step 7. Write Data to the TB Command .................................................................................. 182
5.3.8.1. Batch Method .............................................................................................................. 182
5.3.8.2. GUI Method ................................................................................................................. 182
5.3.9. Tips For Curve Fitting Creep Models ...................................................................................... 182
5.4. Chaboche Material Curve Fitting ................................................................................................... 184
5.4.1. Understanding the Chaboche Material Curve-Fitting Process ................................................ 184
5.4.2. Step 1. Prepare Experimental Data ........................................................................................ 185
5.4.3. Step 2. Input the Experimental Data ...................................................................................... 186
5.4.3.1. Batch Method .............................................................................................................. 186
5.4.3.2. GUI Method ................................................................................................................. 186
5.4.4. Step 3. Select a Material Model Option .................................................................................. 186
5.4.4.1. Batch Method .............................................................................................................. 187
5.4.4.2. GUI Method ................................................................................................................. 187
5.4.5. Step 4. Initialize the Coefficients ............................................................................................ 187
5.4.5.1. Including Isotropic Hardening Models with Chaboche Kinematic Hardening ................. 187
5.4.5.2. General Process for Initializing MISO Option Coefficients .............................................. 188
5.4.5.2.1. Batch Method ..................................................................................................... 189
5.4.5.2.2. GUI Method ........................................................................................................ 189
5.4.6. Step 5. Specify Control Parameters and Solve ........................................................................ 190
viii
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Reference
5.4.6.1. Temperature-Dependent Solutions .............................................................................. 190
5.4.6.2. Batch Method .............................................................................................................. 190
5.4.6.3. GUI Method ................................................................................................................. 191
5.4.7. Step 6. Plot the Experimental Data and Analyze ..................................................................... 191
5.4.7.1. Analyzing Your Curves for Proper Fit ............................................................................. 191
5.4.8. Step 7. Write Data to the TB Command .................................................................................. 191
6. Material Model Combinations ............................................................................................................ 193
7. Understanding Field Variables ............................................................................................................ 199
7.1. User-Defined Field Variables .......................................................................................................... 199
7.1.1. Subroutine for Editing Field Variables .................................................................................... 200
7.2. Data Processing ............................................................................................................................ 200
7.3. Logarithmic Interpolation and Scaling ........................................................................................... 202
7.4. Example: One-Dimensional Interpolation ....................................................................................... 202
7.5. Example: Two-Dimensional Interpolation ....................................................................................... 203
7.6. Example: Multi-Dimensional Interpolation ..................................................................................... 203
8. GUI-Inaccessible Material Properties .................................................................................................. 205
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
ix
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
List of Figures
3.1. Stress-Strain Curve for an Elastic-Plastic Material .................................................................................... 22
3.2.Yield Surface in Principal Stress Space .................................................................................................... 25
3.3. Plastic Strain Flow Rule .......................................................................................................................... 25
3.4. Isotropic Hardening of the Yield Surface ................................................................................................ 26
3.5. Kinematic Hardening of the Yield Surface .............................................................................................. 27
3.6.Yield Surface for von Mises Yield Criterion .............................................................................................. 31
3.7. Stress vs. Total Strain for Bilinear Isotropic Hardening ............................................................................. 34
3.8. Stress vs. Total Strain for Multilinear Isotropic Hardening ........................................................................ 35
3.9. Stress vs. Plastic Strain for Voce Hardening ............................................................................................. 37
3.10. Stress vs. Total Strain for Bilinear Kinematic Hardening ......................................................................... 40
3.11. Stress vs. Total Strain for Multilinear Kinematic Hardening .................................................................... 42
3.12. Power Law Criterion in the Meridian Plane ........................................................................................... 50
3.13. Hyperbolic and Linear Criterion in the Meridian Plane .......................................................................... 51
3.14. Yield Surface for the Cap Criterion ....................................................................................................... 55
3.15. Growth, Nucleation, and Coalescence of Voids at Microscopic Scale ...................................................... 58
3.16. Cast Iron Yield Surfaces for Compression and Tension .......................................................................... 63
3.17. Generalized Maxwell Solid in One Dimension ...................................................................................... 90
3.18. Sphere Discretization by 42 Microplanes ............................................................................................ 101
3.19. Damage Parameter d Depending on the Equivalent Strain Energy ...................................................... 103
3.20. Stress-strain Behavior at Uniaxial Tension ........................................................................................... 103
3.21. Pseudoelasticity (PE) and Shape Memory Effect (SME) ........................................................................ 112
3.22. Typical Superelasticity Behavior ......................................................................................................... 113
3.23. Idealized Stress-Strain Diagram of Superelastic Behavior .................................................................... 114
3.24. Admissible Paths for Elastic Behavior and Phase Transformations ....................................................... 118
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
xi
xii
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
List of Tables
3.1. Linear Material Property Descriptions .................................................................................................... 15
3.2. Implicit Creep Equations ....................................................................................................................... 68
3.3. Superelastic Option Constants ............................................................................................................ 115
3.4. Shape Memory Effect Option Constants .............................................................................................. 119
5.1. Experimental Details for Case 1 and 2 Models and Blatz-Ko .................................................................. 158
5.2. Experimental Details for Case 3 Models ............................................................................................... 158
5.3. Hyperelastic Curve-Fitting Model Types ............................................................................................... 160
5.4. Viscoelastic Data Types and Abbreviations ........................................................................................... 166
5.5. Creep Data Types and Abbreviations ................................................................................................... 176
5.6. Creep Model and Data/Type Attribute ................................................................................................. 177
5.7. Creep Models and Abbreviations ......................................................................................................... 178
6.1. Material Model Combination Possibilities ............................................................................................ 193
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
xiii
xiv
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Behavior
Application
Linear elastic
Plastic and
elastic-plastic
Plastic deformation is
observed in many materials
such as metals, alloys, soils,
rocks, concrete, and
ceramics.
Hyperelastic
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Behavior
Application
deformation, such as
polymers and biological
materials.
Expansion and
swelling
Radiation environments,
bonded materials with
thermal strain mismatch,
and soils that absorb water.
Interaction
These models produce a response based on the Gasket and joint materials
interaction of structures.
and also models of bonded
and separating surfaces
along interfaces or material
cleavage.
Shape memory
alloy
Description
Magnetic
Conductivity
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Parameters
Type
Description
Permittivity and
Permeability
Loss
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Model
Elements
AHYPER
Anisotropic
hyperelasticity
ANEL
ANISO
BB
Bergstrom-Boyce
BH
Magnetic
BISO
Bilinear isotropic
hardening
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Model
Elements
SOLSH190, SHELL208, SHELL209, PLANE223,
SOLID226, SOLID227, REINF263, REINF264, REINF265,
SOLID272, SOLID273, SHELL281, SOLID285, PIPE288,
PIPE289, ELBOW290
BKIN
Bilinear kinematic
hardening
CAST
Cast iron
CDM
Mullins effect
CGCR
Crack growth
fracture criterion
PLANE182, SOLID185
COMP
Composite damage
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Label
(Lab)
Material Model
Elements
CONCR
Concrete or
concrete damage
SOLID65
Concrete damage model using explicit dynamic elements
SOLID164, SOLID168
CREEP
Creep
CTE
Coefficient of
thermal expansion
CZM
Cohesive zone
DISCRETE
Explicit
spring-damper
(discrete)
COMBI165
DMGE
Damage evolution
law
DMGI
Damage initiation
criteria
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Model
Elements
DP
Drucker-Prager
plasticity
SOLID65
DPER
Anisotropic electric
permittivity
EDP
Extended
Drucker-Prager
ELASTIC
Elasticity
EOS
Equation of state
EVISC
Viscoelasticity
EXPE
Experimental data
FCON
Fluid conductance
data
FLUID116
FCLI
Failure criteria
material strength
limits
FLUID
Fluid
HSFLD241, HSFLD242
FOAM
Foam
FRIC
Coefficient of
friction
GASKET
Gasket
GCAP
Geological cap
GURSON
Gurson
pressure-dependent
plasticity
HFLM
FLUID116
HILL
Hill anisotropy
HONEY
Honeycomb
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Label
(Lab)
Material Model
Elements
HYPER
Hyperelasticity
INTER
Contact interaction
JOIN
MPC184
KINH
Multilinear
kinematic hardening
MELAS
Multilinear elasticity
SOLID65
Nonlinear legacy elements only
MISO
Multilinear isotropic
hardening
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Model
Elements
Hill plasticity:
SOLID65, LINK180, SHELL181, PLANE182,
PLANE183, SOLID185, SOLID186, SOLID187,
BEAM188, BEAM189, SOLSH190, SHELL208,
SHELL209, PLANE223, SOLID226, SOLID227,
REINF263, REINF264, REINF265, SOLID272,
SOLID273, SHELL281, SOLID285, ELBOW290
MKIN
Multilinear
kinematic hardening
MOONEY
Mooney-Rivlin
hyperelasticity
MPLANE
Microplane
NLISO
Voce isotropic
hardening law
PERF
Johnson-Champoux-Allard
FLUID30, FLUID220, FLUID221
Equivalent Fluid
Model of a Porous
Media
PIEZ
Piezoelectric matrix
PLASTIC
Plasticity
10
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Label
(Lab)
Material Model
Elements
SOLID226, SOLID227, REINF263, REINF264, REINF265, SOLID272,
SOLID273, SHELL281, SOLID285
PLAW
Plasticity laws
PM
Coupled Pore-Fluid
Diffusion and
Structural Model of
Porous Media
PRONY
Prony series
constants for
viscoelastic
materials
PZRS
Piezoresistivity
RATE
Rate-dependent
plasticity
(viscoplasticity)
SDAMP
Material structural
damping
SHIFT
SMA
STATE
State variables
(user-defined)
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
11
Material Model
Elements
Also, user-defined plasticity or viscoplasticity: PLANE183
SWELL
Swelling
UNIAXIAL
Uniaxial stress-strain
relation
USER
User-defined
WEAR
12
Contact surface
wear
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
13
Material Models
If the form applies to more than one element, it is described here and referenced in the element description.
14
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Units
EX
EY
Description
Elastic modulus, element x direction
Force/Area
EZ
PRXY
PRYZ
PRXZ
NUXY
None
NUYZ
NUXZ
GXY
GYZ
Force/Area
GXZ
ALPX
ALPY
Strain/Temp
ALPZ
CTEX
CTEY
Strain/Temp
CTEZ
THSX
THSY
Strain
THSZ
REFT
Temp
MU
None
ALPD
None
BETD
None
15
Material Models
MP,
Lab
Value
Units
Description
DMPR
None
DENS
Mass/Vol
Mass density
KXX
KYY
KZZ
Heat/Mass*Temp
Specific heat
ENTH
Heat/Vol
HF
Heat /
(Time*Area*Temp)
EMIS
None
Emissivity
Heat/Time
None
VISC
Force*Time/ Length2
Viscosity
SONC
Length/Time
QRATE
MURX
MURY
MURZ
MGXX
MGYY
Current/Length
MGZZ
RSVX
RSVY
Resistance*Area/Length
RSVZ
PERX
PERY
None
PERZ
LSST
SBKX
SBKY
SBKZ
DXX
DYY
Length /Time
DZZ
CREF
CSAT
16
Mass/Length
Saturated concentration
Mass/Length
Reference concentration
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Units
BETX
BETY
Description
Coefficient of diffusion expansion, element x direction
Length /Mass
BETZ
17
Material Models
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:
Constant
Meaning
C1-C6
C7-C12
C13-C18
C19-C21
For a list of the elements that support this material model, see Material Model Element Support (p. 5).
3.3.4. 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,
18
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
19
Material Models
T = element evaluation temperature
TREF = temperature at which zero thermal strains exist (TREF command or REFT )
se(T) = secant coefficient of thermal expansion, with respect to a definition temperature (in
this case, same as TREF) (ALPX )
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:
Tn = temperature at which an se value is being evaluated
To = definition temperature at which the se values are defined (in this case, same as TREF)
in(T) = instantaneous coefficient of thermal expansion at temperature T (CTEX )
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:
where:
ith(T) = thermal strain at temperature T (THSX)
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.
If the se values are based upon a definition temperature other than TREF, then you need to convert
those values to TREF (MPAMOD). Also see the Mechanical APDL Theory Reference.
3.3.6. 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.
20
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
21
Material Models
Figure 3.1: Stress-Strain Curve for an Elastic-Plastic Material
Unloading recovers the elastic portion of the total strain, and if the load is completely removed, a permanent deformation due to the plastic strain remains in the material. Evolution of the plastic strain
depends on the load history such as temperature, stress, and strain rate, as well as internal variables
such as yield strength, back stress, and damage.
To simulate elastic-plastic material behavior, several constitutive models for plasticity are provided. The
models range from simple to complex. The choice of constitutive model generally depends on the experimental data available to fit the material constants.
The following rate-independent plasticity material model topics are available:
3.4.1. Understanding the Plasticity Models
3.4.2. Isotropic Hardening
3.4.3. Kinematic Hardening
3.4.4. Generalized Hill
3.4.5. Drucker-Prager
3.4.6. Gurson
3.4.7. Cast Iron
22
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
3.4.1.2. Strain Decomposition
3.4.1.3.Yield Criterion
3.4.1.4. Flow Rule
3.4.1.5. Hardening
3.4.1.6. Large Deformation
3.4.1.7. Output
3.4.1.8. Resources
3.4.1.1. Nomenclature
Following are the common symbols used in the rate-independent plasticity theory documentation:
Symbol
Definition
Symbol
Definition
Identity tensor
Strain
Young's Modulus
Elastic strain
Elasto-plastic tangent
Plastic strain
Plastic tangent
Stress
Stress components
Principal stresses
Yield stress
Plastic work
Mohr-Coulomb cohesion
Yield function
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
23
Material Models
Symbol
Definition
Symbol
Definition
Plastic potential
For a general model of plasticity that includes arbitrary load paths, the flow theory of plasticity decomposes the incremental strain tensor into elastic and plastic strain increments:
The increment in stress is then proportional to the increment in elastic strain, and the plastic constitutive
model gives the incremental plastic strain as a function of the material state and load increment.
Equation 3.1 (p. 24) is a general function representing the specific form of the yield criterion for each
of the plasticity models. The function is a surface in stress space and an example, plotted in principal
stress space, as shown in this figure:
24
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
Figure 3.2: Yield Surface in Principal Stress Space
where
When the plastic potential is the yield surface from Equation 3.1 (p. 24), the plastic strain increment is
normal to the yield surface and the model has an associated flow rule, as shown in this figure:
Figure 3.3: Plastic Strain Flow Rule
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
25
Material Models
These flow rules are typically used to model metals and give a plastic strain increment that is proportional to the stress increment. If the plastic potential is not proportional to the yield surface, the model
has a non-associated flow rule, typically used to model soils and granular materials that plastically deform
due to internal frictional sliding. For non-associated flow rules, the plastic strain increment is not in the
same direction as the stress increment.
Non-associated flow rules result in an unsymmetric material stiffness tensor. Unsymmetric analysis can
be set via the NROPT command. For a plastic potential that is similar to the yield surface, the plastic
strain direction is not significantly different from the yield surface normal, and the degree of asymmetry
in the material stiffness is small. In this case, a symmetric analysis can be used, and a symmetric material stiffness tensor will be formed without significantly affecting the convergence of the solution.
3.4.1.5. Hardening
The yield criterion for many materials depends on the history of loading and evolution of plastic strain.
The change in the yield criterion due to loading is called hardening and is defined by the hardening
rule. Hardening behavior results in an increase in yield stress upon further loading from a state on the
yield surface so that for a plastically deforming material, an increase in stress is accompanied by an increase in plastic strain.
Two common types of hardening rules are isotropic and kinematic hardening. For isotropic hardening,
the yield surface given by Equation 3.1 (p. 24) has the form:
where
This type of hardening can model the behavior of materials under monotonic loading and elastic unloading, but often does not give good results for structures that experience plastic deformation after
a load reversal from a plastic state.
For kinematic hardening, the yield surface has the form:
26
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
where
As shown in the following figure, the back stress tensor is the center (or origin) of the yield surface,
and plastic loading from
to
results in a change in the back stress and therefore a shift in
the yield surface:
Figure 3.5: Kinematic Hardening of the Yield Surface
Kinematic hardening is observed in cyclic loading of metals. It can be used to model behavior such as
the Bauschinger effect, where the compressive yield strength reduces in response to tensile yielding.
It can also be used to model plastic ratcheting, which is the buildup of plastic strain during cyclic
loading.
Many materials exhibit both isotropic and kinematic hardening behavior, and these hardening rules
can be used together to give the combined hardening model. Other hardening behaviors include
changes in the shape of the yield surface in which the hardening rule affects only a local region of the
yield surface, and softening behavior in which the yield stress decreases with plastic loading.
3.4.1.7. Output
Output quantities specific to the plastic constitutive models are available for use in the POST1 database
postprocessor (/POST1) and in the POST26 time-history results postprocessor (/POST26).
The equivalent stress (label SEPL) is the current value of the yield stress evaluated from the hardening
model.
The accumulated plastic strain (label EPEQ) is a path-dependent summation of the plastic strain rate
over the history of the deformation:
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
27
Material Models
where
The stress ratio (label SRAT) is the ratio of the elastic trial stress to the current yield stress and is an indicator of plastic deformation during an increment. If the stress ratio is:
>1
A plastic deformation occurred during the increment.
<1
An elastic deformation occurred during the increment.
1
The stress state is on the yield surface.
Alternatively, the output quantities can have specialized meanings specific to the given material model.
For example:
For the kinematic hardening models, the equivalent stress is determined by evaluating the kinematic
hardening rule for stress vs. plastic strain using the accumulated plastic strain. The stress ratio for multilinear
kinematic hardening uses the trial stress and yield stress of the first subvolume.
For the generalized Hill model, the equivalent plastic strain is given by:
3.4.1.8. Resources
The following list of resources offers more information about plasticity:
1. Hill, R. The Mathematical Theory of Plasticity. New York: Oxford University Press, 1983.
2. Prager, W.The Theory of Plasticity: A Survey of Recent Achievements. Proceedings of the Institution of
Mechanical Engineers. 169.1 (1955): 41-57.
3. Besseling, J. F.A Theory of Elastic, Plastic, and Creep Deformations of an Initially Isotropic Material Showing
Anisotropic Strain-Hardening, Creep Recovery, and Secondary Creep. ASME Journal of Applied Mechanics.
25 (1958): 529-536.
4. Owen, D. R. J., A.Prakash, O. C. Zienkiewicz.Finite Element Analysis of Non-Linear Composite Materials by
Use of Overlay Systems. Computers and Structures. 4.6 (1974): 1251-1267.
5. Rice, J. R.Continuum mechanics and thermodynamics of plasticity in relation to microscale deformation
mechanisms. Constitutive Equations in Plasticity. Ed. A. Argon. Cambridge, MA: MIT Press, 1975. 23-79.
28
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
6. Chaboche, J. L.Constitutive Equations for Cyclic Plasticity and Cyclic Viscoplasticity. International Journal
of Plasticity. 5.3 (1989): 247-302.
7. Chaboche, J. L.On Some Modifications of Kinematic Hardening to Improve the Description of Ratchetting
Effects. International Journal of Plasticity. 7.7 (1991): 661-678.
8. Shih, C. F., D. Lee.Further Developments in Anisotropic Plasticity. Journal of Engineering Materials and
Technology. 100.3 (1978): 294-302.
9. Valliappan, S., P. Boonlaulohr, I. K. Lee.Non-Linear Analysis for Anisotropic Materials. International Journal
for Numerical Methods in Engineering. 10.3 (1976): 597-606.
10. Drucker, D. C., W. Prager.Soil Mechanics and Plastic Analysis or Limit Design. Quarterly of Applied Mathematics. 10.2 (1952): 157-165.
11. Sandler, I. S, F. L. DiMaggio, G. Y. Baladi.A Generalized Cap Model for Geological Materials. Journal of the
Geotechnical Engineering Division. 102.7 (1976): 683-699.
12. Schwer, L. E., Y. D. Murray.A Three-Invariant Smooth Cap Model with Mixed Hardening. International
Journal for Numerical and Analytical Methods in Geomechanics. 18.10 (1994): 657-688.
13. Foster, C., R. Regueiro, A. Fossum, R. Borja.Implicit Numerical Integration of a Three-Invariant, Isotropic/Kinematic Hardening Cap Plasticity Model for Geomaterials. Computer Methods in Applied Mechanics and Engineering. 194.50-52 (2005): 5109-5138.
14. Pelessone, D.A Modified Formulation of the Cap Model. Technical Report GA-C19579. San Diego: Gulf
Atomics, 1989.
15. Fossum, A.F., J. T. Fredrich.Cap Plasticity Models and Compactive and Dilatant Pre-Failure Deformation.
Pacific Rocks 2000: Rock Around the Rim. Proceedings of the Fourth North American Rock Mechanics Symposium. Eds. J. Girard, M. Liebman, C. Breeds, T. Doe, A. A. Balkema. Rotterdam, 2000: 1169-1176.
16. Gurson, A. L.Continuum Theory of Ductile Rupture by Void Nucleation and Growth: Part I--Yield Criteria
and Flow Rules for Porous Ductile Media. Journal of Engineering Materials and Technology. 99.1 (1977): 215.
17. Needleman, A. V. Tvergaard.An Analysis of Ductile Rupture in Notched Bars. Journal of the Mechanics and
Physics of Solids. 32.6 (1984): 461-490.
18. Arndt, S., B. Svendsen, D. Klingbeil. Modellierung der Eigenspannungen an der Rispitze mit einem
Schdigungsmodell. Technische Mechanik. 4.17 (1997): 323-332.
19. Besson, J., C. Guillemer-Neel.An Extension of the Green and Gurson Models to Kinematic Hardening.
Mechanics of Materials. 35.1-2 (2003): 1-18.
20. Hjelm, H. E.Yield Surface for Grey Cast Iron Under Biaxial Stress. Journal of Engineering Materials and
Technology. 116.2 (1994): 148-154.
21. Chen, W. F., D. J. Han. Plasticity for Structural Engineers. New York: Springer-Verlag, 1988.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
29
Material Models
where
is a scalar function of stress and
is the yield stress that evolves as a function of the
set of material internal variables . This type of hardening can model the behavior of materials under
monotonic loading and elastic unloading, but often does not give good results for structures that experience additional plastic deformation after a load reversal from a plastic state.
Three general classes of isotropic hardening models are available: bilinear, multilinear, and nonlinear.
Each of the hardening models assumes a von Mises yield criterion, unless an anisotropic Hill yield criterion
is defined, and includes an associated flow rule.
Isotropic hardening can also be combined with kinematic hardening and the Extended Drucker-Prager
and Gurson models to provide an evolution of the yield stress. For more information, see Material
Model Combinations (p. 193).
The following topics related to the isotropic hardening material model are available:
3.4.2.1.Yield Criteria and Plastic Potentials
3.4.2.2. General Isotropic Hardening Classes
and
is the von Mises effective stress, also known as the von Mises equivalent stress,
is the yield strength and corresponds to the yield in uniaxial stress loading.
In principal stress space, the yield surface is a cylinder with the axis along the hydrostatic line
and gives a yield criterion that is independent of the hydrostatic stress, as shown in the following figure:
30
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
Figure 3.6: Yield Surface for von Mises Yield Criterion
For an associated flow rule, the plastic potential is the yield criterion in Equation 3.2 (p. 30) and the
plastic strain increment is proportional to the deviatoric stress
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
31
Material Models
where the directional yield stress ratios are the user input parameters and are related to the isotropic
yield stress parameter by:
where is the yield stress in the direction indicated by the value of subscript i. The stress directions
are in the anisotropy coordinate system which is aligned with the element coordinate system (ESYS).
The isotropic yield stress
is entered in the constants for the hardening model.
The Hill yield criterion defines a surface in six-dimensional stress space and the flow direction is normal
to the surface. The plastic strain increments in the anisotropy coordinate system are:
The Hill surface, used with a hardening model, replaces the default von Mises yield surface.
After defining the material data table (TB,HILL), input the required constants (TBDATA):
Constant
Meaning
Property
C1
R11
C2
R22
C3
R33
C4
R12
C5
R23
C6
R13
The constants can be defined as a function of temperature (NTEMP on the TB command), with temperatures specified for the table entries (TBTEMP).
Example 3.1: Hill Surface
/prep7
MP,EX,1,20.0E5
MP,NUXY,1,0.3
32
! ELASTIC CONSTANTS
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
TB,HILL,1
! HILL TABLE
TBDATA,1,1.0,1.1,0.9,0.85,0.9,0.80
! ELASTIC CONSTANTS
TB,HILL,1,1,,PC
! HILL data table for plasticity and creep
TBDATA,1,1.0,1.1,0.9,0.85,0.9,0.80 ! Plasticity Hill parameters
TBDATA,7,1.0,1.0,1.0,1.0 ,1.0,1.0 ! Creep Hill parameters
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
33
Material Models
Figure 3.7: Stress vs. Total Strain for Bilinear Isotropic Hardening
Meaning
Property
C1
Yield stress
C2
Tangent modulus
The constants can be defined as a function of temperature (NTEMP on the TB command), with temperatures specified for the table entries (TBTEMP).
Example 3.3: Bilinear Isotropic Hardening
/prep7
MPTEMP,1,0,500
MPDATA,EX,1,,14E6,12e6
MPDATA,PRXY,1,,0.3,0.3
TB,BISO,1,2
TBTEMP,0.0
TBDATA,1,44E3,1.2E6
TBTEMP,500
TBDATA,1,29.33E3,0.8E6
!
!
!
!
!
34
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
The multilinear hardening behavior is described by a piece-wise linear stress-total strain curve, starting
at the origin and defined by sets of positive stress and strain values, as shown in this figure:
Figure 3.8: Stress vs. Total Strain for Multilinear Isotropic Hardening
The first stress-strain point corresponds to the yield stress. Subsequent points define the elastic plastic
response of the material.
Meaning
Property
Strain value
Stress value
The stress-plastic strain data points are entered into the table via the TBPT command.
Temperature-dependent data can be defined (NTEMP on the TB command), with temperatures specified
for the table entries (TBTEMP). Interpolation between temperatures occurs via stress-vs.-plastic-strain.
Example 3.4: Multilinear Hardening with Plastic Strain
/prep7
MPTEMP,1,0,500
! Define temperature-dependent EX,
MPDATA,EX,1,,14.665E6,12.423e6
MPDATA,PRXY,1,,0.3
TB,PLASTIC,1,2,5,MISO
TBTEMP,0.0
TBPT,DEFI,0,29.33E3
TBPT,DEFI,1.59E-3,50E3
TBPT,DEFI,3.25E-3,55E3
TBPT,DEFI,5.91E-3,60E3
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
35
Material Models
TBPT,DEFI,1.06E-2,65E3
TBTEMP,500
TBPT,DEFI,0,27.33E3
TBPT,DEFI,2.02E-3,37E3
TBPT,DEFI,3.76E-3,40.3E3
TBPT,DEFI,6.48E-3,43.7E3
TBPT,DEFI,1.12E-2,47E3
! Temperature = 500
! Plastic strain, stress at temperature = 500
where G is the shear modulus determined from the user defined elastic constants and
mulated equivalent plastic strain.
is the accu-
Meaning
C1
C2
Property
Initial yield stress
Exponent
36
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
Figure 3.9: Stress vs. Plastic Strain for Voce Hardening
The evolution of the yield stress for this model is specified by the following equation:
Meaning
Property
C1
C2
Linear coefficient
C3
Exponential coefficient
C4
Exponential saturation
parameter
The constants can be defined as a function of temperature (NTEMP on the TB command), with temperatures specified for the subsequent set of constants (TBTEMP).
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
37
Material Models
Example 3.6: Voce Nonlinear Isotropic Hardening
/PREP7
TB,NLISO,1,2,,VOCE
TBTEMP,40
TBDATA,1,280,7e3,155,7e2
TBTEMP,60
TBDATA,1,250,5e3,120,3e2
!
!
!
!
!
where the first term on the right side of the equation represents the rate of work hardening, and the
second defines the rate of static recovery where , and
are material parameters, is the temperature, and
is the current yield stress. A lower limit yield stress can be defined for static recovery such
that:
where
An approximate solution for the static recovery of the yield stress over a time step is to integrate the
static recovery rate and incorporate the work hardening rate into the initial condition. The yield stress
is then given by the solution to:
where
is the yield stress at the end of the previous increment and
stress due to work hardening over the current time step.
38
Meaning
Property
C1
Coefficient
C2
Temperature Coefficient
C3
Exponential
C4
Recovery threshold
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
The constants can be defined as a function of temperature (NTEMP on the TB command), with temperatures specified for the table entries (TBTEMP).
Example 3.7: Isotropic Hardening Static Recovery
/prep7
YOUNGS = 30e3
! Youngs Modulus
NU
= 0.3
! Poissons ratio
SIGMA0 = 18.0E0 ! Initial yield stress
mp,ex,1,YOUNGS
mp,prxy,1,NU
TB,CHAB,1,1,1
! Chaboche kinematic hardening
TBDATA,1,sigma0,
TBDATA,2,1e+3,0
TB,CREEP,1,1, ,1,1
tbdata,1,
! creep model
! No creep strain
where
and
is the shift in the position of the yield surface in stress space and evolves during
Three general classes of kinematic hardening models are available: bilinear, multilinear, and nonlinear.
Each of the hardening models assumes a von Mises yield criterion, unless an anisotropic Hill yield criterion
is defined, and includes an associated flow rule.
Kinematic hardening can also be combined with isotropic hardening and the Gurson model to provide
an evolution of the yield stress. For more information, see Material Model Combinations (p. 193).
The following topics related to the kinematic hardening material model are available:
3.4.3.1.Yield Criteria and Plastic Potentials
3.4.3.2. General Kinematic Hardening Classes
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
39
Material Models
40
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
where G is the elastic shear modulus and the shift strain is numerically integrated from the incremental
shift strain which is proportional to the incremental plastic strain:
where
Meaning
Property
C1
C2
Tangent modulus
The constants can be defined as a function of temperature (NTEMP on the TB command), with temperatures specified for the table entries (TBTEMP).
This model can be used with the TB commands TBOPT option:
For TBOPT 1, no stress relaxation occurs with an increase in temperature. This option is not recommended for non-isothermal problems.
For TBOPT = 1, Rice's hardening rule [5] is applied, which takes stress relaxation with temperature increase into account.
Example 3.8: Bilinear Kinematic Hardening
/prep7
MPTEMP,1,0,500
! Define temperatures for Young's modulus
MPDATA,EX,1,,14E6,12e6
MPDATA,PRXY,1,,0.3,0.3
TB,BKIN,1,2,2,1
! Activate a data table with TBOPT=1
! stress relaxation with temperature
TBTEMP,0.0
! Temperature = 0.0
TBDATA,1,44E3,1.2E6
! Yield = 44,000; Tangent modulus = 1.2E6
TBTEMP,500
! Temperature = 500
TBDATA,1,29.33E3,0.8E6
! Yield = 29,330; Tangent modulus = 0.8E6
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
41
Material Models
Figure 3.11: Stress vs. Total Strain for Multilinear Kinematic Hardening
The model formulation is the sublayer or overlay model of Besselling [3] and Owen, Prakash and Zienkiewicz [4] in which the material is assumed to be composed of a number of sublayers or subvolumes,
all subjected to the same total strain. The number of subvolumes is the same as the number of input
stress-strain points, and the overall behavior is weighted for each subvolume where the weight is given
by:
where
The behavior of each subvolume is elastic-perfectly plastic, with the uniaxial yield stress for each subvolume given by:
where
42
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
The subvolumes undergo kinematic hardening with an associated flow rule and the plastic strain increment for each subvolume is the same as that for bilinear kinematic hardening. The total plastic strain
is given by:
where
Meaning
Property
P1
Strain value
P2
Stress value
The constants can be defined as a function of temperature (NTEMP on the TB command), with temperatures specified for the table entries (TBTEMP).
When entering temperature-dependent stress-strain points, the set of data at each temperature must
have the same number of points. Thermal softening for the multilinear kinematic hardening model is
the same as that for bilinear kinematic hardening with Rice's hardening rule [5].
!
!
!
!
!
!
!
!
!
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
43
Material Models
Example 3.10: Multilinear Kinematic Hardening with Stress vs. Plastic Strain
/prep7
TB,PLASTIC,1,2,3,KINH
TBTEMP,20.0
TBPT,,0.0,1.0
TBPT,,0.1,1.2
TBPT,,0.2,1.3
TBTEMP,40.0
TBPT,,0.0,0.9
TBPT,,0.0900,1.0
TBPT,,0.129,1.05
!
!
!
!
!
!
!
!
!
Stress = 1.0
Stress = 1.2
Stress = 1.3
Stress = 0.9
Stress = 1.0
Stress = 1.05
where
and are user-input material parameters,
of the plastic strain rate.
is the magnitude
During a solution, if there is a change in temperature over an increment (non-isothermal loading), the
back stress terms are scaled in a manner similar to that of bilinear kinematic hardening with Rice's
hardening rule [5].
44
Meaning
Property
C1
C2
C3
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
Constant
Meaning
Property
C4
C5
...
...
...
C(2n)
C(1+2n)
Temperature-dependent data can be defined (NTEMP on the TB command), with temperatures specified
for the table entries (TBTEMP).
Set the TB commands NPTS value equal to n, the number of superimposed kinematic hardening
models.
Example 3.11: Nonlinear Kinematic Hardening
/prep7
TB,CHABOCHE,1,1,3
where the last term on the right side of the equation is the rate of static recovery of the kinematic back
stress component with
and
the material parameters for kinematic static recovery, and
is the
von Mises effective back stress component.
Meaning
Property
C1
Coefficient
C2
Temperature Coefficient
The constants can be defined as a function of temperature (NTEMP on the TB command), with temperatures specified for the table entries (TBTEMP).
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
45
Material Models
Example 3.12: Kinematic Hardening Static Recovery
/prep7
YOUNGS = 30e3
! Youngs Modulus
NU
= 0.3
! Poissons ratio
SIGMA0 = 18.0E0 ! Initial yield stress
mp,ex,1,YOUNGS
mp,prxy,1,NU
TB,CHAB,1,1,1
! Chaboche kinematic hardening
TBDATA,1,sigma0,
TBDATA,2,1e+4,0
TB,CREEP,1,1, ,1,1
tbdata,1,
! creep model
! No creep strain
TB,PLASTIC,1,,1,KSR
TBDATA,1,1e-3,2.0,
where
and
are the user-defined magnitudes of the tension and compression yield strength, respectively. The subscripts on the tension and compression yield stresses correspond to the Voigt notation
coordinate directions
.
From the assumption of incompressible plastic deformation, the mixed subscript coefficients are given
by:
(3.7)
are:
(3.8)
46
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
, then:
and the coefficients in the yield criterion from Equation 3.5 (p. 46) are determined from Equation 3.6 (p. 46) through Equation 3.8 (p. 46), and the user-input tension and compression yield stresses.
Due to the incompressibility assumption,
(3.9)
and for a closed yield surface,
(3.10)
Equation 3.9 (p. 47) and Equation 3.10 (p. 47) must be satisfied throughout the evolution of yield stresses
that result from plastic deformation. The program checks these conditions through 20 percent equivalent
plastic strain, but you must ensure that conditions are satisfied if the deformation exceeds that range.
A bilinear anisotropic work hardening rule is used to evolve the individual components of tension and
compression yield stresses. For a general state of deformation with a bilinear hardening law, the plastic
work is:
where
stress
by:
Equation 3.11 (p. 47) is then the isotropic hardening evolution equation for the tensile and compressive
yield stress components.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
47
Material Models
Meaning
Property
C1-C3
C4-C6
C7-C9
C10-C12
C13-C15
C16-C18
3.4.5. Drucker-Prager
The following topics related to Drucker-Prager plasticity are available:
3.4.5.1. Classic Drucker-Prager
3.4.5.2. Extended Drucker-Prager (EDP)
3.4.5.3. Extended Drucker-Prager Cap
48
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
Angle of internal friction
Dilatancy angle
The amount of dilatancy (the increase in material volume due to yielding) can be controlled via 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 non-associated.
For more information about this material model, see Classic Drucker-Prager Model in the Mechanical
APDL Theory Reference.
Meaning
Property
C1
Force/Area
Cohesion value
C2
Internal friction
C3
Dilatancy
49
Material Models
Meaning
Property
C1
Pressure sensitivity
C2
The constants can be defined as a function of temperature (NTEMP on the TB command), with temperatures specified for the table entries (TBTEMP).
50
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
Constant
Meaning
Property
C1
Pressure sensitivity
C2
Exponent
C3
The constants can be defined as a function of temperature (NTEMP on the TB command), with temperatures specified for the table entries (TBTEMP).
In the following figure, the hyperbolic yield criterion is plotted and compared to the linear yield criterion
shown in the dashed line:
Figure 3.13: Hyperbolic and Linear Criterion in the Meridian Plane
Meaning
Property
C1
Pressure sensitivity
C2
Material parameter
C3
The constants can be defined as a function of temperature (NTEMP on the TB command), with temperatures specified for the table entries (TBTEMP).
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
51
Material Models
where
Meaning
C1
Property
Pressure sensitivity
The material behavior can be defined as a function of temperature (NTEMP on the TB command), with
temperatures specified for the table entries (TBTEMP).
Meaning
Property
C1
Pressure sensitivity
C2
Exponent
The material behavior can be defined as a function of temperature (NTEMP on the TB command), with
temperatures specified for the table entries (TBTEMP).
52
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
the constant
Meaning
Property
C1
Pressure sensitivity
C2
Material parameter
The material behavior can be defined as a function of temperature (NTEMP on the TB command), with
temperatures specified for the table entries (TBTEMP).
where
Associated flow is obtained if the plastic potential form and parameters are set equal to the yield criterion.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
53
Material Models
tb,edp,1,1,2,LFPOT
tbdata,1,0.566206
Example 3.16: EDP -- Power Law Yield Criterion and Flow Potential
/prep7
!!! Define linear elasticity constants
mp,ex,1,2.1e4
mp,nuxy,1,0.45
! Extended Drucker-Prager Material Model Definition
! Power Law Yield Function
tb,edp,1,1,3,PYFUN
tbdata,1,8.33,1.5
! Power Law Plastic Flow Potential
tb,edp,1,1,2,PFPOT
tbdata,1,8.33,1.5
where
, and
, given by:
Three functions define the surfaces that make up the yield criterion. The shear envelope function is
given by:
where
is the cohesion related yield parameter and is a user defined material parameter along with
, , and . This function reduces to the Drucker-Prager criterion for
. For positive values of ,
the shear failure envelope is evaluated at
54
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
The compaction function is itself a function of the shear envelope function and is given by:
where
is the user-defined value of
in the following figure:
, as shown
The expansion function is a function of the shear envelope function and is given by:
where
when
is a user-input material parameter. The expansion function defines the material yield surface
. The expansion cap function reaches peak value at
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
55
Material Models
where is a user-defined material parameter, a ratio of the extension strength to compression strength
in triaxial loading.
Two methods of isotropic hardening can be used to evolve the yield criterion due to plastic deformation.
Hardening of the compaction cap is due to evolution of , which is the intersection of the cap surface
with
shown in Figure 3.14: Yield Surface for the Cap Criterion (p. 55). This value evolves due to
plastic volume strain
where
, and
. The restriction
Evolution of the yield surface at the intersection of the shear envelope with the expansion cap occurs
by combining the cap model with an isotropic hardening model to evolve the value of . The bilinear,
multilinear, or nonlinear isotropic hardening function can be used, and the yield stress from the isotropic
hardening model must be consistent with the value of
calculated from the cap material parameters
given by
.
The following topics related to defining the EDP Cap material model are available:
3.4.5.3.1. Defining the EDP Cap Yield Criterion and Hardening
3.4.5.3.2. Defining the EDP Cap Plastic Potential
3.4.5.3.3. Example EDP Cap Material Model Definition
Material
Property
C1
C2
C3
C4
C5
C6
C7
C8
C9
C10
Hardening parameter
C11
Hardening parameter
The yield criterion and hardening behavior can be defined as a function of temperature (NTEMP on the
TB command), with temperatures specified for the table entries (TBTEMP).
56
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
Material
Property
C1
C2
C3
C4
The plastic flow potential can be defined as a function of temperature (NTEMP on the TB command),
with temperatures specified for the table entries (TBTEMP).
If the plastic flow potential is not defined, the yield surface is used as the flow potential, resulting in
an associated flow model.
3.4.6. Gurson
The Gurson model is used to represent plasticity and damage in ductile porous metals [16][17]. When
plasticity and damage occur, ductile metal undergoes a process of void growth, nucleation, and coalescence. The model incorporates these microscopic material behaviors into macroscopic plasticity behavior based on changes in the void volume fraction, also known as porosity, and pressure. A porosity
increase corresponds to an increase in material damage, resulting in a diminished load-carrying capacity.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
57
Material Models
The yield criterion and flow potential for the Gurson model is:
where
, and
(a): Existing voids grow when the solid matrix is in a hydrostatic-tension state. The solid
matrix is assumed to be incompressible in plasticity so that any material volume growth
is due to the void volume expansion.
(b): Void nucleation occurs, where new voids are created during plastic deformation due
to debonding of the inclusion-matrix or particle-matrix interface, or from the fracture of
the inclusions or particles themselves.
(c): Voids coalesce. In this process, the isolated voids establish connections. Although
coalescence may not discernibly affect the void volume, the load-carrying capacity of
the material begins to decay more rapidly at this stage.
The void volume fraction is the ratio of void volume to the total volume. A volume fraction of 0 indicates
no voids and the yield criterion reduces to the von Mises criterion. A volume fraction of 1 indicates all
58
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
the material is void. The initial void volume fraction,
change of void volume fraction is a combination of the rate of growth and the rate of nucleation:
From the assumption of isochoric plasticity and conservation of mass, the rate of change of void volume
fraction due to growth is proportional to the rate of volumetric plastic strain:
Void nucleation is controlled by either the plastic strain or the stress, and is assumed to follow a normal
distribution of statistics.
In the case of strain-controlled nucleation, the distribution is described by the mean strain,
deviation, . The void nucleation rate due to strain control is given by:
where
, and
In the case of stress-controlled nucleation, void nucleation is determined by the distribution of maximum
normal stress on the interfaces between inclusions and the matrix, equal to
. Stress-controlled
nucleation takes into account the effect of triaxial loading on the rate of void nucleation. The voidnucleation rate for stress control is given by:
and deviation,
The modified void volume fraction, , is used to model the loss of material load carrying capacity associated with void coalescence. When the current void volume fraction reaches a critical value , the
material load carrying capacity decreases rapidly due to coalescence. When the void volume fraction
reaches , the load-carrying capacity of the material is lost completely. The modified void volume
fraction is given by:
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
59
Material Models
3.4.6.2. Hardening
The Gurson model can be combined with one of the isotropic hardening models to incorporate isotropic
hardening of the yield stress in the Gurson yield criterion.
To combine the Gurson model with Chaboche kinematic hardening, the yield criterion is modified to:
where
is the von Mises equivalent modified relative stress, and
are functions of the modified back stress given by:
where
is the kinematic hardening back stress. Then, the modified relative stress is:
Meaning
Property
C1
C2
Initial porosity
C3
C4
C5
60
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
To define stress-controlled nucleation, initialize the material data table (TB,GURSON,,,,SSNU), then input
the following constants (TBDATA):
Constant
Meaning
Property
C1
Nucleation porosity
C2
Mean stress
C3
To define strain-controlled nucleation, initialize the material data table (TB,GURSON,,,,SNNU), then input
the following constants (TBDATA):
Constant
Meaning
Property
C1
Nucleation porosity
C2
Mean strain
C3
Meaning
Property
C1
Critical porosity
C2
Failure porosity
!
!
!
!
!
!
!
!
!
!
!
!
initial porosity
first Tvergaard constant
second Tvergaard constant
third Tvergaard constant = q1^2
critical porosity
failure porosity
nucleation porosity
standard deviation of mean strain
mean strain
initial yielding strength (MPa)
power value for nonlinear isotropic
hardening power law
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
61
Material Models
!base model
tb,gurson,1,,5,base
tbdata,1,sigma_Y,f_0,q1,q2,q3
! Strain-controlled nucleation
tb,gurson,1,,3,snnu
tbdata,1,f_N,strain_N,s_N
! Coalescence
tb,gurson,1,,2,coal
tbdata,1,f_c,f_F
! Power law isotropic hardening
tb,nliso,1,,2,POWER
tbdata,1,sigma_Y,power_N
where
where
and
where
The following figure shows the yield surfaces in compression and tension for perfectly plastic behavior:
62
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rate-Independent Plasticity
Figure 3.16: Cast Iron Yield Surfaces for Compression and Tension
The yield surfaces are plotted in the meridian plane in which the ordinate and abscissa are von Mises
equivalent stress and pressure, respectively.
The evolution of the yield stress in tension and compression follows the user input piecewise linear
stress-strain curves for compression and tension. The tension yield stress evolves as a function of the
equivalent uniaxial plastic strain, . The evolution of the equivalent uniaxial plastic strain is defined
by equating the uniaxial plastic work increment to the total plastic work increment:
The compression yield stress evolves as a function of the equivalent plastic strain, , which is calculated
from the increment in plastic strain determined by consistency with the yield criterion and the flow
potential.
The plastic flow potential is defined by the von Mises yield criterion in compression and results in an
associated flow rule. The flow potential in compression is:
In tension, the Rankine cap yield surface is replaced by an ellipsoidal surface defined by:
where
The plastic Poisson's ratio determines the amount of volumetric expansion during tensile plastic deformation. The tensile flow potential gives a nonassociated flow model and results in an unsymmetric material stiffness tensor.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
63
Material Models
Enter the tensile multilinear hardening stress-strain points into a data table (TB,UNIAXIAL,,,,TENSION).
Do the same for the compressive multilinear hardening stress-strain points (TB,UNIAXIAL,,,,COMPRESSION).
Enter tension and compression stress-strain points into their respective tables (TBPT), with the compression points entered as positive values.
The plastic Poisson's ratio and stress-strain points can be defined as a function of temperature (NTEMP
value on the TB command), with individual temperatures specified for the table entries (TBTEMP).
Example 3.20: Cast Iron
/prep7
mp, ex, 1,14.773E6
mp,nuxy, 1,0.2273
! Define cast iron model
TB,CAST,1,,,ISOTROPIC
TBDATA,1,0.04
TB,UNIAXIAL,1,1,5,TENSION
TBTEMP,10
TBPT,,0.550E-03,0.813E+04
TBPT,,0.100E-02,0.131E+05
TBPT,,0.250E-02,0.241E+05
TBPT,,0.350E-02,0.288E+05
TBPT,,0.450E-02,0.322E+05
TB,UNIAXIAL,1,1,5,COMPRESSION
TBTEMP,10
TBPT,,0.203E-02,0.300E+05
TBPT,,0.500E-02,0.500E+05
TBPT,,0.800E-02,0.581E+05
TBPT,,0.110E-01,0.656E+05
TBPT,,0.140E-01,0.700E+05
64
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
is very small, the Peirce model has less difficulty converging as compared to the Perzyna
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:
Constant
Meaning
C1
C2
where
. Here,
The EVH option can be combined with nonlinear (Chaboche) kinematic hardening using von Mises or
Hill yield criterion.
The six material constants in the EVH option are input (TBDATA) in the order shown:
Constant
Meaning
C1
C2
C3
C4
C5
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
65
Material Models
Constant
Meaning
C6
s0
Units
Stress
Q = Activation energy
Energy / Volume
Energy / Volume
temperature
C2
Q/R
C3
Pre-exponential factor
1 / Time
C4
xi
Stress multiplier
Dimensionless
C5
Dimensionless
C6
h0
Stress
Stress
C7
C8
Dimensionless
C9
Dimensionless
66
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
3.5.5. Creep
The creep strain rate, , 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. 68) and Explicit Creep
Equations (p. 69) 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 Material Model Element Support (p. 5).
Three types of creep equations are available:
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. 193)
for further information. Because 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.)
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-simultaneous 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 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
67
Material Models
Equation
Type
(TBOPT)
1
Strain
Hardening
C1>0
Primary
Time
Hardening
C1>0
Primary
Generalized
Exponential
Generalized
Graham
C1>0
Generalized
Blackburn
C1>0,
C3>0, Primary
C6>0
Modified Time
Hardening
C1>0
Primary
Modified Strain
Hardening
C1>0
Primary
Generalized
Garofalo
C1>0
Secondary
Exponential
form
C1>0
Secondary
10
Norton
C1>0
Secondary
11
Combined
Time
Hardening
Primary
C1>0,
+
C5>0
Secondary
12
Rational
polynomial
C2>0
13
Generalized
Time
Hardening
100
---
C1>0,
Primary
C5>0
Primary
+
Secondary
Primary
User Creep
---
where:
68
Primary
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
---
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 Material Model Element Support (p. 5).
See Creep in the Structural Analysis Guide for more information on this material option.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
69
Material Models
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.
Specify irradiation induced creep with constant C66. Irradiation Induced Explicit Creep Equation for C66
= 5 (p. 76) 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 indicated 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
= ( )(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.0e6, 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 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
TBDATA,1,c1,c2,c3,c4,,1
The explicit creep constants that you enter with the TBDATA are:
Constant
Meaning
C1-CN
where:
= change in equivalent strain with respect to time
= equivalent stress
70
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
where:
where:
x = 0 for C2
x = G + H for C2 < C3
C2 = 6000 psi (default), C3 = 25000 psi (default)
s, r, , G, and H = functions of temperature and stress as described in the reference.
This double exponential equation is valid for Annealed 304 Stainless Steel over a temperature range
from 800 to 1100F. 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 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
71
Material Models
where:
c = limiting value of primary creep strain
p = primary creep time factor
= secondary (minimum) creep strain rate
This standard rational polynomial creep equation is valid for Annealed 304 SS over a temperature range
from 427C to 704C. The equation is described completely in the High Alloy Steels. The first term describes 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.
72
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
A, B, and
This modified rational polynomial equation is valid for Annealed 2 1/4 Cr -1 Mo Low Alloy steel over a
temperature range of 700 - 1100F. 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.
where:
c = limiting value of primary creep strain
p = primary creep time factor
= secondary (minimum) creep strain rate
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
73
Material Models
This standard rational polynomial creep equation is valid for Annealed 2 1/4 Cr - 1 Mo Low Alloy Steel
over a temperature range from 371C to 593C. 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 modification 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.
where:
C1 = Scaling constant
M, N, K = Function of temperature (determined by linear interpolation within table) as follows:
Constant
Meaning
C5
C49
C50
C48 + C5
C48 + C5 + 1
First M value
...
C48 + 2C5
C5th M value
C48 + 2C5
C5th M value
...
C48 + 2C5
C5th M value
C48 + 2C5 +
1
First N value
...
74
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Meaning
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.
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.
where:
c = cpt/(1+pt) +
ln c = -1.350 - 5620/T - 50.6 x 10-6 + 1.918 ln (/1000)
ln p = 31.0 - 67310/T + 330.6 x 10-6 - 1885.0 x 10-12 2
ln
This creep law is valid for Annealed 316 SS over a temperature range from 800F to 1300F. 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 temperature range: 800 - 1300F. Maximum stress allowed for ec calculation: 45,000 psi; minimum stress: 0.0
psi. If T + Toffset < 1160, no creep is computed.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
75
Material Models
where:
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.
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
where:
B = FG + C63
= equivalent stress
76
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelasticity
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 1300F. 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.
3.6. Hyperelasticity
Hyperelastic material behavior is supported by current-technology shell, plane, and solid elements. For
a list of elements that can be used with hyperelastic material models, see Material Model Element
Support (p. 5). You can specify options to describe the hyperelastic material behavior for these elements.
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.
The following hyperelastic material model topics are available:
3.6.1. Arruda-Boyce Hyperelasticity
3.6.2. Blatz-Ko Foam Hyperelasticity
3.6.3. Extended Tube Hyperelasticity
3.6.4. Gent Hyperelasticity
3.6.5. Mooney-Rivlin Hyperelasticity
3.6.6. Neo-Hookean Hyperelasticity
3.6.7. Ogden Hyperelasticity
3.6.8. Ogden Compressible Foam Hyperelasticity
3.6.9. Polynomial Form Hyperelasticity
3.6.10. Response Function Hyperelasticity
3.6.11.Yeoh Hyperelasticity
3.6.12. Special Hyperelasticity
For information about other hyperelastic material models, see Special Hyperelasticity (p. 86).
where:
W = strain energy per unit reference volume
J = determinant of the elastic deformation gradient F
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
77
Material Models
= initial shear modulus of materials
L = limiting network stretch
d = material incompressibility parameter
The initial bulk modulus is defined as:
where:
W = strain energy per unit reference volume
= initial strain shear modulus
I2 and I3= second and third strain invariants
The initial bulk modulus k is defined as:
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.
78
TBOPT
Constants
C1
Gc
Purpose
Crosslinked network
modulus
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelasticity
TBOPT
Constants
Purpose
C2
Ge
Constraint network
modulus
C3
Empirical parameter (0
1)
C4
Extensibility parameter
C5
Incompressibility parameter
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.
where:
W = strain energy per unit reference volume
= initial shear modulus of material
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
79
Material Models
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:
where:
W = strain energy potential
c10, c01 = material constants characterizing the deviatoric deformation of the material
d = material incompressibility parameter
The initial shear modulus is defined as:
where:
d = (1 - 2*) / (C10 + C01)
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:
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:
80
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelasticity
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:
W = strain energy per unit reference volume
= initial shear modulus of the material
d = material incompressibility parameter.
J = determinant of the elastic deformation gradient F
The initial bulk modulus is defined by:
where:
W = strain energy potential
p = principal stretches of the left Cauchy-Green tensor
J = determinant of the elastic deformation gradient
N, p, p and dp = material constants
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
81
Material Models
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 material constants. For this reason, very high values of N are not recommended.
The initial shear modulus is defined by:
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.
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:
1, 1, 2, 2, d1, d2
For N (NPTS) = 3:
1, 1, 2, 2, 3, 3, d1, d2, d3
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.
where:
W = strain energy potential
82
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelasticity
J = determinant of the elastic deformation gradient
N, i, i and k = material constants
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 material constants. For this reason, very high values of N are not recommended.
The initial shear modulus is defined by:
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:
1, 1, 1
For N (NPTS) = 2:
1, 1, 2, 2, 1, 2
For N (NPTS) = 3:
1, 1, 2, 2, 3, 3, 1, 2, 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.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
83
Material Models
where:
W = strain energy potential
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:
c10, c01, c20, c11, c02, d1, d2
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.
84
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelasticity
where N is the NPTS value (TB,HYPER,,,,RESPONSE) and dk represents the material constants incompressibility 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.
where:
W = strain energy potential
J = determinant of the elastic deformation gradient F
N, ci0, and dk = 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 material 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 initial shear modulus is defined by:
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
85
Material Models
and the initial bulk modulus K is defined as:
86
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Hyperelasticity
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.
TBOPT
Constants
Purpose
Input Format
POLY
C1 to C31
Anisotropic strain
energy potential
TB,AHYPER,,,POLY
TBDATA,,A1,A2,A3,B1....
EXP
C1 to C10
Exponential
anisotropic strain
energy potential
TB,AHYPER,,,EXPO
TBDATA,1,A1,A2,A3,B1,B2,B3
TBDATA,7,C1,C2,E1,E2
AVEC
C1 to C3
Material direction
constants
TB,AHYPER,,,AVEC
TBDATA,,A1,A2,A3
BVEC
C1 to C3
Material direction
constants
TB,AHYPER,,,BVEC
TBDATA,,B1,B2, B3
PVOL
C1
Volumetric potential
TB,AHYPER,,,PVOL
TBDATA,,D
You can enter temperature-dependent data for anisotropic hyperelastic material via the TBTEMP command. For the first temperature curve, issue TB, AHYPER,,,TBOPT, then input the first temperature (TBTEMP). The subsequent TBDATA command inputs the data.
The program interpolates the temperature data to the material points automatically using linear interpolation. 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.
87
Material Models
shear distortional behavior defined by an isochoric strain energy potential, while the response to volumetric deformations is still purely elastic and characterized by a volumetric strain energy potential.
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
Constant
Meaning
C1
C2
N0
C3
Pa
C4
N1
( Block )2
Dimensionless
Material constant
s-1(Pa)-m
C5
Property
Units
Pa
C6
Material constant
Dimensionless
C7
Material constant
Dimensionless
C8
Dimensionless
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
Constant
Meaning
Property
C1
Units
1 / Pa
88
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelasticity
TB,CDM,MAT,NTEMPS,NPTS,TBOPT
The material constants for each valid TBOPT value follow:
Modified Ogden-Roxburgh Pseudo-Elastic
TBOPT = PSE2
Constant
Meaning
Property
C1
Damage variable
parameter
C2
Damage variable
parameter
C3
Damage variable
parameter
3.7. Viscoelasticity
Viscoelastic materials are characterized by a combination of elastic behavior, which stores energy during
deformation, and viscous behavior, which dissipates energy during deformation.
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 largedeformation formulations. Also presented is time-temperature superposition for thermorheologically
simple materials and a harmonic domain viscoelastic model.
3.7.1. Viscoelastic Formulation
3.7.2.Time-Temperature Superposition
3.7.3. Harmonic Viscoelasticity
For additional information, see Viscoelasticity in the Structural Analysis Guide.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
89
Material Models
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:
(3.13)
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:
(3.14)
90
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelasticity
(3.15)
where:
G0, K0 = relaxation moduli at t = 0
nG, nK = number of Prony terms
iG, iK = relative moduli
iG, iK = relaxation time
For use in the incremental finite element procedure, the solution for Equation 3.13 (p. 90) at t1 = t0 +
t is:
(3.16)
(3.17)
where si and pi are the deviatoric and pressure components, respectively, of the Cauchy stress for each
Maxwell element.
By default, the midpoint rule is used to approximate the integrals:
(3.18)
(3.19)
An alternative stress integration method is to assume a constant strain rate over the time increment.
Then the stress update is:
(3.20)
(3.21)
The model requires input of the parameters in Equation 3.14 (p. 90) and Equation 3.15 (p. 91). The relaxation 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:
Table Location
Constant
1G
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
91
Material Models
Table Location
Constant
1G
...
...
2(NPTS) - 1
nGG
2(NPTS)
nGG
For the bulk Prony data table, TBOPT = BULK, NPTS = nK, and the constants in the data table follow
this pattern:
Table Location
Constant
1K
1K
...
...
2(NPTS) - 1
nKK
2(NPTS)
nKK
Use TBOPT = INTEGRATION with the Prony table (TB,PRONY) to select the stress update algorithm. If
the table is not defined, or the value of the first table location is equal to 1, then the default midpoint
formula from Equation 3.18 (p. 91) and Equation 3.19 (p. 91) are used. If the value of the first table
location is equal to 2, then the constant strain rate formula from Equation 3.20 (p. 91) and Equation 3.21 (p. 91) are used.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelasticity
(3.23)
where:
right Cauchy-Green deformation tensor
isochoric part of C
determinant of the deformation gradient
The second Piola-Kirchhoff stress in the Maxwell element springs is then:
(3.24)
and the large-deformation stress update for the Maxwell element stresses is given by:
(3.25)
(3.26)
where:
deviatoric component of Si
pressure component of Si
An anisotropic hyperelastic model can also be used for Equation 3.23 (p. 93) , 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. 90)). The hyperelastic parameters for this model are input via a hyperelastic data
table (TB,HYPER). For more information, see Hyperelasticity (p. 77).
3.7.1.4. Dissipation
For a physical interpretation of the Prony series formulation, the dissipated energy in the viscoelastic
material is the energy used to deform the dashpots in the Maxwell elements. The increment of energy
used by the dashpots over a time increment is:
is the pressure,
For an additive decomposition of the Maxwell element strains, the dashpot strain increment is given
by:
where
is the strain increment, is the Maxwell spring compliance, and
is the stress increment.
Defining
as the instantaneous compliance, the Maxwell spring compliance is:
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
93
Material Models
The dashpot strain increment and compliance can be applied to both deviatoric and volumetric components, and then
corresponds to either the deviatoric or volumetric relative modulus for Maxwell
element i.
For large deformation, the incremental dissipation uses strain increments calculated from the rate of
deformation tensor and a Maxwell compliance tensor calculated from the instantaneous hyperelastic
stiffness and the relative moduli.
94
Table Location
Constant
Tr
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelasticity
Table Location
Constant
C1
C2
Constant
Tr
2
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 temperature that is different from the ambient temperature of the material. The shift function is given by:
(3.30)
where:
X = weight parameter
TF = absolute fictive temperature.
The partial fictive temperatures are calculated in the relative temperature scale defined by the input
material parameters. The evolving fictive temperature is given by:
(3.31)
where:
nf = number of partial fictive temperatures
Cfi = fictive temperature relaxation coefficient
Tfi = partial fictive temperature
The evolution of the partial fictive temperature is given by:
(3.32)
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
95
Material Models
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:
(3.33)
where:
T = temperature increment over the time step. The temperature increment in the first increment
is the body temperature at the end of the increment minus the fictive thermal strain reference
temperature, Tref, defined in the shift function table. If Tref is 0 or undefined in the shift function
table, the shift function reference temperature, Tr, is used to calculate the temperature increment
in the first time step.
g and l = glass and liquid coefficients, respectively, of thermal expansion given by:
(3.34)
(3.35)
Equation 3.34 (p. 96) and Equation 3.35 (p. 96) are evaluated at the relative current and fictive temperatures.
The parameters are input in a shift function data table (TB,SHIFT).
For the Tool-Narayanaswamy with fictive temperature shift function, TBOPT = FICT, NPTS = nf, and the
required input constants are:
Table Location
Constant
Tr
H/R
4 to 3(NPTS + 1)
3(NPTS + 1) + 1 to 3(NPTS + 1) + 5
3(NPTS + 1) + 6 to 3(NPTS + 1) +
10
3(NPTS + 1) + 11
Tref
96
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelasticity
(3.38)
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.
97
Material Models
Position
Value
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
2
3
4
Input complex tensile modulus on an experimental data table (TB,EXPE) with TBOPT = EMODULUS.
The data points are defined by:
Position
Value
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
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 experimental 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.
98
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Microplane
3.7.3.4. Stress
The magnitude of the real and imaginary stress components are obtained from expanding Equation 3.36 (p. 97) and using the storage and loss moduli from either the Prony series parameters or the
experimental data:
(3.39)
(3.40)
where:
Re() = real stress magnitude
Im() = imaginary stress magnitude
3.8. Microplane
The microplane model (TB,MPLANE) is based on research by Bazant and Gambarova [1][2] in which the
material behavior is modeled through uniaxial stress-strain laws on various planes.
Directional-dependent stiffness degradation is modeled through uniaxial damage laws on individual
potential failure planes, leading to a macroscopic anisotropic damage formulation.
The model is well suited for simulating engineering materials consisting of various aggregate compositions with differing properties (for example, concrete modeling, in which rock and sand are embedded
in a weak matrix of cements).
The microplane model cannot be combined with other material models.
The following topics concerning the microplane material model are available:
3.8.1. Microplane Modeling
3.8.2. Material Models with Degradation and Damage
3.8.3. Material Parameters Definition and Example Input
3.8.4. Learning More About Microplane Material Modeling
Also see Material Model Element Support (p. 5) for microplane.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
99
Material Models
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], expressed as:
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.
The strain split is expressed as:
where V is the second-order volumetric projection tensor and 1 the second-order identity tensor.
The deviatoric microplane strain vector D is calculated as:
where is the fourth-order identity tensor and the vector n describes the normal on the microsphere
(microplane).
The macroscopic strain is expressed as:
where v and D are the scalar volumetric stress and the deviatoric stress tensor on the microsphere,
and
.
Assume isotropic elasticity:
100
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Microplane
and
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 (p. 100) equation and the derived stresses (p. 100) equation are
solved via numerical integration:
3.8.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]
The following figure illustrates the discretization process:
Figure 3.18: Sphere Discretization by 42 Microplanes
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
101
Material Models
The damage parameter dmic is the normalized damage variable
where
and
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 (p. 102) 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.
The damage evolution is modeled by the following function:
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:
102
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Microplane
Figure 3.19: Damage Parameter d Depending on the Equivalent Strain Energy
, and mic.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
103
Material Models
E is Youngs modulus and is Poissons ratio. Both are microplane elastic properties and are defined
via the MP command.
The parameters k0, k1, k2,
Meaning
Property
C1
k0
C2
k1
C3
k2
C4
Critical equivalent-strain-energy
density
C5
mic
C6
mic
104
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Meaning
Property
Units
C1
kx
Permeability coefficient
Length4/Force/Time
C2
ky
Permeability coefficient
Length4/Force/Time
C3
kz
Permeability coefficient
Length4/Force/Time
Meaning
Property
Units
C1
Biot coefficient
Dimensionless
C2
km
Biot modulus
Force/Length2
3.10.1. Piezoelectricity
Piezoelectric capability (TB,PIEZ) is available with the coupled-field elements. (See Material Model Element
Support (p. 5) for piezoelectricity.) Material properties required for the piezoelectric effects include
the dielectric (relative permittivity) constants, the elastic coefficient matrix, and the piezoelectric matrix.
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 Anisotropy.
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:
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
105
Material Models
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 materials 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.
You can define up to 18 constants (C1-C18) with TBDATA commands (6 per command):
Constant
Meaning
C1-C6
C7-C12
C13-C18
See Piezoelectric Analysis in the Coupled-Field Analysis Guide for more information on this material
model.
3.10.2. Piezoresistivity
Elements with piezoresistive capabilities use the TB,PZRS command to calculate the change in electric
resistivity produced by elastic stress or strain. Material properties required to model piezoresistive materials are electrical resistivity, the elastic coefficient matrix, and the piezoresistive matrix.
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 information.
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:
Constant
Meaning
C1-C6
C7-C12
C13-C18
106
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Meaning
C19-C24
C25-C30
C31-C36
For 2-D problems, a 4x4 matrix relates terms ordered x, y, z, xy via 16 constants.
Constant
Meaning
C1-C4
C7-C10
C13-C16
C19-C22
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.
3.10.3. Magnetism
Elements with magnetic capability use the TB table to input points characterizing B-H curves. Temperature-dependent curves cannot be input.
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:
Constant
Meaning
Property
Magnetomotive
force/length
Flux/Area
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 establish 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 Material Properties and Real Constants in the Low-Frequency Electromagnetic Analysis Guide
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
107
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:
Constant
Meaning
C1-C6
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
3.11. Gasket
The gasket model (TB,GASKET) enables simulating the gasket joints with the interface elements. The
gasket material is usually under compression and is highly nonlinear. The material also exhibits quite
complicated unloading behavior when compression is released.
You can define some general parameters including the initial gap, stable stiffness for numerical stabilization, 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 unloading. 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 accurately model the gasket unloading behavior. When no unloading curves are defined, the material behavior 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:
Gasket
Data Type
General
parameters
TBOPT
Constants
Meaning
C1
C2
PARA
C3
108
Input Format
TB,GASKET,,,,PARA
TBDATA,1,C1,C2,C3
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Gasket
Gasket
Data Type
Compression
load closure
curve
TBOPT
COMP
Linear
unloading
data
LUNL
Nonlinear
unloading
data [2]
NUNL
Transverse
shear
TSS
Constants
Meaning
Input Format
Xi
Closure value.
Yi
Pressure value.
TB,GASKET,,,2,COMP
TBPT,,X1,Y1
TBPT,,X2,Y2
Xi
Closure value on
compression curve where
unloading started.
TB,GASKET,,,2,LUNL
TBPT,,X1,Y1
TBPT,,X2,Y2
Yi
Xi
Closure value.
Yi
Pressure value.
XY, XZ
TB,GASKET,,,2,NUNL
TBPT,,X1,Y1
TBPT,,X2,Y2
TB,GASKET,,,2,TSS
TBDATA,1,TSSXY,TSSXZ
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 temperature
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 interpolation. When the temperature is out of the specified range, the closest temperature point is used.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
109
Material Models
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.
3.12. Swelling
Swelling (TB,SWELL) is a material enlargement (volume expansion) caused by neutron bombardment
or other effects (such as moisture). The swelling strain rate is generally nonlinear and is a function of
factors such as temperature, time, neutron flux level, stress, and moisture content.
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.
Swelling strain is modeled using additive decomposition of strains, expressed as:
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
cannot 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:
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.
The following options for modeling swelling are available:
Linear swelling defines swelling strain rate as a function of fluence rate, expressed as:
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 User-Programmable Features.
Swelling equations are material-specific and are empirical in nature.
110
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Constant
Description
LINE
C1
Linear swelling
TBDATA,1,C1
EXPT
Exponential swelling
TBDATA,1,C1,C2,C3,C4
USER
C1, ..., Cn
User-defined
TBDATA,1,C1,C2,
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 Material Model Element Support (p. 5)
For more information about this material model, see Swelling in the Structural Analysis Guide.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
111
Material Models
Figure 3.21: Pseudoelasticity (PE) and Shape Memory Effect (SME)
112
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
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:
where is the material parameter, is the stress, and 1 is the identity tensor.
The evolution of the martensite fraction, S, is then defined as follows:
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
113
Material Models
where:
where
AS
AS
s
SA
s
SA
where
are the material parameters shown in Figure 3.23: Idealized Stress-Strain Diagram
of Superelastic Behavior (p. 114).
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
(
, respectively) as:
114
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Meaning
Property
C1
C2
C3
C4
C5
C6
115
Material Models
3.13.2.1.The Constitutive Model for Shape Memory Effect
3.13.2.2. Material Parameters for the Shape Memory Effect Option
where:
D = material elastic stiffness tensor
= total strain
= total transformation strain
= deviatoric transformation strain
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 austenite is observed in a stress-free state
h = material parameter related to the hardening of the material during the phase transformation
= indicator function introduced to satisfy the constraint on the transformation norm [1] in
which
116
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
where is defined by
where
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:
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.
J2 and J3 are defined as follows:
where is an internal variable and is called as transformation strain multiplier. and F(Xtr) must satisfy
the classical Kuhn-Tucker conditions, as follows:
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
117
Material Models
where DA is the elastic stiffness tensor of austenite phase, and DS is the elastic stiffness tensor of
martensite phase. The Poissons 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.
The following figure illustrates a number of the mechanical model features:
Figure 3.24: 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
. 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. Because
the material tangent stiffness matrix is generally unsymmetric, use the unsymmetric Newton-Raphson
option (NROPT,UNSYM) to avoid convergence problems.
118
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Hardening parameter
C2
To
Reference temperature
C3
Elastic limit
C4
C5
C6
Em
Martensite modulus
C7
Example 3.23: 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
, is available
Elastic strain energy density is available as part of the strain energy density record SEND (ELASTIC).
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
119
Material Models
120
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
MPC184 Joint
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:
Joint Element
Constant
Meaning
C16
Term D44
C21
Term D66
Universal joint
Slot joint
C1
Term D11
Point-in-plane joint
Translational joint
C1
Term D11
General joint
---
Screw joint
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.
121
Material Models
The appropriate TBOPT labels for each joint element type are shown in the following tables. For a description of each TBOPT label, see Joint Element Specifications (JOINT) in the TB command documentation.
Nonlinear Stiffness Behavior
Joint Element
TBOPT on TB command
JNSA, JNS4
JNSA, JNS6
Universal joint
Slot joint
Point-in-plane joint
Translational joint
General joint
Screw joint
Joint Element
TBOPT on TB command
JNDA, JND4
JNDA, JND6
Universal joint
Slot joint
Point-in-plane joint
Translational joint
General joint
Screw joint
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
122
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
MPC184 Joint
TBPT,,rotation_value_1,moment_value_1
TBPT,,rotation_value_2,moment_value_2
Consider a function where the damping force varies with temperature and relative velocity:
For more information about the Function Tool utility, see Using the Function Tool in the Basic Analysis
Guide.
123
Material Models
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 consideration.
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 necessary.
The default value for the translational joint and the slot joint is set to 0.005*h, where h is a characteristic
length value computed from overall dimensions of the model. The value of h defaults to 1.0 if a characteristic length cannot be computed properly. The default value for the revolute joint is set to 0.001 radians.
The frictional behavior is implemented using a penalty method. Thus, there will be relative elastic slip even
when sticking conditions prevail. The amount of elastic slip depends on the value specified for elastic slip.
In some cases, the default values may result in large elastic slip. Therefore, you should specify an amount
of elastic slip that is appropriate for your model.
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 stickstiffness 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).
Interference Fit Force/Moment
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:
TBOPT Labels for Elements Supporting Coulomb Friction
Friction
Parameter
x-axis Revolute
Joint
z-axis Revolute
Joint
Slot Joint
Translational
Joint
Static Friction
MUS4
MUS6
MUS1
MUS1
Exponential
Friction Law
EXP4
EXP6
EXP1
EXP1
Max. Allowable
Shear
Force/Moment
TMX4
TMX6
TMX1
TMX1
124
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Contact Friction
Elastic Slip
SL4
SL6
SL1
SL1
Interference Fit
Force/Moment
FI4
FI6
FI1
FI1
Stick-Stiffness
SK4
SK6
SK1
SK1
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 xaxis revolute joint.
TB, JOIN, 1, , , MUS4
TBDATA, 1, 0.1
TB, JOIN, 1, , , SK4
TBDATA, 1, 3.0E4
TB, JOIN, 1, , , FI4
TBDATA, 1, 10000.00
!
!
!
!
!
!
Label
Value
Label
Value
Label
Value
Example 2 Specifying temperature dependent friction coefficient and other friction parameters for
a z-axis revolution joint.
TB, JOIN, 1,2 , 1, MUS6
TBTEMP, 10
TBDATA, 1, 0.15
TBTEMP, 20
TBDATA, 1, 0.1
!
TB, JOIN, 1, , , SK4
TBDATA, 1, 3.0E4
TB, JOIN, 1, , , FI4
TBDATA, 1, 10000.00
Example 3
joint.
!
!
!
!
Label
Value
Label
Value
for
for
for
for
stick-stiffness
stick-stiffness
interference fit force
interference fit force
Specifying the exponential law for friction and other friction parameters for a z-axis revolute
Example 4
!
!
!
!
!
!
!
!
!
!
!
!
!
Label
Value
Label
Value
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
125
Material Models
See Understanding Field Variables (p. 199) 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
TBTEMP,100.0
TBDATA,1,MU
TBTEMP,200.0
TBDATA,1,MU
!
!
!
!
!
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.
Note that if the coefficient of friction is defined as a function of temperature, the program always uses
the contact surface temperature as the primary variable (not the average temperature from the contact
and target surfaces).
126
Contact Friction
TBFIELD,SLRV,20.0
TBFIELD,NPRE,150.0
TBDATA,1,MU1,MU2
TBFIELD,NPRE,300.0
TBDATA,1,MU1,MU2
!
!
!
!
!
Define
Define
Define
Define
Define
See Understanding Field Variables (p. 199) 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
TBTEMP,100.0
TBDATA,1,MU1,MU2
TBTEMP,200.0
TBDATA,1,MU1,MU2
!
!
!
!
!
Note that if the coefficient of friction is defined as a function of temperature, the program always uses
the contact surface temperature as the primary variable (not the average temperature from the contact
and target surfaces).
!* LOAD STEP 2
...
TIME,2
SOLVE
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
127
Material Models
Field variables specified with the TBFIELD command are not available for TB,FRIC,,,,USER.
For detailed information on using the USERFRIC subroutine, see Writing Your Own Friction Law
(USERFRIC) in the Contact Technology Guide.
Interaction
Behavior
Description
STANDARD Standard
unilateral contact
ROUGH
NOSEPE
No separation
No separation contact in which the target and contact surfaces are tied
(sliding permitted) for the remainder of the analysis once contact is established (although
sliding is permitted). (Similar to KEYOPT(12) = 2 for pair-based contact
elements.)
BONDED
Bonded contact
(no separation, no
sliding)
Bonded contact in which the target and contact surfaces are bonded in
all directions (once contact is established) for the remainder of the
analysis. (Similar to KEYOPT(12) = 3 for pair-based contact elements.)
ANOSEP
No separation
(always)
128
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Contact Interaction
TBOPT
Label
Interaction
Behavior
Description
ABOND
Bonded (always)
Bonded contact in which contact detection points that are either initially
inside the pinball region or that once involve contact always attach to
the target surface along the normal and tangent directions to the contact
surface (fully bonded). (Similar to KEYOPT(12) = 5 for pair-based contact
elements.)
IBOND
Bonded (initial
contact)
Bonded contact in which the contact detection points that are initially
in a closed state will remain attached to the target surface, and the
contact detection points that are initially in an open state will remain
open throughout the analysis. (Similar to KEYOPT(12) = 6 for pair-based
contact elements.)
All of the above options use one material constant (C1) on the TBDATA command. The value of C1
defines the effect of initial penetration or gap, as described in the table below.
C1
Value
Description
[1]
Include both initial geometrical penetration or gap and offset. (Similar to KEYOPT(9) = 0 for
pair-based contact elements.)
[2]
Exclude both initial geometrical penetration or gap and offset. (Similar to KEYOPT(9) = 1 for
pair-based contact elements.)
Include both initial geometrical penetration or gap and offset, but with ramped effects. (Similar
to KEYOPT(9) = 2 for pair-based contact elements.)
Include offset only (exclude initial geometrical penetration or gap). (Similar to KEYOPT(9) = 3
for pair-based contact elements.)
Include offset only (exclude initial geometrical penetration or gap), but with ramped effects.
(Similar to KEYOPT(9) = 4 for pair-based contact elements.)
Include offset only (exclude initial geometrical penetration or gap) regardless of the initial
contact status (near-field or closed). (Similar to KEYOPT(9) = 5 for pair-based contact elements.)
Include offset only (exclude initial geometrical penetration or gap), but with ramped effects
regardless of the initial contact status (near-field or closed). (Similar to KEYOPT(9) = 6 for
pair-based contact elements.)
1. C1 = 0 is the default for all general contact interaction types if TBDATA is issued but C1 is not specified.
2. C1 = 1 is the default for all general contact interaction types if TBDATA is not issued.
Note that the effects of C1 input are dependent on which TBOPT label was previously defined. The indicated initial gap effect is considered only if TBOPT = ANOSE or TBOPT = ABOND is defined.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
129
Material Models
After writing a USERINTER subroutine to program your interaction model, you incorporate the model
in your analysis by using the command TB,INTER with TBOPT = USER and specifying the interaction
properties on the TBDATA command as shown below.
TB,INTER,1,,2,USER
TBDATA,1,PROP1,PROP2
For detailed information on using the USERINTER subroutine, see Defining Your Own Contact Interaction
(USERINTER) in the Contact Technology Guide.
3.17.1. Exponential Cohesive Zone Material for Interface Elements and Contact
Elements
Interface elements and contact elements allow exponential cohesive zone materials to be used for
simulating interface delamination and other fracture phenomena. To define exponential material behavior, issue the TB,CZM,,,,EXPO command, then specify the following material constants via the TBDATA
command:
Constant
Meaning
Property
C1
max
C2
C3
To define a temperature dependent material, use the TBTEMP command as shown below:
TB,CZM,1,2,,EXPO
TBTEMP,100.0
TBDATA,1,max,n,t
TBTEMP,200.0
TBDATA,1,max,n,t
!
!
!
!
!
3.17.2. Bilinear Cohesive Zone Material for Interface Elements and Contact
Elements
Interface elements and contact elements allow bilinear cohesive zone materials to be used for simulating
interface delamination and other fracture phenomena. To define bilinear material behavior, issue the
TB,CZM,,,,BILI command, then specify the following material constants via the TBDATA command:
130
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Meaning
Property
C1
max
C2
C3
C4
C5
C6 [1]
Ratio of
to
, or ratio of
to
,-max,
Mode I Dominated
Mode II Dominated
Mixed-Mode
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
131
Material Models
Constant
Meaning
Property
C1
Damping coefficient
To define a temperature-dependent material, use the TBTEMP command as shown in the following
example input fragment:
! define first temperature
TBTEMP,100.0
!define damping coefficient at temp 100.0
TBDATA,1,c1
!define second temperature
TBTEMP,200.0
!define damping coefficient at temp 200.0
TBDATA,1,c1
For more information, see Viscous Regularization in the Mechanical APDL Theory Reference.
Meaning
Property
C1
max
C2
C3
C4
C5
C6
1. For contact elements using the force-based model (see the description of KEYOPT(3) for
CONTA175, CONTA176, and CONTA177), input a contact force value for this quantity.
To define a temperature dependent material, use the TBTEMP command as shown below:
TB,CZM,1,2,,CBDD
TBTEMP,100.0
TBDATA,1,max,
TBTEMP,200.0
,max,
,, ! Define material constants at temp 100.0
! Define second temperature
TBDATA,1,max,
,max,
,,
132
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Meaning
Property
C1
max
C2
Gcn
C3
max
C4
Gct
C5
C6
1. For contact elements using the force-based model (see the description of KEYOPT(3) for
CONTA175, CONTA176, and CONTA177), input a contact force value for this quantity.
2. For contact elements using the force-based model (see the description of KEYOPT(3) for
CONTA175, CONTA176, and CONTA177), this quantity is critical fracture energy.
To define a temperature dependent material, use the TBTEMP command as shown in the following
example input fragment:
TB,CZM,1,2,,CBDE
TBTEMP,100.0
TBDATA,1,max,Gcn,max,Gct,,
TBTEMP,200.0
TBDATA,1,max,Gcn,max,Gct,,
!
!
!
!
!
!
Debonding Modes
Debonding involves separation of surfaces forming an interface. The direction of separation determines
the debonding mode. The program detects the debonding mode based on material data that you input
for normal and tangential directions:
Mode I debonding involves separation normal to the interface. It is activated by inputting data items
C1, C2, and C5 on the TBDATA command.
Mode II debonding involves slip tangent to the interface. It is activated by inputting data items C3,
C4, and C5 on the TBDATA command.
Mixed mode debonding involves both normal separation and tangential slip. It is activated by inputting
data items C1, C2, C3, C4, C5, and C6 on the TBDATA command.
133
Material Models
the material properties for any number of healing cycles, but be sure to start with a cycle number of
zero.
For example, the following commands specify healing of the CZM interface if the contact surfaces come
into contact after they are completely debonded:
TB,CZM,1,,,CBDE
TBFIELD,CYCLE,0
TBDATA,1,max,Gcn,max,Gct,,
TBFIELD,CYCLE,1
TBDATA,1, ...
!
!
!
!
!
When the contact interface is completely debonded and the surfaces come into contact again, the debonding parameter is set to 0 thus effectively healing the CZM. The healing cycle is incremented by
one and the appropriate material data is interpolated for this healing cycle.
This healing option is only available when one of the supported cohesive zone materials is used with
contact elements. It is not available when a cohesive zone material is used with interface elements.
Wear coefficient, K
C2
Material hardness, H
C3
Pressure exponent, m
C4
Velocity exponent, n
C5
Optional flag to control how the wear increment is calculated (see below)
134
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Acoustics
Set C5 to -99 to calculate wear for postprocessing purposes only. The Archard model is used to calculate
the wear, but the contact nodes are not moved; thus wear is just a postprocessing variable and does
not affect the solution.
Use the TBFIELD command to define the constants as a function of temperature and/or time. A sample
input with constants as a function of time is shown below:
TB,WEAR,1,,,ARCD
TBFIELD,TIME,0
TBDATA,1,K,H,m,n
TBFIELD,TIME,1
TBDATA,1,K,H,m,n
!
!
!
!
!
Alternatively, you can use the TBTEMP command to define the constants as a function of temperature
alone. A sample input is shown below:
TB,WEAR,1,,,ARCD
TBTEMP,100
TBDATA,1,K,H,m,n
TBTEMP,200
TBDATA,1,K,H,m,n
!
!
!
!
!
3.19. Acoustics
The following topics related to acoustic materials are available:
3.19.1. Equivalent Fluid Model of Perforated Media
3.19.2. Acoustic Frequency-Dependent Materials
3.19.3. Low Reduced Frequency (LRF) Model of Acoustic Viscous-Thermal Media
135
Material Models
TB,PERF,,,,JCA
The effective density is given by:
(3.41)
where:
0 = density of fluid
= fluid resistivity
= porosity
= tortuosity
= viscous characteristic length
= dynamic viscosity
The effective bulk modulus is given by:
(3.42)
where:
= specific heat ratio
P0 = static reference pressure
Prt = Prandtl number
' = thermal characteristic length
The constants C1 through C5 (entered via the TBDATA command) are:
Constant
Meaning
C1
C2
Porosity (defaults to 1)
C3
Tortuosity (defaults to 1)
C4
C5
Additional material parameters are input with the MP and R commands. For more information, see
Equivalent Fluid of Perforated Materials in the Mechanical APDL Theory Reference.
136
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Acoustics
The propagating constant is given by:
where:
= fluid resitivity
f = frequency
= angular frequency
The constant C1 (entered via the TBDATA command) is:
Constant
Meaning
C1
where:
= fluid resitivity
f = frequency
= angular frequency
The constant C1 (entered via the TBDATA command) is:
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
137
Material Models
Constant
Meaning
C1
where:
R = Resistance
X=
= attenuation constant
= phase constant
The constants C1 through C4 (entered via the TBDATA command) are:
Constant
Meaning
C1
C2
C3
C4
3.19.1.5. Complex Density and Velocity Equivalent Fluid Model of Perforated Media
To define a complex impedance-propagating constant equivalent fluid model of a porous medium in
an acoustic full harmonic analysis, issue this command:
TB,PERF,,,,CDV
The complex density is given by:
= r + ji
The complex sound speed is given by:
c = cr + jci
where:
138
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Acoustics
r = real part of complex density (kg/m3)
i = imaginary part of complex density (kg/m3)
cr = real part of complex sound speed (m/s)
ci = imaginary part of complex sound speed (m/s)
The constants C1 through C4 (entered via the TBDATA command) are:
Constant
Meaning
C1
C2
C3
C4
where:
n1 = normal velocity at port 1
1 = pressure at port 1
n2 = normal velocity at port 2
2 = pressure at port 2
Y11, Y12, Y13 = complex admittance elements
1 = internal source related to port 1 (usually zero in acoustic applications)
2 = internal source related to port 2 (usually zero in acoustic applications)
The constants C1 through C12 (entered via the TBDATA command) are:
Constant
Meaning
C1
C2
C3
C4
C5
C6
C7
C8
139
Material Models
Constant
Meaning
C9
C10
C11
C12
For an acoustic 2 x 2 transfer admittance matrix, the port number (SF,Nlist,PORT) can be any positive
integer.
If the two ports of the transfer admittance matrix are connecting to the fluid, the smaller port number
corresponds to port 1 of the 2 x 2 transfer admittance matrix and the greater port number corresponds
to port 2.
If one port of the transfer admittance matrix is connecting to the acoustic-structural interaction interface
and another port is connecting to the fluid, the FSI interface (SF,Nlist,FSI) corresponds to port 1 and
the defined port number (SF,Nlist,PORT) corresponds to port 2 of the transfer admittance matrix.
A pair of ports of the 2 x 2 transfer admittance matrix must be defined in the same element.
where:
Y = Complex admittance elements determined by geometric dimension and material
= Ratio of inner and outer radius for cylindrical structure (default = 1)
The constants C1 through C6 (entered via the TBDATA command) are:
Constant
Meaning
C1
C2
C3
C4
C5
140
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Acoustics
Constant
Meaning
C6
Meaning
C1
C2
C3
C4
C5
C6
C7
Meaning
C1
141
Material Models
Constant
Meaning
C1
C2
Meaning
C1
3.20. Fluids
Fluid material models can be used with hydrostatic fluid elements to model compressible fluids. For
theoretical background on these materials, see Fluid Material Models in the Mechanical APDL Theory
Reference. For more information on using these fluid material models with the hydrostatic fluid elements,
see Modeling Hydrostatic Fluids in the Structural Analysis Guide.
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:
Constant
Meaning
Property
C1
Bulk modulus
C2
C3
0f
Initial density
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
TBTEMP,100.0
TBDATA,1,K,,0f
TBTEMP,200.0
TBDATA,1,K,,0f
!
!
!
!
!
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
142
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Meaning
Property
C1
0f
Initial density
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
TBTEMP,100.0
TBDATA,1,0f
TBTEMP,200.0
TBDATA,1,0f
!
!
!
!
!
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, followed by TBPT commands for up to 100 pressure-volume data points. The data points (X, Y) entered
on TBPT are:
Constant
Meaning
Pressure value
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.
143
Material Models
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.
Define temperatures via TBTEMP commands.
Example 3.24: Input for a User-Defined Material
TB,USER,1,2,4
TBTEMP,1.0
TBDATA,1,19e5,0.3,1e3,100,
TBTEMP,2.0
TBDATA,1,21e5,0.3,2e3,100,
!
!
!
!
!
!
!
!
!
!
For information about state variable support, see Using State Variables with the UserMat Subroutine (p. 145).
TBTEMP,1.0
TBDATA,1,19e5,0.3,1e3,100,
TBTEMP,2.0
TBDATA,1,21e5,0.3,2e3,100,
!
!
!
!
!
!
!
!
!
!
For information about state variable support, see Using State Variables with the UserMatTh Subroutine (p. 145).
144
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
TBTEMP,1.0
TBDATA,1,19e5,0.3,1e3,100,
TBTEMP,2.0
TBDATA,1,21e5,0.3,2e3,100,
!
!
!
!
!
!
!
!
!
!
Contact elements do not support a user-defined cohesive material. However, a similar capability can
be achieved by defining your own contact interface behavior via the userinter subroutine. For more
information, see Defining Your Own Contact Interaction (USERINTER) in the Contact Technology Guide.
For information about state variable support, see Using State Variables with the UserCZM Subroutine (p. 146).
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
145
Material Models
146
TBOPT = 1
TBOPT = 2
C1
C2
C3
C4
C5
C6
C7
C8
C9
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
TBOPT = 1
TBOPT = 2
C10
--
C11
--
C12
--
C13
--
C14
XZIC -- XZ compressive
inclination parameter for Puck
failure index (default = 0.0)
--
C15
--
C16
YZIC -- YZ compressive
inclination parameter for Puck
failure index (default = 0.0)
--
C17
--
C18
--
C19
--
C20
--
C2
C3
C4
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
147
Material Models
Strength
Max
Max
Tsai-Wu
Puck
Hashin
LaRc03/04
Limit
Strain
Stress
Strength
User-Defined
Criterion Criterion Criterion
Constants Criterion Criterion
Ratio
C5
--
--
--
C6
--
--
--
C7
C8
--
--
C9
--
--
--
C10
--
--
--
--
--
C11
--
--
--
--
--
C12
--
--
--
--
--
C13
--
--
--
--
--
C14
--
--
--
--
--
C15
--
--
--
--
--
C16
--
--
--
--
--
C17
--
--
--
--
--
C18
--
--
--
--
--
C19
--
--
--
--
--
C20
--
--
--
--
--
148
TBOPT = 1 or FCRT
NPTS = 4
C1
C2
C3
C4
Material Damage
1. The available failure criteria are as follows:
1 -- Maximum strain
2 -- Maximum stress
3 -- Puck
4 -- Hashin
5 -- LaRc03
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:
TB,DGMI Input
TBOPT
value
Compatible TB,DMGE,,,,TBOPT
Option
1, 2, ..., 19
1 or MPDG
4 (Hashin failure
criteria only)
2 or CDM
TBOPT = 1 or MPDG
NPTS = 4
C1
C2
C3
C4
Constants
C1
TBOPT = 2 or CDM
NPTS = 4
Energy dissipated per unit area from tensile
fiber damage [1]
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
149
Material Models
Constants
TBOPT = 2 or CDM
NPTS = 4
C2
C3
C4
C5
C6
C7
C8
For complex stress state, the equivalent stresses and strains are calculated based on Hashin
failure criteria.
150
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Damage
2. Viscous damping coefficients are also specified respectively for all four damage modes.
For a specific damage mode, the damage evolution is regularized as follows:
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:
TB,DGME,,,,TBOPT option
1 or MPDG
1, 2, ..., 19
2 or CDM
4 (Hashin failure
criteria only)
For a general orthotropic material, the damaged elastic matrix is defined as:
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
151
Material Models
For a transversely isotropic material with plane stress state, primarily adopted for thin fiber-reinforced
composite structures, the damaged elastic matrix can be expressed as:
Four damage modes (fiber tension [rupture], fiber compression [kinking], matrix tension [cracking], and
matrix compression [crushing]) are accounted for. Four damage variables (one for each mode) are used
to measure damage. The damage variables for calculating the damaged elasticity matrix are determined
as follows:
152
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Damage
When a damage mode is initiated, the damage progresses immediately, indicated by the increasing
damage variable for the mode.
Two damage evolution methods are available:
Material property degradation method (TB,DMGE,,,,MPDG)
The material stiffness is instantly reduced based on the damage variables, explicitly specified via the
TB,DMGE command. Any physical failure criteria can be used to detect the onset of the damage.
Continuum damage mechanics method (TB,DMGE,,,,CDM)
Damage variables increase gradually based on the energy amounts dissipated for the various damage
modes. To achieve an objective response, the dissipated energy for each damage mode is regularized
as follows:
The characteristic length Le is calculated from the element area A via the following expressions:
With this characteristic length, the constitutive relation is converted from stress-strain to stress-displacement relation. The damage evolution function is derived from Hashin failure criteria; therefore,
only Hashin failure criteria are allowed via TB,DMGI for detecting damage onset. For a plane stress
state, the equivalent displacements
and equivalent stresses
are given below for all damage
modes.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
153
Material Models
In expressions above,
154
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
155
156
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
157
Column 1
Column 2
Uniaxial Test
Engineering Strain
Engineering Stress
Biaxial Test
Engineering Strain
Engineering Stress
Planar/Shear Test
Engineering Strain
(in loading
direction)
Engineering Stress
Engineering Shear
Strain
Engineering Shear
Stress
Volumetric Test
Hydrostatic Pressure
Column 3
(Optional) Engineering
Normal Stress (normal to
the edge of shear stress)
Column 1
Column 2
Column 3
Uniaxial Test
Engineering Strain
Engineering
Stress
158
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Column 1
Column 2
Column 3
Biaxial Test
Engineering Strain
Engineering
Stress
Shear Test
Engineering
Stress
Volumetric Test
(Optional)
Engineering
Normal Stress
(normal to the
edge of shear
stress)
Hydrostatic Pressure
60.00
118.2
175.2
231.1
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 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
159
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.).
Abbreviation
Order/Options
No. of
Coefficients [1]
Linear/Nonlinear
Fitting
Mooney-Rivlin
moon
2, 3, 5, 9
2/3/5/9+1
Linear
Polynomial
poly
1 to N
Linear
Yeoh
yeoh
1 to N
N+N
Linear
Neo-Hookean
neoh
1+1
Linear
Ogden
ogde
1 to N
2*N+N
Nonlinear
Arruda-Boyce
boyc
2+1
Nonlinear
Gent
gent
2+1
Nonlinear
Blatz-Ko
blat
Nonlinear
Ogden
Hyper-foam
foam
1 to N
2*N+N
Nonlinear
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.
Blatz-Ko and Ogden hyper-foam are compressible models. For Ogden hyper-foam, the experimental
data you supply will require additional fields.
160
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
where:
ID = Index corresponding to the material number
Option2 = model name, as specified in Table 5.3: Hyperelastic Curve-Fitting Model Types (p. 160)
(above)
Option3 = order or number of coefficients, where applicable. Table 5.3: Hyperelastic Curve-Fitting
Model Types (p. 160) (above) specifies the number and type of coefficient(s) necessary for each hyperelastic model type.
161
5.1.5.1. Batch
TBFT,SET,ID,HYPER,Option2,Option3,Option4,Option5
! initialize coefficients
where:
ID = Index corresponding to the material number
Option2 = model name (see Table 5.3: Hyperelastic Curve-Fitting Model Types (p. 160) (above) for
the available models)
Option3 = order, if applicable
Option4 = index of coefficient
Option5 = value of the coefficient
To set Option4 for a reference temperature, or for temperature dependency:
TBFT,SET,ID,HYPER,Option2,Option3,Option4,Option5
where:
ID = Index corresponding to the material number
Option2 = model name. See Table 5.3: Hyperelastic Curve-Fitting Model Types (p. 160) (above) for
the available models.
Option3 = order, if applicable
Option4 = tdep or tref
Option5 = If Option4 = tdep, Option5 = 1 turns on temperature dependency, and Option5
= 0 turns it off. If Option4 = tref, Option5 will be either all, or the reference temperature.
TBFT,FIX,ID,HYPER,Option2,Option3,Option4,Option5
where:
ID = Index corresponding to the material number
Option2 = model name See Table 5.3: Hyperelastic Curve-Fitting Model Types (p. 160) (above) for
the available models.
Option3 = (Blank--not applicable)
Option4 = index of coefficient
Option5 = 1 for fixed and 0 to vary
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.
162
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
5.1.6.1. Batch
TBFT,SOLVE,ID,HYPER,Option2,Option3,Option4, ..., Option7
where:
ID = Index corresponding to the material number
Option2 = model name. See Table 5.3: Hyperelastic Curve-Fitting Model Types (p. 160) (above) for
the models available.
Option3 = order or number of your coefficients. See Table 5.3: Hyperelastic Curve-Fitting Model
Types (p. 160) (above) for possible values.
Option4 = curve-fitting procedure: 0 = non-normalized least squares, 1 = normalized least squares
Option5 = maximum number of iterations
Option6 = tolerance of residual changes
Option7 = tolerance of coefficient changes
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.
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 experiments are plotted in the GUI window. To view specific data and its corresponding fitting result, you
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
163
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 experimental data, repeating step 3 through step 7 until you get an acceptable solution.
5.1.8.1. Batch
TBFT,FSET,ID,HYPER,Option2,Option3
! write data to TB
where:
ID = Index corresponding to the material number
Option2 = model name. See Table 5.3: Hyperelastic Curve-Fitting Model Types (p. 160) (above) for
the models available.
Option3 = order or number of your coefficients. See Table 5.3: Hyperelastic Curve-Fitting Model
Types (p. 160) (above) for possible values.
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.
164
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
165
time
Shear modulus
smod
Bulk modulus
bmod
Temperature
temp
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
166
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
5.2.3.1. Batch
TBFT,EADD,ID,Option1,Option2,Option3,Option4
! input data
where:
ID = index corresponding to the material number
Option1 = experimental data type, either sdec or bdec
Option2 = name of file containing experimental data
Option3 = file name extension
Option4 = file directory
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)
where:
ID = Index corresponding to the material number
Option2 = PVHE (refers to Prony Viscohypoelastic)
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
167
where:
ID = Index corresponding to the material number
CATEGORY = VISCO
Option2 = pshear or pbulk or shift
Option3 = dependent on Option2 as follows:
Option2 = pshear or bulk, Option3 = NONE, or 1 to N
Option2 = shift, Option3 = NONE or TN or WLF
TBFT,FCASE,ID,FINI
168
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
where:
ID = Index corresponding to the material number
Option2 = Case name
Option3 = (Blank--not applicable)
Option4 = Index of coefficient
Option5 = Value of coefficient
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
169
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:
ID = Index corresponding to the material number
Option2 = Case name
Option3 = (Blank--not applicable)
Option4 = Index of coefficient
Option5 = 1 to fix, 0 to vary
Example 5.2: Fixing Coefficient Values
TBFT,FIX,1, myvisco1,,1,1.! Fix the first coefficient to a value set via TBFT,SET
TBFT,FIX,1, myvisco1,,2,1 ! Fix the second coefficient to a value set via TBFT,SET
Temperature dependency uses Option4, and references data files you entered with the appropriate
temp header.
TBFT,SET,ID,VISCO,Option2,Option3,Option4,Option5
where:
ID = Index corresponding to the material number
Option2 = Model name. See Table 5.4: Viscoelastic Data Types and Abbreviations (p. 166) (above)
for the models available.
Option3 = (Blank--not applicable)
Option4 = tdep or tref
Option5 = If Option 4 = tdep, then 1 turns temperature dependency on and 0 turns it off. If
Option 4 = tref, this value will be a specific temperature, or ALL.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
171
where:
ID = Index corresponding to the material number
Option2 = Case name
Option3 = (Blank -- not applicable)
Option4 = comp
Option5 = pshea (for Shear only) or pbulk (for bulk only)
The SOLVE option allows you to specify procedure types, tolerances, and number of iterations.
TBFT,SOLVE,ID,CASE,Option2,Option3,Option4, ..., Option7
172
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
173
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:
Conversion Procedure:
where:
Option2 = Case
Option3 = Case name
174
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Comments
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
175
Comments
they are unacceptable, repeat steps 3 through 5 to
obtain a new curve-fitting solution.
time
creq
dcreq
Equivalent Stress
seqv
Temperature
temp
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. 176).
Following is a typical data input file:
/1,seqv ! indicates first column is stress
/2,creq ! indicates second column is creep strain
/3,temp ! indicates third column is temperature
/4,dcreq ! indicates fourth column is creep strain rate
4000
0.00215869 100
0.000203055
4000
0.00406109 100
0.000181314
4000
0.00664691 100
0.000165303
4000
0.0102068 100
0.000152217
4000
0.0151416 100
0.000140946
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
176
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Thirteen model types are available for creep curve fitting. The model you choose determines the experimental data required for the curve-fitting process. The following table describes the creep data required
to perform curve fitting for each model type. For strain hardening and modified strain hardening, you
must input both creep strain and creep strain rate in the experimental data.
Table 5.6: Creep Model and Data/Type Attribute
Creep Model
Strain Hardening
creq
dcreq
time
seqv
temp
Time Hardening
Generalized Exponential
Generalized Graham
Generalized Blackburn
Generalized Garofalo
Exponential Form
Norton
x
x
where:
ID = Index corresponding to the material number
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
177
=
=
=
=
creep
name of file containing experimental data
file name extension
file directory
where:
ID = Index corresponding to the material number
Option2 = creep model abbreviation. See Table 5.7: Creep Models and Abbreviations (p. 178).
Option3 = not used for creep curve fitting.
The following table describes the creep models available and their abbreviated names for Option2
(above).
Table 5.7: Creep Models and Abbreviations
Model
Number
Name
Fitting
Name/Option2
Strain Hardening
shar
Time Hardening
thar
Generalized Exponential
gexp
Generalized Graham
ggra
Generalized Blackburn
gbla
mtha
msha
Generalized Garofalo
ggar
Exponential Form
expo
10
Norton
nort
11
psth
178
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Name
Fitting
Name/Option2
12
psrp
13
gtha
The experimental data must be consistent with the creep model you choose. See Table 5.6: Creep
Model and Data/Type Attribute (p. 177) for the data types required for each creep model.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
179
where:
ID = Index corresponding to the material number
Option2 = creep model name
Option3 = order, if applicable
Option4 = index of coefficient
Option5 = value of coefficient
To modify the coefficients, use the FIX option of the TBFT command.
TBFT,FIX,ID,CREEP,Option2,Option3,Option4,Option5
where:
ID = Index corresponding to the material number
Option2 = creep model name
Option3 = order, if applicable
Option4 = index of coefficient
Option5 = 0 - variable, 1 - fixed
TBFT,SET,ID,CREEP,Option2,Option3,Option4,Option5
where:
ID = Index corresponding to the material number
Option2 = creep model name
Option3 = order, if applicable
Option4 = tdep or tref
Option5 = For Option4 = tdep, 1 to activate, 0 to deactivate. For Option4 = tdep, a specific
temp value, or all.
180
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
where:
ID = Index corresponding to the material number
Option2 = creep function name (See Table 5.7: Creep Models and Abbreviations (p. 178))
Option3 = ((Blank--not applicable for creep models)
Option4 = Error norm: 0 = non-normalized, 1 = normalized (default)
Option5 = Maximum number of iterations
Option6 = Tolerance of residual changes
Option7 = Tolerance of coefficient changes
Other solving parameters are available. See the TBFT command for details.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
181
where:
ID = Index corresponding to the material number
Option2 = Creep model abbreviation
Option3 = (Blank--not applicable)
Time Hardening
182
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Generalized Graham
Generalized Blackburn
Modified Time
Hardening
Modified Strain
Hardening
Generalized Garofalo
Exponential Form
Norton
Prim+Sec Time
Hardening
Prim+Sec Rational
Polynomial
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
183
Comments
184
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Comments
185
3.16E+02
3.23E+02
7.87E+01
-1.64E+02
-2.31E+02
-2.66E+02
-2.84E+02
-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.
where:
ID = Index corresponding to the material number
Option1 = Experimental data type UNIA (uniaxial test data)
Option2 = Experimental data file name
Option3 = File name extension
Option4 = File directory
186
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
! define case
where:
ID = Index corresponding to the material number
Option2 = CPLA
Option3 = Your specified case name
TBFT,FADD,ID,Category,Option2,Option3
where:
ID = Index corresponding to the material number
Category = PLAS
Option2 = One of the following: CHABOCHE, BISO, MISO, VOCE, or POWER
Option3 = Dependent on Option2, as follows:
When Option2 = CHABOCHE, Option3 = 1 to N
When Option2 = MISO, Option3 = 1 to Niso
When Option2 = BISO, VOCE, or POWER, Option3 is not used
TBFT,FCASE,ID,FINISH
! create case
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
187
188
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
where:
ID = Index corresponding to the material number
Option2 = Case name
Option3 = (Blank--not applicable)
Option4 = Index of coefficient
Option5 = Value of coefficient
Example 5.3: Initialize Coefficients
TBFT,SET,1,case,case1,,1,1.2
TBFT,SET,1,case,case1,,2,1.5
By default, coefficients are not fixed. To fix a coefficient to a value (p. 188) 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
where:
ID = Index corresponding to the material number
Option2 = Case name
Option3 = (Blank--not applicable)
Option4 = Index of coefficient
Option5 = 1 to fix, 0 to vary (default)
Temperature dependency uses Option4 and references your specified data files with the appropriate
temp header:
TBFT,SET,ID,CASE,Option2,Option3,Option4,Option5
where:
ID = Index corresponding to the material number
Option2 = Case name
Option3 = (Blank--not applicable)
Option4 = tdep or tref
Option5 = If Option4 = tdep, then 1 activates temperature dependency 0 deactivates it. If Option4 = tref, this value is either a specific temperature or all temperatures (ALL (p. 188)).
Example 5.4: Fixing Coefficient Values
TBFT,FIX,1,case,case1,,1,1
TBFT,FIX,1,case,case1,,2,1
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
189
where:
ID = Index corresponding to the material number
Option2 = Case name
Option3 = (Blank--not applicable)
Option4 = tdep
Option5 = 1 to activate temperature dependency, 0 to deactivate (default)
TBFT,SET,ID,CASE,Option2,Option3,Option4,Option5
where:
ID = Index corresponding to the material number
Option2 = Case name
Option3 = (Blank--not applicable)
Option4 = tref
Option5 = Valid temperature values found in the experimental data
The SOLVE option allows you to specify procedure types, tolerances, and the number of iterations:
TBFT,SOLVE,ID,CASE,Option2,Option3,Option4, ... , Option7
190
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
191
192
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
With ...
Combination
Type
Command, Label
Link to Example
Plasticity
Combined
Hardening
Bilinear
TB,BISO + TB,CHAB
Plasticity
Combined
Hardening
Multilinear
TB,MISO + TB,CHAB
Plasticity
Combined
Hardening
Multilinear
TB, PLAS,,,,MISO +
TB,CHAB
Plasticity
Combined
Hardening
Nonlinear
TB,NLISO + TB,CHAB
Plasticity
Combined
Hardening
Multilinear
TB,PLAS,,,,MISO +
TB,EDP
Plasticity
Combined
Hardening
Multilinear
TB,MISO+ TB,EDP
Viscoplasticity
Isotropic
Hardening
Bilinear
TB,PLAS,,,,BISO +
TB,EDP
Viscoplasticity
Isotropic
Hardening
Multilinear
TB,MISO + TB,RATE
Viscoplasticity
Isotropic
Hardening
Multilinear
TB,PLAS,,,,MISO +
TB,RATE
Viscoplasticity
Isotropic
Hardening
Nonlinear
TB,NLISO + TB,RATE
Viscoplasticity
Combined
Hardening
Nonlinear
TB,CHAB+
TB,RATE+TB,BISO
Viscoplasticity
Combined
Hardening
Nonlinear
TB,CHAB+
TB,RATE+TB,MISO
Viscoplasticity
Combined
Hardening
Nonlinear
TB,CHAB+
TB,RATE+TB,PLASTIC
Viscoplasticity
Combined
Hardening
Nonlinear
TB,CHAB+
TB,RATE+TB,NLISO
Gurson Plasticity
Isotropic
Hardening
Bilinear
TB,GURS + TB,BISO
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
193
With ...
Combination
Type
Command, Label
Link to Example
Gurson Plasticity
Isotropic
Hardening
Multilinear
TB,GURS + TB,MISO
Gurson Plasticity
Isotropic
Hardening
Multilinear
TB,GURS + TB,PLAS,,,
MISO
Gurson Plasticity
Isotropic
Hardening
Nonlinear
TB,GURS + TB,NLISO
Gurson Plasticity
Isotropic
Hardening
Nonlinear
TB,GURS + TB,CHAB
Gurson Plasticity
Isotropic
Hardening
Nonlinear
TB,GURS + TB,CHAB
+ TB,BISO or MISO or
NLISO, or
TB,PLAS,,,,MISO
Plasticity and
Creep (Implicit)
Isotropic
Hardening
Bilinear
TB,BISO + TB,CREEP
Plasticity and
Creep (Implicit)
Isotropic
Hardening
Multilinear
TB,MISO + TB,CREEP
Plasticity and
Creep (Implicit)
Isotropic
Hardening
Multilinear
TBPLAS,,,,MISO +
TB,CREEP
Plasticity and
Creep (Implicit)
Isotropic
Hardening
Nonlinear
TB,NLISO + TB,CREEP
Plasticity and
Creep (Implicit)
Kinematic
Hardening
Bilinear
TB,BKIN + TB,CREEP
Plasticity and
Creep (Implicit)
Kinematic
Hardening
Nonlinear
TB,CHAB + TB,CREEP
Plasticity and
Creep (Implicit)
Combined
Hardening
Nonlinear
TB,CHAB + TB,CREEP
+ TB,BISO or TB,MISO
or TB,NLISO or
TB,PLAS,,,,MISO
Anisotropic
Plasticity
Isotropic
Hardening
Bilinear
TB,HILL + TB,BISO
Anisotropic
Plasticity
Isotropic
Hardening
Multilinear
TB,HILL + TB,MISO
Anisotropic
Plasticity
Isotropic
Hardening
Multilinear
TB,HILL +
TBPLAS,,,,MISO
Anisotropic
Plasticity
Isotropic
Hardening
Nonlinear
TB,HILL + TB,NLSIO
Anisotropic
Plasticity
Kinematic
Hardening
Bilinear
TB,HILL + TB,BKIN
Anisotropic
Plasticity
Kinematic
Hardening
Multilinear
TB,HILL + TB,MKIN
Anisotropic
Plasticity
Kinematic
Hardening
Multilinear
Anisotropic
Plasticity
Kinematic
Hardening
Multilinear
TB,HILL + TBPLAS,,,,
KINH
194
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Model
With ...
Anisotropic
Plasticity
Kinematic
Hardening
Chaboche
TB,HILL + TB,CHAB
Anisotropic
Plasticity
Combined
Hardening
Bilinear
Isotropic
and
Chaboche
TB,HILL + TB,BISO +
TB,CHAB
Anisotropic
Plasticity
Combined
Hardening
Multilinear
Isotropic
and
Chaboche
TB,HILL + TB,MISO +
TB,CHAB
Anisotropic
Plasticity
Combined
Hardening
Multilinear
Isotropic
and
Chaboche
TB,HILL +
TB,PLAS,,,,MISO +
TB,CHAB
Anisotropic
Plasticity
Combined
Hardening
Nonlinear
Isotropic
and
Chaboche
TB,HILL + TB,NLISO +
TB,CHAB
Anisotropic
Viscoplasticity
Isotropic
Hardening
Bilinear
TB,HILL + TB,RATE +
TB,BISO
Anisotropic
Viscoplasticity
Isotropic
Hardening
Multilinear
TB,HILL + TB,RATE +
TB,MISO
Anisotropic
Viscoplasticity
Isotropic
Hardening
Nonlinear
TB,HILL + TB,RATE +
TB,NLISO
TB,HILL + TB,CREEP
Anisotropic
Creep (Implicit)
Anisotropic
Creep and
Plasticity
(Implicit)
--Isotropic
Hardening
Anisotropic
Creep and
Plasticity
(Implicit)
Isotropic
Hardening
Anisotropic
Creep and
Plasticity
(Implicit)
Isotropic
Hardening
Anisotropic
Creep and
Plasticity
(Implicit)
Isotropic
Hardening
Anisotropic
Creep and
Plasticity
(Implicit)
Combination
Type
Kinematic
Hardening
Bilinear
Command, Label
TB,HILL + TB,CREEP
+ TB,BISO
Multilinear
TB,HILL + TB,CREEP
+ TB,MISO
Multilinear
TB,HILL + TB,CREEP
+ TB,MISO
Nonlinear
TB,HILL + TB,CREEP
+ TB,NLISO
Bilinear
TB,HILL + TB,CREEP
+ TB,BKIN
Link to Example
HILL and CHAB
Example
HILL and BISO and
CHAB Example
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
195
With ...
Anisotropic
Creep and
Plasticity
(Implicit)
Kinematic
Hardening
Anisotropic
Creep and
Plasticity
(Implicit)
Combined
Hardening
Combination
Type
Command, Label
Link to Example
Chaboche
TB,HILL + TB,CHAB +
TB,CREEP
Chaboche
TB,HILL + TB,CHAB +
TB,CREEP
TB,BISO or TB,MISO or
TB,NLISO or
TB,PLAS,,,,MISO
Anisotropic
Hyperelasticity
and
Viscoelasticity
Finite Strain
Nonlinear
Visco-AnisotropicAnisotropic
Hyperelasticity Elasticity
TB,AHYPER +
TB,PRONY
Hyperelasticity
and
Viscoelasticity
Nonlinear
Finite Strain
Isotropic
Visco-Hyperelasticity
Elasticity
TB,HYPER +
TB,PRONY
Extended
Drucker-Prager
(EDP) and Creep
(Implicit)
Isotropic
Hardening
Bilinear,
Multilinear,
or Nonlinear
TB,EDP + TB,CREEP +
TB,BISO or TB,MISO
or TB,NLISO
Isotropic
Hardening
Bilinear,
Multilinear,
or Nonlinear
TB,EDP + TB,CREEP +
TB,BISO or TB,MISO
or TB,NLISO
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.
Bilinear Isotropic Hardening (TB,BISO)
Bilinear Kinematic Hardening (TB,BKIN)
Chaboche Nonlinear Kinematic Hardening (TB,CHAB)
Creep (Implicit) (TB,CREEP) -- Also see Creep in the Structural Analysis Guide.
Hill Anisotropy (TB,HILL]
Multilinear Isotropic Hardening (TB,MISO)
Multilinear Kinematic Hardening (TB,MKIN or KINH)
Nonlinear Isotropic Hardening (TB,NLISO).
Rate-Dependent Plasticity (TB,RATE)
Hyperelasticity (TB,HYPER)
Anisotropic Hyperelasticity (TB, AHYPER) -- Also see Anisotropic Hyperelasticity in the Mechanical APDL
Theory Reference.
196
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Viscoelasticity
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
197
198
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
199
For more information, see Subroutine UserFld (Update User-Defined Field Variables) in the Programmer's
Reference
You can use any appropriate API subroutine in the Programmer's Reference to query and update the
user-defined field variables.
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Data Processing
TB,FRIC,1,2, ,ISO
TBFIELD,TEMP,100.0
TBFIELD,SLDI,0.1
TBDATA,1,0.3
TBFIELD,SLDI,0.5
TBDATA,1,0.5
TBFIELD,TEMP,200.0
TBFIELD,SLDI,0.2
TBDATA,1,0.2
TBFIELD,SLDI,0.7
TBDATA,1,0.1
= 100.00
= 0.10000
FRICTION DIR.
1
TEMPERATURE
SLIDE DIST
FRICTION COEFF.
0.30000
= 100.00
= 0.50000
FRICTION DIR.
1
TEMPERATURE
SLIDE DIST
FRICTION COEFF.
0.50000
= 200.00
= 0.20000
FRICTION DIR.
1
TEMPERATURE
SLIDE DIST
FRICTION COEFF.
0.20000
= 200.00
= 0.70000
FRICTION DIR.
1
FRICTION COEFF.
0.10000
0.3
0.5
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):
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
201
200
x0 = 0.2,
y0 = 0.35
y1 = 0.5
x1 = 0.5
and solving for the interpolated values using Equation (1), we obtain
202
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
x0 = 0.2,
y0 = 0.2
y1 = 0.14
x1 = 0.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.
t = 180,
t0 = 100,
y0 = 0.45
y1 = 0.16
t1 = 200
and solving for the interpolated values using Equation (1), we obtain
Equation (8)
203
The figure shows how the data points are projected first in one dimension and then in the second dimension. This method was extrapolated and implemented for multiple dimensions.
Implementation of the radial basis function is expressed as:
where N is the number of data points and O is the number of free variables (or the order of the interpolation). Input data is (xj,1, xj,2, , zj) where j varies from 1 to the N.
The unknown values are ai (where i varies from 1 to N ) and c. The equation is evaluated for all data
points provided in the input to calculate the ai and c values.
Reference
For further information about multidimensional field-variable interpolation, consult this reference:
1. Amidror, Isaac.Scattered Data Interpolation Methods for Electronic Imaging Systems: A Survey . Journal
of Electronic Imaging 11:2 (2002: 157-176.
204
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Material Property
AFDM
Acoustic frequency
dependent material
AHYPER
Anisotropic hyperelasticity
BB
Bergstrom-Boyce
CDM
Mullins effect
CGCR
CTE
Coefficient of thermal
expansion
CZM
Cohesive zone
DMGE
DMGI
EDP
Extended Drucker-Prager
ELASTIC
Elasticity
EXPE
Experimental data
FCLI
FLUID
Fluid
FRIC
Coefficient of friction
GURSON
Gurson pressure-dependent
plasticity
INTER
User-defined contact
interaction
MPLANE
Microplane
PERF
PM
Porous media
WEAR
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
205
206
Release 16.0 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.