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

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

2、刀具是刚体且锋利，只考虑刀具的温度传导,不考虑刀具、工件的振动,假 设 条 件,3,由于刀具和工件的切削厚度方向上，切削工程中层厚不变，所以按平面应变来模拟,论文主要内容,1.铝合金A357切削加工有限元模型,4,1.1材料模型 1.1.1 A357的Johnson-Cook本构模型 材料本构模型用来描述材料的力学性质，表征材料变形过程中的动态响应。在材料微观组织结构一定的情况下，流动应力受到变形程度、变形速度、及变形温度等因素的影响非常显著。这些因素的任何变化都会引起流动应力较大的变动。因此材料本构模型一般表示为流动应力与应变、应变率、温度等变形参数之间的数学函数关系。建立材料本构模型，无论

3、是在制定合理的加工工艺方面，还是在金属塑性变形理论的研究方面都是极其重要的。在以塑性有限元为代表的现代塑性加工力学中，材料的流动应力作为输入时的重要参数，其精确度也是提高理论分析可靠度的关键。在本课题研究中，材料本构模型是切削加工数值模拟的必要前提，是预测零件铣削加工变形的重要基础，只有建立了大变形情况下随应变率和温度变化的应力应变关系，才能够准确描述材料在切削加工过程的塑性变形规律，继而才能在确定的边界条件和切削载荷下预测零件的变形大小及趋势,论文主要内容,1.铝合金A357切削加工有限元模型,4,1.1.1 A357的Johnson-Cook本构模型 在切削过程中，工件在高温、大应变下发生

4、弹塑性变形，被切削材料在刀具的作用下变成切屑时的时间很短，而且被切削层中各处的应变、应变速率和温度并不均匀分布且梯度变化很大。因此能反映出应变、应变速率、温度对材料的流动应力影响的本构方程，在切削仿真中极其关键。当前常用的塑性材料本构模型主要有：Bodner-Paton、Follansbee-Kocks、Johnson-Cook、 Zerrilli-Armstrong等模型，而只有Johnson-Cook模型描述材料高应变速率下热粘塑性变形行为。JohnsonCook模型认为材料在高应变速率下表现为应变硬化、应变速率硬化和热软化效应，JohnsonCook模型如下所示,式中第一项描述了材料的应

5、变强化效应，第二项反映了流动应力随对数应变速率增加的关系，第三项反映了流动应力随温度升高指数降低的关系。 、Tr分别表示参考应变速率和参考温度，Tm为材料熔点。式中A、B、n、C、m、D、k是7 个待定参数；A、B、n表征材料应变强化项系数；C表征材料应变速率强化项系数；m表征材料热软化系数； , 分别为常温材料熔点,论文主要内容,1.铝合金A357切削加工有限元模型,4,1.1.2 材料失效准则 实现切屑从工件分离，本文采用的是剪切失效模型。剪切失效模型是基于等效塑性应变在积分点的值，当损伤参数达到1时，单元即失效，失效参数定义如下： 式中： 为失效参数， 为等效塑性应变初始值， 为等效塑性

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

7、应力，q为Mises应力),温度，预定义域变量。这里采用JohnsonCook模型定义失效应变,论文主要内容,1.铝合金A357切削加工有限元模型,4,1.1.2 材料失效准则,式中： 为低于转变温度的条件下测得的实效常数。 为参考应变率， 为塑性应变率。 由下式确定,是当前温度， 是熔点， 是室温,论文主要内容,1.铝合金A357切削加工有限元模型,4,1.1.2 材料失效准则,图1.1中 和 为材料开始损伤时的屈服应力和等效塑性应变。 是材料失效时即图中D=1时的等效塑性应变。材料失效时的等效塑性应变 依赖于单元的特征长度，不能作为描述材料损伤演化的准则。相反，材料损伤演化的准则又等效塑性

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

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

10、=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 摩擦模型 金属切削过程中，刀具前刀面的摩擦状态非常复杂，通常把前刀面得摩擦区分为粘结区和滑动区，粘结区的摩擦状态与材料的临界剪应力有关，滑动区可近似认为摩擦系数为常值可以用下式来表示： 式中：

11、为接触面的滑动剪切应力； 为摩擦系数； 为接触面上的压力； 为材料的临界屈服压力,论文主要内容,1.铝合金A357切削加工有限元模型,4,1.2.1 质量放大,论文主要内容,1.研究意义,修正引信,零件加工难度大,阻力环面积小，弹道修正能力受限制,修正模块微型化技术限制,迫弹修正机构,结构简单，加工容易,经济性好,位于尾翼,较大的弹道修正能力,5,论文主要内容,2.修正机构设计,1.天线座 2.GPS天线 3.弹尾 4.控制组件底座 5.惯性开关 6.DSP组件 7.控制组件压盖 8.GPS接收机 9.电池,2.1修正机构作用原理,修正机构作用原理图,控制系统结构图,6,论文主要内容,2.修正

