223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

fixed RGC homogenization in case dt==0 

removed misplaced parameter statement in IO.f90
This commit is contained in:
Franz Roters 2013-02-13 09:36:06 +00:00
parent 0856a7aa32
commit 35346b70b5
2 changed files with 6 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
code/IO.f90
View File

@ -226,7 +226,7 @@ subroutine IO_open_inputFile(myUnit,model)
integer(pInt) :: myStat integer(pInt) :: myStat
character(len=1024) :: path character(len=1024) :: path
character(len=4), parameter :: InputFileExtension2 = '.pes' character(len=4) :: InputFileExtension2 = '.pes'
#ifdef Abaqus #ifdef Abaqus
path = trim(getSolverWorkingDirectoryName())//trim(model)//InputFileExtension2 ! attempt .pes, if it exists: it should be used path = trim(getSolverWorkingDirectoryName())//trim(model)//InputFileExtension2 ! attempt .pes, if it exists: it should be used

5
code/homogenization_RGC.f90
View File

@ -385,6 +385,11 @@ function homogenization_RGC_updateState(&
real(pReal), dimension(:,:), allocatable :: tract,jmatrix,jnverse,smatrix,pmatrix,rmatrix real(pReal), dimension(:,:), allocatable :: tract,jmatrix,jnverse,smatrix,pmatrix,rmatrix
real(pReal), dimension(:), allocatable :: resid,relax,p_relax,p_resid,drelax real(pReal), dimension(:), allocatable :: resid,relax,p_relax,p_resid,drelax
if (dt==0.0_pReal) then ! no need to do anything if time increment is zero (thgis is usually the first call in FEM)
homogenization_RGC_updateState = .true. ! pretend everything is fine
return
endif
!* ------------------------------------------------------------------------------------------------------------- !* -------------------------------------------------------------------------------------------------------------
!*** Initialization of RGC update state calculation !*** Initialization of RGC update state calculation
!* Get the dimension of the cluster (grains and interfaces) !* Get the dimension of the cluster (grains and interfaces)