第六讲：非线性分析.ppt

1、 Innovating through simulation北京怡格明思工程技术有限公司 第六讲 非线性分析 王慎平 北京怡格明思工程技术有限公司 Innovating through simulation北京怡格明思工程技术有限公司非线性有限元分析结构的非线性问题就是指结构的刚度随其变形而改变。所有的物理结构都是非线性的，而线性分析只是一种方便的近似，这对一些简单设计来说通常是精确的，但显然对许多结构模拟是不够的，诸如加工过程的模拟，锻造过程，冲压，压溃分析和橡胶问题等的分析。由于刚度依赖位移，所以不能再用初始柔度（将刚度阵集成并求逆即可得到柔度阵）乘以所加载荷的方法来计算任何载荷作用下的位

2、移。在非线性分析中，结构的刚度阵在分析过程中必须进行多次的集成和求逆，这就使得非线性分析求解比线性分析要昂贵得多。Innovating through simulation北京怡格明思工程技术有限公司非线性的来源与一般解法非线性的来源与一般解法1.材料非线性材料非线性非线性弹性非线性弹性弹塑性弹塑性超弹性超弹性粘弹性与粘塑性粘弹性与粘塑性 Innovating through simulation北京怡格明思工程技术有限公司2.几何非线性大偏转或变形；大扭曲；结构不稳定性(屈曲)预紧力效应 Innovating through simulation北京怡格明思工程技术有限公司3.边界非线性两个

3、物体的接触边界随加载和变形而两个物体的接触边界随加载和变形而改变引起的接触非线性（其中包含有改变引起的接触非线性（其中包含有摩擦接触和无摩擦接触）；摩擦接触和无摩擦接触）；Innovating through simulation北京怡格明思工程技术有限公司 求解非线性问题主要有两种方法隐式方法能够求解静态和动态方程；需要求解一组矩阵方程以便获得增量步结束时的状态；需要进行多次迭代；ABAQUS/standard应用该方法求解非线性问题；显式方法只能求解动态平衡方程；可以用来求解准静态问题；下一步的分析只依赖于上一步分析结束时的结果；不需要进行迭代求解；ABAQUS/Explicit应用显式方

4、法求解；Innovating through simulation北京怡格明思工程技术有限公司求解平衡方程典型的非线性问题具有所有三种形式的非线性。在方程中必须包括非线性项。一般的，每个自由度的非线性方程是耦合的。静态平衡的基本表达式为：由单元应力引起的加在节点上的内力，I，与外力，P，必须平衡，即：(Eq.3.1)Innovating through simulation北京怡格明思工程技术有限公司牛顿-拉普森求解技术第一次迭代(i=1)假定前面收敛增量步的解u0，P0，为已知的。在当前增量步中，将一个小的增量P，载荷施加到结构上。ABAQUS基于u0处的切线刚度K0 确定位移修正c1；AB

5、AQUS更新模型的状态为u1，形成K1 并计算 I1。总载荷PTOTAL与内力 I1的差称为残差，R1=PTOTAL-I1.如果 R1 在模型的每个自由度上都非常小(在容差范围之内)，结构就是平衡的。l默认的容差R1必须小于在整个时间段上作用于结构上的平均力的0.5%。lABAQUS自动计算这个在空间和时间上的平均力。如果迭代不能得到收敛的解，ABAQUS执行另外的迭代，以找到收敛的解。Innovating through simulation北京怡格明思工程技术有限公司第二次迭代(i=2)基于更新的刚度K1，计算新的位移纠正 c2，并且把新的残差R2与容差进行比较，察看在u2 处是否得到收敛

6、解。Innovating through simulation北京怡格明思工程技术有限公司该过程将一直重复，直到力的残差在允许的容差之内。每次迭代i需要：1.形成切线刚度Ki。2.求解系统方程组，得到位移修正ci+1。修正位移的估计值：ui+1=ui+ci+1。3.基于ui+1 计算内力向量Ii+1。4.进行平衡收敛判断：是否Ri+1 在容差之内?是否 Innovating through simulation北京怡格明思工程技术有限公司一般的，每个分析步(*STEP)需要几个增量步。Innovating through simulation北京怡格明思工程技术有限公司非线性输入文件 Inno

7、vating through simulation北京怡格明思工程技术有限公司非线性输入文件*HEADING CANTILEVER BEAM EXAMPLE-LARGE DISPLACEMENT*NODE 1,0.,0.11,200.,0.*NGEN1,11,1*ELEMENT,TYPE=B211,1,3*ELGEN,ELSET=BEAMS1,5,2,1*BEAM SECTION,SECTION=RECT,ELSET=BEAMS,MATERIAL=MAT150.,5.*MATERIAL,NAME=MAT1*ELASTIC2.E5,.3*BOUNDARY1,1,6*AMPLITUDE,NAME=

8、RAMP0.0,0.0,0.5,0.3,1.0,1.0 Innovating through simulation北京怡格明思工程技术有限公司*RESTART,WRITE,FREQ=3*STEP,NLGEOM,INC=25APPLY POINT LOAD*STATIC0.1,1.0,0.001,1.0*CLOAD,AMPLITUDE=RAMP11,2,-1200.*NODE PRINT,FREQ=1U,RF*EL PRINT,FREQ=10S,E*NODE FILE,FREQ=5U*END STEP 和线性输入的主要不同最小时间增量最大时间增量分析步的时间周期建议的初始时间增量先前定义的载荷幅

