非线性问题解法及其在Duncan—Chang模型中的应用.pdfVIP

张 波：非线性问题解法及其在Duncan—Chang模型中的应用 95 非线性 问题解法及其在 Duncan—Chang 模型中的应用 张 波 (陕西理工学院土木工程与建筑系． 陕西 汉中 723001) 【摘 要】 本文主要探讨了非线性方程的求解方法，指出增量法是Duncan—Chang模型数值计算的一种 比 较合适方法 ，这种方法不仅可得到其加载过程的历史结果 ，便于绘制其变形过程的应力应变曲线，而且验证总子 步数为45时，即可满足精度要求。 【关键词】 非线性；增量法；Duncan—Chang模型 【中图分类号】 TU431 【文献标识码】 B 【文章编号】 1001—6864(2011)08—0095—03 非线性问题一般可分为几何非线性、状态非线性 为凹性关系时的单变量问题可能会导致发散。 及材料非线性。一般将因大变形引起的非线性称为几 P 何非线性问题 ；接触、碰撞等则属于状态非线性 问题 ； 材料非线性问题一般是指材料的应力 一应变关系呈 现非线性特性，如岩土、土壤等材料，其应力应变关系 在一定范围内呈现出非线性弹性或弹塑性性质；实际 应用中，有些问题可能兼具以上的几种情况。另外，有 些材料还呈现出蠕变性质，这里不加讨论。计算中，对 qbo ’ 。妒 这些非线性问题都可归结到求解一个联立方程组，这 图1直接迭代法迭代过程 一 方程组的系数依赖于一些基本变量的某种函数，具 有如下形式： 1．2 Newton—Raphson法(简称N—R法) H4,+ =0 (1) 若令 ({}+。)H({咖}，){}+{f}=0 (3) 由于系数矩阵 依赖于未知数 或其导数，显 则在已获得方程 (1)的第r次的近似解 {咖}后 ， 然 ，式(1)就成为非线性方程组。这时，通常不能直接 为求得改进 的近似解可利用仅保 留线性项的泰勒 求解该方程组，必须采用迭代法进行求解 ，本文将对各 (Taylor)级数展开式 ： 种迭代方法作简单介绍，并讨论其在岩土 Duncan— Chang非线性弹性模型计算中的应用。 t4(t4,t)} …r)}+()rA…=0 1 非线性 问题的数值方法 (4) 1．1 直接迭代法(逐次逼近法) 则有：{咖}+。={咖}+△{}=0 (5) 这一方法是通过逐次求解来完成的，而在每次求 前式中d{}／d{6}为切线矩阵，即： 解时，总用前一次解得的未知量 的值求出系数矩阵 的当前值。其第次迭代的值可近似为： [酬 ] (6) 西 =一H[(西)] 厂 (2) 于是式(5)中△{}=一[ ] { }=一 如果迭代过程是收敛的，当r趋于无穷大时，西 的 [ ]I_(日({咖}){咖}+{) (7) 极限值就趋于真实解。为计算 日，要求先假设未知数