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

哈尔滨建筑大学王焕定教授制作第1页，共28页，星期日，2025年，2月5日非线性问题可分为三类：材料非线性、几何非线性和边界非线性。我们只讨论前两类问题。不管那类非线性问题，最终都归结为一组非线性方程Ψ(a)=0，a为待求的未知量。Ψ(a)=0可写成平衡方程的形式Ψ(a)=P(a)-R=K(a)a-R=0对非线性方程Ψ(a)=0，一般只能用数值方法求近似解答。其实质是，用一系列线性方程组的解去逼近所讨论非线性方程组的解。分段线性法第2页，共28页，星期日，2025年，2月5日1.1直接迭代法1.2牛顿法和修正牛顿法1.3增量方法1.4增量弧长法第3页，共28页，星期日，2025年，2月5日1.1直接迭代法Ψ(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)视为不平衡力并作为衡量收敛的标准第4页，共28页，星期日，2025年，2月5日对于单变量问题的非线性方程，直接迭代法的计算过程如图1所示。图上给出的是和之间的关系，而不是和，之间的关系对单变量情况，直接迭代实质是“割线”法，一定条件下这种迭代过程是收敛的，但对多自由度情况，由于未知量通过矩阵K(an)的元素互相耦合，在迭代过程中可能会出现不稳定现象。a1=K(a0)-1R直至Δan=an+1-an满足收敛条件第5页，共28页，星期日，2025年，2月5日收敛性第6页，共28页，星期日，2025年，2月5日1.2牛顿法和修正牛顿法如果将非线性方程Ψ(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)-1Pn，an+1=an+Δan第7页，共28页，星期日，2025年，2月5日Δan≈-KT(an)-1Pn，an+1=an+Δan直至Δan满足收敛性第8页，共28页，星期日，2025年，2月5日1.2牛顿法和修正牛顿法如果在迭代计算的每一步内，矩阵KT都用初始近似解KT0计算，在这种情况下，仅第一步迭代需要完全求解一个线性方程组，如果将KT0三角分解并存储起来，而以后各步迭代中采用迭代公式则只需对上式右端项中的进行回代就行了。这种方法称为修正的牛顿法。返首页第9页，共28页，星期日，2025年，2月5日使用修正的牛顿法求解非线性方程组，虽然每一步迭代所花费的计算时间减少了，但迭代过程的收敛速度也降低了。为了提高修正牛顿法的收敛速度可采用某些过量修正技术。第10页，共28页，星期日，2025年，2月5日加速技术：修正牛顿法的过量修正技术－搜索办法返首页在算出后，新的近似解由下式给出（i=1,2,···,N）其中是大于1的正数，它称为过量修正因子。确定的一维搜索办法。将看做N维空间中的搜索方向，我们希望在该方向上找到一个更好的近似值，即找到一个式中的最好的值。虽然沿这一方向，不能期望求得精确解，但我们可以迭择因子(在搜索问题中称为步长因子)，使在搜索方向上的分量为零，即上式是一个关于的单变量非线性方程。通常可用一些比较简单的方法来估算出的大小。第11页，共28页，星期日，2025年，2月5日加速技术：Aitken加速技术返首页在算出后，新的近似解由下式给出其中是大于1的正数，它称为加速因子。第12页，共28页，星期日，2025年，2月5日第13页，共28页，星期日，2025年，2月5日第14页，共28页，星期日，2025年，2月5日1.2牛顿法和修正牛顿法其中μn的作用是改变切线