8.2. Material Model Combination Examples

You can combine several material model options to simulate complex material behaviors. Material Model Combinations presents the model options you can combine along with the associated TB command labels and links to example input listings.

The following example input listings are presented in sections identified via the TB command labels (Lab values).

 

8.2.1. RATE and CHABOCHE and PLASTIC (BISO) Example

Combining rate-dependent plasticity (viscoplasticity), Chaboche nonlinear kinematic hardening plasticity, and bilinear isotropic hardening plasticity:

MP,EX,1,185.0E3     ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,RATE,1,,,PERZYNA              ! RATE TABLE
TBDATA,1,0.5,1 

TB,CHAB,1                        ! CHABOCHE TABLE
TBDATA,1,180,100,3

TB,PLAS,1,,,BISO                 ! PLASTIC(BISO) TABLE
TBDATA,1,180,200

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

 

8.2.2. RATE and CHABOCHE and PLASTIC (MISO) and PLASTIC (KSR) Example

Combining multilinear isotropic hardening plasticity, rate-dependent plasticity (viscoplasticity), and Chaboche nonlinear kinematic hardening plasticity with kinematic static recovery: 

MP,EX,1,79650e6
MP,NUXY,1,0.33

TB,CHAB,1,,2,TRATE              ! DEFINE CHABOCHE MATERIAL DATA
TBDATA,1,1.5e8
TBDATA,2,62511e7,2000
TBDATA,4,62511e6,1000

TB,RATE,1,,,PERZYNA             ! RATE TABLE
TBDATA,1,0.1,2
 
TB,PLAS,1, , ,MISO              ! PLASTIC data table
TBPT,DEFI,0,1.5e8     
TBPT,DEFI,0.01,2.8e8
TBPT,DEFI,0.05,4e8
TBPT,DEFI,0.1,4.1e8

TB,PLASTIC,1,,2,KSR2            ! Kinematic hardening static recovery
TBDATA,1,1e-12,2.1,
TBDATA,3,1e-12,2.15,

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

For information about the PLASTIC (MISO and KSR) options, see Multilinear Isotropic Hardening and Kinematic Hardening Static Recovery.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.3. RATE and CHABOCHE and PLASTIC (MISO) Example

Combining multilinear isotropic hardening plasticity, rate-dependent plasticity (viscoplasticity), and Chaboche nonlinear kinematic hardening plasticity:

MP,EX,1,185E3                    ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,RATE,1,,,PERZYNA              ! RATE TABLE
TBDATA,1,0.5,1 

TB,CHAB,1                        ! CHABOCHE TABLE
TBDATA,1,180,100,3               ! THIS EXAMPLE ISOTHERMAL

TB,PLAS,,,,MISO                  ! PLASTIC (MISO) TABLE
TBPT,,0.0,180
TBPT,,0.99795,380

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.4. RATE and CHABOCHE and NLISO Example

Combining viscoplasticity and Chaboche nonlinear kinematic hardening plasticity and nonlinear isotropic hardening plasticity:

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,RATE,1,,,PERZYNA              ! RATE TABLE
TBDATA,1,0.5,1 

TB,CHAB,1,3,5                    ! CHABOCHE TABLE
TBTEMP,20                        ! THIS EXAMPLE TEMPERATURE DEPENDENT
TBDATA,1,500,20000,100,40000,200,10000
TBDATA,7,1000,200,100,100,0
TBTEMP,40
TBDATA,1,880,204000,200,43800,500,10200
TBDATA,7,1000,2600,2000,500,0
TBTEMP,60
TBDATA,1,1080,244000,400,45800,700,12200
TBDATA,7,1400,3000,2800,900,0

TB,NLISO,1,2                     ! NLISO TABLE
TBTEMP,40
TBDATA,1,880,0.0,80.0,3
TBTEMP,60
TBDATA,1,1080,0.0,120.0,7

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

For information about the NLISO option, see Nonlinear Isotropic Hardening.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.5. PLASTIC (BISO) and CHABOCHE Example

Combining bilinear isotropic hardening plasticity with Chaboche nonlinear kinematic hardening plasticity:

MP,EX,1,185.0E3                  ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,PLAS,1,,,BISO                 ! PLASTIC (BISO) TABLE
TBDATA,1,180,200

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.6. PLASTIC (MISO) and CHABOCHE Example

Combining multilinear isotropic hardening and Chaboche nonlinear kinematic hardening:

MP,EX,1,185E3                     ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,CHAB,1                         ! CHABOCHE TABLE
TBDATA,1,180,100,3                ! THIS EXAMPLE ISOTHERMAL

TB,PLAS,,,,MISO                   ! PLASTIC (MISO) TABLE
TBPT,,0.0,180
TBPT,,0.99795,380

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.7. NLISO and CHABOCHE Example

Combining nonlinear isotropic hardening plasticity with Chaboche nonlinear kinematic hardening plasticity:

MP,EX,1,20.0E5                    ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,CHAB,1,3,5                     ! CHABOCHE TABLE
TBTEMP,20                         ! THIS EXAMPLE TEMPERATURE DEPENDENT
TBDATA,1,500,20000,100,40000,200,10000
TBDATA,7,1000,200,100,100,0
TBTEMP,40
TBDATA,1,880,204000,200,43800,500,10200
TBDATA,7,1000,2600,2000,500,0
TBTEMP,60
TBDATA,1,1080,244000,400,45800,700,12200
TBDATA,7,1400,3000,2800,900,0

TB,NLISO,1,2                       ! NLISO TABLE
TBTEMP,40
TBDATA,1,880,0.0,80.0,3
TBTEMP,60
TBDATA,1,1080,0.0,120.0,7

For information about the NLISO option, see Nonlinear Isotropic Hardening.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.8. PLASTIC (MISO) and EDP Example

Combining multilinear isotropic hardening Extended Drucker-Prager plasticity:

/prep7
mp,ex,1,2.1e4			! Elastic Properties
mp,nuxy,1,0.1

ys=7.894657
sl=1000.0

tb,edp,1,,,LYFUN
tbdata,1,2.2526,ys

tb,edp,1,,,LFPOT
tbdata,1,0.566206

tb,plas,1,1,2,miso
tbpt,defi,0.0,7.894
tbpt,defi,1,1007.894

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

For information about the EDP option, see Extended Drucker-Prager (EDP).

 

8.2.9. GURSON and PLASTIC (BISO) Example

Combining bilinear isotropic hardening with Gurson plasticity: 

q1=1.5
q2=1
q3=q1*q1
sigma_Y=E/300.0

f_0= 0.04
f_N= 0.04
S_N=0.1
strain_N=0.3

tb,GURS,1,,5,BASE            ! Gurson's BASE model
tbdata,1,sigma_Y,f_0,q1,q2,q3

tb,GURS,1,,3,SNNU            ! Gurson's SNNU model
tbdata,1,f_N,strain_N,S_N

TB,PLAS,1,,,BISO             ! PLASTIC (BISO) TABLE
TBDATA,1,sigma_Y,0.1         ! Isotropic hardening yield stress definition overrides

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

For information about the GURSON option, see Gurson.

 

8.2.10. GURSON and PLASTIC (MISO) Example

Combining multilinear isotropic hardening with Gurson plasticity: 

q1=1.5
q2=1
q3=q1*q1
sigma_Y=E/300.0
f_0= 0.04
f_N= 0.04
S_N=0.1
strain_N=0.3

tb,GURS,1,,5,BASE                ! Gurson's BASE model
tbdata,1,sigma_Y,f_0,q1,q2,q3

tb,GURS,1,,3,SNNU                ! Gurson's SNNU model
tbdata,1,f_N,strain_N,S_N

tb,plas,1,,4,miso
tbpt, defi, 0.0, sigma_Y
tbpt, defi, 1, 10.0* sigma_Y
! Isotropic hardening yield stress definition overrides

For information about the GURSON option, see Gurson.

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

 

