【讨论】接触单元技术----摘自ANSYS用户区

edodo

以下摘自ANSYS公司的官方网站，就其中我认为唯一的方法就是我们必须试用不同的值直到找到正确的值。嗬嗬！不就等于说计算出来的结果肯定不是真正的解么。只是一个我们通过不断调整而得到的一个近似解么！那不也就是说接肯定不正确！：）
看来以后接触问题还是不要用ANSYS来做了！：）
  
接触单元技术
  
    在有限元分析中，接触单元通常用来描述两物体相互接触或滑动的界面。近年来，ANSYS开发了一系列的接触单元。刚开始有节点对节点单元CONTAC12和CONTAC52，接着有节点对地单元CONTAC26，然后有节点对面单元CONTAC48和CONTAC49。最近几年，我们引入一类面对面接触单元CONTA169和CONTA174，同时还有一种新的节点对节点单元CONTA178。
  
   虽然接触单元的参数具有多样性，但我们在使用他们时可谨记重要的一点，他们具有一个共同的特点，即除了CONTA178的KEYOPT（2）=0或1外，所有的接触单元都有接触刚度。在现实中实际上相邻结构之间只是一种空隙，但在有限元分析中，这种空隙是一带有刚度的接触单元，这是因为通过刚度矩阵来实现接触算法的。一些接触单元要求使用者输入刚度值，同时另外的接触单元若没有输入则使用缺省值。分析工程师所面对的问题就是针对给定的条件确定一个合理的刚度值。如果过高，问题将会不收敛，如果过低，可能得到错误的结果。那么我们所面对的问题是怎样才能找到一个正确的刚度值？
    
    我认为唯一的方法就是我们必须试用不同的值直到找到正确的值。也就是刚开始我们应该使用一个较小的值，然后稳步的增加直到分析的结果不再有什么变化。那么对于我们这一特定分析的问题，这一点就是我们所想要的合适值。
   
    我们可举例说明，如图1所示，平行放置两个悬臂梁，并有少许的交迭，下面的左边固支，上面的右边固支，当在上面梁的自由端施加一个向下位移时，梁变形弯曲并接触下面的梁，然后一起向下运动。用SOLID45单元划分梁，用TARGE170和CONTA174面面接触单元来描述相互作用。在此基础上，把CONTA174单元的刚度从非常低变到非常高，从而来观察它对结果的影响和收敛的迭代次数。图2说明了下梁自由端的偏移随接触单元刚度的变化情况，当刚度增加时，偏移量接近一个常数值（我们可以假定它是一个"正确"的结果。）图3说明求解所需的迭代次数，当接触单元刚度增加时，求解所需的迭代次数也是增加的，并服从指数关系。如果刚度过高，问题很有可能根本就不收敛。图4说明在上梁自由端接触单元的渗透量，当刚度增加时，渗透量降低。
    从这些图可知，当接触单元的刚度为10e6时，可获得合理精确的结果。任何大于该值的刚度对下梁的偏移量没有什么影响，而求解所需的迭代数却显著的增加。对于这个题目，10e6的刚度是很适合的。但是，如果改变边界条件、网格密度、两梁之间的相对位置、材料特性或梁的几何形状，能获得满意结果的接触刚度值将是不同的。比如，如果网格密度增加，则接触单元数将增加，每一个单元上的载荷将降低。如果接触单元数增加两倍，一个合适的接触单元刚度值应为原来的一半。
由于每个题目都是不一样的，所以在求解之前并没有通用的方法来确定接触单元刚度的最佳值。我们不得不试算一个我们认为合适的值然后查看计算结果。一个有经验的分析工程师可能只查看一个计算结果来判定所取值的合适度，但对于大多数情况而言，最好用一个合理而不过度精确的刚度值进行第一次求解，然后用10倍于该值的刚度进行第二次求解，如果两者结果相差很小，而迭代数增加很多，那么我们则正好取得了曲线上的突变点，从而获得相当好的结果。
接触单元刚度问题仅仅是一个例子，即对于分析工程师来说，总是置疑于分析结果的正确与否是非常重要的，并要意识到数值仿真的局限性和潜在的假设及他们怎样影响所分析问题的结果。  
  
  
  

  

edodo

现在没有了！
对不起！

杜晓帆

没有图啊，不过也谢谢，新手