1、湖南大学硕士学位论文大型薄板局部特征成形过程中的凸包开裂与翘曲回弹控制技术研究姓名：毛丽臣申请学位级别：硕士专业：车辆工程指导教师：李光耀；汪晨20120524大型薄板局部特征成形过程中的凸包开裂与翘曲回弹控制技术研究摘要大型薄板局部特征类冲压件作为底板、背板、基体类零件，在汽车工业和家用电器业中应用十分普遍。这类零件的特点是板料厚度薄( O 6 o 8 m r n ) 、零件面积大、整体刚度较差，且由分布不均、深浅不一的局部特征所构成。对于较复杂的大型薄板类冲压件，由于不均匀的局部特征成形，使其在生产制造过程中经常会出现局部特征开裂、打鼓、翘曲回弹等问题，而开裂和翘曲回弹尤为严重，影响了产品
2、质量和装配精度，延长了产品的更新换代周期。对于开裂问题，从上世纪9 0 年代起，就有大量国内外学者进行了相关研究，一般通过调节压边力和设置拉延筋就能得到很好的解决。但是对于一些特殊的零件，比如背板零件，成形过程中既不能调节压边力又不能设置拉延筋。对于这种零件的开裂问题，必须寻求新的解决方法。对于回弹问题，近年来国内外学者在回弹预测与回弹控制方面做了大量的研究工作，但是对于由于局部特征成形导致的翘曲回弹问题研究的比较少。对于翘曲回弹的控制，目前也缺乏实际可行的一般方法和流程。本文基于工程实际，利用有限元数值模拟方法对大型薄板局部特征成形类零件成形过程中出现的凸包开裂和翘曲回弹控制技术进行了系统的
3、研究。本文研究的具体内容包括：( I ) 对背板零件本身的结构成形特点和成形后的应力进行了分析，并对背板零件出现翘曲回弹问题的原因和机理进行了阐述。提出背板零件成形脱模后，交替分布的拉压应力的释放以及残余弯矩使零件未成形区域产生失稳现象是导致其出现翘曲问题的根本原因。( 2 ) 提出了多步成形危险截面转移法来解决背板零件出现的矩形凸包底部开裂问题。为了使局部特征有更好的成形质量，建立了材料S E C C N 5 的质量函数曲线，并利用自适应响应面法对中间构型的几何参数凹模入口圆角R l 、凸模圆角R 2 以及成形深度H 进行了优化。( 3 ) 研究了基于局部特征的翘曲回弹控制技术。对背板零件进
4、行了特征提取，并对提取的局部特征进行了分析。提出了“弯矩抵消”方法和“波动剪切”方法，有效地控制了单特征成形过程中的翘曲问题，并利用多特征进行了验证，结果表明，这两种方法对于多特征成形过程中的翘曲回弹控制依然有效，从而可以推广到背板零件的翘曲回弹控制中。( 4 ) 研究了基于整体刚度的翘曲回弹控制技术。提出将形貌优化方法用于背板零件强化压痕的布置。对一个典型的背板零件在其未变形区域利用形貌优化方法布置了强化压痕，并对形貌优化结果进行了工程可制造化处理。计算结果表塞型璧丝堡二一明，利用形貌优化方法布置强化压痕后，背板零件的各阶非刚体模态固有频率都有较大程度的提高，背板零件的整体刚度明显增强，从而
5、可以对背板零件的翘曲起到抑制作用。关键词：模组背板；局部特征成形；翘曲回弹；多步成形：形貌优化A b s t r a c tL a r g es c a l e s h e e tm e t a ls t a m p i n gp a r t s ，s u c ha sb a s e p l a t e ，b a c kp a n e l ，m a t r i xp a r t s ，a r ev e r yw i d e l yu s e di nt h ea u t o m o t i v ei n d u s t r ya n dh o u s e h o l de l e c t r
6、i c a la p p l i a n c ei 1 1 d u s t r i e s T h ec h a l 麓c t e r i S t i co ft h e s ep a r t si st h i nt h i c k n e s s ( 0 6 O 8 r a m ) ，l a r g es i z ea n dD o o rs t i f f n e s s A n dt h e ya r ec o m p o s e do ft h el o c a lf e a t u r e sw h i c ha r en o n u n i f o r md i s t r i
7、b u t i o na n dd i f f e r e n td e p t h s D u et ot h eu n e v e nf o r m i n go ft h el o c a lf e a t u r e s ，t h ec o m p l e xl a r g es h e e tm e t a ls t a m p i n gp a r t so f t e no c c u rc o n v e xc r a c k i n g ，l i g h tw r i n k l i n go fl a r g ea r e aa n do v e r a l lw a r
8、p i n gs p r i n g b a c k T h e s ed e f e c t ss e r i o u s l ya f f e c tt h ep r o d u c tq u a l i t ya n dt h ea c c u r a c yo ft h ea s s e m b l y , a n de x t e n dt h eu p g r a d ec y c l eo ft h ep r o d u c t S i n c et h e19 9 0 s ，al a r g en u m b e ro fr e s e a r c h e r sd i dm
9、u c hr e s e a r c ha b o u tt h ec r a c l ( i n gp r o b l e m G e n e r a l l y , c r a c k i n gp r o b l e mc a nb es o l v e dw e l lb ya d j u s t i n gt h eB H Fa n ds e t t i n gd r a w b e a d s B u tf o rs o m es p e c i a lp a r t s ，s u c ha sb a c kp a n e l ，t h eB H Fc a nn o tb ea d
10、 j u s t e dd u r i n gt h es t a m p i n gp r o c e s s F o rt h ec r a c k i n gp r o b l e mo ft h i sk i n do fp a r t s ，n e ws o l u t i o n sn e e dt ob ef o u n d A b o u tt h es p r i n g b a c kp r e d i c t i o na n dc o n t r o lD r o b l e m ，m a n yr e s e a r c h e r sd i dm u c hw o
11、r k si nt h er e c e n ty e a r s B u tt h es t u d yw o r ka b o u tt h ew a r p i n gs p r i n g b a c kc a u s e db yl o c a lf e a t u r e sf o r m i n gi sr e l a t i v e l yf e w A n da tp r e s e n t ，i ti sv e r yl a c ko fp r a c t i c a lm e t h o d sa n dp r o c e s s e sf o rt h ew a r p
12、 i n gs p r i n g b a c kc o n t r 0 1 B a s e do ne n g i n e e r i n gp r a c t i c e ，t h ef i n i t ee l e m e n tn u m e r i c a ls i m u l a t i o nm e t h o di su s e di nt h i sp a p e rt of i n dt h es o l u t i o n st Ot h ec o n v e xc r a c k i n ga n dw a r p i n gs p r i n g b a c kp r
13、 o b l e mc a u s e db yl o c a lf e a t u r e sf o r m i n gd u r i n gt h el a r g es h e e tp a r t ss t a m p l n gp r o c e s s T h em a i nc o n t e n to ft h i sp a p e ri sa sb e l o w ：( 11B a s e do nt h ea n a l y s i so ft h es t r u c t u r ea n df o r m i n gc h a r a c t e r i s t i c
14、 so ft h eb a c kp a n e lp a r t s ，t h er e a s o n sa n dm e c h a n i s m so ft h ew a r p i n gs p r i n g b a c kp r o b l e ma r ee x p l a i n e d T h i sp a p e rc o n s i d e r st h a tt h eu n s t a b i l i t yc a u s e db yt h er e l e a s i n go ft h ea l t e n l a t i v ed i s t r i b
15、u t i o no ft h et e n s i l ea n dc o m p r e s s i v es t r e s sa n dt h er e s i d u a lm o m e n ti st h er o o tr e a s o no ft h ew a r p i n gs p r i n g b a c kp r o b l e m ( 2 ) T h em u l t i s t e pf o r m i n gm e t h o di sp r o p o s e dt os o l v et h eb o t t o mc r a c k i n gp r
16、o b l e mo ft h er e c t a n g u l a rc o n v e xd u r i n gt h eb a c kp a n e ls t a m p i n gp r o c e s s I no r d e rt og e tb e t t e rf o r m i n gq u a l i t y , t h eq u a l i t yf u n c t i o no ft h eS E C C N 5i se s t a b l i s h e da n dt h ea d a p t i v er e s p o n s es u r f a c em
17、 e t h o di su s e dt oo p t i m i z et h eg e o m e t r i cp a r a m e t e r s 【d i eV工程硕士学位论文e n t r a n c ef i l l e tr a d i u sR 1 ，p u n c hf i l l e tr a d i u sR 2 ，f o r m i n gd e p t hH ) o ft h ei n t e r m e d i a t ec o n f i g u r a t i o n ( 3 ) T h ew a r p i n gs p r i n g b a c kc
18、 o n t r o lm e t h o db a s e do nt h el o c a lf e a t u r e si si n v e s t i g a t e di nt h i sp a p e r T h em e t h o do ff e a t u r ee x t r a c t i o ni su s e dt oe x t r a c tas m a l lLt y p ef e a t u r ef r o mt h eb a c kp a n e lp a r t s T h e “m o m e m c o u n t e r b a l a n c e
19、 d m e t h o da n d“w a v e c u t t i n g m e t h o da r ep r o p o s e dt os o l v et h ew a r p i n gs p r i n g b a c kp r o b l e ma n dw o r kv e r yw e l li nt h ec o n t r o lo ft h ew a r p i n gs p r i n g b a c kp r o b l e mw h i c ho c c u r si nt h es i n g l ef e a t u r ea n dm u l t
20、i p l ef e a t u r e sf o r m i n gp r o c e s s ，t h u sc a nb eu s e dt oc o n t r o lt h ew a r p i n gs p r i n g b a c kp r o b l e md u r i n gt h eb a c kp a n e ls t a m p i n gp r o c e s s ( 4 ) T h ew a r p i n gs p r i n g b a c kc o n t r o lm e t h o db a s e do nt h ei n t e g r a ls
21、t i f f n e s si ss t u d i e di nt h i sp a p e r T h et o p o g r a p h yo p t i m i z a t i o nm e t h o di su s e dt od e s i g nt h el a y o u to ft h es t r e n g t h e n e di m p r i n t T h es t r e n g t h e n e di m p r i n tl a y o u ti sd e s i g n e di nt h eu n d e f o r m e dr e g i o
22、 no fat y p i c a lb a c kp a n e lp a r t sb yt h et o p o g r a p h yo p t i m i z a t i o nm e t h o da n dt h ee n g i n e e r i n gm a n u f a c t u r i n gt r e a t m e n ti sa p p l i e do nt h eo p t i m i z a t i o nr e s u l t T h ec o m p u t a t i o nr e s u l ts h o w st h a te v e r yn
23、 o n - r i g i db o a yo r d e rm o d a ln a t u r a lf r e q u e n c yh a sag r e a ti m p r o v e m e n t ，a n dt h ei n t e g r a ls t i f f n e s so ft h eb a c kp a n e li ss i g n i f i c a n t l ye n h a n c e d T h u s ，t h i sm e t h o dC a nb eu s e dt oc o n t r o lt h ew a r p i n gs p r
24、 i n g b a c kp r o b l e m K e yW o r d s ：m o d u l eb a c kp a n e l ；l o c a lf e a t u r ef o r m i n g ；w a r p i n gs p r i n g b a e k ；m u l t i s t e pf o r m i n g ；t o p o g r a p h yo p t i m i z a t i o nV：。：三呈兰塑兰兰垒兰丝童三：。：。：。：一! = ! = ! E ! ! ! ! ! ! ! ! ! = ! = ! = ! = = = ! = ! ! ! !
25、 = ! = ! ! ! = ! = = = = = = = = = = 2 5 一。插图索引图1 1 数值模拟仿真技术在冲压模具设计应用中的流程2图1 2 弹性恢复示意图5图1 3 角度回弹和卷曲回弹示意图6图1 4 扭曲回弹示意图7图1 5 翘曲回弹示意图7图2 1 模组背板C A D 模型一13图2 2 模组背板打凸模具结构图1 4图2 3 成形后的模组背板零件1 4图2 4 实际生产模具1 4图2 5 打凸工序的零件C A D 模型15图2 6 凸包成形深度示意图1 5图2 7 模组背板冲压工序流程图15图2 8 打凸过程16图2 9 上折弯过程16图2 1 0 凸包开裂的主要形式1
26、6图2 ，11 模组背板打鼓1 7图2 1 2 模组背板翘曲1 7图2 1 3 材料S E C C N 5 的K r u p k o w s k y 硬化曲线1 8图2 1 4 模组背板成形后X 方向应力分布一2 4图2 1 5 模组背板成形后Y 方向应力分布一2 4图2 1 6 背板零件脱模后的状态2 4图2 1 7 背板零件脱模后X 方向残余应力状态分布一2 5图2 18 背板零件脱模后Y 方向残余应力状态分布一2 5图2 1 9 某一凸包及周围X 方向应力分布一2 5图2 2 0 某一凸包及周围Y 方向应力分布2 5图2 2 1 模组背板成形后X X 弯矩分布2 6图2 2 2 模组背板
27、成形后Y - Y 弯矩分布2 6图2 2 3 某一局部特征x X 弯矩分布2 6图2 2 4 某一局部特征Y - Y 弯矩分布2 6图3 1 背板零件开裂的主要形式2 8图3 2 矩形局部特征C A D 模型2 9X大型薄板局部特征成形过程中的凸包开裂与翘曲回弹控制技术研究图3 3P A M S T A M P2 G 仿真分析模型2 9图3 4 矩形局部特征一步成形后F L D 图3 0图3 5 两步成形示意图及相关尺寸3 0图3 6 第一步构型C A D 模型3 0图3 7 第二步构型C A D 模型31图3 8 第一步构型仿真模型31图3 9 第二步构型仿真模型31图3 1 0 中间构型F
28、 L D 图3 2图3 1 1 最终构型F L D 图3 2图3 1 2 一步成形和两步成形二维F L D 图比较3 2图3 1 3 中间构型几何参数组合一3 3图3 1 4 中间构型几何参数组合二3 3图3 15r a d i o s s o n e s t e p 一步法仿真模型3 4图3 1 6P A M S T A M P2 G 增量法仿真模型3 4图3 17r a d i o s s o n e s t e p 一步法仿真结果( 厚度减薄率) 3 5图3 18P A M S T A M P2 G 增量法仿真结果( 厚度减薄率) 3 5图3 1 9 材料S C E E N 5 质量函数
29、曲线3 6图3 2 0 自适应响应面方法迭代示意图3 8图3 2 1 目标函数F ( x ) 的迭代历程3 9图3 2 2 约束条件G 1 ( x ) 的迭代历程3 9图3 2 3 设计变量圆角R 1 的迭代历程3 9图3 2 4 设计变量圆角R 2 的迭代历程4 0图3 2 5 设计变量成形深度H 的迭代历程4 0图3 2 6 第一步成形优化前与优化后成形极限图对比4 l图3 2 7 第二步成形优化前与优化后成形极限图对比4 1图3 2 8 最终成形后零件优化前与优化后的对比4 1图4 1L 型特征C A D 模型4 4图4 2L 型特征P A M S T A M P2 G 仿真分析模型4
30、4图4 3L 型特征成形后X 方向应力分布4 5图4 4L 型特征成形后Y 方向应力分布4 5图4 5 单个L 型特征回弹分析结果4 5图4 6L 型特征回弹结束之后的状态4 6图4 7 反对称L 型特征C A D 模型4 6图4 8 反对称特征P A M S T A M P2 G 仿真分析模型4 7图4 9 反对称特征成形后X 方向应力分布4 7X：! ：至堡至型呈兰篁兰鎏圣：，：。：。一图4 1 0 反对称特征成形后Y 方向应力分布4 7图4 11 反对称特征回弹分析结果4 8图4 1 2 反对称特征回弹结束之后的状态4 8图4 1 3 对称特征C A D 模型4 9图4 1 4 对称特征
31、P A M S T A M P2 G 仿真分析模型4 9图4 1 5 对称特征成形后X 方向应力分布状态4 9图4 1 6 对称特征成形后X 方向应力分布状态5 0图4 1 7 对称特征回弹分析结果5 0图4 1 8 对称特征回弹结束之后的状态5 0图4 1 9 单个L 型特征成形之后X X 弯矩云图一5 1图4 2 0 单个L 型特征成形之后Y - Y 弯矩云图5 2图4 2 1 “弯矩抵消”方法分析模型5 2图4 2 2 特征筋与局部特征下凹方向示意图5 2图4 2 3 特征筋与局部特征同向C A D 模型5 3图4 2 4 特征筋与局部特征同向P A M S T A M P2 G 分析模
32、型5 3图4 2 5 “弯矩抵消”正向模型成形之后X - X 弯矩云图5 4图4 2 6 “弯矩抵消 正向模型成形之后Y - Y 弯矩云图5 4图4 2 7 “弯矩抵消”正向模型回弹分析结果5 4图4 2 8 原始模型与“弯矩抵消”正向模型脱模之后状态的比较5 5图4 2 9 特征筋与局部特征反向C A D 模型5 5图4 3 0 特征筋与局部特征反向P A M S T A M P2 G 分析模型5 6图4 3 1 “弯矩抵消”反向模型成形之后X - X 弯矩云图5 6图4 3 2 “弯矩抵消”反向模型成形之后Y Y 弯矩云图5 6图4 3 3 “弯矩抵消”反向模型回弹分析结果5 7图4 3
33、4 原始模型与“弯矩抵消”反向模型脱模之后状态比较5 7图4 3 5 “波动剪切”方法假想分析模型5 9图4 3 6 特征筋与局部特征下凹方向示意图5 9图4 3 7 特征筋与局部特征同向C A D 模型5 9图4 3 8 特征筋与局部特征同向P A M S T A M P2 G 分析模型6 0图4 3 9 原始模型与“波动剪切”正向模型X 方向应力分布比较6 0图4 4 0 原始模型与“波动剪切”正向模型Y 方向应力分布比较6 0图4 4 1 “波动剪切”正向模型回弹分析结果6 1图4 4 2 原始模型与“波动剪切”正向模型脱模之后状态比较6 1图4 4 3 特征筋与局部特征反向C A D
34、模型6 2图4 4 4 特征筋与局部特征反向P A M S T A M P2 G 分析模型6 2X I大型薄板局部特征成形过程中的凸包开裂与翘曲I 旦I 弹控制技术研究图4 4 5 原始模型与“波动剪切”反向模型x 方向应力分布比较6 2图4 4 6 原始模型与“波动剪切”反向模型Y 方向应力分布比较6 3图4 4 7 “波动剪切”反向模型回弹分析结果6 3图4 4 8 原始模型与“波动剪切”正向模型脱模之后状态比较6 3图4 4 9 各模型最大回弹量6 4图4 5 0 验证算例一的四种模型的C A D 模型6 5图4 5 1 原始模型回弹分析结果6 6图4 5 2 “弯矩抵消”模型回弹分析结
35、果6 6图4 5 3 “波动剪切”模型回弹分析结果6 6图4 5 4 混合模型回弹分析结果6 7图4 5 5 各模型最大回弹量6 7图4 5 6 验证算例二的各种模型的C A D 模型6 8图4 5 7 原始模型回弹分析结果6 9图4 5 8 “弯矩抵消”模型回弹分析结果6 9图4 5 9 “波动剪切”模型回弹分析结果6 9图4 6 0 混合模型回弹分析结果7 0图4 6 1 各模型最大回弹量7 0图5 1 背板零件简化C A D 模型7 3图5 2 背板零件模态分析有限元网格7 3图5 3 各阶模态固有频率计算结果7 4图5 4 第一阶非刚体模态振型7 4图5 5 第二阶非刚体模态振型7 4
36、图5 6 第三阶非刚体模态振型7 5图5 7 第四阶非刚体模态振型7 5图5 8 第五阶非刚体模态振型7 5图5 9 第六阶非刚体模态振型7 5图5 1 0 带强化压痕的实际背板零件一7 6图5 11 带强化压痕的实际背板零件二7 6图5 1 2 形貌优化的设计流程7 7图5 1 3 设计区域与非设计区域示意图7 8图5 1 4 设计区域与非设计区域有限元网格7 8图5 1 5 背板零件局部特征成形脱模后的状态7 8图5 16 第2 阶非刚体模态振型7 9图5 17 第3 阶非刚体模态振型7 9图5 18 背板零件形貌优化加筋对称平面8 0工程硕士学位论文图5 1 9 强化压痕参数示意图8 0
37、图5 2 0 优化各阶模态固有频率计算结果8 l图5 2 l 优化二各阶模态固有频率计算结果8 l图5 2 2 优化三各阶模态固有频率计算结果8 2图5 2 3 强化压痕位置分布优化结果8 2图5 2 4 优化四各阶模态固有频率计算结果8 2图5 2 5 强化压痕位置分布优化结果8 3图5 2 6 背板零件形貌优化后强化压痕分布8 4图5 2 7 背板零件形貌优化后等值线云图8 4图5 2 8 强化压痕工程可制造化处理后的结构8 4图5 2 9 重构后结构各阶模态固有频率计算结果8 5X大型薄板局部特征成形过程中的凸包开裂与翘曲回弹控制技术研究附表索引表1 1 破裂产生的原因以及解决方法4表2
38、 1 材料S E C C N 5 成分参数18表2 2S E C C N 5 材料参数一18表3 1 1 0 个测量点厚度减薄率对比与相对误差3 5表3 2 设计变量的初始值与优化值( 单位：m m ) 3 8表4 1 原始模型、“弯矩抵消”模型和“波动剪切”模型比较6 4表4 2 “弯矩抵消”模型、“波动剪切”模型、混合模型与原始模型比较6 7表4 3 “弯矩抵消”模型、“波动剪切”模型、混合模型与原始模型比较一7 0表5 1 背板零件无约束非刚体模态频率与振型描述7 4表5 2 优化过程描述7 9表5 3 原始模型与最终优化方案对比8 3表5 4 原始模型与工程可制造化后模型对比8 5工程
39、硕士学位论文第1 章绪论1 1 引言近几年来，在“汽车下乡”、“家电下乡”和“以旧换新 等政策的刺激下，汽车以及家电类产品需求不断扩大，销售量急剧上升。据相关统计资料显示，我国已经成为全球最大的汽车产销国，2 0 1 0 年中国汽车产量已超过日本和美国，位居世界第一，达到1 8 0 0 万辆。汽车工业已经成为我国经济的重要支柱产业之一，它有力促进了化工、材料、制造、电子等多个产业的发展。薄板冲压成形作为一种极其重要的金属塑性成形方法，在汽车、航空、电器和国防等工业中得到了广泛的应用。在最近十几年，由于汽车和家用电器等工业的崛起，薄板成形技术得到了空前的发展。薄板冲压成形技术得以持续的快速发展，
40、主要归功于它的许多优点，如材料利用率高，产品尺寸精度稳定，生产操作简单，容易实现自动化；大规模生产加工成本低、生产效率高等特点。薄板类零件通过优化设计可达到高强度、高刚度和质量轻的要求。有人说汽车技术的竞争就是薄板冲压成形技术的竞争，这从一定意义上讲是有道理的，因为汽车车身的设计和制造是影响整个汽车产品效益的最关键因素，而汽车车身零件绝大部分是薄板冲压件。模组背板是汽车电器和家用电器上的一个重要零件，其生产过程是一种典型的薄板冲压成形过程，其制造难度和产品要求非常高。模组背板平面度要求很高，不仅要求性能的波动范围小，而且不同冲压方式( 连续冲压和单冲) 、不同模具、不同冲压道次、不同变形量等对
41、钢板的机械性能要求都不一样。破裂、起皱、和回弹是薄板冲压成形过程中最常见的三大问题。在模组背板零件的冲压成形过程中，同样存在这三大问题，分别对应：局部特征开裂、打鼓和翘曲回弹。这三大缺陷严重影响了背板零件产品质量和装配精度，延长了产品更新换代周期，所以研究背板等大型薄板局部特征成形类零件的凸包开裂、打鼓和翘曲回弹具有重要的工程意义和实用价值。1 2 数值模拟技术在薄板冲压成形中的应用薄板冲压成形技术有着悠久的历史，在汽车工业、家电行业等领域应用非常广泛。在冲压产品的整个设计制造过程中，模具的设计制造水平是制约产品开发周期和产品品质的关键因素。传统的工艺制定和模具设计大都依赖于工程技术人员的经验
42、，通过不断反复的试模、修模、调模来保证零件的质量。这样的方法用大型薄板局部特征成形过程中的凸包开裂与翘曲回弹控制技术研究于产品设计，不仅周期长、费用高，产品品质往往难以得到保证。随着计算机技术、数值计算方法等相关学科的迅速发展，冲压成形计算机模拟仿真技术不断成熟，并在冲压工艺与模具设计中发挥着重要的作用。运用冲压成形过程模拟仿真技术可以帮助工程师解决传统方法难以解决的冲压工艺和模具设计难题。通过给定模具和相应的工艺方案，定义材料的变形过程，就可以预测材料的流动趋势以及成形缺陷( 起皱、拉裂、回弹) ，从而可以判断模具和工艺方案的合理性，模拟结果可以用于指导模具的修改以及冲压工艺的优化。每一次模
43、拟仿真计算过程就相当于一次试模过程( 如图1 1 所示) 。因此成熟的模拟仿真技术不仅可以降低试模次数，缩短产品开发周期，降低成本，还可以提高产品质量以及市场竞争力l J J 。图1 1 数值模拟仿真技术在冲压模具设计应用中的流程2 0 世纪8 0 年代后期，随着计算机技术飞速发展和有限元方法的曰益成熟，在世界汽车工业应用需求的推动下，板料成形过程的计算机模拟迎来了蓬勃发展的时期。为了促进板料成形模拟技术的研究和应用，国际上发起了定期( 每三年一届) 召开的国际板料成形三维数值模拟会议N U M I g H E E T ( I n t e r n a t i o n a lC o n f e
44、r e n c eo n N u m e r i c a lS i m u l m i o no f 3 DS h e e tF o r m i n g ) ，从1 9 9 3 年开始至2 0 1 1年，已经连续召开8 届。N U M I S H E E T 上论文的发表以及标准考题的提出及应答，有力地推动了板料成形模拟技术的研究【2 】。随着冲压成形仿真技术的迅速发展，以及对有限元方法模拟冲压过程认识的不断深入，人们对车身零部件如行李箱盖板、油底壳、车门外板、横梁等的成形工程硕士学位论文过程成功的进行了模拟，仿真结果被用于指导模具和冲压工艺的设计，并且涌现出了一批优秀的商业软件。目前国际上应
45、用较广泛的专门针对板料成形的C A E软件有法国E S I 集团的P A M S T A M P ，美国E T A 和L S T C 公司联合开发的D Y N A F O R M 、A l t a i r 公司开发的H Y P E R F O R M ，瑞士A u t o F o r m 工程有限公司A U T O F O R M 以及德国G N S 公司的I N D E E D 等。国内也有不少研究机构开发了具有自主知识产权的板料成形模拟仿真系统，比如吉林大学的K M A S 、华中科技大学的F A S T A M P 、湖南大学的C A D E M 以及北京航空航天大学的S H E E T
46、 F O R M等软件。利用这些成熟的商业软件，许多学者在板料的拉裂预测与消除、起皱预测与消除、回弹计算与控制、压边力优化、成形力预测等方面做了大量的工作p J ，如D o n g 和L i n l 4 l 禾l J 用数值模拟技术分析了桑塔纳2 0 0 0 的某外板的成形性，模拟结果与实验结果吻合得相当好。M a k i n o u c h i e l 5 1 利用冲压成形仿真软件模拟了若干冲压件的成形过程，计算了它们的坯料形状、减薄量、残余应力、回弹量，并预测了成形后的一些常见缺陷。T a k u d a 6 等利用有限元分析软件研究了在温成形条件下铝合金圆柱件深拉伸的成形情况，并与实验结
47、果进行了对比，结果表明，仿真分析可以成功预测颈缩发生的位置和成形极限。L i u 7 1 等利用有限元分析软件研究了变压力对U 形件回弹的影响，并且提出了一套有效的变压边力行程方法，该方法在保证零件不出现拉裂和起皱缺陷的情况下，能最大程度的减少回弹量。1 3 薄板冲压成形中破裂问题研究现状1 3 1 破裂产生的原因及解决方法在薄板冲压成形中有四种典型的成形方式：拉深、弯曲、翻边( 包括扩孔) 、胀形。在这四种典型的成形中，出现的破裂大致可分为a 破裂、破裂和弯曲破裂这三种形式【8 J o 这些破裂产生的力学原因如下：1 1 口破裂：由于板料所受拉应力超过材料强度极限所引起的破裂。2 )破裂：由
48、于板料的伸长变形超过材料的局部延伸率引起的破裂。3 1 弯曲破裂：由于弯曲变形区的外层材料中拉应力过大引起的破裂。以上是这三种破裂问题产生的力学原因，从整个冲压过程成形过程分析，产生破裂问题的主要因素有：1 ) 材料参数与材料性能有关的参数，如厚向异性系数r 和硬化指数n 等对零件成形质量有很大的影响。表现在成形极限图( F L D ) 上为：硬化指数n 越大，破裂极限曲线越高。厚向异性系数r 值越大，应变路径与破裂极限曲线的交点值越大，越不容易出现破裂问题。当r 值 1 时，厚度方向的变形比面内困难，拉深不易变薄，从而不易出现破裂问题。大型薄板局部特征成形过程中的凸包开裂与翘曲回弹控制技术研
49、究2 ) 模具型面参数模具型面的圆角半径取值也是影响零件是否会出现破裂的重要原因之一。若圆角过小，材料流动困难，圆角处变薄剧烈，就会出现圆角开裂缺陷。另外工艺补充的几何参数取值不合理也会使零件出现开裂问题。3 ) 工艺过程参数冲压过程的工艺参数，比如压边力、拉延阻力、摩擦系数、成形速度等都对零件的成形质量影响很大，其中压边力和拉延阻力的调节是改善零件成形质量的一项重要措施。表1 1 所示为冲压成形过程中产生破裂的原因以及解决方法的总结。表1 1 破裂产生的原因以及解决方法1 3 2 破裂问题的研究现状在过去几十年里，许多学者对于薄板冲压成形过程中破裂问题的预测以及消除，做了大量研究。如K e
50、e l e r 【9 1 和G o o d w i n 0 0 1 利用网格应变分析技术提出了成形极限图( F L D ) 的概念，并指出根据应变点坐标所处区域可以预测零件的开裂。D o e g e 和G r o c h e l l l l 提出将有限元模拟与成形极限图F L D 相结合可以用于预测零件的开裂，并且可以将拉深引起的开裂和胀形引起的开裂区分开来。A h m e t o g l u t l 2 】等人利用有限元软件模拟了深盒形件铝合金的成形，并通过控制压边力和润滑状况来避免出现破裂和起皱缺陷。U t h a i s a n g s u k 【13 J 等人使用成形极限应力图( F
