系列2-原创#cohesive+隐式分析+通用接触--报错node numbering原因
本帖最后由 weizhiyuanmali 于 2019-1-20 18:33 编辑********************************************************************************************
首先声明（做了模型分析）：
批量插入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-18ADJUSTED 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-18ADJUSTED 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+、专门定制视频、全程亲自教学、各种模型调试及解答问题等等，倾囊相教）
最近也遇到这个问题，隐式分析+cohesive+通用接触，提交提示零/负厚度单元。感觉只要用了通用接触，abaqus就会检查整个模型单元，并且把零厚度cohesive当作实体单元检查，导致报错。不去的是不是有这个问题。谢谢你的经验分享 隐式分析只能用face to face？
页:
[1]