ansys高级非线性分析—几何不稳定性

1、几何不稳定性几何不稳定性September 30, 2001Inventory #0014919-2几何不稳定性几何不稳定性本章综述本章综述本章阐述几何不稳定性问题本章阐述几何不稳定性问题, 即关于屈曲的问题，将介绍以下技术即关于屈曲的问题，将介绍以下技术: 特征值屈曲特征值屈曲 载荷控制载荷控制 位移控制位移控制 弧长法弧长法September 30, 2001Inventory #0014919-3几何不稳定性几何不稳定性 . 本章综述本章综述本章包括以下主题本章包括以下主题:A. 结构稳定性背景结构稳定性背景B. 线性线性(特征值特征值)屈曲过程屈曲过程 C. 非线性屈曲技术背景非线性屈

2、曲技术背景D. 非线性前屈曲过程非线性前屈曲过程E. 非线性后屈曲过程非线性后屈曲过程September 30, 2001Inventory #0014919-4几何不稳定性几何不稳定性A. 结构稳定性背景结构稳定性背景很多结构需要评价它们的结构稳定性，细柱体、压杆和真空罐都是很多结构需要评价它们的结构稳定性，细柱体、压杆和真空罐都是稳定性非常重要的结构的例子。稳定性非常重要的结构的例子。在不稳定性在不稳定性(屈曲屈曲)的开始的开始, 在载荷没有实质性变化的情况下在载荷没有实质性变化的情况下(除了一除了一个小的载荷扰动个小的载荷扰动), 结构的位移将有一个非常大的变化结构的位移将有一个非常大的

3、变化 u。FF稳定稳定不稳定不稳定September 30, 2001Inventory #0014919-5几何不稳定性几何不稳定性 . 结构稳定性背景结构稳定性背景当增加轴向载荷当增加轴向载荷(F)时时, 一个理想化的端部固定的柱体将呈现下述行一个理想化的端部固定的柱体将呈现下述行为。为。uF 分叉点分叉点稳定平衡稳定平衡中性平衡中性平衡不稳定平衡不稳定平衡FcrFFuSeptember 30, 2001Inventory #0014919-6几何不稳定性几何不稳定性 . 结构稳定性背景结构稳定性背景分叉点分叉点分叉点分叉点 是载荷历程中的一点是载荷历程中的一点, 该点可能存在两个分支解。

4、该点可能存在两个分支解。在理想化的端部固定柱体的情况下在理想化的端部固定柱体的情况下, 在临界载荷在临界载荷(Fcr)下下, 柱体可向左或向右屈曲，因此可能存在两个载荷路径。柱体可向左或向右屈曲，因此可能存在两个载荷路径。 在实际结构中在实际结构中, 几何缺陷的存在或力的扰动几何缺陷的存在或力的扰动(P 0) 将决定载荷路径的方向。将决定载荷路径的方向。FFuPSeptember 30, 2001Inventory #0014919-7几何不稳定性几何不稳定性 . 结构稳定性背景结构稳定性背景稳定、不稳定及中性平衡稳定、不稳定及中性平衡考虑下图所示球的平衡，若表面向上凹考虑下图所示球的平衡，若

5、表面向上凹, 平衡是稳定的平衡是稳定的, 扰动时扰动时, 球返球返回初始位置。若表面向下凹回初始位置。若表面向下凹, 平衡是不稳定的平衡是不稳定的, 扰动时扰动时, 球将滚开。若球将滚开。若表面是平的表面是平的, 球处于中性平衡球处于中性平衡, 扰动时扰动时, 钢球将保持在新的位置。钢球将保持在新的位置。 稳定稳定 不稳定不稳定 中性中性September 30, 2001Inventory #0014919-8几何不稳定性几何不稳定性 . 结构稳定性背景结构稳定性背景临界载荷临界载荷当当 F Fcr 时时, 柱体处柱体处于不稳定平衡状态于不稳定平衡状态, 任何扰动力将引起坍塌。任何扰动力将引

