sign-transfer done right

This commit is contained in:
Martin Diehl 2021-06-21 23:26:17 +02:00
parent 4b82cb4c14
commit c49ba5419b
2 changed files with 3 additions and 4 deletions

View File

@ -282,8 +282,7 @@ module subroutine isotropic_dotState(Mp,ph,en)
endif
dot%xi(en) = dot_gamma &
* ( prm%h_0 + prm%h_ln * log(dot_gamma) ) &
* abs( 1.0_pReal - stt%xi(en)/xi_inf_star )**prm%a *prm%h &
* sign(1.0_pReal, 1.0_pReal - stt%xi(en)/xi_inf_star)
* sign(abs(1.0_pReal - stt%xi(en)/xi_inf_star)**prm%a *prm%h, 1.0_pReal-stt%xi(en)/xi_inf_star)
else
dot%xi(en) = 0.0_pReal
endif

View File

@ -364,8 +364,8 @@ module subroutine phenopowerlaw_dotState(Mp,ph,en)
! calculate left and right vectors
left_SlipSlip = 1.0_pReal + prm%h_int
xi_slip_sat_offset = prm%f_sat_sl_tw*sqrt(sumF)
right_SlipSlip = abs(1.0_pReal-stt%xi_slip(:,en) / (prm%xi_inf_sl+xi_slip_sat_offset)) **prm%a_sl &
* sign(1.0_pReal,1.0_pReal-stt%xi_slip(:,en) / (prm%xi_inf_sl+xi_slip_sat_offset))
right_SlipSlip = sign(abs(1.0_pReal-stt%xi_slip(:,en) / (prm%xi_inf_sl+xi_slip_sat_offset)) **prm%a_sl, &
1.0_pReal-stt%xi_slip(:,en) / (prm%xi_inf_sl+xi_slip_sat_offset))
!--------------------------------------------------------------------------------------------------
! shear rates