abaqus-铝合金A357切削加工有限元模拟

1、铝合金A357切削加工有限元模拟 Design and Research of General Trajectory Corrector For an Aircraft汇报人：张怡雯 By：ABAQUS.1.abaqus商用仿真软件中限元模型建立铝合金A357切削加工有限元模型论文主要内容模拟结果分析论文主要内容123.2.论文主要内容 1.铝合金A357切削加工有限元模型装夹条件；装夹条件；切削参；切削参； 刀具的儿何参数；刀具的儿何参数；切削路径切削路径 影影 响响 因因 素素 忽略加工过程中，由于温度变化引起的金相组忽略加工过程中，由于温度变化引起的金相组织及其它的化学变化织及其它的化学

2、变化;被加工对象的材料是各向同性的；被加工对象的材料是各向同性的； 刀具是刚体且锋利，只考虑刀具的温度传导刀具是刚体且锋利，只考虑刀具的温度传导;不考虑刀具、工件的振动；不考虑刀具、工件的振动； 假假 设设 条条 件件.3.由于刀具和工件的切削厚度方向上，切削工程由于刀具和工件的切削厚度方向上，切削工程中层厚不变，所以按平面应变来模拟；中层厚不变，所以按平面应变来模拟；论文主要内容 1.铝合金A357切削加工有限元模型.4.1.1材料模型材料模型1.1.1 A357的Johnson-Cook本构模型 材料本构模型用来描述材料的力学性质，表征材料变形过程中的动态响应。在材料微观组织结构一定的情况

3、下，流动应力受到变形程度、变形速度、及变形温度等因素的影响非常显著。这些因素的任何变化都会引起流动应力较大的变动。因此材料本构模型一般表示为流动应力与应变、应变率、温度等变形参数之间的数学函数关系。建立材料本构模型，无论是在制定合理的加工工艺方面，还是在金属塑性变形理论的研究方面都是极其重要的。在以塑性有限元为代表的现代塑性加工力学中，材料的流动应力作为输入时的重要参数，其精确度也是提高理论分析可靠度的关键。在本课题研究中，材料本构模型是切削加工数值模拟的必要前提，是预测零件铣削加工变形的重要基础，只有建立了大变形情况下随应变率和温度变化的应力应变关系，才能够准确描述材料在切削加工过程的塑性变

4、形规律，继而才能在确定的边界条件和切削载荷下预测零件的变形大小及趋势。论文主要内容 1.铝合金A357切削加工有限元模型.4.1.1.1 A357的Johnson-Cook本构模型 在切削过程中，工件在高温、大应变下发生弹塑性变形，被切削材料在刀具的作用下变成切屑时的时间很短，而且被切削层中各处的应变、应变速率和温度并不均匀分布且梯度变化很大。因此能反映出应变、应变速率、温度对材料的流动应力影响的本构方程，在切削仿真中极其关键。当前常用的塑性材料本构模型主要有：Bodner-Paton、Follansbee-Kocks、Johnson-Cook、 Zerrilli-Armstrong等模型，而

5、只有Johnson-Cook模型描述材料高应变速率下热粘塑性变形行为。JohnsonCook模型认为材料在高应变速率下表现为应变硬化、应变速率硬化和热软化效应，JohnsonCook模型如下所示： 01ln1mnrmrT TA BcTT 式中第一项描述了材料的应变强化效应，第二项反映了流动应力随对数应变速率增加的关系，第三项反映了流动应力随温度升高指数降低的关系。 、Tr分别表示参考应变速率和参考温度，Tm为材料熔点。式中A、B、n、C、m、D、k是7 个待定参数；A、B、n表征材料应变强化项系数；C表征材料应变速率强化项系数；m表征材料热软化系数； , 分别为常温材料熔点。论文主要内容 1.

6、铝合金A357切削加工有限元模型.4.1.1.2 材料失效准则 实现切屑从工件分离，本文采用的是剪切失效模型。剪切失效模型是基于等效塑性应变在积分点的值，当损伤参数达到1时，单元即失效，失效参数定义如下： 式中： 为失效参数， 为等效塑性应变初始值， 为等效塑性应变增量， 为失效应变。失效应变 设定以来于以下几个方面：依据塑性应变率，无量纲压应力与偏应力之比p/q(p为压应力，q为Mises应力),温度，预定义域变量。这里采用JohnsonCook模型定义失效应变。 0plplplf论文主要内容 1.铝合金A357切削加工有限元模型.4.1.1.2 材料失效准则 实现切屑从工件分离，本文采用的

