铝合金预拉伸板残余应力数值模拟方法研究_陈树梁.pdf
《铝合金预拉伸板残余应力数值模拟方法研究_陈树梁.pdf》由会员分享,可在线阅读,更多相关《铝合金预拉伸板残余应力数值模拟方法研究_陈树梁.pdf(6页珍藏版)》请在咨信网上搜索。
1、基金项目:四川省重点研发计划(2021YFG0060)收稿日期:20210901修回日期:20210912第 40 卷第 2 期计算机仿真2023 年 2 月文章编号:10069348(2023)02029205铝合金预拉伸板残余应力数值模拟方法研究陈树梁,徐雷*,胡元昊,高涵(四川大学机械工程学院,四川 成都 610065)摘要:针对铝合金预拉伸板残余应力仿真结果误差较大的问题,提出一种考虑铝合金材料参数变化的数值模拟方法。铝合金预拉伸板淬火和拉伸过程中屈服强度和弹性模量的变化考虑进来,其中屈服强度的变化通过设置温度阈值实现,弹性模量的变化通过将每个单元积分点弹性模量表示成温度和塑性变形的函
2、数来实现。通过验证发现,上述方法有效地降低了仿真结果与实验数据之间的误差,能很好地预测铝合金板残余应力的大小和变化趋势。通过对照仿真发现淬火与拉伸工艺之间的自然时效时间越长,所提方法的效果越显著。关键词:铝合金预拉伸板;残余应力;数值模拟;自然时效中图分类号:TP391.9文献标识码:BStudy on Numerical Simulation Method of esidualStress of Aluminum Alloy PreStretching Board Based onChange of Material ParametersCHEN Shuliang,XU Lei*,HU Yu
3、anhao,GAO Han(College of Mechanical Engineering,Sichuan University,Chengdu Sichuan 610065,China)ABSTACT:Aiming at the problem of large error of in the simulation results of residual stress of simulation resultsof aluminum alloy prestretching board,a numerical simulation method considering the change
4、 of aluminum alloy ma-terial parameters is proposed When establishing the simulation model,the changes of yield strength and elastic mod-ulus during quenching and stretching of aluminum alloy prestretching board were are taken into account The changeof yield strength wais realized by setting the tem
5、perature threshold,and the change of elastic modulus wais realized byexpressing the elastic modulus of each element integral point as a function of temperature and plastic deformationThrough The verification shows,it is found that this method can effectively reduce the error between the simulationre
6、sults and the experimental data,thusand can well predict the magnitude and change trend of residual stress in alu-minum alloy plate Throughrough the comparative simulation experiment,it is found,it is found that the longer thenatural aging time between quenching and drawing process,the more signific
7、ant the effect of this method isKEYWODS:Aluminum alloy prestretching board;esidual stress;Numerical simulation;Natural aging1引言铝合金尤其是 7xxx 铝合金因其密度小、强度高、硬度高等优良性能,广泛应用在航空航天产品的薄壁结构件中。其刚度低、壁薄、弹性模量小的特点使得残余应力导致的加工变形成为了影响铝合金厚板加工精度的重要因素。而且为了增加薄壁构件的强度及其轻量化,薄壁构件大型化和整体化的趋势也使得加工变形问题更为突出。残余应力主要包括初始残余应力和加工产生残余应力,
8、随着高速铣削工艺的发展和加工参数的优化,加工产生的残余应力越来越小,因此对铝合金预拉伸板初始残余应力的精确数值模拟对薄壁结构件变形的预测和加工工艺调整尤为重要。针对铝合金厚板初始残余应力的研究,曹海龙等研究了7055 铝合金淬火条件对残余应力分布规律的影响1。张园园等对 7075 铝合金做了大量研究,包括淬火的换热系数反求,进行淬火以及预拉伸仿真得到铝合金厚板的初始残余应力分布规律以及残余应力与拉升量和拉伸速率的关系等2。李雪春等认为,一般工程应用中都把弹性模量作为常数,实292际上金属的弹性模量在塑性变形过程中,随着塑性变形程度的增加会发生变化。另外李春雪等研究了铝合金的弹性模量值随塑性变形
9、程度及加载方式变化的规律3。刘瑶琼等,林峰等分析了塑性变形对铝合金弹性模量的影响,通过对模型分层建立了考虑弹性模量差异的 7050 铝合金厚板“淬火预拉伸”有限元模型,分析了弹性模量差异对 7050 铝合金厚板残余应力的影响4,5。实际上,铝合金预拉伸板在淬火与预拉伸整个过程中,材料的一些参数是变化的,其中对预拉伸板初始残余应力影响较大的主要是弹性模量与屈服极限。现有的仿真很少见到考虑屈服极限变化的情况,而考虑弹性模量差异的仿真基本都是采用对模型分层的方法,这种方法将材料进行了离散化处理,虽然减少了仿真的误差,但同样导致仿真结果不连续,与实际不够吻合。针对铝合金材料淬火预拉伸过程中参数变化这个
10、问题,本文以 7075 铝合金预制板为对象进行仿真,仿真模型不仅考虑了铝合金屈服极限的变化,而且将仿真模型的每个积分点的材料弹性模量与其温度、等效塑性变形联系起来,使得仿真结果更加精确。2仿真原理及参数确定铝合金预拉伸厚板的工艺路线如图 1 所示图 1铝合金预拉伸厚板的工艺路线图固溶的时候铝合金厚板内部应力几乎为零,因此铝合金厚板的初始残余应力主要跟淬火、预拉伸和时效过程相关。因此本文主要探究淬火和预拉伸过程。3淬火过程数学模型零件的淬火过程是应该热交换过程,根据能量守恒定律和 Fourier 传热定律,零件的瞬态温度场应满足传热分析的过程控制方程,即xkxtx+ykyty+zkztz+Q c
11、Tt=0(1)在式(1)中,kx、ky、kz分别为 X、Y、Z 三个方向的热传导系数。为密度,Q 为材料内部相变潜热,c 为比热,t 为时间。铝合金板从固溶到淬火转移时间很短,而且淬火温度远低于相变温度,因此可以忽略淬火过程的相变,即可令 Q=0。另外淬火过程工件的热辐射传热也远小于工件与介质之间的对流换热,因此传热过程的边界条件可以视为第三类边界条件,即 kTn=H(TW TC)(2)在式(2)中,k 为边界方向的热传导率;H 为边界的换热系数;n 为边界外法线方向向量;Tw,Tc分别为工件表面温度与介质温度2。3.17075 铝合金热力学参数7075 铝合金材料的热力学参数如密度、弹性模量
12、、导热系数、比热容以及热膨胀系数都与温度 T 相关。具体参数如表 1 所示 68。表 1铝合金 7075 的热力学基本参数表温度密度kg/m3弹性模量 GPa热导率W/(m)线膨胀系数106/比热J/(kg)泊松比2028007100015521.685010027756519316123.490020027505626217524.39700.3330027253798218525.2102040027003150019330.7112050026752500019731.41320另外淬火过程中的换热系数也是非常重要的参数,这里采用浸没淬火的方式,换热系数见表 2。表 27075 铝合金浸
13、没淬火(水淬)换热系数表温度()换热系数温度()换热系数20130012100133509.515021.54007.520022450725016.55006.753.2应力应变计算原理2.1.1弹性应变阶段弹性阶段应力应变关系采用胡克定律进行描述,即=E(3)在式(3)中,E 为杨氏模量,一般计算都将材料的弹性模量看成常量或者随温度变化的变量进行计算,而实际上,杨氏模量 E 还是关于塑性应变和应变率的变量,由于拉伸过程为低应变率,因此忽略应变率的影响。另外,由塑性变形产生的杨氏模量变化会随着时间由快到慢逐渐恢复9,而且当多次加载和单次加载到达相同塑性变形值时,多次加载明显造成弹性模量下降更
14、多10,11,另外考虑到淬火过程产生的塑性变形比拉伸过程要小得多,因此忽略淬火阶段与预拉伸阶段之间的弹性模量恢复。查阅相关文献,得到了 7075 铝合金杨氏模量与塑性变形的函数关系式5,12。E=5k=1ak5k(GPa)(4)式中:系数 k(1k5)分别为 4057350、797878.143、55548.186、1596.196、70.952。根据表 1 中杨氏模量与温度的数据拟合得到弹性模量与温度的关系392ET=7 104T3 0.522T2 9.673T+71325(5)弹性模量微观上主要取决于原子间的结合力,而温度变化会引起原子间距变化从而影响弹性模量变化13。对于面心立方结构的铝
15、合金,塑性变形造成的弹性模量变化主要是由变形中形成的主要织构成分含量的变化引起的14。而温度也会影响织构的变化和恢复,但考虑到淬火时间短,且淬火阶段整体上产生的塑性变形只占整个过程的小部分。因此为了方便研究,假设温度与塑性变形都独立作用于杨氏模量,且塑性变形对不同温度下杨氏模量的影响是等比列的,可以将弹性模量表示成式(6)。E=E70.95ET(6)2.1.2塑性应变阶段当材料达到屈服极限进入塑性应变阶段的应力主要与应变、温度相关。实际上,铝合金固溶淬火处理之后的铝合金材料屈服强度在自然时效下会得到提高,其变化趋势15 如图 2 所示。图 2铝合金板淬火后的时效特征图根据图 2 可以看出:室温
16、情况下,在淬火后约 0.3 小时内,屈服强度几乎维持在 135MPa 不变,而在淬火 0.3 小时之后提高的速率由快速到缓慢缓慢转变,经过一个小时的自然时效,屈服强度达到 180MPa 左右,而淬火经过较长的自然时效(100h)屈服强度可以达到 300MPa 左右。但是大部分铝合金仿真研究都是使用固定屈服极限来代替整个过程进行仿真。比如采用淬火拉伸之后材料的屈服极限进行仿真,甚至采用进行人工时效后的屈服极限进行仿真,这样仿真出来的结果跟实际情况有一定的差异。因此对不同阶段采用相对应的屈服极限进行仿真是非常有必要的,为了得到不同阶段铝合金材料的屈服极限特性,查阅文献得到 7075T651 铝合金
17、的屈服极限随温度变化如表 4 所示16。7075 铝合金在屈服强度为 225MPa 时的应力应变数据如表 5 所示4。表 47075T651 铝合金屈服极限的温度特性表温度20100200300400500620屈服极限456389275.747.135.516.30表 57075 铝合金应力应变数据表塑性应变00.0060.0160.0260.034应力 MPa225260300330350因为铝合金材料淬火后不同时效时间应力应变曲线不方便测量,为了方便研究,假设不同自然时效时间的铝合金的软化系数与 7075T651 相同以及不同自然时效时间的铝合金的塑性应变的硬化系数相同17。然后根据前面
18、提到的数据可以推导出不同自然时效铝合金在不同温度下的拉伸应力应变曲线。虽然与实际参数存在误差,但是该研究主要是为了探究改进仿真的可行性和有效性,在误差允许范围内是可以接受的。4仿真模型的构建4.1基本假设影响铝合金预拉伸板材残余应力的因素非常多。为了突出主要影响因素,仿真过程做出如下假定条件:1)材料视为各向同性的弹塑性材料。2)淬火介质温度始终保持恒定。3)忽略淬火转移时间,即不考虑淬火工件专员过程中与空气的热量交换。4)忽略淬火工件进入淬火介质的过程。5)拉伸速率很小,忽略其对残余应力消除的影响。6)视不同厚度铝合金板时效效果相同。4.2模型设置使用 800mm220mm50mm 的 70
19、75 铝合金厚板加热到475固溶处理,然后在 20的水中淬火 300s,淬火后放置不同的时间(屈服极限分别强化到 180MPa,300MPa)再进行拉伸消除残余应力处理,拉伸量分别为 1.5%、2.0%、2.5%。设置四组对比仿真,第一组是考虑屈服极限和弹性模量变化的改进仿真组;第二组是不考虑弹性模量变化以及以拉伸时屈服极限作为整个仿真时间的屈服极限的仿真;第三组是只考虑屈服极限变化的仿真组;第四组是只考虑弹性模量变化的仿真组。为了简化模型,这里假设厚板各个表面同时淬火,温度场,应力场对称分布,所以使用 1/8 模型仿真以减少计算量。仿真采用 4 个分析步,分别对应淬火、拉伸、卸载和剪裁。淬火
20、过程的热力耦合计算的常见方式有顺序耦合和完全耦合,因为 ABAQUS 子程序 HETVAL 需要在温度位移网格单元中起作用,因此这里采用完全耦合的方式,分析步选用温度位移耦合、相互作用选择表面热交换。考虑到模型边角淬火过程存在应力集中现象,所以网格采用二次完全积分网格。拉伸过程通过对厚板一端拉伸需要距离实现。卸载过程通过取消激活拉伸过程的加载实现。实际生产过程中,需要对淬火拉伸后的铝合金厚板进行边缘剪裁,这里在厚板中心截取尺寸为 120 mm120 mm50 mm 的区域进行应力测量。在剪裁分析步中通过生死单元法实现这一492步骤。4.3改进的仿真方法4.3.1屈服强度变化铝合金板淬火后在自然
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 铝合金 拉伸 残余 应力 数值 模拟 方法 研究 陈树梁
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。