8.2.11. NLISO and GURSON Example

Combining nonlinear isotropic hardening with Gurson plasticity: 

q1=1.5
q2=1
q3=q1*q1
sigma_Y=E/300.0

f_0= 0.04
f_N= 0.04
S_N=0.1
strain_N=0.3
Power_N=0.1

tb,GURS,1,,5,BASE                ! Gurson's BASE model
tbdata,1,sigma_Y,f_0,q1,q2,q3

tb,GURS,1,,3,SNNU                ! Gurson's SNNU model
tbdata,1,f_N,strain_N,S_N

! Isotropic hardening yield stress definition overrides
tb,nliso,1,1,2,POWER
tbdata,1,sigma_Y,power_N

For information about the NLISO option, see Nonlinear Isotropic Hardening.

For information about the GURSON option, see Gurson.

 

8.2.12. GURSON and CHABOCHE Example

Combining Gurson plasticity with Chaboche nonlinear kinematic hardening plasticity:

! Gurson coefficients
Q1=1.5e-16                   ! FIRST TVERGAARD CONSTANT
Q2=1                         ! SECOND TVERGAARD CONSTANT
Q3=Q1*Q1                     ! THIRD TVERGAARD CONSTANT
F_0=0.04                     ! INITIAL POROSITY
F_N=0.005e-12                ! VOL. FRAC. OF VOID NUCL. PARTICLES
S_N=0.05                     ! STAND. DEV. OF MEAN STRN FOR NUCLEA.
STRAIN_N=0.1e2               ! MEAN STRAIN FOR NUCLEATIONS

TB,GURS,1,,5,BASE            ! BASE DEFINED
TBDATA,1,320.,F_0,Q1,Q2,Q3

TB,GURS,1,,3,SNNU            ! SNNU DEFINED
TBDATA,1,F_N,STRAIN_N,S_N

! Chaboche kinematic hardening
TB,CHAB,1,,2
TBDATA,1,320.,150000.,300.,1500.,0.1

For information about the GURSON option, see Gurson.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.13. GURSON and CHABOCHE and PLASTIC (BISO) Example

Combining Gurson plasticity with Chaboche nonlinear kinematic hardening plasticity and bilinear isotropic hardening:

! Gurson coefficients
Q1=1.5e-16                   ! FIRST TVERGAARD CONSTANT
Q2=1                         ! SECOND TVERGAARD CONSTANT
Q3=Q1*Q1                     ! THIRD TVERGAARD CONSTANT
F_0=0.04                     ! INITIAL POROSITY
F_N=0.005e-12                ! VOL. FRAC. OF VOID NUCL. PARTICLES
S_N=0.05                     ! STAND. DEV. OF MEAN STRN FOR NUCLEA.
STRAIN_N=0.1e2               ! MEAN STRAIN FOR NUCLEATIONS

TB,GURS,1,,5,BASE            ! BASE DEFINED
TBDATA,1,320.,F_0,Q1,Q2,Q3

TB,GURS,1,,3,SNNU            ! SNNU DEFINED
TBDATA,1,F_N,STRAIN_N,S_N

! Chaboche kinematic hardening
TB,CHABOCHE,1,,2
TBDATA,1,320.,150000.,300.,1500.,0.1

! Bilinear isotropic hardening
TB,PLAS,1,,,BISO
TBDATA,1,320.,0.1

For information about the GURSON option, see Gurson.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

 

8.2.14. RATE and PLASTIC (BISO) Example

Combining bilinear isotropic hardening plasticity with rate-dependent plasticity (viscoplasticity):

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,PLAS,1,,,BISO                 ! PLASTIC (BISO) TABLE
TBDATA,1,9000,10000

TB,RATE,1,,,PERZYNA              ! RATE TABLE
TBDATA,1,0.5,1

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

 

8.2.15. RATE and PLASTIC (MISO) Example

Combining multilinear isotropic hardening plasticity with rate-dependent plasticity (viscoplasticity):

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,PLAS,,,,MISO                  ! MISO TABLE
TBPT,,0.00000,30000
TBPT,,4.00E-3,32000
TBPT,,8.10E-3,33800
TBPT,,1.25E-2,35000
TBPT,,2.18E-2,36500
TBPT,,3.10E-2,38000
TBPT,,4.05E-2,39000

TB,RATE,1,,,PERZYNA              ! RATE TABLE
TBDATA,1,0.5,1

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

 

8.2.16. RATE and NLISO Example

Combining nonlinear isotropic hardening plasticity with rate-dependent plasticity (viscoplasticity):

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,NLISO,1                       ! NLISO TABLE
TBDATA,1,30000,100000,5200,172

TB,RATE,1,,,PERZYNA              ! RATE TABLE
TBDATA,1,0.5,1

For information about the NLISO option, see Nonlinear Isotropic Hardening.

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

 

8.2.17. PLASTIC (BISO) and CREEP Example

Combining bilinear isotropic hardening plasticity with implicit creep:

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,PLAS,1,,,BISO                 ! PLASTIC (BISO) TABLE
TBDATA,1,9000,10000

TB,CREEP,1,,,2                   ! CREEP TABLE
TBDATA,1,1.5625E-14,5.0,-0.5,0.0

Combining bilinear isotropic hardening plasticity with explicit creep:

MP,EX,1,20000
MP,NUXY,1,0.3

TB,PLAS,1,,,BISO
TBDATA,1,1000,6500

TB,CREEP,1,,,0
TBDATA,1,3.125e-13,5,2,100,1e-5,2

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

For information about the CREEP option, see Creep Model and Creep Option.

 

8.2.18. PLASTIC (MISO) and CREEP Example

Combining the multilinear isotropic hardening option with implicit creep:

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,PLAS,,,,MISO                  ! MISO TABLE
TBPT,,0.00000,30000
TBPT,,4.00E-3,32000
TBPT,,8.10E-3,33800
TBPT,,1.25E-2,35000
TBPT,,2.18E-2,36500
TBPT,,3.10E-2,38000
TBPT,,4.05E-2,39000

TB,CREEP,1,,,2                   ! CREEP TABLE
TBDATA,1,1.5625E-14,5.0,-0.5,0.0

Combining the multilinear isotropic hardening option with explicit creep:

MP,EX,1,20000
MP,NUXY,1,0.3

TB,PLAS,1,,,MISO
TBPT,DEFINE,0,1000
TBPT,DEFINE,0.2,2000

TB,CREEP,1,,,0
TBDATA,1,3.125e-20,5,1,0,,1

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

For information about the CREEP option, see Creep Model and Creep Option.

 

8.2.19. PLASTIC (KINH) and CREEP Example

Combining the multilinear kinematic hardening option with explicit creep:

MP,EX,1,20000
MP,NUXY,1,0.3

TB,PLAS,1,,,KINH
TBPT,DEFINE,0,1000
TBPT,DEFINE,0.2,2000

TB,CREEP,1,,,0
TBDATA,1,3.125e-20,5,1,0,,1

For information about the PLASTIC (KINH) option, see Multilinear Kinematic Hardening.

For information about the CREEP option, see Creep Model and Creep Option.

 

8.2.20. NLISO and CREEP Example

Combining nonlinear isotropic hardening plasticity with implicit creep:

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,NLISO,1                       ! NLISO TABLE
TBDATA,1,30000,100000,5200,172

TB,CREEP,1,,,2                   ! CREEP TABLE
TBDATA,1,1.5625E-14,5.0,-0.5,0.0

For information about the NLISO option, see Nonlinear Isotropic Hardening.

For information about the CREEP option, see Creep Model and Creep Option.

 

8.2.21. PLASTIC (BKIN) and CREEP Example

Combining bilinear kinematic hardening plasticity with implicit creep:

MP,EX,1,1e7                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.32 

TB,PLAS,1,,,BKIN              ! PLASTIC (BKIN) TABLE
TBDATA,1,42000,1000

