铝合金厚板时效成形回弹补偿算法

  • 格式:pdf
  • 大小:479.36 KB
  • 文档页数:8

and sphere-part are reduced below 0.4 mm,after four iterations.And after five iterations,the maximal shape error of the saddle-part is reduced below 0.5 mm。except for the four comers of it.This algorithm’S advanta-
1 回弹预测
1.1蠕变本构模型
1994年,Kowalewski等‘171给出了单轴受力 条件下的蠕变损伤统一本构方程。该方程通过 引入状态变量,可以描述铝合金材料在蠕变初 期直至蠕变第3阶段的位错强化、高温下的蠕 变软化现象和晶界处空穴形核长大而引起的
第5期
黄霖等:铝合金厚板时效成形回弹补偿算法
软化。 时效成形中材料不会进入蠕变第3阶段,只
算法具有计算精度高、收敛速度快的优点。
关键词:铝合金;时效成形;有限元;回弹;模具设计
中图分类号:TB31
文献标识码:A
Abstract:The age forming technique,characterized by huge springback,has been developed tO manufacture in—
后回弹量的准确预测,提出了有限元计算驱动的 适用于时效成形的模具型面修正算法。
[11]Karafillis A P,Boyce M C.Tooling and binder design for
sh舱‘”伯¨o”…8 P眦∞58∞mP“5叶…g 5pn嘞g”k
(2)通过计算3个实例,即单曲率的圆柱面:兰::峨Int。e。rn。。at,io。n。a(1。J)o:裂。£.№曲in。Tod8“8
航空学报 设定初始模具型面S牛=so作为迭代起点, 当max(E。)≤ERRToL时停止迭代,ERRTOL 为可接受的最大型面误差。
回弹后的板料曲面
第29卷
Fig.2
图2偏差调节法流程图 Flowchart of deviation adjusting(DA)method
(b)型面误差分布
3计算实例
(a)型面误差分布
图5马鞍形零件设计型面
Fig.5 Design surface of saddle-part
万方数据
迭代次数 (”型面最大误差的收敛情况 图6 马鞍面零件计算结果 Fig.6 Saddle-part calculation results
4结论
(1)基于有限元模型对铝合金厚板时效成形
(b)掣嘶误差分御

量 划 嗤 X 堪 嚏 尉4球面零件计算结果
Fig.4 Sphere-part calculation results
3.3马鞍面零件


圆柱面能简化为二维问题,而球面具有轴对


称性。为了更好地验证补偿修正算法,针对更具
恒 副
有一般性的马鞍形曲面零件进行了分析。图5为


X/mm (a)零件与其设计型面的截面线
万方数据
(a)零件与其设计型面的截面线
第5期
黄霖等:铝合金厚板时效成形回弹补偿算法
所设计的马鞍形曲面零件,板料长780 mm,宽 580 iTlm,厚度为25 mm。迭代计算中,调整控制 系数k取0.9。
图6为模具迭代调整的计算结果。可见,经过 5次补偿修正,零件型面最大误差从12.360 mm减 小到0.712 mm,且最大误差分布于零件的4个 角部,在零件其他部位误差分布均匀,且均小于 0.500 mm,如图6(a)所示。由于零件的复杂性 增加,使收敛速度减慢,且所能达到的成形精度降 低,但是仍然控制在可接受的范围内。
tegral wing-skin panel components.The die surface should be modified,based on accurate prediction of spring。
back.tO eliminate the effect of springback on forming precision.A finite element method(FEM)drove itera—
襄1 155℃下蠕变本构方程中材料常数的值
Table 1 Values of material constants for creep constitutive equations at 155℃
材料常数 A/h一1
数值
6 ×"
B/MPa一1
仉叫
踟/MPa
3∞
h/MPa
2∞
H’
Oo

1_
m1
OJ
1.2有限元模型的建立
2补偿修模算法
2.1型面误差定义
考虑到时效成形所适用的均为大曲率半径的 曲面零件,将实际零件曲面与其设计型面进行配 准后放置在模型空间中,并保持曲面质心点法矢 与模型空间的Z轴重合。定义两曲面间Z方向 的距离为型面误差。
采用离散的三角网格模型来描述零件的设计 型面。如图1所示,从零件的有限元网格模型上 一节点P沿Z轴向设计型面投影,与三角面片 ABC相交于点P 7,则p-p-,即为点P处的误差。 遍历零件有限元模型上的每一个节点,使得到完 整的型面误差分布。
万方数据
z l
火P
图1型面误差定义
Fig.1 Definition of surface shape error
2.2模具型面调整算法
根据有限元模型预测的型面成形误差向设计 模具型面投影,据此得到模具上各点的偏移量并 对模具型面进行调整,消除回弹对零件形状误差 的影响,如图2所示。为便于描述,使用二维模型 对算法迸行说明。
第29卷第5期
2008年
9月
航空学报
ACTA AERONAUTICA ET ASTRONAUTICA SINICA
V01.29 No.5 Sept. 2008
文章编号:1000一6893(2008)05—1406—05
铝合金厚板时效成形回弹补偿算法
黄霖,万敏
(北京航空航天大学飞行器制造工程系,北京 100191) Compensation Algorithm for Springback in Age Forming for Aluminum Alloy Thick Plate
ges of high precision and fast convergence rate are confirmed.
Key words:aluminum alloy;age forming;FEM;springback;die design
回弹是板料成形中造成形状误差的主要因 素,使成形后的零件与其设计形状之间产生偏差, 并将影响后续的装配。时效成形中的回弹是由应 力松弛后的残余应力引起的。由于时效成形是在 小于材料屈服强度的应力状态下进行的,应力松 弛的温度与时间受到时效工艺的限制,所以成形 后的大回弹量成为该工艺的一大特点和需要解决 的关键技术问题。
传统的控制回弹的方法多基于回弹机理,通 过增加板料成形过程的拉伸量减小回弹对零件形 状的影响阻21。随着有限元技术及其他数值方法 在板料成形回弹仿真及预测中的广泛应用[3{],并 发展出基于回弹精确预测的非解析回弹补偿算 法。目前,主要的非解析回弹补偿算法为:回弹向 前法(K&B)[乳1川,响应面法(RSM)[12143和偏差调
图3为模具迭代调整的计算结果。可见,随 着修模次数的增加,零件与其设计型面误差逐步 减少;经过4次补偿修模后,零件型面最大误差从 6.778 0 mm减小到0.315 7 mm。
3.2球面零件
零件设计型面为半径是6 000 mm的球面, 板料尺寸760 mm×760 mm,厚度是35 mm。迭
迭代次数 (c)型面最人误差的收敛情况
(1)根据经过i次修正的模具型面S午,分析 铝合金厚板的时效成形,并预测卸载后的回弹。
(2)计算回弹后零件曲面与零件设计曲面 SD之间的误差,即计算各个对应点之间的偏移量
Ef。
(3)调整定义模具型面的节点,生成新的模 具型面S铲1一S}+kEi,返回第(1)步进行下一次 迭代计算。其中五为调整控制系数,取值一般介 于0.8至1.2之间,可根据回弹率来确定,回弹率 高志取较大值,反之则取较小值。
编写ABAQUS子程序CREEP将材料的蠕 变本构模型应用到有限元分析中,通过数值方法 模拟铝合金厚板时效成形及回弹。
根据时效成形工艺过程,建立如下分析步骤: (1)在板料背面施加足够的载荷,使其发生 变形并与模具表面完全紧密接触。 (2)保持载荷及板料的变形量相当一段时 间,在恒定温度下,使材料发生应力松弛,部分弹 性变形转变为蠕变应变,材料内部应力下降。 (3)卸载,板料回弹,测量回弹量。 对有限元模型预测的回弹量进行了实验验 证,验证了时效成形有限元模型的准确性口…。
Huang Lin.Wan Min
(Department of Aircraft Manufacturing Engineering,Beijing University of Aeronautics and Astronautics,Beijing 100191,China)
摘要:时效成形是一种适用于飞机整体壁板零件制造的成形工艺,具有回弹量大的特点,因此需要在准确
图3圆柱面零件计算结果
Fig.3 Arc-part caculation results
代计算中,调整控制系数k取0.8。 图4为模具迭代调整的计算结果。可见,对
于球形零件,第1次修模就表现出显著的效果,零 件与其设计型面误差显著减少;经过4次回弹补 偿修正,零件型面最大误差从18.310 0 mm减小 到0.390 4 mm。
需描述材料在蠕变第1、第2阶段的特性,因此不 需要考虑该本构方程中的损伤变量。通过简化和 必要的修正,提出如下蠕变本构方程n引: