Skip to content
Snippets Groups Projects
Commit 4c2c4a84 authored by Dave Whipp's avatar Dave Whipp
Browse files

Fixed even more typos

parent 164b80d8
No related branches found
No related tags found
No related merge requests found
...@@ -41,7 +41,8 @@ type (parameters) params ...@@ -41,7 +41,8 @@ type (parameters) params
double precision :: J2d,J3d,viscosity,pressure,compressibility,straintot double precision :: J2d,J3d,viscosity,pressure,compressibility,straintot
double precision :: friction_angle double precision :: friction_angle
character (len=8) plasticity_type character (len=8) plasticity_type
character (len=12) plasticity_ss_type character (len=12) plasticity_ss_type_coh
character (len=12) plasticity_ss_type_phi
double precision plasticity_parameters(9) double precision plasticity_parameters(9)
logical :: is_plastic,flag_vrm_pb logical :: is_plastic,flag_vrm_pb
...@@ -369,8 +370,8 @@ logical :: flag_vrm_pb ...@@ -369,8 +370,8 @@ logical :: flag_vrm_pb
!------------------------------------------------------------------------------| !------------------------------------------------------------------------------|
double precision :: strain_soft_phi_in,strain_soft_phi_out,strain_soft_phi double precision :: strain_soft_phi_in,strain_soft_phi_out,strain_soft_phi
double precision :: strain_soft_c_in,strain_soft_c_out,strain_soft_phi double precision :: strain_soft_c_in,strain_soft_c_out,strain_soft_c
double precision :: fact double precision :: fact,strain_soft_ref
! Set strain softening onset/termination and min values ! Set strain softening onset/termination and min values
strain_soft_c_in = plasticity_parameters(6) strain_soft_c_in = plasticity_parameters(6)
...@@ -394,8 +395,8 @@ if (strain_soft_c_in.gt.0.d0) then ...@@ -394,8 +395,8 @@ if (strain_soft_c_in.gt.0.d0) then
flag_vrm_pb = .true. flag_vrm_pb = .true.
end select ss_type_c end select ss_type_c
! Soften cohesion if beyond onset value ! Soften cohesion if beyond onset value
if (strain_soft_ref.gt.strain_soft_in_c) then if (strain_soft_ref.gt.strain_soft_c_in) then
fact=(strain_soft_ref-strain_soft_in_c)/(strain_soft_out_c-strain_soft_in_c) fact=(strain_soft_ref-strain_soft_c_in)/(strain_soft_c_out-strain_soft_c_in)
fact=min(fact,1.d0) fact=min(fact,1.d0)
c=c+(strain_soft_c-c)*fact c=c+(strain_soft_c-c)*fact
endif endif
...@@ -415,8 +416,8 @@ if (strain_soft_phi_in.gt.0.d0) then ...@@ -415,8 +416,8 @@ if (strain_soft_phi_in.gt.0.d0) then
flag_vrm_pb = .true. flag_vrm_pb = .true.
end select ss_type_phi end select ss_type_phi
! Soften phi if beyond onset value ! Soften phi if beyond onset value
if (strain_soft_ref.gt.strain_soft_in_phi) then if (strain_soft_ref.gt.strain_soft_phi_in) then
fact=(strain_soft_ref-strain_soft_in_phi)/(strain_soft_out_phi-strain_soft_in_phi) fact=(strain_soft_ref-strain_soft_phi_in)/(strain_soft_phi_out-strain_soft_phi_in)
fact=min(fact,1.d0) fact=min(fact,1.d0)
phi=phi+(strain_soft_phi-phi)*fact phi=phi+(strain_soft_phi-phi)*fact
endif endif
...@@ -425,4 +426,4 @@ endif ...@@ -425,4 +426,4 @@ endif
end subroutine strain_soften end subroutine strain_soften
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment