VUMAT子程序计算报错

CADTJ

笔者采用自编的梁单元单轴本构模型计算轴向受压梁，计算中出现如下报错：
Abaqus/Explicit 6.11-1                           DATE 05-Jul-2012  TIME 17:18:37
-------------------------------------------------------------------------------
PREPROCESSOR WARNING MESSAGES
-------------------------------------------------------------------------------

***WARNING: There are 1 warning messages in the data (.dat) file.  Please
            check the data file for possible errors in the input file.

-------------------------------------------------------------------------------
  USER SUBROUTINE VUMAT
-------------------------------------------------------------------------------

***WARNING: A material defined in user subroutine VUMAT must be defined as
            purely elastic (using the initial elastic modulus) at the
            beginning of the analysis (stepTime=0). This is an informative
            message. It does not necessarily indicate that user subroutine
            VUMAT is incorrectly defined.

-------------------------------------------------------------------------------
MODEL INFORMATION (IN GLOBAL X-Y COORDINATES)
-------------------------------------------------------------------------------
   Total mass (unscaled) in model = 540.00
   Center of mass of model (based on unscaled mass) = ( 0.000000E+00, 6.750000E-01, 0.000000E+00)
    Moments of Inertia :
                 About Center of Mass              About Origin
      I(XX)          9.085275E+01                  3.368903E+02
      I(YY)          1.440000E+01                  1.440000E+01
      I(ZZ)          9.085275E+01                  3.368903E+02
      I(XY)          0.000000E+00                  0.000000E+00
      I(YZ)          0.000000E+00                  0.000000E+00
      I(ZX)         -2.286989E-17                 -2.286989E-17
-------------------------------------------------------------------------------
STABLE TIME INCREMENT INFORMATION
-------------------------------------------------------------------------------

  The stable time increment estimate for each element is based on
  linearization about the initial state. The stable time increment estimates
  do not account for mass scaling.

   Initial time increment = 3.55067E-05
   Statistics for all elements:
      Mean = 3.55067E-05
      Standard deviation = 0.0000
   Most critical elements:
    Element number   Rank    Time increment   Increment ratio
    (Instance name)
   ----------------------------------------------------------
           1          1       3.550669E-05      1.000000E+00
    BEAM-1
           3          2       3.550669E-05      1.000000E+00
    BEAM-1
           5          3       3.550669E-05      1.000000E+00
    BEAM-1
           6          4       3.550669E-05      1.000000E+00
    BEAM-1
           7          5       3.550669E-05      1.000000E+00
    BEAM-1
           9          6       3.550669E-05      1.000000E+00
    BEAM-1
          10          7       3.550669E-05      1.000000E+00
    BEAM-1
           2          8       3.550670E-05      9.999997E-01
    BEAM-1
           4          9       3.550670E-05      9.999997E-01
    BEAM-1
           8         10       3.550671E-05      9.999995E-01
    BEAM-1

  Instance name for the most critical element will be printed if the critical
  element changes or at the beginning of a new page.

Abaqus/Explicit 6.11-1                           DATE 05-Jul-2012  TIME 17:18:38
The double precision Abaqus/Explicit executable will be used in this analysis.

-------------------------------------------------------------------------------
SOLUTION PROGRESS
-------------------------------------------------------------------------------
STEP 1  ORIGIN 0.0000
  Total memory used for step 1 is approximately 368.9 kilobytes.
  Direct time incrementation control will be used.
  Variable mass scaling factor at zero increment: 1.0000
***WARNING: The direct user control is used in the current analysis step. In
            this case, Abaqus/Explicit will not check that the computed
            response is stable during the step. Please check carefully the
            energy history and other response variables to ensure that a valid
            and stable response has been obtained.
              STEP     TOTAL       CPU      STABLE    CRITICAL    KINETIC      TOTAL    PERCENT
INCREMENT     TIME      TIME      TIME   INCREMENT     ELEMENT     ENERGY     ENERGY  CHNG MASS
        0  0.000E+00 0.000E+00  00:00:00 1.000E-05          10  3.037E-03  3.037E-03   0.000E+00
INSTANCE WITH CRITICAL ELEMENT: BEAM-1
ODB Field Frame Number      0 of   1000 requested intervals at increment zero.
***ERROR: Floating Point Exception detected in increment 200.  Job exiting.
      200  2.000E-03 2.000E-03  00:00:00 1.000E-05           0 NaN        NaN          0.000E+00
ODB Field Frame Number      1 of   1000 requested intervals at  2.000010E-03
Nodes involved in floating point error:

The following checklist may be helpful in diagnosing the error:
    1. Check contact definitions for problems such as excessive initial
       overclosure or unrealistic tied definition between contact pairs.
       A vector plot of velocities or accelerations will usually help to
       identify contact problems.
    2. Check stiffness (elastic modulus) and mass (density) definitions for
       consistent units and verify that the combination is reasonable.
    3. Check for poor mesh definition.
    4. Check the boundary conditions for an excessive loading rate. The
       *DIAGNOSTICS, DEFORMATION SPEED CHECK=DETAIL option may be used to obtain
       detailed diagnostics information.
    5. Check the current status of the structure to see if it has
       totally failed.
    6. A dashpot or a very stiff spring may cause the analysis to go unstable.
       The *DYNAMIC, DIRECT option may be used to control the time increment
       directly.

***ERROR: The fatal error will generate a new field output in the last
          increment. All variables applicable to the current procedure and
          material type(s) will be written to the output database.

Nodes involved in floating point error:

***ERROR: Floating Point Exception detected in increment 202.  Job exiting.
      202  2.010E-03 2.010E-03  00:00:00 1.000E-05           0 NaN        NaN          0.000E+00
A special ODB Field Frame with all relevant variables was written at     2.010E-03
  THE ANALYSIS HAS NOT BEEN COMPLETED

对于错误警告信息中Floating Point Exception这是什么意思？有什么可以解决的方法吗？