Lecture3-Nonlinearity

Abaqus/Standard中的非线性分析第三讲L3.2IntroductiontoAbaqus/StandardandAbaqus/Explicit概述•非线性结构力学•求解平衡方程•非线性输入文件•非线性悬臂梁分析的输出•习题3：非线性静力学非线性结构力学L3.4IntroductiontoAbaqus/StandardandAbaqus/Explicit非线性结构力学一些材料非线性的例子•非线性的来源•材料非线性：•非线性弹性塑性•材料损伤•失效机制•等等。•注意：如果预定义了温度或场变量，则材料与温度或场变量的相关性不引入非线性。L3.5IntroductiontoAbaqus/StandardandAbaqus/Explicit•边界条件非线性：•接触问题•在分析过程中边界条件变化。•严重不连续形式的非线性.非线性结构力学自接触的例子：例子手册1.1.17，压缩减震器L3.6IntroductiontoAbaqus/StandardandAbaqus/Explicit非线性结构力学•几何非线性：•大绕度和大变形•大旋转•结构不稳定（屈曲）•预载荷效应几何非线性例子：聚合物键盘罩求解平衡方程L3.8IntroductiontoAbaqus/StandardandAbaqus/Explicit求解平衡方程•典型的非线性问题具有所有三种形式的非线性。•在方程中必须包括非线性项。•一般的，每个自由度的非线性方程是耦合的。•静态平衡的基本表达式为：由单元应力引起的加在节点上的内力I，与外力P，必须平衡，即：0)()(uIuP(Eq.3.1)L3.9IntroductiontoAbaqus/StandardandAbaqus/Explicit求解平衡方程L3.10IntroductiontoAbaqus/StandardandAbaqus/Explicit求解平衡方程•为了求解非线性平衡问题，在Abaqus/Standard中使用基于牛顿-拉普森技术的增量迭代法。•假定前一步载荷增量的解，u0，为已知。•假定在第i次迭代之后，得到近似解ui。设ci+1为离散平衡方程的精确解和当前解之差（Equation3.1），所以有：•将方程3.2的左手边在近似解ui附近，以泰勒级数方式展开，可以得到(Eq.3.2)(Eq.3.3).0)()(11iiiicuIcuP1()()()()....0iiiiiPuIuPuIucuu.L3.11IntroductiontoAbaqus/StandardandAbaqus/Explicit•略去高阶项，方程可以写为其中为切线刚度。•解的下一次近似为•注意：如果载荷与位移相关（比如，旋转的表面压力），刚度矩阵中包含载荷刚度的贡献。求解平衡方程,)()(1iiiiuIuPcKuuPuuIKiii)()(.11iiicuuL3.12IntroductiontoAbaqus/StandardandAbaqus/Explicit求解平衡方程•网格中的静力平衡1.施加规定的“载荷增量”。2.迭代，直到每个节点的所有节点力的和非常小。3.满足平衡后，更新模型状态。4.回到第1步，施加下一个载荷增量。L3.13IntroductiontoAbaqus/StandardandAbaqus/Explicit•单自由度例子•非线性弹簧•求解u(P)或P(u)—典型的，在分析步中，载荷从0递增到PFINAL。•时间经常从0变到1。求解平衡方程L3.14IntroductiontoAbaqus/StandardandAbaqus/Explicit求解平衡方程•牛顿-拉普森求解技术•第一次迭代(i=1)•假定前面收敛增量步的解u0，P0，为已知的。•在当前增量步中，将一个小的增量P，载荷施加到结构上。•Abaqus基于u0处的切线刚度K0确定位移修正c1；前一增量步结束时，总载荷PTOTAL和内力间的关系为：.:0100IPcKITOTALL3.15IntroductiontoAbaqus/StandardandAbaqus/Explicit•Abaqus更新模型的状态为u1，形成K1并计算I1。•总载荷PTOTAL与内力I1的差称为残差，R1:R1=PTOTALI1.•如果R1在模型的每个自由度上都非常小(在容差范围之内)，结构就是平衡的。•默认的容差R1必须小于结构对时间平均力的0.5%。•Abaqus自动计算时间平均力。•如果迭代不能得到收敛的解，Abaqus执行另外的迭代，以找到收敛的解。求解平衡方程L3.16IntroductiontoAbaqus/StandardandAbaqus/Explicit•第二次迭代(i=2)•基于更新的刚度K1，计算新的位移纠正c2，并且•把新的残差R2与容差进行比较，察看在u2处是否得到收敛解，。求解平衡方程.1211:IPcKITOTALL3.17IntroductiontoAbaqus/StandardandAbaqus/Explicit求解平衡方程•该过程将一直重复，直到力的残差在允许的容差之内。每次迭代i需要：1.形成切线刚度Ki。2.求解系统方程组，得到位移修正ci+1。•修正位移的估计值：ui+1=ui+ci+1。3.基于ui+1计算内力向量Ii+1。4.进行平衡收敛判断：•是否Ri+1在容差之内?•是否#11iterijjccu?L3.18IntroductiontoAbaqus/StandardandAbaqus/Explicit求解平衡方程•一般的，每个分析步(*STEP)需要几个增量步。非线性输入文件L3.20IntroductiontoAbaqus/StandardandAbaqus/Explicit非线性输入文件*HEADINGCANTILEVERBEAMEXAMPLE--LARGEDISPLACEMENT*NODE1,0.,0.11,200.,0.*NGEN1,11,1*ELEMENT,TYPE=B211,1,3*ELGEN,ELSET=BEAMS1,5,2,1*BEAMSECTION,SECTION=RECT,ELSET=BEAMS,MATERIAL=MAT150.,5.*MATERIAL,NAME=MAT1*ELASTIC2.E5,.3*BOUNDARY1,1,6*AMPLITUDE,NAME=RAMP0.0,0.0,0.5,0.3,1.0,1.0L3.21IntroductiontoAbaqus/StandardandAbaqus/Explicit*RESTART,WRITE,FREQ=3*STEP,NLGEOM,INC=25APPLYPOINTLOAD*STATIC0.1,1.0,0.001,1.0*CLOAD,AMPLITUDE=RAMP11,2,-1200.*NODEPRINT,FREQ=1U,RF*ELPRINT,FREQ=10S,E*NODEFILE,FREQ=5U*ENDSTEP和线性输入的主要不同最小时间增量最大时间增量分析步的时间周期建议的初始时间增量先前定义的载荷幅值函数和线性输入主要的不同非线性输入文件L3.22IntroductiontoAbaqus/StandardandAbaqus/Explicit•分析步和过程输入*STEP,NLGEOM,INC=25•NLGEOM:包括所有由以下原因引起的几何非线性效应：•大挠度、大旋转、大变形。•预载荷（初始应力）。•载荷刚度。•如果上面列出的项不重要，应用NLGEOM选项得到的结果同没有应用NLGEOM选项得到的结果类似，但是求解的费用更高。•INC=25:在本例中允许的最大增量为25：•如果在施加全部载荷之前达到了最大增量数量，程序将会中止。•保证程序的运行时间不会太长—用户可以重新启动分析。•默认值为100。非线性输入文件L3.23IntroductiontoAbaqus/StandardandAbaqus/Explicit非线性输入文件•对于所有的瞬态过程，都有类似的时间增量，它们是*STATIC*DYNAMIC*HEATTRANSFER*VISCO*COUPLEDTEMPERATURE-DISPLACEMENT*SOILS*MODALDYNAMIC(只允许固定时间增量)*COUPLEDTHERMAL-ELECTRIC•物理（真实）时间或名义时间的度量取决于过程和时间相关过程或率相关行为。非线性悬臂梁分析的输出L3.25IntroductiontoAbaqus/StandardandAbaqus/Explicit非线性悬臂梁分析的输出•状态(.sta)文件•总结分析的过程—显示计算过程中使用的自动时间增量。•在作业运行的同时，可以检查状态文件。•在每个成功的增量之后，记录一行。SUMMARYOFJOBINFORMATION:STEPINCATTSEVEREEQUILTOTALTOTALSTEPINCOFDOFIFDISCONITERSITERSTIME/TIME/LPFTIME/LPFMONITORRIKSITERSFREQ1110330.1000.1000.10001210220.2000.2000.10001310220.3500.3500.15001410220.5750.5750.22501510440.9130.9130.33751610221.001.000.08750THEANALYSISHASCOMPLETEDSUCCESSFULLYL3.26IntroductiontoAbaqus/StandardandAbaqus/Explicit•自动的时间增量•试探算法（基于多年的经验）控制时间积分的精度。•在静力学分析中，基于迭代次数达到收敛。•容易得到收敛解（比最大允许迭代数量少很多）：增加迭代步长•不容易收敛或发散：减小增量步长•否则：保持同样的增量步长•自动的时间增量步长控制工作很好。如果没有特殊的原因，用户不要改变它。提示：对于高度非线性问题，推荐初始的时间增量为总时间增量的一小部分（比如10%）。非线性悬臂梁分析的输出L3.27IntroductiontoAbaqus/StandardandAbaqus/Explicit•信息(.msg)文件•包括：•所有的收敛控制：•利用*CONTROLS选项覆盖默认值—不是经常需要•关于具体模型特征的细节：•非默认模型特征•使用NLGEOM参数•重启动文件的写出频率•所有的迭代细节非线性悬臂梁分析的输出L3.28IntroductiontoAbaqus/StandardandAbaqus/Explicit非线性悬臂梁分析的输出•求解器信息：•数值奇异：这表明在线性方程组的求解过程中造成了大量数位的丢失，以至于结果是不可靠的。最常见的原因是静力应力分析过程中没有约束的刚体模式。•零主元：在线性方程组求解过程中，当没有相应的刚度抵抗所施加的力时，会出现这样的信息。最常见的原因是没有约束的刚体模式和过度约束的自由度。•负特征值：负特征值表明刚度矩阵是非正定的；比如，超过结构的屈曲载荷。L3.29IntroductiontoAbaqus/StandardandAbaqus/Explicit非线性悬臂梁分析的输出•有用的修理故障信息：•最大残差的位置•过度变形的位置•接触条件变化的位置L3.30IntroductiontoAbaqus/StandardandAbaqus/ExplicitSTEP1STATICANALYSISSMALLDISPLACEMENTANALYSISAUTOMATICTIMECONTROLWITH-ASUGGESTEDINITIALTIMEINCREMENTOF0.100ANDATOTALTIMEPERIODOF1.00THEMINIMUMTIMEINCREMENTALLOWEDIS1.000E-03THEMAXIMUMTIMEINCREMENTALLOWEDIS1.00CONVERGE