03 非线性分析.docVIP

第三部分 非线性分析 非线性有限元概述 1.1非线性行为 非线性结构的基本特征是结构刚度随载荷的改变而变化。如果绘制一个非线性结构的载荷－位移曲线，则力与位移的关系是非线性函数。 引起结构非线性的原因几何非线性大应变，大位移，大旋转 b 材料非线性 塑性，超弹性，粘弹性，蠕变 状态改变非线性 接触，单元死活 非线性行为——不能使用叠加原理! 结构响应与路径有关，也就是说加载的顺序可能是重要的。 结构响应与施加的载荷可能不成比例。 非线性分析的应用 一些典型的非线性分析的应用包括： 非线性屈曲失稳分析 金属成形研究 碰撞与冲击分析 制造过程分析（ 装配、部件接触等） 材料非线性分析 （塑性材料、聚合物） 橡胶底密封 一个包含几何非线性（大应变与大变形），材料非线性（橡胶），及状态非线性（接触的例子。 求解一个结构的平衡问题通常等于求解结构的总位能的驻值问题。结构总位能 增量法：就是将荷载分成一系列的荷载增量，即 A 要点：在每一个荷载增量求解完成后，继续进行下一个荷载增量之前，调整刚度矩阵以反映结构刚度的变化。 B 增量法的优点：可以追踪结构变形历程，这对于材料或几何非线性（特别是极限值屈曲分析）十分有用。 增量法的缺点：随着荷载步增量的增加而产生积累误差，导致荷载－位移曲线飘移。 对飘移进行平衡修正，可以大大提高增量法的精度。应用最广的就是在每一级载荷增量上用 3、迭代法割线刚度法收敛性差，因此很少应用 切线刚度法 Newton-Raphson 2.2 Newton-Raphson迭代法 优点：对于一致的切向刚度矩阵有二次收敛速度。 Newton-Raphson 法迭代求解使用下列方程 [KT] = 切向刚度矩阵 {u} = 位移增量 {Fa} = 施加的载荷矢量 {Fnr} = 内力矢量 Newton-Raphson 法需要一个收敛的度量以决定何时结束迭代。 给定外部载荷（Fa），内部载荷（ Fnr ）（由单元应力产生并作用于节点），在一个体中，外部载荷必须与内力相平衡。 收敛是平衡的度量。 迭代中的Newton-Raphson 不平衡量是： 不平衡量 R 实际上从未真正等于零。当不平衡量小到误差允许范围内时，可中止迭代，得到平衡解。 虽然使用一致切向刚度的Newton-Raphson法具有平方的收敛速度，但它不能保证一定收敛！只有初始构形在收敛半径以内， Newton-Raphson 才可以保证收敛。Newton-Raphson迭代法： A F.N.R: 缺点：消耗机时！ 优点：适用于高度非线性问题。 I.N.R: 缺点：收敛性较差。 优点：消耗机时少。 M.N.R: 优缺点介于前两者之间 计算效率： 2.3 收敛准则 1、Newton-Raphson迭代法的收敛判据： 位移/旋转增量准则 位移收敛准则有时不可靠！ 力/力矩平衡准则 不平衡力收敛准则有时也不可靠！ 能量准则(ansys中不使用，用户可自己定义) 2.4 载荷步、子步与平衡迭代步 1、载荷步： 在ANSYS中，结构上施加的载荷及边界条件由一系列定义的载荷步来描述。 子步： 给定载荷步中的载荷是逐步施加上去的 ，载荷的每个增量称之为子步。 平衡迭代步：每个子步中为得到收敛解而进行的迭代步。 每个载荷步与子步都与“时间”相关联。子步也叫时间步。 在率相关分析（蠕变，粘塑性）与瞬态分析中，“时间 ”代表真实的时间。 对于率无关的静态分析，“时间” 表示加载次序。在静态分析中，“时间”可设置为任何适当的值。 建模技巧: 在静态分析中，“ 时间 ”可设置为给定载荷的大小，这样将易于绘制载荷－位移曲线。 2.5 对于非线性分析，有三种方程求解器供选择：Sparse 波前求解器 PCG 2、梁/壳模型，或梁/壳/实体模型：使用sparse 求解器。 三维实体模型（ Solid92 或 Solid45 ），自由度数相对较大（ 100,000 ）：使用PCG 求解器 。 病态问题或单元刚度矩阵带宽大（包含在输出文件中）：使用sparse 求解器。 非对称矩阵：使用sparse 求解器。 注意: 如果可使用并行处理，波前求解器可能比sparse 求解器速度快，因为波前求解器对并行计算进行了优化 收敛范数L1范数，L2范数，无限范数 L1范数:用不平衡量绝对值的和与收敛判据作对比。 L2范数:用力不平衡量的 SRSS （平方和的平方根）。 无限范数:检查所有自由度的最大不平衡量。（此选项的作用是为收敛独立检查模型的每个自由度。） 若响应平稳且时间步长适当小，预测可以加速收敛。 若响应不平稳或分析中涉及大旋转，预测可能导致发散！ 对大旋转分析不要使用预测。 2.6