工程有限元方法材料非线性要点PPT(41页)_详细.ppt

弹塑性本构关系 (c.3) 平面应力问题的应力应变关系 平面应力问题的应力应变分量 平面应力问题的弹性、塑性矩阵 其中 弹塑性本构关系 (3) 其它硬化材料的应力应变关系 (a) 理想弹塑性材料 (b) 运动硬化材料 (c) 混合硬化材料 屈服条件 应力应变关系的获得 弹塑性本构关系 2.4 弹塑性全量应力应变关系 (略) 全量理论的基本假设 应力主方向和应变主方向重合，且在整个加载过程中主方向保持不变 全量理论的适用性及特点 一般情况下很难满足全量理论基本假设 若实际分析中为单参数加载且变形满足小应变条件时，上述假设近似成立 全量理论使得整个分析包括本构关系的建立得到简化 弹塑性有限元分析 3.1 弹塑性有限元分析概述 对于小变形的弹塑性问题，几何方程与弹性分析相同，不同之处是单元的应力－应变关系可能是线性的或非线性的，因此集合单元得到的总刚度矩阵是与应力水平相关而出现的非线性方程组。 求解弹塑性问题，一般采用荷载增量法。当出现屈服点后，每次增加的荷载应适当减小。 弹塑性有限元分析 3.2 弹塑性问题的增量方程 由于材料和结构的弹塑性行为与加载及变形历史相关，在进行结构的弹塑性分析时，通常将载荷分成若干个增量，然后对于每一载荷增量，将弹塑性方程线性化，从而使得弹塑性非线性分析转化为一系列线性问题 增量方程 (1) 平衡方程 (2) 几何关系 (3) 物理方程 (4) 边界条件 其中 弹塑性有限元分析 3.3 弹塑性问题的增量有限元格式 虚位移原理（增量形式） 引入应力应变关系 矩阵形式 弹塑性有限元分析 有限元离散 单元内位移插值 应变位移关系 系统平衡方程 单元节点增量位移向量 型函数矩阵 应变矩阵 弹塑性刚度矩阵 不平衡力向量 其中 外载荷向量 内力向量 避免解的漂移 弹塑性有限元分析 增量有限元实施步骤 (1) 线性化弹塑性本构关系 (2) 形成增量有限元方程 (3) 求解增量有限元方程，获得增量位移 (4) 积分本构方程确定新的应力状态 (5) 检查平衡条件，若不满足平衡条件，则继续迭代 即检查外载荷向量和内力向量是否满足平衡 此处迭代可选择非线性方程组求解的欧拉法、N-R法或mN-R法 3 弹塑性有限元分析 每一增量步中基于N-R方法的迭代步骤 (1) 计算刚度阵和不平衡力向量 (2) 形成增量有限元方程 (3) 求解增量有限元方程获得增量位移，更新位移 弹塑性有限元分析 (4) 计算单元应变增量和应力增量修正量 (5) 检查平衡条件，若不满足平衡条件，则返回(1)继续迭代 常用的收敛准则包括位移、平衡力或能量收敛准则： 弹塑性有限元分析 上一增量步或上一迭代步结束时的应力 弹塑性有限元分析 克罗内克符号 * 布拉格；蔡格勒 * * 有限元方法与应用材料非线性 张有为 工程力学系 材料非线性要点 材料非线性问题分类 非线性弹性问题 弹塑性问题 弹塑性有限元分析 弹塑性本构关系 弹塑性有限元分析 材料非线性问题的来源 材料非线性的来源：体系的非线性是由于材料的应力与应变关系的非线性引起 增压器涡轮机轮盘叶片组件弹塑性变形 分析类型 特点 描述方法 应力和应变 仅材料 非线性 位移和应变无限小，应力应变关系是非线性的 仅材料非线性 工程应力 工程应变 材料非线性问题的来源 引例： 解：由于仅考虑材料非线性，因此a和b两段的应变可表示为 力平衡关系可表示为 此外本构关系可表示为 弹性加载 塑性加载 弹性卸载 材料非线性问题的来源 (1) 在初始加载阶段，由于外载荷较小，a和b两段将均处于弹性变形范围，有 (2) 随着载荷的增加，a段均处于弹性，b段处于塑性，即 由加载曲线可知a段将一直处于弹性变形阶段 材料非线性问题的来源 (3) 载荷达到峰值后，系统将处于弹性卸载，即有 外载荷与交界处位移之间的关系曲线如下图所示 材料非线性问题分类 非线性弹性问题 应力与应变的关系为非线性，且加载与卸载应力应变间的对应关系相同。 弹塑性问题 加载的过程中同时产生可恢复的弹性变形和不可恢复的塑性变形。弹塑性应力和应变间不再保持一一对应的关系，应变不仅依赖于当时的应力状态，而且还依赖于整个的加载历史。 弹塑性本构关系 2.1 材料的弹塑性性质 (1) 单轴拉伸实验 理想弹塑性材料 应变硬化材料 低碳钢拉伸应力应变曲线 颈缩阶段 强化阶段 屈服阶段 弹性阶段 弹塑性本构关系 (2) 几点说明 由于弹塑性的应力－应变关系不是一一对应的，因此研究弹塑性问题时，只有在确定的加载（或卸载）条件下才有明确的意义。 为了避免应力应变间的多值性带来的困难，不宜追求全应力与全应变之间的全量本构关系，应建立在一定加载路线条件下的增量关系。 弹塑性本构关系 为简化分析，结构工程中可采用理想弹塑性模型和弹性线性