[工学]材料加工力学基础-3.pptVIP

u R u1 u2 △R1 方程的求解 一般不能直接求解，因为在刚度方程中含有 迭代求解方法 第五章 刚塑性有限元法：基于罚函数法的刚塑性有限元法 将变形过程分成多个增量步 对于第i个增量步 设速度场初值 由线性化的整体平衡方程 求解速度场修正量 修正构形，求解其它物理量 Newton-Raphson迭代法 减速因子 引入减速因子的原则： 该因子必须小于1； 迭代求解过程的收敛判据 依据速度范数的判据 依据节点力不平衡量的判据 采用迭代法时必须设定迭代收敛判据，以决定迭代过程何时结束。 定义 过大则精度差；过小则时间长。 若解答是精确的，则： 依据单元变形能变化率与节点力作功率的关系： 定义收敛判据： 第五章 刚塑性有限元法：基于罚函数法的刚塑性有限元法 应变速率、应变、应力、构形的求解 第五章 刚塑性有限元法：基于罚函数法的刚塑性有限元法 速度收敛后，得到变形体的速度场分布 单元应变速率： 单元应变增量： 单元应变总量： 单元应变总量： 单元应力： 构形的求解： 在t 时刻，构形坐标为 ，在 时刻： 粘塑性材料定义 第五章 刚塑性有限元法：刚粘塑性有限元 流动应力与应变速率有关： ，甚至 对于刚粘塑性材料， 如设流动应力 代入 和 的相应表达式，即得到粘塑性材料的平衡方程。 单元平衡方程 第五章 刚塑性有限元法：基于Lagrange乘子法的刚塑性有限元 在单元层次上，将 作为待定常数 思考问题：该方程在形式上有什么特殊之处？对角线元素不占优的方程组如何求解？ 单元平衡方程的线性化 第五章 刚塑性有限元法：基于Lagrange乘子法的刚塑性有限元 为速度场的一个近似解，与之对应的 为 , 是速度 设 场的一个修正量，并认为 为常数，则平衡方程的第一式成为 （忽略了高阶微量） 也即： 第五章 刚塑性有限元法：基于Lagrange乘子法的刚塑性有限元 再一次忽略高阶微量后，得到 第五章 刚塑性有限元法：基于Lagrange乘子法的刚塑性有限元 令 得到线性化后的单元平衡方程： 该方程是包含Lagrange乘子的平衡方程，由单元刚度矩阵组装成整体平衡方程以后，即可求解。 为什么要假设初始速度场？ 第1步和第n步初始速度场选取方法的区别。 假设初始速度场的几种方法 初等方法：依据滑移线法或上限法 近似泛函法： Lagrange 乘子法 罚函数法 修正的罚函数法 正负号的选取规则 第五章 刚塑性有限元法：初始速度场的假设 通过“近似泛函”和“简化”的有限元法求解初始速度场。 基于罚函数的初始速度场 第五章 刚塑性有限元法：初始速度场的假设 对于单元e： 组装 引入速度边界条件 求解初始速度场。 由此可见：泛函的构造思路是，(1) 使有限元方程直接线性化；(2) 使初始速度场近似满足体积不变条件；(3) 使外载荷作功率为正。 直接迭代法 基本思想： 基本方法： 第五章 刚塑性有限元法：初始速度场的假设 将非线性的材料本构关系近似线性化，从而得到线性的整体有限元平衡方程，并将其解答作为刚塑性有限元的初始速度场。 定义 为材料的粘性系数。将“率”形式的塑性流动方程改写成 在寻求初始解答时，将 作为常数，于是 和 成线性关系，塑性变形功率应该有因子1/2。 Lagrange乘子法： 罚函数法： 修正的罚函数法： 第五章 刚塑性有限元法：初始速度场的假设 由此得到三种不同的方法： 与原泛函的差别仅在第一项！ 令近似泛函的一次变分等于零，即得到满足体积不变条件和速度边界条件的速度场近似解，该近似解可作为刚塑性有限元的初始速度场。 基于罚函数的有限元方程 第五章 刚塑性有限元法：初始速度场的假设 由此可见：(1) 有限元方程已直接线性化；(2) 初始速度场近似满足体积不变条件。 粘性系数取法： 初始值： 为提高精度而迭代。 摩擦力作用特点 与相对运动方向有关（分流点和分流线的概念） 摩擦力不能超过剪切屈服应力 影响摩擦力的因素复杂 没有分流线（点）的情况 摩擦力方向：与相对运动方向相反 摩擦力大小： 库仑摩擦定律： 剪切摩擦模型： 第五章 刚塑性有限元法：摩擦边界的处理 1. 求解过程需要迭代； 2. 摩擦剪应力不能超过K。 m是摩擦因子 M=1表示粘着 时，计算有可能发生不稳定。 有分流线（点）的情况 分流点（线）两侧摩擦力方向相反； 分流点处摩擦力为零； 摩擦力大小仍是按剪切摩擦或库仑摩擦计算。 第五章