7、是剪切失效模型。剪切失效模型是基于等效塑性应变在积分点的值，当损伤参数达到1时，单元即失效，失效参数定义如下： 式中： 为失效参数， 为等效塑性应变初始值， 为等效塑性应变增量， 为失效应变。失效应变 设定以来于以下几个方面：依据塑性应变率，无量纲压应力与偏应力之比p/q(p为压应力，q为Mises应力),温度，预定义域变量。这里采用JohnsonCook模型定义失效应变。 0plplplf论文主要内容 1.铝合金A357切削加工有限元模型.4.1.1.2 材料失效准则 12345exp1ln1plopdddddq p lf= 式中： 为低于转变温度的条件下测得的实效常数。 为参考应变率， 为

8、塑性应变率。 由下式确定： 0()/()1transitiontransitonmelttransitiontransitionmeltmeltforforfor 是当前温度， 是熔点， 是室温。 melt transition论文主要内容 1.铝合金A357切削加工有限元模型.4.1.1.2 材料失效准则 plyo0和 图1.1中 和 为材料开始损伤时的屈服应力和等效塑性应变。 是材料失效时即图中D=1时的等效塑性应变。材料失效时的等效塑性应变 依赖于单元的特征长度，不能作为描述材料损伤演化的准则。相反，材料损伤演化的准则又等效塑性位移 或者断裂耗散能量 决定。 yo 0p l 0p l 0

9、p l plu fG论文主要内容 1.铝合金A357切削加工有限元模型.4.1.1.2 材料失效准则 当材料开始损伤破坏时，应力应变曲线已经不能准确的描述材料的行为。继续应用该应力应变曲线会导致应变集中，变化过于依赖建模时所画的网格，以致当网格变密后耗散能量反而降低。Hillerborg能量失效法被提出用应力位移响应曲线来表征破坏过程减少了分析对网格的依赖性。利用脆性断裂概念定义一个使单元破坏的能量Gf作为材料的参数。通过这种方法，损伤开始的软化效应是一种应力位移响应而不是应力应变响应。破坏能量有下式表示： 表达式中的 为等效塑性位移，它描述了当损伤开始之后裂纹变化的屈服应力，在损伤开始之前

10、=0.在损伤开始之后 =L ，L为与积分点相关的单元特征长度，单元特征长度的定义基于单元的集合形状，平面单元长度为积分点区域面积的平方根，而立方体单元长度为积分点区域体积的立方根。 0fGplplffplpluplyyoL ddu plu plu plu pl论文主要内容 1.铝合金A357切削加工有限元模型.4.1.1.2 材料失效准则 基于有效塑性位移定义损伤演化用Linear方法定义即如下图所示： 该准则使有效塑性位移达到 = 时，材料的刚度完全丧失，模型的失效网格被自动删除，也就是材料此时发生断裂，切屑开始形成。 plu plfu论文主要内容 1.铝合金A357切削加工有限元模型.4.

11、1.1.3 A357与刀具材料参数 A357铝合金，密度=2680Kg/m3,弹性模量E=79GP，泊松比=0.33其他参数如下表： 论文主要内容 1.铝合金A357切削加工有限元模型.4.1.1.3 A357与刀具材料参数 A357铝合金，密度=2680Kg/m3,弹性模量E=79GP，泊松比=0.33其他参数如下表： 刀具使用的是硬质合金，密度=15000Kg/m3,弹性模量E=210GP，泊松比=0.22其其它参数如下表： 论文主要内容 1.铝合金A357切削加工有限元模型.4.1.2 摩擦模型摩擦模型 金属切削过程中，刀具前刀面的摩擦状态非常复杂，通常把前刀面得摩擦区分为粘结区和滑动区

12、，粘结区的摩擦状态与材料的临界剪应力有关，滑动区可近似认为摩擦系数为常值可以用下式来表示： 式中： 为接触面的滑动剪切应力； 为摩擦系数； 为接触面上的压力； 为材料的临界屈服压力。 cnsm in （ ，） c n s论文主要内容 1.铝合金A357切削加工有限元模型.4.1.2.1 质量放大 论文主要内容 1.研究意义修正引信修正引信零件加工难度大；零件加工难度大；阻力环面积小，弹道修正能力受限制阻力环面积小，弹道修正能力受限制 修正模块微型化技术限制；修正模块微型化技术限制；迫弹修正机构迫弹修正机构结构简单，加工容易；结构简单，加工容易；经济性好；经济性好； 位于尾翼；位于尾翼；较大的弹

13、道修正能力较大的弹道修正能力.5.论文主要内容 2.修正机构设计1.天线座 2.GPS天线 3.弹尾 4.控制组件底座 5.惯性开关 6.DSP组件7.控制组件压盖 8.GPS接收机 9.电池2.1修正机构作用原理修正机构作用原理修正机构作用原理图 控制系统结构图.6.论文主要内容 2.修正机构设计整体效果图解锁前解锁后1.垫片 2.螺钉 3.阻尼片 4.弹簧 5.滑块.7.论文主要内容 2.修正机构设计2.2发射强度校核发射强度校核 基于修正弹发射的安全性，需进行发射强度校核，本文采用瞬态动力学数值方法进行强度校核。膛压曲线图 模型剖视图模型网格图弹丸施加载荷图 .8.论文主要内容 2.修正