6、起坍塌。 当当 F = Fcr 时时, 柱体柱体处于中性平衡状态，把这个力定义为临界载荷。处于中性平衡状态，把这个力定义为临界载荷。September 30, 2001Inventory #0014919-9几何不稳定性几何不稳定性 . 结构稳定性背景结构稳定性背景极限载荷极限载荷在实际结构中在实际结构中, 很难达到临界载荷，因为扰动和非线性行为很难达到临界载荷，因为扰动和非线性行为, 低于低于临界载荷时结构通常变得不稳定。临界载荷时结构通常变得不稳定。uF 分叉点分叉点Fcr实际的结构响应实际的结构响应, 低于临界载低于临界载荷时出现不稳定性。荷时出现不稳定性。September 30, 2

7、001Inventory #0014919-10几何不稳定性几何不稳定性B. 线性特征值屈曲线性特征值屈曲前屈曲和坍塌载荷分析的分析技术包括前屈曲和坍塌载荷分析的分析技术包括: 线性特征值屈曲线性特征值屈曲 非线性屈曲分析非线性屈曲分析本节主要讨论第一种方法本节主要讨论第一种方法-线性特征值屈曲。线性特征值屈曲。Fu理想载荷路径理想载荷路径有缺陷结构的载荷路径有缺陷结构的载荷路径 前屈曲前屈曲线性特征线性特征 值屈曲值屈曲非线性屈曲非线性屈曲September 30, 2001Inventory #0014919-11几何不稳定性几何不稳定性 . 线性特征值屈曲线性特征值屈曲特征值屈曲分析特征

8、值屈曲分析 预测一个预测一个理想线弹性理想线弹性 结构的理论屈曲强度结构的理论屈曲强度(分叉点分叉点)特征值公式决定结构的分叉点，该方法与线弹性屈曲分析的教科书特征值公式决定结构的分叉点，该方法与线弹性屈曲分析的教科书所述方法一致。所述方法一致。Euler 柱体的特征值屈曲解与经典柱体的特征值屈曲解与经典Euler 解吻合。解吻合。September 30, 2001Inventory #0014919-12几何不稳定性几何不稳定性 . 线性特征值屈曲线性特征值屈曲然而然而, 缺陷和非线性行为阻止大多数实际结构达到理想的弹性屈曲缺陷和非线性行为阻止大多数实际结构达到理想的弹性屈曲强度，特征值屈

9、曲一般产生强度，特征值屈曲一般产生非保守非保守 解解, 使用时应谨慎。使用时应谨慎。理想载荷路径理想载荷路径有缺陷结构的载荷路径有缺陷结构的载荷路径Fu 前屈曲前屈曲 分叉点分叉点极限载荷极限载荷September 30, 2001Inventory #0014919-13几何不稳定性几何不稳定性 . 线性特征值屈曲线性特征值屈曲尽管特征值屈曲一般产生非保守的结果尽管特征值屈曲一般产生非保守的结果, 线性屈曲分析仍有两个优线性屈曲分析仍有两个优点点: 相对不费时相对不费时(快捷快捷)的分析。的分析。 为了提供更真实的结果为了提供更真实的结果, 屈曲模态形状可用作非线性屈曲分析的初始屈曲模态形状

10、可用作非线性屈曲分析的初始几何缺陷。几何缺陷。September 30, 2001Inventory #0014919-14几何不稳定性几何不稳定性 . 线性特征值屈曲线性特征值屈曲线性屈曲分析基于经典的特征值问题。为推导特征值问题线性屈曲分析基于经典的特征值问题。为推导特征值问题, 首先求首先求解线弹性前屈曲载荷状态解线弹性前屈曲载荷状态 P0 的载荷的载荷-位移关系，即给定位移关系，即给定 P0 求解求解P0 = Keu0得到得到u0 =施加载荷施加载荷 P0 的位移结果的位移结果s s =与与u0对应的应力对应的应力September 30, 2001Inventory #0014919

