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

minor changes to multiphysicsOut

This commit is contained in:
Pratheek Shanthraj 2014-08-10 11:14:43 +00:00
parent 0a01d5f3a1
commit 8f444b0695
2 changed files with 19 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
code/crystallite.f90
View File

@ -3945,11 +3945,17 @@ function crystallite_postResults(ipc, ip, el)
ip, & !< integration point index
ipc !< grain index
real(pReal), dimension(1+crystallite_sizePostResults(microstructure_crystallite(mesh_element(4,el)))+ &
1+plasticState(material_phase(ipc,ip,el))%sizePostResults) :: &! + &
! 1+damageState(material_phase(ipc,ip,el))%sizePostResults + &
! 1+thermalState(material_phase(ipc,ip,el))%sizePostResults) :: &
#ifdef multiphysicsOut
real(pReal), dimension(1+crystallite_sizePostResults(microstructure_crystallite(mesh_element(4,el))) + &
1+plasticState(material_phase(ipc,ip,el))%sizePostResults + &
1+damageState(material_phase(ipc,ip,el))%sizePostResults + &
1+thermalState(material_phase(ipc,ip,el))%sizePostResults) :: &
crystallite_postResults
#else
real(pReal), dimension(1+crystallite_sizePostResults(microstructure_crystallite(mesh_element(4,el)))+ &
1+plasticState(material_phase(ipc,ip,el))%sizePostResults) ::
crystallite_postResults
#endif
real(pReal), dimension(3,3) :: &
Ee
real(pReal), dimension(4) :: &
@ -4073,6 +4079,7 @@ function crystallite_postResults(ipc, ip, el)
constitutive_postResults(crystallite_Tstar_v(1:6,ipc,ip,el), crystallite_Fe, &
crystallite_temperature(ip,el), ipc, ip, el)
c = c + plasticState(material_phase(ipc,ip,el))%sizePostResults
#ifdef multiphysicsOut
crystallite_postResults(c+1) = real(damageState(material_phase(ipc,ip,el))%sizePostResults,pReal) ! size of constitutive results
c = c + 1_pInt

11
code/homogenization.f90
View File

@ -646,10 +646,15 @@ subroutine materialpoint_postResults
thePos = thePos + 1_pInt
grainLooping :do g = 1,myNgrains
#ifdef multiphysicsOut
theSize = (1 + crystallite_sizePostResults(myCrystallite)) + &
(1 + plasticState(material_phase(g,i,e))%sizePostResults)! + & !ToDo
! (1 + damageState(material_phase(g,i,e))%sizePostResults) + &
! (1 + thermalState(material_phase(g,i,e))%sizePostResults)
(1 + plasticState(material_phase(g,i,e))%sizePostResults) + & !ToDo
(1 + damageState(material_phase(g,i,e))%sizePostResults) + &
(1 + thermalState(material_phase(g,i,e))%sizePostResults)
#else
theSize = (1 + crystallite_sizePostResults(myCrystallite)) + &
(1 + plasticState(material_phase(g,i,e))%sizePostResults)
#endif
materialpoint_results(thePos+1:thePos+theSize,i,e) = crystallite_postResults(g,i,e) ! tell crystallite results
thePos = thePos + theSize
enddo grainLooping