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

minor changes (syntax polishing)

This commit is contained in:
Philip Eisenlohr 2007-04-11 10:03:43 +00:00
parent e1f2daf753
commit 6095ce0972
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
trunk/constitutive.f90
View File

@ -868,7 +868,7 @@ enddo ! cp_element
end subroutine end subroutine
function constitutive_HomogenizedC(ipc,ip,el) function constitutive_homogenizedC(ipc,ip,el)
!********************************************************************* !*********************************************************************
!* This function returns the homogenized elacticity matrix * !* This function returns the homogenized elacticity matrix *
!* INPUT: * !* INPUT: *
@ -910,9 +910,9 @@ implicit none
!* Definition of variables !* Definition of variables
integer(pInt) ipc,ip,el integer(pInt) ipc,ip,el
integer(pInt) matID,i,k,l,m,n integer(pInt) matID,i,k,l,m,n
real(pReal) Tstar_v(6) real(pReal), dimension(6) :: Tstar_v
real(pReal) Lp(3,3) real(pReal), dimension(3,3) :: Lp
real(pReal) dLp_dTstar(3,3,3,3) real(pReal), dimension(3,3,3,3) :: dLp_dTstar
real(pReal), dimension(constitutive_Nstatevars(ipc,ip,el)) :: state,gdot_slip,dgdot_dtauslip,tau_slip real(pReal), dimension(constitutive_Nstatevars(ipc,ip,el)) :: state,gdot_slip,dgdot_dtauslip,tau_slip
!* Get the material-ID from the triplet(ipc,ip,el) !* Get the material-ID from the triplet(ipc,ip,el)
@ -942,7 +942,7 @@ return
end subroutine end subroutine
function constitutive_DotState(Tstar_v,state,ipc,ip,el) function constitutive_dotState(Tstar_v,state,ipc,ip,el)
!********************************************************************* !*********************************************************************
!* This subroutine contains the constitutive equation for * !* This subroutine contains the constitutive equation for *
!* calculating the velocity gradient * !* calculating the velocity gradient *
@ -962,7 +962,7 @@ implicit none
integer(pInt) ipc,ip,el integer(pInt) ipc,ip,el
integer(pInt) matID,i integer(pInt) matID,i
real(pReal), dimension(6) :: Tstar_v real(pReal), dimension(6) :: Tstar_v
real(pReal), dimension(constitutive_Nstatevars(ipc,ip,el)) :: constitutive_DotState,& real(pReal), dimension(constitutive_Nstatevars(ipc,ip,el)) :: constitutive_dotState,&
state,gdot_slip,tau_slip,self_hardening state,gdot_slip,tau_slip,self_hardening
!* Get the material-ID from the triplet(ipc,ip,el) !* Get the material-ID from the triplet(ipc,ip,el)
@ -978,7 +978,7 @@ do i=1,constitutive_Nstatevars(ipc,ip,el)
enddo enddo
!* Hardening for all systems !* Hardening for all systems
constitutive_DotState=matmul(constitutive_HardeningMatrix(1:material_Nslip(matID),1:material_Nslip(matID),& constitutive_dotState=matmul(constitutive_HardeningMatrix(1:material_Nslip(matID),1:material_Nslip(matID),&
material_CrystalStructure(matID)),self_hardening) material_CrystalStructure(matID)),self_hardening)
return return