CN106334970B - 一种断续冲击载荷作用下的高速铣刀波动力学损伤识别方法 - Google Patents
一种断续冲击载荷作用下的高速铣刀波动力学损伤识别方法 Download PDFInfo
- Publication number
- CN106334970B CN106334970B CN201610966431.2A CN201610966431A CN106334970B CN 106334970 B CN106334970 B CN 106334970B CN 201610966431 A CN201610966431 A CN 201610966431A CN 106334970 B CN106334970 B CN 106334970B
- Authority
- CN
- China
- Prior art keywords
- mrow
- milling cutter
- msup
- mfrac
- wave
- 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
- 238000003801 milling Methods 0.000 title claims abstract description 172
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000005520 cutting process Methods 0.000 claims abstract description 54
- 230000006378 damage Effects 0.000 claims abstract description 41
- 239000002245 particle Substances 0.000 claims abstract description 37
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 238000010219 correlation analysis Methods 0.000 claims abstract description 5
- 239000000463 material Substances 0.000 claims description 41
- 238000000429 assembly Methods 0.000 claims description 23
- 230000000712 assembly Effects 0.000 claims description 23
- 230000035939 shock Effects 0.000 claims description 15
- 238000010276 construction Methods 0.000 claims description 13
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 5
- 230000005284 excitation Effects 0.000 claims description 4
- 230000000644 propagated effect Effects 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 238000009527 percussion Methods 0.000 claims description 2
- 230000006735 deficit Effects 0.000 abstract description 5
- 230000004069 differentiation Effects 0.000 abstract description 4
- 238000012360 testing method Methods 0.000 description 22
- 238000004458 analytical method Methods 0.000 description 10
- 230000008859 change Effects 0.000 description 7
- 239000013078 crystal Substances 0.000 description 6
- 230000001066 destructive effect Effects 0.000 description 5
- 206010061619 Deformity Diseases 0.000 description 4
- 230000001815 facial effect Effects 0.000 description 4
- 230000005012 migration Effects 0.000 description 4
- 238000013508 migration Methods 0.000 description 4
- 238000012098 association analyses Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000003776 cleavage reaction Methods 0.000 description 2
- 238000000329 molecular dynamics simulation Methods 0.000 description 2
- 230000007017 scission Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002149 energy-dispersive X-ray emission spectroscopy Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000008653 root damage Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q17/00—Arrangements for observing, indicating or measuring on machine tools
- B23Q17/09—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool
- B23Q17/0952—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool during machining
- B23Q17/0957—Detection of tool breakage
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q2717/00—Arrangements for indicating or measuring
- B23Q2717/006—Arrangements for indicating or measuring in milling machines
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
- Milling Processes (AREA)
Abstract
本发明公开了一种断续冲击载荷作用下的高速铣刀波动力学损伤识别方法,其技术要点在于由以下步骤构成:一、建立高速铣刀瞬态切削力模型和质点动力学方程;二、建立高速铣刀应力波的波动方程;三、求解高速铣刀应力波衰减应力值和衰减率;四、识别高速铣刀波动力学损伤。本发明可揭示出高速铣刀宏介观结构域之间的内在联系,识别出铣刀局部宏介观损伤,可确定高速铣刀的高频应力波汇集区域,并对其进行跨尺度关联分析,结合分子动力学仿真,识别出高速铣刀易发生损伤部位及介观尺度下损伤的类型,重点解决了高速、断续冲击载荷作用下铣刀组件潜在的损伤位置与损伤类型识别问题,同时也解决了高速铣刀波动力学损伤多尺度演变问题。
Description
技术领域:
本发明涉及一种高速铣刀损伤的识别方法,具体涉及一种断续冲击载荷作用下的高速铣刀波动力学损伤识别方法。
背景技术:
在高速、断续切削条件下,随着铣刀结构、材料和切削参数的改变,铣刀内部的能量和质量分布也随之发生改变,在获得高切削效率的同时,铣刀应力场分布的不均匀性显著增加,其内部能量不断积聚。在此条件下,铣刀一旦发生损伤或存在潜在缺陷,其破坏程度将迅速发展,若不能及时发现,产生的多是不可修复的严重事故。
在高速铣刀安全性研究中发现,采用铣刀强度失效判据和铣刀整体振动失稳判据,无法正确描述和评定铣刀损伤的产生和发展,铣刀最先发生的局部性损伤和破坏与作用在刀齿上的切削力冲击所引起的应力脉冲密切相关。分析发现,切削力冲击引起的应力脉冲,不会立即引起铣刀结构内所有质点的振动,而是按照组件介质的波速在铣刀内传播,形成具有不同速率过程的应力波。进一步研究距离切削力作用点不同位置处的应力波波形,发现铣刀内部的应力波具有衰减特性和在不同铣刀组件介质内传播时的折射、反射等波动特性。
该结果表明,在ISO15641标准和铣刀安全性判据的双重约束下,尽管切削力引起的应力波不会使高速铣刀产生整体性破坏,但因高速切削冲击产生的诸如刀片破损、刀体及联接元件结合面微小破坏等铣刀局部性的波动力学损伤现象却经常发生,使得铣刀实际失效寿命值远小于预先规定的最小失效寿命值。受此问题影响,相关企业高速加工设备的利用率仅为40%左右,波动力学损伤引发的高速铣刀局部破坏性失效问题,已成为制约高速铣削加工能效大幅度提高的瓶颈。基于此,提出了一种断续冲击载荷作用下的高速铣刀波动力学损伤识别方法是非常必要的。
已有的铣刀损伤识别的方法主要有损伤力学和分子动力学。连续损伤力学忽略了介观和微观上的控制机制,因而不能合理地处理分布式微损伤非平衡演化导致的失效问题。分子运动学能够刻画出微损伤形成细节,但微损伤演变过程中的微裂纹、微孔洞等特征,都很难在这种计算方法中体现出来,从而难以推导出宏观材料真实物理特征。因此,上述方法均无法单独用于解决高速铣刀波动力学损伤多尺度演变问题。
发明内容:
本发明针对波动力学损伤引发的高速铣刀局部破坏性失效问题,提供了一种断续冲击载荷作用下的高速铣刀波动力学损伤识别方法。通过本发明可确定高速铣刀的高频应力波汇集区域,并对其进行跨尺度关联分析,结合分子动力学仿真,识别出高速铣刀易发生损伤部位及介观尺度下损伤的类型,重点解决了高速、断续冲击载荷作用下铣刀组件潜在的损伤位置与损伤类型识别问题,同时也解决了高速铣刀波动力学损伤多尺度演变问题。
本发明的断续冲击作用下的高速铣刀波动力学损伤识别方法,为实现上述目的所采用的技术方案由以下步骤构成:
一、建立高速铣刀瞬态切削力模型和质点动力学方程
根据铣刀切削工艺参数、铣刀结构参数、铣刀组件材料参数,获取铣刀瞬态切削力波形满足式(1):
式中,p为作用于刀齿的单位切削层面积上的切削力,fz为每齿进给量,ap为轴向铣削深度,Fi(t)为第i个刀齿瞬态切削力,为随时间变化的刀齿瞬时位置角,m为参与切削的刀齿数量,F(t)为m个刀齿的瞬态切削力,t为时间,
利用铣刀质点运动微分方程表征其局部结构的动力响应,质点动力学方程满足式(2):
式中,u为组件内质点相对位移,T为弦内张力,ρ为弦线密度,x为波阵面随时间变化的空间位置,Fi(x,t)为刀齿瞬态切削力在该横向振动的激励;式(2)中,质点运动特征值为质点的位置u(x,t)和质点运动速度v的方程为:
二、建立高速铣刀应力波的波动方程
利用铣刀结构变化矩阵、铣刀组件材料本构关系、铣刀组件质点相对运动,构建铣刀组件内应力波的波动方程,
首先根据铣刀结构将其简化为轴对称结构,获得运动平衡方程为:
为了方便表明高速铣刀波动方程的推导过程,将材料本构关系定义为最简单的形式:
σ=Eε (5)
采用式(4)、式(5),并根据牛顿第二定律,建立高速铣刀波动方程为
c0 2=E/ρ0,ζ=-dA(x)/A(x)dx (7)
波速是描述铣刀应力波的重要参数满足式(8):
上述各式中:ρ0为材料密度,A(x)为刀具结构变化矩阵,v为刀具内的质点运动速度,σ为应力值;E为弹性模量,ε为应变;C为应力波的波速,x为波阵面随时间变化的空间位置;c0 2为材料内部应力波的传播速度,ζ为铣刀组件结构的修正系数;
三、求解高速铣刀应力波衰减应力值和衰减率
定义衰减率η为单位应力经过单位位移的衰减数值:
η=-(1/σx)dσx/dx (9)
对式(9)两边积分可得:
σx=σ0e-ηx (10)
应力波在铣刀内部经过多次反射后的波速C与初始波速c有如下关系:
τ=l/C=λ/c N=(C/c)(λ/l) (11)
高速面铣刀切削时承受的切削力冲击载荷所形成的冲击波长为冲击传递距离L的2倍,因此,上式可改写为下式:
N=2CL/(cl) (12)
对上式两边做偏微分变化,其中等式右端对L与l进行偏微分变化,得出下公式(13):
q=2ρ2C/(ρ1c+ρ2C) (13)
瞬态切削力载荷作用下,高速铣刀内部形成的应力波反射次数趋于无穷大的时候,铣刀刀齿区域内应力波衰减后的应力值σx为:
对上式进行整理得衰减率的求解方程为:
上式中,η为衰减率,σx为衰减后应力值,σ0为x=0处的应力值(Mpa),x为应力波传播距离;τ为冲击作用时间,λ为波长,c为应力波初始速度,N为反射次数,L为应力波传播距离,l为应力波在铣刀中传播的长度,q为应力波的反射效率,ρ1和ρ2分别为应力波传播时所经过的不同介质材料密度,B为应力传递比,σi为应力波单次反射后的应力,ηxAB为应力波单次反射的衰减率;
四、识别高速铣刀波动力学损伤
依据高速铣刀应力波的传播特性,确定应力波汇集区域,利用铣刀该区域跨尺度关联分析结果,结合分子动力学仿真,识别出铣刀局部宏观和介观损伤。
本发明的有益效果:本发明首先根据铣刀切削工艺参数、铣刀结构参数、铣刀组件材料参数来获取铣刀瞬态切削力波形,推导出反映铣刀局部结构响应的动力学方程,该方程的初始条件为质点位置与质点运动速度;再依据铣刀结构变化矩阵、铣刀组件材料本构关系和质点相对运动,构建铣刀应力波的波动方程,定量描述瞬态切削力激励下铣刀内部应力应变变化的波动性;之后推导出铣刀应力波衰减后的应力值和衰减率的求解方程,铣刀内部应力波的衰减会改变应力场分布状态,造成局部应力累积,进而引发细微小尺度损伤萌发;最后依据高速铣刀应力波的传播特性,确定应力波汇集区域,利用铣刀该区域跨尺度关联分析结果,结合分子动力学仿真,识别出铣刀局部宏介观损伤。
与现有铣刀损伤识别的方法相比,本发明可揭示出高速铣刀宏介观结构域之间的内在联系,有效解决了高速铣刀波动力学损伤多尺度演变问题,识别出铣刀局部宏介观损伤,为后续跨尺度设计更为高效和稳定的高能效铣刀提供理论依据。采用本发明的方法,可确定高速铣刀的高频应力波汇集区域,并对其进行跨尺度关联分析,结合分子动力学仿真,识别出高速铣刀易发生损伤部位及介观尺度下损伤的类型,重点解决了高速、断续冲击载荷作用下铣刀组件潜在的损伤位置与损伤类型识别问题。
附图说明:
图1为高速铣刀波动力学损伤识别方法的流程图;
图2为高速铣刀铣削过程图;
图3为多齿瞬态切削力波形图;
图4为霍普金森压杆实验中35CrMo的应力-应变关系的实验结果图;
图5为霍普金森压杆实验中42CrMo的应力-应变关系的实验结果图;
图6为冲击载荷0.4MPa下铣刀组件材料35CrMo的试件在超景深显微镜下的试件表面变形图;
图7为冲击载荷0.8MPa下铣刀组件材料35CrMo的试件在超景深显微镜下的试件表面变形图;
图8为冲击载荷0.4MPa下铣刀组件材料35CrMo的试件在金相显微镜下的试件表面变形图;
图9为冲击载荷0.8MPa下铣刀组件材料35CrMo的试件在金相显微镜下的试件表面变形图;
图10为冲击载荷0.4MPa下铣刀组件材料42CrMo的试件在超景深显微镜下的试件表面变形图;
图11为冲击载荷0.8MPa下铣刀组件材料42CrMo的试件在超景深显微镜下的试件表面变形图;
图12为冲击载荷0.4MPa下铣刀组件材料42CrMo的试件在金相显微镜下的试件表面变形图;
图13为冲击载荷0.8MPa下铣刀组件材料42CrMo的试件在金相显微镜下的试件表面变形图;
图14为42CrMo原子群运动图;
图15为晶格尺寸仿真结果图;
图16为原子群能量仿真结果图;
图17为高速铣刀损伤跨尺度分析流程图;
图18为刀体齿根损伤位置图;
图19为螺钉与刀片定位面损伤位置图;
图20为刀体齿根处损伤中位错运动的介观分析图;
图21为刀体齿根处损伤中位错塞积的介观分析图;
图22为刀体齿根处损伤中晶界迁移和微裂纹扩展的介观分析图;
图23为螺钉与刀片定位面损伤中位错运动的介观分析图;
图24为螺钉与刀片定位面损伤中位错塞积和微裂纹扩展的介观分析图;
图25为螺钉与刀片定位面损伤中晶界迁移和晶面解理的介观分析图。
具体实施方式:
以下结合各附图对本发明的具体实施方式进行说明:
参照图1,本发明的断续冲击载荷作用下的高速铣刀波动力学损伤识别方法由以下步骤构成:
一、建立高速铣刀瞬态切削力模型和质点动力学方程
高速铣刀切削时,由于其刀齿切削的断续性和切削厚度的变化,使得波形呈周期性变化的切削力成为引起刀具系统振动的主要激振力之一,其铣削过程及瞬态切削力波形如图2、图3所示,图2中,vf为铣刀进给速度,n为铣刀转速,kc,x、kc,y、kc,z分别为沿三个方向的模态刚度,cc,x、cc,y、cc,z分别为三个方向的模态阻尼,图3中,tm-2、tm-1、tm分别为第m-2、第m-1、第m个刀齿的切削时间。
多个刀齿同时参与切削时,高速铣刀瞬态切削力为:
式中,p为作用于刀齿的单位切削层面积上的切削力,fz为每齿进给量,ap为轴向铣削深度,Fi(t)为第i个刀齿瞬时切削力,φi(t)为随时间变化的刀齿瞬时位置角,m为参与切削的刀齿数量,F(t)为m个刀齿的瞬态切削力,t为时间。
高速铣刀瞬态切削力的冲击不仅引起铣刀整体振动,而且会引起铣刀局部振动,为揭示铣刀局部结构变形和损伤,利用铣刀质点运动微分方程表征其局部结构的动力响应。质点动力学方程满足式(2):
式中,u为组件内质点相对位移,T为弦内张力,ρ为弦线密度,x为波阵面随时间变化的空间位置,Fi(x,t)为刀齿瞬态切削力在该横向振动的激励;式(2)即为刀齿瞬态切削力激励下的铣刀质点受迫振动方程。其中,质点运动特征值为质点的位置u(x,t)和质点运动速度v
由此可见,利用铣刀质点运动特征值可以推导出反映铣刀局部结构响应的动力学方程,质点位置与质点运动速度为该动力学方程的初始条件。
二、建立高速铣刀应力波的波动方程
应力波作为铣刀内部应力载荷的一种扰动形式和质点的运动结果,有着特定特征参数对其描述,波动方程是研究应力波问题的关键控制方程,为此,考虑高速铣刀应力波的特点,对其波动方程进行简单推导,
根据铣刀结构将其简化为轴对称结构,运动平衡方程满足式(4):
将铣刀组件材料的本构关系定义为最简单的形式如式(5):
σ=Eε (5)
采用式(5)、式(6),根据牛顿第二定律,建立高速铣刀波动方程如式(6):
c0 2=E/ρ0,ζ=-dA(x)/A(x)dx (7)
波速是描述铣刀应力波的重要参数,在忽略介质转变效应的前提下,应力波的波速可以利用波阵面空间位置对时间的微商来获得,
上述各式中:ρ0为材料密度,A(x)为刀具结构变化矩阵,v为刀具内的质点运动速度,σ为应力值,E为弹性模量,ε为应变,u为组件内质点的相对位移,C为应力波的波速,x为波阵面随时间变化的空间位置;c0 2为材料内部应力波的传播速度,ζ为铣刀组件结构的修正系数。
三、求解高速铣刀应力波衰减率和衰减应力值
定义衰减率η为单位应力经过单位位移的衰减数值如式(9):
η=-(1/σx)dσx/dx (9)
对式(9)两边积分可得式(10):
σx=σ0e-ηx (10)
将应力波近似为矩形波,则应力波在铣刀内部经过多次反射次后的波速C与初始波速c关系如式(11):
τ=l/C=λ/c N=(C/c)(λ/l) (11)
高速面铣刀切削时承受的切削力冲击载荷所形成的冲击波长为冲击传递距离L的2倍,因此,式(11)可改写为式(12):
N=2CL/(cl) (12)
对式(12)两边做偏微分变化,其中等式右端对L与l进行偏微分变化,如式(13):
q=2ρ2C/(ρ1c+ρ2C) (13)
瞬态切削力载荷作用下,高速铣刀内部形成的应力波反射次数趋于无穷大的时候,铣刀刀齿区域内应力波衰减后的应力值σx为:
对式(14)整理得衰减率的求解方程为:
上述各式中,η为衰减率,σx为衰减后应力值,σ0为x=0处的应力值(Mpa),x为应力波传播距离,τ为冲击作用时间,λ为波长,c初始速度,N为反射次数,L为应力波传播距离,l为应力波在铣刀中传播的长度,q为应力波的反射效率,ρ1和ρ2分别为应力波传播时所经过的不同介质材料密度,B为应力传递比,σi为应力波单次反射后的应力,ηxAB为应力波单次反射的衰减率。
四、识别高速铣刀波动力学损伤
选取高速铣刀组件材料35CrMo、42CrMo,进行霍普金森压杆实验,获得铣刀组件材料在高速冲击载荷作用下的应力-应变关系和试件变形检测结果,如图4至图13所示。
实验结果表明,冲击载荷加载后试件呈波状界面,这是冲击加载后的铣刀组件材料的基本组织特征,高速、断续冲击载荷作用并没有使铣刀组件材料发生大规模的结构变化,但在一定程度上引起其内部结构变形和合并。
提出高速铣刀介观损伤分析方法,首先利用铣刀组件能谱分析结果,构建高速铣刀刀体原子群模型,然后采用第一性原理求解原子群能量最低的稳定构型,利用力连接方法进行铣刀组件分子动力学仿真,结果如图14至图16所示。图中,铣刀组件材料为42CrMo,加载应力0.97GPa。
由图14至图16可知,铣刀组件材料42CrMo在加载应力作用下,晶格结构尺寸呈现出波动性的变化,总体能量也呈现出明显的周期性波动。原子群构型受力中央有较为明显的位错带出现,发生了的位错滑移与位错攀移运动;在受力部位的两端交界带处由于位错密度的不断增加,位错间的相互作用力增大,发生严重的位错塞积,引起了应力集中。
依据上述分析结果,高速铣刀波动力学损伤识别首先根据铣刀切削工艺参数、铣刀结构参数、铣刀组件材料参数,获取铣刀瞬态切削力波形,推导出反映铣刀局部结构响应的动力学方程;依据铣刀结构变化矩阵、铣刀组件材料本构关系和质点相对运动,构建铣刀应力波的波动方程,定量描述瞬态切削力激励下铣刀内部应力应变变化的波动性;推导出铣刀应力波衰减后的应力值和衰减率的求解方程;依据高速铣刀应力波的传播特性,确定应力波汇集区域,利用铣刀该区域跨尺度关联分析结果,结合分子动力学仿真,识别出铣刀局部宏介观损伤,如图17所示。图中,FG1为铣削工艺参数,FG2为铣刀结构参数,FG2为铣刀组件材料参数,S(x,ε)为铣刀原子群能量,f(x,ε)为铣刀原子群运动特征,Mt为铣刀原子群构型,Wt为元素数百分比,At为元素质量百分比,FRm为介观损伤特征,FRs为宏观损伤特征。
采用上述方法分析直径63mm四齿高速铣刀的高频应力波汇集区域,并对其进行跨尺度关联分析,结果如图18至图25所示。图中,铣刀刀体材料为42CrMo,螺钉材料为35CrMo,刀片材料为硬质合金。
结果表明,根据铣刀组件内部应力波传递特点,其能量主要集中在刀体的齿根处,宏观上损伤主要表现形式是以剪切断裂为主的组合型损伤,螺钉与刀片定位面组合变形和螺钉头上下部分均受到不同程度的拉伸作用,沿轴向切削力方向受到剪切载荷。介观上刀体齿根处的损伤类型为位错滑移、位错攀移、位错塞积、微裂纹扩展、晶界迁移,螺钉与刀片定位面的损伤类型分为位错滑移、位错攀移、位错塞积、微裂纹扩展、晶界迁移、晶面解理。
Claims (1)
1.一种断续冲击作用下的高速铣刀波动力学损伤识别方法,其特征在于由以下步骤构成:
一、建立高速铣刀瞬态切削力模型和质点动力学方程
根据铣刀切削工艺参数、铣刀结构参数、铣刀组件材料参数,获取铣刀瞬态切削力波形满足式(1):
式(1)中,p为作用于刀齿的单位切削层面积上的切削力,fz为每齿进给量,ap为轴向铣削深度,Fi(t)为第i个刀齿瞬态切削力,为随时间变化的刀齿瞬时位置角,m为参与切削的刀齿数量,F(t)为m个刀齿的瞬态切削力,t为时间,
利用铣刀质点运动微分方程表征其局部结构的动力响应,质点动力学方程满足式(2):
<mrow>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mi>T</mi>
<mi>&rho;</mi>
</mfrac>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>&rho;</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式(2)中,u为组件内质点相对位移,T为弦内张力,ρ为弦线密度,x为波阵面随时间变化的空间位置,Fi(x,t)为刀齿瞬态切削力在横向振动的激励;式(2)中,质点运动特征值为质点的位置u(x,t)和质点运动速度v的方程为:
<mrow>
<mi>v</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
二、建立高速铣刀应力波的波动方程
利用铣刀结构变化矩阵、铣刀组件材料本构关系、铣刀组件质点相对运动,构建铣刀组件内应力波的波动方程,
首先根据铣刀结构将其简化为轴对称结构,获得运动平衡方程为:
<mrow>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>v</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&sigma;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>&sigma;</mi>
<mfrac>
<mrow>
<mi>d</mi>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
将材料本构关系定义为:
σ=Eε (5)
采用式(4)、式(5),并根据牛顿第二定律,建立高速铣刀波动方程为
<mrow>
<msubsup>
<mi>c</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mi>d</mi>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
<msubsup>
<mi>c</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<msubsup>
<mi>&zeta;c</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
c0 2=E/ρ0,ζ=-dA(x)/A(x)dx (7)
波速是描述铣刀应力波的重要参数满足式(8):
<mrow>
<mi>C</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
上述式(4)至式(8)中:ρ0为材料密度,A(x)为刀具结构变化矩阵,v为刀具内的质点运动速度,σ为应力值;E为弹性模量,ε为应变;C为应力波的波速,x为波阵面随时间变化的空间位置;c0 2为材料内部应力波的传播速度,ζ为铣刀组件结构的修正系数;
三、求解高速铣刀应力波衰减应力值和衰减率
定义衰减率η为单位应力经过单位位移的衰减数值:
η=-(1/σx)dσx/dx (9)
对式(9)两边积分可得:
σx=σ0e-ηx (10)
应力波在铣刀内部经过多次反射后的波速C与初始波速c有如下关系:
τ=l/C=λ/c N=(C/c)(λ/l) (11)
高速面铣刀切削时承受的切削力冲击载荷所形成的冲击波长为冲击传递距离L的2倍,因此,式(11)可改写为下式:
N=2CL/(cl) (12)
对式(12)两边做偏微分变化,其中等式右端对L与l进行偏微分变化,得出下公式(13):
q=2ρ2C/(ρ1c+ρ2C) (13)
瞬态切削力载荷作用下,高速铣刀内部形成的应力波反射次数趋于无穷大的时候,铣刀刀齿区域内应力波衰减后的应力值σx为:
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>q</mi>
<mo>&lsqb;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&eta;x</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
</mrow>
</msup>
<mo>+</mo>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>&eta;x</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>q</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>n&eta;x</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
</mrow>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>q</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>n&eta;x</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
对式(14)进行整理得衰减率的求解方程为:
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>x</mi>
</mfrac>
<mi>ln</mi>
<mo>&lsqb;</mo>
<msup>
<mi>q</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>-</mo>
<mi>q</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>4</mn>
<msup>
<mi>B</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>q</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>-</mo>
<mi>q</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>B</mi>
<mo>=</mo>
<msub>
<mi>&sigma;</mi>
<mi>x</mi>
</msub>
<mo>/</mo>
<msub>
<mi>&sigma;</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
式(9)至式(15)中,η为衰减率,σx为衰减后应力值,σ0为x=0处的应力值(Mpa),x为应力波传播距离;τ为冲击作用时间,λ为波长,c为应力波初始速度,N为反射次数,L为应力波传播距离,l为应力波在铣刀中传播的长度,q为应力波的反射效率,ρ1和ρ2分别为应力波传播时所经过的不同介质材料密度,B为应力传递比,σi为应力波单次反射后的应力,ηxAB为应力波单次反射的衰减率;
四、识别高速铣刀波动力学损伤
依据高速铣刀应力波的传播特性,确定应力波汇集区域,利用铣刀该区域跨尺度关联分析结果,结合分子动力学仿真,识别出铣刀局部宏观和介观损伤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610966431.2A CN106334970B (zh) | 2016-10-28 | 2016-10-28 | 一种断续冲击载荷作用下的高速铣刀波动力学损伤识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610966431.2A CN106334970B (zh) | 2016-10-28 | 2016-10-28 | 一种断续冲击载荷作用下的高速铣刀波动力学损伤识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106334970A CN106334970A (zh) | 2017-01-18 |
CN106334970B true CN106334970B (zh) | 2018-06-01 |
Family
ID=57841017
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610966431.2A Expired - Fee Related CN106334970B (zh) | 2016-10-28 | 2016-10-28 | 一种断续冲击载荷作用下的高速铣刀波动力学损伤识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106334970B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108436596B (zh) * | 2017-06-06 | 2019-10-25 | 哈尔滨理工大学 | 一种基于高速铣刀组件原子群构型的铣刀损伤预后方法 |
CN109940459B (zh) * | 2019-04-10 | 2020-06-05 | 哈尔滨理工大学 | 一种高效铣刀损伤的多尺度识别方法 |
CN110988135B (zh) * | 2019-12-23 | 2022-03-29 | 北京工业大学 | 一种无约束振动处理下疲劳损伤修复振动参数的确定方法 |
CN114861395B (zh) * | 2022-03-30 | 2024-08-16 | 西北核技术研究所 | 一维杆中弹塑性应力波参量和能量解析方法及模型 |
CN114769692B (zh) * | 2022-04-19 | 2023-01-03 | 哈尔滨理工大学 | 高效铣刀后刀面摩擦磨损边界跨尺度识别方法 |
CN116738730B (zh) * | 2023-06-17 | 2023-12-05 | 哈尔滨理工大学 | 高效铣刀后刀面摩擦应力波传播与衰减特性的解算方法 |
CN116956681B (zh) * | 2023-07-26 | 2024-01-30 | 哈尔滨理工大学 | 铣削钛合金应力波叠加特性与刀齿应力波能量解算方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS55106751A (en) * | 1979-01-31 | 1980-08-15 | Yoshiaki Kakino | Tool damage detector |
US4694686A (en) * | 1984-06-18 | 1987-09-22 | Borg-Warner Corporation | Cutting tool wear monitor |
CN102564314B (zh) * | 2011-12-06 | 2014-04-16 | 上海交通大学 | 用于检测立铣刀磨损状态的正交视觉检测系统 |
CN202607422U (zh) * | 2012-05-30 | 2012-12-19 | 襄樊学院 | 采用dlp的铣刀磨损在机三维测量装置 |
CN104647138A (zh) * | 2013-11-19 | 2015-05-27 | 湖北文理学院 | 用于铣床铣刀磨损状态监测的微机监控装置 |
CN105619180B (zh) * | 2016-04-05 | 2017-10-17 | 哈尔滨理工大学 | 一种振动作用下的高速铣刀刀齿磨损差异性检测方法 |
CN105945651A (zh) * | 2016-06-22 | 2016-09-21 | 哈尔滨理工大学 | 球头铣刀精密铣削用的刀具磨损在线检测装置及检测方法 |
-
2016
- 2016-10-28 CN CN201610966431.2A patent/CN106334970B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN106334970A (zh) | 2017-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106334970B (zh) | 一种断续冲击载荷作用下的高速铣刀波动力学损伤识别方法 | |
Dalguer et al. | The importance of the dynamic source effects on strong ground motion during the 1999 Chi-Chi, Taiwan, earthquake: Brief interpretation of the damage distribution on buildings | |
Argon et al. | Cavity formation from inclusions in ductile fracture | |
Lim et al. | Wave motion in auxetic solids | |
Chizari et al. | Experimental and numerical study of water jet spot welding | |
He et al. | Experimental and numerical analysis of water jet peening on 6061 aluminum alloy | |
Correig et al. | Aftershock series of event February 18, 1996: An interpretation in terms of self‐organized criticality | |
Dali et al. | Comparison between dynamic and non-dynamic cutting tool option in FEM simulation for producing dimple structure | |
Xu et al. | Rock breaking model under dynamic load with the application of torsional and axial percussion hammer | |
Yuen et al. | Sandwich panels subjected to blast loading | |
CN108436596B (zh) | 一种基于高速铣刀组件原子群构型的铣刀损伤预后方法 | |
Rizzo et al. | Analysis of the geometric parameters of a solitary waves-based harvester to enhance its power output | |
Uenishi et al. | Dynamic fragmentation of ice spheres: two specific fracture patterns | |
Navarro et al. | Waves in periodic dissipative laminate metamaterial generated by plate impact | |
Choi et al. | Uncertainty assessment of structural modeling in the seismic response analysis of nuclear facilities | |
Liu et al. | Analytical and experimental study of determining the optimal number of wedge shape cutting teeth in coring bits used in percussive drilling | |
Lovely et al. | Regions of reduced static stress drop near fault tips for large strike-slip earthquakes | |
Ji et al. | Average static stress drop of shallow crustal earthquakes inferred from earthquake source spectra | |
Harkness et al. | Optimization of ultrasonic horns for momentum transfer and survivability in high-frequency/low frequency planetary drill tools | |
Radchenko | Influence of Reinforcement and Attack Angle on the Destruction of a Concrete Barrier Under a High-Speed Impact | |
Mehtiyev et al. | Review of basic equations for evaluating drilling efficiency | |
Paglialunga et al. | Origin of the temporal evolution of elastic properties during laboratory seismic cycle. | |
Stone et al. | Sensitivity of modeled topographic effects to kinematic source parameters in 3D simulations of M6. 5-7.0 earthquakes | |
Rajagopal et al. | Numerical studies of ultrasonic waves on the surface of curved objects | |
Pankow et al. | Bayesian SPAC Tomography with a Dense Nodal Array |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180601 |
|
CF01 | Termination of patent right due to non-payment of annual fee |