液压成形模拟计算中断

wdjsc

我在dynaform中模拟液压成形过程(Dynaform的求解器是ls-dyna)，但不知怎么的，Ls-dyna分析计算时总是莫名其妙地中断，下面是我的message文件。大家看看是怎么回事？
      Date: 05/17/2003      Time: 22:10:18   
  
     ___________________________________________________
     |                                                 |
     |  Livermore  Software  Technology  Corporation   |
     |                                                 |
     |  7374 Las Positas Road                          |
     |  Livermore, CA 94550                            |
     |  Tel: (925) 449-2500  Fax: (925) 449-2507       |
     |  www.lstc.com                                   |
     |_________________________________________________|
     |                                                 |
     |  LS-DYNA, A Program for Nonlinear Dynamic       |
     |  Analysis of Structures in Three Dimensions     |
     |  Version: 950d         Date: 03/23/2000         |
     |                                                 |
     |  Platform   : Windows NT/98/95 (Intel)          |
     |  Precision  : Single precision (32 bits)        |
     |                                                 |
     |  Unauthorized use infringes LSTC copyrights     |
     |_________________________________________________|
  
  Input file: D:\OILBOX~1\INPUTD~1.DYN                                                
expanding   memory to  26214400
  
********  notice  ********  notice  ********  notice  ********
*                                                            *
* This is the LS-DYNA Explicit Finite Element code.          *
* Neither LSTC nor the authors assume any responsibility for *
* the validity, accuracy, or applicability of any results    *
* obtained from this system. The user must verify his own    *
* results.                                                   *
*                                                            *
* LSTC endeavors to make the LS-DYNA code as complete,       *
* accurate and easy to use as possible.                      *
* Suggestions and comments are welcomed.  Please report any  *
* errors encountered in either the documentation or results  *
* immediately to LSTC through your site focus.               *
*                                                            *
* Copyright (C) 1990-1999                                    *
* by Livermore Software Technology Corp. (LSTC)              *
* All rights reserved                                        *
*                                                            *
********  notice  ********  notice  ********  notice  ********
contracting memory to     24000
expanding   memory to     83540
expanding   memory to    173200
expanding   memory to    235496
expanding   memory to    259996
expanding   memory to    280033
expanding   memory to    302423
contracting memory to         2
  
  LS-DYNA will perform a structural only analysis  
  
1WWW                                                                     
                         ls-dyna (version 950d      )     date 03/23/2000
expanding   memory to     15851
contracting memory to     11754
expanding   memory to    150497
expanding   memory to    152032
expanding   memory to    218376
expanding   memory to    220422
expanding   memory to    222502
expanding   memory to    223906
expanding   memory to    226907
expanding   memory to    230214
expanding   memory to    230216
expanding   memory to    253957
expanding   memory to    254917
expanding   memory to    267462
expanding   memory to    324310
expanding   memory to    329429
expanding   memory to    375633
expanding   memory to    400525
contracting memory to    400525
contracting memory to    400525
expanding   memory to    400605
expanding   memory to    405375
expanding   memory to    408431
       0 t 0.0000E+00 dt 0.00E+00 flush i/o buffers
expanding   memory to    414651
  master surface has   373 segments,     0 renumbered
  
  slave surface has   732 segments,     0 renumbered
  
  connectivity of improperly numbered master segments
     8708     8707     8925     8929
     8925     8924     8927     8929
     8723     8695     8708     8929
     8927     8724     8723     8929
     8925     8707     8706     8928
  additional segments are being renumbered but not printed
  master surface has   325 segments,    72 renumbered
  
  slave surface has   732 segments,     0 renumbered
  
  master surface has   125 segments,     0 renumbered
  
  slave surface has   732 segments,     0 renumbered
  
  The LS-DYNA time step size should not exceed     0.159E-04
  to avoid contact instabilities.  If the step size is
  bigger then scale the penalty of the offending surface.
  
contracting memory to    408431
       0 t 0.0000E+00 dt 0.00E+00 flush i/o buffers
  
Slave  side interface #   1
# nodes on + segment side =   154 max dist = 3.15843E+01 min dist = 4.01054E-01
# nodes on - segment side =     0 max dist =-1.00000E+10 min dist = 1.00000E+10
  
Slave  side interface #   2
# nodes on + segment side =     0 max dist =-1.00000E+10 min dist = 1.00000E+10
# nodes on - segment side =   495 max dist = 3.55199E+01 min dist = 4.01031E-01
  
