金属刚粘塑性变形的能量泛函与动力显式算法有限元分析.pdfVIP

栅 搀 Vl June 1998D1．19No．2 金属刚粘塑性变形的能量泛函 与动力显式算法有限元分析 弋 ＼l3 25 解春雷 李尚健 黄树槐 (华中理工买军 模报复覆互 家重点实验室，武汉，43007) D』 摘 要 设计出刚粘塑性动态变形场的能量泛函．采用罚函数法和Lagrangian乘子法取消 对运动容许速度场满足体积不可压缩的约束条件，运用虚功原理和广义变分原理，推导出刚粘塑 性变形的动力分析显式算法有限元方程， 及速率形式的中心差分球解的时间积分显式．给出了 关t词 苎兰 鬯 竖有限元／： 鸳 锯 1 引言 金属成形过程的有限元理论和模拟技术已用于分析各类金属成形问题，刚塑性有限元 分析方法是解决体积成形问题的最有效的方法_1 ]．近年来，动力显式有限元方法的研究和 应用取得了较大的进展口 ]，动力显式方法具有 算法简单}计算机开销小；不需要迭代求解 高次非线性方程组，避免了臆式算法中的迭代求解不收敛的问题；若采用一点积分方法，又 可太大减小计算工作量．现有的动力显式有限元算法，是针对一般变形情况的，适宜于弹性 和弹塑性变形问题，如果采用弹塑性有限元模拟高温条件下的金属体积成形过程，就需要金 属在高温下的弹性模量和泊松比等材料特性．然而，对于体积成形问题，特别是热体积成形 问题，一般地说，弹性变形和回弹是可以忽略的，采用刚粘塑性有限元方法时就不需要弹性 模量和泊橙比参数．为此，本文探讨了金属体积成形的刚粘塑性变形动力分析显式有限元算 法． 2 刚粘塑性动态变形场的能■泛函 2．1 问题的定义 设刚塑性变形体 ，在f时刻的现时构形上，其表面积S—S +S+＆， 是外力 已知 表面}S 是与另一物体接触表面．此表面力为 ；S 是已知边界速度 表面}如图1所示．质 点的空间位置坐标记为 ，则位移、速度和加速度分别为 *一 一 ，如一主， 一 ( 一 1，2，3) 变形体 上．在塑性变形医满足刚塑性方程和边界条件“]． 2．2 能量泛函 运用不完全广义变分原理，设计刚粘塑性动态变形场的能量泛函 ，在变形场满足几 何方程和应力边界条件下，当她 一0使泛函取得驻值时，所得变分方程的解也为真实解．下 · 国家自拣科学基金资助项目． 1997—03—12收到第1稿．1997—10—23收到修改稿 第2期 解春雷等：盘属刚粘塑性变形的能量泛函与动力显式算法有限元分折 ·1 57· 面分别采用Lagrang[an乘子法和罚函数法，取消对 运动容许速度场 的体积不可压缩的约束条件． (1)能量泛函 采用Lagrangian乘子法建立能量泛函 ，在 满足变形几何方程和位移边界条件的完全速度场 中，使泛函 ’ 一 』 { 如d +』 {砘 』vf；,~dV — J 玉ds+J (南)d +J (1) 取驻值，即嘶 —o的速度场 ．及其所对应的应力 - 场 一定是刚粘塑性变形动力问题的真实解． 图1问题定义 (1)式中右端第一项的物理意义是完全速度场引起的物体的动能变化率，{；表示物质导数； 为质点速度，下标i表示i方向分量；P为材料密度；y为阻尼系数； 为体积力；E( )为功 函数，满足 一 ／ ／为Cauchy应力偏张量；毫J， 为应变速率张量和体积应变速率；^为 Lagrangian乘子．