Modification of /MAT/LAW190 to improve stability in compression

laurent-altr · Aug 22, 2024 · 548b965 · 548b965
1 parent c7f0712
commit 548b965
Show file tree

Hide file tree

Showing 9 changed files with 481 additions and 535 deletions.
diff --git a/engine/source/materials/mat/mat090/sigeps90.F b/engine/source/materials/mat/mat090/sigeps90.F
@@ -391,7 +391,7 @@ SUBROUTINE SIGEPS90(
               END DO
 
             END DO  ! K=1,3
-              
+
             E(1:NEL) = MAX(E0,EMAX(1:NEL))
             E(1:NEL) = EMAX(1:NEL)
             ET(1:NEL) = EMIN(1:NEL)/E0  ! Not used
@@ -417,6 +417,7 @@ SUBROUTINE SIGEPS90(
             ILEN(1:NEL) = NPF(IFUNC(1)+1) / 2 - IAD(1:NEL) - VARTMP(1:NEL,KK)
 
             CALL VINTER2(TF,IAD,VARTMP(1,KK),ILEN,NEL,STRAIN(1,K),DYDX,SQSTAT(1,K))
+            SQSTAT(1:NEL,K) = YFAC1 * SQSTAT(1:NEL,K)
 
             DO I=1,NEL
               IF (TFLAG == 2 .AND. STRAIN(I,K) < ZERO) SQSTAT(I,K) = E0*STRAIN(I,K)  

diff --git a/engine/source/materials/mat/mat190/conversion.F b/engine/source/materials/mat/mat190/conversion.F