11、-15几何不稳定性几何不稳定性 . 线性特征值屈曲线性特征值屈曲假设假设前屈曲位移很小前屈曲位移很小, 在在任意任意 状态下状态下(P, u, s s) 增量平衡方程增量平衡方程由下式给出由下式给出 P = Ke + Ks s(s s) u式中式中Ke = 弹性刚度矩阵弹性刚度矩阵Ks s(s s)=某应力状态某应力状态 s s 下计算的下计算的初始应力矩阵初始应力矩阵September 30, 2001Inventory #0014919-16几何不稳定性几何不稳定性 . 线性特征值屈曲线性特征值屈曲假设前屈曲行为是一个外加载荷假设前屈曲行为是一个外加载荷 P0 的线性函数的线性函数,P =

12、 l lP0u = l lu0s s = l ls s0 则可得则可得Ks s(s s) = llKs s(s s0 0)因此因此, 整个前屈曲范围整个前屈曲范围 内的增量平衡方程变为内的增量平衡方程变为 P = Ke + l lKs s(s s0 0) uSeptember 30, 2001Inventory #0014919-17几何不稳定性几何不稳定性 . 线性特征值屈曲线性特征值屈曲在在不稳定性开始不稳定性开始 (屈曲载荷屈曲载荷Pcr) 时时, 在在 P 0 的情况下的情况下, 结构会结构会出现一个变形出现一个变形 u。把上述表达式把上述表达式 ( P 0) 代入前面的代入前面的前屈

13、曲范围内前屈曲范围内 的增量平衡方程的增量平衡方程, 则有则有Ke + l lKs s(s s0 0) u = 0上述关系代表上述关系代表经典的特征值问题。经典的特征值问题。September 30, 2001Inventory #0014919-18几何不稳定性几何不稳定性 . 线性特征值屈曲线性特征值屈曲为了满足前面的关系为了满足前面的关系, 必须有必须有:detKe + l lKs s(s s0 0) = 0在在 n 个自由度的有限元模型中个自由度的有限元模型中, 上述方程产生上述方程产生 l l (特征值特征值) 的的 n阶多阶多项式，这种情况下特征向量项式，这种情况下特征向量 un

14、表示屈曲时叠加到系统上的变形表示屈曲时叠加到系统上的变形，由计算出的，由计算出的 l l 最小值给定弹性临界载荷最小值给定弹性临界载荷Pcr。September 30, 2001Inventory #0014919-19几何不稳定性几何不稳定性 . 特征值屈曲过程特征值屈曲过程特征值屈曲分析包括以下四个主要步骤特征值屈曲分析包括以下四个主要步骤: 1. 建模建模2. 获得带有预应力的静力解获得带有预应力的静力解3. 获得特征值屈曲解获得特征值屈曲解4. 查看结果查看结果September 30, 2001Inventory #0014919-20几何不稳定性几何不稳定性 . 特征值屈曲过程特征

15、值屈曲过程建模建模该任务与大多数其它分析类似该任务与大多数其它分析类似, 除了下面两点除了下面两点: 只有线性行为有效，非线性单元处理为线性，它们的刚度基于初始状只有线性行为有效，非线性单元处理为线性，它们的刚度基于初始状态态, 且不能改变。且不能改变。 必须定义杨氏模量，材料特性可能是线性、各向同性或各向异性必须定义杨氏模量，材料特性可能是线性、各向同性或各向异性, 忽忽略非线性特性。略非线性特性。September 30, 2001Inventory #0014919-21几何不稳定性几何不稳定性 . 特征值屈曲过程特征值屈曲过程获得带有预应力的静力解获得带有预应力的静力解当获得静力解时当

16、获得静力解时, 必须必须 设置预应力标识设置预应力标识, 以进行后面的特征值屈曲以进行后面的特征值屈曲分析。分析。 Main Menu Preprocessor Loads Analysis Options 或键入命令或键入命令: PSTRES,ONSeptember 30, 2001Inventory #0014919-22几何不稳定性几何不稳定性 . 特征值屈曲过程特征值屈曲过程获得带有预应力的静力解获得带有预应力的静力解通常单位载荷就足够了，计算出的特征值代表施加载荷上的屈曲载通常单位载荷就足够了，计算出的特征值代表施加载荷上的屈曲载荷因子。荷因子。注意特征值代表所有载荷的比例因子，若某