12、机构设计,整体效果图,解锁前,解锁后,1.垫片 2.螺钉 3.阻尼片 4.弹簧 5.滑块,7,论文主要内容,2.修正机构设计,2.2发射强度校核,基于修正弹发射的安全性，需进行发射强度校核，本文采用瞬态动力学数值方法进行强度校核,膛压曲线图,模型剖视图,模型网格图,弹丸施加载荷图,8,论文主要内容,2.修正机构设计,总体等效应力云图,总体等效应变云图,可以看出，弹体最大等效应力未超过材料的屈服强度，弹体等效应变量也较小，所以弹丸满足强度要求,9,论文主要内容,2.修正机构设计,2.3机构理论分析,扭簧的工作转矩为,空气阻力矩为,阻尼片的转动角速度,从而有,10,论文主要内容,2.修正机构设计,

13、设y1=，y2=d/dt，则可得微分状态方程,张开角度随时间变化 转动角速度随时间变化,11,论文主要内容,2.修正机构设计,2.4样机收星试验,修正机构样机及控制元件,样机收星试验过程,12,论文主要内容,3.流场数值模拟分析,3.1计算模型,修正前弹丸表面网格,修正后弹丸表面网格,13,论文主要内容,3.流场数值模拟分析,3.2阻力系数曲线,FLUENT软件中有关选项的设置为,1) 选择密度基求解器,2) 流体为理想气体,3) 来流压力为一个大气压,4) 采用多重网格法，层数设置为5层,5) 设置马赫数范围为0.62.0,14,论文主要内容,3.流场数值模拟分析,阻力系数计算结果,弹丸阻力

14、系数曲线图,阻力系数曲线基本满足阻力系数随马赫数变化的规律,15,论文主要内容,3.流场数值模拟分析,阻力系数曲线拟合结果,阻力系数拟合曲线,16,论文主要内容,3.流场数值模拟分析,3.3气动特性分析,M=0.6速度等值线云图,M=0.6压力等值线云图,17,论文主要内容,3.流场数值模拟分析,M=0.9压力等值线云图,M=0.9速度等值线云图,18,论文主要内容,3.流场数值模拟分析,M=1.2速度等值线云图,M=1.2压力等值线云图,19,论文主要内容,4.一维弹道仿真分析,4.1迫击炮弹弹道分析,采用四阶龙格-库塔法求解，得,迫弹外弹道模型,迫弹不同射角弹道曲线,20,论文主要内容,4

15、.一维弹道仿真分析,速度变化曲线 弹道倾角变化曲线,弹丸在相同初速下，随着射角增大，弹丸的飞行时间越长，落地速度越大，落地时的弹道倾角增大,21,论文主要内容,4.一维弹道仿真分析,4.2弹道修正能力分析,阻尼片张开前后弹丸阻力系数,迫弹全装药极限修正弹道曲线,T为阻尼片张开时刻,22,论文主要内容,4.一维弹道仿真分析,4.3阻尼片张开时刻算法,阻尼片张开时刻算法流程图,23,论文主要内容,4.一维弹道仿真分析,迫弹弹道修正仿真结果,24,论文主要内容,4.一维弹道仿真分析,阻尼片延时展开对修正精度的影响,弹道上升段 弹道下降段,25,论文主要内容,5.基于地理信息一维弹道修正的研究,高坡地

16、形 低洼地形,5.1不同地形弹道修正分析,不同地形弹道修正,26,论文主要内容,5.基于地理信息一维弹道修正的研究,5.2程序实现,基于不同地形进行弹道修正，其确定阻尼片张开时刻的算法,1)先采用水平面上进行弹道修正确定阻尼片张开时刻的算法求解出水平距离满足修正精度的张开时刻作为初始值； (2) 代入弹道方程，求解地形曲线上的弹道落点，若满足精度，则结束计算，否则采用在初始值附近一定时间范围内对时间变量进行循环，直至满足精度要求,高坡地形弹道修正仿真,27,论文主要内容,5.基于地理信息一维弹道修正的研究,低洼地形弹道修正仿真,确定不同地形弹道修正阻尼片张开时刻的算法总体有效，基于地形匹配进行

17、弹道修正比水平面上修正的误差要大，这是因为基于地形匹配弹道修正的误差是水平射程和竖直高度上两方面误差的累积,28,论文主要内容,5.基于地理信息一维弹道修正的研究,5.3精度影响因素分析,45射角 60射角,高坡地形,低洼地形,45射角 60射角,目标高度因素的影响,高坡地形弹道修正落点误差随目标点高度的增高呈不确定变化趋势；低洼地形其落点误差随目标点高度的降低(图中高度绝对值增大)而增大,29,论文主要内容,5.基于地理信息一维弹道修正的研究,5.3精度影响因素分析,45射角 60射角,高坡地形,低洼地形,45射角 60射角,地形坡度因素的影响,高坡地形和低洼地形下，弹道修正误差都随坡度斜率的增大而增大。这是因为随着坡度斜率的增大，阻尼片张开时刻调节过程中，弹丸落点在高度上调节的速度增大，导致弹道修正关于竖直高度这一因素的误差增大,30,论文主要内容,6.总结与进一步研究方向,6.1内容总结,1）对迫击炮弹一维弹道修正增阻机构进行了设计,2）进行了合理性分析，如发射强度校核、收星试验等,