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

隱式分析通用接觸的節(jié)點(diǎn)被自動調(diào)整導(dǎo)致計(jì)算提交出錯.jpg
隱式分析通用接觸的節(jié)點(diǎn)被自動調(diào)整.jpg
有無通用接觸的pre文件差別--隱式分析通用接觸會自動調(diào)整節(jié)點(diǎn).jpg
有無通用接觸的inp文件內(nèi)容對比.jpg

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

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

批量插入cohesive單元+顯示分析              可以順利提交計(jì)算;

批量插入cohesive單元+顯示分析+通用接觸     可以順利提交計(jì)算;

批量插入cohesive單元+隱式分析              可以順利提交計(jì)算;

批量插入cohesive單元+隱式分析+通用接觸   不可以順利提交計(jì)算;

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

報錯信息:

批量插入cohesive單元+隱式分析+通用接觸,模型提交分析時,出現(xiàn)報錯,如下:

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,不能計(jì)算

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

問題分析:

 1 既然報錯inp文件,那么我們檢查了批量插入cohesive單元+隱式分析在有無通用接觸下的inp文件

   發(fā)現(xiàn)除了通用接觸沒有設(shè)置外,其他的完全相同,那么,我們現(xiàn)在可以確定是由于通用接觸引起的

   報錯而不能進(jìn)行提交計(jì)算(做了大量的模型對比)。

 2 那么,既然知道是由于通用接觸的加入導(dǎo)致的模型不能計(jì)算,那么我們就需要分析為什么通用接觸的

   加入會導(dǎo)致計(jì)算不能提交呢,另外,現(xiàn)在除了報錯節(jié)點(diǎn)編號錯誤外,沒有其他錯誤拋出,那么我們首先

   檢查aba拋出的節(jié)點(diǎn)編號錯誤,經(jīng)過仔細(xì)檢查分析inp文件中報錯的單元(報錯的單元是cohesive單元),

   分析結(jié)果表明:節(jié)點(diǎn)編號并沒有問題(編號正確,而且是逆時針的),而且,對于含通用接觸的inp(不

   能計(jì)算)和不含通用接觸的inp(能計(jì)算)兩者的檢驗(yàn)結(jié)果是一樣的。

3  既然檢查了abaqus拋出的節(jié)點(diǎn)編號問題后發(fā)現(xiàn),節(jié)點(diǎn)編號是沒有問題的,那么接下來怎么辦,...,陷入

   死循環(huán)了,研究了好多天不能找到問題所在...

4  等等,除了這個拋出的錯誤外,還有警告呢,警告有兩個,第二個沒什么用,第一個很重要,它說有15416個

   節(jié)點(diǎn)被調(diào)整了,后續(xù)的調(diào)整可能引起這些及誒單不在他們的主面上,這些節(jié)點(diǎn)已經(jīng)被定義在了節(jié)點(diǎn)集合WarnNodeAdjust中,

   那么我們看一下這個調(diào)整的節(jié)點(diǎn)集合吧,進(jìn)入后處理,發(fā)現(xiàn)這個集合就是我們插入cohesive區(qū)域內(nèi)的一些節(jié)點(diǎn),

   也就是說在使用批量插入cohesive單元+隱式分析+通用接觸時,節(jié)點(diǎn)被abaqus自動進(jìn)行了調(diào)整,然后導(dǎo)致了不能計(jì)算。

5  因?yàn)閕np文件根本不能正確提交計(jì)算,所以我們只能看.pre前處理文件了,對比有無通用接觸的兩個前處理文件發(fā)現(xiàn),

   包含通用接觸的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      

   好吧,這應(yīng)該就是癥結(jié)所在了,我們有分析了隱式分析有無通用接觸情況下的pre文件,發(fā)現(xiàn)只有批量插入cohesive單元+隱式分析+通用接觸

   的時候會出現(xiàn)節(jié)點(diǎn)被調(diào)整的問題,總的來說就是因?yàn)閍baqus自動調(diào)整了節(jié)點(diǎn),導(dǎo)致了我們的inp文件不能計(jì)算,而出現(xiàn)的節(jié)點(diǎn)編號錯誤的提示。

   因?yàn)橹笥?jì)算的時候才會自動調(diào)整這些節(jié)點(diǎn),所以有無通用接觸的兩種情況下輸出的inp都是相同的,因此只分析inp文件我們并不能找到根本

   原因。

解決辦法:

  1 刪除報錯的單元(治標(biāo)不治本,而且如果報錯有很多單元怎么辦呢?)

  2 不使用通用接觸(但是如果是模型中開裂后的裂縫又在壓應(yīng)力作用下閉合了怎么辦?)

  3 把模型的所有外表面建立為一個surface集合,通用接觸處使用selected surface pair,可以成功提交計(jì)算,但是計(jì)算特別慢,而且計(jì)算結(jié)果總

    感覺不對,大家可以嘗試。

  4 徹底的解決辦法暫時沒找到(或者根本沒有,因?yàn)榕坎迦隿ohesive單元,隱式分析只要使用批量插入cohesive單元,它就會調(diào)整節(jié)點(diǎn)),如果

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

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

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

登錄后免費(fèi)查看全文
立即登錄
App下載
技術(shù)鄰APP
工程師必備
  • 項(xiàng)目客服
  • 培訓(xùn)客服
  • 平臺客服

TOP

16
5
6