17、载荷是常数，而其它载注意特征值代表所有载荷的比例因子，若某载荷是常数，而其它载荷是变量荷是变量, 则需确保常载荷的应力刚度矩阵没有被乘以因子则需确保常载荷的应力刚度矩阵没有被乘以因子(后面讨后面讨论论)求解模型求解模型 Main Menu Solution -Solve- Current LS 或键入命令或键入命令: SOLVESeptember 30, 2001Inventory #0014919-23几何不稳定性几何不稳定性 . 特征值屈曲过程特征值屈曲过程获得特征值屈曲解获得特征值屈曲解完成静态求解后完成静态求解后, 退出并重新进入求解器退出并重新进入求解器, 并指定分析类型为特征并指定

18、分析类型为特征值屈曲值屈曲: Main Menu Finish Main Menu Solution -Analysis Type- New Analysis 或键入命令或键入命令: FINISH/SOLUANTYPE,BUCKLESeptember 30, 2001Inventory #0014919-24几何不稳定性几何不稳定性 . 特征值屈曲过程特征值屈曲过程获得特征值屈曲解获得特征值屈曲解指定特征值提取方法和要提取的屈曲模态数目指定特征值提取方法和要提取的屈曲模态数目: Main Menu Solution Analysis Options 或键入命令或键入命令: BUCOPT,LAN

19、B,3,0Block Lanczos 是推荐的特是推荐的特征值提取方法。本例中征值提取方法。本例中, 要要求求 3 个模态。个模态。September 30, 2001Inventory #0014919-25几何不稳定性几何不稳定性 . 特征值屈曲过程特征值屈曲过程获得特征值屈曲解获得特征值屈曲解指定要写入结果文件模态数。指定要写入结果文件模态数。 Main Menu Solution -Load Step Opts - Expansion Pass Expand Modes . 或键入命令或键入命令: BUCOPT,LANB,3,0也可计算出相应的应力也可计算出相应的应力分布。分布。Sep

20、tember 30, 2001Inventory #0014919-26几何不稳定性几何不稳定性 . 特征值屈曲过程特征值屈曲过程常量和变量载荷的注释常量和变量载荷的注释可以对特征值进行迭代可以对特征值进行迭代, 调整变量载荷直到特征值变为调整变量载荷直到特征值变为1.0 或接近或接近于于 1.0。考虑一个自重为。考虑一个自重为 WO 和外加载荷和外加载荷 A 的杆的例子，可以迭代的杆的例子，可以迭代，调整，调整 A 的值直到的值直到 l l = 1.0。September 30, 2001Inventory #0014919-27几何不稳定性几何不稳定性 . 特征值屈曲过程特征值屈曲过程查看

21、结果查看结果可以在通用后处理器中查看特征值屈曲分析的结果，结果包括载荷可以在通用后处理器中查看特征值屈曲分析的结果，结果包括载荷因子、屈曲模态和相对应力分布。因子、屈曲模态和相对应力分布。 Main Menu General Postproc Results Summary . 或键入命令或键入命令: SET,LIST“Set” 列表明屈曲模态列表明屈曲模态数数, “Time” 值表示相应值表示相应的载荷因子。的载荷因子。September 30, 2001Inventory #0014919-28几何不稳定性几何不稳定性 . 特征值屈曲过程特征值屈曲过程查看结果查看结果屈曲模态的最大位移归一

22、化为屈曲模态的最大位移归一化为1.0，因此，因此, 位移不能代表真实的变形位移不能代表真实的变形，且应力是相对于屈曲模态。，且应力是相对于屈曲模态。通常查看最初少数的屈曲模态是有益的，在随后的非线性屈曲分析通常查看最初少数的屈曲模态是有益的，在随后的非线性屈曲分析中中, 结构的高阶屈曲模态可能是重要的。结构的高阶屈曲模态可能是重要的。若存在密排的特征值若存在密排的特征值, 这表明该结构对缺陷敏感，应执行具有适当这表明该结构对缺陷敏感，应执行具有适当的缺陷或扰动的非线性屈曲分析。的缺陷或扰动的非线性屈曲分析。September 30, 2001Inventory #0014919-29几何不稳定

