接触问题的有限元分析.ppt

6.3 接触对生成与参数设置 （2）接触分析的参数设置 5）点─点接触分析 3.设置平衡方程个数为适合时间步长大小的数值，在“SOLCONTROL”打开时，依据不同的问题，其缺省值为15至26 迭代。 4.打开预测—校正选项。 5.设置N—S选项为FULL。 6.对于涉及到的摩擦分析，其中法向和切向运动是强耦合时，“NROPT，UNSYM”是很有用的。 7.许多收敛困难是由于使用了大的接触刚度值。 8.相反，如果出现了过渗透，则可能是KN的值太小。 6.4 赫兹接触理论 （1）两球接触或球与平面接触 6.4 赫兹接触理论 （1）两球接触或球与平面接触 （2）两圆柱接触或圆柱与平面接触 6.4 赫兹接触理论 （2）两圆柱接触或圆柱与平面接触 6.5 典型实例分析 （1）平面接触分析 问题描述 理论分析 两平板为均匀的直线接触，因而两平板的应力、变形相等，Y向应力为100 MPa。 ANSYS 分析的操作步骤 (1) 选取单元和输入材料属性 ① 选取单元：Main MenuPreprocessorElement typeAdd/Edit/Delete，出现一个对话框，单击[Add]，又出现一个“Element Type”对话框，在“Library of Element Type”左面的列表栏中选择“Solid”后面列表栏中选择“Quad 4 node 182”，单击[OK]；再单击[Close]。 6.5 典型实例分析 （1）平面接触分析 ② 输入材料参数：Main MenuPreprocessorMaterial Props Material Models，在“Material Model Available”下面的对话框中，双击打开“StructuralLinearElastic Isotropic”，又出现一个对话框，输入“EX=1000, PRXY=0.3”，单击[OK]，单击“Material Exit”，完成材料属性的设置。 (2) 建立网格模型 ①建立几何模型：MainMenuPreprocessorModelingCreateAreasRectangleBy Dimensions，弹出一个“Create Rectangle by Dimensions”的对话框，在“X1，X2”后面依次输入“0”和“5”，在“Y1，Y2”后面依次输入“0”和“2.5”，单击[Apply] 6.5 典型实例分析 （1）平面接触分析 在弹出的同一个对话框中，仅在“Y1，Y2”后面依次输入“2.5”和“5”，单击[OK]，输出窗口上显示2个矩形。 ②设置单元大小：MainMenuPreprocessorMeshingSize CntrlsManualSizeGlobal Size，在“SIZE Element edge length”后面输入“0.25”，单元[OK]。 ③对面划分单元：Main MenuPreprocessorMeshingMeshAreasMapped3 or 4 sided，弹出一个“Mesh Areas”的拾取框，单击[Pick All]，则完成面网格的划分，生成的网格结果 6.3 接触对生成与参数设置 （2）接触分析的参数设置 1）面─面接触分析——决定接触刚度 ②用最终载荷的一部分来运行分析(刚好使接触完全建立)。 ③检查每个子步中的渗透和平衡迭代的次数。 如果整体不收敛是由于渗透量过大引起，而不是由残余力和位移增量引起，这时FKN的值太大或者FTOL的值过小； 如果整体收敛需要许多平衡迭代步才能满足残余力和位移的收敛，则表明FKN和FTOL的值高估了。 ④根据需要调整FKN、FKT、FTOLN或SLOT的值，开始一个完整的分析过程，如果在整体平衡迭代中，用户能够支配渗透量，这时可以增大FTOLN的值来获得更大的许可渗透量，或者增大FKN的值。 1）面─面接触分析——选择摩擦模型 粘和状态(stick）：在基本的库仑摩擦模型中，两个接触面在开始相互滑动之前，在界面上会有某一大小的剪应力产生。 6.3 接触对生成与参数设置 （2）接触分析的参数设置 1）面─面接触分析——选择摩擦模型 粘聚力 滑动状态(sliding): 一旦剪应力超过此值后，两个表面将开始相互滑动。 粘和/滑动计算决定什么时候一个点从粘和状态到滑动状态，或从滑动状态到粘和状态。 对于另一种系统支持的摩擦模型，可以使用USERFRIC子程序来定义用户需要的摩擦模型。 对于无摩擦、粗糙和绑定接触，接触单元的刚度矩阵是对称的。 大多数摩擦接触问题采用对称求解器来求解，但如果