CN110597180B - 一种基于切削过程仿真的单晶铜微铣削力预测方法 - Google Patents
一种基于切削过程仿真的单晶铜微铣削力预测方法 Download PDFInfo
- Publication number
- CN110597180B CN110597180B CN201910850915.4A CN201910850915A CN110597180B CN 110597180 B CN110597180 B CN 110597180B CN 201910850915 A CN201910850915 A CN 201910850915A CN 110597180 B CN110597180 B CN 110597180B
- Authority
- CN
- China
- Prior art keywords
- slip
- micro
- tensor
- shear
- single crystal
- 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.)
- Expired - Fee Related
Links
- 239000013078 crystal Substances 0.000 title claims abstract description 145
- 238000003801 milling Methods 0.000 title claims abstract description 98
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 title claims abstract description 73
- 229910052802 copper Inorganic materials 0.000 title claims abstract description 71
- 239000010949 copper Substances 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000004088 simulation Methods 0.000 title claims abstract description 33
- 238000005520 cutting process Methods 0.000 title claims abstract description 30
- 239000000463 material Substances 0.000 claims abstract description 39
- 238000012545 processing Methods 0.000 claims abstract description 11
- 239000013598 vector Substances 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 11
- 230000001186 cumulative effect Effects 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000005516 engineering process Methods 0.000 claims description 8
- 230000005489 elastic deformation Effects 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 6
- 238000012217 deletion Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 4
- 230000037430 deletion Effects 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 239000002178 crystalline material Substances 0.000 claims description 2
- 238000009795 derivation Methods 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 238000009826 distribution Methods 0.000 claims description 2
- 230000008571 general function Effects 0.000 claims description 2
- 230000003993 interaction Effects 0.000 claims description 2
- 230000002452 interceptive effect Effects 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 claims description 2
- 238000005728 strengthening Methods 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims 1
- 238000003754 machining Methods 0.000 abstract description 5
- 238000005094 computer simulation Methods 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 4
- 241000755266 Kathetostoma giganteum Species 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- UONOETXJSWQNOL-UHFFFAOYSA-N tungsten carbide Chemical compound [W+]#[C-] UONOETXJSWQNOL-UHFFFAOYSA-N 0.000 description 1
- 239000011882 ultra-fine particle Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/401—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by control arrangements for measuring, e.g. calibration and initialisation, measuring workpiece for machining purposes
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/34—Director, elements to supervisory
- G05B2219/34242—For measurement only
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明属于微小零件精密高效加工领域，特别涉及一种基于切削过程仿真的单晶铜微铣削力预测方法。由于单晶铜材料的特殊性，本发明首先综合考虑了晶体变形运动学、晶体塑性本构关系，率相关晶体的硬化规律等因素对于微铣削单晶铜的影响；然后对本构方程的增量形式进行了推导，采用Fortran语言编写了单晶铜VUMAT用户材料子程序，将晶体塑性本构引入到有限元仿真中；最后采用ABAQUS软件对单晶铜微铣削过程进行仿真建模，实现了单晶铜微铣削力预测。本发明以有限元切削仿真为基础，实现单晶铜微铣削力的准确预测，为研究单晶铜微铣削过程提供技术支撑，提高单晶铜微铣削加工精度与效率，具有实际应用价值。
Description
技术领域
本发明属于微小零件精密高效加工领域，特别涉及一种基于切削过程仿真的单晶铜微铣削力预测方法。
技术背景
随着科学技术的进步，航空航天等领域微小结构/零件除了要满足加工工艺上的要求，还要满足一些特殊要求，如需要有优质的导电性能。单晶铜由于具有一致的结晶取向，且只有一个晶粒，内部没有晶界，其导电性能超过一般的金属且具有良好的热、疲劳和蠕变性能，使其在航空航天等领域应用广泛。微铣削技术是加工单晶铜微小零件的有效手段。研究微小零件微铣削技术，对提高其加工精度及效率具有重要意义。微铣削过程中，微铣削力作为一个重要过程物理参数，研究微铣削力建模方法对于优化微小零件微铣削加工工艺、提高加工质量具有重要指导作用。
目前，微铣削力建模方法主要有力学解析法、有限元仿真法及智能算法建模法等。力学解析法考虑了刀具参数、工件材料屈服强度和切削条件等因素的影响，但建模过程中做出的假设及简化会使得预测精度降低且建模过程复杂。有限元仿真法节省成本，可以综合考虑刀具及工件材料的影响，但计算成本高，预测精度差。智能算法建模法需要大量试验样本，且无法考虑实际加工过程情况，导致微铣削力预测精度低。尽管微铣削力建模方法较多，但由于单晶铜材料的各向异性，不同晶体取向的单晶铜微铣削加工机理也不同，目前还缺乏切削过程仿真的单晶铜微铣削力预测方法。
卢晓红等人在《一种基于刀具磨损效应的微铣削力建模方法》中提出了一种基于刀具磨损效应的微铣削力建模方法。应用DEFORM有限元软件建立微铣削三维仿真模型。采用Johnson-Cook材料本构模型表征材料属性，并将刀具磨损效应引入微铣削力建模过程，实现一种模型对于微铣削力的精准预测，提高模型鲁棒性。但文献并没有开展有关单晶铜这类各向异性材料的研究。
尽管微铣削力建模方法较多，但由于单晶铜材料的各向异性，不同晶体取向的单晶铜微铣削加工机理也不同，目前还缺乏切削过程仿真的单晶铜微铣削力预测方法。
发明内容
为解决上述问题，本发明提出一种基于切削过程仿真的单晶铜微铣削力预测方法。以有限元切削仿真为基础，实现对于单晶铜微铣削力的准确预测，为单晶铜微铣削过程提供技术支撑，提高单晶铜微铣削加工精度与效率，具有实际应用价值。
本发明采用的技术方案是：
一种基于切削过程仿真的单晶铜微铣削力预测方法，由于单晶铜材料的特殊性，首先综合考虑了晶体变形运动学、晶体塑性本构关系、率相关晶体的硬化规律等因素对于微铣削单晶铜的影响；然后对本构方程的增量形式进行推导，采用Fortran语言编写单晶铜VUMAT用户材料子程序，将晶体塑性本构引入到有限元仿真中；最后采用ABAQUS软件对单晶铜微铣削过程进行仿真建模，实现单晶铜微铣削力预测，方法的具体步骤如下：
步骤1：构建晶体变形运动学模型
晶体塑性理论认为晶体的变形是由晶格畸变和位错滑移组成。其中晶格畸变可以通过连续介质力学中的弹性部分来描述，位错滑移可以通过连续介质力学来表示。如图1所示，晶体的总变形梯度F被分解为：
F＝Fe·Fp (1)
其中，Fe为弹性变形梯度，Fp为塑性的变形梯度。
滑移方向：
m*(α)＝Fem(α) (2)
其中，m(α)和m*(α)分别是晶体的晶格发生畸变前的第α滑移系滑移方向的单位向量和晶体的晶格发生畸变后的第α滑移系滑移方向的单位向量。
滑移面法向：
n*(α)＝(Fe -1)Tn(α) (3)
其中，n(α)和n*(α)分别是晶体的晶格发生畸变前的第α滑移面滑移方向的单位法向量和晶体的晶格发生畸变后的第α滑移面滑移方向的单位法向量。
当前状态的速度梯度张量L为：
其中，
和
分别为F、Fe和Fp的导数；Le和Lp分别是弹性速度梯度张量和塑性速度梯度张量。
建立各滑移系的滑移剪应变与晶体塑性变形的关系：
其中，
为在滑移系α上的滑移剪切应变率，对所有激活的滑移系进行求和。
另外，速度梯度张量L还可以分解为对称部分D和非对称部分W，D为变形率张量，W为旋率张量；变形率张量D和旋率张量W也可以按弹性部分和塑性部分，即：
L＝D+W，D＝De+Dp，W＝We+Wp (6)
其中，De为弹性变形率张量，Dp为塑性变形率张量，We为弹性旋转张量，Wp为塑性旋转张量。
并满足：
因此，变形率张量D与旋转张量W的表达式如下：
这些晶体变形学的公式，将晶体变形中的滑移剪切率与宏观上的变形率相联系起来。
步骤2：晶体塑性本构关系的建立
由于弹性势的存在，确定了晶格的弹性变形率张量De和柯西应力焦曼率
之间存在着如下的关系：
其中，σ为当前的柯西应力；Z表示单位二阶张量；L表示一组弹性模量张量，包括：Lijkl，Ljikl，Lijlk，Lklij；它们具有完全对称的特性；
为晶格旋转轴上的共旋转应力率，与材料旋转轴上的共旋转应力率
有如下的关系：
其中：
由此，建立了应力率、变形率和滑移剪切应变率的关系，而应变率的得出先求出各滑移系的剪切应变率，而它是根据硬化规律来计算的。
步骤3：率相关晶体的硬化规律分析
单晶体滑移的程度，取决于外加载荷所产生的切应力的大小、晶体结构的类型以及开动滑移面相对于切应力的取向。不同取向的单晶体产生滑移时所需的拉伸载荷是不同的。可以认为滑移开始时，这一数值与单晶体一般应力应变曲线上的屈服应力相当。
当滑移面上存在临界分切应力(CRSS)，当滑移方向的切应力达到临界值，就会出现滑移现象。临界分切应力是晶体的固有属性，由晶体构型、材料组成和温度决定，与外力无关。晶体在滑移系α上的滑移剪切应变率
由分解剪应力τα决定：
其中，为了简化
表达形式将
简化为
表示滑移系α的参考应变率，gα表示滑移系的强度，fα是一个无量纲通用函数，描述应变率与应力的依赖关系。
采用幂分布来描述多晶蠕变fα(x)：
fα(x)＝x|x|n-1 (14)
其中，n为率敏感指数；x为函数自变量。
由式(13)所示，还要得到滑移系的强度gα的值，其演化方程可由下式确定：
其中，
为gα的导数；hαβ为滑移硬化模量，表示滑移系中的滑移剪应变对滑移系所产生的硬化。当α＝β时，hαα为自硬化模量，即相同滑移系自身变形对自身的硬化；当α≠β时，hαβ为潜硬化模量，即不同滑移系间产生的硬化；
为晶体在滑移系β上的滑移剪切应变率。
用幂函数来描述自硬化模量hαα：
其中，h0为初始硬化模量，τ0为初始屈服应力，τs为剪切屈服强度，γ为滑移系累计剪切应变，其中γ表达如下：
其中，t是时间变量。
潜硬化模量为：
hαβ＝qh(γ)(α≠β) (18)
其中，q是常数。
晶体材料中三个强化阶段的滑移硬化模量：
hαβ＝qhαα(β≠α) (20)
其中，hs为易滑移时的硬化模量；γα为在滑移系α上的累计剪切应变；γβ为在滑移系β上的累计剪切应变；函数G与交互式硬化有关：
其中，γ0为滑移系间相互作用达到峰值强度后的滑移量，每个分量fαβ表示特定滑移作用的强度值。
步骤4：本构模型的增量形式推导
有限元计算采用增量法，因此首先应将本构关系写成增量形式。
根据率相关的切线系数方法，定义第α个滑移系上的剪切应变增量Δγ(α)为：
其中，
和
分别为t+Δt时刻的剪切应变量和t时刻的剪切应变量。
写成线性插值方式为：
其中，θ为积分参数，取值范围为[0,1]，当θ为零时，对应为简单欧拉时间积分格式，一般为了保证具有较高的精度和稳定性，推荐在[0.5,1]的范围内取；
和
分别为t+Δt时刻的剪切率和t时刻的剪切率。
滑移剪切应变率是关于分解剪应力和当前强度的函数，将公式(23)进行泰勒展开，如下：
由公式(22-24)可推知：
引入各滑移系的施密特因子
和张量
由晶体滑移的硬化方程，可知当前滑移系的强度增量Δgα为：
由此可知，剪应力的增量△τα为：
其中，
和
为各滑移系的施密特张量。
可得共回转应力增量Δσij为：
其中，
和
为各滑移系参考基向量上的位移分向量；σjk、σij和σik为各滑移系参考基向量上的柯西分应力；△εkk和Δεij各滑移系参考基向量上的应变增量。
由给定的应变增量得出剪切应变增量：
其中，δαβ为克罗内克尔记号。
求解方程组(30)可解的该增量步中的滑移剪切应力的增量，带入式(27)、(28)、(29)，其他未知量均可求得。
步骤5：ABAQUS/VUMAT用户材料子程序编写，输入步骤1-4得到的参数，构建晶体弹塑性本构模型
采用Fortran语言编写单晶铜VUMAT用户材料子程序，VUMAT用户材料子程序流程如图2所示。第一步，从ABAQUS中读取材料参数PROPS(n)，定义积分点序号I并且令I＝1。第二步，定义滑移系组数、积分系数和旋转增量；计算滑移方向滑移面法向、滑移变形张量、状态变量、剪切应变和累计剪切应变；计算滑移旋转张量、剪切应变率和其导数；求出自潜硬化模量；更新应力和应变；进而判断I是否是最后一个积分点序号，若不满足，则执行I＝I+1，循环第二步，若满足条件，则继续向下输出；最后，将更新好的应力和应变输出到ABAQUS中，结束程序。
步骤6：基于有限元方法获得单晶铜微铣削力预测
测绘微铣刀几何结构，建立微铣刀模型，导入ABAQUS软件中，将其设置为刚体。建立三维加工工件模型，设置为弹塑性体。对模型进行网格划分，刀具模型网格类型采用R3D3型网格，工件模型网格类型采用C3D8R，并将刀具及工件接触处切削区域的网格局部加密。
有限元仿真中，为了保证仿真结果的准确性，需要选择合适的摩擦模型来反映实际加工中的摩擦。微铣削加工过程中切屑与前刀面的接触区域有着两种摩擦类型，分别为粘结区和滑移区。
摩擦特性表达如下：
τf＝τs(τf≥τs) (31)
τf＝μσn(τf＜τs) (32)
其中，式(31)表示粘结区，式(32)表示滑动区，其中τf是摩擦剪切力，τs是工件的剪切屈服强度，μ是摩擦系数，σn是接触区域的正应力。根据摩擦剪切力与材料剪切屈服强度的比较，判断出当前处于的摩擦区。
采用单元删除的建模技术模拟材料去除。为实现单元删除，需要设置材料去除的准则。采用一种基于剪切应变的标准来解释材料去除过程中的各向异性：
max(γ-γcr,γsl,min-γsl,cr)≥0，γsl,min＝min(γα)，α＝1,2...N(33)
其中，γsl,cr和γcr分别是单滑移系统上剪切应变的临界值和所有滑移系统上的累积剪切应变，监测单个滑移系统上的剪切以及由于滑移系统引起的整体滑移；如果达到单个系统的临界值或累积滑移的临界值就删除该单元；γsl，min为滑移系α上累计剪切应变最小值。
设置主轴转速、每齿进给量、轴向切深等参数，即可对单晶铜微铣削力进行预测。
单晶铜微铣削仿真过程如图3所示，该参数下ABAQUS软件输出的三向切削力如图4所示，图中结果为采用ABAQUS后处理软件滤波之后的结果。
本发明的有益效果：本发明所述的基于切削过程仿真的单晶铜微铣削力预测方法无需大量实验而且对单晶铜这种各向异性材料适应性强，可实现单晶铜微铣削力预测，提高微铣削力的预测精度和效率，具有实际应用价值。
附图说明
图1是晶体变形运动示意图。
图2是VUMAT子程序流程图。
图3是单晶铜微铣削仿真过程。
图4是ABAQUS软件输出的切削力。
图5(a)是<100>晶向单晶铜微铣削Fx方向仿真值与实验值的对比图。
图5(b)是<100>晶向单晶铜微铣削Fy方向仿真值与实验值的对比图。
图5(c)是<100>晶向单晶铜微铣削Fz方向仿真值与实验值的对比图。
图6(a)是<110>晶向单晶铜微铣削Fx方向仿真值与实验值的对比图。
图6(b)是<110>晶向单晶铜微铣削Fy方向仿真值与实验值的对比图。
图6(c)是<110>晶向单晶铜微铣削Fz方向仿真值与实验值的对比图。
图7(a)是<111>晶向单晶铜微铣削Fx方向仿真值与实验值的对比图。
图7(b)是<111>晶向单晶铜微铣削Fy方向仿真值与实验值的对比图。
图7(c)是<111>晶向单晶铜微铣削Fz方向仿真值与实验值的对比图。
具体实施方式
下面结合附图和技术方案详细说明本发明的具体实施方式。
考虑到单晶铜微小零件微铣削过程中，微铣削力是一个重要的过程参量，因而研究单晶铜微铣削力建模方法对优化微小零件微铣削加工工艺、提高加工质量具有重要指导作用。此外，由于单晶铜材料的各向异性对于微铣削过程中的切削力变化有一定的影响。据此，针对单晶铜微小零件微铣削预测力难题，发明了一种基于切削过程仿真的单晶铜微铣削力预测方法。
以两刃平头铣刀微铣削直槽为例，采用的刀具是CrN涂层的超微粒子碳化钨两刃平头立铣刀，刀具的直径为0.4mm、刀具刃长为1.2mm、刃口圆弧半径1.9μm螺旋角20°、周刃前角9°、周刃后角12°，借助ABAQUS软件进行仿真，详细说明本发明实施过程。
首先确定选用的工件材料为单晶铜，其材料参数性能为密度：8,960kg/m3；弹性常数C11:168,150Mpa；弹性常数C12:121,400Mpa；弹性常数C44:75,400Mpa；比热容8,960J/kg·℃；热膨胀系数0.000016/℃；热传导系数：385W/m·K；熔化温度1,356K。
将通过扫描电镜得到的图片导入到AutoCAD中进行精确描摹，得到其底面轮廓，将之导入Creo软件中，建立微铣刀的三维几何模型，采用有限元分析前处理软件Hypermesh对微铣刀模型进行网格划分，网格类型采用R3D3型网格，为了保证微铣刀的形状精度，在微铣刀的切削刃部分划分了较为密集的网格，刀具网格划分了23967个，最后导入至ABAQUS软件中设置为刚体。模型工件为阶梯状半环形体，设置为弹塑性材料。网格类型采用C3D8R，在与刀具接触的切削区域使用网格局部加密，网格数量为180540，网格最小尺寸为1.5μm。利用Fortran语言编写单晶铜VUMAT用户材料子程序，将晶体塑性本构引入到有限元仿真中。定义材料参数时带入上述材料参数以及率敏感指数n＝50，初始硬化模量h0＝180MPa，易滑移时的硬化模量hs＝24MPa，剪切屈服强度τs＝52MPa，初始屈服应力τ0＝4MPa，常数q＝1。
ABAQUS软件中：在Property模块中，剪切摩擦系数设置为0.48，设置材料去除的准则，公式(33)中单滑移系统上剪切应变的临界值γsl,cr值为6.0，所有滑移系统上的累积剪切应变γcr的值为0.068。在Assembly模块中，调入工件和刀具模型。在模型中，通过调整铣刀和工件的相对位置确定切削深度和进给距离。定义边界条件，将工件侧面及底面自由度严格约束。在Step模块中，建立初始步后依次插入微铣削加工分析步、退刀分析步、约束转换分析步，过程类型选Dynamics，Explicit，场变量输出中勾选力，以便读取切削力值。选用了三种单晶铜(<100>,<110>,<111>)，上表面分别是(100)、(110)和(111)晶面，进给方向分别为<001>，<001>，<110>。在Job模块中创建数据，检查任务，设置刀具主轴转速为40,000rpm；设置每齿进给量为3μm/z；轴向切削深度为40μm。数据检查无误后，递交任务，进行有限元分析。在仿真完成后，输出X、Y、Z三个方向的切削力。
表1<100>晶向单晶铜微铣削实验与仿真的切削力对比
从图5(a)，5(b)，5(c)中可以看出，<100>晶向单晶铜仿真得到的微铣削力与实验得到的微铣削力基本吻合，但仍然存在偏差，如表1所示为<100>晶向单晶铜微铣削力的实验值和仿真值的对比，从表中可以看出最大相对误差为28.6％，平均相对误差为14.8％。因此验证了所建立的<100>晶向的单晶铜微铣削仿真模型的准确性。
表2<110>晶向单晶铜微铣削实验与仿真的切削力对比
如图6(a)，6(b)，6(c)所示，<110>晶向单晶铜在进给方向力Fx，径向力Fy和轴向力Fz三个方向力的仿真值与实验值对比。结果表明，仿真得到的微铣削力与实验得到的微铣削力基本吻合，但仍然存在偏差，如表2所示为<110>晶向单晶铜微铣削力的实验值和仿真值的对比，从表中可以看出最大相对误差为29.4％，平均相对误差为19.5％。因此验证了所建立的<110>晶向的单晶铜微铣削仿真模型的准确性。
表3<111>晶向单晶铜微铣削实验与仿真的切削力对比
如图7(a)，7(b)，7(c)所示，<111>晶向单晶铜在进给方向力Fx，径向力Fy和轴向力Fz三个方向力的仿真值与实验值对比。结果表明，仿真得到的微铣削力与实验得到的微铣削力基本吻合，但仍然存在偏差，如表3所示为<111>晶向单晶铜微铣削力的实验值和仿真值的对比，从表中可以看出最大相对误差为28.0％，平均相对误差为15.8％。因此验证了所建立的<111>晶向的单晶铜微铣削仿真模型的准确性。
综上，验证了所建立的<100>，<110>，<111>晶向的单晶铜微铣削仿真模型的准确性。
Claims (1)
1.一种基于切削过程仿真的单晶铜微铣削力预测方法，其特征在于，具体步骤如下：
步骤1：构建晶体变形运动学模型
晶体的总变形梯度F被分解为：
F＝Fe·Fp (1)
其中，Fe为弹性变形梯度，Fp为塑性的变形梯度；
滑移方向：
m*(α)＝Fem(α) (2)
其中，m(α)和m*(α)分别是晶体的晶格发生畸变前的第α滑移系滑移方向的单位向量和晶体的晶格发生畸变后的第α滑移系滑移方向的单位向量；
滑移面法向：
n*(α)＝(Fe -1)Tn(α) (3)
其中，n(α)和n*(α)分别是晶体的晶格发生畸变前的第α滑移面滑移方向的单位法向量和晶体的晶格发生畸变后的第α滑移面滑移方向的单位法向量；
当前状态的速度梯度张量L为：
其中，
和
分别为F、Fe和Fp的导数；Le和Lp分别是弹性速度梯度张量和塑性速度梯度张量；
建立各滑移系的滑移剪应变与晶体塑性变形的关系：
其中，
为在滑移系α上的滑移剪切应变率，对所有激活的滑移系进行求和；
另外，速度梯度张量L可分解为对称部分D和非对称部分W，D为变形率张量，W为旋率张量；变形率张量D和旋率张量W可按弹性部分和塑性部分，即：
L＝D+W，D＝De+Dp，W＝We+Wp (6)
其中，De为弹性变形率张量，Dp为塑性变形率张量，We为弹性旋转张量，Wp为塑性旋转张量；
并满足：
因此，变形率张量D与旋转张量W的表达式如下：
步骤2：晶体塑性本构关系的建立
由于弹性势的存在，确定晶格的弹性变形率张量De和柯西应力焦曼率
之间存在着如下的关系：
其中，σ为当前的柯西应力；Z表示单位二阶张量；L表示一组弹性模量张量，包括：Lijkl，Ljikl，Lijlk，Lklij，具有完全对称的特性；
为晶格旋转轴上的共旋转应力率，与材料旋转轴上的共旋转应力率
有如下的关系：
其中：
步骤3：率相关晶体的硬化规律分析
晶体在滑移系α上的滑移剪切应变率
由分解剪应力τα决定：
其中，将
简化为
表示滑移系α的参考应变率，gα表示滑移系的强度，fα是一个无量纲通用函数，描述应变率与应力的依赖关系；
采用幂分布来描述多晶蠕变fα(x)：
fα(x)＝x|x|n-1 (14)
其中，n为率敏感指数；x为函数自变量；
滑移系的强度gα，其演化方程由下式确定：
其中，
为gα的导数；hαβ为滑移硬化模量，表示滑移系中的滑移剪应变对滑移系所产生的硬化；当α＝β时，hαα为自硬化模量，即相同滑移系自身变形对自身的硬化；当α≠β时，hαβ为潜硬化模量，即不同滑移系间产生的硬化；
为晶体在滑移系β上的滑移剪切应变率；
用幂函数来描述自硬化模量hαα：
其中，h0为初始硬化模量，τ0为初始屈服应力，τs为剪切屈服强度，γ为滑移系累计剪切应变，其中γ表达如下：
其中，t是时间变量；
潜硬化模量为：
hαβ＝qh(γ) (18)
其中，α≠β，q是常数；
晶体材料中三个强化阶段的滑移硬化模量：
hαβ＝qhαα (20)
其中，β≠α，hs为易滑移时的硬化模量；γα为在滑移系α上的累计剪切应变；γβ为在滑移系β上的累计剪切应变；函数G与交互式硬化有关：
其中，γ0为滑移系间相互作用达到峰值强度后的滑移量，每个分量fαβ表示特定滑移作用的强度值；
步骤4：本构模型的增量形式推导
根据率相关的切线系数方法，定义第α个滑移系上的剪切应变增量Δγ(α)为：
其中，
和
分别为t+Δt时刻的剪切应变量和t时刻的剪切应变量；
写成线性插值方式为：
其中，θ为积分参数，取值范围为[0,1]；
和
分别为t+Δt时刻的剪切率和t时刻的剪切率；
滑移剪切应变率是关于分解剪应力和当前强度的函数，将公式(23)进行泰勒展开，如下：
由公式(22-24)可推知：
引入各滑移系的施密特因子
和张量
由晶体滑移的硬化方程，可知当前滑移系的强度增量Δgα为：
则，剪应力的增量Δτα为：
其中，
和
为各滑移系的施密特张量；
得到共回转应力增量Δσij为：
其中，
和
为各滑移系参考基向量上的位移分向量；σjk、σij和σik为各滑移系参考基向量上的柯西分应力；Δεkk和Δεij各滑移系参考基向量上的应变增量；
由给定的应变增量得出剪切应变增量：
其中，δαβ为克罗内克尔记号；
步骤5：ABAQUS/VUMAT用户材料子程序编写，输入步骤1-4得到的参数，构建晶体弹塑性本构模型
第一步，从ABAQUS中读取材料参数PROPS(n)，定义积分点序号I并且令I＝1；第二步，定义滑移系组数、积分系数和旋转增量；计算滑移方向滑移面法向、滑移变形张量、状态变量、剪切应变和累计剪切应变；计算滑移旋转张量、剪切应变率和其导数；求出自潜硬化模量；更新应力和应变；进而判断I是否是最后一个积分点序号，若不满足，则执行I＝I+1，循环第二步，若满足条件，则继续向下输出；最后，将更新好的应力和应变输出到ABAQUS中，结束程序；
步骤6：基于有限元方法获得单晶铜微铣削力预测
测绘微铣刀几何结构，建立微铣刀模型，导入ABAQUS软件中，将其设置为刚体；建立三维加工工件模型，设置为弹塑性体；对模型进行网格划分，刀具模型网格类型采用R3D3型网格，工件模型网格类型采用C3D8R，并将刀具及工件接触处切削区域的网格局部加密；
微铣削加工过程中切屑与前刀面的接触区域有着两种摩擦类型，分别为粘结区和滑移区；
摩擦特性表达如下：
τf＝τs (31)
其中，τf≥τs，式(31)表示粘结区；
τf＝μσn (32)
其中，τf＜τs，式(32)表示滑动区，其中τf是摩擦剪切力，τs是工件的剪切屈服强度，μ是摩擦系数，σn是接触区域的正应力；
根据摩擦剪切力与材料剪切屈服强度的比较，判断出当前处于的摩擦区；
采用单元删除的建模技术模拟材料去除；为实现单元删除，需要设置材料去除的准则，采用一种基于剪切应变的标准来解释材料去除过程中的各向异性：
max(γ-γcr,γsl,min-γsl,cr)≥0，γsl,min＝min(γα)，α＝1,2...N (33)
其中，γsl,cr和γcr分别是单滑移系统上剪切应变的临界值和所有滑移系统上的累积剪切应变，监测单个滑移系统上的剪切以及由于滑移系统引起的整体滑移；如果达到单个系统的临界值或累积滑移的临界值就删除该单元；γsl，min为滑移系α上累计剪切应变最小值；
设置主轴转速、每齿进给量、轴向切深参数，即可对单晶铜微铣削力进行预测。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910850915.4A CN110597180B (zh) | 2019-09-10 | 2019-09-10 | 一种基于切削过程仿真的单晶铜微铣削力预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910850915.4A CN110597180B (zh) | 2019-09-10 | 2019-09-10 | 一种基于切削过程仿真的单晶铜微铣削力预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110597180A CN110597180A (zh) | 2019-12-20 |
| CN110597180B true CN110597180B (zh) | 2020-08-14 |
Family
ID=68858426
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910850915.4A Expired - Fee Related CN110597180B (zh) | 2019-09-10 | 2019-09-10 | 一种基于切削过程仿真的单晶铜微铣削力预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110597180B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111523268B (zh) * | 2020-04-22 | 2021-06-01 | 四川大学 | 一种基于机器学习的材料抗疲劳优化设计方法 |
| CN112683715B (zh) * | 2020-11-27 | 2022-03-29 | 华南理工大学 | 一种韧性金属材料临界切削条件的预测方法 |
| CN112989667B (zh) * | 2021-03-29 | 2022-04-19 | 山东大学 | 晶体材料零件三维晶体塑性有限元建模方法及系统 |
| CN114974478B (zh) * | 2022-06-10 | 2024-05-31 | 山东大学 | 考虑应变率的晶体金属材料直角微切削建模方法及系统 |
| CN115035964B (zh) * | 2022-06-10 | 2024-04-02 | 山东大学 | 晶体金属材料微铣削性能预测方法及系统 |
| CN115422797B (zh) * | 2022-08-26 | 2024-08-09 | 中国矿业大学 | 一种考虑应变率影响的大应变挤出切削带材织构预测方法 |
| CN118246190A (zh) * | 2023-03-21 | 2024-06-25 | 哈尔滨理工大学 | 一种双相不锈钢加工工艺 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101853324A (zh) * | 2010-06-24 | 2010-10-06 | 西北工业大学 | 圆周铣过程中铣削力建模方法 |
| US7823428B1 (en) * | 2006-10-23 | 2010-11-02 | Wright State University | Analytical method for use in optimizing dimensional quality in hot and cold rolling mills |
| CN102873381A (zh) * | 2012-09-29 | 2013-01-16 | 西安交通大学 | 一种基于动力学模型的高速铣削工艺参数优化方法 |
| CN103397255A (zh) * | 2013-08-09 | 2013-11-20 | 武汉钢铁(集团)公司 | 一种各向异性小的高性能易切削钢 |
| CN105117547A (zh) * | 2015-08-25 | 2015-12-02 | 大连理工大学 | 微铣削镍基高温合金预测残余应力的仿真方法 |
| CN105243195A (zh) * | 2015-09-16 | 2016-01-13 | 大连理工大学 | 一种微铣削镍基高温合金加工硬化的预测方法 |
| CN106156430A (zh) * | 2016-07-06 | 2016-11-23 | 大连理工大学 | 一种基于刀具磨损效应的微铣削力建模方法 |
| CN106407526A (zh) * | 2016-09-05 | 2017-02-15 | 大连理工大学 | 一种微铣削过程刀具后刀面磨损预测方法 |
| CN108416084A (zh) * | 2018-01-23 | 2018-08-17 | 南京理工大学 | 考虑复合材料弹塑性与损伤耦合的弹塑性损伤有限元算法 |
-
2019
- 2019-09-10 CN CN201910850915.4A patent/CN110597180B/zh not_active Expired - Fee Related
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7823428B1 (en) * | 2006-10-23 | 2010-11-02 | Wright State University | Analytical method for use in optimizing dimensional quality in hot and cold rolling mills |
| CN101853324A (zh) * | 2010-06-24 | 2010-10-06 | 西北工业大学 | 圆周铣过程中铣削力建模方法 |
| CN102873381A (zh) * | 2012-09-29 | 2013-01-16 | 西安交通大学 | 一种基于动力学模型的高速铣削工艺参数优化方法 |
| CN103397255A (zh) * | 2013-08-09 | 2013-11-20 | 武汉钢铁(集团)公司 | 一种各向异性小的高性能易切削钢 |
| CN105117547A (zh) * | 2015-08-25 | 2015-12-02 | 大连理工大学 | 微铣削镍基高温合金预测残余应力的仿真方法 |
| CN105243195A (zh) * | 2015-09-16 | 2016-01-13 | 大连理工大学 | 一种微铣削镍基高温合金加工硬化的预测方法 |
| CN106156430A (zh) * | 2016-07-06 | 2016-11-23 | 大连理工大学 | 一种基于刀具磨损效应的微铣削力建模方法 |
| CN106407526A (zh) * | 2016-09-05 | 2017-02-15 | 大连理工大学 | 一种微铣削过程刀具后刀面磨损预测方法 |
| CN108416084A (zh) * | 2018-01-23 | 2018-08-17 | 南京理工大学 | 考虑复合材料弹塑性与损伤耦合的弹塑性损伤有限元算法 |
Non-Patent Citations (3)
| Title |
|---|
| Cutting parameters optimization for MRR under the constraints of surface roughness and cutter breakage in micro-milling process;Xiaohong Lu;《Journal of Mechanical Science and Technology》;20180723;第3379-3388页 * |
| Inconel 718微铣削加工表面残余应力有限元仿真;卢晓红等;《东北大学学报(自然科学版)》;20170215(第02期);第254-259页 * |
| 精密超声振动切削单晶铜的计算机仿真研究;杨亮等;《系统仿真学报》;20070420(第04期);第738-741,753页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110597180A (zh) | 2019-12-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110597180B (zh) | 一种基于切削过程仿真的单晶铜微铣削力预测方法 | |
| Limido et al. | SPH method applied to high speed cutting modelling | |
| Uhlmann et al. | Cutting simulation with the meshfree finite pointset method | |
| Cheng et al. | Modeling and simulation of material removal rates and profile accuracy control in abrasive flow machining of the integrally bladed rotor blade and experimental perspectives | |
| Ayhan | Simulation of three-dimensional fatigue crack propagation using enriched finite elements | |
| Linjiang et al. | Inverse identification of constitutive parameters of Ti2AlNb intermetallic alloys based on cooperative particle swarm optimization | |
| Gu et al. | MPM simulations of high-speed and ultra high-speed machining of titanium alloy (Ti–6Al–4V) based on fracture energy approach | |
| Gylienė et al. | The numerical analysis of cutting forces in high feed face milling, assuming the milling tool geometry | |
| Muhamad et al. | Design improvement of steering knuckle component using shape optimization | |
| Constantin et al. | 3D FEM analysis of cutting processes | |
| Li et al. | Effect of cutting process adjustment on crystallographic texture of machined surface layer of titanium alloy | |
| Tang et al. | Finite-element simulation of conventional and high-speed peripheral milling of hardened mold steel | |
| Zong et al. | FEM optimization of tool geometry based on the machined near surface's residual stresses generated in diamond turning | |
| Ren et al. | Recent developments on thermo-mechanical simulations of ductile failure by meshfree method | |
| Hashmi et al. | A critical review of modeling and simulation techniques for loose abrasive based machining processes | |
| Marscher | A phenomenological model of abradable wear in high performance turbomachinery | |
| Kumar et al. | Finite element model based on abaqus/explicit to analyze the temperature effects of turning | |
| Lan et al. | Research on the influence of cutting parameters on chips in the process of diamond bead cutting pipeline steel | |
| Zhou et al. | Influence of cutting and clamping forces on machining distortion of diesel engine connecting rod | |
| Yang et al. | The simulation of cutting force and temperature field in turning of Inconel 718 | |
| Xue et al. | Finite element analysis of revolving tip-based cutting process | |
| Boubakri et al. | Micromachining simulation using a crystal plasticity model: ALE and CEL approaches | |
| Grissa et al. | A numerical–analytical approach to predict high cycle fatigue performance of finish machined AISI 316L steel | |
| Peng et al. | Analysis of micro/mesoscale sheet forming process by strain gradient plasticity and its characterization of tool feature size effects | |
| Childerhouse et al. | Rapid acquisition of digital fingerprints of Ti-6Al-4V macrotexture from machining force measurement data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200814 |