接触分析定稿.ppt

8-* ., 刚性面 下面是对于一个刚性面, 低次单元和高次单元两种划分方式都用了的例子. 低次单元的计算量少, 但是需要更多的单元. 高次单元的计算量大, 但是, 需要的单元少. 低次 高次 8-* ., 改变刚性面的形状 对于具有低次或高次目标元素(不是基本元素)的已划分网格的目标面, 目标节点能够旋转到局部坐标系. 这允许目标面在一个分析中改变形状. 例如, 希望用刚性圆柱面进行扩径. 此圆柱面可以划分为高次四边形, 将目标节点从直角坐标系旋转到圆柱坐标系, 并在刚性目标面上施加径向位移, 就可以对圆柱扩径. 刚性面 8-* ., 接触算法 选择一个接触算法 关键字选项(2) 增强的拉格朗日法(关键字选项(2)=0)是缺省选项, 推荐于一般应用. 它对罚刚度不太敏感, 但是也要求给出一个侵入容差. 能够用罚函数法(关键字选项 (2)=1)这个选项. 它推荐应用于单元非常扭曲、大摩擦系数和/或用增大的拉格朗日法收敛行为不好的问题. 8-* ., 确定罚刚度 对于面对面接触单元, ANSYS基于单元类型、材料性质和它下面的单元尺寸确定接触刚度. 可以用实常数FKN给出接触刚度的一个比例系数或绝对值. 惩罚刚度(FKN)应该足够大, 使接触侵入量小. 同时也应足够小, 使问题没有病态矩阵. FKN值通常在 0.01~10之间. 对于体积变形问题, 用值1.0 (默认), 对于弯曲问题, 用值 0.1. 8-* ., 确定侵入容差 侵入容差(FTOLN)是与接触单元下面的实体单元深度(h)相乘的比例系数. 若此值太小, 会引起收敛困难. 绝对不要用太小的容差! 增大罚刚度(FKN)将减少侵入. 将FKN值增大100倍会相应地减少侵入, 但是接 触压力只改变5%. 8-* ., 确定罚刚度 确定一个好的刚度值需要做些实验. 下面的处理可以作为指南: 1. 开始用一个较小的 FKN值. 2. 运行最终载荷的一部分的分析. 3. 检查侵入和每一子步的平衡迭代数. 如果收敛受侵入容差的驱使, 可能是FKN 值估计不足或 FTOLN 值太小. 如果需要多次迭代才能使残余力收敛, 而不是侵入, FKN值可能估计得太高. 4. 调整 FKN 或 FTOLN值, 运行整个分析. 注意:从一个载荷步到另一个载荷步, FKN 和 FTOLN 值都能修改, 并且能在重起动分析中调整. 8-* ., 摩擦 在接触模型中包括摩擦, 就允许在两个接触体间产生剪力. 在库仑摩擦模型中, 在两个表面开始相对滑动之前, 这两个表面承受的剪应力要上升到一个定值. 此阶段称为粘着. 当滑动发生时, 库仑摩擦模型定义平衡的摩擦应力? : 式中, p是单位接触压力, ? 是摩擦系数(用材料性质MU定义). 一旦剪应力超过此定值, 这两个表面将相对滑动. 此阶段就是滑动. 摩擦系数 ?可以是任意非负值. 8-* ., 摩擦 |?| ?max ? p 滑动 粘着 ANSYS 提供一个实常数TAUMAX, 用于定义等效剪应力的最大值. 这样, 不管接触压力值是多大, 只要等效剪应力达到最大值TAUMAX, 就会发生滑动. 该剪应力极限值通常用于接触压力会变得非常大的情况. TAUMAX的一个合理的上限估计值是: 式中 ?y 是表面附近材料的范?米赛斯屈服应力. TAUMAX的最好值通常来源于经验数据. 8-* ., 摩擦 包含摩擦的接触问题会产生不对称的刚度矩阵. 然而, 使用不对称的方程求解器比对称求解器的计算费用更高. 由于这个原因, ANSYS用一种对称算法进行计算. 利用此算法, 可以解决多数包含摩擦的接触问题. 如果遇到了一个收敛缓慢的问题, 可以用不对称求解选项 (关键选项(6)=1). 在这种情况下, 应该用稀疏矩阵求解器 (优选的)或波前法求解器. 8-* ., 接触检查点 ANSYS面对面单元默认用高斯积分点作为接触检查点. 高斯点 接触面 接触面 ANSYS 允许通过设置关键字选项(4)=1将接触检查点移到节点. 为了保证收敛和精度, 建议保持默认设置. 8-* ., 高斯点的优点 对于高次单元, 节点力连续. 高斯点不能“ 滑离”目标面的边缘. 垂直于高斯点只定义了唯一的侵入, 不要求目标面的平滑程度. 接触检查点更多. Gauss point 8-* ., 接触检查点 如果有一个包括角部接触的问题, 使用高斯点作为接触检查点会导致角部的过分侵入. 这时也不用将接触检查点移到节点, 可以将面对面接触单元与节点对面接触单元混合起来用. 面对面单元 节点对面单元用于在角部建立接触模型 8-* ., 对称接触与不对称接触 在