TB,CREEP,1,,,6                ! CREEP TABLE
TBDATA,1,7.4e-21,3.5,0,0,0,0

Combining bilinear kinematic hardening plasticity with explicit creep:

MP,EX,1,20000 
MP,NUXY,1,0.3 

TB,PLAS,1,,,BKIN 
TBDATA,1,1000,6500 

TB,CREEP,1,,,0 
TBDATA,1,3.125e-13,5,2,100,1e-5,2

For information about the PLASTIC (BKIN) option, see Bilinear Kinematic Hardening.

For information about the CREEP option, see Creep Model and Creep Option.

 

8.2.22. CHABOCHE and PLASTIC (KSR2) and CREEP Example

Combining Chaboche nonlinear kinematic hardening with static recovery and implicit creep:

YOUNGS = 30e3       ! Young’s Modulus
NU     = 0.3        ! Poisson’s ratio
SIGMA0 = 18.0E0     ! Initial yield stress
				                       
mp,ex,1,YOUNGS
mp,prxy,1,NU

TB,CHAB,1,1,1,TRATE  ! Chaboche kinematic hardening
TBDATA,1,sigma0,
TBDATA,2,1e+4,0  

TB,CREEP,1,1, ,1,1   ! Creep model
tbdata,1,            ! No creep strain

TB,PLASTIC,1,,1,KSR2 ! Kinematic hardening static recovery
TBDATA,1,1e-3,2.0,

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

For more information about the PLASTIC (KSR2) option, see Kinematic Hardening Static Recovery.

For information about the CREEP option, see Creep Model and Creep Option.

 

8.2.23. HILL and PLASTIC (BISO) Example

Combining anisotropic plasticity with bilinear isotropic hardening.

MP,EX,1,20.0E5                        ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,HILL,1,2                           ! HILL TABLE
TBTEMP,100
TBDATA,1,1,1.0402,1.24897,1.07895,1,1
TBTEMP,200
TBDATA,1,0.9,0.94,1.124,0.97,0.9,0.9

TB,PLAS,1,,,BISO                      ! PLASTIC (BISO) TABLE
TBTEMP,100
TBDATA,1,461.0,374.586
TBTEMP,200
TBDATA,1,400.0,325.0

For information about the HILL option, see Hill Yield Criterion.

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

 

8.2.24. HILL and PLASTIC (MISO) Example

Combining multilinear isotropic hardening with Hill anisotropic plasticity:

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,PLAS,,,,MISO                  ! MISO TABLE
TBPT,,0.00000,30000
TBPT,,4.00E-3,32000
TBPT,,8.10E-3,33800
TBPT,,1.25E-2,35000
TBPT,,2.18E-2,36500
TBPT,,3.10E-2,38000
TBPT,,4.05E-2,39000

TB,HILL,1                        ! HILL TABLE
TBDATA,1,1.0,1.1,0.9,0.85,0.9,0.80

For information about the HILL option, see Hill Yield Criterion.

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

 

8.2.25. HILL and NLISO Example

Combining anisotropic plasticity with nonlinear isotropic hardening:

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,NLISO,1                       ! NLISO TABLE
TBDATA,1,30000,100000,5200,172

TB,HILL,1                        ! HILL TABLE
TBDATA,1,1.0,1.1,0.9,0.85,0.9,0.80

For information about the HILL option, see Hill Yield Criterion.

For information about the NLISO option, see Nonlinear Isotropic Hardening.

 

8.2.26. HILL and PLASTIC (BKIN) Example

Combining anisotropic plasticity with bilinear kinematic hardening:

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,PLAS,1,,,BKIN                 ! PLASTIC (BKIN) TABLE
TBDATA,1,9000,10000

TB,HILL,1                        ! HILL TABLE
TBDATA,1,1.0,1.1,0.9,0.85,0.9,0.80

For information about the HILL option, see Hill Yield Criterion.

For information about the PLASTIC (BKIN) option, see Bilinear Kinematic Hardening.

 

8.2.27. HILL and CHABOCHE Example

Combining anisotropic plasticity with Chaboche nonlinear kinematic hardening:

MP,EX,1,185E3                    ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,CHAB,1                        ! CHABOCHE TABLE
TBDATA,1,180,400,3,0

TB,HILL,1                        ! HILL TABLE
TBDATA,1,1.0,1.1,0.9,0.85,0.9,0.80

For information about the HILL option, see Hill Yield Criterion.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.28. HILL and PLASTIC (BISO) and CHABOCHE Example

Combining anisotropic plasticity with bilinear isotropic hardening and Chaboche nonlinear kinematic hardening:

MP,EX,1,185E3                       ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,CHAB,1                           ! CHABOCHE TABLE
TBDATA,1,180,100,3

TB,PLAS,1,,,BISO                    ! PLASTIC (BISO) TABLE
TBDATA,1,180,200

TB,HILL,1                           ! HILL TABLE
TBDATA,1,1.0,1.1,0.9,0.85,0.9,0.80

For information about the HILL option, see Hill Yield Criterion.

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.29. HILL and PLASTIC (MISO) and CHABOCHE Example

Combining multilinear isotropic hardening option with Hill anisotropic plasticity and Chaboche nonlinear kinematic hardening:

MP,EX,1,185E3                       ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,CHAB,1                           ! CHABOCHE TABLE
TBDATA,1,185,100,3

TB,PLAS,,,,MISO                     ! MISO TABLE
TBPT,,0.001,185
TBPT,,0.998,380

TB,HILL,1                           ! HILL TABLE
TBDATA,1,1.0,1.1,0.9,0.85,0.9,0.80

For information about the HILL option, see Hill Yield Criterion.

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.30. HILL and NLISO and CHABOCHE Example

Combining anisotropic plasticity with nonlinear isotropic hardening and Chaboche nonlinear kinematic hardening:

MPTEMP,1,20,200,400,550,600,650            ! ELASTIC CONSTANTS
MPTEMP,,700,750,800,850,900,950
!
MPDATA,EX,1,1,1.250E4,1.210E4,1.140E4,1.090E4,1.070E4,1.050E4 
MPDATA,EX,1,,1.020E4,0.995E4,0.963E4,0.932E4,0.890E4,0.865E4 
!
MPDATA,EY,1,1,1.250E4,1.210E4,1.140E4,1.090E4,1.070E4,1.050E4 
MPDATA,EY,1,,1.020E4,0.995E4,0.963E4,0.932E4,0.890E4,0.865E4 
!
MPDATA,EZ,1,1,1.250E4,1.210E4,1.140E4,1.090E4,1.070E4,1.050E4
MPDATA,EZ,1,,1.020E4,0.995E4,0.963E4,0.932E4,0.890E4,0.865E4
!
MPDATA,PRXY,1,1,0.351,0.359,0.368,0.375,0.377,0.380 
MPDATA,PRXY,1,,0.382,0.384,0.386,0.389,0.391,0.393
!
MPDATA,PRYZ,1,1,0.351,0.359,0.368,0.375,0.377,0.380 
MPDATA,PRYZ,1,,0.382,0.384,0.386,0.389,0.391,0.393
!
MPDATA,PRXZ,1,1,0.351,0.359,0.368,0.375,0.377,0.380
MPDATA,PRXZ,1,,0.382,0.384,0.386,0.389,0.391,0.393
!
MPDATA,GXY,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4,1.040E4
MPDATA,GXY,1,,1.020E4,1.000E4,0.973E4,0.946E4,0.908E4,0.887E4
!	      
MPDATA,GYZ,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4,1.040E4
MPDATA,GYZ,1,,1.020E4,1.000E4,0.973E4,0.946E4,0.908E4,0.887E4
!								      
MPDATA,GXZ,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4,1.040E4
MPDATA,GXZ,1,,1.020E4,1.000E4,0.973E4,0.946E4,0.908E4,0.887E4


