清华大学计算固体力学第六次课件_求解方法与稳定性.ppt

5 稳定性和连续方法 平衡解答的分支 为了更好地理解一个系统，必须确定它的平衡路径或者分支，以及它们的稳定性。在结构力学中，通过直接获得动态解答可以回避与不稳定行为相关的困难，如关于位移摄动的系数行列式为零的非平凡解答。当结构加载超过它的临界点或在动态模拟中的分支点时，结构动态地通过了最接近的稳定分支。然而，不稳定是不容易出现的，并且缺陷敏感性的可能性也是不清晰的。因此，为了理解结构的行为或者整个过程，必须小心地检查它的平衡行为。 5 稳定性和连续方法 平衡解答的分支 一个具有反对称分支的系统对于缺陷是非常敏感的，在图(e)中，最大荷载随着缺陷变化。在实际结构中缺陷是不可避免的，一个完好结构的理论分岔点不是强度的真正度量；实际结构可以在一个远远低于理论值的荷载下屈曲。 许多结构行为的奇怪表现可能被动态模拟掩盖了。例如在承受轴向荷载的圆柱壳中，交叉分支是反对称的，如图(e)所示。 5 稳定性和连续方法 平衡解答的分支 为了达到确定非线性系统性质的目的，通常将平衡解答分成若干分支，当一个参数变化时，它是连续地线性描述系统的响应。称这些分支为平衡分支。追踪模型的平衡分支作为参数 的函数。 寻找 使得残数为零（即系统是处于平衡的） 非线性系统显示了三种分岔行为： 1 转折点：结构分析中的临界点，在该点上分支斜率改变符号。 2 静止分岔：简单分岔，在这些点上两个平衡分支交叉。 3 Hopf分岔：在这些点上平衡分支与一个周期运动的分支交叉。 5 稳定性和连续方法 平衡解答的分支 图示浅桁架行为展示了一个转折点（或者临界点）；A和B是转折点。在转折点之后，一个分支可以是稳定或者是不稳定的。在这种情况下，在例6.4会看到，第一个转折点A 后面的分支是不稳定的，而第二个转折点B 后面的分支是稳定的。(例6.4中字母有区别) 5 稳定性和连续方法 平衡解答的分支 在图a中显示的梁是分岔的典型例题。两个分支的交叉点B是分岔点。在分岔点后，基本分支的继续部分BC成为不稳定的。分岔点B对应于Euler梁的屈曲荷载。这种类型的分支通常称为音叉。 Hopf分岔在率无关结构中不常见，发生在率相关材料或者主动控制的问题。在Hopf分岔中，在一个分支的终点不可能得到稳定平衡的结果。代之是这个平衡分支分成为两个分支，其结果对时间是周期性变化的。 5 稳定性和连续方法 连续方法和弧长法 跟踪平衡分支的方法称为连续方法。平衡分支的跟踪是相当困难的；对于连续性，仍然没有实现强健和自动化的过程。下面描述基于参数化的连续方法-弧长法。 如果采用载荷增量求解，希望在A附近收敛，有可能跳到另一个替代平衡点A’；同样情况也适用于位移增量求解，有可能从B跳到B’。如何跟踪真正的平衡轨迹，采用弧长法。 5 稳定性和连续方法 连续方法和弧长法 在跟踪分支时，荷载参数通常从零开始并逐渐增加。对于参数 ? 的每一个增量，计算平衡方程解答，即找到关于方程的解。 或者 在弧长法中，除了逐渐增加荷载参数 ? 之外，在位移—荷载参数空间中的弧长度量也需要逐渐增加，这是通过在平衡方程中增加参数化方程实现的。 弧长增量 载荷增量 首先讨论如何应用荷载参数跟踪分支。 位置变化 5 稳定性和连续方法 连续方法和弧长法 在步骤中弧长被横切的近似值 载荷增量比例因数 一个可能的比例因数与材料刚度矩阵的对角线单元有关，即 产生许多其它类型的参数化方程。方程的整个系统包括平衡方程和 弧长方程 参数化方程 (1) (2) 5 稳定性和连续方法 连续方法和弧长法 这样，有了一个附加的方程和一个附加的未知数 使得弧长s(圆的半径)是逐渐增加的。弧长法既不是力加载，也不是位移加载，而是载荷－位移混合加载。 假设节点外部荷载的分布不随模型的变形而改变；当荷载因子是从属荷载的系数时，如压力，根据荷载几何的影响，必须修改方法以考虑到荷载分布的变化。上述不包括惯性项，因为连续方法仅适用于平衡问题。 在结构力学中，弧长法经常称为Riks方法（Riks,1972） 对于一个自由度的问题，这一过程是最容易解释的。图示浅桁架，假设在点n得到平衡解答。在荷载—位移平面中，弧长方程(1)是一个围绕点n的圆。参数方程(2)的解答是与该圆交叉的平衡分支。在点n增量荷载参数是无效的，因为它又带回到分支上。 在弧长法中，以沿分支的弧长形式重新表示问题，因此一个分支可能跟随一个下降的荷载。在步骤中不需要增加荷载，而仅需要增加弧长参数，这是跟踪分支的完美的自然方式。 在两个自由度的问题中，弧长方程可以围绕平衡点定义一个球面或球，以及分支与球交叉的一个