correct traction in linkslip

githubgrasp · githubgrasp · commit e017811fe03d · 2025-09-08T15:32:05.000+01:00
diff --git a/src/sm/Materials/InterfaceMaterials/linkslip.C b/src/sm/Materials/InterfaceMaterials/linkslip.C
@@ -168,7 +168,7 @@ LinkSlip :: giveEngTraction_3d(const FloatArrayF<3> &jump, GaussPoint *gp, TimeS
 
     if ( f > 0 ) { //plastic response.
         //Reduced stress by increasing plastic strain.
-        traction.at(1) = evaluateBondStress(tempKappa);;
+      traction.at(1) = sgn(traction.at(1)) * evaluateBondStress(tempKappa);
     }
 
     //Compute the lateral stress components

Original file line number	Diff line number	Diff line change
`@@ -168,7 +168,7 @@ LinkSlip :: giveEngTraction_3d(const FloatArrayF<3> &jump, GaussPoint *gp, TimeS`
`168`	`168`
`169`	`169`	`if ( f > 0 ) { //plastic response.`
`170`	`170`	`//Reduced stress by increasing plastic strain.`
`171`		`- traction.at(1) = evaluateBondStress(tempKappa);;`
	`171`	`+ traction.at(1) = sgn(traction.at(1)) * evaluateBondStress(tempKappa);`
`172`	`172`	`}`
`173`	`173`
`174`	`174`	`//Compute the lateral stress components`