ANSYS高级非线性培训手册第09章课程.ppt

July 30, 1999 第九章 几何不稳定性 几何不稳定性本章综述 本章阐述几何不稳定性问题, 即关于屈曲的问题，将介绍以下技术: 特征值屈曲 载荷控制 位移控制 弧长法 几何不稳定性 ... 本章综述 本章包括以下主题: A. 结构稳定性背景 B. 线性(特征值)屈曲过程 C. 非线性屈曲技术背景 D. 非线性前屈曲过程 E. 非线性后屈曲过程 几何不稳定性A. 结构稳定性背景 很多结构需要评价它们的结构稳定性，细柱体、压杆和真空罐都是稳定性非常重要的结构的例子。 在不稳定性(屈曲)的开始, 在载荷没有实质性变化的情况下(除了一个小的载荷扰动), 结构的位移将有一个非常大的变化{?u}。 几何不稳定性 ... 结构稳定性背景 当增加轴向载荷(F)时, 一个理想化的端部固定的柱体将呈现下述行为。 几何不稳定性 ... 结构稳定性背景 分叉点 分叉点 是载荷历程中的一点, 该点可能存在两个分支解。 在理想化的端部固定柱体的情况下, 在临界载荷(Fcr)下, 柱体可向左或向右屈曲，因此可能存在两个载荷路径。在实际结构中, 几何缺陷的存在或力的扰动(P ? 0) 将决定载荷路径的方向。 几何不稳定性 ... 结构稳定性背景 稳定、不稳定及中性平衡 考虑下图所示球的平衡，若表面向上凹, 平衡是稳定的, 扰动时, 球返回初始位置。若表面向下凹, 平衡是不稳定的, 扰动时, 球将滚开。若表面是平的, 球处于中性平衡, 扰动时, 钢球将保持在新的位置。 几何不稳定性 ... 结构稳定性背景 临界载荷 当 F Fcr 时, 柱体处于稳定平衡状态，若引入一个小的扰动力(P ? 0), 然后卸载, 柱体将返回到它的初始位置。当 F Fcr 时, 柱体处于不稳定平衡状态, 任何扰动力将引起坍塌。 当 F = Fcr 时, 柱体处于中性平衡状态，把这个力定义为临界载荷。 几何不稳定性 ... 结构稳定性背景 极限载荷 在实际结构中, 很难达到临界载荷，因为扰动和非线性行为, 低于临界载荷时结构通常变得不稳定。 几何不稳定性B. 线性特征值屈曲 前屈曲和坍塌载荷分析的分析技术包括: 线性特征值屈曲 非线性屈曲分析 本节主要讨论第一种方法--线性特征值屈曲。 几何不稳定性 ... 线性特征值屈曲 特征值屈曲分析 预测一个理想线弹性 结构的理论屈曲强度(分叉点) 特征值公式决定结构的分叉点，该方法与线弹性屈曲分析的教科书所述方法一致。Euler 柱体的特征值屈曲解与经典Euler 解吻合。 几何不稳定性 ... 线性特征值屈曲 然而, 缺陷和非线性行为阻止大多数实际结构达到理想的弹性屈曲强度，特征值屈曲一般产生非保守 解, 使用时应谨慎。 几何不稳定性 ... 线性特征值屈曲 尽管特征值屈曲一般产生非保守的结果, 线性屈曲分析仍有两个优点: 相对不费时(快捷)的分析。 为了提供更真实的结果, 屈曲模态形状可用作非线性屈曲分析的初始几何缺陷。 几何不稳定性 ... 线性特征值屈曲 线性屈曲分析基于经典的特征值问题。为推导特征值问题, 首先求解线弹性前屈曲载荷状态 {P0} 的载荷-位移关系，即给定 {P0} 求解 {P0} = [Ke]{u0} 得到 {u0} = 施加载荷 {P0} 的位移结果 {s} = 与{u0}对应的应力 几何不稳定性 ... 线性特征值屈曲 假设前屈曲位移很小, 在任意 状态下({P}, {u}, {s}) 增量平衡方程由下式给出 {DP} = [[Ke] + [Ks(s)]]{Du} 式中 [Ke] = 弹性刚度矩阵 [Ks(s)] = 某应力状态 {s} 下计算的初始应力矩阵 几何不稳定性 ... 线性特征值屈曲 假设前屈曲行为是一个外加载荷 {P0} 的线性函数, {P} = l{P0} {u} = l{u0} {s} = l{s0} 则可得 [Ks(s)] = l[Ks(s0)] 因此, 整个前屈曲范围 内的增量平衡方程变为 {DP} = [[Ke] + l[Ks(s0)]]{Du} 几何不稳定性 ... 线性特征值屈曲 在不稳定性开始 (屈曲载荷{Pcr}) 时, 在 {DP} ? 0 的情况下, 结构会出现一个变形 {Du}。 [[Ke] + l[Ks(s0)]]{Du} = {DP} 把上述表达式 ({DP} ? 0) 代入前面的前屈曲范围内 的增量平衡方程, 则有 [[Ke] + l[Ks(s0)]]{Du} = {0} 上述关系代表经典的特征值问题。 几何不稳定性 ... 线性特征值屈曲 为了满足前面的关系, 必须有: det[[Ke] + l[Ks(s0)]]