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

reordered calculation of current stiffness to enable restarting with different resolution

This commit is contained in:
Martin Diehl 2012-06-22 09:20:23 +00:00
parent 64823fe6ea
commit 0606845f15
1 changed files with 11 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
code/DAMASK_spectral.f90
View File

@ -97,7 +97,6 @@ program DAMASK_spectral
itmax,& itmax,&
itmin, & itmin, &
memory_efficient, & memory_efficient, &
update_gamma, &
divergence_correction, & divergence_correction, &
DAMASK_NumThreadsInt, & DAMASK_NumThreadsInt, &
fftw_planner_flag, & fftw_planner_flag, &
@ -549,13 +548,11 @@ program DAMASK_spectral
0.0_pReal,ielem,1_pInt,sigma,dsde,P_real(i,j,k,1:3,1:3),dPdF) 0.0_pReal,ielem,1_pInt,sigma,dsde,P_real(i,j,k,1:3,1:3),dPdF)
C = C + dPdF C = C + dPdF
enddo; enddo; enddo enddo; enddo; enddo
C_ref = C * wgt C = C * wgt
C_ref = C
call IO_write_jobBinaryFile(777,'C_ref',size(C_ref)) call IO_write_jobBinaryFile(777,'C_ref',size(C_ref))
write (777,rec=1) C_ref write (777,rec=1) C_ref
close(777) close(777)
call IO_write_jobBinaryFile(777,'C',size(C))
write (777,rec=1) C
close(777)
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! restore deformation gradient and stiffness from saved state ! restore deformation gradient and stiffness from saved state
@ -624,9 +621,9 @@ program DAMASK_spectral
write(538) 'increments', bc(1:N_Loadcases)%incs ! one entry per loadcase write(538) 'increments', bc(1:N_Loadcases)%incs ! one entry per loadcase
write(538) 'startingIncrement', restartInc - 1_pInt ! start with writing out the previous inc write(538) 'startingIncrement', restartInc - 1_pInt ! start with writing out the previous inc
write(538) 'eoh' ! end of header write(538) 'eoh' ! end of header
write(538) materialpoint_results(1_pInt:materialpoint_sizeResults,1,1_pInt:Npoints) ! initial (non-deformed or read-in) results
if (debugGeneral) write(6,'(a)') 'Header of result file written out' if (debugGeneral) write(6,'(a)') 'Header of result file written out'
endif endif
write(538) materialpoint_results(1_pInt:materialpoint_sizeResults,1,1_pInt:Npoints) ! initial (non-deformed or read-in) results
flush(538) flush(538)
!################################################################################################## !##################################################################################################
@ -709,7 +706,7 @@ program DAMASK_spectral
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! calculate reduced compliance ! calculate reduced compliance
if(size_reduced > 0_pInt) then ! calculate compliance in case stress BC is applied if(size_reduced > 0_pInt) then ! calculate compliance in case stress BC is applied
C_lastInc = math_rotate_forward3333(C*wgt,bc(loadcase)%rotation) ! calculate stiffness from former inc C_lastInc = math_rotate_forward3333(C,bc(loadcase)%rotation) ! calculate stiffness from former inc
temp99_Real = math_Plain3333to99(C_lastInc) temp99_Real = math_Plain3333to99(C_lastInc)
k = 0_pInt ! build reduced stiffness k = 0_pInt ! build reduced stiffness
do n = 1_pInt,9_pInt do n = 1_pInt,9_pInt
@ -790,8 +787,9 @@ program DAMASK_spectral
enddo; enddo; enddo enddo; enddo; enddo
call debug_info() call debug_info()
! for test of regridding ! for test of regridding
!if(mod(inc-1,bc(loadcase)%restartFrequency) == 0_pInt .and. restartInc/=inc) & ! if( bc(loadcase)%restartFrequency > 0_pInt .and. &
! call quit(-1*(restartInc+1)) ! trigger exit to regrid ! mod(inc-1,bc(loadcase)%restartFrequency) == 0_pInt .and. &
! restartInc/=inc) call quit(-1*(restartInc+1)) ! trigger exit to regrid
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! copy one component of the stress field to to a single FT and check for mismatch ! copy one component of the stress field to to a single FT and check for mismatch
@ -1024,7 +1022,6 @@ program DAMASK_spectral
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! updated deformation gradient ! updated deformation gradient
do k = 1_pInt, res(3); do j = 1_pInt, res(2); do i = 1_pInt, res(1) do k = 1_pInt, res(3); do j = 1_pInt, res(2); do i = 1_pInt, res(1)
F(i,j,k,1:3,1:3) = F(i,j,k,1:3,1:3) - deltaF_real(i,j,k,1:3,1:3)*wgt ! F(x)^(n+1) = F(x)^(n) + correction; *wgt: correcting for missing normalization F(i,j,k,1:3,1:3) = F(i,j,k,1:3,1:3) - deltaF_real(i,j,k,1:3,1:3)*wgt ! F(x)^(n+1) = F(x)^(n) + correction; *wgt: correcting for missing normalization
enddo; enddo; enddo enddo; enddo; enddo
@ -1048,6 +1045,7 @@ program DAMASK_spectral
enddo ! end looping when convergency is achieved enddo ! end looping when convergency is achieved
CPFEM_mode = 1_pInt ! winding forward CPFEM_mode = 1_pInt ! winding forward
C = C * wgt
write(6,'(a)') '' write(6,'(a)') ''
write(6,'(a)') '==================================================================' write(6,'(a)') '=================================================================='
if(err_div > err_div_tol .or. err_stress > err_stress_tol) then if(err_div > err_div_tol .or. err_stress > err_stress_tol) then
@ -1073,13 +1071,11 @@ program DAMASK_spectral
write (777,rec=1) F write (777,rec=1) F
close (777) close (777)
restartInc=totalIncsCounter restartInc=totalIncsCounter
call IO_write_jobBinaryFile(777,'C',size(C))
write (777,rec=1) C
close(777)
endif endif
if (update_gamma) then
write(6,'(a)') 'update C_ref '
C_ref = C*wgt
endif
endif ! end calculation/forwarding endif ! end calculation/forwarding
enddo ! end looping over incs in current loadcase enddo ! end looping over incs in current loadcase
deallocate(c_reduced) deallocate(c_reduced)