计算固体计算力学 - 第二章 非线性方程组的解法汇总.pptVIP

* * * * * * * * * * * * * * * * * 求解非线性方程组的另一类方法是增量方法。使用这种方法需要知道“荷载”项(R)为零时问题的解(a)0。在实际问题中，(R)经常代表真实荷载，(a)0 代表结构位移。在问题的初始状态，它们均为零。这种从问题的初值开始，随着荷载列阵(R)按增量形式逐渐增大，研究(a)i的变化规律的方法，称为增量方法。 当问题的性质与加载的历史有关时，如弹塑性问题，则必须采用增量方法。 2.4 增量方法 设“荷载”(R)在任一增量步的值为λ(R)，λ 为荷载增量因子，(R)为标准荷载列阵,则非线性方程Ψ(a)= 0可写为 引入切线矩阵且略去高阶小量后可改写为 Ψ(a,λ)=P(a)-λR=0 若λ+Δλ时的解答为a+Δa，象牛顿法一样，将Ψ (a+Δ(a)),λ+Δλ) 按Taylor级数展开，则可得 设荷载增量因子λ分别取如下值 am+1=am+Δam 则荷载(R)可分成M级，第m级荷载为λmR，其增量为(λm+1-λm)R=ΔλmR。 由此可得：Δam=[KT(am,λm)]-1ΔλmR 在以上各式中，下标m表示增量步的步数，而λm=1的解对应于Ψ(a)=P(a)-R=0的解。 在一个自变量的情况下，求解非线性方程组的过程如下图所示。如果Δλm足够小，则认为所得解即为方程组的合理的