2007-03-20 19:25:22 +05:30
|
|
|
|
|
|
|
!##############################################################
|
|
|
|
MODULE prec
|
|
|
|
!##############################################################
|
|
|
|
|
|
|
|
implicit none
|
2007-03-29 12:45:12 +05:30
|
|
|
! *** Precision of real and integer variables ***
|
2007-03-20 19:25:22 +05:30
|
|
|
integer, parameter :: pReal = 8
|
2007-03-20 21:28:53 +05:30
|
|
|
integer, parameter :: pInt = 4
|
2007-03-29 12:45:12 +05:30
|
|
|
! *** Numerical parameters ***
|
|
|
|
! *** How frequently the jacobian is recalculated ***
|
|
|
|
integer (pInt), parameter :: ijaco = 5_pInt
|
2007-04-10 16:52:06 +05:30
|
|
|
! *** Maximum number of internal cutbacks in time step ***
|
2007-04-11 15:35:08 +05:30
|
|
|
integer(pInt), parameter :: nCutback = 7_pInt
|
2007-04-10 16:52:06 +05:30
|
|
|
! *** Maximum number of regularization attempts for Jacobi inversion ***
|
2007-04-11 15:35:08 +05:30
|
|
|
integer(pInt), parameter :: nReg = 1_pInt
|
2007-03-29 12:45:12 +05:30
|
|
|
! *** Perturbation of strain array for numerical calculation of FEM Jacobi matrix ***
|
|
|
|
real(pReal), parameter :: pert_e=1.0e-5_pReal
|
2007-04-11 15:35:08 +05:30
|
|
|
! *** Maximum number of iterations in outer (state variables) loop ***
|
|
|
|
integer(pInt), parameter :: nState = 50_pInt
|
|
|
|
! *** Convergence criteria for outer (state variables) loop ***
|
|
|
|
real(pReal), parameter :: tol_State = 1.0e-4_pReal
|
2007-03-29 12:45:12 +05:30
|
|
|
! *** Maximum number of iterations in inner (stress) loop ***
|
2007-04-11 15:35:08 +05:30
|
|
|
integer(pInt), parameter :: nStress = 2000_pInt
|
2007-03-29 12:45:12 +05:30
|
|
|
! *** Convergence criteria for inner (stress) loop ***
|
2007-04-11 15:35:08 +05:30
|
|
|
real(pReal), parameter :: tol_Stress = 1.0e-3_pReal
|
2007-03-29 12:45:12 +05:30
|
|
|
! *** Factor for maximum stress correction in inner (stress) loop ***
|
|
|
|
real(pReal), parameter :: crite = 1.0e-1_pReal
|
2007-03-20 19:25:22 +05:30
|
|
|
|
|
|
|
END MODULE prec
|