系列2-原創#cohesive+隱式分析+通用接觸--報錯node numbering原因

隱式分析通用接觸的節點被自動調整導致計算提交出錯.jpg
隱式分析通用接觸的節點被自動調整.jpg
有無通用接觸的pre文件差別--隱式分析通用接觸會自動調整節點.jpg
有無通用接觸的inp文件內容對比.jpg

********************************************************************************************

首先聲明(做了模型分析):

批量插入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單元,它就會調整節點),如果

    有找到徹底解決辦法的,可以相互探討,

********************************************************************************************

ABAQUS斷裂模擬收徒 ,保證快速學會各種ABAQUS斷裂模擬方法    1200/人(將享有各種插件以及程序,價值3000+、專門定制視頻、全程親自教學、各種模型調試及解答問題等等,傾囊相教)

登錄后免費查看全文
立即登錄
App下載
技術鄰APP
工程師必備
  • 項目客服
  • 培訓客服
  • 平臺客服

TOP

16
5
6