TB,NLISO,1                                ! NLISO TABLE
TBDATA,1,180,0.0,100.0,5

!

TB,CHAB,1                                 ! CHABOCHE TABLE
TBDATA,1,180,100,3


TB,HILL,1,5                               ! HILL TABLE
TBTEMP,750.0
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,800.0				 
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,850.0				 
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,900.0				 
TBDATA,1,1.0,1.0,1.0,1.00,1.00,1.00
TBTEMP,950.0				 
TBDATA,1,1.0,1.0,1.0,1.00,1.00,1.00

For information about the HILL option, see Hill Yield Criterion.

For information about the NLISO option, see Nonlinear Isotropic Hardening.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.31. HILL and RATE and PLASTIC (BISO) Example

Combining anisotropic viscoplasticity with bilinear isotropic hardening plasticity:

MPTEMP,1,20,400,650,800,950      ! ELASTIC CONSTANTS
!
MPDATA,EX,1,1,30.00E6,27.36E6,25.20E6,23.11E6,20.76E6
!
MPDATA,EY,1,1,30.00E6,27.36E6,25.20E6,23.11E6,20.76E6
!
MPDATA,EZ,1,1,30.00E6,27.36E6,25.20E6,23.11E6,20.76E6
!
MPDATA,PRXY,1,1,0.351,0.359,0.368,0.375,0.377 
!
MPDATA,PRYZ,1,1,0.351,0.359,0.368,0.375,0.377
!
MPDATA,PRXZ,1,1,0.351,0.359,0.368,0.375,0.377
!
MPDATA,GXY,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4
!             
MPDATA,GYZ,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4
!                                                                     
MPDATA,GXZ,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4


TB,PLAS,1,,,BISO                 ! PLASTIC (BISO) TABLE
TBDATA,1,45000,775000


TB,RATE,1,2,,PERZYNA             ! RATE TABLE
TBTEMP,20
TBDATA,1,0.1,0.3
TBTEMP,950
TBDATA,1,0.3,0.5


TB,HILL,1,5                      ! HILL TABLE
TBTEMP,750.0
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,800.0                             
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,850.0                             
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,900.0                             
TBDATA,1,1.0,1.0,1.0,1.00,1.00,1.00
TBTEMP,950.0                             
TBDATA,1,1.0,1.0,1.0,1.00,1.00,1.00

For information about the HILL option, see Hill Yield Criterion.

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

 

8.2.32. HILL and RATE and NLISO Example

Combining anisotropic viscoplasticity with nonlinear isotropic hardening plasticity:

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,NLISO,1                       ! NLISO TABLE
TBDATA,1,30000,100000,5200,172

TB,HILL,1                        ! HILL TABLE
TBDATA,1,1.0,1.1,0.9,0.85,0.9,0.80

TB,RATE,1,,,PERZYNA              ! RATE TABLE
TBDATA,1,0.5,1

For information about the HILL option, see Hill Yield Criterion.

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

For information about the NLISO option, see Nonlinear Isotropic Hardening.

 

8.2.33. HILL and CREEP Example

Combining anisotropic viscoplasticity with implicit creep:

MPTEMP,1,20,200,400,550,600,650           ! ELASTIC CONSTANTS
MPTEMP,,700,750,800,850,900,950
!
MPDATA,EX,1,1,1.250E4,1.210E4,1.140E4,1.090E4,1.070E4,1.050E4 
MPDATA,EX,1,,1.020E4,0.995E4,0.963E4,0.932E4,0.890E4,0.865E4 
!
MPDATA,EY,1,1,1.250E4,1.210E4,1.140E4,1.090E4,1.070E4,1.050E4 
MPDATA,EY,1,,1.020E4,0.995E4,0.963E4,0.932E4,0.890E4,0.865E4 
!
MPDATA,EZ,1,1,1.250E4,1.210E4,1.140E4,1.090E4,1.070E4,1.050E4
MPDATA,EZ,1,,1.020E4,0.995E4,0.963E4,0.932E4,0.890E4,0.865E4
!
MPDATA,PRXY,1,1,0.351,0.359,0.368,0.375,0.377,0.380 
MPDATA,PRXY,1,,0.382,0.384,0.386,0.389,0.391,0.393
!
MPDATA,PRYZ,1,1,0.351,0.359,0.368,0.375,0.377,0.380 
MPDATA,PRYZ,1,,0.382,0.384,0.386,0.389,0.391,0.393
!
MPDATA,PRXZ,1,1,0.351,0.359,0.368,0.375,0.377,0.380
MPDATA,PRXZ,1,,0.382,0.384,0.386,0.389,0.391,0.393
!
MPDATA,GXY,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4,1.040E4
MPDATA,GXY,1,,1.020E4,1.000E4,0.973E4,0.946E4,0.908E4,0.887E4
!	      
MPDATA,GYZ,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4,1.040E4
MPDATA,GYZ,1,,1.020E4,1.000E4,0.973E4,0.946E4,0.908E4,0.887E4
!								      
MPDATA,GXZ,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4,1.040E4
MPDATA,GXZ,1,,1.020E4,1.000E4,0.973E4,0.946E4,0.908E4,0.887E4


TB,CREEP,1,,,2                           ! CREEP TABLE
TBDATA,1,5.911E-34,6.25,-0.25


TB,HILL,1,5                              ! HILL TABLE
TBTEMP,750.0
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,800.0				 
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,850.0				 
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,900.0				 
TBDATA,1,1.0,1.0,1.0,1.00,1.00,1.00
TBTEMP,950.0				 
TBDATA,1,1.0,1.0,1.0,1.00,1.00,1.00

For information about the HILL option, see Hill Yield Criterion.

For information about the CREEP option, see Creep Model and Creep Option.

 

8.2.34. ANISO and CREEP Example

Combining anisotropic viscoplasticity with explicit creep:

MP,EX,1,30E6
MP,NUXY,1,0.3

TB,ANISO,1
TBDATA,1,3E4,3E4,3E4,3E6,3E6,3E6
TBDATA,7,3E4,3E4,3E4,3E6,3E6,3E6
TBDATA,13,17321,17321,17321
TBDATA,16,1013514,1013514,1013514

TB,CREEP,1,,,0
TBDATA,1,1e-18,2,1,0,,1

For information about the ANISO option, see Generalized Hill Anisotropy.

For information about the CREEP option, see Creep Model and Creep Option.

 

8.2.35. HILL and CREEP and PLASTIC (BISO) Example

Combining anisotropic implicit creep with bilinear isotropic hardening plasticity:

MPTEMP,1,20,200,400,550,600,650            ! ELASTIC CONSTANTS
MPTEMP,,700,750,800,850,900,950
!
MPDATA,EX,1,1,1.250E4,1.210E4,1.140E4,1.090E4,1.070E4,1.050E4 
MPDATA,EX,1,,1.020E4,0.995E4,0.963E4,0.932E4,0.890E4,0.865E4 
!
MPDATA,EY,1,1,1.250E4,1.210E4,1.140E4,1.090E4,1.070E4,1.050E4 
MPDATA,EY,1,,1.020E4,0.995E4,0.963E4,0.932E4,0.890E4,0.865E4 
!
MPDATA,EZ,1,1,1.250E4,1.210E4,1.140E4,1.090E4,1.070E4,1.050E4
MPDATA,EZ,1,,1.020E4,0.995E4,0.963E4,0.932E4,0.890E4,0.865E4
!
MPDATA,PRXY,1,1,0.351,0.359,0.368,0.375,0.377,0.380 
MPDATA,PRXY,1,,0.382,0.384,0.386,0.389,0.391,0.393
!
MPDATA,PRYZ,1,1,0.351,0.359,0.368,0.375,0.377,0.380 
MPDATA,PRYZ,1,,0.382,0.384,0.386,0.389,0.391,0.393
!
MPDATA,PRXZ,1,1,0.351,0.359,0.368,0.375,0.377,0.380
MPDATA,PRXZ,1,,0.382,0.384,0.386,0.389,0.391,0.393
!
MPDATA,GXY,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4,1.040E4
MPDATA,GXY,1,,1.020E4,1.000E4,0.973E4,0.946E4,0.908E4,0.887E4
!	      
MPDATA,GYZ,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4,1.040E4
MPDATA,GYZ,1,,1.020E4,1.000E4,0.973E4,0.946E4,0.908E4,0.887E4
!								      
MPDATA,GXZ,1,1,1.190E4,1.160E4,1.110E4,1.080E4,1.060E4,1.040E4
MPDATA,GXZ,1,,1.020E4,1.000E4,0.973E4,0.946E4,0.908E4,0.887E4


TB,PLAS,1,,,BISO                         ! PLASTIC (BISO) TABLE
TBDATA,1,180,200


TB,CREEP,1,,,2                           ! CREEP TABLE
TBDATA,1,5.911E-34,6.25,-0.25


TB,HILL,1,5                              ! HILL TABLE
TBTEMP,750.0
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,800.0				 
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,850.0				 
TBDATA,1,1.0,1.0,1.0,0.93,0.93,0.93
TBTEMP,900.0				 
TBDATA,1,1.0,1.0,1.0,1.00,1.00,1.00
TBTEMP,950.0				 
TBDATA,1,1.0,1.0,1.0,1.00,1.00,1.00

For information about the HILL option, see Hill Yield Criterion.

For information about the CREEP option, see Creep Model and Creep Option.

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

 

8.2.36. HILL and CREEP and PLASTIC (MISO) Example

Combining multilinear isotropic hardening with Hill anisotropic plasticity and implicit creep:

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,PLAS,1,,7,MISO                ! MISO TABLE
TBPT,,0.00000,30000
TBPT,,4.00E-3,32000
TBPT,,8.10E-3,33800
TBPT,,1.25E-2,35000
TBPT,,2.18E-2,36500
TBPT,,3.10E-2,38000
TBPT,,4.05E-2,39000

TB,HILL,1                        ! HILL TABLE
TBDATA,1,1.0,1.1,0.9,0.85,0.9,0.80

TB,CREEP,1,,,2                   ! CREEP TABLE
TBDATA,1,1.5625E-14,5.0,-0.5,0.0

For information about the HILL option, see Hill Yield Criterion.

For information about the CREEP option, see Creep Model and Creep Option.

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

 

8.2.37. HILL and CREEP and NLISO Example

Combining anisotropic implicit creep with nonlinear isotropic hardening plasticity:

MP,EX,1,20.0E5                   ! ELASTIC CONSTANTS
MP,NUXY,1,0.3

TB,NLISO,1                       ! NLISO TABLE
TBDATA,1,30000,100000,5200,172

TB,HILL,1                        ! HILL TABLE
TBDATA,1,1.0,1.1,0.9,0.85,0.9,0.80

TB,CREEP,1,,,2                   ! CREEP TABLE
TBDATA,1,1.5625E-14,5.0,-0.5,0.0

For information about the HILL option, see Hill Yield Criterion.

For information about the CREEP option, see Creep Model and Creep Option.

For information about the NLISO option, see Nonlinear Isotropic Hardening.

 

8.2.38. HILL and CREEP and PLASTIC (BKIN) Example

Combining anisotropic implicit creep with bilinear kinematic hardening plasticity:

MP,EX,1,1e7                     ! ELASTIC CONSTANTS
MP,NUXY,1,0.32 

TB,PLAS,1,,,BKIN                ! PLASTIC (BKIN) TABLE
TBDATA,1,42000,1000   

TB,CREEP,1,,,6                  ! CREEP TABLES
TBDATA,1,7.4e-21,3.5,0,0,0,0  

TB,HILL,1                       ! HILL TABLE
TBDATA,1,1.15,1.05,1.0,1.0,1.0,1.0

For information about the HILL option, see Hill Yield Criterion.

For information about the CREEP option, see Creep Model and Creep Option.

For information about the PLASTIC (BKIN) option, see Bilinear Kinematic Hardening.

 

8.2.39. HYPER and PRONY Example

Combining implicit hyperelasticity and viscoelasticity: 

c10=293   
c01=177   
TB,HYPER,1,,,MOON             !!!! type 1 is Mooney-Rivlin
TBDATA,1,c10,c01
a1=0.1
a2=0.2
a3=0.3
t1=10
t2=100
t3=1000
tb,prony,1,,3,shear          ! define Prony constants
tbdata,1,a1,t1,a2,t2,a3,t3

For information about hyperelasticity, see Hyperelasticity and Hyperelasticity Model.

For information about the viscoelasticity, see Viscoelasticity and Viscoelasticity Model.

 

8.2.40. AHYPER and PRONY Example

Combining anisotropic hyperelasticity and viscoelasticity: 

! defining material constants for anisotropic hyperelastic option with TB,AHYPER command
tb,ahyper,1,1,31,poly
! a1,a2,a3
tbdata,1,10,2,0.1
! b1,b2,b3
tbdata,4,5,1,0.1
! c2,c3,c4,c5,c6
tbdata,7,1,0.02,0.002,0.001,0.0005
! d2,d3,d4,d5,d6
tbdata,12,1,0.02,0.002,0.001,0.0005
! e2,e3,e4,e5,e6
tbdata,17,1,0.02,0.002,0.001,0.0005
! f2,f3,f4,f5,f6
tbdata,22,1,0.02,0.002,0.001,0.0005
! g2,g3,g4,g5,g6
tbdata,27,1,0.02,0.002,0.001,0.0005

!compressibility parameter d
tb,ahyper,1,1,1,pvol
tbdata,1,1e-3

!orientation vector A=A(x,y,z)
tb,ahyper,1,1,3,avec
tbdata,1,1,0,0
!orientation vector B=B(x,y,z)
tb,ahyper,1,1,3,bvec
tbdata,1,1/sqrt(2),1/sqrt(2),0

! defining material constants for Prony series with TB,PRONY command
      a1=0.1
      a2=0.2
      a3=0.3
      t1=10
      t2=100
      t3=1000
      tb,prony,1,,3,shear          ! define Prony constants
      tbdata,1,a1,t1,a2,t2,a3,t3

For information about anisotropic hyperelasticity, see Anisotropic Hyperelasticity (TB,AHYPER) and Anisotropic Hyperelasticity Model.

Viscoelastic behavior is assumed to be isotropic. For information about the viscoelasticity, see Viscoelasticity and Viscoelasticity Model.

 

8.2.41. HYPER and PRONY and CDM Example

Combining hyperelasticity, viscoelasticity and Mullins effect:

! Ogden hyperelastic potential
a1 = 2.7971
m1 = 0.77817
a2 = -2.7188
m2 = -0.011229
a3 = 10.505
m3 = 1.269e-7
a4 = 0.33382
m4 = 16.169

tb,hyper,1,,4,ogden
tbdata,1,m1,a1
tbdata,3,m2,a2
tbdata,5,m3,a3
tbdata,7,m4,a4
tbdata,9,0.0,0.,0.,0.

! Viscoelasticity
TB,PRONY,1,1,4,SHEAR
TBDATA,1, 0.1 , 1E+2
TBDATA,3, 0.1 , 1E+1
TBDATA,5, 0.1 , 1E+0
TBDATA,7, 0.1 , 1E-1

! Ogden Roxburgh Mullins effect
TB,CDM,1,,3,PSE2         
TBDATA,1,1.1,50.0,0.2
 

8.2.42. HYPER and CDM Example

Combining hyperelasticity and Mullins effect damage:
! Ogden hyperelastic potential
a1 = 2.7971
m1 = 0.77817
a2 = -2.7188
m2 = -0.011229
a3 = 10.505
m3 = 1.269e-7
a4 = 0.33382
m4 = 16.169

