原创系列2：Cohesive隐式分析与通用接触报错解析

首先声明（做了模型分析）：

批量插入cohesive单元+显示分析              可以顺利提交计算；

批量插入cohesive单元+显示分析+通用接触     可以顺利提交计算；

批量插入cohesive单元+隐式分析              可以顺利提交计算；

批量插入cohesive单元+隐式分析+通用接触   不可以顺利提交计算；

报错信息：

批量插入cohesive单元+隐式分析+通用接触，模型提交分析时，出现报错，如下：

The node numbering might not be correct for element 8133 instance part-1-1

The node numbering might not be correct for element 16211 instance part-1-1

Analysis Input File Processor exited with an error.

警告信息：

Not all the nodes that have been adjusted were printed. Specify *preprint,contact=yes for complete printout.

Not all the nodes that have been adjusted were printed. Specify *preprint,contact=yes for complete printout.

15416 nodes have been adjusted more than once. The subsequent adjustments may cause these nodes not to lie on their master surface. The nodes have been identified in node set WarnNodeAdjust.

MPCS (EXTERNAL or INTERNAL, including those generated from rigid body definitions), KINEMATIC COUPLINGS, AND/OR EQUATIONS WILL ACTIVATE ADDITIONAL DEGREES OF FREEDOM

直接aborted，不能计算

问题分析：

1 既然报错inp文件，那么我们检查了批量插入cohesive单元+隐式分析在有无通用接触下的inp文件

  发现除了通用接触没有设置外，其他的完全相同，那么，我们现在可以确定是由于通用接触引起的

  报错而不能进行提交计算（做了大量的模型对比）。

2 那么，既然知道是由于通用接触的加入导致的模型不能计算，那么我们就需要分析为什么通用接触的

  加入会导致计算不能提交呢，另外，现在除了报错节点编号错误外，没有其他错误抛出，那么我们首先

  检查aba抛出的节点编号错误，经过仔细检查分析inp文件中报错的单元（报错的单元是cohesive单元），

  分析结果表明：节点编号并没有问题（编号正确，而且是逆时针的），而且，对于含通用接触的inp（不

  能计算）和不含通用接触的inp（能计算）两者的检验结果是一样的。

3  既然检查了abaqus抛出的节点编号问题后发现，节点编号是没有问题的，那么接下来怎么办，...，陷入

  死循环了，研究了好多天不能找到问题所在...

4  等等，除了这个抛出的错误外，还有警告呢，警告有两个，第二个没什么用，第一个很重要，它说有15416个

  节点被调整了，后续的调整可能引起这些及诶单不在他们的主面上，这些节点已经被定义在了节点集合WarnNodeAdjust中，

  那么我们看一下这个调整的节点集合吧，进入后处理，发现这个集合就是我们插入cohesive区域内的一些节点，

  也就是说在使用批量插入cohesive单元+隐式分析+通用接触时，节点被abaqus自动进行了调整，然后导致了不能计算。

5  因为inp文件根本不能正确提交计算，所以我们只能看.pre前处理文件了，对比有无通用接触的两个前处理文件发现，

  包含通用接触的pre文件中多了好多下面的信息（打开inp文件，修改前面的 *Preprint, echo=NO, model=NO, history=NO, contact=YES）：

  TRACKING LOAD BALANCE

      General_Contact_Faces / General_Contact_Faces

     PROCESSOR 1:  5784 FACETS

     PROCESSOR 2:  4864 FACETS

     PROCESSOR 3:  5632 FACETS

     PROCESSOR 4:  5376 FACETS

     PROCESSOR 5:  5632 FACETS

     PROCESSOR 6:  5632 FACETS

 SLAVE : General_Contact_Faces                                                          

 MASTER: General_Contact_Faces                                                          

             DISTANCE

  NODE       ADJUSTED    REMARKS/WARNINGS

        3   1.73472E-18  ADJUSTED TO ( 8.86807E-03, 1.70870E-02,  0.0000    ).

        4    0.0000      ADJUSTED TO ( 8.58252E-03, 1.68079E-02,  0.0000    ).

        6    0.0000      ADJUSTED TO ( 9.13970E-03, 1.67108E-02,  0.0000    ).

        7    0.0000      ADJUSTED TO ( 8.86807E-03, 1.70870E-02, 1.00000E-03).

        8    0.0000      ADJUSTED TO ( 8.58252E-03, 1.68079E-02, 1.00000E-03).

       10   1.73472E-18  ADJUSTED TO ( 9.13970E-03, 1.67108E-02, 1.00000E-03).

.

.

.

.

.

.

--------------------------------------------------------------------------------

Elements with large aspect ratio

     Element       Adjusted nodes

------------------ --------------

     PART-1-1.4466       NO      

     PART-1-1.4825       NO      

     PART-1-1.4826       NO      

     PART-1-1.4827       NO      

     PART-1-1.4829       NO      

     PART-1-1.4830       NO      

     PART-1-1.5021       NO      

     PART-1-1.5084       NO      

     PART-1-1.5137       NO      

     PART-1-1.5179       NO      

     PART-1-1.5311       NO      

     PART-1-1.5469       NO      

     PART-1-1.5471       NO      

     PART-1-1.5472       NO      

     PART-1-1.5488       NO      

.

.

.

    LABELS REFERENCED IN THE ABOVE TABLE

     -19-     INTPROP-1                                                                      

     -20-     General_Contact_Faces                                                          

                                I N I T I A L   C O N T A C T   S T A T U S

***NOTE: The following initial contact status information will reflect some

         but not all types of automated nodal position adjustments. Output

         variables such as CSTATUS and COPEN written to the output database

         (.odb) at the initial output frame (if the CSTATUS and CDISP output

         variables are requested) will reflect all types of automated nodal

         position adjustments. Viewing CSTATUS and COPEN in a post-processor

         is the recommended way to determine the initial contact status.

         Initial contact status information may not be available in the .dat

         file in a future version.

   SLAVE SURFACE:         General_Contact_Faces                                                          

   MASTER SURFACE:        General_Contact_Faces                                                          

      NODE        STATUS            CLEARANCE            NOTE

      NUMBER

         3       CLOSED            6.9228E-19                        

         4       CLOSED            6.9228E-19                        

         5        OPEN             1.9240E-03                        

         6       CLOSED             0.000                            

         7       CLOSED             0.000                            

         8       CLOSED            1.1027E-18                        

         9        OPEN             1.9240E-03                        

        10       CLOSED             0.000                            

        11       CLOSED             0.000                            

        12       CLOSED             0.000                            

        13       CLOSED             0.000                            

        14       CLOSED             0.000                            

        15       CLOSED            5.7674E-19      

  好吧，这应该就是症结所在了，我们有分析了隐式分析有无通用接触情况下的pre文件，发现只有批量插入cohesive单元+隐式分析+通用接触

  的时候会出现节点被调整的问题，总的来说就是因为abaqus自动调整了节点，导致了我们的inp文件不能计算，而出现的节点编号错误的提示。

  因为之后计算的时候才会自动调整这些节点，所以有无通用接触的两种情况下输出的inp都是相同的，因此只分析inp文件我们并不能找到根本

  原因。

解决办法：

 1 删除报错的单元（治标不治本，而且如果报错有很多单元怎么办呢？）

 2 不使用通用接触（但是如果是模型中开裂后的裂缝又在压应力作用下闭合了怎么办？）

 3 把模型的所有外表面建立为一个surface集合，通用接触处使用selected surface pair，可以成功提交计算，但是计算特别慢，而且计算结果总

   感觉不对，大家可以尝试。

 4 彻底的解决办法暂时没找到（或者根本没有，因为批量插入cohesive单元，隐式分析只要使用批量插入cohesive单元，它就会调整节点），如果

   有找到彻底解决办法的，可以相互探讨，

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删