CN105046031B - 一种预测铝合金薄壁构件机加工产生的形变量的方法 - Google Patents
一种预测铝合金薄壁构件机加工产生的形变量的方法 Download PDFInfo
- Publication number
- CN105046031B CN105046031B CN201510574547.7A CN201510574547A CN105046031B CN 105046031 B CN105046031 B CN 105046031B CN 201510574547 A CN201510574547 A CN 201510574547A CN 105046031 B CN105046031 B CN 105046031B
- Authority
- CN
- China
- Prior art keywords
- mrow
- residual stress
- msup
- machining
- wall member
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Investigating And Analyzing Materials By Characteristic Methods (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种预测铝合金薄壁构件机加工产生的形变量的方法,涉及预测铝合金薄壁构件机加工产生的形变量的技术,属于机械加工领域。它为了解决现有技术对铝合金薄壁构件的整个机加工过程进行模拟的计算量很大的问题,测试坯材的残余应力和铝合金薄壁构件机加工的表层残余应力,用等效表层残余应力的分布代替表层残余应力的分布,将坯材的残余应力输入到ABAQUS有限元软件,去除机加工过程中需要去除的部位,将等效表层残余应力施加到ABAQUS分析模型,计算由于坯材的残余应力释放和机加工表层残余应力共同作用下的铝合金薄壁构件形变量,从而可以预测铝合金薄壁构件机加工产生的形变量。本发明可用于预测铝合金薄壁构件机加工产生的形变量。
Description
技术领域
本发明涉及预测铝合金薄壁构件机加工变形的技术,属于机械加工领域。
背景技术
随着航空、航天技术的发展,航空航天已成为高精尖技术密集的行业。铝合金薄壁或框架类整体构件被越来越广泛的应用于航空航天、武器装备、精密仪器等行业,零件精度的要求也不断提高。铝合金薄壁构件受坯材残余应力以及机加工残余应力影响,易产生很大的加工变形,加工变形是航空航天结构件制造技术所面对的最突出问题之一。随着计算机技术的发展,利用计算机计算铝合金薄壁构件加工过程中的变形成为目前科研技术开发和研究的重要途径,但由于机加工过程复杂,很难通过对铝合金薄壁构件的整个机加工过程进行模拟的方式来预测机加工产生的形变量,且计算量很大。
发明内容
本发明是为了解决由于机加工过程复杂且计算量大,导致难以用计算机对铝合金薄壁构件的整个机加工过程进行模拟来预测机加工形变量的问题,从而提供一种预测铝合金薄壁构件机加工产生的形变量的方法。
本发明所述的一种预测铝合金薄壁构件机加工产生的形变量的方法,该方法包括以下步骤:
步骤一:采用裂纹柔度法对铝合金坯材的残余应力进行测试,获取铝合金坯材的残余应力沿深度方向的分布,对所述残余应力的分布进行拟合;
步骤二:采用XRD(X射线衍射)残余应力测试法获得铝合金薄壁构件机加工的表层残余应力沿深度方向的分布,对表层残余应力的分布进行拟合;
步骤三:采用新的表层残余应力σ即等效表层残余应力的分布代替步骤二得到的拟合后的表层残余应力的分布,所述新的表层残余应力的分布与步骤二得到的拟合后的表层残余应力对中性面的弯矩相同,
等效表层残余应力的分布为:
其中
d为深度值,z为深度变量,式中的数值大小等于机加工的表层残余应力对中性面的等效作用力,σmachine为步骤二中拟合后的表层残余应力;
步骤四:建立铝合金薄壁构件的ABAQUS有限元模型,将步骤一得到的拟合后的铝合金坯材的残余应力通过SIGINI子程序文件输入到ABAQUS有限元软件,去除机加工过程中需要去除的部位;
步骤五:将步骤三中的等效表层残余应力σ施加到ABAQUS有限元模型中,分布在层深为d的表层单元,计算由于铝合金坯材的残余应力释放和机加工表层残余应力共同作用下的构件的形变量。
上述步骤一中对残余应力分布的拟合采用勒让德函数实现;
勒让德函数为:
式中z'为归一化的深度;
铝合金坯材的残余应力σ'沿深度方向的分布拟合为:
式中Ai为拟合系数。
上述步骤二中对表层残余应力分布的拟合采用拉格朗日多项式函数实现。
上述步骤四中去除机加工过程中需要去除的部位采用生死单元法实现。
本发明所述的一种预测铝合金薄壁构件机加工产生的形变量的方法,测试坯材的残余应力和铝合金薄壁构件机加工的表层残余应力,用等效表层残余应力的分布代替表层残余应力的分布,将坯材的残余应力输入到ABAQUS有限元软件,去除机加工过程中需要去除的部位,将等效表层残余应力σ施加到ABAQUS分析模型,计算由于坯材的残余应力释放和机加工表层残余应力共同作用下的铝合金薄壁构件的形变量,从而可以预测铝合金薄壁构件机加工的形变量,本发明可以准确预测铝合金薄壁构件机加工产生的形变量,很大程度的减小预测的计算量,可以针对复杂构件的变形规律进行分析,优化机加工参数,采用本方法进行预测的误差最大仅为13%。
附图说明
图1是具体实施方式一所述的一种预测铝合金薄壁构件机加工产生的形变量的方法的原理框图。
图2是具体实施方式五中的210x100x20mm的2A14-T6铝合金坯材的残余应力及拟合的坯材的残余应力分布图。
图3为具体实施方式五中的两种铣削方式的示意图,通常构件在铣削过程中的铣削区域可分为两种,一种是侧壁的铣削,A代表侧壁区域,一种是平面区域的铣削,B代表平面区域,分别测试这两种铣削方式的表层残余应力。
图4为具体实施方式五中的A的机加工表层残余应力分布及其拟合曲线图。
图5为具体实施方式五中的B的机加工表层残余应力分布及其拟合曲线图。
图6为具体实施方式五中的预测的铝合金薄壁构件机加工变形的分布图,C代表底面。
图7为具体实施方式五中的底面机加工变形的实验结果与模拟结果的对比图。
具体实施方式
具体实施方式一:参照图1具体说明本实施方式,一种预测铝合金薄壁构件机加工产生的形变量的方法,该方法包括以下步骤:
步骤一:采用裂纹柔度法对铝合金坯材的残余应力进行测试,获取铝合金坯材的残余应力沿深度方向的分布,对所述残余应力的分布进行拟合;
步骤二:采用XRD(X射线衍射)残余应力测试法获得铝合金薄壁构件机加工的表层残余应力沿深度方向的分布,对表层残余应力的分布进行拟合;
步骤三:采用新的表层残余应力σ即等效表层残余应力的分布代替步骤二得到的拟合后的表层残余应力的分布,所述新的表层残余应力的分布与步骤二得到的拟合后的表层残余应力对中性面的弯矩相同,
等效表层残余应力的分布为:
其中
式中的数值大小等于机加工的表层残余应力对中性面的等效作用力,σmachine为步骤二中拟合后的表层残余应力;
步骤四:建立铝合金薄壁构件的ABAQUS有限元模型,将步骤一得到的拟合后的铝合金坯材的残余应力通过SIGINI子程序文件输入到ABAQUS有限元软件,去除机加工过程中需要去除的部位;
步骤五:将步骤三中的等效表层残余应力σ施加到ABAQUS有限元模型中,分布在层深为d的表层单元,计算由于铝合金坯材的残余应力释放和机加工表层残余应力共同作用下的构件的形变量。
测试深度接近表层残余应力几乎零的深度。
具体实施方式二:本实施方式是对具体实施方式一所述的一种预测铝合金薄壁构件机加工产生的形变量的方法作进一步说明,本实施方式中,步骤一中对残余应力分布的拟合采用勒让德函数实现;
勒让德函数为:
式中z'为归一化的深度;
铝合金坯材的残余应力σ'沿深度方向的分布拟合为:
式中Ai为拟合系数。
具体实施方式三:本实施方式是对具体实施方式一所述的一种预测铝合金薄壁构件机加工产生的形变量的方法作进一步说明,本实施方式中,步骤二中对表层残余应力分布的拟合采用拉格朗日多项式函数实现。
具体实施方式四:本实施方式是对具体实施方式一所述的一种预测铝合金薄壁构件机加工产生的形变量的方法作进一步说明,本实施方式中,步骤四去除机加工过程中需要去除的部位采用生死单元法实现。
具体实施方式五:本实施方式是对具体实施方式一所述的一种预测铝合金薄壁构件机加工产生的形变量的方法的验证。
图2是采用裂纹柔度法获得的210×100×20mm的2A14-T6铝合金坯材的残余应力及拟合的坯材的残余应力分布图,将此应力作为初始残余应力施加到坯材的ABAQUS有限元模型中;图3为铣削示意图,分别测试这两种铣削方式的表层残余应力,作为构件机加工引入的表层残余应力;图4为A的机加工表层残余应力分布及其拟合曲线图,根据拟合的曲线,通过公式(2)计算A的等效表层残余应力,图5为B的机加工表层残余应力分布及其拟合曲线图,根据拟合的曲线,通过公式(2)计算B的等效表层残余应力,图6为预测的铝合金薄壁构件机加工变形的分布图,每个区域的形变量如表1所示,采用三坐标测量仪测量了底面的变形,通过测量得到的实验结果与采用本发明的方法得到的模拟结果对比图为图7,从图中可以看出,模拟结果和实验结果吻合较好,模拟结果和实验结果最大偏差为13%,横坐标中的长向代表图6中的x方向。
表1
区域 | 形变量(mm) |
1 | 0.3187 |
2 | 0.1932 |
3 | 0.6774 |
4 | -0.5777 |
5 | -0.1833 |
6 | -0.3088 |
7 | -0.4343 |
Claims (4)
1.一种预测铝合金薄壁构件机加工产生的形变量的方法,其特征在于,该方法包括以下步骤:
步骤一:采用裂纹柔度法对铝合金坯材的残余应力进行测试,获取铝合金坯材的残余应力沿深度方向的分布,对所述残余应力的分布进行拟合;
步骤二:采用XRD残余应力测试法获得铝合金薄壁构件机加工的表层残余应力沿深度方向的分布,对表层残余应力的分布进行拟合;
步骤三:采用新的表层残余应力σ即等效表层残余应力的分布代替步骤二得到的拟合后的表层残余应力的分布,所述新的表层残余应力的分布与步骤二得到的拟合后的表层残余应力对中性面的弯矩相同,
等效表层残余应力的分布为:
<mrow>
<mi>&sigma;</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mover>
<mi>&sigma;</mi>
<mo>&OverBar;</mo>
</mover>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>&le;</mo>
<mi>z</mi>
<mo>&le;</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>&GreaterEqual;</mo>
<mi>d</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中
<mrow>
<mover>
<mi>&sigma;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>d</mi>
</msubsup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>c</mi>
<mi>h</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
</mrow>
</msub>
<mi>d</mi>
<mi>z</mi>
</mrow>
<mi>d</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
d为深度值,z为深度变量,式中的数值大小等于机加工的表层残余应力对中性面的等效作用力,σmachine为步骤二中拟合后的表层残余应力;
步骤四:建立铝合金薄壁构件的ABAQUS有限元模型,将步骤一得到的拟合后的铝合金坯材的残余应力通过SIGINI子程序文件输入到ABAQUS有限元软件,去除机加工过程中需要去除的部位;
步骤五:将步骤三中的等效表层残余应力σ施加到ABAQUS有限元模型中,分布在层深为d的表层单元,计算由于铝合金坯材的残余应力释放和机加工表层残余应力共同作用下的构件的形变量。
2.根据权利要求1所述的一种预测铝合金薄壁构件机加工产生的形变量的方法,其特征在于,步骤一中,对残余应力分布的拟合采用勒让德函数实现;
勒让德函数为:
<mrow>
<msub>
<mi>L</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>z</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>d</mi>
<mrow>
<msup>
<mi>dz</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>z</mi>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>z</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>n</mi>
<mo>!</mo>
</mrow>
</mfrac>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>&le;</mo>
<msup>
<mi>z</mi>
<mo>&prime;</mo>
</msup>
<mo>&le;</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>&le;</mo>
<mi>n</mi>
<mo>&le;</mo>
<mn>11</mn>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中z'为归一化的深度;
铝合金坯材的残余应力σ'沿深度方向的分布拟合为:
<mrow>
<msup>
<mi>&sigma;</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msup>
<mi>z</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mn>11</mn>
</munderover>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<msub>
<mi>L</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>z</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中Ai为拟合系数。
3.根据权利要求1所述的一种预测铝合金薄壁构件机加工产生的形变量的方法,其特征在于,步骤二中,对表层残余应力分布的拟合采用拉格朗日多项式函数实现。
4.根据权利要求1所述的一种预测铝合金薄壁构件机加工产生的形变量的方法,其特征在于,步骤四中去除机加工过程中需要去除的部位采用生死单元法实现。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510574547.7A CN105046031B (zh) | 2015-09-10 | 2015-09-10 | 一种预测铝合金薄壁构件机加工产生的形变量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510574547.7A CN105046031B (zh) | 2015-09-10 | 2015-09-10 | 一种预测铝合金薄壁构件机加工产生的形变量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105046031A CN105046031A (zh) | 2015-11-11 |
CN105046031B true CN105046031B (zh) | 2018-01-30 |
Family
ID=54452572
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510574547.7A Active CN105046031B (zh) | 2015-09-10 | 2015-09-10 | 一种预测铝合金薄壁构件机加工产生的形变量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105046031B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113237583B (zh) * | 2021-05-13 | 2022-03-15 | 中南大学 | 一种镁合金筒形件残余应力评估预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887472A (zh) * | 2009-05-12 | 2010-11-17 | 通用汽车环球科技运作公司 | 预测淬火铝铸件中残余应力和变形的方法 |
CN102831274A (zh) * | 2012-08-30 | 2012-12-19 | 中南林业科技大学 | 铝合金厚板淬火应力场分布的修正方法 |
CN103294862A (zh) * | 2013-05-29 | 2013-09-11 | 上海飞机制造有限公司 | 碳纤维增强树脂基复合材料加工变形的仿真模拟方法 |
CN104077442A (zh) * | 2014-06-24 | 2014-10-01 | 南京晨光集团有限责任公司 | 基于有限元分析的大型整体薄壁件加工精度控制方法 |
-
2015
- 2015-09-10 CN CN201510574547.7A patent/CN105046031B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887472A (zh) * | 2009-05-12 | 2010-11-17 | 通用汽车环球科技运作公司 | 预测淬火铝铸件中残余应力和变形的方法 |
CN102831274A (zh) * | 2012-08-30 | 2012-12-19 | 中南林业科技大学 | 铝合金厚板淬火应力场分布的修正方法 |
CN103294862A (zh) * | 2013-05-29 | 2013-09-11 | 上海飞机制造有限公司 | 碳纤维增强树脂基复合材料加工变形的仿真模拟方法 |
CN104077442A (zh) * | 2014-06-24 | 2014-10-01 | 南京晨光集团有限责任公司 | 基于有限元分析的大型整体薄壁件加工精度控制方法 |
Non-Patent Citations (3)
Title |
---|
Minimization of Residual Stress in an Al-Cu Alloy Forged Plate b;Ya-Bo Dong etal.;《Journal of Materials Engineering and Performance》;20150421;第2256-2265页 * |
Reducing residual stress in 2014 aluminium alloy die forgings;D.A.Tanner eatl.;《Materials and Design》;20070718;第1489-1496页 * |
基于裂纹柔度法的铝合金预拉伸板;唐志涛 等;《中国有色金属学报》;20070930;第17卷(第9期);第1404-1409页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105046031A (zh) | 2015-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108304657B (zh) | 基于有限元的机床关重件残余应力多工序连续建模仿真方法 | |
CN108182325B (zh) | 一种薄壁结构件切削加工变形预测分析方法 | |
Gang | Study on deformation of titanium thin-walled part in milling process | |
D’Alvise et al. | Modelling of part distortion due to residual stresses relaxation: an aeronautical case study | |
Denkena et al. | Experimental investigation and modeling of thermal and mechanical influences on shape deviations in machining structural parts | |
CN107657129B (zh) | 基于装夹力监测的薄壁件残余应力变形感知预测方法 | |
Liu et al. | Micro-cutting of single-crystal metal: Finite-element analysis of deformation and material removal | |
Thabet et al. | CFD simulation of the air flow around a car model (Ahmed body) | |
Wang et al. | A multi-region model for numerical simulation of micro bulk forming | |
Haichao et al. | A machining position optimization approach to workpiece deformation control for aeronautical monolithic components | |
CN104182585A (zh) | 硅胶弹性体手感分析及疲劳寿命预测的有限元模拟方法 | |
Cerutti et al. | Prediction of machining quality due to the initial residual stress redistribution of aerospace structural parts made of low-density aluminium alloy rolled plates | |
Li et al. | Modeling and simulation of machining distortion of pre-bent aluminum alloy plate | |
Zhan et al. | A high efficient surface-based method for predicting part distortions in machining and shot peening | |
Denkena et al. | Simulation of residual stress related part distortion | |
CN114878046A (zh) | 一种测量厚板焊接件内部残余应力的方法 | |
Al-Athel et al. | Eulerian volume of solid (VOS) approach in solid mechanics and metal forming | |
CN105046031B (zh) | 一种预测铝合金薄壁构件机加工产生的形变量的方法 | |
Neugebauer et al. | Roller Hemming Simulation: State of the art and application limits | |
CN111967146B (zh) | 基于变形测试的钛合金薄壁件残余应力分布反向识别方法 | |
Zheng et al. | Study of distortion on milled thin-wall aluminum parts influenced by initial residual stress and toolpath strategy | |
Yazdani | Steady and unsteady numerical analysis of the DrivAer model | |
Wang et al. | Machining distortion in the milling of multi-frame components | |
Yadav et al. | Controlling deformations of thin-walled Al 6061-T6 components by adaptive clamping | |
Li et al. | Identification of material parameters from punch stretch test |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |