接触报错，接触从有到无接触

xiaaiyin

*** WARNING ***                         CP =    3160.128   TIME= 18:48:54
Contact element 13203 (real ID 6) status changes abruptly from contact  
(with target element 13079) -> no-contact.                              
*** WARNING ***                         CP =    3366.767   TIME= 18:50:56
Contact element 12905 (real ID 5) status changes abruptly from contact  
(with target element 12647) -> no-contact.                              
我后来查看了这两个接触对，图片如下

这张是Contact element 13203和target element 13079，图中选中的那一点即为这两个单元。（是螺栓帽的小环和钢板的接触）

这张是Contact element 12905和target element 12647，图中所选的那一小块即为这两个单元。（这是上下两块钢板所形成的接触）
中间的孔均为螺栓孔。
然后就出现*** ERROR ***                           CP =    3571.658   TIME= 18:52:50
One or more elements have become highly distorted.  Excessive           
distortion of elements is usually a symptom indicating the need for     
corrective action elsewhere.  Try incrementing the load more slowly     
(increase the number of substeps or decrease the time step size).  You  
may need to improve your mesh to obtain elements with better aspect     
ratios.  Also consider the behavior of materials, contact pairs,        
and/or constraint equations.  If this message appears in the first      
iteration of first substep, be sure to perform element shape checking.  
*** WARNING ***                         CP =    3571.846   TIME= 19:06:11
The unconverged solution (identified as time 5 substep 999999) is      
output for analysis debug purposes.  Results should not be used for     
any other purpose.                                                      
*** ERROR ***                           CP =    3573.265   TIME= 19:06:12
LSSOLVE looping is stopped because solution is unconvergent and NCNV,1  
set.      
不收敛了！                                         
下面使用cncheck检查的结果，希望各位高手施于援手~~
并且用12.0总是算不了几步，用10.0却可以计算很多荷载步，不知何故？
*** NOTE ***                            CP =       5.366   TIME= 19:58:35
The default contact stiffness starting from Release 12.0 is no longer   
affected by defined plasticities input by the TB commands.  You should  
confirm that the appropriate contact stiffness was used.               
*** NOTE ***                            CP =       5.569   TIME= 19:58:35
Deformable-deformable contact pair identified by real constant set 1   
and contact element type 4 has been set up.                             
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN          1.0000   
The resulting contact stiffness              0.41200E+06
Default penetration tolerance factor FTOLN   0.10000   
The resulting penetration tolerance           1.0000   
Max. initial friction coefficient MU         0.40000   
Default tangent stiffness factor FKT          1.0000   
Default Max. friction stress TAUMAX          0.10000E+21
Average contact surface length                14.669   
Average contact pair depth                    10.000   
Default pinball region factor PINB            2.0000   
The resulting pinball region                  20.000   
*WARNING*: Initial penetration is included.
*** NOTE ***                            CP =       5.569   TIME= 19:58:35
Max.  Initial penetration 0 was detected between contact element 11583  
and target element 11512.                                               
****************************************
  
*** NOTE ***                            CP =       5.569   TIME= 19:58:35
Deformable-deformable contact pair identified by real constant set 2   
and contact element type 4 has been set up.                             
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN          1.0000   
The resulting contact stiffness              0.41200E+06
Default penetration tolerance factor FTOLN   0.10000   
The resulting penetration tolerance           1.0000   
Max. initial friction coefficient MU         0.40000   
Default tangent stiffness factor FKT          1.0000   
Default Max. friction stress TAUMAX          0.10000E+21
Average contact surface length                15.444   
Average contact pair depth                    10.000   
Default pinball region factor PINB            2.0000   
The resulting pinball region                  20.000   
*WARNING*: Initial penetration is included.
*** NOTE ***                            CP =       5.569   TIME= 19:58:35
Max.  Initial penetration 2.842170943E-14 was detected between contact  
element 12017 and target element 11759.                                 
****************************************
  
