材料成形有限元法基础.ppt

板料成形有限元法—薄膜单元 局部坐标系以1节点为坐标原点，x轴与单元1、2边重合并指向2节点，z轴与单元法向量nz平行。单位法向量nz为 板料成形有限元法—薄膜单元 随体局部坐标系oxyz与空间整体坐标系OXYZ之间坐标转换矩阵 ?为 板料成形有限元法—薄膜单元 单元局部坐标自由度向量ue与整体坐标自由度向量Ue的变换关系为 板料成形有限元法—薄膜单元 四边形薄膜单元 板料成形有限元法—薄膜单元 局部坐标系以1节点为坐标原点，x轴与单元1、2边重合并指向2节点，z轴与单元法向量nz平行。单位法向量nz为 板料成形有限元法—薄膜单元 随体局部坐标系oxyz与空间整体坐标系OXYZ之间坐标转换矩阵 ?为 板料成形有限元法—薄膜单元 单元局部坐标自由度向量ue与整体坐标自由度向量Ue的变换关系为 板料成形有限元法—单元平衡方程 单元平衡方程 单元插值关系(局部坐标系) 单元几何关系(局部坐标系) 单元本构关系(局部坐标系) 最小势能原理(局部坐标系) 坐标变换关系 板料成形有限元法—单元平衡方程 代入局部坐标系的单元3个关系式, 得 代入坐标转换关系式, 得 板料成形有限元法—单元平衡方程 ke 为局部坐标系的单元刚度矩阵 f e 为局部坐标系的单元载荷向量 Ke 为整体坐标系的单元刚度矩阵 F e 为整体坐标系的单元载荷向量 板料成形有限元法—本构方程 在金属塑性大变形有限元分析时经常采用流动理论本构方程，其它本构方程很少采用 基于形变理论或非经典的角点理论本构方程虽然可以比较准确板料失稳后的局部化变形过程，但是板料成形属于强约束过程，对角点本构方程不敏感，而且板料成形也并不十分关心板料失稳后的局部化变形过程 目前板料成形有限元模拟的精度也不高，也没有必要十分强调本构方程的影响 板料成形有限元法—本构方程 考虑具有光滑屈服面屈服函数的弹塑性体，假设温度对变形速度的影响很小，可以忽略不计。这样全应变率可以分解为弹性应变率和塑性应变率之和，即 板料成形有限元法—本构方程 采用Hooke定律，弹性应变率 Piola Kirchhoff应力的Jaumann速率 用第二 表示为 塑性应变率 用流动法则和屈服函数 f 表示为 板料成形有限元法—本构方程 nij是屈服面（应力空间f = 0的曲面）的单位法向量 ? = 0，应力处于弹性状态, 应力点位于屈服面以内 ? = 1，应力处于塑性加载状态,应力点位于屈服面上 h 表示当前状态的加工硬化率 板料成形有限元法—本构方程 式中 应变率 板料成形有限元法—本构方程 由 与Cauchy应力的Jaumann速率 的关系 如果材料不可压缩 材料本构方程为 板料成形有限元法—J2流动理论 J2流动理论 Mises屈服函数为 表示应力 的偏量 板料成形有限元法—J2流动理论 材料的J2流动本构关系 式中 h 可以由单向拉伸实验确定 板料成形有限元法—J2流动理论 E 为弹性模量 Et 为单向拉伸真应力—对数应变曲线的切线模量 K 为材料强化系数; n 为材料强化指数 板料成形有限元法—屈服函数 Hill正交各向异性函数 一般的，若把各向异性主轴作为随体坐标系的x,y,z轴，则Hill屈服函数可以表示成 F、G、H、L、M、N为各向异性参数，由实验确定 板料成形有限元法—屈服函数 在平面应力状态下 Hill各向异性函数 板料成形有限元法—迭代解法 非线性方程组迭代解法 板材冲压成形数值模拟是一个强非线性问题，涉及到几何、材料和边界三重非线性。如果采用隐式有限元法就要求解非线性有限元方程组。非线性方程组一般是采用线性化方法，通过一系列线性解逼近非线性解。但是这种方法是有局限性的，而且有时解的漂移误差很大。因此，一般都采用迭代法求解非线性有限元方程组，即Newton-Raphson法 有限元基本方法—总结 ? 为弹性体的应变 ? 为弹性体的应力 u为弹性体的可容位移 弹性体处于平衡状态时，其势能应为最小 P 为表面分布力 G 为弹性体的体力 A 为弹性体的表面积 V 为弹性体的体积 有限元基本方法—总结 单元插值关系 单元几何关系 单元本构关系 N为单元形函数矩阵 L为单元几何微分算子 为单元弹性矩阵 单元节点自由度向量 有限元基本方法—总结 B 称为应变矩阵 单元平衡方程或单元刚度方程 k 称为单元刚度矩阵 f 称为单元载荷向量 单元刚度矩阵的特性 对称性 奇异性 主元恒正且对角占优 有限元基本方法—总结 单元模型构造 整体坐标系法 局部坐标系法 采用多项式函数 易于进行微分和积分运算 随着多项式函数阶次的增加，可以提高有限元法的计算精