计算固体计算力学第二章非线性方程组的解法.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级数展开，则可得 设荷载增量