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

write output information only for active homogenisations. removed unused variables

This commit is contained in:
Pratheek Shanthraj 2015-04-21 14:10:34 +00:00
parent d049eadce7
commit bec17127d0
1 changed files with 41 additions and 47 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
code/homogenization.f90
View File

@ -46,8 +46,6 @@ module homogenization
materialpoint_subFrac, & materialpoint_subFrac, &
materialpoint_subStep, & materialpoint_subStep, &
materialpoint_subdt materialpoint_subdt
integer(pInt), private :: &
homogenization_maxSizeState
logical, dimension(:,:), allocatable, private :: & logical, dimension(:,:), allocatable, private :: &
materialpoint_requested, & materialpoint_requested, &
materialpoint_converged materialpoint_converged
@ -126,18 +124,18 @@ subroutine homogenization_init()
mesh_element, & mesh_element, &
FE_Nips, & FE_Nips, &
FE_geomtype FE_geomtype
use lattice, only: & #ifdef FEM
lattice_referenceTemperature use crystallite, only: &
crystallite_sizePostResults
#else
use constitutive, only: & use constitutive, only: &
constitutive_maxSizePostResults, & constitutive_maxSizePostResults, &
constitutive_damage_maxSizePostResults, & constitutive_damage_maxSizePostResults, &
constitutive_thermal_maxSizePostResults, & constitutive_thermal_maxSizePostResults, &
constitutive_vacancy_maxSizePostResults constitutive_vacancy_maxSizePostResults
use crystallite, only: & use crystallite, only: &
#ifdef FEM
crystallite_sizePostResults, &
#endif
crystallite_maxSizePostResults crystallite_maxSizePostResults
#endif
use material use material
use homogenization_none use homogenization_none
use homogenization_isostrain use homogenization_isostrain
@ -245,6 +243,7 @@ subroutine homogenization_init()
if (worldrank == 0_pInt) then if (worldrank == 0_pInt) then
call IO_write_jobFile(FILEUNIT,'outputHomogenization') call IO_write_jobFile(FILEUNIT,'outputHomogenization')
do p = 1,material_Nhomogenization do p = 1,material_Nhomogenization
if (any(material_homog == p)) then
i = homogenization_typeInstance(p) ! which instance of this homogenization type i = homogenization_typeInstance(p) ! which instance of this homogenization type
knownHomogenization = .true. ! assume valid knownHomogenization = .true. ! assume valid
select case(homogenization_type(p)) ! split per homogenization type select case(homogenization_type(p)) ! split per homogenization type
@ -279,6 +278,7 @@ subroutine homogenization_init()
do e = 1_pInt,field_Noutput(p) do e = 1_pInt,field_Noutput(p)
write(FILEUNIT,'(a,i4)') trim(field_output(e,p))//char(9),field_sizePostResult(e,p) write(FILEUNIT,'(a,i4)') trim(field_output(e,p))//char(9),field_sizePostResult(e,p)
enddo enddo
endif
enddo enddo
close(FILEUNIT) close(FILEUNIT)
endif endif
@ -412,9 +412,7 @@ subroutine materialpoint_stressAndItsTangent(updateJaco,dt)
FEsolving_execIP, & FEsolving_execIP, &
terminallyIll terminallyIll
use mesh, only: & use mesh, only: &
mesh_element, & mesh_element
mesh_NcpElems, &
mesh_maxNips
use material, only: & use material, only: &
plasticState, & plasticState, &
damageState, & damageState, &
@ -745,9 +743,6 @@ subroutine materialpoint_postResults
FEsolving_execElem, & FEsolving_execElem, &
FEsolving_execIP FEsolving_execIP
use mesh, only: & use mesh, only: &
#ifdef FEM
mesh_maxNips, &
#endif
mesh_element mesh_element
use material, only: & use material, only: &
mappingHomogenization, & mappingHomogenization, &
@ -902,7 +897,6 @@ subroutine homogenization_partitionDeformation(ip,el)
HOMOGENIZATION_ISOSTRAIN_ID, & HOMOGENIZATION_ISOSTRAIN_ID, &
HOMOGENIZATION_RGC_ID HOMOGENIZATION_RGC_ID
use crystallite, only: & use crystallite, only: &
crystallite_partionedF0, &
crystallite_partionedF crystallite_partionedF
use homogenization_isostrain, only: & use homogenization_isostrain, only: &
homogenization_isostrain_partitionDeformation homogenization_isostrain_partitionDeformation