14、机构设计总体等效应力云图 总体等效应变云图 可以看出，弹体最大等效应力未超过材料的屈服强度，弹体等效应变量也较小，所以弹丸满足强度要求。 .9.论文主要内容 2.修正机构设计2.3机构理论分析机构理论分析扭簧的工作转矩为： 4111023667E dMD n空气阻力矩为： 2221sin4xWMSLC阻尼片的转动角速度： 12dJMMdt2122MMdddtdtJ从而有： .10.论文主要内容 2.修正机构设计设y1=，y2=d/dt，则可得微分状态方程：12224112112012123sin3()644(0)0,(0)0 xwyyS v CyMME dyyJD nmLyy 00.010.0

15、20.030.040.050102030405060708090100时间（s）张开角度（）00.010.020.030.040.05050100150200250300350时间（s）转动角速度（rad/s） 张开角度随时间变化 转动角速度随时间变化.11.论文主要内容 2.修正机构设计2.4样机收星试验样机收星试验 修正机构样机及控制元件 样机收星试验过程 .12.论文主要内容 3.流场数值模拟分析3.1计算模型计算模型修正前弹丸表面网格 修正后弹丸表面网格 .13.论文主要内容 3.流场数值模拟分析3.2阻力系数曲线阻力系数曲线FLUENT软件中有关选项的设置为：软件中有关选项的设置为：

16、(1) 选择密度基求解器；(2) 流体为理想气体；(3) 来流压力为一个大气压；(4) 采用多重网格法，层数设置为5层；(5) 设置马赫数范围为0.62.0。.14.论文主要内容 3.流场数值模拟分析阻力系数计算结果：阻力系数计算结果：弹丸阻力系数曲线图 阻力系数曲线基本满足阻力系数随马赫数变化的规律。阻力系数曲线基本满足阻力系数随马赫数变化的规律。.15.论文主要内容 3.流场数值模拟分析阻力系数曲线拟合结果：阻力系数曲线拟合结果：阻力系数拟合曲线 状态状态马赫数马赫数/Ma拟合拟合方法方法参数结果参数结果最大最大误差误差修正前修正前Ma1.2Logistic0.1470.5448.77e7

17、-190.0281.2Ma2Logistic0.5750.1161.362e-68.1070.007修正后修正后Ma0.9Logistic1.1741.365873-8.50.0160.910347y=1.5x-3375(2450，300)18.38(2448.1，297.1)3.5284y=x-100(300，200)9.42(296.9，196.9)4.4182y=x-200(280，80)8.0(281，81)1.4高坡地形弹道修正仿真 .27.论文主要内容 5.基于地理信息一维弹道修正的研究低洼地形弹道修正仿真 装药号射角()地形函数目标点 (m)张开时刻(s)落点 (m)误差(m)6

18、46y=-x+2900(3000，-100)14.54(2994，-94)8.5684y=-x+700(850，-150)33.74(848.6，-148.3)2545y=-4x+14000(3620，-480)/10585y=-x+440(620，-180)30.52(616.8，-176.8)4.5448y=-4x+10560(2800，-640)17.96(2799.8，-639.3)0.7484y=-5x+2300(600，-700)24.68(598.2，-691)9.1383y=-3x+900(480，-540)13.04(482.8，-548.4)8.9284y=-x+100(3

19、00，-200)9.12(297，-197)4.3146y=-x+1100(1200，-100)10(1200.6，-100.6)0.8 确定不同地形弹道修正阻尼片张开时刻的算法总体有效，基于地形匹配进行弹道修正比水平面上修正的误差要大，这是因为基于地形匹配弹道修正的误差是水平射程和竖直高度上两方面误差的累积。 .28.论文主要内容 5.基于地理信息一维弹道修正的研究5.3精度影响因素分析精度影响因素分析45射角 60射角 高坡地形低洼地形45射角 60射角 目标高度因素的影响：目标高度因素的影响： 高坡地形弹道修正落点误差随目标点高度的增高呈不确定变化趋势；低洼地形其落点误差随目标点高度的降低(图中高度绝对值增大)而增大。 .29.论文主要内容 5.基于地理信息一维弹道修正的研究5.3精度影响因素分析精度影响因素分析45射角 60射角 高坡地形低洼地形45射角 60射角 地形坡度因素的影响：地形坡度因素的影响： 高坡地形和低高坡地形和低洼地形下，弹道修洼地形下，弹道修正误差都随坡度斜正误差都随坡度斜率的增大而增大。率的增大而增大。这是因为随着坡度这是因为随着坡度斜率的增大，阻尼斜率的增大，阻尼片张开时刻调节过片张开时刻调节过程中，弹丸落点在程中，弹丸落点在高度上调节的速度高度上调节的速度增大，导致弹道修