Merge branch 'main' of https://gitlab.uliege.be/moammm/moammmPublic/code

MaterialModels/FiniteStrain/Finite_VEVP/1e_strain_cyc.inp

@@ -35,6 +35,22 @@
 3380000.0,33800000.0
 *DEPVAR
 108
+*INITIAL CONDITIONS, TYPE=SOLUTION
+1.0,1.0,1.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0
 **
 *STEP, name=Step-1, inc=1000, NLGEOM=YES
 *Static

MaterialModels/FiniteStrain/Finite_VEVP/1e_strain_mono.inp

@@ -35,6 +35,22 @@
 3380000.0,33800000.0
 *DEPVAR
 108
+*INITIAL CONDITIONS, TYPE=SOLUTION
+1.0,1.0,1.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0
 **
 *STEP, name=Step-1, inc=1000, NLGEOM=YES
 *Static

MaterialModels/FiniteStrain/Finite_VEVP/README.md

@@ -67,6 +67,11 @@ The UMAT implementation uses the material parameters **_in the following order_*
 |49 | Retardation Time Deviatoric - Maxwell Branch 7  |      $g_7$                    |    s        |
 |50 | Retardation Time Deviatoric - Maxwell Branch 8  |      $g_8$                    |    s        |
 
+## Remark
+
+Lines 164-168 have to be uncommented for use in FFTMAD but have to be commented for ABAQUS or cm3Libraries.
+If you comment out, then in ABAQUS, use the inp file to set initial values for F_vp.
+
 ## Disclaimer
 
 Code related to MOAMMM project

MaterialModels/FiniteStrain/Finite_VEVP/dogbone.inp

@@ -4990,6 +4990,22 @@
 3380000.0,33800000.0
 *DEPVAR
 108
+*INITIAL CONDITIONS, TYPE=SOLUTION
+1.0,1.0,1.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,
+0.0
 ** -----------------------------------------------------------------------------
 **
 ** BOUNDARY CONDITIONS

MaterialModels/FiniteStrain/Finite_VEVP/umat.f

@@ -161,11 +161,11 @@ C     !--------------------------------------------------------------
         ! for the first time step (initial condition) in FFTMAD.
         ! Can be commented out for ABAQUS, If you comment out, then
         ! in abaqus use the inp file to set initial values for F_vp.
-         IF ((KINC .EQ. 1) .AND. (KSTEP .EQ. 1)) THEN
-          STATEV(1) = 1.D0
-          STATEV(2) = 1.D0
-          STATEV(3) = 1.D0
-         END IF
+C         IF ((KINC .EQ. 1) .AND. (KSTEP .EQ. 1)) THEN
+C          STATEV(1) = 1.D0
+C          STATEV(2) = 1.D0
+C          STATEV(3) = 1.D0
+C         END IF
 
          ! State variables at previous time step
          CALL voit2mat(STATEV(1:9), F_vp_n(:,:))