板材多点成形过程中成形力的数值模拟
- 格式:pdf
- 大小:368.52 KB
- 文档页数:8
" ( ( ! !! ! #! " ! !# " % % #! %) ", #)
$ ’ !
( ) ’
" 为壳单元中面对应点的坐标; / , 在 [ G’ , ] 区间内变 式中: ! ! & 为单元的厚度, ’ !F& ! ! 是参数坐标, 化; 下标% 为单元节点的编号; ( 为节点% 的形函数矩阵。 # " % 为单元节点处法向矢量, % ", #)
厚度的圆柱面、 球面及马鞍面三种典型形状的成形力进行了数值模拟, 分析了成形力的变化过 程及影响因素。结果表明: 在成形过程中, 随着基本体与板料接触数目的增多, 成形力会越来 越大。厚板成形、 硬质材料成形时都将需要较大的成形力。由于变形的复杂性, 成形双向曲率 的目标面时所需的成形力明显高于单向曲率。这些结果对于选择多点成形设备、 更好地控制 板材多点成形过程以减少成形缺陷具有重要的指导作用。 关键词: 板材; 多点成形; 成形力; 数值模拟 中图分类号: > ? * ) ! 文献标识码: .
・
・
[
]
相反。成形力作用于工具和板材相互接触的区域内, 要计算成形力就要计算接触界面的作用力。采用 有限元分析板材成形时, 工具作用到板材上的接触力是与板材上的有限元节点内力平衡的, 因此, 通过 计算出各单元上的节点内力的合力, 就可得到总成形力。 在有限元分析中, 单元. 的节点内力由如下公式计算:
文章编号: ( ) ! " # ! $ % & ’ # ( ) ) & ) ( $ ) ( ( " $ ) "
板材多点成形过程中成形力的数值模拟
宋雪松,蔡中义,李明哲
(吉林大学 辊锻工艺研究所, 吉林 长春 ! ) * ) ) ( (
摘
要: 基于动力显式有限元格式, 给出了多点成形中成形力的计算公式。对不同材料、 不同
第!期
宋雪松, 等: 板材多点成形过程中成形力的数值模拟
・! ! $・
点成形压力机实现的, 压力机工作时, 油缸的成形力通过活动横梁施加到上 (或下) 基本体群, 然后作用 到板材上。许多学者对板材多点成形数值模拟进行了研究。如采用静力隐式算法的多点成形过程分
[ ] [ , ] [ , ] ! ! " # $ % & 析 、 起皱与压痕缺陷的数值模拟 、 不同工艺过程的数值模拟 等等。但对于多点成形中的成
O 7 4 , * & P 7 , ( . E , ( ) ) & C
! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! ! !
・
) ( !
* 、 * 由节点法向矢 H H 式中: $ $ % 的速度分量; ) + + ( %) ( %、 %分别为局部坐标系下的节点 ( %、 %为" % 的变化速率; % , , , 量分量定义为: * * * * H H H -+ H H -+ H H -+ H % ’ ) + ( ’ ) ( ! ) ( 1 ) ( E ! E + 1 ’ + E ! + ’ -+ 1 , , , , , , , , ( ) 1 ! * * * * H H H H H H H H & ( + ) ( ! -+ ( E + ( 1 -+ ( ’ + ( E -+ ( ! + ( ’ -+ ( 1 ’ ) ! ) 1 ) E , , , , ! " 4 成形力的计算 板材成形过程中的成形力就是工具对板材的作用力, 其值与板材对工具的反作用力大小相等、 方向
(7 , * 8 89 * : 3 + / + ) ; 3 ; & ; ’, < 3 8 3 += + 3 > ’ : ) 3 ; 1 ? + @ 1 & +! * ) ) ( ( 1 3 + ?) , , 2, ,
: , 3 4 * + % ( ’ + 2 5 @ A7 8 B = @ A 8 2 1 3 : @ 6 4 3 : 3 B D 3 8 3 B @ @ 4 @ 1 @ 8 B 5 : = @ 1 @ B = @ D 7 E 1 3 8 D 7 E : @ 5 3 8 B = @ 5 = @ @ B1 @ B 2 41 F 4 B 3 + < C 9 ( ) , , H IC , @ , : 4 3 8 A E 3 : 2 4 7 3 8 B D 7 E 1 3 8 E 7 : @ 5 5J @ E @ D 7 E 1 F 4 2 B @ A , > = @ D 7 E 1 3 8 D 7 E : @ 5 D 7 E B = E @ @ B 3 : 2 4 5 = 2 @ 53 < C 9 G 9 < C C , 5 = @ E 3 : 2 4 2 8 A5 2 A A 4 @5 F E D 2 : @7 DA 3 D D @ E @ 8 B1 2 B @ E 3 2 4 5J 3 B =A 3 D D @ E @ 8 B B = 3 : K 8 @ 5 5J @ E @5 3 1 F 4 2 B @ A , > = @L 2 E 3 2 B 3 7 8 C = 3 5 B 7 E 3 @ 5 7 D B = @ D 7 E 1 3 8 D 7 E : @J @ E @ 2 8 2 4 ; @ A , > = @ E @ 5 F 4 B 5 3 8 A 3 : 2 B @ B = 2 B B = @ D 7 E 1 3 8 D 7 E : @M @ : 7 1 @ 5 4 2 E @ E 2 8 A 9 < 9 9 4 2 E @ E 2 4 7 8 3 B =B = @3 8 : E @ 2 5 @: 7 8 B 2 : B 3 8 F 8 : =8 F 1 M @ EA F E 3 8 = @D 7 E 1 3 8 E 7 : @ 5 5 , > = @D 7 E 1 3 8 7 E : @ 9 9J 9C 9B 9C 9D 3 8 : E @ 2 5 @ 5 @ L 3 A @ 8 B 4 3 B = B = @ 5 = @ @ B B = 3 : K 8 @ 5 5 2 8 A B = @ E @ 4 2 B 3 L @ 4 4 2 E @ D 7 E 1 3 8 D 7 E : @ 3 5 8 @ @ A @ A 3 D B = @1 2 B @ E 3 2 4 3 5 <J < 9 9 5 B 3 D D @ E , > = @ D 7 E 1 3 8 7 E : @ D 7 E B = @7 M @ : BJ 3 B =B = @5 F E D 2 : @7 DA 7 F M 4 @: F E L 2 B F E @ 3 57 M L 3 7 F 5 4 4 2 E @ E B = 2 8B = 2 B 9D N < 9 J 3 B = 5 3 8 4 @ : F E L 2 B F E @A F @ B 7 B = @A @ D 7 E 1 2 B 3 7 87 D B = @ D 7 E 1 @ E 3 51 7 E @: 7 1 4 @ 6B = 2 8B = 2 B 7 D B = @ 4 2 B B @ E , > = @ 5 @ 9 C E @ 5 F 4 B 5 2 E @ 5 3 8 3 D 3 : 2 8 B D 7 E : = 7 7 5 3 8 H IC E @ 5 5 2 8 A: 7 8 B E 7 4 4 3 8 H IC E 7 : @ 5 5 B 7 2 4 4 @ L 3 2 B @ D 7 E 1 3 8 @ D @ : B 5 , 9 9G 9G 9A : ; ; ; 5 $ , % 8 * 5 = @ @ B1 @ B 2 4 1 F 4 B 3 + 7 3 8 B D 7 E 1 3 8 D 7 E 1 3 8 D 7 E : @ 8 F 1 @ E 3 : 2 4 5 3 1 F 4 2 B 3 7 8 C 9 9 67
! " # $ % & ’ ( ) * & # " ) ( + & , -, . . , % # & . , % ’ $ * & -* 0 $ $ +# $ + ( ) / # " ) + & 1 , & + . , % # & % , ’ $ * * 2 /2
! " # $% & ’ ( ) * + . /0 1 * + ( 3, 4 /5 3 + ( 6 1 ’ ,, , 2 ,
形力目前还没有人进行系统的数值模拟研究。作者应用动力显式算法, 研究钢、 铜、 铝三种材料, 在圆柱 面、 球面及马鞍面三种典型形状多点成形过程中成形力的变化及有关影响因素。这些结果对于选择多 点成形设备、 更好地控制板材多点成形过程以减少成形缺陷具有重要的指导作用。
’ 成形力的计算方法
! " ! # $ % ’ ( ) * + , . , / 1 * 2 3 / & 0 0壳单元 (2 壳单元是基 ( ) * + , . 的) / 0 ( ( ’ # 1薄壳单元提供了多种不同的算法。其中 2 3 4 6 7 8 9 : ; * < 7 = <) 5 5 于完全平面的几何形状, 没有考虑翘曲问题, 因而尽管 2 但分析扭曲问题时 < 是一种非常高效的单元, 却导致不正确的结果。/ (/ ) 壳单元考虑了单元非平面几何形状的问题, 能给出比较好的扭 > 9 3 7 * ( @ > ( ? 曲计算结果, 但计算的效率比较低。2 (2 壳单元是对 2 (2 3 4 6 7 8 9 : ; * A ; B * C 9 @ = B AC) 3 4 6 7 8 9 : ; * < 7 = <) 5 ? ? 5 5 壳单元的一种改进, 可以进行扭曲分析, 而且, 增加的计算量也不大, 其计 2 AC 单元中引入了翘曲刚度, 算速度只相当于 2 < 单元的’ D ! "倍。2 AC 壳单元比较适用于复杂形状零件的板材成形数值模拟。 在2 AC 壳单元中壳面几何形状由下式插值得到: