ansys有限元分析重点关于接触刚度.pdfVIP

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