非线性有限元.pptVIP

非线性有限元 * 非线性问题 大多数实际问题属于非线性问题，根据产生非线性的原因，非线性问题主要有三种类型： ●材料非线性(物理非线性) ●几何非线性 ●接触非线性 材料非线性 应力与应变之间为非线性关系，通常与加载历史有关，加载和卸载不是同一途径，因而其物理方程中的弹性矩阵是应变的函数。但材料非线性问题仍 属于小变形问题，位移和 应变是微量，其几何方程 是线性的。土、岩石、混 凝土等具有典型的材料非 线性性质，应当按材料非 线性问题处理。 几何非线性 几何非线性属于大变形问题，位移和应变或者它们中一个是有限量。可能会有三种情况:大位移(包括线位移和角位移)、小应变；小位移、大应变；大位移、大应变。此时反映应变和位移关系的几何方程是非线性方程， 如果应力和应变之间的关系也是非线性的，就变成了更复杂的双重非线性问题。在几何非线性问题中一般都认为应力在弹性范围内，应力与应变之间呈线性关系。工程中的实体结构和板壳结构都存在几何非线性问题，例如弹性薄壳的大挠度分析，压杆或板壳在弹性屈曲后的稳定性问题。 接触非线性 由于接触体的变形和接触边界的摩擦作用，使得部分边界条件随加载过程而变化，且不可恢复。这种由边界条件的可变性和不可逆性产生的非线性问题，称为接触非线性。 材科非线性有限元法 材料非线性是由本构关系的非线性引起的。但它和线弹性有限元一样，都属于小变形问题，因而关于形函数的选取、应变矩阵、应力矩阵及刚度矩阵的形式都是相同的，不同的仅在刚度矩阵是按非线性弹性或弹塑性矩阵计算的，这是材料非线性有限元的基本内容。 全量形式的非线性弹性本构方程： 全量形式非线性有限元 D为弹性矩阵 平均等效应力 平均等效应变 和 之间存在单值函数关系 关系由试验确定，对于简单拉伸，就是单轴的关系。 单元切线刚度矩阵 可得整体平衡方程 整体切线刚度矩阵 与位移 有关， 这是一个非线性方程组。一般用迭代 法求解。 整体平衡方程 （一）全量形式非线性有限元求解方法 直接迭代法 初应力法 初应变法 (1) 直接迭代法 对非线性方程组 设其初始的近似解为 ，由此确定近似的 矩阵 可得出改进的近似解 重复这一过程，以第i次近似解求出第i＋1次近似解的迭代公式为直接迭代法 对非线性方程组 直到 变得充分小，即近似解收敛时，终止迭代。 在迭代过程中,得到的近似解一般不会满足 作为对平衡偏离的一种度量，称为失衡力。 直接迭代法的计算过程 设 为具有一阶导数的连续函数， 是方程的第i 次近似解。若 在 附近按一阶Taylor级数展开 Newton—Raphson方法 引入记号 Newton—Raphson迭代公式为 Newton—Raphson迭代法的计算过程 但 Newton-Raphson 法不能保证在所有情况下都收敛! 仅当初始构形在收敛半径 内时 Newton-Raphson 才收敛. F δ 位移 载荷 收敛半径 如果 δ初始 在收敛半径内, 解将收敛; 否则解发散. δ初始 ? F δ 位移 载荷 δ初始 发散! F δ 位移 载荷 δ初始 收敛 初始点在收敛半径外部 初始点在收敛半径内部 F δ δ start 如果初始构形在收敛半径外部, 有两种技术可帮助获得收敛解: 递增加载使目标更接近初始点 F δ δ start F1 用收敛增强工具扩大收敛半径 通常结合两种 策略获得收敛. 一般的规律是系统任何方面的突变会导致收敛困难. 刚度突变. 载荷突变. 最佳收敛行为是把突变分成一系列很多小的递增的变化. 采用渐变加载. 采用小的时间步. m-Newton-Raphson法 m-Newton—Raphson迭代公式为 m-Newton—Raphson迭代法的计算过程 q-Newton-Raphson法 每次迭代后用一个简单的方法修正 的修正要满足以下的拟牛顿方程 DFP(Davidon-Fletcher-Powell)公式 BFS（Broyden-Fletcher-shanno）公式 *