*** NOTE ***                            CP =       5.569   TIME= 19:58:35
Deformable-deformable contact pair identified by real constant set 3   
and contact element type 4 has been set up.                             
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN          1.0000   
The resulting contact stiffness              0.51500E+06
Default penetration tolerance factor FTOLN   0.10000   
The resulting penetration tolerance          0.80000   
Max. initial friction coefficient MU         0.40000   
Default tangent stiffness factor FKT          1.0000   
Default Max. friction stress TAUMAX          0.10000E+21
Average contact surface length                14.669   
Average contact pair depth                    8.0000   
Default pinball region factor PINB            2.0000   
The resulting pinball region                  16.000   
*WARNING*: Initial penetration is included.
*** NOTE ***                            CP =       5.569   TIME= 19:58:35
Max.  Initial penetration 1.421085472E-14 was detected between contact  
element 12292 and target element 12192.                                 
****************************************
  
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Deformable-deformable contact pair identified by real constant set 4   
and contact element type 4 has been set up.                             
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN          1.0000   
The resulting contact stiffness              0.51500E+06
Default penetration tolerance factor FTOLN   0.10000   
The resulting penetration tolerance          0.80000   
Max. initial friction coefficient MU         0.40000   
Default tangent stiffness factor FKT          1.0000   
Default Max. friction stress TAUMAX          0.10000E+21
Average contact surface length                14.669   
Average contact pair depth                    8.0000   
Default pinball region factor PINB            2.0000   
The resulting pinball region                  16.000   
*WARNING*: Initial penetration is included.
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Max.  Initial penetration 1.421085472E-14 was detected between contact  
element 12485 and target element 12361.                                 
****************************************
  
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Deformable-deformable contact pair identified by real constant set 5   
and contact element type 4 has been set up.                             
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN          1.0000   
The resulting contact stiffness              0.41200E+06
Default penetration tolerance factor FTOLN   0.10000   
The resulting penetration tolerance           1.0000   
Max. initial friction coefficient MU         0.40000   
Default tangent stiffness factor FKT          1.0000   
Default Max. friction stress TAUMAX          0.10000E+21
Average contact surface length                15.444   
Average contact pair depth                    10.000   
Default pinball region factor PINB            2.0000   
The resulting pinball region                  20.000   
*WARNING*: Initial penetration is included.
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Max.  Initial penetration 2.842170943E-14 was detected between contact  
element 12830 and target element 12574.                                 
****************************************
  
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Deformable-deformable contact pair identified by real constant set 6   
and contact element type 4 has been set up.                             
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN          1.0000   
The resulting contact stiffness              0.41200E+06
Default penetration tolerance factor FTOLN   0.10000   
The resulting penetration tolerance           1.0000   
Max. initial friction coefficient MU         0.40000   
Default tangent stiffness factor FKT          1.0000   
Default Max. friction stress TAUMAX          0.10000E+21
Average contact surface length                14.669   
Average contact pair depth                    10.000   
Default pinball region factor PINB            2.0000   
The resulting pinball region                  20.000   
*WARNING*: Initial penetration is included.
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Max.  Initial penetration 0 was detected between contact element 13191  
and target element 13122.                                               
****************************************
  
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Deformable-deformable contact pair identified by real constant set 7   
and contact element type 4 has been set up.                             
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN          1.0000   
The resulting contact stiffness              0.29429E+06
Default penetration tolerance factor FTOLN   0.10000   
The resulting penetration tolerance           1.4000   
Max. initial friction coefficient MU         0.40000   
Default tangent stiffness factor FKT          1.0000   
Default Max. friction stress TAUMAX          0.10000E+21
Average contact surface length                14.669   
Average contact pair depth                    14.000   
Default pinball region factor PINB            2.0000   
The resulting pinball region                  28.000   
*WARNING*: Initial penetration is included.
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Max.  Initial penetration 3.552713679E-15 was detected between contact  
element 13329 and target element 13276.                                 
****************************************
  
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Deformable-deformable contact pair identified by real constant set 8   
and contact element type 4 has been set up.                             
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN          1.0000   
The resulting contact stiffness              0.29429E+06
Default penetration tolerance factor FTOLN   0.10000   
The resulting penetration tolerance           1.4000   
Max. initial friction coefficient MU         0.40000   
Default tangent stiffness factor FKT          1.0000   
Default Max. friction stress TAUMAX          0.10000E+21
Average contact surface length                14.448   
Average contact pair depth                    14.000   
Default pinball region factor PINB            2.0000   
The resulting pinball region                  28.000   
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
One of the contact searching regions contains at least 68 target        
elements.  You may reduce the pinball radius.                           
*WARNING*: Initial penetration is included.
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Max.  Initial penetration 7.105427358E-15 was detected between contact  
element 13623 and target element 13465.                                 
****************************************
  
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Deformable-deformable contact pair identified by real constant set 9   
and contact element type 4 has been set up.                             
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN          1.0000   
The resulting contact stiffness              0.29429E+06
Default penetration tolerance factor FTOLN   0.10000   
The resulting penetration tolerance           1.4000   
Max. initial friction coefficient MU         0.40000   
Default tangent stiffness factor FKT          1.0000   
Default Max. friction stress TAUMAX          0.10000E+21
Average contact surface length                14.669   
Average contact pair depth                    14.000   
Default pinball region factor PINB            2.0000   
The resulting pinball region                  28.000   
*WARNING*: Initial penetration is included.
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Max.  Initial penetration 0 was detected between contact element 13938  
and target element 13875.                                               
****************************************
  
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Deformable-deformable contact pair identified by real constant set 10   
and contact element type 4 has been set up.                             
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN          1.0000   
The resulting contact stiffness              0.29429E+06
Default penetration tolerance factor FTOLN   0.10000   
The resulting penetration tolerance           1.4000   
Max. initial friction coefficient MU         0.40000   
Default tangent stiffness factor FKT          1.0000   
Default Max. friction stress TAUMAX          0.10000E+21
Average contact surface length                14.537   
Average contact pair depth                    14.000   
Default pinball region factor PINB            2.0000   
The resulting pinball region                  28.000   
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
One of the contact searching regions contains at least 68 target        
elements.  You may reduce the pinball radius.                           
*WARNING*: Initial penetration is included.
*** NOTE ***                            CP =       5.585   TIME= 19:58:35
Max.  Initial penetration 2.842170943E-14 was detected between contact  
element 14242 and target element 14084.                                 
****************************************
  