23、性几何不稳定性 . 特征值屈曲过程特征值屈曲过程其它考虑事项其它考虑事项:有些情况下有些情况下, 在特征值屈曲分析中计算出负的特征值，在特征值提在特征值屈曲分析中计算出负的特征值，在特征值提取过程中遇到数值困难时会发生这种情况。取过程中遇到数值困难时会发生这种情况。 在这种情况下，可指在这种情况下，可指定特征值提取的偏移点定特征值提取的偏移点 (BUCOPT)， 在偏移点附近提取特征值最在偏移点附近提取特征值最精确，这需要对临界载荷值有一定的了解。精确，这需要对临界载荷值有一定的了解。September 30, 2001Inventory #0014919-30几何不稳定性几何不稳定性 . 特

24、征值屈曲过程特征值屈曲过程其它考虑事项其它考虑事项:在屈曲分析中在屈曲分析中, 压力压力-载荷刚度矩阵对精确地计算载荷因子通常是重载荷刚度矩阵对精确地计算载荷因子通常是重要的。缺省时要的。缺省时, 对特征值屈曲分析对特征值屈曲分析 ANSYS自动自动 包括压力包括压力-载荷刚度载荷刚度矩阵。尽管不是推荐的矩阵。尽管不是推荐的, 用户仍可手动激活或停用压力用户仍可手动激活或停用压力-载荷刚度的载荷刚度的使用使用, 通过通过: Main Menu Solution Unabridged Menu Main Menu Solution -Load Step Opts- Solution Ctrl 或

25、键入命令或键入命令: SOLCONTROL,INCPSeptember 30, 2001Inventory #0014919-31几何不稳定性几何不稳定性C. 非线性屈曲背景非线性屈曲背景下图为一般的非线性载荷变形曲线，该图说明理想载荷路径、有缺下图为一般的非线性载荷变形曲线，该图说明理想载荷路径、有缺陷结构的载荷路径和该结构的实际动态响应。陷结构的载荷路径和该结构的实际动态响应。Fu理想载荷路径理想载荷路径有缺陷结构的载荷路径有缺陷结构的载荷路径 实际动态响应实际动态响应前屈曲前屈曲后屈曲后屈曲理想静态行理想静态行为为 分叉点分叉点极限点极限点September 30, 2001Invent

26、ory #0014919-32几何不稳定性几何不稳定性 . 非线性屈曲背景非线性屈曲背景前面前面 B 节中讨论了节中讨论了线性线性 特征值屈曲过程特征值屈曲过程(前面幻灯片中的理想载荷前面幻灯片中的理想载荷路径路径)。有几种分析技术用于计算结构的有几种分析技术用于计算结构的非线性非线性 静力变形响应，这些技术静力变形响应，这些技术包括包括: 载荷控制载荷控制 位移控制位移控制 弧长法弧长法September 30, 2001Inventory #0014919-33几何不稳定性几何不稳定性 . 非线性屈曲背景非线性屈曲背景载荷控制载荷控制:如下图所示如下图所示, 考虑浅拱的快速通过分析，当以增

27、量载荷考虑浅拱的快速通过分析，当以增量载荷 (F) 求解该求解该问题时问题时, 求解采用求解采用载荷控制载荷控制来完成。来完成。FFFFappu用载荷控制能达到用载荷控制能达到Fapp吗？吗？September 30, 2001Inventory #0014919-34几何不稳定性几何不稳定性 . 非线性屈曲背景非线性屈曲背景载荷控制载荷控制:使用使用 Newton-Raphson 载荷控制的困难是求解不能通过不稳定点载荷控制的困难是求解不能通过不稳定点。在不稳定点。在不稳定点 (Fcr), 切线刚度矩阵切线刚度矩阵 KT 是奇异的，使用载荷控制是奇异的，使用载荷控制, Newton-Raph

