type('61.inp')
*HEADING
 VERIFICATION OF ABAQUS ELEMENT STATUS OUTPUT TO MATLAB (STATUS, RECORD KEY 61)
*PART,NAME=PART-1
*NODE
 1,  0.0, 0.0
 4,  1.0, 0.0
 13, 0.0, 1.0
 16, 1.0, 1.0
*NGEN, NSET=BOT
 1,4,1
*NGEN, NSET=TOP
 13,16,1
*NFILL, NSET=NALL
 BOT, TOP, 3, 4
*ELEMENT,TYPE=CPS4R
 1,1,2,6,5
*ELGEN, ELSET=EA
 1, 3, 1, 1, 3, 4, 3
*SOLID SECTION,ELSET=EA,MATERIAL=GLASS_EPOXY,ORIENT=RECT
 1.,
*NSET,NSET=LEFT
 1,5,9,13
*NSET, NSET=RIGHT
 4, 8, 12, 16
*ORIENTATION,NAME=RECT
 1.0, 0.0, 0.0, 0.0, 1.0, 0.0
 3,0.0
*END PART
*ASSEMBLY,NAME=ASSEMBLY-1
*INSTANCE,NAME=PART-1-1,PART=PART-1
*END INSTANCE
*END ASSEMBLY
*MATERIAL,NAME=GLASS_EPOXY
*ELASTIC,TYPE=LAMINA
 53.8E9,17.9E9,0.25,8.96E9,8.96E9,6.88E9
*DAMAGE INITIATION,CRITERION=HASHIN,ALPHA=0.0
 1034E6,1034E6,27.6E6,138E6,41.4E6,50E6
*DAMAGE EVOLUTION, TYPE=ENERGY, SOFTENING=LINEAR
 10E6,10E6,5E6,5E6
*DAMAGE STABILIZATION
 0.0001, 0.0001, 0.0001, 0.0001
*STEP,INC=200,NLGEOM
*STATIC
 0.001,0.01,,0.001
*BOUNDARY
 ASSEMBLY-1.PART-1-1.RIGHT,1,1,0.001
 ASSEMBLY-1.PART-1-1.LEFT,1,1
 ASSEMBLY-1.PART-1-1.BOT, 2,2
*FILE FORMAT,ASCII
*EL FILE
 STATUS
*END STEP
*STEP,INC=200,NLGEOM
*STATIC
 0.01,0.01,,0.001
*MODEL CHANGE, TYPE=ELEMENT, REMOVE
 PART-1-1.5,
*EL FILE
 STATUS
*END STEP
*STEP,INC=200,NLGEOM
*STATIC
 0.001,0.01,,0.001
*BOUNDARY
 ASSEMBLY-1.PART-1-1.RIGHT,1,1,0.1
*EL FILE
 STATUS
*END STEP
*STEP,INC=200,NLGEOM
*STATIC
 0.01,0.01,,0.001
*MODEL CHANGE, TYPE=ELEMENT, ADD=STRAIN FREE
 PART-1-1.5,
*EL FILE
 STATUS
*END STEP
*STEP,INC=200,NLGEOM
*STATIC
 0.001,0.01,,0.001
*BOUNDARY
 ASSEMBLY-1.PART-1-1.RIGHT,1,1,0.2
*EL FILE
 STATUS
*END STEP