这种停止计算算正常吗

lij27

我做了个球自由落体的计算，能够实现与地面的两次碰撞（大概1秒钟），而且感觉很好，在我还想接着算下去看看后面的运动情况（与地面再次碰撞）的时候，计算停止了。但是我的最小时间步长已经设得很小了，并且模型并未发生什么变形，计算也未达到我所设的终止步，终止时间没有设，我曾试着在restart中把最小时间步长再设小些，还是不行，请问这可能是什么原因造成的，属于正常情况吗？dytran中除上以上一些准则外，还有什么情况会终止计算。谢谢

worker

“计算也未达到我所设的终止步”说明属于非正常停止计算。检查一下ERROR文件看看有没有出错提示。
最有可能的情况是计算步长太小，超过设定的最小时间步长。
在拉格朗日有限元计算中，时间步长变小大都因为单元畸变引起应力波在一个单元内传播时间过短造成。但根据你所说“模型并未发生什么变形”与此矛盾。
所以最好贴出你的输入文件，或出错信息。

lij27

以下就是我的输入文件，球能发生三次碰撞，就停止计算了，请指正。
START
CEND
ENDSTEP=300000
CHECK=NO
TITLE= Jobname is: 01
TLOAD=1
TIC=1
SPC=1
$ Output result for request: ball
TYPE (ball) = ARCHIVE
ELEMENTS (ball) = 1
SET 1 = 1 THRU 1636
ELOUT (ball) = AREA
STEPS (ball) = 0 THRU END BY 100
SAVE (ball) = 10000
$ Output result for request: rst
TYPE (rst) = RESTART
STEPS (rst) = 0 THRU END BY 30000
SAVE (rst) = 1
$------- Parameter Section ------
PARAM,CONTACT,THICK,0.0
PARAM,INISTEP,5.0e-6
PARAM,MINSTEP,1.0e-14
$------- BULK DATA SECTION -------
BEGIN BULK
INCLUDE 05.bdf
$ ------- GRAVITATION -----
TLOAD1         1     444               0
GRAV         444              10       0      -1       0
$
$ ========== PROPERTY SETS ==========
$
$            * ball *
$
PSHELL         1       1    .003
$
$            * rigid *
$
PSHELL         2       2    .005
$
$
$ ========= MATERIAL DEFINITIONS ==========
$
$
$ -------- Material ball id =1
MAT1           1   5e+08              .3    3000
$
$ -------- Material rigid id =2
MATRIG         2    7850                        
$
$ ======== Load Cases ========================
$
$
$ ------- Rigid Body: bodyfor
$
BODYFOR        3LAGRANGE  SHQUAD                                        +      
+              0CONSTANT     9.8       0      -1       0
$
$ -------- Contact : contace
$
CONTACT        4    ELEM    ELEM       2       3                        +      
+             V4    BOTH             YES                                
$
$ Adaptive Slave Contact Set for contace
$
SET1           2    1062    THRU    1536     481    THRU    1061       1+      
+           THRU     480
$
$ Adaptive Master Contact Set for contace
$
SET1           3    1537    THRU    1636
$
$ ------- Rigid Body Object rigid_constraint -----
$ ---- Reference Node for Rigid body is 1720
TLOAD1         1       5              12
FORCE          5     MR2               1       0       0       0
TLOAD1         1    1005              12
MOMENT      1005     MR2               1       0       0       0
$
$
ENDDATA

lij27

另外出错信息也没有什么提示呀，就是以下这些
%W-P2082603-P2_CHECK_DMATEP_FAILURE,,,                                         
        Shell elements with elastic material 1 do not have the                  
        strains and effective stress available as extended output.              
        Only shell elements with elasto-plastic materials do have               
        access to the additional data for output purposes.                     
                                                                                
%I-P3044801-P3_WRITE_RBE2_SUMMARY,,,                                          
        For a comlete MATRIG-RBE2 summary including a list of grid points for   
        each assembly, please use PARAM,RBE2INFO,GRIDON

lij27

再一次认真看了看变形的单元，发现确实有几个单元发生了较大的变形，问题极有可能就出在这了，之前没有检查出来。但是有一点想不太通，就是为什么第一次碰撞的速度是最快的，都没有产生大的变形，上升时应该有一个恢复的过程，之后碰撞的速度应该越来越小变形也应该越来越小才对。会不会是材料定义有问题。因为这是我随便想到的一个例子，材料也是随便定义的。