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
|
|
|
|
|
! *** Maximum number of internal cutbacks in time step ***
|
|
|
|
|
integer(pInt), parameter :: ncut=7_pInt
|
|
|
|
|
! *** Perturbation of strain array for numerical calculation of FEM Jacobi matrix ***
|
|
|
|
|
real(pReal), parameter :: pert_e=1.0e-5_pReal
|
|
|
|
|
! *** Maximum number of iterations in outer (statevariables) loop ***
|
|
|
|
|
integer(pInt), parameter :: nouter = 50_pInt
|
|
|
|
|
! *** Convergence criteria for outer (statevariables) loop ***
|
|
|
|
|
real(pReal), parameter :: tol_outer = 1.0e-4_pReal
|
|
|
|
|
! *** Maximum number of iterations in inner (stress) loop ***
|
|
|
|
|
integer(pInt), parameter :: ninner = 2000_pInt
|
|
|
|
|
! *** Convergence criteria for inner (stress) loop ***
|
|
|
|
|
real(pReal), parameter :: tol_inner = 1.0e-3_pReal
|
|
|
|
|
! *** 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
|
