diff --git a/Modules/Abaqus/Explicit/abaqus_test.sub b/Modules/Abaqus/Explicit/abaqus_test.sub
new file mode 100644
index 0000000..f4b0a70
--- /dev/null
+++ b/Modules/Abaqus/Explicit/abaqus_test.sub
@@ -0,0 +1,35 @@
+#!/bin/bash
+#SBATCH --job-name Abaqus
+#SBATCH --nodes 2
+#SBATCH --ntasks-per-node 2
+#SBATCH --cpus-per-task 4
+#SBATCH --time 00:10:00
+ 
+module purge
+module load intel/14.0.1 intelmpi/4.1.3 abaqus/6.14-1/intel-14.0.1
+
+# Generating the node list from the SLURM variables
+cpus_per_node=$(echo $SLURM_TASKS_PER_NODE | sed -e 's/(x[0-9])//')
+if [ "x${SLURM_CPUS_PER_TASK}" != "x" ]; then
+   cpus_per_node=$(( ${cpus_per_node} * ${SLURM_CPUS_PER_TASK} ))
+fi
+
+nodes_list=""
+for i in $(scontrol show hostnames | sort -u); do
+   if [ "x${nodes_list}" != "x" ]; then
+      nodes_list="${nodes_list},"
+   fi
+   nodes_list="${nodes_list}['$i',${cpus_per_node}]"
+done
+
+mp_host_list="[${nodes_list}]"
+
+cat <<EOF >> abaqus_v6.env
+mp_host_list=${mp_host_list}
+EOF
+
+total_nb_cpus=$(( ${cpus_per_node} * ${SLURM_NNODES} ))
+
+echo CPUS=$total_nb_cpus NODES=$mp_host_list
+
+abaqus inter cpus=$total_nb_cpus j=discbrake_3d_xpl
\ No newline at end of file
diff --git a/Modules/Abaqus/Explicit/discbrake_3d_xpl.inp b/Modules/Abaqus/Explicit/discbrake_3d_xpl.inp
new file mode 100644
index 0000000..c237f9e
--- /dev/null
+++ b/Modules/Abaqus/Explicit/discbrake_3d_xpl.inp
@@ -0,0 +1,224 @@
+*HEADING
+ THERMALLY COUPLED ANALYSIS OF A SYMMETRIC DISC
+ BRAKE (EXPLICIT)
+ RIGID BODY (PIN NSET)
+*NODE
+   1,  90.0000000E-3,   0.00000000E-3, 0.0E-3
+  36,  88.6326978E-3, -15.62833599E-3, 0.0E-3
+ 101, 100.0000000E-3,   0.00000000E-3, 0.0E-3
+ 136,  98.4807753E-3, -17.36481777E-3, 0.0E-3
+ 401, 135.0000000E-3,   0.00000000E-3, 0.0E-3
+ 436, 132.9490467E-3, -23.44250399E-3, 0.0E-3
+1001,  90.0000000E-3,   0.00000000E-3, 3.0E-3
+1036,  88.6326978E-3, -15.62833599E-3, 3.0E-3
+1101, 100.0000000E-3,   0.00000000E-3, 3.0E-3
+1136,  98.4807753E-3, -17.36481777E-3, 3.0E-3
+1401, 135.0000000E-3,   0.00000000E-3, 3.0E-3
+1436, 132.9490467E-3, -23.44250399E-3, 3.0E-3
+2101, 100.0000000E-3,   0.00000000E-3, 4.5E-3
+2136,  98.4807753E-3, -17.36481777E-3, 4.5E-3
+2401, 135.0000000E-3,   0.00000000E-3, 4.5E-3
+2436, 132.9490467E-3, -23.44250399E-3, 4.5E-3
+3101, 100.0000000E-3,   0.00000000E-3, 5.0E-3
+3136,  98.4807753E-3, -17.36481777E-3, 5.0E-3
+3401, 135.0000000E-3,   0.00000000E-3, 5.0E-3
+3436, 132.9490467E-3, -23.44250399E-3, 5.0E-3
+*NODE, NSET=CENTER
+1000, 0., 0., 0., 
+*NGEN, LINE=C, NSET=SUPPORT
+   1,   36, 1,  , 0., 0., 0.0E-3, 0., 0., 100.E-3
+1001, 1036, 1,  , 0., 0., 3.0E-3, 0., 0., 100.E-3
+*NGEN, LINE=C, NSET=INNER1
+ 101,  136, 1,  , 0., 0., 0.0E-3, 0., 0., 100.E-3
+1101, 1136, 1,  , 0., 0., 3.0E-3, 0., 0., 100.E-3
+*NGEN, LINE=C, NSET=INNER2
+2101, 2136, 1,  , 0., 0., 4.5E-3, 0., 0., 100.E-3
+3101, 3136, 1,  , 0., 0., 5.0E-3, 0., 0., 100.E-3
+*NSET, NSET=INNER
+INNER1, INNER2
+*NGEN, LINE=C, NSET=OUTER
+ 401,  436, 1,  , 0., 0., 0.0E-3, 0., 0., 100.E-3
+1401, 1436, 1,  , 0., 0., 3.0E-3, 0., 0., 100.E-3
+2401, 2436, 1,  , 0., 0., 4.5E-3, 0., 0., 100.E-3
+3401, 3436, 1,  , 0., 0., 5.0E-3, 0., 0., 100.E-3
+*NFILL, NSET=DISK
+INNER, OUTER, 3, 100
+*NSET, NSET=SYMM, GENERATE
+101, 436
+*ELEMENT, TYPE=C3D8RT, ELSET=SUPPORT
+  1,   1, 101, 102,   2, 1001, 1101, 1102, 1002
+ 36,  36, 136, 101,   1, 1036, 1136, 1101, 1001
+*ELEMENT, TYPE=C3D8RT, ELSET=DISK
+101, 101, 201, 202, 102, 1101, 1201, 1202, 1102
+136, 136, 236, 201, 101, 1136, 1236, 1201, 1101
+*ELGEN, ELSET=SUPPORT
+  1, 35, 1, 1
+*ELGEN, ELSET=DISK
+101, 35, 1, 1, 3, 100, 100, 3, 1000, 1000
+136,  1, 1, 1, 3, 100, 100, 3, 1000, 1000
+*ELSET, ELSET=FULLDISK
+SUPPORT, DISK
+*ELSET, ELSET=DISKSURF, GENERATE
+2101, 2436
+*SOLID SECTION,ELSET=FULLDISK,MATERIAL=STEEL,
+ORIENTATION=CYL
+*MATERIAL, NAME=STEEL
+*ELASTIC
+209.E9, .3
+*EXPANSION
+11.E-6, 
+*DENSITY
+7800., 
+*SPECIFIC HEAT
+452., 
+*CONDUCTIVITY
+48., 
+*ORIENTATION, SYSTEM=CYLINDRICAL, NAME=CYL
+0., 0., 0., 0., 0., 1.
+3, 0.
+*RIGID BODY, REF NODE=1000, PIN NSET=SUPPORT
+*NODE
+10001,  29.E-3, 97.26895702E-3,  5.0E-3
+10005, -29.E-3, 97.26895702E-3,  5.0E-3
+10301,  38.E-3, 127.4558747E-3,  5.0E-3
+10305, -38.E-3, 127.4558747E-3,  5.0E-3
+11001,  29.E-3, 97.26895702E-3,  5.5E-3
+11005, -29.E-3, 97.26895702E-3,  5.5E-3
+11301,  38.E-3, 127.4558747E-3,  5.5E-3
+11305, -38.E-3, 127.4558747E-3,  5.5E-3
+12001,  29.E-3, 97.26895702E-3,  7.0E-3
+12005, -29.E-3, 97.26895702E-3,  7.0E-3
+12301,  38.E-3, 127.4558747E-3,  7.0E-3
+12305, -38.E-3, 127.4558747E-3,  7.0E-3
+13001,  29.E-3, 97.26895702E-3, 10.0E-3
+13005, -29.E-3, 97.26895702E-3, 10.0E-3
+13301,  38.E-3, 127.4558747E-3, 10.0E-3
+13305, -38.E-3, 127.4558747E-3, 10.0E-3
+14001,  29.E-3, 97.26895702E-3, 15.0E-3
+14005, -29.E-3, 97.26895702E-3, 15.0E-3
+14301,  38.E-3, 127.4558747E-3, 15.0E-3
+14305, -38.E-3, 127.4558747E-3, 15.0E-3
+*NGEN, LINE=C, NSET=INPAD
+10001,10005,1,,0.,0.,5.0E-3,0.,0.,100.E-3
+11001,11005,1,,0.,0.,5.5E-3,0.,0.,100.E-3
+12001,12005,1,,0.,0.,7.0E-3,0.,0.,100.E-3
+13001,13005,1,,0.,0.,10.0E-3,0.,0.,100.E-3
+14001,14005,1,,0.,0.,15.0E-3,0.,0.,100.E-3
+*NGEN, LINE=C, NSET=OUTPAD
+10301,10305,1,,0.,0.,5.0E-3,0.,0.,100.E-3
+11301,11305,1,,0.,0.,5.5E-3,0.,0.,100.E-3
+12301,12305,1,,0.,0.,7.0E-3,0.,0.,100.E-3
+13301,13305,1,,0.,0.,10.0E-3,0.,0.,100.E-3
+14301,14305,1,,0.,0.,15.0E-3,0.,0.,100.E-3
+*NFILL, NSET=NPAD
+INPAD, OUTPAD, 3, 100
+*NSET, NSET=NPADBACK, GENERATE
+14001, 14305
+*NSET, NSET=ALL
+DISK, NPAD, SUPPORT
+*ELEMENT, TYPE=C3D8RT
+10001, 10001, 10101, 10102, 10002, 11001, 11101,
+11102, 11002
+*ELGEN, ELSET=EPAD
+10001, 4, 1, 1, 3, 100, 100, 4, 1000, 1000
+*ELSET, ELSET=EPADBACK, GENERATE
+13001, 13204
+*SOLID SECTION, ELSET=EPAD, MATERIAL=PADMAT,
+ORIENTATION=CYL
+*MATERIAL, NAME=PADMAT
+*ELASTIC
+2200.E6, .25,  20.0
+1300.E6, .25, 100.0
+ 530.E6, .25, 200.0
+ 320.E6, .25, 300.0
+*EXPANSION
+10.E-6, 20.0
+30.E-6,200.0
+*DENSITY
+1550., 
+*SPECIFIC HEAT
+1200., 
+*CONDUCTIVITY
+0.9, 
+*INITIAL CONDITIONS, TYPE=TEMPERATURE
+ALL,20.
+*ELSET,ELSET=EALL
+EPAD,FULLDISK
+** DEFINE SURFACES 
+*SURFACE, NAME=MASTER
+DISKSURF, S2
+*SURFACE, NAME=SLAVE
+10001, S1
+10002, S1
+10003, S1
+10004, S1
+10101, S1
+10102, S1
+10103, S1
+10104, S1
+10201, S1
+10202, S1
+10203, S1
+10204, S1
+*ELSET, ELSET=QA_TEST, GEN
+2201, 2209, 1
+*NSET, NSET=QA_TEST, GEN 
+3201, 3210, 1
+*STEP
+ PRESS THE PAD AGAINST THE DISK
+*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
+ , 0.001
+*BOUNDARY
+SYMM, 3, 3
+CENTER, 1, 6
+NPADBACK, 1, 2
+** DEFINE CONTACT PAIR
+*CONTACT PAIR, INTERACTION=INTERF,WEIGHT=0.0
+SLAVE,MASTER
+*SURFACE INTERACTION, NAME=INTERF
+*FRICTION
+0.37,,,20.0
+0.38,,,100.0
+0.41,,,200.0
+0.39,,,300.0
+0.24,,,400.0
+*GAP CONDUCTANCE, PRESSURE
+5.0E7,0.0
+5.0E7,1000.0
+*GAP HEAT GENERATION
+1.0, 
+*AMPLITUDE,NAME=RAMP
+0.0,0.0,0.001,1.0
+*DLOAD,AMP=RAMP
+EPADBACK, P2, 1.7E6
+*FILM,AMP=RAMP
+DISKSURF, F2, 20., 100.
+****
+****    OUTPUT REQUESTS
+****
+*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=1 
+*ELEMENT OUTPUT, ELSET=EALL
+LE, 
+*NODE OUTPUT, NSET=ALL
+NT, 
+****    OUTPUT FOR QA PURPOSES
+*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=1 
+*NODE OUTPUT, NSET=QA_TEST
+NT, 
+*OUTPUT, HIST, TIME INTERVAL=8E-5
+*ENERGY OUTPUT, VAR=PRESELECT 
+*FILE OUTPUT,NUM=1
+*NODE FILE, NSET=QA_TEST
+NT, 
+*EL FILE, ELSET=QA_TEST
+LE,
+*END STEP
+*STEP
+ ROTATE THE DISK BY 60 DEGREES
+*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
+ , 0.015
+*AMPLITUDE,NAME=RAMP2,DEFINITION=SMOOTH STEP
+0.,0.,0.015,1.0
+*BOUNDARY, TYPE=DISPLACEMENT,AMP=RAMP2
+CENTER, 6, 6, 1.04717
+*END STEP
diff --git a/Modules/Abaqus/Standard/abaqus_test.sub b/Modules/Abaqus/Standard/abaqus_test.sub
new file mode 100644
index 0000000..d645553
--- /dev/null
+++ b/Modules/Abaqus/Standard/abaqus_test.sub
@@ -0,0 +1,12 @@
+#!/bin/bash 
+#SBATCH --job-name Abaqus
+#SBATCH --nodes 1
+#SBATCH --ntasks 1
+#SBATCH --cpus-per-task 8
+#SBATCH --time 00:10:00
+ 
+module purge
+module load intel/14.0.1 intelmpi/4.1.3 abaqus/6.14-1/intel-14.0.1
+
+# Abaqus/standard uses thread parallelism
+abaqus inter cpus=$SLURM_CPUS_PER_TASK j=umatmst3 user=umatmst3.f
\ No newline at end of file
diff --git a/Modules/Abaqus/Standard/umatmst3.f b/Modules/Abaqus/Standard/umatmst3.f
new file mode 100644
index 0000000..2eca1ad
--- /dev/null
+++ b/Modules/Abaqus/Standard/umatmst3.f
@@ -0,0 +1,220 @@
+      SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
+     1 RPL,DDSDDT,DRPLDE,DRPLDT,STRAN,DSTRAN,
+     2 TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,MATERL,NDI,NSHR,NTENS,
+     3 NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,CELENT,
+     4 DFGRD0,DFGRD1,NOEL,NPT,KSLAY,KSPT,KSTEP,KINC)
+C
+      INCLUDE 'ABA_PARAM.INC'
+C
+      CHARACTER*80 MATERL
+      DIMENSION STRESS(NTENS),STATEV(NSTATV),
+     1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),
+     2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),
+     3 PROPS(NPROPS),COORDS(3),DROT(3,3),
+     4 DFGRD0(3,3),DFGRD1(3,3)
+C
+      DIMENSION EELAS(6),EPLAS(6),FLOW(6)
+      PARAMETER (ONE=1.0D0,TWO=2.0D0,THREE=3.0D0,SIX=6.0D0)
+      DATA NEWTON,TOLER/10,1.D-6/
+C
+C -----------------------------------------------------------
+C     UMAT FOR ISOTROPIC ELASTICITY AND ISOTROPIC PLASTICITY
+C     J2 FLOW THEORY
+C     CAN NOT BE USED FOR PLANE STRESS
+C -----------------------------------------------------------
+C     PROPS(1) - E
+C     PROPS(2) - NU
+C     PROPS(3) - SYIELD
+C     CALLS AHARD FOR CURVE OF SYIELD VS. PEEQ
+C -----------------------------------------------------------
+C
+      IF (NDI.NE.3) THEN
+         WRITE(6,1)
+ 1       FORMAT(//,30X,'***ERROR - THIS UMAT MAY ONLY BE USED FOR ',
+     1          'ELEMENTS WITH THREE DIRECT STRESS COMPONENTS')
+      ENDIF
+C
+C     ELASTIC PROPERTIES
+C
+      EMOD=PROPS(1)
+      ENU=PROPS(2)
+      IF(ENU.GT.0.4999.AND.ENU.LT.0.5001) ENU=0.499
+      EBULK3=EMOD/(ONE-TWO*ENU)
+      EG2=EMOD/(ONE+ENU)
+      EG=EG2/TWO
+      EG3=THREE*EG
+      ELAM=(EBULK3-EG2)/THREE
+C
+C     ELASTIC STIFFNESS
+C
+      DO 20 K1=1,NTENS
+        DO 10 K2=1,NTENS
+           DDSDDE(K2,K1)=0.0
+ 10     CONTINUE
+ 20   CONTINUE
+C
+      DO 40 K1=1,NDI
+        DO 30 K2=1,NDI
+           DDSDDE(K2,K1)=ELAM
+ 30     CONTINUE
+        DDSDDE(K1,K1)=EG2+ELAM
+ 40   CONTINUE
+      DO 50 K1=NDI+1,NTENS
+        DDSDDE(K1,K1)=EG
+ 50   CONTINUE
+C
+C    CALCULATE STRESS FROM ELASTIC STRAINS
+C
+      DO 70 K1=1,NTENS
+        DO 60 K2=1,NTENS
+           STRESS(K2)=STRESS(K2)+DDSDDE(K2,K1)*DSTRAN(K1)
+ 60     CONTINUE
+ 70   CONTINUE
+C
+C    RECOVER ELASTIC AND PLASTIC STRAINS
+C
+      DO 80 K1=1,NTENS
+         EELAS(K1)=STATEV(K1)+DSTRAN(K1)
+         EPLAS(K1)=STATEV(K1+NTENS)
+ 80   CONTINUE
+      EQPLAS=STATEV(1+2*NTENS)
+C
+C    IF NO YIELD STRESS IS GIVEN, MATERIAL IS TAKEN TO BE ELASTIC
+C
+      IF(NPROPS.GT.2.AND.PROPS(3).GT.0.0) THEN
+C
+C       MISES STRESS
+C
+        SMISES=(STRESS(1)-STRESS(2))*(STRESS(1)-STRESS(2)) +
+     1          (STRESS(2)-STRESS(3))*(STRESS(2)-STRESS(3)) +
+     1          (STRESS(3)-STRESS(1))*(STRESS(3)-STRESS(1))
+        DO 90 K1=NDI+1,NTENS
+              SMISES=SMISES+SIX*STRESS(K1)*STRESS(K1)
+ 90     CONTINUE
+        SMISES=SQRT(SMISES/TWO)
+C
+C       HARDENING CURVE, GET YIELD STRESS
+C
+        NVALUE=NPROPS/2-1
+        CALL AHARD(SYIEL0,HARD,EQPLAS,PROPS(3),NVALUE)
+C
+C       DETERMINE IF ACTIVELY YIELDING
+C
+        IF (SMISES.GT.(1.0+TOLER)*SYIEL0) THEN
+C
+C         FLOW DIRECTION
+C
+          SHYDRO=(STRESS(1)+STRESS(2)+STRESS(3))/THREE
+          ONESY=ONE/SMISES
+          DO 110 K1=1,NDI
+             FLOW(K1)=ONESY*(STRESS(K1)-SHYDRO)
+ 110      CONTINUE
+          DO 120 K1=NDI+1,NTENS
+             FLOW(K1)=STRESS(K1)*ONESY
+ 120      CONTINUE
+C
+C       SOLVE FOR EQUIV STRESS, NEWTON ITERATION
+C
+          SYIELD=SYIEL0
+          DEQPL=0.0
+          DO 130 KEWTON=1,NEWTON
+             RHS=SMISES-EG3*DEQPL-SYIELD
+             DEQPL=DEQPL+RHS/(EG3+HARD)
+             CALL AHARD(SYIELD,HARD,EQPLAS+DEQPL,PROPS(3),NVALUE)
+             IF(ABS(RHS).LT.TOLER*SYIEL0) GOTO 140
+ 130      CONTINUE
+          WRITE(6,2) NEWTON
+ 2        FORMAT(//,30X,'***WARNING - PLASTICITY ALGORITHM DID NOT ',
+     1        'CONVERGE AFTER ',I3,' ITERATIONS')
+ 140      CONTINUE
+          EFFHRD=EG3*HARD/(EG3+HARD)
+C
+C       CALC STRESS AND UPDATE STRAINS
+C
+          DO 150 K1=1,NDI
+             STRESS(K1)=FLOW(K1)*SYIELD+SHYDRO
+             EPLAS(K1)=EPLAS(K1)+THREE*FLOW(K1)*DEQPL/TWO
+             EELAS(K1)=EELAS(K1)-THREE*FLOW(K1)*DEQPL/TWO
+ 150      CONTINUE
+          DO 160 K1=NDI+1,NTENS
+             STRESS(K1)=FLOW(K1)*SYIELD
+             EPLAS(K1)=EPLAS(K1)+THREE*FLOW(K1)*DEQPL
+             EELAS(K1)=EELAS(K1)-THREE*FLOW(K1)*DEQPL
+ 160      CONTINUE
+          EQPLAS=EQPLAS+DEQPL
+          SPD=DEQPL*(SYIEL0+SYIELD)/TWO
+C
+C       JACOBIAN
+C
+          EFFG=EG*SYIELD/SMISES
+          EFFG2=TWO*EFFG
+          EFFG3=THREE*EFFG2/TWO
+          EFFLAM=(EBULK3-EFFG2)/THREE
+          DO 220 K1=1,NDI
+             DO 210 K2=1,NDI
+                DDSDDE(K2,K1)=EFFLAM
+ 210         CONTINUE
+             DDSDDE(K1,K1)=EFFG2+EFFLAM
+ 220      CONTINUE
+          DO 230 K1=NDI+1,NTENS
+             DDSDDE(K1,K1)=EFFG
+ 230      CONTINUE
+          DO 250 K1=1,NTENS
+             DO 240 K2=1,NTENS
+                DDSDDE(K2,K1)=DDSDDE(K2,K1)+FLOW(K2)*FLOW(K1)
+     1                                       *(EFFHRD-EFFG3)
+ 240         CONTINUE
+ 250      CONTINUE
+        ENDIF
+      ENDIF
+C
+C STORE STRAINS IN STATE VARIABLE ARRAY
+C
+      DO 310 K1=1,NTENS
+        STATEV(K1)=EELAS(K1)
+        STATEV(K1+NTENS)=EPLAS(K1)
+ 310  CONTINUE
+      STATEV(1+2*NTENS)=EQPLAS
+C
+      RETURN
+      END
+C
+C
+      SUBROUTINE AHARD(SYIELD,HARD,EQPLAS,TABLE,NVALUE)
+C
+      INCLUDE 'ABA_PARAM.INC'
+      DIMENSION TABLE(2,NVALUE)
+C
+C    SET YIELD STRESS TO LAST VALUE OF TABLE, HARDENING TO ZERO
+      SYIELD=TABLE(1,NVALUE)
+      HARD=0.0
+C
+C   IF MORE THAN ONE ENTRY, SEARCH TABLE
+C
+      IF(NVALUE.GT.1) THEN
+        DO 10 K1=1,NVALUE-1
+           EQPL1=TABLE(2,K1+1)
+           IF(EQPLAS.LT.EQPL1) THEN
+             EQPL0=TABLE(2,K1)
+             IF(EQPL1.LE.EQPL0) THEN
+                WRITE(6,1)
+ 1              FORMAT(//,30X,'***ERROR - PLASTIC STRAIN MUST BE ',
+     1                 'ENTERED IN ASCENDING ORDER')
+                CALL XIT
+              ENDIF
+C
+C           CURRENT YIELD STRESS AND HARDENING
+C
+            DEQPL=EQPL1-EQPL0
+            SYIEL0=TABLE(1,K1)
+            SYIEL1=TABLE(1,K1+1)
+            DSYIEL=SYIEL1-SYIEL0
+            HARD=DSYIEL/DEQPL
+            SYIELD=SYIEL0+(EQPLAS-EQPL0)*HARD
+            GOTO 20
+            ENDIF
+ 10         CONTINUE
+ 20         CONTINUE
+      ENDIF
+      RETURN
+      END
diff --git a/Modules/Abaqus/Standard/umatmst3.inp b/Modules/Abaqus/Standard/umatmst3.inp
new file mode 100644
index 0000000..1008f06
--- /dev/null
+++ b/Modules/Abaqus/Standard/umatmst3.inp
@@ -0,0 +1,85 @@
+*HEADING
+ UMAT - MISES PLASTICITY, UNIAXIAL TENSION TEST, C3D8 --- 
+*NODE,NSET=ALLN
+1,0.,0.,0.
+2,1.,0.,0.
+3,1.,1.,0.
+4,0.,1.,0.
+5,0.,0.,1.
+6,1.,0.,1.
+7,1.,1.,1.
+8,0.,1.,1.
+*NODE,NSET=ALLN1
+11,0.,0.,0.
+12,1.,0.,0.
+13,1.,1.,0.
+14,0.,1.,0.
+15,0.,0.,1.
+16,1.,0.,1.
+17,1.,1.,1.
+18,0.,1.,1.
+*ELEMENT,TYPE=C3D8,ELSET=ALLE
+1,1,2,3,4,5,6,7,8
+*ELEMENT,TYPE=C3D8,ELSET=ALLE1
+10,11,12,13,14,15,16,17,18
+*SOLID SECTION,ELSET=ALLE,MATERIAL=ALLE
+*SOLID SECTION,ELSET=ALLE1,MATERIAL=ALLE1
+*MATERIAL,NAME=ALLE
+*USER MATERIAL,CONSTANTS=8
+200.E3,.3,200.,0.,220.,.0009,220.,.0029
+*DEPVAR
+13, 
+*MATERIAL,NAME=ALLE1
+*USER MATERIAL,CONSTANTS=8,UNSYMM,TYPE=MECHANICAL
+200.E3,.3,200.,0.,220.,.0009,220.,.0029
+*DEPVAR
+13, 
+*BOUNDARY
+1,PINNED
+2,2
+5,2
+6,2
+4,1
+5,1
+8,1
+2,3
+3,3
+4,3
+11,PINNED
+12,2
+15,2
+16,2
+14,1
+15,1
+18,1
+12,3
+13,3
+14,3
+*STEP,NLGEOM,INC=20
+*STATIC,DIRECT
+1.,20.
+*BOUNDARY
+7,3,,.004
+5,3,,.004
+6,3,,.004
+8,3,,.004
+*EL PRINT
+S, 
+SINV, 
+E, 
+EE, 
+SDV, 
+*NODE PRINT
+U,RF
+*EL FILE,FREQ=10
+S,E,SDV
+*END STEP
+*STEP,NLGEOM,INC=20,UNSYMM=YES
+*STATIC,DIRECT
+1.,20.
+*BOUNDARY
+17,3,,.004
+15,3,,.004
+16,3,,.004
+18,3,,.004
+*END STEP