*preprint, model=yes, history=yes
*restart, write
**
*node, nset=nodes_tfv0
 1001, -1.,-1.,-1.
 1002,  1.,-1.,-1.
 1003,  1., 1.,-1.
 1004, -1., 1.,-1.
 1005, -1.,-1., 1.
 1006,  1.,-1., 1.
 1007,  1., 1., 1.
 1008, -1., 1., 1.
 1011, -1.,-1.,-2.
 1012,  1.,-1.,-2.
 1013,  1., 1.,-2.
 1014, -1., 1.,-2.
 1015, -1.,-1.,-1.
 1016,  1.,-1.,-1.
 1017,  1., 1.,-1.
 1018, -1., 1.,-1.
 1021, -1.,-1.,-4.
 1022,  1.,-1.,-4.
 1023,  1., 1.,-4.
 1024, -1., 1.,-4.
 1025, -1.,-1.,-2.
 1026,  1.,-1.,-2.
 1027,  1., 1.,-2.
 1028, -1., 1.,-2.
 1031, -1.,-1.,-5.
 1032,  1.,-1.,-5.
 1033,  1., 1.,-5.
 1034, -1., 1.,-5.
 1035, -1.,-1.,-4.
 1036,  1.,-1.,-4.
 1037,  1., 1.,-4.
 1038, -1., 1.,-4.
**
*element, type=c3d8r, elset=cube_tfv0
 121, 1021,1022,1023,1024,1011,1012,1013,1014
 111, 1011,1012,1013,1014,1001,1002,1003,1004
 101, 1001,1002,1003,1004,1005,1006,1007,1008
 131, 1031,1032,1033,1034, 1021,1022,1023,1024
**
*solid section, material=mat_tfv0, elset=cube_tfv0
 1.
**
*material, name=mat_tfv0
*density
 0.01
*elastic
 1.0e3, 0.30
**
*elset, elset=ehist_tfv0
 101, 131
**
*nset, nset=nhist_tfv0
 1001, 1011, 1021, 1031
*nset, nset=nfixed_tfv0
 1031, 1032, 1033, 1034
*nset, nset=ndisp_tfv0
 1005, 1006, 1007, 1008
**
*boundary
 nfixed_tfv0, 1, 3
**
*Amplitude, Name=ramp1
 0., 0., 1., 1.
**
** Step 1
**
*step, nlgeom=YES, inc=100 
*static
 1., 1.
**
*boundary, amp=ramp1
 ndisp_tfv0, 3, 3, -0.001
**
*output, field, number interval=4
*element output
 mises
*node output
 u, rf
*output, history, frequency=1
*node output, nset=nhist_tfv0
 u
*element output, elset=ehist_tfv0
 mises
*end step