*** WARNING ***                         CP =       5.585   TIME= 19:58:35
Maximum friction coefficient 0.4 is defined in the model.  Switch to   
the unsymmetric solver (NROP,UNSYM) instead if convergence difficulty   
is encountered.                                                         
      10 CONTACT PAIRS ARE SELECTED
     CONTACT PAIR HAVING REAL ID =     1 IS INITIALLY CLOSED
     CONTACT PAIR HAVING REAL ID =     2 IS INITIALLY CLOSED
     CONTACT PAIR HAVING REAL ID =     3 IS INITIALLY CLOSED
     CONTACT PAIR HAVING REAL ID =     4 IS INITIALLY CLOSED
     CONTACT PAIR HAVING REAL ID =     5 IS INITIALLY CLOSED
     CONTACT PAIR HAVING REAL ID =     6 IS INITIALLY CLOSED
     CONTACT PAIR HAVING REAL ID =     7 IS INITIALLY CLOSED
     CONTACT PAIR HAVING REAL ID =     8 IS INITIALLY CLOSED
     CONTACT PAIR HAVING REAL ID =     9 IS INITIALLY CLOSED
     CONTACT PAIR HAVING REAL ID =    10 IS INITIALLY CLOSED

请问这些接触都是哪错了呢。
下面是面面接触的生成命令流


/prep7
allsel
!生成接触单元
*do,k,3,21,2                      !10组循环
ET,3,targe170                  !定义三维目标单元为targe170
ET,4,conta174             !定义三维接触单元为conta17
j=(k-1)/2
r,j                                   !定义接触实常数j
real,j              
!前面已经定义了mp,mu,8,0.3，材料8，在这就不定义了
mat,8                !指定材料8
cmsel,s,tar%j%       !target
type,3                               !指定目标单元类型k
nsla,s,1                              !选择目标面上的所有节点，包括边界点
esln,s,0                              !选择依附于这些节点的单元
esurf,all                             !构建目标单元 ESURF, XNODE, Tlab, Shape
cmsel,s,con%j%  !contact
type,4                               !指定接触单元类型k+1
nsla,s,1                               !选择接触面上的所有节点，包括边界点
esln,s,0                                  !选择依附于这些节点的单元
esurf,all
*enddo
allsel,all

孙武

顶起！

孙武

Contact element 13203 (real ID 6) status changes abruptly from contact  
(with target element 13079) -> no-contact.
求高手告诉如何解决这个错误

hg_boy

这不是错误。
只是一个警告信息。

这个信息仅仅是告诉你在结构受到外载荷过程中，该处的接触状态发生了变化。
本来是接触的两个单元，因为外力作用导致他们脱开，由接触状态变为非接触状态。

xiaaiyin

hg_boy 发表于 2012-4-26 20:40
这不是错误。
只是一个警告信息。

这个接触警告后，后面很快就出现单元扭曲错误了，请教版主，这两者之间有什么联系吗

hg_boy

这只是一个辅助诊断信息，你可以根据这个提示信息，查看是哪里脱开了，为什么这里的接触会脱开。

有些情况，脱开是正常范围内的。比如载荷的方向本身就是使两个部件朝着分离的方向的，如果载荷增大了，就有可能使局部接触状态发生改变而脱开。

有时候，脱开意味着结构状态发生变化，譬如如果载荷过大，接触界面上的摩擦力不足以阻止相对滑动趋势，两个部件要发生相对移动了，此时结构处于欠约束状态了，计算就没法收敛。

这个和你后面单元过度扭曲，没有必然的联系，只是有可能有一定的关联。

单元过度扭曲，可能是你施加的载荷过大。
也有可能你的载荷并不大，但是材料参数的单位搞错了，使得材料很弱，显得你施加的载荷很大。
还有可能是你是接触相关设置不当，导致局部过大的载荷，使得单元变形过度，从而扭曲。
还有可能和网格有关。

还有可能就是因为接触状态发生改变，使得载荷的分配发生变化，从而局部单元受力过大而扭曲。

可能性很多，具体要结合模型来判断，没有看到模型，只能泛泛而谈。

xiaaiyin

hg_boy 发表于 2012-4-26 22:33  这只是一个辅助诊断信息，你可以根据这个提示信息，查看是哪里脱开了，为什么这里的接触会脱开。   有些情况 ...

感谢版主的专业回答，为我指引了光亮之路

xxh060113

hg_boy 发表于 2012-4-26 22:33
这只是一个辅助诊断信息，你可以根据这个提示信息，查看是哪里脱开了，为什么这里的接触会脱开。

有些情况 ...

对螺栓施加预紧力时后进行非线性求解常会出现一个节点位移超限，超过默认值，这又是怎么回事？常常会导致求解不收敛，难道是给的荷载步太小的缘故？还是没有解决问题。。。

xuanzhuan

不加任何载荷子步试试。