9、值函数和线性输入主要的不同非线性输入文件 Innovating through simulation北京怡格明思工程技术有限公司分析步和过程输入*STEP,NLGEOM,INC=25NLGEOM:包括所有由以下原因引起的几何非线性效应：大挠度、大旋转、大变形。预载荷（初始应力）。载荷刚度。如果上面列出的项不重要，应用NLGEOM选项得到的结果同没有应用NLGEOM选项得到的结果类似，但是求解的费用更高。INC=25:在本例中允许的最大增量为25：如果在施加全部载荷之前达到了最大增量数量，程序将会中止。保证程序的运行时间不会太长用户可以重新启动分析。默认值为100。非线性输入文件 Innovat

10、ing through simulation北京怡格明思工程技术有限公司非线性悬臂梁分析的输出状态(.sta)文件总结分析的过程显示计算过程中使用的自动时间增量。在作业运行的同时，可以检查状态文件。在每个成功的增量之后，记录一行。SUMMARY OF JOB INFORMATION: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 0 3 3 0.100 0.100 0.1000 1 2

11、1 0 2 2 0.200 0.200 0.1000 1 3 1 0 2 2 0.350 0.350 0.1500 1 4 1 0 2 2 0.575 0.575 0.2250 1 5 1 0 3 3 0.913 0.913 0.3375 1 6 1 0 2 2 1.00 1.00 0.08750 Innovating through simulation北京怡格明思工程技术有限公司自动的时间增量试探算法（基于多年的经验）控制时间积分的精度。在静力学分析中，基于迭代次数达到收敛。容易得到收敛解（比最大允许迭代数量少很多）：增加迭代步长不容易收敛或发散：减小增量步长否则：保持同样的增量步长自动的

12、时间增量步长控制工作很好。如果没有特殊的原因，用户不要改变它。提示：提示：对于高度非线性问题，推荐初始的时间增量为总时间增量的一小部分（比如10%）。非线性悬臂梁分析的输出 Innovating through simulation北京怡格明思工程技术有限公司信息(.msg)文件包括：所有的收敛控制：利用*CONTROLS选项覆盖默认值不是经常需要关于具体模型特征的细节：非默认模型特征使用NLGEOM参数重启动文件的写出频率所有的迭代细节非线性悬臂梁分析的输出 Innovating through simulation北京怡格明思工程技术有限公司常见的收敛性问题单元扭曲（Element dis

13、tortion）当一个单元在某一积分点的体积变为负值时，会 提示该信息。检查属性、载荷、边界条件，一般来说需要对网格进行细化；过度屈服（Excessive yielding）在一个大变形问题中，如果材料属性包含一定 的塑性，会出现该信息。检查载荷（是否过载？），材料定义或者对网 格进行细化；沙漏问题（Hourglassing）在应用一次减缩积分单元时产生的非物理变形。如果伪沙漏刚度过大，可以采用细化网格或者应用其他类型的单元类型 （可以应用一次全积分单元或者二次单元）；Innovating through simulation北京怡格明思工程技术有限公司某些在ABAQUS/Standard出现

14、的收敛性问题数值奇异（Numerical singularities）这个信息通常说明在一个线性方程求解中丢失了太多的数据信息，所 得到的分析结果是不可靠的。最可能的原因是在进行静态应力分析时 对刚体没有进行约束；零 pivots（Zero pivots）最可能的原因是存在未约束的刚体和过约束的自由度；负特征值（Negative eigenvalues）负特征值一般说明刚度矩阵不一定是正值。Innovating through simulation北京怡格明思工程技术有限公司某些在ABAQUS/Explicit中出现的收敛性问题超过波前速度（Excessive wave speed）这个错误信

15、息说明模型中的某种因素导致了一个或多个单元在一个增量 步中极度变形。这个信息很有可能与上一条刚出现的的警告信息相关。可能产生错误的原因是分析过程中包含了接触，不稳定的时间增量（用 户调整了时间增量或直接对时间增量进行控制），网格扭曲和沙漏；网格扭曲（Mesh distortion）重新细分网格或者应用自适应网格；沙漏 （Hourglassing）细分网格或者应用非缺省沙漏控制模式；ABAQUS/Explicit 提供了几种形式的沙漏控制模式；Innovating through simulation北京怡格明思工程技术有限公司 S T E P 1 S T A T I C A N A L Y S

16、 I S SMALL DISPLACEMENT ANALYSIS AUTOMATIC TIME CONTROL WITH-A SUGGESTED INITIAL TIME INCREMENT OF 0.100 AND A TOTAL TIME PERIOD OF 1.00 THE MINIMUM TIME INCREMENT ALLOWED IS 1.000E-03 THE MAXIMUM TIME INCREMENT ALLOWED IS 1.00 CONVERGENCE TOLERANCE PARAMETERS FOR FORCE CRITERION FOR RESIDUAL FORCE FOR A NONLINEAR PROBLEM 5.000E-03 CRITERION FOR DISP.CORRECTION IN A NONLINEAR PROBLEM 1.000E-02 INITIAL VALUE OF TIME AVERAGE FORCE 1.000E-02 AVERAGE FORCE IS TIME AVERAGE FORCE ALTERNATE CRIT.FOR RE