铝合金预拉伸板残余应力数值模拟方法研究_陈树梁.pdf

基金项目:四川省重点研发计划(2021YFG0060)收稿日期:20210901修回日期:20210912第 40 卷第 2 期计算机仿真2023 年 2 月文章编号:10069348(2023)02029205铝合金预拉伸板残余应力数值模拟方法研究陈树梁，徐雷*，胡元昊，高涵(四川大学机械工程学院，四川 成都 610065)摘要:针对铝合金预拉伸板残余应力仿真结果误差较大的问题，提出一种考虑铝合金材料参数变化的数值模拟方法。铝合金预拉伸板淬火和拉伸过程中屈服强度和弹性模量的变化考虑进来，其中屈服强度的变化通过设置温度阈值实现，弹性模量的变化通过将每个单元积分点弹性模量表示成温度和塑性变形的函数来实现。通过验证发现，上述方法有效地降低了仿真结果与实验数据之间的误差，能很好地预测铝合金板残余应力的大小和变化趋势。通过对照仿真发现淬火与拉伸工艺之间的自然时效时间越长，所提方法的效果越显著。关键词:铝合金预拉伸板;残余应力;数值模拟;自然时效中图分类号:TP391.9文献标识码:BStudy on Numerical Simulation Method of esidualStress of Aluminum Alloy PreStretching Board Based onChange of Material ParametersCHEN Shuliang，XU Lei*，HU Yuanhao，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 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 temperature 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 simulationresults 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 significant the effect of this method isKEYWODS:Aluminum alloy prestretching board;esidual stress;Numerical simulation;Natural aging1引言铝合金尤其是 7xxx 铝合金因其密度小、强度高、硬度高等优良性能，广泛应用在航空航天产品的薄壁结构件中。其刚度低、壁薄、弹性模量小的特点使得残余应力导致的加工变形成为了影响铝合金厚板加工精度的重要因素。而且为了增加薄壁构件的强度及其轻量化，薄壁构件大型化和整体化的趋势也使得加工变形问题更为突出。残余应力主要包括初始残余应力和加工产生残余应力，随着高速铣削工艺的发展和加工参数的优化，加工产生的残余应力越来越小，因此对铝合金预拉伸板初始残余应力的精确数值模拟对薄壁结构件变形的预测和加工工艺调整尤为重要。针对铝合金厚板初始残余应力的研究，曹海龙等研究了7055 铝合金淬火条件对残余应力分布规律的影响1。张园园等对 7075 铝合金做了大量研究，包括淬火的换热系数反求，进行淬火以及预拉伸仿真得到铝合金厚板的初始残余应力分布规律以及残余应力与拉升量和拉伸速率的关系等2。李雪春等认为，一般工程应用中都把弹性模量作为常数，实292际上金属的弹性模量在塑性变形过程中，随着塑性变形程度的增加会发生变化。另外李春雪等研究了铝合金的弹性模量值随塑性变形程度及加载方式变化的规律3。刘瑶琼等，林峰等分析了塑性变形对铝合金弹性模量的影响，通过对模型分层建立了考虑弹性模量差异的 7050 铝合金厚板“淬火预拉伸”有限元模型，分析了弹性模量差异对 7050 铝合金厚板残余应力的影响4，5。实际上，铝合金预拉伸板在淬火与预拉伸整个过程中，材料的一些参数是变化的，其中对预拉伸板初始残余应力影响较大的主要是弹性模量与屈服极限。现有的仿真很少见到考虑屈服极限变化的情况，而考虑弹性模量差异的仿真基本都是采用对模型分层的方法，这种方法将材料进行了离散化处理，虽然减少了仿真的误差，但同样导致仿真结果不连续，与实际不够吻合。针对铝合金材料淬火预拉伸过程中参数变化这个问题，本文以 7075 铝合金预制板为对象进行仿真，仿真模型不仅考虑了铝合金屈服极限的变化，而且将仿真模型的每个积分点的材料弹性模量与其温度、等效塑性变形联系起来，使得仿真结果更加精确。2仿真原理及参数确定铝合金预拉伸厚板的工艺路线如图 1 所示图 1铝合金预拉伸厚板的工艺路线图固溶的时候铝合金厚板内部应力几乎为零，因此铝合金厚板的初始残余应力主要跟淬火、预拉伸和时效过程相关。因此本文主要探究淬火和预拉伸过程。3淬火过程数学模型零件的淬火过程是应该热交换过程，根据能量守恒定律和 Fourier 传热定律，零件的瞬态温度场应满足传热分析的过程控制方程，即xkxtx+ykyty+zkztz+Q cTt=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 铝合金材料的热力学参数如密度、弹性模量、导热系数、比热容以及热膨胀系数都与温度 T 相关。具体参数如表 1 所示 68。表 1铝合金 7075 的热力学基本参数表温度密度kg/m3弹性模量 GPa热导率W/(m)线膨胀系数106/比热J/(kg)泊松比2028007100015521.685010027756519316123.490020027505626217524.39700.3330027253798218525.2102040027003150019330.7112050026752500019731.41320另外淬火过程中的换热系数也是非常重要的参数，这里采用浸没淬火的方式，换热系数见表 2。表 27075 铝合金浸没淬火(水淬)换热系数表温度()换热系数温度()换热系数20130012100133509.515021.54007.520022450725016.55006.753.2应力应变计算原理2.1.1弹性应变阶段弹性阶段应力应变关系采用胡克定律进行描述，即=E(3)在式(3)中，E 为杨氏模量，一般计算都将材料的弹性模量看成常量或者随温度变化的变量进行计算，而实际上，杨氏模量 E 还是关于塑性应变和应变率的变量，由于拉伸过程为低应变率，因此忽略应变率的影响。另外，由塑性变形产生的杨氏模量变化会随着时间由快到慢逐渐恢复9，而且当多次加载和单次加载到达相同塑性变形值时，多次加载明显造成弹性模量下降更多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。对于面心立方结构的铝合金，塑性变形造成的弹性模量变化主要是由变形中形成的主要织构成分含量的变化引起的14。而温度也会影响织构的变化和恢复，但考虑到淬火时间短，且淬火阶段整体上产生的塑性变形只占整个过程的小部分。因此为了方便研究，假设温度与塑性变形都独立作用于杨氏模量，且塑性变形对不同温度下杨氏模量的影响是等比列的，可以将弹性模量表示成式(6)。E=E70.95ET(6)2.1.2塑性应变阶段当材料达到屈服极限进入塑性应变阶段的应力主要与应变、温度相关。实际上，铝合金固溶淬火处理之后的铝合金材料屈服强度在自然时效下会得到提高，其变化趋势15 如图 2 所示。图 2铝合金板淬火后的时效特征图根据图 2 可以看出:室温情况下，在淬火后约 0.3 小时内，屈服强度几乎维持在 135MPa 不变，而在淬火 0.3 小时之后提高的速率由快速到缓慢缓慢转变，经过一个小时的自然时效，屈服强度达到 180MPa 左右，而淬火经过较长的自然时效(100h)屈服强度可以达到 300MPa 左右。但是大部分铝合金仿真研究都是使用固定屈服极限来代替整个过程进行仿真。比如采用淬火拉伸之后材料的屈服极限进行仿真，甚至采用进行人工时效后的屈服极限进行仿真，这样仿真出来的结果跟实际情况有一定的差异。因此对不同阶段采用相对应的屈服极限进行仿真是非常有必要的，为了得到不同阶段铝合金材料的屈服极限特性，查阅文献得到 7075T651 铝合金的屈服极限随温度变化如表 4 所示16。7075 铝合金在屈服强度为 225MPa 时的应力应变数据如表 5 所示4。表 47075T651 铝合金屈服极限的温度特性表温度20100200300400500620屈服极限456389275.747.135.516.30表 57075 铝合金应力应变数据表塑性应变00.0060.0160.0260.034应力 MPa225260300330350因为铝合金材料淬火后不同时效时间应力应变曲线不方便测量，为了方便研究，假设不同自然时效时间的铝合金的软化系数与 7075T651 相同以及不同自然时效时间的铝合金的塑性应变的硬化系数相同17。然后根据前面提到的数据可以推导出不同自然时效铝合金在不同温度下的拉伸应力应变曲线。虽然与实际参数存在误差，但是该研究主要是为了探究改进仿真的可行性和有效性，在误差允许范围内是可以接受的。4仿真模型的构建4.1基本假设影响铝合金预拉伸板材残余应力的因素非常多。为了突出主要影响因素，仿真过程做出如下假定条件:1)材料视为各向同性的弹塑性材料。2)淬火介质温度始终保持恒定。3)忽略淬火转移时间，即不考虑淬火工件专员过程中与空气的热量交换。4)忽略淬火工件进入淬火介质的过程。5)拉伸速率很小，忽略其对残余应力消除的影响。6)视不同厚度铝合金板时效效果相同。4.2模型设置使用 800mm220mm50mm 的 7075 铝合金厚板加热到475固溶处理，然后在 20的水中淬火 300s，淬火后放置不同的时间(屈服极限分别强化到 180MPa，300MPa)再进行拉伸消除残余应力处理，拉伸量分别为 1.5%、2.0%、2.5%。设置四组对比仿真，第一组是考虑屈服极限和弹性模量变化的改进仿真组;第二组是不考虑弹性模量变化以及以拉伸时屈服极限作为整个仿真时间的屈服极限的仿真;第三组是只考虑屈服极限变化的仿真组;第四组是只考虑弹性模量变化的仿真组。为了简化模型，这里假设厚板各个表面同时淬火，温度场，应力场对称分布，所以使用 1/8 模型仿真以减少计算量。仿真采用 4 个分析步，分别对应淬火、拉伸、卸载和剪裁。淬火过程的热力耦合计算的常见方式有顺序耦合和完全耦合，因为 ABAQUS 子程序 HETVAL 需要在温度位移网格单元中起作用，因此这里采用完全耦合的方式，分析步选用温度位移耦合、相互作用选择表面热交换。考虑到模型边角淬火过程存在应力集中现象，所以网格采用二次完全积分网格。拉伸过程通过对厚板一端拉伸需要距离实现。卸载过程通过取消激活拉伸过程的加载实现。实际生产过程中，需要对淬火拉伸后的铝合金厚板进行边缘剪裁，这里在厚板中心截取尺寸为 120 mm120 mm50 mm 的区域进行应力测量。在剪裁分析步中通过生死单元法实现这一492步骤。4.3改进的仿真方法4.3.1屈服强度变化铝合金板淬火后在自然时效下屈服强度会得到提高。淬火后零件内部残余应力都低于屈服强度，所以屈服强度提高过程中不会直接引起应力场的变化，但是提高后的屈服强度会影响拉伸阶段的应力场分布。所以屈服强度提高的快慢不会对仿真结果产生影响，即提高的过程不是所关心的。另外，当铝合金板整体温度与水温很接近时，铝合金板初始残余应力基本不再变化，因此，为了简单有效地实现铝合金板淬火与拉伸过程中屈服极限的变化，设置了一个接近水温的阈值(20.5)，当材料温度低于这个阈值时，屈服极限便转变为时效之后(拉伸时)的屈服极限。以此实现拉伸过程仿真按照实际的屈服极限进行。4.3.2弹性模量变化因为弹性模量不仅仅与温度有关，也与塑性变形有关，在仿真软件已有的功能中不能直接定量地描述弹性模量与温度和塑性变形的关系。但是 ABAQUS 为用户提供了灵活的用户子程序接口，通过这些接口可以实现一些用户需要的功能。通过 HETVAL 与 USDFLD 子程序来获得单元积分点的温度与塑性变形，并将用户定义的场变量与温度和塑性变形定量地联系起来，在仿真属性设置里面将弹性模量定义成场变量的函数，然后通过场变量来实现弹性模量的变化。另外，HETVAL 与 USDFLD 子程序在每个增量步结束时传入积分点的参数，因此在仿真过程中控制增量步足够小，即可认为弹性模量是实时地更新和连续地变化。5仿真结果分析5.1改进仿真有效性验证为了验证改进仿真的有效性，参考了张园园论文中的实验数据，该数据源于“973 课题组”残余应力测试实验。该实验方法中不同拉伸量实验的拉伸过程都在淬火后的一个小时内分别完成2，但没有准确的时间。因此首先对铝合金板淬火后经过一个小时的自然时效，屈服强度达到 180MPa 时进行预拉伸的情况进行了仿真。为了保证结果对比的严谨性，仿真取值的位置与试验测量位置相同。将采用改进仿真与未改进仿真的应力峰值结果与实际的试验测量结果进行对比。以轧向应力为例，厚板的轧向和横向对称中心位置沿厚度方向的最大拉压力与最大压应力试验测量结果与仿真结果见表 6。表 6仿真与试验应力峰值数据对比表拉伸量应力来源最大拉应力 MPa最大压应力 MPa1.5%实验16.6418.85未改进仿真27.7127.91改进仿真22.0624.362.0%实验14.2915.98拉伸量应力来源最大拉应力 MPa最大压应力 MPa未改进仿真19.6321.48改进仿真16.1019.212.5%实验13.9515.93未改进仿真15.5317.71改进仿真12.4416.24通过上表求得未改进仿真结果与改进的仿真结果误差，其中正数表示仿真值大于试验值，负数表示仿真值小于于试验值，结果如表 7 所示。表 7改进与未改进仿真应力峰值误差表拉伸量应力来源最大拉应力误差 最大压应力误差1.5%未改进仿真66.53%48.06%改进仿真32.57%29.23%2.0%未改进仿真37.37%34.42%改进仿真12.67%20.21%2.5%未改进仿真11.33%11.17%改进仿真10.82%1.95%从表 7 中可以看出，改进的仿真比未改进的仿真峰值应力误差平均消除了 50%左右，说明改进的仿真方法对残余应力的预测更加准确。不过三种不同的拉伸量误差存在较大的差异，考虑到仿真是在淬火经过一个小时自然时效后进行拉伸，而实验是在一个小时内分别完成，拉伸的顺序有先后。而且，淬火后的一个小时内。铝合金材料的屈服极限变化较大。因此猜测误差的差异一部分原因是由于自然时效时间不同导致拉伸时屈服极限存在较大差异造成的。为了验证这一猜测，假设淬火后屈服极限未变化之前进行拉伸试验，对此进行了仿真，其中改进仿真结果与试验结果误差见表 8。表 8改进仿真应力峰值误差表拉伸量应力类型数值 MPa误差1.5%最大拉压力18.5511.48%最大压应力21.5414.27%2.0%最大拉压力14.021.89%最大压应力17.8511.70%2.5%最大拉压力11.4318.06%最大压应力15.194.65%从表 8 中结果可以看出拉伸量 1.5%和拉伸量 2.0%的误差有了明显的降低，且误差相近。而拉升量 2.5%的误差反而更大。因此结果验证了猜测，即拉伸量 1.5%、2%的试验发生在拉伸量 2.5%之前。综合来看，虽然仿真与实验的时效时间存在差异，但是改进仿真相比未改进仿真能有效地减少误差。考虑到这只是应力峰值的对比，为了验证厚度方向应力变化趋势的准确性，以拉伸量 2.5%轧制方向应力为例，将改进的仿真结果与试验结果及采用的分层模型仿真进行比592较12，结果如图 3 所示。图 3轧制方向应力沿厚度方向变化趋势对比图从图 3 中可以看出改进之后的仿真与试验结果趋势基本一致，而且在厚度 010mm 厚度区间，改进的仿真比采用分层模型的仿真能更好预测应力的趋势。通过应力峰值与应力趋势的对比，可以认为改进的仿真能够更好的预测铝合金预拉伸板的初始残余应力。5.2不同自然时效时间两种仿真差异对比前面是在自然时效 1 小时进行拉伸的仿真。为了探究改进的仿真和未改进仿真在不同自然时效时间后进行拉伸的差异。将淬火后自然时效 1 小时(屈服强度 180MPa)进行预拉伸与自然时效 100 小时(屈服强度 300MPa)进行预拉伸的仿真结果进行对比，以拉伸量 2.5%为例，其中心位置轧制方向的应力沿厚度方向趋势比较如图 4 所示。图 4不同自然时效时间两种仿真差异对比图从图中可以看出:1)随着自然时效时间的增加(屈服极限愈大)，应力值总体上变大。所以，拉伸应该在淬火后越快进行越好。2)随着自然时效时间的增加，改进仿真与未改进仿真的结果差异也越大，尤其是靠近表面处的应力差异巨大。3)随着自然时效时间的增加，改进仿真和为改进仿真压应力峰值位置点也出现较大的差异。因此淬火与拉伸之间的间隔越久，改进仿真与未改进仿真差异越大，使用改进的仿真就越有必要。5.3两种材料参数变化分别对仿真结果影响分析改进的仿真考虑了弹性模量和屈服强度的变化，为了探究改进仿真中引入的弹性模量变化与屈服强度变化分别对误差消除的贡献，因为自然时效越长，改进仿真与未改进仿真差 异 越 明 显，因 此 以 自 然 时 效 100 小 时(屈 服 强 度300MPa)为例，采用控制变量法分别做了只改进弹性模量变化和只改进屈服强度变化的对比仿真，不同拉伸量的效果类似，以拉伸量 2.5%为例，将两组对照仿真结果与改进仿真和未改进仿真结果进行对比，结果如图 5 所示。图 5不同改进仿真与未改进仿真结果对比图从图 5 可以看出只改进屈服强度与只改进弹性模量的仿真都能使总体应力得到大幅减小，并且使得压应力峰值右移。而且只改进弹性模量的仿真比只改进屈服极限的仿真效果更明显。对于应力降低的原因，只改进弹性模量的仿真是因为改进之后的弹性模量受塑性变形的影响而降低，从而导致了应力总体下降。而只改进屈服极限的仿真是因为淬火过程屈服强度降低导致更小的弹性变形和更大的塑性变形，而残余应力产生主要是因为弹性变形不均匀产生的，因此更小的弹性变形降低了残余应力。另外，从图 5 可以看出相比于只改进弹性模量，改进仿真结果与只改进屈服强度结果更接近，这是因为改进仿真考虑了屈服强度的变化，使得弹性变形相比只改进弹性模量仿真更小，虽然塑性变形更大导致弹性模量降低。但是自然时效较长的情况下弹性变形变化的影响更大，从而导致改进仿真结果与只改进屈服强度结果更接近。6结论本文针对铝合金预拉伸板残余应力仿真误差较大的问题，提出一种考虑铝合金材料参数变化的数值模拟方法，将铝合金预拉伸板淬火和拉伸过程中屈服强度和弹性模量的变化考虑进来。并得到以下结论:1)改进的仿真有效地降低了误差，能很好地预测铝合金板残余应力的大小和变化趋势。对于铝合金厚板淬火拉伸工艺以及铝合金板加工变形预测都具有指导意义。2)淬火与拉伸工艺之间的自然时效时间越长，应力值总体也越大。所以，拉伸应该在淬火后越快进行越好。(下转第 516 页)692Computer Vision and Pattern ecognition Workshops(CVPW)IEEE Computer Society，2017 17Li B，en W，Fu D，et al Benchmarking Single Image Dehazingand Beyond J IEEE Transactions on Image Processing，2018 18X Liu，Y Ma，Z Shi，et al GridDehazeNet:AttentionBased MultiScale Network for Image DehazingC 2019 IEEE/CVF Inter-national Conference on Computer Vision(ICCV)，2019:73137322 作者简介韦照川(1973)，男(瑶族)，四川省成都市人，副教授，主要研究领域为信号处理。王皓坤(1995)，男(汉族)，四川省西昌市人，硕士生，主要研究领域为图像处理。纪元法(1975)，男(汉族)，山东省莘县人，教授，硕士研究生导师，主要研究领域为信号处理、卫星导航。(上接第 296 页)3)淬火与拉伸工艺之间的自然时效时间越长，改进仿真与未改进仿真差异越大，使用改进的仿真就越有必要。4)弹性模量和屈服强度的变化对误差的消除都起到重要的作用。参考文献:1Hai Long Cao，et al Numerical Simulation of Quenching and PreStretching esidual Stress in 7085 Aluminum Alloy PlateJ Ma-terials Science Forum，2016，3960:211217 2张园园 铝合金厚板淬火过程及预拉伸热力仿真与实验研究 D 中南大学，2008 3李雪春，杨玉英，王永志，包军 塑性变形对铝合金弹性模量的影响J 中国有色金属学报，2002，(4):701705 4刘瑶琼，龚海，孙燕杰 7050 铝合金厚板弹性模量分布差异及其对预拉伸后残余应力的影响J 热加工工艺，2019，48(24):2529，36 5Feng Lin，et al Investigation on esidual Stresses in Pre Stretching Process of Aluminum Alloy Thick Plate Based onMaterial NonUniformityJ Key Engineering Materials，2016，4283:821829 6朱才朝，罗家元，钟渝 考虑夹持影响的铝合金板拉伸模拟及试验 J 材料科学与工艺，2011，19(6):1621 7柯映林，董辉跃 7075 铝合金厚板预拉伸模拟分析及其在淬火残余应力消除中的应用 J中国有色金属学报，2004，14(4):639645 8Kethireddy Narender，et al Temperature Dependence of Densityand Thermal Expansion of Wrought Aluminum Alloys 7041，7075and 7095 by Gamma ay Attenuation MethodJ Journal ofModern Physics，2013，4(3):331336 9Villuendas，Aranzazu;Jorba，Jordi;oca，Antoni The ole of Pre-cipitates in the Behavior of Young s Modulus in Aluminum Alloys J Metallurgical and Materials Transactions A，2014，45(9):38573865 10Liu Y L，Zhu Y X，Dong W Q，Yang H Springback predictionmodel considering the variable Y oung s modulus for the bendingrectangular 3A21 tubeJ Journal of Materials Engineering andPerformance，2013，22(1):916 11Morestin F，Boivin M On the necessity of taking into account thevariation in the Y oung modulus with plastic strain in elasticplastic softwareJ Nuclear Engineering and Design，1996，162(1):107116 12秦国华，林锋，左敦稳 基于材料非均匀性铝厚板预拉伸残余应力的预测方法J 中国有色金属学报，2015，25(3):626633 13徐志东，范子亮 金属材料的弹性模量随温度变化规律的唯象解释J 西南交通大学学报，1993，(2):8792 14李雪春，杨玉英，王永志，于忠奇 塑性变形中织构的演化规律及其对弹性模量的影响C 中国机械工程学会、北京机械工程学会制造业与未来中国 2002 年中国机械工程学会年会论文集，中国机械工程学会、北京机械工程学会:中国机械工程学会，2002:1 15梁世斌 铝合金挤压加工技术手册M 长沙:中南大学出版社，2019:632 16李大峰，丁华锋，刘立斌，王孟，蔡奎 7075 铝合金板淬火残余应力模拟及实验研究J 机械研究与应用，2012，25(3):9295 17Zhihua Chen，et al Experimental investigation on the postfiremechanical properties of structural aluminum alloys 6061T6 and7075T73J ThinWalled Structures，2016，106:187200 作者简介陈树梁(1997)，男(汉族)，湖南衡阳人，硕士研究生，主要研究方向:智能制造、计算机仿真。徐雷(1973)，女(汉族)，四川成都人，副教授，博士，主要研究方向:智能制造、工业机器人。胡元昊(1997)，男(汉族)，四川绵阳人，硕士研究生，主要研究方向:智能制造、产线平衡。高涵(1996)，男(汉族)，湖北孝感人，硕士研究生，主要研究方向:智能制造、结构优化。615