上心盘热模压成形过程的模拟分析
. h) P" w6 t7 b5 f& V7 Y摘要:对铁路用上心盘热模压成形过程进行了数值模拟分析,得到了上心盘成形过程的应力场、应变场。对在实验中工艺参数选择及模具修正起到了重要指导作用。
! c* I6 f+ x" P( L; W关键词:模压成形;数值模拟
) E' J1 {: X0 @3 F" |
h$ A$ z; _4 v. k( |% t一、引言 & M1 M! r6 n* x1 G6 w1 z$ w- g
$ a: s E+ o8 w3 C; L1 s9 f
随着我国国民经济的发展,铁路运输能力也在不断增强,这就要求进一步提高铁路运行的可靠性、安全性和快速性。 5 ~' O1 E% x8 }* S
' j% |- B& L9 s4 i+ H$ n- D铁路货车上心盘是货车车体底部一个重要的支撑零件,在铁路运行中一直承受车体的整体重量,要求其具有较高的抗疲劳性能和整体综合性能。
( I+ x; ? C* {- `$ L3 g: F- M! @3 s; f( F4 [( K6 Z0 v2 K- y
某机车厂针对上心盘生产和使用中存在的问题将生产工艺从铸造生产改为水压机上的热模压成形。为深入研究该成形方式是否可行,本文采用有限元法对上心盘热模压成形过程进行了数值模拟分析。 ) u* r$ w- W G% c& X5 i
$ t( `5 f8 t; X+ I* G2 U' W! w: Y$ } X o4 P二、模拟分析模型的建立
1 t/ \) |' ^8 o; p3 Q r- Y3 c( y' V
) c) H+ [. u( l; K% E: h G上心盘从结构上具有对称性,以两个平面坐标轴对称,在长度方向上相差很小,热模压成形过程中金属流动也以两个坐标轴对称,金属成形过程中对中间环形部分具有几乎完全一致的流动规律。为提高计算效率、节省计算时间,取其轴向截面的二分之一来建立几何模型,该模型可以较完整地描述出金属在成形过程中的流动情况,如图1所示。 $ F/ {' ^* `0 g- R& D
& f$ N- n+ K \* r K从数值模拟分析的角度出发,借助基于修正的拉格朗日描述的弹塑性大变形有限元方法,忽略了变形过程中温度对变形的影响(变形速度较快),应用接触单元建立包括上、下模及坯料在内的整体模型,对整个成形过程进行数值模拟。。
+ |9 @' y/ j3 t" ?! t计算所用坯料尺寸为420×500×45mm的厚板,材料为25号钢,由于高温热态成形,为简化计算过程,材料特性取为各向同性双线性材料模式,此时材料的性能参数如表1。 表125号钢的材料性能参数
: J: @) ]' ?* r; d! G7 x4 f 温度(℃) | σs(MPa) | σb(MPa) | G(MPa) | E(GPa) | 20 | 260 | 470 | 195 | 210.8 | 1200 | 52 | 94 | 17 | 16.3 |
6 q* ~* n# s+ [取泊松比μ=0.3,摩擦系数f=0.12。
$ v6 c6 d% m7 U' R: f( ~7 t
. \0 n% b0 [; W; `考虑到坯料的尺寸较大,变形体采用平面四边形四节点单元作为基本单元,以保证坯料按所要求的形状成形。在整个模型中,共有二维实体单元336个,上模具接触单元13个,坯料上表面接触单元45个;下模具表面单元13个,坯料下表面与之对应的接触单元46个。整个模型共有453个单元,427个节点。 $ x- W" Q) {5 O6 L( y
/ N2 P; p+ i. N4 ]' Y由于取整个模具与坯料对称面的一半为模型,因此对于模具及坯料的左侧,只允许有上下方向的位移,不允许其左右移动,同时在下模的表面加上位移约束使其固定,这与热模压成形过程中下模具不动是完全一致的。 - C6 p+ H; i1 R8 l$ h+ g
0 T5 N" a& p- C2 h. F; V! B对所建立的有限元模型进行求解,其迭代步数为335步,迭代总次数为1952次,总的计算时间为6h42min。
7 m9 B7 T% y$ a. r9 U. e5 n1 M7 t+ ]& j
三、计算结果分析
% z/ r$ J6 B' s' B* Y
/ t+ D0 ~( Y0 F! i2 R0 O9 o1.应力场
* n0 u& x; A' Y: P8 U
8 ^+ J1 e. f3 E; _2 l: c图3~图5是在上述条件下计算得到的X方向应力、Y方向应力和等效应力分布规律。从应力场分布规律可以看出,坯料在成形过程中能较好地充填模具型腔,成形性较好,没有充不满、折叠等缺陷生成。在上心盘230mm圆环面外端R30mm圆角处,X、Y两个方向上受到较大的压应力,应力平均值为100MPa;在125mm圆环面上R20mm圆角处沿X方向出现了拉应力,但未达到破坏的程度。对于疲劳实验易产生裂纹的300mm圆环面底部R20mm处,在两个方向上均为压应力,不会产生裂纹。 $ \( Z/ @8 X9 b1 E
8 w1 I+ u$ `. ]( R6 I
2.应变场
L; _1 q, t/ e2 f! A( @2 x
% G! y& Q7 z+ I3 p2 C图6~图8是计算得到的X方向应变、Y方向应变及等效应变分布规律。从应变分布可以看出,各处应变差别较大,上心盘中心部位受X方向应变较大,即沿X方向伸长较多,而在140mm圆环面处出现了压缩应变。在Y方向上,上心盘中心部位受压缩变形最大,但由于中心部位受三向压应力状态,不存在产生裂纹的危险性。 9 T( B/ E. _# D1 I3 a
4 M/ l5 Q% E1 f; ~3 r7 o$ V
# t& t! d) o4 S
从图8中可看出,整个锻件等效应力分布很均匀,差值很小,因此卸载后不会产生很大的应力集中。对25号钢上心盘成形过程,由于材料本身具有很好的成形性能,且是热模压成形,因此成形后锻件残余应力很小,不会影响到上心盘的使用寿命。 1 W! t3 z0 |& F( p/ }$ ~; V! H3 R
5 @' B0 j. p5 y# D
四、结论 $ C& k0 g: F# W" k; w9 q3 }* a0 d
( ^. ~. o* F7 s# ]" F; N; |
对上心盘热模压成形过程进行了数值模拟分析,明确得到了整个成形过程中应力场、应变场的变化情况。对上心盘热模压成形时模具设计及工艺参数选择具有重要指导意义,生产的产品完全满足性能要求。 |