非线性分析中弧长法的简介.docVIP

非线性分析弧长法的读书报告 1.弧长法的原理 弧长法属于双重目标控制方法，即在求解过程中同时控制荷载因子和位移增量的步长。其基本的控制方程为： （1） 式中： 为荷载因子增量数值 为荷载比例系数，用于控制弧长法中荷载因子增量所占的比重 为固定的半径 在求解过程中，荷载因子增量在迭代中是变化的，下列非线性静力平衡的迭代求解公式中存在个未知数，即 （2） 这样，在弧长法中一共存在个未知数，根据约束方程： 即为附加的控制方程，问题才能得到解答，此时，可以根据值的取值分为两种弧长法，其中，时的弧长法称为球面弧长法，时的弧长法称为柱面弧长法。 1.1球面弧长法 如下图所示，根据图中所示来说明弧长法的求解策略： 在第步荷载增量第次迭代分析中，结构的位移增量可由式（2）来计算，即 （3） 由于刚度矩阵不对称以及带宽被改变，通常情况下直接联立求解式（1）和式（3）中个变量相当困难。通常，将式(3)中分解为两个部分，即 （4） 式中， （5） 说明的第一项为采用荷载控制的标准切线刚度迭代求解的 结果； （6） 说明的第二项为考虑荷载下按当前结构刚度计算出来的位移增量。 至此，式(4)中仍有为未知量，现在只有借助控制方程（1）来求解。第次迭代结束后相应于迭代初始点（上图中点）的位移增量为 (7) 将（7）式代入控制方程（1）式，并考虑迭代前后弧长保持不变，故得到： (8) 式中： 与分别为第次迭代前后荷载因子相对于迭代初始点（上图中 点）的增量，它们有如下关系： (9) 将(4)代入(9)中，可得： （10） 式中： 由此，可通过解上述一元二次方程得到何在因子的增量，从而进一步确定当前的荷载水平和位移向量，即 (11) (12) 1.2柱面弧长法 经过研究发现，荷载比例系数对于最终分析结果的影响是有限的，尤其是结构非线性成都较高时，这种影响甚微。令，从而减少程序中的未知数和提高求解效率。于是把从原先所谓的“球面弧长法”简化为“柱面弧长法”，结构的控制方程则简化成： （13） 同理可得：(10)中关于的一元二次方程式中的系数可以简化成: 由上式求解出两个根，从而得到两个位移增量，即： （14） （15） 为了保证迭代方向尽量保持一致，可以根据两次迭代结束时的位移增量和的夹角最小准则来判断到底哪个更合适。 两个向量之间的夹角可以按照下式确定： （16） 将一元二次方程求出的两个解分别代入上式求出各自的夹角，夹角小的的荷载因子增量为所求的值。可按照式（11）和（12）更新结构的荷载水平和位移向量。 1.3弧长法的简化形式 上述弧长法的求解过程，需要求解一元二次方程，计算量大，因此，为简化计算，提出了另一种控制方程，用垂直于迭代向量的平面代替圆弧，把弧长不变的条件改为向量与向量始终保持正交，即满足下列控制方程： （17） 写成矩阵形式为 与前面的解法相同，可求解上述一元二次方程得 相比之下，用上式求解容易多了，其余步骤完全同上述相关内容。 1.4弧长法的求解步骤 对于第1个增量步()第1次迭代()分析,选定参考荷载，即 确定了初始弧长增量。 输入期望迭代次数；如果采用球面弧长法，则输入荷载参与比例系 数。 存储结构初始切线刚度。 在第次增量步分析中，迭代流程如下。 = 1 \* GB3 ① 求解出。 = 2 \* GB3 ② 记录迭代次数，对结构刚度矩阵进行三角分解或计算“当前刚度系数”以判别矩阵是否正定。 = 3 \* GB3 ③ 更新结构的变形向量，计算结构的恢复力箱梁和非平衡力向量。 = 4 \* GB3 ④ 如果采用切线刚度迭代技术，则要根据当前结构的变形向量更新结构刚度矩阵；如果采用初始刚度迭代技术，则只需在每次增量分析中的初始迭代中根据上一次增量结束时的结构位移向量来更新结构刚度即可，在增量步中不用更新。 = 5 \* GB3 ⑤ 计算和，若采用初始刚度迭代技术，在整个增量步迭代中为定值，不必重复计算。 = 6 \* GB3 ⑥ 求解荷载因子增量。 = 7 \* GB3 ⑦ 按式（4）计算，由式（1