非线性方程组的解法.pptVIP

钢筋混凝土结构的非线性有限元分析 莫 婵 主要内容 非线性方程组的求解方法 1.1 结构分析的非线性问题 1.2 非线性方程组的求解方法 欧拉折线法 增量法 修正的欧拉折线法 割线刚度迭代法 迭代法 切线刚度迭代法 等刚度迭代法 1.3 收敛准则 1.4 实例比较 1.1 结构分析的非线性问题 1.2 非线性方程组的求解方法 不论是材料非线性问题还是几何非线性问题，经过有限元离散后，最终都归结为求解一个n个变量，n个方程组的非线性方程组： —— 平衡条件 其中[K]为总刚度矩阵；[ ]为节点位移列阵；[P]为节点荷载列阵。 总体刚度矩阵可由单元刚度矩阵按标准方法集合而成： —— 几何关系 —— 本构关系 在线弹性材料中， 是常量；在非线性问题中， 是应力的状态函数，即： 1 欧拉折线法 将荷载P分为m个增量， 每一个荷载增量产生一个位移 ，因而 施加n个荷载增量后，总荷载为 。 计算第n个位移增量公式： 计算总的位移与应力： 2 修正的欧拉折线法 利用每一步荷载增量的始，末刚度的某种加权平均值替代起 始刚度来计算本步荷载增量的位移。 其中， 为加权参数， ，一般可取 。然后利用 来推求 ，利用这一刚度来求本步荷载增量下的位移，即 当取 时，就相当于求常微分方程数值解的中点龙格-库 塔法，全部公式为： 1 割线刚度迭代法 在某级荷载P作用下，用初始刚度 矩阵[K0]求得位移的第一次近似值： 然后利用 求得单元的应变，进而求得应力。根据应 力状态确定即时的本构矩阵，根据这一本构关系矩阵即 可求得新的割线刚度矩阵[K1]可求得位移的第二次近似值： 重复上述步骤，每次可由下列公式求得进一步的近似值： 直到 与 充分接近为止。 2 切线刚度迭代法 首先取初始刚度矩阵[K0],求得位移的 第一次近似值 由初始位移可以求得单元应变，进而 求得单元应力。由单元应力可以求得相 应的节点荷载[P1]。第二步，用相应于 时的即时切线模量[K1] 在荷载 作用下求得位移增量 ，即 进而求得位移的第二次近似值为： 重复上述步骤，即 3 等刚度迭代法 计算过程如右图所示： 具体步骤如下： （1）用初始刚度[K0],求出位移的第一 近似值 （2）按 求出单元应变 ，由单元应变求得单元应力 ，由应力可以求得相当的节点力为 这样 与原加荷载的差为 （3）将 再加于结构，仍用初始刚度[K0]求得附加位移 ，从而求得第二次位移的近似值 （4）重复以上的步骤，直到 与 之差达到足够小或者 足够小为止。 1.3 收敛标准 常用向量的三个范数为： （1）采用结点位移序列收敛标准时，若所求得的位移序列为： 则收敛标准可表达为： （2）采用不平衡结点力序列收敛标准时，若所求得的不平衡结点力序列为： 则收敛标准可表达为： 1.4 实例比较 【例1】 材料本构关系为 试用增量法求解位移？ 解： （1）弹性解 （2）非线性精确解 由 得 （3）欧拉折线法求解 用有限元法求解时，刚度矩阵只有一个元素，即为 由 ，可得切线弹性模量为： 割线模量为： 相应地，全量刚度（割线刚度）为： 增量刚度（切线刚度）为： 用欧拉折线时，将荷载分为3个增量， 为第一级荷载增量时， ，因而有 将计算结果与精确解相比较，发现结果变小。 改用中点刚度法计算