Slave  side interface #   3
# nodes on + segment side =   399 max dist = 4.01040E-01 min dist = 4.00995E-01
# nodes on - segment side =     0 max dist =-1.00000E+10 min dist = 1.00000E+10
  
  calculation with mass scaling for minimum dt
    added mass   =   8.6114E+01
    physical mass=   4.4982E-03
    ratio        =   1.9144E+04
  
       1 t 0.0000E+00 dt 5.05E-04 flush i/o buffers
       1 t 0.0000E+00 dt 5.05E-04 write d3plot file
  
  ******** termination due to mass increase *******
          added mass       =     8.6114E+01
          percent increase =     1.9144E+06
  
            nodes with largest mass increase
  
   node id      old mass        new mass        increase
   -------    ------------    ------------    ------------
      8243     1.07488E-06     1.57180E-01     1.57179E-01
      7375     1.07488E-06     1.57180E-01     1.57179E-01
      7943     1.07488E-06     1.57180E-01     1.57178E-01
      7695     1.07488E-06     1.57180E-01     1.57178E-01
      7667     1.07488E-06     1.57179E-01     1.57178E-01
      8235     1.07488E-06     1.57179E-01     1.57178E-01
      7915     1.07488E-06     1.57179E-01     1.57178E-01
      7367     1.07488E-06     1.57179E-01     1.57178E-01
      8186     1.02815E-06     1.52601E-01     1.52600E-01
      7538     1.02815E-06     1.52601E-01     1.52600E-01
      8242     1.25989E-06     1.52371E-01     1.52370E-01
      7912     1.25988E-06     1.52371E-01     1.52370E-01
      7726     1.25989E-06     1.52371E-01     1.52370E-01
      8236     1.25989E-06     1.52371E-01     1.52370E-01
      7374     1.25989E-06     1.52371E-01     1.52370E-01
      7698     1.25989E-06     1.52371E-01     1.52370E-01
      7884     1.25989E-06     1.52371E-01     1.52370E-01
      7368     1.25989E-06     1.52371E-01     1.52370E-01
      8072     9.81682E-07     1.43926E-01     1.43925E-01
      7424     9.81682E-07     1.43926E-01     1.43925E-01
   total added mass for these 20 nodes =  3.06944E+00
  
       1 t 5.0490E-04 dt 5.05E-04 write d3dump01 file
       1 t 5.0490E-04 dt 5.05E-04 write d3plot file
       1 t 5.0490E-04 dt 5.05E-04 write d3plot file
  
  E r r o r   t e r m i n a t i o n
  
  T i m i n g   i n f o r m a t i o n
                        CPU(seconds)   %CPU  Clock(seconds) %Clock
  ----------------------------------------------------------------
  Initialization ....... 2.0000E+00   50.00     1.9800E+00   50.77
  Element processing ... 1.0000E+00   25.00     8.3000E-01   21.28
  Binary databases ..... 0.0000E+00    0.00     5.0000E-02    1.28
  ASCII database ....... 0.0000E+00    0.00     0.0000E+00    0.00
  Contact algorithm .... 1.0000E+00   25.00     1.0400E+00   26.67
  Contact entities ..... 0.0000E+00    0.00     0.0000E+00    0.00
  Rigid bodies ......... 0.0000E+00    0.00     0.0000E+00    0.00
  ----------------------------------------------------------------
  T o t a l s            4.0000E+00  100.00     3.9000E+00  100.00
  
  Problem time       =    5.0490E-04
  Problem cycle      =         1
  Total CPU time     =         4 seconds (   0 hours  0 minutes  4 seconds)
  CPU time per zone cycle  =   1286173 nanoseconds
  Clock time per zone cycle=   1234726 nanoseconds
  Start time   05/17/2003 22:10:18   
  End time     05/17/2003 22:10:23   
  Elapsed time        5 seconds (   0 hours  0 minutes  5 seconds)

caeguy

時間步長巨大

belldingding

看着像是接触设置不对呀

wise

因为你的质量缩放太大造成的，减小*control_timestep中的DT2MS值。一般质量增加量要严格控制。

wdjsc

有劳各位了。能否再告诉我具体该怎么改，也就是说这些值该怎么确定？谢谢！

belldingding

修改K文件，最好自己试试，用二分缩减法。

wyons

added mass = 8.6114E+01  
          percent increase = 1.9144E+06  
  
以上信息说明程序停止是因为质量增加太大，达到 percent 1.9144E+06 。
  
The LS-DYNA time step size should not exceed 0.159E-04  
to avoid contact instabilities. If the step size is  
bigger then scale the penalty of the offending surface.  
  
这部分说明，根据你的网格和材料模型计算出来的最大时间步长为0.159E-04 。
  
解决的办法有三个：
１、不使用质量缩放，采用自动时间步长
２、质量缩放时，指定时间步长不要超过程序给的时间步长太多，并且使用负值，使质量缩放只限于特征长度较小的单元。
３、修改模型或者材料，使单元网格均匀，时间步长增大。