tb,hyper,1,,4,ogden
tbdata,1,m1,a1
tbdata,3,m2,a2
tbdata,5,m3,a3
tbdata,7,m4,a4
tbdata,9,0.0,0.,0.,0.

! Ogden Roxburgh Mullins effect
TB,CDM,1,,3,PSE2         
TBDATA,1,1.1,50.0,0.2
 

8.2.43. HYPER with Embedded Fibers Example

! ISOTROPIC MATRIX
TB,HYPER,1,,,NEO
TBDATA,1,1E2,1E-6

! TWO EMBEDDED FIBER DIRECTIONS
TB,HYPER,1,1,,EXF1
TBDATA,1,1,0,0
TBDATA,4,1/SQRT(2),1/SQRT(2),0

! FIBER TENSION PROPERTIES
TB,HYPER,1,,,EX1
TBDATA,1,1000,2         ! FIBER 1
TBDATA,3,2000,2         ! FIBER 2

! FIBER COMPRESSION PROPERTIES
TB,HYPER,1,,,EXA1
TBDATA,1,10,2         ! FIBER 1
TBDATA,3,0,0          ! FIBER 2, TENSION ONLY
 

8.2.44. HYPER and CDM with Embedded Fibers Example

! MOONEY RIVLIN MATRIX MATERIAL
TB,HYPER,1,1,2,MOONEY
TBDATA,1,1E2,120,1E-5 

! TWO EMBEDDED FIBER DIRECTIONS   
TB,HYPER,1,1,,EXF1 
TBDATA,1,1,0,0          ! FIBER 1  
TBDATA,4,0,1,0          ! FIBER 2

! FIBER TENSION PROPERTIES  
TB,HYPER,1,1,,EX1
TBDATA,1,8,2         ! FIBER 1  
TBDATA,3,4,2         ! FIBER 2  

! FIBER COMPRESSION PROPERTIES  
TB,HYPER,1,1,,EXA1  
TBDATA,1,0.1,0         ! FIBER 1 
TBDATA,3,0.1,0         ! FIBER 2

! OGDEN ROXBURGH MULLINS EFFECT
TB,CDM,1,,3,PSE2         
TBDATA,1,1.1,50.0,0.2
 

8.2.45. HYPER, CDM, and PRONY with Embedded Fibers Example

! MOONEY RIVLIN MATRIX MATERIAL
TB,HYPER,1,1,2,MOONEY
TBDATA,1,1E2,120,1E-5 

! TWO EMBEDDED FIBER DIRECTIONS   
TB,HYPER,1,1,,EXF1 
TBDATA,1,1,0,0          ! FIBER 1  
TBDATA,4,0,1,0          ! FIBER 2

! FIBER TENSION PROPERTIES  
TB,HYPER,1,1,,EX1
TBDATA,1,8,2         ! FIBER 1  
TBDATA,3,4,2         ! FIBER 2  

! FIBER COMPRESSION PROPERTIES  
TB,HYPER,1,1,,EXA1  
TBDATA,1,0.1,0         ! FIBER 1 
TBDATA,3,0.1,0         ! FIBER 2

! VISCOELASTICITY
TB,PRONY,1,1,4,SHEAR
TBDATA,1, 0.1 , 1E+2
TBDATA,3, 0.1 , 1E+1
TBDATA,5, 0.1 , 1E+0
TBDATA,7, 0.1 , 1E-1

! OGDEN ROXBURGH MULLINS EFFECT
TB,CDM,1,,3,PSE2         
TBDATA,1,1.1,50.0,0.2
 

8.2.46. EDP and CREEP and PLASTIC (MISO) Example

Combining Extended Drucker-Prager with implicit creep and multilinear hardening:

ys=100.0
alpha=0.1
!
!define edp for material 1
!
tb,edp,1,,,LYFUN
tbdata,1,alpha,ys
tb,edp,1,,,LFPOT
tbdata,1,alpha
!
!define miso hardening for material 1
!
tb,plastic,1,,2,miso
tbpt,defi,0.0,ys
tbpt,defi,1,1000+ys
!
!define implicit creep for material 1
!
tb,creep,1,,4,1
tbdata,1,1.0e-2,0.5,0.5,0.0

/solu
KBC,0
nlgeom,on
cnvtol,F,1.0,1.0e-10
rate,on
outres,all,all
time,5 
nsub,100,1000,10
solv

For information about the EDP option, see:

For information about the CREEP option, see Creep Model and Creep Option.

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

 

8.2.47. CAP and CREEP and PLASTIC (MISO) Example

Combining geomaterial cap with implicit creep and multilinear hardening: 

TB,EDP,1,,11,CYFUN
tbdata, 	1,     1.0
tbdata, 	2,     1.0
tbdata, 	3,     -80
tbdata, 	4,     10
tbdata, 	5,    0.001
tbdata, 	6,     2
tbdata, 	7,     0.05
tbdata, 	8,     1.0
tbdata, 	9,     0.6
tbdata,        10,     3.0/1000
tbdata,        11,     0.0
tb,plastic,1,,2,miso
tbpt,defi,0.0,8.0
tbpt,defi,1.0,100.0
tb,creep,1,,4,1
tbeo,capc,shea
tbdata,1,1.0e-4,0.6,0.4,0.0
tb,creep,1,,4,1
tbeo,capc,comp
tbdata,1,2.0e-2,0.5,0.5,0.0nlgeom,on
cnvtol,F,1.0,1.0e-10
rate,on
outres,all,all
time,5 
nsub,100,1000,10
solv

For information about the cap model, see:

For information about the CREEP option, see Creep Model and Creep Option.

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

 

8.2.48. CHABOCHE and CREEP Example

Combining Chaboche nonlinear kinematic hardening with implicit creep:

YOUNGS= 30e3
NU    =0.3
SIGMA0= 18.0

MP,EX,1,YOUNGS
MP,PRXY,1,NU

TB,CHAB,1,1,3
TBDATA,1,SIGMA0,
TBDATA,2,5174000,4607500,17155,1040,895.18,9

TB,CREEP,1,1, ,1
TBDATA,1, 2.0E-10,0.01,0.1,0

Combining Chaboche nonlinear kinematic hardening with explicit creep:

MP,EX,1,200000
MP,NUXY,1,0.3

TB,CHAB,1,,
TBDATA,1,980,22400,0

TB,CREEP,1,,,0
TBDATA,1,3.125e-20,5,1,0,,1

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

For information about the CREEP option, see Creep Model and Creep Option.

 

8.2.49. CHABOCHE and CREEP and NLISO Example

Combining Chaboche nonlinear kinematic hardening with implicit creep and nonlinear power law isotropic hardening:

YOUNGS= 30e3
NU    =0.3
SIGMA0= 18.0

MP,EX,1,YOUNGS
MP,PRXY,1,NU

TB,CHAB,1,1,3
TBDATA,1,SIGMA0,
TBDATA,2,5174000,4607500,17155,1040,895.18,9

TB,CREEP,1,1, ,1
TBDATA,1, 2.0E-10,0.01,0.1,0

TB,NLISO,1,,,POWER
TBDATA,1,SIGMA0
TBDATA,2,0.5

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

For information about the CREEP option, see Creep Model and Creep Option.

For information about the NLISO option, see Nonlinear Isotropic Hardening.

 

8.2.50. CHABOCHE and CREEP and HILL Example

Combining anisotropic plasticity with Chaboche nonlinear kinematic hardening and implicit creep:

YOUNGS= 30e3
NU    =0.3
SIGMA0= 18.0

MP,EX,1,YOUNGS
MP,PRXY,1,NU

TB,CHAB,1,1,3
TBDATA,1,SIGMA0,
TBDATA,2,5174000,4607500,17155,1040,895.18,9

TB,CREEP,1,1, ,1
TBDATA,1, 2.0E-10,0.01,0.1,0

