马上注册，结识高手，享用更多资源，轻松玩转三维网社区。
您需要 登录 才可以下载或查看，没有帐号？注册
x
汽车厚板料零件冲压成形分析及回弹计算
7 ~' ^, B8 E4 t( W' [% U ( f8 e: h3 ]5 ~0 i* l9 w
0 Z7 K/ t/ @* J6 ^2 P4 s$ S, j! c7 ^) s作者：中国第一汽车集团 富壮 王广盛
7 t; M+ |0 n/ }" v3 A
. n9 B1 \- G7 ^5 O1 p) b
2 k3 g/ _3 U; Z2 U# x! E' m W / R+ d+ d! X* `7 D
摘 要：汽车上板厚大于5mm 的厚板料零件的冲压成形CAE技术在材料、工艺、计算和评估等方面都与薄板料零件有所不同，基于MSC.Marc 软件并结合作者在厚板料零件冲压成形CAE 分析方面的实际工作，对计算模型建立时需注意的问题如单元选择、单元划分、屈服准则、硬化曲线、工况设定和回弹计算等进行了详细说明，并对厚板料零件上的伸长类翻边结构的成形极限问题进行了探讨。 , @9 [1 `5 u8 d1 V0 d& n
关键词：厚板料；冲压成形；成形极限；CAE ; I- t4 x8 s M2 L. E, D: A
0 S7 l! U& x! T+ Y. H7 ]
引言
/ {/ X6 B2 X, U2 {' H, D& I# U/ S+ J- b! ^; G. J9 {( k) g
随着我国汽车板料零件设计、制造水平的不断提高，薄板料零件冲压成形CAE 技术的应用已日趋成熟，相关产品的设计和制造部门针对不同软件及计算方法建立起了对应的材料、工艺、计算和评估方面的标准和规范。这些标准和规范经过实践的检验和修正，目前在产品设计和生产制造环节中得到了广泛应用。 4 F2 l1 ]/ |8 `" ]# ~, J
7 u3 g4 P$ n% _与薄板料零件不同，对于板厚大于5mm 的厚板料零件，例如商用车车架横梁、纵梁和加强板类零件，其在冲压成形、失效判定和回弹计算方面还没有一个明确的计算方法和分析思路，应用也远不如薄板料零件冲压成形CAE 技术广泛和成熟，这是与厚板料零件冲压成形的特点及其CAE 技术有关的。
" X2 k! H9 V$ r) F& M2 f9 Y$ r" @4 D- i% u7 c
目前国内针对这方面的研究相对少，这部分工作也有进一步研究和完善的必要，为此作者将近年关于厚板料零件冲压成形CAE 技术方面的工作进行了总结，并对其中一些具体问题进行了深入探讨。当然由于个人能力有限并且所面对问题又是行业内公认的“顽疾”，因此所做的工作远没有达到解决精确回弹计算的程度。
$ X ?% T+ s* H+ q9 T- P6 _1 V
9 N3 a7 c0 m4 i1 X0 w3 F+ z本文所讨论的相关内容都是基于MSC.Marc 平台的，选择MSC.Marc 软件除了非线性计算功能方面的考量外，更主要的是作者有十年以上该软件的使用经验，对于成形和回弹计算模型的精度和效率的控制有一定把握。
/ X( d$ @( v8 F* p& o3 u5 G& [: f" T2 M% O1 {. ?4 x
1 厚板料零件冲压成形及其CAE 技术的特点
, Z6 |0 h6 w3 s" p/ d8 l
2 h4 X. M c/ a: w与薄板料零件相比厚板料零件在冲压成形及其CAE 技术方面有如下特点： 3 }1 W# @% l0 W4 N0 \; {2 C6 N
' W; W3 N" `3 x; r w1)从材料方面讲，厚板料零件的材料一般采用热轧碳素钢板或热轧低合金高强度钢板。与冷轧薄板料相比，热轧厚板料的表面质量差、厚度公差大、材料力学性能不稳定，并且由于材料的延伸率较低，获得有效硬化曲线数据的应变范围较窄(工程应力达到抗拉强度之前的部分有效)，硬化曲线一般不采用真实应力应变曲线，而是采用幂指函数来替代。
; @; |0 B& F- y( ^! d- f8 r. s" ^9 h0 ^0 n' y h
2)厚板料零件一般是以梁类结构件为主，其成形方式主要为弯曲，并伴有局部的伸长类或压缩类翻边，少数情况下还包括胀形。成形时一般不采用压边圈，但都有背压垫，冲压工艺设计相对简单。单从判定零件冲压结构工艺性角度而言，通过经验数据和简单的一步法计算即可实现，因此一般不需要经过冲压成形CAE 评估。图1 为一汽某商用车车架纵梁的端部结构。 ) b* Q0 S5 o& d- M* A
' n8 y3 L3 f* a; d/ N
" g4 r3 c" }2 c4 q3 Q3 V
! g' J( A8 r) w
; O6 f2 r: B4 I& Z2 [) E虽然厚板料零件冲压成形CAE 评估的需求相对于薄板料零件要少得多，但随着厚板料零件材料强度级别的不断提高，回弹所造成零件尺寸精度的问题越来越严重，模具设计和后期的工艺调试都要求对零件回弹的性质及大小有所了解，以便采取相应的对策和补救方案。因此，厚板料零件冲压成形CAE 技术的工作重点在于对零件回弹趋势的准确判断和对其量值的合理估算。图2 为一汽某商用车车架纵梁的回弹计算结果。
, u5 k1 R0 d% `0 A$ y+ e* J0 Q4 Q
; i0 w# ~( X; W7 d
) |% Y* ?3 V: ~8 ^, X; T7 J
3)对于厚板料零件，其弯曲半径与板厚之比一般都很小，板厚方向的应力及其应力变化不容忽视，即弯曲部位的应力状态为三向应力状态。因回弹计算时弯曲部位的应力状态对回弹变形影响极大，故不宜采用壳单元进行厚板料零件成形及回弹计算。
, K m1 i# k' l) O: d5 ^8 l4 Q0 x( E5 v9 E( K! C8 l" k! t
另外，采用实体单元将会给动力显式算法带来麻烦。因为动力显式算法中的临界时间步长与单元上两节点间的最小距离成正比，采用实体单元后(一般料厚方向上的单元要分3 层以上)临界时间步长将大大减小，从而使计算费用大幅提升。并且由于弯曲成形时边缘较大面积的坯料处于自由状态，为了有效控制动力响应问题，凸模虚拟速度一般要小于2m/s，这将进一步增加计算时间，增大计算成本。 + e& d5 h2 n# z+ F' x. ]. \* E: ^
* |; q. x; C) A# m/ E1 Y7 {( }4)成形失效开裂判定准则是冲压零件能否通过成形性评估的主要依据，从目前情况看厚板料零件成形失效开裂判定准则还不十分完善，而对于薄板料零件成形的FLD 图则很难直接用于厚板料零件的成形性分析，尤其是冲裁后坯料边缘存在硬化和撕裂情况下的成形极限判定。
3 v: D; E0 J2 x& n$ P0 }) Q9 }; Z+ I( G2 A. T3 z
2 计算模型建立时需要注意的问题
: v/ H+ p7 V- W4 R8 C0 D" R7 C6 F" s- a6 Y4 D# }0 D. z p6 y
2.1 单元选择
" p' j' c, f9 H; H' f
3 a1 E6 ~" z$ ~4 s# J# @MSC.Marc 软件有两种单元供厚板料零件冲压成形CAE 分析时选用。一种是普通实体单元，另一种是实体厚壳单元。 ( I3 h9 t& M) k% r8 S, f5 v* ]; O$ G
0 y2 P/ M/ L @( d8 E5 V实体厚壳单元是基于选择性缩减积分和假定应变技术的[1]，由于是八节点拓扑结构，因此可以实现板料内外两侧的双面接触。
7 G, i% q8 z0 [$ W m [ V: \# l( A1 k& X; \# y2 b
实体厚壳单元采用平面内一点积分，而厚向的积分点数可由用户来指定。实体厚壳单元能够较好的分析弯曲圆角部位的塑性变形，但对于伸长类或压缩类翻边部位如果单元密度不够，则无法准确追踪塑性应变的梯度，这将直接影响到零件回弹计算中扭转变形的趋势和量值。另外，实体厚壳单元在使用Hill1948 和Barlat1991 屈服准则时可能会存在一些问题[2]。图3 和图4 为采用同等平面单元密度、不同单元类型时横梁伸长类翻边部位的应变分布和回弹计算结果。
" Q( ?3 |/ O' n4 R
8 B, ?7 q9 h9 {( y# F: u5 i' R4 k/ k
9 l- l9 T1 x1 q3 {8 G) @4 W+ L: U8 K8 ~- ]1 W
注：为了便于说明回弹变形，在显示设定上作了相应放大处理
7 k6 k/ S ?& z6 Y, q图4 采用同等单元密度不同单元类型时厚板料零件的回弹计算结果
4 m" J* z) D3 G% k$ O2 d2 J
针对上述情况，对于那些以弯曲变形为主、伸长类或压缩类翻边变形程度比较弱的情况，推荐采用实体厚壳单元。对于那些伸长类或压缩类翻边变形程度比较大的厚板料零件，分析时还是建议采用普通实体单元。 : Q. {" U' d# O! x& u
4 S f; t( h" z! ]5 ]9 w {/ |
通常情况下实体厚壳单元偏“软”，普通实体单元偏“硬”，实际回弹变形可能介于二个计算结果之间。 % j' m: A& P. l; Y
; l5 Y6 T+ S$ q
2.2 单元划分
; K0 V. t6 U' D( w* v! \9 ^8 b7 E% v+ o6 W5 l3 A0 o3 z9 }
计算模型建立时首先遇到的问题是坯料单元划分。对于厚板料零件成形用坯料单元的划分通常应保证：
& B" ], E0 k# A/ y/ o
; a! t' u1 T- E+ K. F1)弯曲圆角部位的单元其位向应与弯曲轴线垂直。这是因为三维实体单元为避免刚硬问题而采用了假定应变技术，而假定应变技术对单元的形状比较敏感[3]，为了保证成形计算时的应力结果的准确性，尽可能使单元位向与主应力方向一致。 4 N$ k/ }9 B3 r' y' ~+ h
$ P( k; M. o% |5 P+ E% B
2)弯曲圆角部位应至少保证四个单元，这对于弯曲圆角部位的应力梯度影响非常大。
% D" ?' S; [$ \, o9 z5 E6 f. ^: t* t& ^- {4 f3 ]
3)有拉伸和压缩变形的部位细分单元，即伸长类或压缩类翻边部位的单元应当细化。
8 g, ?; k; a8 y! X- o9 E c* l8 l; m
4)采用普通实体单元时，板料厚度方向分层推荐4 层，最少3 层。采用实体厚壳单元时，厚向积分点数为7。这一要求是为了准确描述回弹后零件厚向残余应力分布状态的，图5 为回弹后零件断面上的残余应力分布。 7 Z& D5 d) e1 E* M/ J L
6 R6 a" n" {* v
$ `$ D& L L: O7 O, K5 V' t% \
图5 回弹后的应力状态 ( H+ ^" P6 c8 E& U- K [4 Y. n
5)单元长边与短边的比值尽可能的小，一般要控制在3 以内[4]。这一要求是为了避免数值计算误差、确保计算精度的。
+ M5 A9 B5 p: @& @! s9 O, V: a" n5 f; L" x: P: x; O* o2 c t
2.3 屈服准则
+ z" g4 F9 H& g- P8 D4 B( H: }- Q3 d3 F8 e
MSC.Marc 软件中用于板料成形的屈服准则有Von Mises、Hill 1948 和Barlat 1991，对于钢板类冲压零件经常使用的是Hill1948 屈服准则。
+ [2 Y$ a* x- |8 b- V' t$ r, {1 A6 Y5 @4 F' w
MSC.Marc 软件需输入的 材料各向异性参数有R0、R45、R90、Y45 /Y0 和Y90 /Y0。其中Y45 /Y0和Y90 /Y0 可按式(1)计算。
9 G! ^8 R3 C3 j6 A1 ]
0 h+ [, w a6 j% S% w* }; q; Z2 R" V9 W; K. G* I' y; U
" m! s7 r: s# J2 \- G定义坯料轧制方向时，对于普通实体单元只能利用3D ANISO 方法进行设定，其中第一矢量方向为轧制方向，第二矢量方向为坯料平面内垂直于第一矢量的方向。对于实体厚壳单元，可以利用XY PLANE、YZ PLANE 或ZX PLANE 方法设定轧制方向，轧制方向为坯料平面与所选中平面的交线。 8 Z) A# R7 F w" a( S8 x& J H. j- p# l
# j% y9 O' ?* s2.4 硬化曲线 0 \* n) ?4 i" w5 w; \( p& C# I
" k, D: o1 A3 c5 q& N需要进行回弹计算的厚板料零件成形分析，材料参数最好输入真实硬化曲线。因为弯曲截面上应力分布与硬化曲线形状相一致，如图6 所示，为保证弯曲截面上的应力分布更加准确，推荐输入材料硬化曲线。 - z( b8 ]2 K0 I* Z9 F5 O5 D
) |4 a/ c8 U1 {6 R2 t6 T3 T
4 v a3 M( V- U4 a5 \; G2 u; N3 J2 D- J% b7 U O3 K3 h; I$ F
但是考虑到材料延伸率较低，获得有效硬化曲线数据的应变范围比较窄，硬化曲线采用真实应力应变曲线有一定困难，因此一般都采用幂指函数来替代。 , E3 J/ y0 R9 N$ v0 t' ?
( p3 Y: h/ K2 E, y8 n& }3 J
式(2)、式(3)和式(4)给出了几个根据σs、σb 和n 值进行材料硬化曲线估算的公式。首先，材料硬化曲线假定为幂指函数形式： * B8 }3 a/ K, g
/ K6 Z+ k$ ~/ W
0 j; w% h$ r* D9 d9 B7 F, v
3 L* C! v8 A% K5 L$ s1 x6 x根据单向拉伸试验σb 的定义和分散性失稳理论可得到： * ]3 j+ S6 x- e1 z+ v! M; r1 w
; m3 P1 @* k3 `9 H
6 S/ D6 m. Q6 ^( k
其中：e = 2.718。
! c" f+ p1 q2 M' U' M
& e2 V/ ?7 c: D+ A3 i+ k通常情况下，材料的实际强度要大于标准给定值30～50MPa，而对于σs 大于345MPa 的普通高强度钢板，一般硬化指数n ≈ 70/σs[5]。 |
发表于 2011-9-12 10:39:18
楼主