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

polished Abaqus Standard, changed inclusion of OMP functionality from "spectral style" (include) to "marc style" (use) to prevent explicit from crashing

This commit is contained in:
Martin Diehl 2013-03-25 18:25:15 +00:00
parent 966ad2826b
commit e6f6c22a30
2 changed files with 199 additions and 158 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

353
code/DAMASK_abaqus_std.f
View File

@ -1,7 +1,7 @@
! Copyright 2011-13 Max-Planck-Institut für Eisenforschung GmbH ! Copyright 2011-13 Max-Planck-Institut für Eisenforschung GmbH
! !
! This file is part of DAMASK, ! This file is part of DAMASK,
! the Düsseldorf Advanced MAterial Simulation Kit. ! the Düsseldorf Advanced Material Simulation Kit.
! !
! DAMASK is free software: you can redistribute it and/or modify ! DAMASK is free software: you can redistribute it and/or modify
! it under the terms of the GNU General Public License as published by ! it under the terms of the GNU General Public License as published by
@ -16,26 +16,18 @@
! You should have received a copy of the GNU General Public License ! You should have received a copy of the GNU General Public License
! along with DAMASK. If not, see <http://www.gnu.org/licenses/>. ! along with DAMASK. If not, see <http://www.gnu.org/licenses/>.
! !
!############################################################## !--------------------------------------------------------------------------------------------------
!* $Id$ ! $Id$
!******************************************************************** !--------------------------------------------------------------------------------------------------
! Material subroutine for Abaqus !> @author Philip Eisenlohr, Max-Planck-Institut für Eisenforschung GmbH
! !> @author Franz Roters, Max-Planck-Institut für Eisenforschung GmbH
! written by P. Eisenlohr, !> @author Koen Janssens, Paul Scherrer Institut
! F. Roters, !> @author Arun Prakash, Fraunhofer IWM
! K. Janssens 2, !> @brief interfaces DAMASK with Abaqus/Standard
! A. Prakash 3 !> @details put the included file abaqus_v6.env in either your home or model directory,
! !> it is a minimum Abaqus environment file containing all changes necessary to use the
! MPI fuer Eisenforschung, Duesseldorf !> DAMASK subroutine (see Abaqus documentation for more information on the use of abaqus_v6.env)
! 2 PSI, Switzerland !--------------------------------------------------------------------------------------------------
! 3 Fraunhofer IWM, Freiburg
!
! REMARK: put the included file abaqus_v6.env in either your home
! or model directory, it is a minimum Abaqus environment file
! containing all changes necessary to use the MPIE subroutine
! (see Abaqus documentation for more information on the use of abaqus_v6.env)
!
!********************************************************************
#ifndef INT #ifndef INT
#define INT 4 #define INT 4
@ -52,51 +44,48 @@
module DAMASK_interface module DAMASK_interface
implicit none implicit none
character(len=4), dimension(2), parameter :: InputFileExtension = ['.pes','.inp'] character(len=4), dimension(2), parameter :: INPUTFILEEXTENSION = ['.pes','.inp']
character(len=4), parameter :: LogFileExtension = '.log' character(len=4), parameter :: LOGFILEEXTENSION = '.log'
contains contains
!-------------------- !--------------------------------------------------------------------------------------------------
!> @brief just reporting
!--------------------------------------------------------------------------------------------------
subroutine DAMASK_interface_init() subroutine DAMASK_interface_init()
!--------------------
write(6,*) write(6,'(/,a)') ' <<<+- DAMASK_abaqus init -+>>>'
write(6,*) '<<<+- DAMASK_abaqus init -+>>>' write(6,'(a)') ' $Id$'
write(6,*) '$Id$'
#include "compilation_info.f90" #include "compilation_info.f90"
write(6,*)
end subroutine DAMASK_interface_init end subroutine DAMASK_interface_init
!--------------------
function getSolverWorkingDirectoryName() !--------------------------------------------------------------------------------------------------
!-------------------- !> @brief using Abaqus/Standard function to get working directory name
use prec !--------------------------------------------------------------------------------------------------
character(1024) function getSolverWorkingDirectoryName()
implicit none implicit none
character(1024) getSolverWorkingDirectoryName integer :: lenOutDir
integer(pInt) LENOUTDIR
getSolverWorkingDirectoryName='' getSolverWorkingDirectoryName=''
CALL GETOUTDIR( getSolverWorkingDirectoryName, LENOUTDIR ) call getoutdir(getSolverWorkingDirectoryName, lenOutDir)
getSolverWorkingDirectoryName=trim(getSolverWorkingDirectoryName)//'/' getSolverWorkingDirectoryName=trim(getSolverWorkingDirectoryName)//'/'
! write(6,*) 'getSolverWorkingDirectoryName', getSolverWorkingDirectoryName
end function getSolverWorkingDirectoryName end function getSolverWorkingDirectoryName
!-------------------- !--------------------------------------------------------------------------------------------------
function getSolverJobName() !> @brief using Abaqus/Standard function to get solver job name
!-------------------- !--------------------------------------------------------------------------------------------------
use prec character(1024) function getSolverJobName()
implicit none implicit none
character(1024) getSolverJobName, JOBNAME integer :: lenJobName
integer(pInt) LENJOBNAME
getSolverJobName='' getSolverJobName=''
CALL GETJOBNAME(getSolverJobName , LENJOBNAME ) call getJobName(getSolverJobName, lenJobName)
! write(6,*) 'getSolverJobName', getSolverJobName
end function getSolverJobName end function getSolverJobName
@ -124,33 +113,39 @@ end module DAMASK_interface
#include "CPFEM.f90" #include "CPFEM.f90"
subroutine UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,& subroutine UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,&
RPL,DDSDDT,DRPLDE,DRPLDT,STRAN,DSTRAN,& RPL,DDSDDT,DRPLDE,DRPLDT,STRAN,DSTRAN,&
TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,NDI,NSHR,NTENS,& TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,NDI,NSHR,NTENS,&
NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,CELENT,& NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,CELENT,&
DFGRD0,DFGRD1,NOEL,NPT,KSLAY,KSPT,KSTEP,KINC) DFGRD0,DFGRD1,NOEL,NPT,KSLAY,KSPT,KSTEP,KINC)
use prec, only: pReal, & use prec, only: &
pInt pReal, &
use numerics, only: numerics_unitlength pInt
use FEsolving, only: cycleCounter, & use numerics, only: &
theInc, & numerics_unitlength
calcMode, & use FEsolving, only: &
lastMode, & cycleCounter, &
theTime, & theInc, &
theDelta, & calcMode, &
lastIncConverged, & lastMode, &
outdatedByNewInc, & theTime, &
outdatedFFN1, & theDelta, &
terminallyIll, & lastIncConverged, &
symmetricSolver outdatedByNewInc, &
use math, only: invnrmMandel outdatedFFN1, &
use debug, only: debug_info, & terminallyIll, &
debug_reset, & symmetricSolver
debug_levelBasic, & use math, only: &
debug_level, & invnrmMandel
debug_abaqus use debug, only: &
use mesh, only: mesh_FEasCP, & debug_info, &
mesh_ipCoordinates debug_reset, &
debug_levelBasic, &
debug_level, &
debug_abaqus
use mesh, only: &
mesh_FEasCP, &
mesh_ipCoordinates
use CPFEM, only: & use CPFEM, only: &
CPFEM_general, & CPFEM_general, &
CPFEM_init_done, & CPFEM_init_done, &
@ -160,122 +155,165 @@ subroutine UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,&
CPFEM_COLLECT, & CPFEM_COLLECT, &
CPFEM_RESTOREJACOBIAN, & CPFEM_RESTOREJACOBIAN, &
CPFEM_BACKUPJACOBIAN CPFEM_BACKUPJACOBIAN
use homogenization, only: &
use homogenization, only: materialpoint_sizeResults, materialpoint_results materialpoint_sizeResults, &
materialpoint_results
implicit none implicit none
CHARACTER(80) CMNAME integer(pInt), intent(in) :: &
integer(pInt) ndi, nshr, ntens, nstatv, nprops, noel, npt,& nDi, & !< Number of direct stress components at this point
kslay, kspt, kstep, kinc nShr, & !< Number of engineering shear stress components at this point
real(pReal) STRESS(NTENS),STATEV(NSTATV),& nTens, & !< Size of the stress or strain component array (NDI + NSHR)
DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),& nStatV, & !< Number of solution-dependent state variables
STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),& nProps, & !< User-defined number of material constants
PROPS(NPROPS),COORDS(3),DROT(3,3),& noEl, & !< element number
DFGRD0(3,3),DFGRD1(3,3) nPt,& !< integration point number
real(pReal) SSE, SPD, SCD, RPL, DRPLDT, DTIME, TEMP,& kSlay, & !< layer number (shell elements etc.)
DTEMP, PNEWDT, CELENT kSpt, & !< section point within the current layer
real(pReal), dimension (3,3) :: pstress ! not used, but needed for call of cpfem_general kStep, & !< step number
real(pReal), dimension (3,3,3,3) :: dPdF ! not used, but needed for call of cpfem_general kInc !< increment number
character(len=80), intent(in) :: &
cmname !< uses-specified material name, left justified
real(pReal), intent(in) :: &
DTIME, &
TEMP, &
DTEMP, &
CELENT
real(pReal), dimension(1), intent(in) :: &
PREDEF, &
DPRED
real(pReal), dimension(2), intent(in) :: &
TIME !< step time/total time at beginning of the current increment
real(pReal), dimension(3), intent(in) :: &
COORDS
real(pReal), dimension(nTens), intent(in) :: &
STRAN, & !< total strains at beginning of the increment
DSTRAN !< strain increments
real(pReal), dimension(nProps), intent(in) :: &
PROPS
real(pReal), dimension(3,3), intent(in) :: &
DROT, & !< rotation increment matrix
DFGRD0, & !< F at beginning of increment
DFGRD1 !< F at end of increment
real(pReal), intent(inout) :: &
PNEWDT, & !< ratio of suggested new time increment
SSE, & !< specific elastic strain engergy
SPD, & !< specific plastic dissipation
SCD, & !< specific creep dissipation
RPL, & !< volumetric heat generation per unit time at the end of the increment
DRPLDT !< varation of RPL with respect to the temperature
real(pReal), dimension(nTens), intent(inout) :: &
STRESS !< stress tensor at the beginning of the increment, needs to be updated
real(pReal), dimension(nStatV), intent(inout) :: &
STATEV !< solution-dependent state variables
real(pReal), dimension(nTens), intent(out) :: &
DDSDDT, &
DRPLDE
real(pReal), dimension(nTens,nTens), intent(out) :: &
DDSDDE !< Jacobian matrix of the constitutive model
! local variables real(pReal) :: temperature ! temp by Abaqus is intent(in)
real(pReal), dimension (3,3) :: pstress ! not used, but needed for call of cpfem_general
real(pReal), dimension (3,3,3,3) :: dPdF ! not used, but needed for call of cpfem_general
real(pReal), dimension(6) :: stress_h real(pReal), dimension(6) :: stress_h
real(pReal), dimension(6,6) :: ddsdde_h real(pReal), dimension(6,6) :: ddsdde_h
integer(pInt) computationMode, i, cp_en integer(pInt) :: computationMode, i, cp_en
logical :: cutBack logical :: cutBack
temperature = temp ! temp is intent(in)
DDSDDT = 0.0_pReal
DRPLDE = 0.0_pReal
if (iand(debug_level(debug_abaqus),debug_levelBasic) /= 0 .and. noel == 1 .and. npt == 1) then if (iand(debug_level(debug_abaqus),debug_levelBasic) /= 0 .and. noel == 1 .and. npt == 1) then
!$OMP CRITICAL (write2out) !$OMP CRITICAL (write2out)
write(6,*) 'el',noel,'ip',npt write(6,*) 'el',noel,'ip',npt
write(6,*) 'got kinc as',kinc write(6,*) 'got kInc as',kInc
write(6,*) 'got dStran',dstran write(6,*) 'got dStran',dStran
call flush(6) flush(6)
!$OMP END CRITICAL (write2out) !$OMP END CRITICAL (write2out)
endif endif
if (.not. CPFEM_init_done) call CPFEM_initAll(temp,noel,npt) if (.not. CPFEM_init_done) call CPFEM_initAll(temperature,noel,npt)
cp_en = mesh_FEasCP('elem',noel) cp_en = mesh_FEasCP('elem',noel)
if (time(2) > theTime .or. kinc /= theInc) then ! reached convergence if (time(2) > theTime .or. kinc /= theInc) then ! reached convergence
terminallyIll = .false. terminallyIll = .false.
cycleCounter = -1 ! first calc step increments this to cycle = 0 cycleCounter = -1 ! first calc step increments this to cycle = 0
if (kinc == 1) then ! >> start of analysis << if (kinc == 1) then ! start of analysis
lastIncConverged = .false. ! no Jacobian backup lastIncConverged = .false. ! no Jacobian backup
outdatedByNewInc = .false. ! no aging of state outdatedByNewInc = .false. ! no aging of state
lastMode = .false. ! pretend last step was collection lastMode = .false. ! pretend last step was collection
calcMode = .false. ! pretend last step was collection calcMode = .false. ! pretend last step was collection
!$OMP CRITICAL (write2out) !$OMP CRITICAL (write2out)
write (6,'(i8,1x,i2,1x,a)') noel,npt,'<< UMAT >> start of analysis..!'; call flush(6) write (6,'(i8,1x,i2,1x,a)') noel,npt,'<< UMAT >> start of analysis..!'; flush(6)
!$OMP END CRITICAL (write2out) !$OMP END CRITICAL (write2out)
else if (kinc - theInc > 1) then ! >> restart of broken analysis << else if (kinc - theInc > 1) then ! restart of broken analysis
lastIncConverged = .false. ! no Jacobian backup lastIncConverged = .false. ! no Jacobian backup
outdatedByNewInc = .false. ! no aging of state outdatedByNewInc = .false. ! no aging of state
lastMode = .true. ! pretend last step was calculation lastMode = .true. ! pretend last step was calculation
calcMode = .true. ! pretend last step was calculation calcMode = .true. ! pretend last step was calculation
!$OMP CRITICAL (write2out) !$OMP CRITICAL (write2out)
write (6,'(i8,1x,i2,1x,a)') noel,npt,'<< UMAT >> restart of analysis..!'; call flush(6) write (6,'(i8,1x,i2,1x,a)') noel,npt,'<< UMAT >> restart of analysis..!'; flush(6)
!$OMP END CRITICAL (write2out) !$OMP END CRITICAL (write2out)
else ! >> just the next inc << else ! just the next inc
lastIncConverged = .true. ! request Jacobian backup lastIncConverged = .true. ! request Jacobian backup
outdatedByNewInc = .true. ! request aging of state outdatedByNewInc = .true. ! request aging of state
lastMode = .true. ! assure last step was calculation lastMode = .true. ! assure last step was calculation
calcMode = .true. ! assure last step was calculation calcMode = .true. ! assure last step was calculation
!$OMP CRITICAL (write2out) !$OMP CRITICAL (write2out)
write (6,'(i8,1x,i2,1x,a)') noel,npt,'<< UMAT >> new increment..!'; call flush(6) write (6,'(i8,1x,i2,1x,a)') noel,npt,'<< UMAT >> new increment..!'; flush(6)
!$OMP END CRITICAL (write2out) !$OMP END CRITICAL (write2out)
endif endif
else if ( dtime < theDelta ) then ! cutBack
else if ( dtime < theDelta ) then ! >> cutBack << cutBack = .true.
cutBack = .true. terminallyIll = .false.
terminallyIll = .false. cycleCounter = -1 ! first calc step increments this to cycle = 0
cycleCounter = -1 ! first calc step increments this to cycle = 0 calcMode = .true. ! pretend last step was calculation
calcMode = .true. ! pretend last step was calculation !$OMP CRITICAL (write2out)
!$OMP CRITICAL (write2out) write(6,'(i8,1x,i2,1x,a)') noel,npt,'<< UMAT >> cutback detected..!'; flush(6)
write(6,'(i8,1x,i2,1x,a)') noel,npt,'<< UMAT >> cutback detected..!'; call flush(6) !$OMP END CRITICAL (write2out)
!$OMP END CRITICAL (write2out) endif ! convergence treatment end
endif ! convergence treatment end calcMode(npt,cp_en) = .not. calcMode(npt,cp_en) ! ping pong (calc <--> collect)
calcMode(npt,cp_en) = .not. calcMode(npt,cp_en) ! ping pong (calc <--> collect) if (calcMode(npt,cp_en)) then ! now calc
if (calcMode(npt,cp_en)) then ! now calc
computationMode = CPFEM_CALCRESULTS computationMode = CPFEM_CALCRESULTS
if ( lastMode .neqv. calcMode(npt,cp_en) ) then ! first after ping pong if ( lastMode .neqv. calcMode(npt,cp_en) ) then ! first after ping pong
call debug_reset() ! resets debugging call debug_reset() ! resets debugging
outdatedFFN1 = .false. outdatedFFN1 = .false.
cycleCounter = cycleCounter + 1 cycleCounter = cycleCounter + 1
endif endif
if(outdatedByNewInc) then if(outdatedByNewInc) then
outdatedByNewInc = .false. outdatedByNewInc = .false.
computationMode = ior(computationMode,CPFEM_AGERESULTS) ! calc and age results computationMode = ior(computationMode,CPFEM_AGERESULTS) ! calc and age results
endif endif
else ! now collect else ! now collect
computationMode = CPFEM_COLLECT computationMode = CPFEM_COLLECT
if(lastMode .neqv. calcMode(npt,cp_en) .and. .not. terminallyIll) then if(lastMode .neqv. calcMode(npt,cp_en) .and. .not. terminallyIll) then
call debug_info() ! first after ping pong reports debugging call debug_info() ! first after ping pong reports debugging
endif endif
if (lastIncConverged) then if (lastIncConverged) then
lastIncConverged = .false. lastIncConverged = .false.
computationMode = ior(computationMode,CPFEM_BACKUPJACOBIAN) ! backup Jacobian after convergence computationMode = ior(computationMode,CPFEM_BACKUPJACOBIAN) ! backup Jacobian after convergence
elseif ( cutBack ) then elseif ( cutBack ) then
cutBack = .false. cutBack = .false.
computationMode = ior(computationMode,CPFEM_RESTOREJACOBIAN) ! restore Jacobian after cutback computationMode = ior(computationMode,CPFEM_RESTOREJACOBIAN) ! restore Jacobian after cutback
endif endif
mesh_ipCoordinates(1:3,npt,cp_en) = numerics_unitlength * COORDS mesh_ipCoordinates(1:3,npt,cp_en) = numerics_unitlength * COORDS
endif endif
theTime = time(2) ! record current starting time theTime = time(2) ! record current starting time
theDelta = dtime ! record current time increment theDelta = dtime ! record current time increment
theInc = kinc ! record current increment number theInc = kinc ! record current increment number
lastMode = calcMode(npt,cp_en) ! record calculationMode lastMode = calcMode(npt,cp_en) ! record calculationMode
if (iand(debug_level(debug_abaqus),debug_levelBasic) /= 0) then if (iand(debug_level(debug_abaqus),debug_levelBasic) /= 0) then
!$OMP CRITICAL (write2out) !$OMP CRITICAL (write2out)
write(6,'(a16,1x,i2,1x,a,i8,a,i8,1x,i5,a)') 'computationMode',computationMode,'(',cp_en,':',noel,npt,')'; call flush(6) write(6,'(a16,1x,i2,1x,a,i8,a,i8,1x,i5,a)') 'computationMode',computationMode,'(',cp_en,':',noel,npt,')'; flush(6)
!$OMP END CRITICAL (write2out) !$OMP END CRITICAL (write2out)
endif endif
call CPFEM_general(computationMode,dfgrd0,dfgrd1,temp,dtime,noel,npt,stress_h,ddsdde_h, pstress, dPdF) call CPFEM_general(computationMode,dfgrd0,dfgrd1,temperature,dtime,noel,npt,stress_h,ddsdde_h, pstress, dPdF)
! Mandel: 11, 22, 33, SQRT(2)*12, SQRT(2)*23, SQRT(2)*13 ! Mandel: 11, 22, 33, SQRT(2)*12, SQRT(2)*23, SQRT(2)*13
! straight: 11, 22, 33, 12, 23, 13 ! straight: 11, 22, 33, 12, 23, 13
@ -300,19 +338,22 @@ subroutine UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,&
statev = materialpoint_results(1:min(nstatv,materialpoint_sizeResults),npt,mesh_FEasCP('elem', noel)) statev = materialpoint_results(1:min(nstatv,materialpoint_sizeResults),npt,mesh_FEasCP('elem', noel))
if ( terminallyIll ) pnewdt = 0.5_pReal ! force cutback directly ? if ( terminallyIll ) pnewdt = 0.5_pReal ! force cutback directly ?
end subroutine UMAT end subroutine UMAT
!********************************************************************
! This subroutine replaces the corresponding Marc subroutine
!********************************************************************
subroutine quit(mpie_error)
use prec, only: pInt !--------------------------------------------------------------------------------------------------
!> @brief calls the exit function of Abaqus/Standard
!--------------------------------------------------------------------------------------------------
subroutine quit(mpie_error)
use prec, only: &
pInt
implicit none implicit none
integer(pInt) mpie_error integer(pInt) mpie_error
flush(6) flush(6)
call xit call xit
end subroutine quit
end subroutine quit

4
code/numerics.f90
View File

@ -139,11 +139,11 @@ subroutine numerics_init
IO_warning, & IO_warning, &
IO_timeStamp IO_timeStamp
#ifndef Marc #ifdef Spectral
!$ use OMP_LIB, only: omp_set_num_threads ! Use the standard conforming module file for omp if not using MSC.Marc !$ use OMP_LIB, only: omp_set_num_threads ! Use the standard conforming module file for omp if not using MSC.Marc
#endif #endif
implicit none implicit none
#ifdef Marc #ifndef Spectral
!$ include "omp_lib.h" ! use the not F90 standard conforming include file to prevent crashes with some versions of MSC.Marc !$ include "omp_lib.h" ! use the not F90 standard conforming include file to prevent crashes with some versions of MSC.Marc
#endif #endif
integer(pInt), parameter :: fileunit = 300_pInt ,& integer(pInt), parameter :: fileunit = 300_pInt ,&