RXX=1.1
RYY=1.1
RZZ=1.0
RXY=1.0
RYZ=1.0
RZX=1.0

TB,HILL,1,1
TBDATA,1,RXX,RYY,RZZ,RXY,RYZ,RZX

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

For information about the CREEP option, see Creep Model and Creep Option.

For information about the HILL option, see Hill Yield Criterion.

 

8.2.51. CHABOCHE and CREEP and HILL

Combining anisotropic plasticity with Chaboche nonlinear kinematic hardening and implicit creep:

YOUNGS= 30e3
NU    =0.3
SIGMA0= 18.0

MP,EX,1,YOUNGS
MP,PRXY,1,NU

TB,CHAB,1,1,3
TBDATA,1,SIGMA0,
TBDATA,2,5174000,4607500,17155,1040,895.18,9

TB,CREEP,1,1, ,1
TBDATA,1, 2.0E-10,0.01,0.1,0

RXX=1.1
RYY=1.1
RZZ=1.0
RXY=1.0
RYZ=1.0
RZX=1.0

TB,HILL,1,1
TBDATA,1,RXX,RYY,RZZ,RXY,RYZ,RZX

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

For information about the CREEP option, see Creep Model and Creep Option.

For information about the HILL option, see Hill Yield Criterion.

 

8.2.52. CREEP and RATE and CHABOCHE and PLASTIC (MISO)

Combining multilinear isotropic hardening plasticity, rate-dependent plasticity (viscoplasticity), creep, and Chaboche nonlinear kinematic hardening plasticity:

MP,EX,1,79650e6
MP,NUXY,1,0.33

TB,CREEP,1,1, ,2                ! DEFINE CREEP MATERIAL DATA
TBDATA,1, 5.0E-8,0.4,0.1,0

TB,CHAB,1,,2                    ! DEFINE CHABOCHE MATERIAL DATA
TBDATA,1,1.5e8
TBDATA,2,62511e7,2000
TBDATA,4,62511e6,1000

TB,RATE,1,,,PERZYNA             ! RATE TABLE
TBDATA,1,0.1,2
 
TB,PLASTIC,1, , ,MISO           ! Activate TB,PLASTIC data table
TBPT,DEFI,0,1.5e8     
TBPT,DEFI,0.01,2.8e8
TBPT,DEFI,0.05,4e8
TBPT,DEFI,0.1,4.1e8

For information about the CREEP option, see Creep Model and Creep Option.

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

 

8.2.53. CREEP and RATE and CHABOCHE and PLASTIC (MISO) and PLASTIC (KSR2) Example

Combining multilinear isotropic hardening plasticity, rate-dependent plasticity (viscoplasticity), creep, and Chaboche nonlinear kinematic hardening plasticity with kinematic static recovery:

MP,EX,1,79650e6
MP,NUXY,1,0.33

TB,CREEP,1,1, ,2                ! DEFINE CREEP MATERIAL DATA
TBDATA,1, 5.0E-8,0.4,0.1,0

TB,CHAB,1,,2,TRATE              ! DEFINE CHABOCHE MATERIAL DATA
TBDATA,1,1.5e8
TBDATA,2,62511e7,2000
TBDATA,4,62511e6,1000

TB,RATE,1,,,PERZYNA             ! RATE TABLE
TBDATA,1,0.1,2
 
TB,PLASTIC,1, , ,MISO           ! Activate TB,PLASTIC data table
TBPT,DEFI,0,1.5e8     
TBPT,DEFI,0.01,2.8e8
TBPT,DEFI,0.05,4e8
TBPT,DEFI,0.1,4.1e8

TB,PLASTIC,1,,2,KSR2            ! Kinematic hardening static recovery
TBDATA,1,1e-12,2.1,
TBDATA,3,1e-12,2.15,

For information about the CREEP option, see Creep Model and Creep Option.

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

 

8.2.54. CAST and CHABOCHE Example

Combining cast iron with Chaboche nonlinear kinematic hardening:

/prep7
mp,  ex, 1,14.773E6
mp,nuxy, 1,0.2273

! Define cast iron model
TB,CAST,1,,,ISOTROPIC
TBDATA,1,0.04

TB,CAST,1,,,ROUNDING
TBDATA,1,0.1

TB,CAST,1,1,5,TENSION
TBPT,,0.000E-00,0.813E+04
TBPT,, 1.13E-04,0.131E+05
TBPT,, 8.69E-04,0.241E+05
TBPT,, 1.55E-03,0.288E+05
TBPT,, 2.32E-03,0.322E+05

TB,CAST,1,1,5,COMPRESSION
TBPT,,0.000E-00,0.300E+05
TBPT,, 1.62E-03,0.500E+05
TBPT,, 4.07E-03,0.581E+05
TBPT,, 6.56E-03,0.656E+05
TBPT,, 9.26E-03,0.700E+05

TB,CHABOCHE,1,1,3,TRATE
TBDATA,1,0,4.0783E3,0.8165
TBDATA,4,3.8973E4,743.0119,3.8973E5,7430.119

For information about the CAST option, see Cast Iron.

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

 

8.2.55. EDP and PELAS and PLASTIC (MISO) Example

Combining multilinear isotropic hardening, extended Drucker-Prager plasticity, and porous elasticity:

/PREP7
! poro-elasticity parameters
YoungsModulus=3000             ! initial Young's modulus
Kappa=0.24e-2                  ! swell index
NU0=0.2                        ! Poisson's ratio
E0=0.34                        ! initial void ratio
k0=YoungsModulus/(3*(1-2*nu0)) ! initial bulk modulus
pt_el=k0*Kappa/(1+E0)          ! elastic tensile strength limit

! strength parameters of the Drucker-Prager model
alpha=0.857     ! EDP yield function parameter – pressure
sig_y=2.424     ! EDP yield function parameter – strength
alpha_bar=0.772 ! EDP plastic potential parameter – pressure

! porous elastic model
TB,PELAS,1,,,POISSON
TBDATA,1,kappa,pt_el,nu0,e0

! linear yield function
TB,EDP,1,,,LYFUN
TBDATA,1,alpha

! multilinear isotropic hardening
TB,PLAS,1,,2,MISO
TBPT,DEFI,0,sig_y
TBPT,DEFI,1E-003,1.05*sig_y

! linear plastic potential
TB,EDP,1,,,LFPOT
TBDATA,1,alpha_bar

For information about the EDP option, see Extended Drucker-Prager (EDP).

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

For information about the PELAS option, see Porous Elasticity.

 

8.2.56. HYPER and PLASTIC (BISO) Example

Combining isotropic hyperelasticity and von Mises plasticity with isotropic bilinear hardening:

! material parameters
YoungsModulus=20E3
PoissonsRatio=0.2
BulkModulus=YoungsModulus/(3.0*(1.0-2.0*PoissonsRatio))
ShearModulus=YoungsModulus/(2.0*(1.0+PoissonsRatio))
YieldStress=1.5
PlasticTangentModulus=YoungsModulus/5000

! isotropic hyperelastic model (Neo-Hookean)
TB,HYPER,1,,,NEO
TBDATA,1,ShearModulus,2.0/BulkModulus

! von Mises plasticity with bilinear isotropic hardening
TB,PLASTIC,1,,,BISO
TBDATA,1,YieldStress,PlasticTangentModulus

For information about the finite-strain plasticity model, see Finite-Strain Plasticity.

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

For information about the HYPER option, see Hyperelasticity.

 

8.2.57. PLASTIC (KINH) and CDM (GDMG) Example

Combining von Mises plasticity with multilinear kinematic hardening, and generalized damage:

! isotropic elasticity
MP,EX,1,200000
MP,NUXY,1,0.3

! von Mises plasticity with multilinear kinematic hardening
sig0=60.0
TB,PLASTIC,1,,,KINH
TBPT,,0.0,sig0
TBPT,,0.1,sig0*1.2
TBPT,,0.2,sig0*1.3
TBPT,,2.0,sig0*1.3

