第二单元：非线性代数方程组的数值解法.pptVIP

2000.3 哈尔滨建筑大学 王焕定教授制作 1.1 直接迭代法 1.2 牛顿法和修正牛顿法 1.3 增量方法 1.4 增量弧长法 * 非线性代数方程组的数值解法 有限单元法II ——2004级硕士生课程 主讲教师：周志勇 (博士) 同济大学土木工程桥梁工程系 同济大学土木学院桥梁工程系 非线性问题可分为三类：材料非线性、几何非线性和边界非线性。我们只讨论前两类问题。 不管那类非线性问题，最终都归结为一组非线性方程Ψ(a)=0，a为待求的未知量。 Ψ(a)=0可写成平衡方程的形式 Ψ(a) =P(a)-R=K(a) a -R=0 对非线性方程Ψ(a)=0，一般只能用数值方法求近似解答。其实质是，用一系列线性方程组的解去逼近所讨论非线性方程组的解。 分段线性法 1.1 直接迭代法 1.2 牛顿法和修正牛顿法 1.3 增量方法 1.4 增量弧长法 Ψ(a) =P(a)-R=K(a) a -R=0 设初始未知量为a0，根据上式有 a1= K(a0)-1R 如果问题是收敛的， a1将比a0有所改善。如此反复迭代可得 an+1= K(an)-1R Δan=an+1- an 当设范数为 或设范数为 收敛条件则为 将Ψ(an)视为不平衡力并作为衡量收敛的标准 对于单变量问题的非线性方程，直接迭代法的计算过程如图1所示。图上给出的是 和 之间的关系，而不是 和 ， 之间的关系 对单变量情况，直接迭代实质是“割线”法，一定条件下这种迭代过程是收敛的，但对多自由度情况，由于未知量通过矩阵K(an)的元素互相耦合，在迭代过程中可能会出现不稳定现象。 a1= K(a0)-1R 直至Δan=an+1- an 满足收敛条件 收敛性 如果将非线性方程Ψ(a) =0在an 附近展开，则 记 KT(an)=[Ψ’(a)]n，Pn =Ψ(an) Δan≈-[Ψ’(a)]n-1Ψ(an) 切线矩阵 不平衡力 如此逐步计算，即可得到非线性方程的解答，这就是牛顿-拉夫森法。 Ψ(a) =Ψ(an)+ [Ψ’(a)]n Δan+。。。 =0 又如果[Ψ’(a)]n的逆存在，则Δan 近似等于 则 Δan≈-KT(an)-1 Pn ， an+1=an+Δan Δan≈-KT(an)-1 Pn ， an+1=an+Δan 直至 Δan 满足收敛性 1.2 牛顿法和修正牛顿法 如果在迭代计算的每一步内，矩阵KT都用初始近似解KT0计算，在这种情况下，仅第一步迭代需要完全求解一个线性方程组，如果将KT0三角分解并存储起来，而以后各步迭代中采用迭代公式 则只需对上式右端项中的 进行回代就行了。这种方法称为修正的牛顿法。 返首页 使用修正的牛顿法求解非线性方程组，虽然 每一步迭代所花费的计算时间减少了，但迭代过 程的收敛速度也降低了。为了提高修正牛顿法的收敛速度可采用某些过量修正技术。 加速技术：修正牛顿法的过量修正技术－搜索办法 返首页 在算出 后，新的近似解由下式给出 （i= 1,2,···,N） 其中 是大于1的正数，它称为过量修正因子。 确定 的一维搜索办法。将 看做N 维空间中的搜 索方向，我们希望在该方向上找到一个更好的近似值， 即找到一个 式中的最好的 值。 虽然沿这一方向，不能期望求得精确解，但我们可以 迭择因子 (在搜索问题中称为步长因子)，使在搜索 方向上 的分量为零，即 上式是一个关于 的单变量非线性方程。通常可用一些 比较简单的方法来估算出 的大小。 加速技术：Aitken加速技术 返首页 在算出 后，新的近似解由下式给出 其中 是大于1的正数， 它称为加速因子。 1.2 牛顿法和修正牛顿法 其中μn的作用是改变切线矩阵KT的主对角元素，使奇异性或病态得到改善。。 此外，在某些非线性问题(如理想塑性和软化塑性问题、塑性卸载)中用牛顿法，迭代过程中切线矩阵可能是奇异的或病态的，为了克服这一现象，可有多种处理方法，其一是按下式来求 使用某种算法的计算效率，除了与收敛速度有 关外，还与每一步迭代所花费的计算量有关。关于 每步的计算量，牛顿法最大，而修正牛顿法最小。 因此在实际问题的计算中判断使用哪种方法效率较 高，往往需要进行数值实验。总的看来