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

output homogogenized F and P 

was disabled for historic reasons only
This commit is contained in:
Martin Diehl 2022-02-19 15:35:38 +01:00
parent b44a862a8a
commit 968e55b0bc
3 changed files with 18 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
PRIVATE

@ -1 +1 @@
Subproject commit 857b994fb7222ab15a2b8c4ded2bba8787d7feb6 Subproject commit 0f9076b0b589f5134b0304f6676e0076b99d887f

19
src/homogenization_mechanical.f90
View File

@ -182,8 +182,10 @@ module subroutine mechanical_results(group_base,ho)
character(len=*), intent(in) :: group_base character(len=*), intent(in) :: group_base
integer, intent(in) :: ho integer, intent(in) :: ho
integer :: ou
character(len=:), allocatable :: group character(len=:), allocatable :: group
group = trim(group_base)//'/mechanical' group = trim(group_base)//'/mechanical'
call results_closeGroup(results_addGroup(group)) call results_closeGroup(results_addGroup(group))
@ -194,12 +196,17 @@ module subroutine mechanical_results(group_base,ho)
end select end select
!temp = reshape(homogenization_F,[3,3,discretization_nIPs*discretization_Nelems]) do ou = 1, size(output_mechanical(1)%label)
!call results_writeDataset(group,temp,'F',&
! 'deformation gradient','1') select case (output_mechanical(ho)%label(ou))
!temp = reshape(homogenization_P,[3,3,discretization_nIPs*discretization_Nelems]) case('F')
!call results_writeDataset(group,temp,'P',& call results_writeDataset(reshape(homogenization_F,[3,3,discretization_nCells]),group,'F', &
! '1st Piola-Kirchhoff stress','Pa') 'deformation gradient','1')
case('P')
call results_writeDataset(reshape(homogenization_P,[3,3,discretization_nCells]),group,'P', &
'deformation gradient','1')
end select
end do
end subroutine mechanical_results end subroutine mechanical_results

8
src/phase_mechanical.f90
View File

@ -331,7 +331,7 @@ module subroutine mechanical_results(group,ph)
integer, intent(in) :: ph integer, intent(in) :: ph
call crystallite_results(group,ph) call results(group,ph)
select case(phase_plasticity(ph)) select case(phase_plasticity(ph))
@ -882,7 +882,7 @@ end function integrateStateRK
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief writes crystallite results to HDF5 output file !> @brief writes crystallite results to HDF5 output file
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine crystallite_results(group,ph) subroutine results(group,ph)
character(len=*), intent(in) :: group character(len=*), intent(in) :: group
integer, intent(in) :: ph integer, intent(in) :: ph
@ -926,7 +926,7 @@ subroutine crystallite_results(group,ph)
if (any(phase_lattice(ph) == ['hP', 'tI'])) & if (any(phase_lattice(ph) == ['hP', 'tI'])) &
call results_addAttribute('c/a',phase_cOverA(ph),group//'/mechanical/'//output_mechanical(ph)%label(ou)) call results_addAttribute('c/a',phase_cOverA(ph),group//'/mechanical/'//output_mechanical(ph)%label(ou))
end select end select
enddo end do
contains contains
@ -948,7 +948,7 @@ subroutine crystallite_results(group,ph)
end function to_quaternion end function to_quaternion
end subroutine crystallite_results end subroutine results
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------