《非线性分析中弧长法的简介》由会员分享，可在线阅读，更多相关《非线性分析中弧长法的简介（7页珍藏版）》请在装配图网上搜索。
1、非线性分析弧长法的读书报告1.弧长法的原理弧长法属于双重目标控制方法，即在求解过程中同时控制荷载因子和位移增量的步长。其基本的控制方程为： （1）式中：为荷载因子增量数值为荷载比例系数，用于控制弧长法中荷载因子增量所占的比重为固定的半径在求解过程中，荷载因子增量在迭代中是变化的，下列非线性静力平衡的迭代求解公式中存在个未知数，即 （2）这样，在弧长法中一共存在个未知数，根据约束方程： 即为附加的控制方程，问题才能得到解答，此时，可以根据值的取值分为两种弧长法，其中，时的弧长法称为球面弧长法，时的弧长法称为柱面弧长法。1.1球面弧长法如下图所示，根据图中所示来说明弧长法的求解策略：在第步荷载增量
2、第次迭代分析中，结构的位移增量可由式（2）来计算，即 （3）由于刚度矩阵不对称以及带宽被改变，通常情况下直接联立求解式（1）和式（3）中个变量相当困难。通常，将式(3)中分解为两个部分，即 （4）式中， （5） 说明的第一项为采用荷载控制的标准切线刚度迭代求解的结果； （6） 说明的第二项为考虑荷载下按当前结构刚度计算出来的位移增量。至此，式(4)中仍有为未知量，现在只有借助控制方程（1）来求解。第次迭代结束后相应于迭代初始点（上图中点）的位移增量为 (7)将（7）式代入控制方程（1）式，并考虑迭代前后弧长保持不变，故得到： (8)式中：与分别为第次迭代前后荷载因子相对于迭代初始点（上图中 点
3、）的增量，它们有如下关系： (9)将(4)代入(9)中，可得： （10）式中：由此，可通过解上述一元二次方程得到何在因子的增量，从而进一步确定当前的荷载水平和位移向量，即(11)(12)1.2柱面弧长法经过研究发现，荷载比例系数对于最终分析结果的影响是有限的，尤其是结构非线性成都较高时，这种影响甚微。令，从而减少程序中的未知数和提高求解效率。于是把从原先所谓的“球面弧长法”简化为“柱面弧长法”，结构的控制方程则简化成： （13）同理可得：(10)中关于的一元二次方程式中的系数可以简化成:由上式求解出两个根，从而得到两个位移增量，即： （14） （15）为了保证迭代方向尽量保持一致，可以根据两次
4、迭代结束时的位移增量和的夹角最小准则来判断到底哪个更合适。两个向量之间的夹角可以按照下式确定： （16）将一元二次方程求出的两个解分别代入上式求出各自的夹角，夹角小的的荷载因子增量为所求的值。可按照式（11）和（12）更新结构的荷载水平和位移向量。1.3弧长法的简化形式上述弧长法的求解过程，需要求解一元二次方程，计算量大，因此，为简化计算，提出了另一种控制方程，用垂直于迭代向量的平面代替圆弧，把弧长不变的条件改为向量与向量始终保持正交，即满足下列控制方程： （17）写成矩阵形式为与前面的解法相同，可求解上述一元二次方程得相比之下，用上式求解容易多了，其余步骤完全同上述相关内容。1.4弧长法的求
5、解步骤（1） 对于第1个增量步()第1次迭代()分析,选定参考荷载，即 确定了初始弧长增量。（2） 输入期望迭代次数；如果采用球面弧长法，则输入荷载参与比例系 数。（3） 存储结构初始切线刚度。（4） 在第次增量步分析中，迭代流程如下。 求解出。 记录迭代次数，对结构刚度矩阵进行三角分解或计算“当前刚度系数”以判别矩阵是否正定。 更新结构的变形向量，计算结构的恢复力箱梁和非平衡力向量。 如果采用切线刚度迭代技术，则要根据当前结构的变形向量更新结构刚度矩阵；如果采用初始刚度迭代技术，则只需在每次增量分析中的初始迭代中根据上一次增量结束时的结构位移向量来更新结构刚度即可，在增量步中不用更新。 计算
6、和，若采用初始刚度迭代技术，在整个增量步迭代中为定值，不必重复计算。 求解荷载因子增量。 按式（4）计算，由式（11）、（12）更新当前的荷载水平和位移向量。 收敛性判别。如果满足收敛准则，则终止当前增量步下的迭代进城，记录迭代次数进入步骤；如果不满足收敛准则，则需要继续迭代，记录迭代次数，令，重复步骤 判别当前荷载水平是否达到期望值或超过一定的增量步数。如是，则分析结束，输出数据；如不是，则令，更新结构切线刚度矩阵，计算当前增量步中的弧长增量，返回步骤。2.弧长法的适用范围弧长法适用性很强，收敛性和稳定性明显好于其他处理负刚度问题的方法，它既可以用于加工软化结构，也可以适用于加工硬化结构，在
7、非线性程度较高的体系应优先考虑采用该方法。但是该方法的计算量很大，对一般非线性问题，还是建议选择其他简单的方法。3.弧长法的评述（1）弧长法属于双重目标控制方法，即在求解过程中同时控制荷载因子和位移增量的步长，从理论上来说，任何方法都应在极值点附近存在刚度奇异的问題，但是控制位移法和弧长法中，迭代点正好落在极值点附近的概率很小，在现实中很难遇见，除非遇到非线性程度很高的结构体系。在分析混凝土开裂的加工软化问題时，偶尔也会在运用弧长法时出现一些问題，如在应用球面孤长法时得不到方程式的实数解，或应用柱面弧长法时在临界点附近发散，找不到交点。（2）弧长法属于自动步长法，只要给出一些控制参数，步长由程序自动计算，此方法为当前的主流计算方法。（3）弧长法中的荷载增量预测因子只是初步控制当前荷载增量步的迭代过程，不能完全准确估计结构迭代结束时的荷载水平，结构最终的荷载水平是迭代结束后的数值。4.参考文献1何政，欧进萍.钢筋混凝土结构非线性分析M.哈尔滨工业大学出版社,2007.2向天宇，赵人达，刘海波. 将弧长法应用于结构的几何非线性有限元分析J.铁道学报,2003.6,(2).3向天宇，赵人，刘海波.混凝土结构全过程非线性分析的弧长法研究J.铁道学报,2002,(3).4江见鲸.钢筋混凝土结构非线性有限元分析M.陕西科学技术出版社,1994.
非线性分析中弧长法的简介
