求救：VUMAT的input文件出错，可能是沙漏问题，好几天都搞不出来，谢谢各位

自己编的最基础的弹性VUMAT，可能是沙漏的设置问题，怎么都不对，请大家伙帮我看看，小弟先谢谢了。
只用一个单元，单元是C3D8R，但是INPUT文件总是有问题，提示是沙漏问题，却不知怎么加关键字*HOURGLASS STIFFNESS ，还有和其他关键字*SOLIDSECTION CONTROL 是怎么配合的。
DATA文件里是这样写的：
PROCESSING PART, INSTANCE, AND ASSEMBLY INFORMATION
*******************************************************
END PROCESSING PART, INSTANCE, AND ASSEMBLY INFORMATION
***********************************************************
OPTIONS BEING PROCESSED
***************************
*Heading
*Node
*Element, type=C3D8R
*Nset, nset=ASSEMBLY_PART-ELASTIC-1__PICKEDSET2
*Elset, elset=ASSEMBLY_PART-ELASTIC-1__PICKEDSET2
*Nset, nset=ASSEMBLY__PICKEDSET5
*Nset, nset=ASSEMBLY__PICKEDSET6
*Nset, nset=ASSEMBLY__PICKEDSET7
*Elset, elset=ASSEMBLY__PICKEDSET5
*Elset, elset=ASSEMBLY__PICKEDSET6
*Elset, elset=ASSEMBLY__PICKEDSET7
*Elset, elset=ASSEMBLY___PICKEDSURF4_S4
*surface, type=ELEMENT, name=ASSEMBLY__PICKEDSURF4
*surface, type=ELEMENT, name=ASSEMBLY__PICKEDSURF4
*material, name=ELASTIC
*density
*usermaterial, constants=2
*solidsection, elset=ASSEMBLY_PART-ELASTIC-1__PICKEDSET2, controls=EC-1, material=ELASTIC
*sectioncontrols, name=EC-1, hourglass=STIFFNESS
***WARNING: THE PARAMETER HOURGLASS = STIFFNESS ON THE SECTION CONTROLS OPTION
IS RELEVANT FOR THESE ELEMENTS: C3D8R, CAX4R, CGAX4R, CPEG4R,
CPE4R, CPS4R, M3D4R, M3D9R, S4R5, S4R, S8R5, S9R5, SC6R, SC8R AND
THEIR HYBRID, THERMAL AND PRESSURE COUNTERPARTS WHEREVER
APPLICABLE. IT IS ALSO RELEVANT FOR ALL TYPES OF MODIFIED
TRIANGULAR AND TETRAHEDRAL ELEMENTS. THIS WARNING CAN BE IGNORED
IF THE FEATURE IS APPLIED TO THESE ELEMENT TYPES ONLY.
*boundary
*boundary
*boundary
*solidsection, elset=ASSEMBLY_PART-ELASTIC-1__PICKEDSET2, controls=EC-1, material=ELASTIC
*hourglassstiffness
***WARNING: THE *HOURGLASS STIFFNESS OPTION IS RELEVANT FOR THESE ELEMENTS
ONLY: C3D4H, C3D8R, C3D8RH, EC3D8R, C3D10M, C3D10MH, CAX4R,
CAX4RH, CAX6M, CAX6MH, CPEG4R, CPEG4RH, CPE4R, CPE4RH, CPE6M,
CPE6MH, CPS4R, CPS6M, M3D4R, M3D9R, S3/S3R, S4/S4R, S4R5,S8R,
S8R5, S9R5, SC6R, SC8R, STRI3, AND SAXA ELEMENTS. FOR THESE
ELEMENTS ONLY THE DRILL SCALE FACTOR IS RELEVANT: STRI3, S3/S3R,
S4, S8R, AND SAXA ELEMENTS.
***WARNING: NULL VALUES SPECIFIED FOR HOURGLASS STIFFNESS HAVE BEEN IGNORED
*surface, type=ELEMENT, name=ASSEMBLY__PICKEDSURF4
*output, field, variable=PRESELECT
*output, history, variable=PRESELECT
*output, field, variable=PRESELECT
*output, history, variable=PRESELECT
*output, field, variable=PRESELECT
*output, history, variable=PRESELECT
***ERROR: 1 elements have been defined with zero hour glass stiffness. You may
use *hourglass stiffness or change the element type. The elements
have been identified in element set ErrElemZeroHourGlassStiffness.
***NOTE: DUE TO AN INPUT ERROR THE ANALYSIS PRE-PROCESSOR HAS BEEN UNABLE TO
INTERPRET SOME DATA. SUBSEQUENT ERRORS MAY BE CAUSED BY THIS OMISSION
*Step, name=Step-compression
*output, field, variable=PRESELECT
*output, history, variable=PRESELECT
*Step, name=Step-compression
*Step, name=Step-compression
*static
*boundary
*boundary
*boundary
*dsload
*output, field, variable=PRESELECT
*output, history, variable=PRESELECT
*endstep
*Step, name=Step-compression
*static
*boundary
*boundary
*boundary
*output, field, variable=PRESELECT
*output, history, variable=PRESELECT
*endstep
***WARNING: USER SUBROUTINE UMAT WILL BE USED WITH THE STAVEV ARRAY
DIMENSIONED TO ZERO SINCE THE *DEPVAR OPTION IS NOT USED WITH THIS
MATERIAL. CONSEQUENTLY, DEFINING STATEV ENTRIES IN SUBROUTINE UMAT
WILL CAUSE CODE EXECUTION ERRORS.
- (RAMP) OR (STEP) - INDICATE USE OF DEFAULT AMPLITUDES ASSOCIATED WITH THE STEP
以下是我的INPUT文件：
*Heading
** Job name: Job-elastic Model name: Model-1
** Generated by: Abaqus/CAE Version 6.8-1
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=part-elastic
*Node
1, -5., -5., 10.
2, -5., 5., 10.
3, -5., -5., 0.
4, -5., 5., 0.
5, 5., -5., 10.
6, 5., 5., 10.
7, 5., -5., 0.
8, 5., 5., 0.
*Element, type=C3D8R
1, 5, 6, 8, 7, 1, 2, 4, 3
*Nset, nset=_PickedSet2, internal, generate
1, 8, 1
*Elset, elset=_PickedSet2, internal
1,
** Section: Section-elastic
*Solid Section, elset=_PickedSet2, controls=EC-1, material=elastic
*hourglass stiffness
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=part-elastic-1, part=part-elastic
*End Instance
**
*Nset, nset=_PickedSet5, internal, instance=part-elastic-1, generate
1, 7, 2
*Elset, elset=_PickedSet5, internal, instance=part-elastic-1
1,
*Nset, nset=_PickedSet6, internal, instance=part-elastic-1, generate
5, 8, 1
*Elset, elset=_PickedSet6, internal, instance=part-elastic-1
1,
*Nset, nset=_PickedSet7, internal, instance=part-elastic-1
1, 2, 5, 6
*Elset, elset=_PickedSet7, internal, instance=part-elastic-1
1,
*Elset, elset=__PickedSurf4_S4, internal, instance=part-elastic-1
1,
*Surface, type=ELEMENT, name=_PickedSurf4, internal
__PickedSurf4_S4, S4
*End Assembly
**
** ELEMENT CONTROLS
**
*Section Controls, name=EC-1, hourglass=STIFFNESS
1., 1., 1.
**
** MATERIALS
**
*Material, name=elastic
*density
2.698E-9
*User Material, constants=2
2.1e+11, 0.33
** ----------------------------------------------------------------
**
** STEP: Step-compression
**
*Step, name=Step-compression
*Static
1., 1., 1e-05, 1.
**
** BOUNDARY CONDITIONS
**
** Name: BC-X Type: Displacement/Rotation
*Boundary
_PickedSet6, 1, 1
** Name: BC-Y Type: Displacement/Rotation
*Boundary
_PickedSet5, 2, 2
** Name: BC-Z Type: Displacement/Rotation
*Boundary
_PickedSet7, 3, 3
**
** LOADS
**
** Name: Load-compression Type: Pressure
*Dsload
_PickedSurf4, P, 1e+10
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step