请教一个关于碰撞的问题

ligezhan

各位大哥，我由于做球穿透板的那个例子老是不对，从论坛上找到启示直接运行官方例子的dat文件，但还是出现同样的错误，百思不得其解，不知道是软件的原因还是例子本身就有错误，特来向各位大哥请教
这是求穿透板的错误文件：
DATE  :30-04-05  TIME  :13:21:28                                                
                                                                                
                                                                                
                                                                                
%I-P2044902-P2XXTXIX_PACKET_ZVARS,,,                                          
        Requested variable EFFPL is unknown and changed into EFFPLS            
                                                                                
%I-P2044902-P2XXTXIX_PACKET_ZVARS,,,                                          
        Requested variable EFFST is unknown and changed into EFFSTS            
                                                                                
%W-P2108906-P2_DEFINE_VARIABLES_FOR_EDIT,,,                                    
        Variable EFFPLS is not legal in input set PLATE and is                  
        therefore ignored.                                                      
        This may result in an empty edit set.                                   
                                                                                
%W-P2108906-P2_DEFINE_VARIABLES_FOR_EDIT,,,                                    
        Variable EFFSTS is not legal in input set PLATE and is                  
        therefore ignored.                                                      
        This may result in an empty edit set.                                   
                                                                                
%W-P2108906-P2_DEFINE_VARIABLES_FOR_EDIT,,,                                    
        Variable TXX is not legal in input set PLATE and is                     
        therefore ignored.                                                      
        This may result in an empty edit set.                                   
                                                                                
%I-P4142610-P4_CONT_OSCNTC_V4,,,                                               
        The suggested number of cubes from the contact cube sorting algoritm is
        5400. LIMCUB is set to 2000                                             
        The contact runs faster if LIMCUB is raised to at least the number      
        that is suggested by the algoritm. This can be achieved by setting      
        PARAM,LIMCUB to 5400 at the beginning of a new analysis.               
        Please note that more memory will be used and it is not allowed in a   
        restart analysis.                                                      
                                                                                
%I-P4142610-P4_CONT_OSCNTC_V4,,,                                               
        The suggested number of cubes from the contact cube sorting algoritm is
        3645. LIMCUB is set to 2000                                             
        The contact runs faster if LIMCUB is raised to at least the number      
        that is suggested by the algoritm. This can be achieved by setting      
        PARAM,LIMCUB to 3645 at the beginning of a new analysis.               
        Please note that more memory will be used and it is not allowed in a   
        restart analysis.                                                      
                                                                                
%I-P4142610-P4_CONT_OSCNTC_V4,,,                                               
        The suggested number of cubes from the contact cube sorting algoritm is
        3125. LIMCUB is set to 2000                                             
        The contact runs faster if LIMCUB is raised to at least the number      
        that is suggested by the algoritm. This can be achieved by setting      
        PARAM,LIMCUB to 3125 at the beginning of a new analysis.               
        Please note that more memory will be used and it is not allowed in a   
        restart analysis.                                                      
                                                                                
%I-P4142610-P4_CONT_OSCNTC_V4,,,                                               
        The suggested number of cubes from the contact cube sorting algoritm is
        2304. LIMCUB is set to 2000                                             
        The contact runs faster if LIMCUB is raised to at least the number      
        that is suggested by the algoritm. This can be achieved by setting      
        PARAM,LIMCUB to 2304 at the beginning of a new analysis.               
        Please note that more memory will be used and it is not allowed in a   
        restart analysis.                                                      
                                                                                
%I-P4142610-P4_CONT_OSCNTC_V4,,,                                               
        The suggested number of cubes from the contact cube sorting algoritm is
        2304. LIMCUB is set to 2000                                             
        The contact runs faster if LIMCUB is raised to at least the number      
        that is suggested by the algoritm. This can be achieved by setting      
        PARAM,LIMCUB to 2304 at the beginning of a new analysis.               
        Please note that more memory will be used and it is not allowed in a   
        restart analysis.                                                      
NO MORE MESSAGES OF THIS KIND WILL BE WRITTEN                                 
                                                                                
                                                                                
                                                                                
JOBNAME IS: RIGID_IMAPCT                                                        
----------------------------------------------------                           
MSC.Dytran 2004 Single-CPU Windows                                             
SUMMARY OF MESSAGES (INCLUDING NON-PRINTED MESSAGES)                           
INFORMATIVE         :      7                                                   
WARNINGS            :      3                                                   
NO ERRORS                                                                       

这是梁撞击墙的错误文件：
DATE  :30-04-05  TIME  :13:57:15                                                
                                                                                
                                                                                
                                                                                
%W-P2108906-P2_DEFINE_VARIABLES_FOR_EDIT,,,                                    
        Variable EFFST-IN is not legal in input set BEAM and is                 
        therefore ignored.                                                      
        This may result in an empty edit set.                                   
                                                                                
%W-P2108906-P2_DEFINE_VARIABLES_FOR_EDIT,,,                                    
        Variable EFFST-MID is not legal in input set BEAM and is               
        therefore ignored.                                                      
        This may result in an empty edit set.                                   
                                                                                
%W-P2108906-P2_DEFINE_VARIABLES_FOR_EDIT,,,                                    
        Variable EFFST-OUT is not legal in input set BEAM and is               
        therefore ignored.                                                      
        This may result in an empty edit set.                                   
                                                                                
%W-P2108906-P2_DEFINE_VARIABLES_FOR_EDIT,,,                                    
        Variable EFFPL-IN is not legal in input set BEAM and is                 
        therefore ignored.                                                      
        This may result in an empty edit set.                                   
                                                                                
%W-P2108906-P2_DEFINE_VARIABLES_FOR_EDIT,,,                                    
        Variable EFFPL-MID is not legal in input set BEAM and is               
        therefore ignored.                                                      
        This may result in an empty edit set.                                   
NO MORE MESSAGES OF THIS KIND WILL BE WRITTEN                                 
                                                                                
%W-P4216001-V4LZT8I1_MOD2_STRAIN_RATE_SPIN,,,                                 
        Time step 1.097E-06 too large for correct stress rotation               
        Use corotational stress formulation                                    
                                                                                
                                                                                
                                                                                
JOBNAME IS: BEAM                                                               
----------------------------------------------------                           
MSC.Dytran 2004 Single-CPU Windows                                             
SUMMARY OF MESSAGES (INCLUDING NON-PRINTED MESSAGES)                           
WARNINGS            :      6                                                   
NO ERRORS

Lcl008

这些是警告和信息.
%%E,%%F才是错误.看来你还要多看书.

ligezhan

Lcd008兄，多谢您的回帖，我知道这是警告，但正是由于它们才没有结果啊，我应该怎么修改呢

worker

Lcl008兄说得对，dytran例子也有同样警告出现但等效应力应变结果还是有的。
ligezhan兄，你的dat文件中有几个地方跟例子不太一样。
ELOUT (plate) = EFFPL-OUT EFFPL-IN EFFST-OUT EFFST-IN FAIL-OUT FAIL-IN
例子中是ELOUT (Plate) = EFFPL-MID,EFFST-MID
因为板单元默认都是三点高斯积分，每个积分点的应力应变结果可以通过积分点输出。你只输出了上下两点，中间点没有。

另外你的刚性球单元是一般板单元的属性
$            * rigid_ball *
$
PSHELL         1       1       1
例子中为DUMMY单元
PSHELL1,200, ,DUMMY
建议你改为
PSHELL1,1, ,DUMMY

再有你的主从接触面使用属性PROP定义的
CONTACT        1    PROP    PROP     1       2                          +A000003
+A000003      V4    BOTH           YES       1.0                        +A000004
+A000004                                                                +A000005
+A000005                            ON
这样你的两个SURFACE就没用了
$ Slave contact surface for contact
$
SURFACE        1            ELEM       3
SET1           3    1568    THRU    2016     987    THRU    1567     406+A000004
+A000004    THRU     986       1    THRU     405
$
$ Master contact surface for contact
$
SURFACE        2            ELEM       4
SET1           4    4341    THRU    4516    3760    THRU    4340    3179+A000005
+A000005    THRU    3759    2598    THRU    3178    2017    THRU    2597
$
建议你要么改用CONTACT,1,SURF,SURF,1,2定义接触，要么把两个SURFACE删掉。
其他的信息
%I-P4142610-P4_CONT_OSCNTC_V4,,,                                               
        The suggested number of cubes from the contact cube sorting algoritm is
        22090. LIMCUB is set to 2000                                            
等一般没啥大问题。还有就是你没有定义节点分离BJOIN这样恐怕作不出例子中破片飞散的效果。如果你不关心这个也无所谓。
再算一下试试还有问题就把dat文件贴来让大家看看。

ligezhan

好的，多谢worker兄和Lcd008兄的指点，我试试看