有限元非线性分析-正式课件-2011-01-06.ppt

* * * * * * * * * * * * * * * * * * * * * * * ANSYS求解非线性问题分析 ANSYS的求解器通过一系列的联立方程来预测工程系统的响应。一种近似的非线性求解是将载荷分成一系列的载荷增量，在每一个增量求解完毕并在进行下一个增量前，程序将调整刚度矩阵以反映结构刚度矩阵的非线性变化，纯粹的增量近似将不可避免的随着载荷增量的累积误差，导致结果失去平衡。 ANSYS使用牛顿-拉普森平衡迭代法(Newton-Raphson)克服这种困难，它迫使每一个载荷增量的末端解达到平衡收敛(在某个容限范围内) Newton-Raphson (N-R)迭代法 迭代求解的几个层次 分为求解层次： 载荷步: 顶层级别在一定时间范围内明确定义的载荷步组成，假定载荷在载荷步内是线性的，在每一个载荷步内，为了逐步加载，可以控制程序进行一系列的平衡迭代获得收敛的解，更多的子步通常会导致较好的精度，然而却需要更长的运行时间。 子步(Substep): 在每一个子步内，程序将进行一系列的平衡迭代以获得收敛的解。 平衡迭代：迭代不成功的时候就通过二分法减小时间步长，直到获得收敛为止。 控制收敛的方法 子步数或者时间步长: 用户可以通过指定实际的子步数也可以通过指定时间步长控制步数 自动时间步长: ANSYS基于结构的特征和系统的响应来调整时间步长 收敛容限：在平衡迭代过程中，用户要设置收敛容限，当迭代误差小于用户的设置值时，终止迭代。 分析问题类型的设置 步长控制选项 非线性控制选项 瞬态动力分析参数 ANSYS的弹塑性分析 前处理:（PREP7） 处理器:（Solution) 后处理: (Post1:通用的后处理器 ) (Post26: 时间-历史后处理器) 材料属性的输入 先输入弹性模量和泊松比，然后再输入塑性应变的表格 材料属性的输入-TC4材料参数 弹性模量:1.16E11Pa 泊松比：0.35 屈服强度：825MPa 屈服后的曲线斜率为1/4的弹性模量 耐压壳体2D模型应力校验结果 耐压壳体的3D端盖模型的应力校验 * * * * * * * * * * * * * * * * * * * 非线性结构分析简介 几何非线性（大应变、屈曲分析等） 材料非线性（弹塑性分析） 接触分析（高度非线性） ANSYS的设置 第一部分 非线性结构分析简介 什么是非线性结构分析 固体力学在分析线性弹性体系时，假设 节点位移无限小， 材料的应力应变满足胡克定律 加载时边界条件的性质保持不变， 若不满足上述条件之一的就是非线性问题。 结构分析的基本变量 在材料确定的情况下，基本的力学变量是位移、应变、应力 结构非线性产生的原因 固体力学问题从本质上讲是非线性的，引起非线性的原因很多，可以分成三种类型： 几何非线性：结构的变形使体系的受力状态发生显著改变，以致于不能用线性分析方法的非线性问题，比如大位移小应变问题、大位移大应变问题，结构变形引起的载荷大小、方向或者边界支撑条件的变化问题。 材料非线性：由于加载历史、环境状况（比如温度）及加载时间总量的影响使得材料的应力与应变关系不符合胡克定律的问题称为材料非线性问题，包括：弹塑分析、蠕变分析、超弹性分析 状态变化（包括接触）：许多普通结构表现出一种与状态相关的非线性行为，例如，轴承套可能是接触的，也可能是不接触的。这些系统的刚度由于系统状态的改变而突然变化，接触是一种很普遍的非线性行为。 非线性问题的求解 非线性分析的仍然三个步骤，但需要反复迭代。求解一般可以分为：增量法、迭代法和混合法，非线性分析时，结构的平衡实际上是在结构发生变形之后达到的。因此，分析的基本问题都是求出当前载荷作用下的平衡状态。 有限元离散系统的平衡方程为 Pt-ft =0 (Pt是t时刻节点外载荷，ft是t时刻节点力), 增量法求解的基本思想是:假设t时刻的解已知，在t+Δt时刻有P t+Δt –ft+Δt =0，由于t时刻的解已知，故ft+Δt =ft +f， 第二部分几何非线性 大应变和大转动 小的变形和小的应变分析假定位移小到足够使所得到的刚度改变无足轻重。这种刚度不变假定意味着使用基于最初几何形状的结构刚度的一次迭代足以计算出小变形分析中的位移。 大应变分析考虑由单元的形状和取向改变导致的刚度改变。因为刚度受位移影响，且反之亦然，所以在大应变分析中需要迭代求解来得到正确的位移。通过发出 NLGEOM ，ON（GUI路径Main MenuSolutionAnalysis Options)，来激活大应变效应。这种效应改变单元的形状和取向，且还随单元转动表面载荷。（集中载荷和惯性载荷保持它们最初的方向。）在大多数实体单元（包括所有的大应变和超弹性单元），以及部分的壳单元中大应变特性是可用的 大应变