From 1ed276ce2d8e65bb5425325557a0fc76991ee02b Mon Sep 17 00:00:00 2001 From: Pratheek Shanthraj Date: Thu, 30 Jul 2015 10:29:25 +0000 Subject: [PATCH] quadratic softening for finite dissipation energy --- code/kinematics_cleavage_opening.f90 | 2 +- code/source_damage_anisoBrittle.f90 | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/code/kinematics_cleavage_opening.f90 b/code/kinematics_cleavage_opening.f90 index 8ff4ca2ba..bc4c0c661 100644 --- a/code/kinematics_cleavage_opening.f90 +++ b/code/kinematics_cleavage_opening.f90 @@ -254,7 +254,7 @@ subroutine kinematics_cleavage_opening_LiAndItsTangent(Ld, dLd_dTstar3333, Tstar traction_t = dot_product(Tstar_v,lattice_Scleavage_v(1:6,2,index_myFamily+i,phase)) traction_n = dot_product(Tstar_v,lattice_Scleavage_v(1:6,3,index_myFamily+i,phase)) traction_crit = kinematics_cleavage_opening_critLoad(f,instance)* & - damage(homog)%p(damageOffset) + damage(homog)%p(damageOffset)*damage(homog)%p(damageOffset) udotd = & sign(1.0_pReal,traction_d)* & kinematics_cleavage_opening_sdot_0(instance)* & diff --git a/code/source_damage_anisoBrittle.f90 b/code/source_damage_anisoBrittle.f90 index 2df036a64..225b73875 100644 --- a/code/source_damage_anisoBrittle.f90 +++ b/code/source_damage_anisoBrittle.f90 @@ -336,7 +336,7 @@ subroutine source_damage_anisoBrittle_dotState(Tstar_v, ipc, ip, el) traction_n = dot_product(Tstar_v,lattice_Scleavage_v(1:6,3,index_myFamily+i,phase)) traction_crit = source_damage_anisoBrittle_critLoad(f,instance)* & - damage(homog)%p(damageOffset) + damage(homog)%p(damageOffset)*damage(homog)%p(damageOffset) sourceState(phase)%p(sourceOffset)%dotState(1,constituent) = & sourceState(phase)%p(sourceOffset)%dotState(1,constituent) + & source_damage_anisoBrittle_sdot_0(instance)* &