《lecture6-cAbaqusStandard中的接触问题.ppt》由会员分享,可在线阅读,更多相关《lecture6-cAbaqusStandard中的接触问题.ppt(19页珍藏版)》请在金锄头文库上搜索。
1、 Dassault Systmes, 2008Abaqus/Standard中的接触问题Lecture 6 Dassault Systmes, 2008Introduction to Abaqus/Standard and Abaqus/ExplicitL6.2简介将接触作将接触作为模型定模型定义的一部分的一部分主主/从方法中网格密度的考从方法中网格密度的考虑Abaqus/Standard中的接触中的接触逻辑习题5:密封接触:密封接触 Dassault Systmes, 2008将接触作为模型定义的一部分 Dassault Systmes, 2008Introduction to Abaqus
2、/Standard and Abaqus/ExplicitL6.4将接触作为模型定义的一部分在Abaqus/Standard 整个contact定义作为模型数据必须在step之前接触对在分析历程中可以使用*MODEL CHANGE选项激活或者不激活*MODEL CHANGE, TYPE=CONTACT PAIR, ADD | REMOVEsurface_1, surface_2 Dassault Systmes, 2008在严格的Master/Slave定义中要考虑网格密度 Dassault Systmes, 2008Introduction to Abaqus/Standard and Ab
3、aqus/ExplicitL6.6在严格的Master/Slave定义中要考虑网格密度从属表面应该是网格划分得更精细的表面.如果网格密度相近似,从属表面应由更为柔软的材料组成.IncorrectMaster surface placed on fine mesh Gross penetration into slave surface.CorrectMaster surface placed on coarse mesh Minimal penetration into slave surface. Dassault Systmes, 2008在 Abaqus/Standard中CONTACT
4、规则 Dassault Systmes, 2008Introduction to Abaqus/Standard and Abaqus/ExplicitL6.8在 Abaqus/Standard中CONTACT规则对于点接触需要施加于点接触需要施加约束束有多种不同的约束施加方法.缺省情况下, Abaqus/Standard使用拉格朗日多项式方法.对于每个潜在的接触点,接触条件通过单个的、非线性的、不等式约束描述:其中h 为“侵彻”, uN 为自由度 Dassault Systmes, 2008Introduction to Abaqus/Standard and Abaqus/Explicit
5、L6.9在 Abaqus/Standard中CONTACT规则一个增量步内状一个增量步内状态的原理的原理图(默(默认状状态) 检查求解是否收敛4形成并求解系统方程2确定初始激活的接触1开始增量步结束增量步Yes(至少一个收敛准则不满足)确定是否趋于收敛5(减少增量步大小并重新求解)确定接触状态的变化3No(在收敛容差范围内)NoYesNewton 迭代 Dassault Systmes, 2008Introduction to Abaqus/Standard and Abaqus/ExplicitL6.10在 Abaqus/Standard中CONTACT规则默默认的接触算法(参考前的接触算法
6、(参考前页的流程的流程图) 确定每个点当前的接触状态(接触或张开)对于一个分析步的第一个增量步,基于初始模型的状态;否则,基于求解的外插 计算刚度,施加相应的接触约束的影响。形成系统方程并求解。 接触压力和间隙是否与接触状态一致。接触状态改变(张开/闭合或粘结/滑移)通常导致系统方程的重大变化。由接触状态改变产生的迭代被标记为严重不连续迭代(SDIs) 123 Dassault Systmes, 2008Introduction to Abaqus/Standard and Abaqus/ExplicitL6.11在 Abaqus/Standard中CONTACT规则是否达到收敛?默认的, A
7、baqus量化与SDIs关联的接触不相容性。那些关于接触状态的不相容性必须十分小才能达到收敛。也必须确保力的剩余误差和求解修正很小才能达到平衡如果接触状态和平衡条件满足各自的收敛准则,这个增量步完成如果没有达到收敛,它可能会收敛吗?Abaqus估算趋势,例如连续迭代中接触状态改变的次数,以确定是否继续迭代或减小增量步长重算如果可能收敛,更新基于 的接触条件和刚度,重新求解系统方程;否则,利用更小的增量步重试。453 Dassault Systmes, 2008Introduction to Abaqus/Standard and Abaqus/ExplicitL6.12在 Abaqus/Sta
8、ndard中CONTACT规则接触打印接触打印输出例子出例子参考: Example Problem 1.3.4, “Deep drawing of a cylindrical cup”状态 (.sta) 文件:SUMMARY OF JOB INFORMATION: MONITOR NODE: 200 DOF: 2 STEP INC ATT SEVERE EQUIL TOTAL TOTAL STEP INC OF DOF IF DISCON ITERS ITERS TIME/ TIME/LPF TIME/LPF MONITOR RIKS ITERS FREQ 1 1 1 1 1 2 1.00 1
9、.00 1.000 0.000 2 1 1 0 1 1 2.00 1.00 1.000 0.000 3 1 1 10 0 10 2.01 0.0100 0.01000 -0.000600 3 2 1 7 1 8 2.02 0.0200 0.01000 -0.00120 3 3 1U 10 0 10 2.02 0.0200 0.01500 -0.00120 3 3 2 5 0 5 2.02 0.0238 0.003750 -0.00142 3 4 1 3 1 4 2.03 0.0294 0.005625 -0.00176 3 5 1 2 3 5 2.04 0.0378 0.008438 -0.0
10、0227 3 6 1 4 1 5 2.05 0.0505 0.01266 -0.00303 3 7 1 4 1 5 2.07 0.0695 0.01898 -0.00417 3 8 1 6 1 7 2.10 0.0979 0.02848 -0.00588 3 9 1U 8 0 8 2.10 0.0979 0.04271 -0.00588 3 9 2 4 3 7 2.11 0.109 0.01068 -0.00652 3 10 1 2 3 5 2.12 0.125 0.01602 -0.00748 3 11 1 4 2 6 2.15 0.149 0.02403 -0.00892 . . . Da
11、ssault Systmes, 2008Introduction to Abaqus/Standard and Abaqus/ExplicitL6.13在 Abaqus/Standard中CONTACT规则Message文件文件, Step 3, Increment 6: INCREMENT 6 STARTS. ATTEMPT NUMBER 1, TIME INCREMENT 1.266E-02 CONTACT PAIR (ASURF,BSURF) NODE 167 IS NOW SLIPPING. CONTACT PAIR (ASURF,BSURF) NODE 171 IS NOW SLIP
12、PING. : : : CONTACT PAIR (ASURF,BSURF) NODE 153 OPENS. CONTACT PRESSURE/FORCE IS -845822. CONTACT PAIR (ASURF,BSURF) NODE 161 OPENS. CONTACT PRESSURE/FORCE IS -1.50656E+006. CONTACT PAIR (ASURF,BSURF) NODE 163 OPENS. CONTACT PRESSURE/FORCE IS -108355. CONTACT PAIR (ASURF,BSURF) NODE 165 OPENS. CONTA
13、CT PRESSURE/FORCE IS -620880. CONTACT PAIR (CSURF,DSURF) NODE 363 OPENS. CONTACT PRESSURE/FORCE IS -3.5893E+006. CONTACT PAIR (ESURF,FSURF) NODE 309 IS NOW SLIPPING. 6 SEVERE DISCONTINUITIES OCCURRED DURING THIS ITERATION. 5 POINTS CHANGED FROM CLOSED TO OPEN 1 POINTS CHANGED FROM STICKING TO SLIPPI
14、NG从节点滑移;如果应用Lagrange摩擦或发生滑移反转,将产生SDIs。假定的接触状态中不相容的被检测到 SDIDue to slip reversal*PRINT, CONTACT=YES,将细节打印输出。(对于错误信息的处理非常有用) Dassault Systmes, 2008Introduction to Abaqus/Standard and Abaqus/ExplicitL6.14在 Abaqus/Standard中CONTACT规则Message文件文件, Step 3, Increment 6 (contd): CONVERGENCE CHECKS FOR SEVERE D
15、ISCONTINUITY ITERATION 1 MAX. PENETRATION ERROR -8.16193E-009 AT NODE 331 OF CONTACT PAIR (ESURF,FSURF) MAX. CONTACT FORCE ERROR -4369.44 AT NODE 363 OF CONTACT PAIR (CSURF,DSURF) THE ESTIMATED CONTACT FORCE ERROR IS LARGER THAN THE TIME-AVERAGED FORCE. AVERAGE FORCE 5.393E+03 TIME AVG. FORCE 3.147E
16、+03 LARGEST RESIDUAL FORCE -1.110E+04 AT NODE 333 DOF 2 LARGEST INCREMENT OF DISP. -7.782E-04 AT NODE 329 DOF 2 LARGEST CORRECTION TO DISP. -1.737E-05 AT NODE 337 DOF 2 FORCE EQUILIBRIUM NOT ACHIEVED WITHIN TOLERANCE. AVERAGE MOMENT 114. TIME AVG. MOMENT 89.8 ALL MOMENT RESIDUALS ARE ZERO LARGEST IN
17、CREMENT OF ROTATION 1.853E-33 AT NODE 100 DOF 6 LARGEST CORRECTION TO ROTATION 6.489E-34 AT NODE 300 DOF 6 THE MOMENT EQUILIBRIUM EQUATIONS HAVE CONVERGED接触状态的收敛性检查平衡的收敛性检查不仅是接触的不相容太大,而且是力的平衡也没有达到 Dassault Systmes, 2008Introduction to Abaqus/Standard and Abaqus/ExplicitL6.15在 Abaqus/Standard中CONTACT
18、规则需要四次需要四次额外的迭代;前三次都是外的迭代;前三次都是SDIs(包含接触不相容)。(包含接触不相容)。 在最后的迭代中接触和平衡在最后的迭代中接触和平衡检查都通都通过,迭代收,迭代收敛。 CONVERGENCE CHECKS FOR SEVERE DISCONTINUITY ITERATION 2 . CONVERGENCE CHECKS FOR SEVERE DISCONTINUITY ITERATION 3 . CONVERGENCE CHECKS FOR SEVERE DISCONTINUITY ITERATION 4 . CONVERGENCE CHECKS FOR EQUIL
19、IBRIUM ITERATION 1 MAX. PENETRATION ERROR -1.38869E-014 AT NODE 331 OF CONTACT PAIR (ESURF,FSURF) MAX. CONTACT FORCE ERROR -0.00111133 AT NODE 331 OF CONTACT PAIR (ESURF,FSURF) THE CONTACT CONSTRAINTS HAVE CONVERGED. AVERAGE FORCE 5.244E+03 TIME AVG. FORCE 3.123E+03 LARGEST RESIDUAL FORCE -9.24 AT N
20、ODE 367 DOF 1 LARGEST INCREMENT OF DISP. -7.809E-04 AT NODE 129 DOF 2 LARGEST CORRECTION TO DISP. 4.229E-08 AT NODE 137 DOF 2 THE FORCE EQUILIBRIUM EQUATIONS HAVE CONVERGED AVERAGE MOMENT 109. TIME AVG. MOMENT 89.0 ALL MOMENT RESIDUALS ARE ZERO LARGEST INCREMENT OF ROTATION 1.925E-33 AT NODE 100 DOF
21、 6 LARGEST CORRECTION TO ROTATION 2.049E-35 AT NODE 100 DOF 6 THE MOMENT EQUILIBRIUM RESPONSE WAS LINEAR IN THIS INCREMENTNo SDIs in this iteration Dassault Systmes, 2008Introduction to Abaqus/Standard and Abaqus/ExplicitL6.16在 Abaqus/Standard中CONTACT规则增量摘要增量摘要: ITERATION SUMMARY FOR THE INCREMENT:
22、5 TOTAL ITERATIONS, OF WHICH 4 ARE SEVERE DISCONTINUITY ITERATIONS AND 1 ARE EQUILIBRIUM ITERATIONS. CURRENT VALUE OF MONITOR NODE 200 D.O.F. 2 IS -3.028E-03 TIME INCREMENT COMPLETED 1.266E-02, FRACTION OF STEP COMPLETED 5.047E-02 STEP TIME COMPLETED 5.047E-02, TOTAL TIME COMPLETED 2.05 Dassault Sys
23、tmes, 2008Introduction to Abaqus/Standard and Abaqus/ExplicitL6.17在 Abaqus/Standard中CONTACT规则Abaqus/Viewer中的接触中的接触诊断断约束的节点将要张开:不一致的接触状态选中这个来看模型中哪里的接触状态改变了 Dassault Systmes, 2008习题5:密封圈接触 Dassault Systmes, 2008Introduction to Abaqus/Standard and Abaqus/ExplicitL6.19习题5:密封圈接触习题任任务1.评估超弹材料模型2.定义面-面接触3.应用边界条件4.完成大位移分析5.显示结果