! generalized damage
! damage profile
c_eta=10
eta_cr1=0
zz=0.6
! fatigue damage
m0=2.8
n0=0.8
p0=0.6E-2
a0=50
TB,CDM,1,,,GDMG
TBDATA,1,eta_cr1,zz,m0,n0,p0,a0
TBDATA,7,c_eta

For information about the PLASTIC (KINH) option, see Multilinear Kinematic Hardening.

For information about the CDM (GDMG) option, see Regularized Generalized Damage for Fatigue and Thermomechanical Fatigue.

 

8.2.58. CREEP and CDM (GDMG) Example

Combining von Norton creep model and generalized damage:

! isotropic elasticity
MP,EX,1,200E+03
MP,NUXY,1,0.3

! Norton creep
TB,CREEP,1,,,10,
TBDATA,1,1.5625E-14,5.0,100.0

! generalized damage
! damage profile
c_eta=10
eta_cr1=1E-5
zz=2
! fatigue damage
m0=1.0
n0=1.0
p0=0.6
a0=600
TB,CDM,1,,,GDMG
TBDATA,1,eta_cr1,zz,m0,n0,p0,a0
TBDATA,7,c_eta

For information about the CREEP option, see Creep Model and Creep Option.

For information about the CDM (GDMG) option, see Regularized Generalized Damage for Fatigue and Thermomechanical Fatigue.

 

8.2.59. PLASTIC (MISO) and CREEP and CDM (GDMG) Example

Combining von Mises plasticity with isotropic multilinear hardening, creep, and generalized damage:

! isotropic elasticity
MP,EX,1,200E+03
MP,NUXY,1,0.3

! von Mises plasticity with multilinear isotropic hardening
TB,PLASTIC,1,,,MISO
TBPT,DEFI,0,100
TBPT,DEFI,0.01,475
TBPT,DEFI,0.1,3850

! time hardening creep
TB,CREEP,1,,,2,
TBDATA,1,1.5625e-14,5.0,0.0,100.0

! generalized damage
! damage profile
c_eta=10
eta_cr1=0
zz=2
! fatigue damage
m0=1.0
n0=1
p0=0.6
a0=600
TB,CDM,1,,,GDMG
TBDATA,1,eta_cr1,zz,m0,n0,p0,a0
TBDATA,7,c_eta

For information about the PLASTIC (MISO) option, see Multilinear Isotropic Hardening.

For information about the CREEP option, see Creep Model and Creep Option.

For information about the CDM (GDMG) option, see Regularized Generalized Damage for Fatigue and Thermomechanical Fatigue.

 

8.2.60. CHABOCHE and RATE (EVH) and CDM (GDMG) Example

Combining von Mises plasticity with Chaboche nonlinear kinematic hardening, exponential viscohardening, and generalized damage:

! isotropic linear elasticity
Em=120e9
nu=0.3
TB,ELAS,1
TBDATA,1,Em,nu

! von Mises plasticity
! Chaboche nonlinear kinematic hardening
s0=6e6
gam1=2e2
gam2=1e2
Ck1=5e6
Ck2=10e6
TB,CHABOCHE,1,,2
TBDATA,1,s0
TBDATA,2,Ck1,gam1
TBDATA,4,Ck2,gam2

! rate dependent plasticity
! exponential visco-hardening (EVH)
Kr=100e6
nn=1.5
br=10
Qr=8e6
TB,RATE,1,,,EVH
TBDATA,1,s0,0,Qr,br,1/nn,Kr

! generalized damage
! damage profile
c_eta=10e-6
eta_cr1=0
zz=4
! fatigue damage
m0=1.5
n0=0.011
p0=0.5
a0=500e6
TB,CDM,1,,,GDMG
TBDATA,1,eta_cr1,zz,m0,n0,p0,a0
TBDATA,7,c_eta

For information about the CHABOCHE option, see Nonlinear Kinematic Hardening.

For information about the RATE option, see Rate-Dependent Plasticity (Viscoplasticity) and Viscoplasticity Model.

For information about the CDM (GDMG) option, see Regularized Generalized Damage for Fatigue and Thermomechanical Fatigue.

 

8.2.61. EDP and CDM (GDMG) Example

Combining von extended Drucker-Prager model and generalized damage:

! linear elasticity
Em=2.35e5
nu=0.3
TB,ELAS,1
TBDATA,1,Em,nu

! linear Drucker-Prager model
TB,EDP,1,,,LYFUN
TBDATA,1,1/3,100

! generalized damage
! damage profile
c_eta=9
eta_cr1=0.001
zz=2
! fatigue damage
m0=1.0
n0=1.0
p0=0.6
a0=110
TB,CDM,1,,,GDMG
TBDATA,1,eta_cr1,zz,m0,n0,p0,a0
TBDATA,7,c_eta

For information about the EDP option, see Extended Drucker-Prager (EDP).

For information about the CDM (GDMG) option, see Regularized Generalized Damage for Fatigue and Thermomechanical Fatigue.

 

8.2.62. CAST and CDM (GDMG) Example

Combining von Cast Iron model and generalized damage:

! linear elasticity
MP,EX,1,14.7E5
MP,NUXY,1,0.23

! cast iron model
TB,CAST,1,,,ISOTROPIC
TBDATA,1,0.04
! tensile hardening
TB,CAST,1,1,5,TENSION
TBPT,,0.00E-00,80
TBPT,,0.11E-03,130
TBPT,,0.87E-03,240
TBPT,,1.55E-03,290
TBPT,,2.32E-03,320
! compressive hardening
TB,CAST,1,1,5,COMPRESSION
TBPT,,0.00E-00,300
TBPT,,1.62E-03,500
TBPT,,4.07E-03,581
TBPT,,6.56E-03,656
TBPT,,9.26E-03,700

! generalized damage
! damage profile
c_eta=100
eta_cr1=0.0
zz=0.6
! fatigue damage
m0=2.8
n0=0.8
p0=0.6e-2
a0=500
TB,CDM,1,,,GDMG
TBDATA,1,eta_cr1,zz,m0,n0,p0,a0
TBDATA,7,c_eta

For information about the CAST option, see Cast Iron.

For information about the CDM (GDMG) option, see Regularized Generalized Damage for Fatigue and Thermomechanical Fatigue.

 

8.2.63. PLASTIC (BISO) and CDM (DUCTILE and EXPDMG) Example

Combining von Mises plasticity with bilinear isotropic hardening and ductile damage with exponential damage evolution:

! isotropic linear elasticity
TB,ELASTIC,1,,,ISOT
TBDATA,1,70E3 ! Young’s modulus [MPa]
TBDATE,2,0.33 ! Poisson’s ratio [-]

! von Mises plasticity model with bilinear isotropic hardening
TB,PLASTIC,1,,,BISO
TBDATA,1,350  ! initial yield strength [MPa]
TBDATA,2,100  ! plastic tangent modulus [MPa]

! ductile damage criterion 
TB,CDM,1,,,DUCTILE
! stress triaxiality [-], damage initiation threshold [-]
TBPT,DEFI,0.00,1.00
TBPT,DEFI,0.11,0.61
TBPT,DEFI,0.22,0.37
TBPT,DEFI,0.33,0.22
TBPT,DEFI,0.44,0.14
TBPT,DEFI,0.56,0.08
TBPT,DEFI,0.67,0.05
TBPT,DEFI,0.78,0.03
TBPT,DEFI,0.89,0.02
TBPT,DEFI,1.00,0.01

! exponential damage evolution law
TB,CDM,1,,,EXPDMG
TBDATA,1,100 ! initial slope of damage function [1/mm]

For information about the PLASTIC (BISO) option, see Bilinear Isotropic Hardening.

For information about the CDM (DUCTILE and EXPDMG) option, see Ductile Damage.