板料成形数值模拟中的弹塑性有限元方法.pdfVIP

板料成形数值模拟中的弹塑性有限元方法 杨玉英徐伟力金朝海挚郭刚 (哈尔滨工业大学材料科学与工程学院 15000n r争深圳高等职业技术学院 518055) 摄耍本文舟绍了板料成形数值模拟中的弹塑性有限元方法的单元选取、平衡方程建立、本构关系选取等技 术问题t并用独立开发的静力隐式的弹塑性有限元程序对几个比较典型的板料成形过程进行了数值模拟，有 助于从量化的角度加深对板料成形规律的认识 ‘ 1引吾 自从山田嘉昭从本世纪60年代推导出小变形情况下的弹塑性本构矩阵以来fl】，弹塑性有限元 方法无论是在关键算法上，还是在本构关系、屈服准则等技术问题上都取得了较大的进展。 尤其是在板料成形领域，它正逐渐从实验室阶段走向实用化阶段Ⅲ．本课题组从90年代初就 开始进行板料成形过程数值模拟软件的开发工作·现已分别开发出两维、三维的弹塑性有限 元模拟软件。本文简要介绍一下板料成形数值模拟中的弹塑性有限元方法的一些技术问题， 并利用独立开发的静力隐式的三维弹塑性有限元模拟软件对几个典型的板料冲压成形过程进 行了数值模拟。 2弹塑性有限元程序的基本理论 2．1单元的选取 在板料成形成形过程的数值模拟中．常用的单元类型有体元、壳元、膜元。本文采用考 虑横向剪切变形的Mindlin壳单元，它既克服了采用体元时若单元细化则计算代价昂贵，著单 元粗化则有局部过剐性的缺点，又克／ltlT膜元忽略弯曲效应，不能正确模拟起皱现象的缺点。 但是，Mindlin壳单元的引用，增加了节点的自由度，并使编程过程变得复杂。 2．2平衡方程的建立及本构矩阵的选取 板料成形过程在非线性领域基本上属于一个大位移、大转动、大应变的问题，也是非线 性领域内最复杂的一种情况。为此本文采用修正的拉格朗日列式方法(简称u．L法)来建立平衡 方程，也称弹塑性有限元方法的控制方程。 修正的拉格朗日列式法在求，+△f时刻物体的平衡解时，选择的参考构形是最后一个己 知平衡状态，即认为t时刻为“来”变形物体，去求f十△，时刻物体的各个未知变量。这时， 虚功方程为 I，“6js，，占”4jE，，dv=占“6‘W (1) 其中 H6tt。口=ttdu+型+6：S?，H6：Eq=o+公+6tE口 (1·a) 矿6忍=圭c券十警+等铮 (1-b) 肛。1=f，弘-17：(cj舻L1西 ：(芷眦]=j，：[BⅣ￡17：【仃l：【口帆j咖 E 弦}=j，．pc姒巧)咖 l(3- ；{盯]为柯西应力矩阵；；{口)为柯西应力矢量：{C)为本构矩阵；fR)为戴荷列阵；{ 黧j 位移增量列阵：IBL]、IBM】为单元形状函数中的线性和非线性部分，具体表达式视 异，可参考文献[3—61：右上标f、i一1为迭代次数。 对于本构矩阵．在小变形情况下．可直接引用标准的弹塑性矩阵，而对于大变形的情测，j 构矩阵需通过焦曼应力率和变形率之间的关系，通过一系列推导才能晟终得到。本文采用} 本构矩阵为 ，旧=，[D口卜。[∥d】 o 式中，，[Dep】为小变形情况下的弹塑陛矩阵，电即由山田嘉昭推导出的弹塑性本构矩阵。 ，f盯，，j为由柯西应力分量构成的矩阵，芭豹具体表达式为 2巧】