28、son 法不收敛。然而法不收敛。然而, 该类型的分析对描述结构的该类型的分析对描述结构的前屈曲前屈曲 行为是有用的。行为是有用的。FappuFcrKT = 0使用载荷控制只有使用载荷控制只有Fcr 可达到。可达到。KT 0September 30, 2001Inventory #0014919-35几何不稳定性几何不稳定性 . 非线性屈曲背景非线性屈曲背景位移控制位移控制:当拱由增量位移加载时当拱由增量位移加载时, 与力相反与力相反, 采用采用位移控制位移控制 进行求解。位移进行求解。位移控制的优点是控制的优点是, 除除 Fcr外外, 它产生一个稳定的解。它产生一个稳定的解。(强加的位移在不强

29、加的位移在不稳定点提供一个附加约束。稳定点提供一个附加约束。)FappuUYUYUY用位移控制用位移控制 能够达到能够达到 Fapp. ( 此时此时 Fapp是强是强加的位移加的位移 UY处的反作处的反作用力。用力。)September 30, 2001Inventory #0014919-36几何不稳定性几何不稳定性 . 非线性屈曲背景非线性屈曲背景位移控制位移控制:位移控制的缺点是只有在知道施加什么位移时才适用位移控制的缺点是只有在知道施加什么位移时才适用! 如果拱上施如果拱上施加压力载荷加压力载荷, 而不是集中力而不是集中力, 位移控制不可能使用。位移控制不可能使用。P对于较复杂的载荷状

30、态对于较复杂的载荷状态, 一般一般也不清楚施加什么位移。也不清楚施加什么位移。September 30, 2001Inventory #0014919-37几何不稳定性几何不稳定性 . 非线性屈曲背景非线性屈曲背景弧长法弧长法:弧长法弧长法 是一种求解方法是一种求解方法, 用于获得不稳定性问题用于获得不稳定性问题 (KT 0) 或负的或负的切线刚度切线刚度 (KT 0) 的数值稳定解。的数值稳定解。弧长法可用于弧长法可用于比例载荷比例载荷 的的静态静态 问题。问题。尽管弧长法能求解复杂的力尽管弧长法能求解复杂的力-变形响应问题变形响应问题, 但它最适合求解没有突但它最适合求解没有突然分叉点的平

31、滑响应问题。然分叉点的平滑响应问题。FuSeptember 30, 2001Inventory #0014919-38几何不稳定性几何不稳定性 . 非线性屈曲背景非线性屈曲背景弧长法弧长法:弧长法同时求解载荷和位移弧长法同时求解载荷和位移, 与与 Newton-Raphson 法相似，然而法相似，然而,引入了一个附加的未知项引入了一个附加的未知项-载荷因子载荷因子l l (-1 l l Preprocessor -Modeling- Update Geom 或键入命令或键入命令: UPGEOM施加到原始几何形状上施加到原始几何形状上的位移的乘子的位移的乘子模态数模态数特征值屈曲分析的结果特征值

32、屈曲分析的结果文件。文件。September 30, 2001Inventory #0014919-48几何不稳定性几何不稳定性 . 非线性前屈曲过程非线性前屈曲过程建模建模 - 初始缺陷初始缺陷初始缺陷的量级将影响非线性屈曲分析的结果，初始缺陷将去除载初始缺陷的量级将影响非线性屈曲分析的结果，初始缺陷将去除载荷荷- 变形响应中的明显不连续性。变形响应中的明显不连续性。缺陷的值相对于结构的总体尺寸是小的，该值应该与实际结构中的缺陷的值相对于结构的总体尺寸是小的，该值应该与实际结构中的缺陷缺陷(真实的或假设的真实的或假设的)的尺寸匹配，制造公差可用于估计不完整性的尺寸匹配，制造公差可用于估计不完

33、整性的量级。的量级。September 30, 2001Inventory #0014919-49几何不稳定性几何不稳定性 . 非线性前屈曲过程非线性前屈曲过程求解求解非线性屈曲分析是具有几何非线性效应的静态分析非线性屈曲分析是具有几何非线性效应的静态分析, 分析延伸至结分析延伸至结构的极限载荷点，确保激活几何非线性构的极限载荷点，确保激活几何非线性(NLGEOM,ON)。推荐使用求解控制推荐使用求解控制(缺省缺省)。使用全使用全 Newton-Raphson 选项选项, 不打开自适应下降不打开自适应下降(求解控制的缺求解控制的缺省项省项)。September 30, 2001Inventor

34、y #0014919-50几何不稳定性几何不稳定性 . 非线性前屈曲过程非线性前屈曲过程求解求解激活自动时间步长激活自动时间步长(求解控制的缺省项求解控制的缺省项)，打开自动时间步长时，打开自动时间步长时, 程程序序自动自动 搜索屈曲载荷，若在给定载荷下求解不收敛搜索屈曲载荷，若在给定载荷下求解不收敛, 则程序二分并则程序二分并在一个更小的载荷下尝试新的求解，同样地在一个更小的载荷下尝试新的求解，同样地, 最小时间步长最小时间步长 将影响将影响结果的精度。结果的精度。uFapp“Time”1,3,524子步数子步数FlimitFapp67September 30, 2001Inventory

35、#0014919-51几何不稳定性几何不稳定性 . 非线性前屈曲过程非线性前屈曲过程求解求解若时间步较大若时间步较大, 则可能（虽然未必）则可能（虽然未必）“跳过跳过”不稳定点获得一个不稳定点获得一个“快速通过快速通过”的求解，务必在时间历程后处理器中画出载荷位移曲线的求解，务必在时间历程后处理器中画出载荷位移曲线。FappuFFF稳定的稳定的不稳定的不稳定的稳定的稳定的September 30, 2001Inventory #0014919-52几何不稳定性几何不稳定性 . 非线性前屈曲过程非线性前屈曲过程求解求解务必设置小的务必设置小的最小时间步长最小时间步长 以允许二分。以允许二分。对于

36、外加载荷值对于外加载荷值, 采用比特征值屈曲载荷高采用比特征值屈曲载荷高 1020% 的值通常是好的值通常是好的选择。前已述及的选择。前已述及, 为便于后处理为便于后处理, 可设定可设定“时间时间”等于外加载荷值等于外加载荷值。务必写出足够的子步数的结果务必写出足够的子步数的结果(OUTRES), 以便于在通用后处理器中以便于在通用后处理器中能查看载荷位移曲线。能查看载荷位移曲线。September 30, 2001Inventory #0014919-53几何不稳定性几何不稳定性 . 非线性前屈曲过程非线性前屈曲过程查看结果查看结果在时间历程后处理器中查看载荷位移曲线，在面内载荷的情况下在时

37、间历程后处理器中查看载荷位移曲线，在面内载荷的情况下, 需要画出面外需要画出面外(侧向的侧向的)位移与载荷。位移与载荷。FFUXxy画出力画出力 (F) 和和 X方向方向(UX)的位移。的位移。September 30, 2001Inventory #0014919-54几何不稳定性几何不稳定性 . 非线性前屈曲过程非线性前屈曲过程其它考虑事项其它考虑事项:载荷位移曲线能帮助确定求解发散的原因是物理不稳定性还是数值载荷位移曲线能帮助确定求解发散的原因是物理不稳定性还是数值不稳定性不稳定性。要。要认识到认识到不收敛的解未必意味着结构达到其最大载荷不收敛的解未必意味着结构达到其最大载荷!Septe

38、mber 30, 2001Inventory #0014919-55几何不稳定性几何不稳定性 . 非线性前屈曲过程非线性前屈曲过程其它考虑事项其它考虑事项:当结构在其屈曲载荷附近时当结构在其屈曲载荷附近时, 切线刚度将接近于切线刚度将接近于 0，从载荷，从载荷-位移曲位移曲线的斜率可以确定是数值或物理不稳定。线的斜率可以确定是数值或物理不稳定。数值不稳定性数值不稳定性物理不稳定性物理不稳定性(屈曲屈曲)uFapp不收敛解不收敛解最后收敛解最后收敛解KT 0uFapp不收敛解不收敛解最后收敛解最后收敛解KT 0September 30, 2001Inventory #0014919-56几何不稳

39、定性几何不稳定性 . 练习练习请参考附加练习题请参考附加练习题:练习练习12: 拱的屈曲拱的屈曲September 30, 2001Inventory #0014919-57几何不稳定性几何不稳定性E. 非线性后屈曲过程非线性后屈曲过程后屈曲分析技术包括后屈曲分析技术包括: 位移控制位移控制 动力学动力学 弧长法弧长法Fu有缺陷结构的静态载荷有缺陷结构的静态载荷路径路径实际动态响应实际动态响应后屈曲后屈曲理想静态行理想静态行为为极限点极限点September 30, 2001Inventory #0014919-58几何不稳定性几何不稳定性 . 非线性后屈曲过程非线性后屈曲过程位移控制位移控制

40、:如前所述如前所述, 位移控制可用于预测简单载荷条件下的后屈曲响应，位位移控制可用于预测简单载荷条件下的后屈曲响应，位移控制的主要缺点是对于复杂载荷移控制的主要缺点是对于复杂载荷, 一般不清楚应施加什么位移。一般不清楚应施加什么位移。FappuUYUYUY用位移控制可达到用位移控制可达到 Fapp. September 30, 2001Inventory #0014919-59几何不稳定性几何不稳定性 . 非线性后屈曲过程非线性后屈曲过程动力学动力学:使用载荷控制的静态稳定性问题能够通过非线性瞬态动力学分析来使用载荷控制的静态稳定性问题能够通过非线性瞬态动力学分析来求解。动力学分析中求解。动力

41、学分析中, 当结构当结构“快速通过快速通过”(动态跃变动态跃变)时时, 不计算不计算软化响应，动力学的主要缺点是它不易衰减不希望的动态效应软化响应，动力学的主要缺点是它不易衰减不希望的动态效应(振振铃铃)。FappuFFF动态响应动态响应September 30, 2001Inventory #0014919-60几何不稳定性几何不稳定性 . 非线性后屈曲过程非线性后屈曲过程弧长法弧长法:弧长法弧长法 是一种求解方法，用于获得在是一种求解方法，用于获得在比例比例 载荷下经历不稳定性载荷下经历不稳定性 (KT 0) 或负的切线刚度或负的切线刚度 (KT Solution -Analysis Ty

42、pe- Soln Control Solution Control Advanced NL tab 或键入命令或键入命令: ARCLENSeptember 30, 2001Inventory #0014919-67几何不稳定性几何不稳定性 . 非线性后屈曲过程非线性后屈曲过程求解求解 - 弧长参数弧长参数因子因子MAXARC 和和MINARC 是弧长半径乘子是弧长半径乘子, 用于定义弧长半径的用于定义弧长半径的极限。极限。(类似于自动时间步长的子步数的设定。类似于自动时间步长的子步数的设定。)MAXARC 缺省值为缺省值为 10MINARC 缺省值为缺省值为 0.001September 30

43、, 2001Inventory #0014919-68几何不稳定性几何不稳定性 . 非线性后屈曲过程非线性后屈曲过程求解求解 - 弧长终止准则弧长终止准则弧长求解通过载荷因子弧长求解通过载荷因子(l l) 达到达到 1.0 或设定求解极限来终止。或设定求解极限来终止。弧长求解可以通过达到第一个极限点或达到指定节点的最大位移准弧长求解可以通过达到第一个极限点或达到指定节点的最大位移准则来终止。则来终止。September 30, 2001Inventory #0014919-69几何不稳定性几何不稳定性 . 非线性后屈曲过程非线性后屈曲过程查看结果查看结果弧长分析中时间与载荷因子有关，因此当后处理弧长分析时弧长分析中时间与载荷因子有关，因此当后处理弧长分析时, 不要不要 参照参照 “time” 值结果