CN112372372B - Efficient milling cutter accumulated friction wear boundary identification and verification method - Google Patents
Efficient milling cutter accumulated friction wear boundary identification and verification method Download PDFInfo
- Publication number
- CN112372372B CN112372372B CN202011439558.1A CN202011439558A CN112372372B CN 112372372 B CN112372372 B CN 112372372B CN 202011439558 A CN202011439558 A CN 202011439558A CN 112372372 B CN112372372 B CN 112372372B
- Authority
- CN
- China
- Prior art keywords
- cutter
- milling
- boundary
- wear
- tooth
- 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
- 238000003801 milling Methods 0.000 title claims abstract description 338
- 238000000034 method Methods 0.000 title claims abstract description 69
- 238000012795 verification Methods 0.000 title claims abstract description 11
- 238000004088 simulation Methods 0.000 claims abstract description 105
- 230000001186 cumulative effect Effects 0.000 claims abstract description 51
- 238000002474 experimental method Methods 0.000 claims abstract description 24
- 230000009471 action Effects 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 14
- 238000010276 construction Methods 0.000 claims abstract description 6
- 238000005520 cutting process Methods 0.000 claims description 57
- 238000009826 distribution Methods 0.000 claims description 45
- 238000005259 measurement Methods 0.000 claims description 22
- 230000008859 change Effects 0.000 claims description 18
- 238000013461 design Methods 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 8
- 238000010168 coupling process Methods 0.000 claims description 8
- 238000005859 coupling reaction Methods 0.000 claims description 8
- 238000009795 derivation Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 238000000691 measurement method Methods 0.000 abstract description 8
- 238000012512 characterization method Methods 0.000 abstract description 4
- 230000007547 defect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 26
- 229910001069 Ti alloy Inorganic materials 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 230000007704 transition Effects 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000003754 machining Methods 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 238000005482 strain hardening Methods 0.000 description 2
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 101150011213 uvs-1 gene Proteins 0.000 description 1
Images
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/0904—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 before or after machining
- B23Q17/0914—Arrangements for measuring or adjusting cutting-tool geometry machine tools
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Numerical Control (AREA)
Abstract
本发明为一种高效铣刀累积摩擦磨损边界识别与验证方法,属于铣刀技术领域,是针对现有铣削实验方法研究刀齿后刀面最大磨损宽度随铣削行程变化,无法揭示铣削振动和刀齿误差作用下刀齿后刀面摩擦磨损边界形成过程的缺陷所提出,其包括高效铣刀刀齿后刀面摩擦磨损边界测量与表征方法、高效铣刀轴向分层铣削仿真模型及边界条件构建方法、刀齿后刀面瞬时摩擦磨损边界识别方法、刀齿后刀面累积摩擦磨损边界解算及验证方法。本发明提出高效铣刀刀齿后刀面摩擦磨损边界测量与表征方法,建立刀齿后刀面瞬时摩擦磨损边界仿真识别方法,提出刀齿后刀面累积摩擦磨损边界解算方法,为揭示刀齿后刀面累积摩擦磨损边界形成过程提供了一种有效的模型。
The invention relates to a high-efficiency milling cutter cumulative friction and wear boundary identification and verification method, which belongs to the technical field of milling cutters, and is aimed at the existing milling experimental method to study the variation of the maximum wear width of the cutter tooth flank with the milling stroke, and cannot reveal the milling vibration and cutter The defects in the formation process of the flank friction and wear boundary of the cutter tooth under the action of tooth error are proposed, which include the measurement and characterization method of the friction and wear boundary of the cutter tooth flank of the high-efficiency milling cutter, the simulation model and boundary conditions of the high-efficiency milling cutter axial layered milling The construction method, the identification method of the instantaneous friction and wear boundary of the tooth flank, the calculation and verification method of the cumulative friction and wear boundary of the tooth flank. The present invention proposes a method for measuring and characterizing the boundary of friction and wear on the flank of a high-efficiency milling cutter, establishes a simulation and identification method for the boundary of instantaneous friction and wear on the flank of a cutter, and proposes a method for calculating the boundary of cumulative friction and wear on the flank of a cutter. The formation process of cumulative friction and wear boundary on flank flank provides an effective model.
Description
技术领域:Technical field:
本发明属于铣刀技术领域,具体涉及一种高效铣刀累积摩擦磨损边界识别与验证方法。The invention belongs to the technical field of milling cutters, and in particular relates to a method for identifying and verifying the accumulated friction and wear boundary of an efficient milling cutter.
背景技术:Background technique:
高效铣刀因其高效的切削性能而广泛应用于航空航天领域大型钛合金结构件加工中,在铣削过程中,受铣削振动和刀齿误差影响,铣刀刀齿后刀面累积摩擦磨损边界的形成与演变存在不确定性,使得高效铣刀的整体使用寿命难以准确的评估和预测,由此导致工件加工质量不稳定,生产成本增加。High-efficiency milling cutters are widely used in the machining of large-scale titanium alloy structural parts in the aerospace field due to their efficient cutting performance. During the milling process, affected by milling vibration and cutter tooth errors, the cumulative friction and wear boundary of the flank of milling cutter teeth is limited. There is uncertainty in the formation and evolution, which makes it difficult to accurately evaluate and predict the overall service life of high-efficiency milling cutters, resulting in unstable workpiece processing quality and increased production costs.
受铣削振动和刀齿误差的影响,铣刀刀齿后刀面的摩擦磨损边界形成过程的不确定性,是预测刀具使用寿命的关键。目前,虽然已有关于铣刀刀齿后刀面摩擦磨损的研究,但是主要集中在通过铣削实验揭示刀齿后刀面摩擦磨损宽度随铣削行程的变化特性,通过对刀齿后刀面磨损宽度的研究却无法揭示随铣削形成变化铣刀刀齿后刀面与工件之间的接触关系,同时由于实验局限性,无法通过铣削实验获取刀齿后刀面瞬时摩擦磨损边界,进而导致无法通过铣削实验的方法揭示出刀齿后刀面瞬时与累积摩擦磨损边界曲线之间的关系和累积摩擦磨损边界曲线的形成过程。Affected by milling vibration and tooth error, the uncertainty of the formation process of the friction and wear boundary of the flank of the milling cutter is the key to predicting the service life of the tool. At present, although there have been researches on the friction and wear of the flank of milling cutter teeth, they mainly focus on revealing the variation characteristics of the friction and wear width of the tooth flank with the milling stroke through milling experiments. However, due to the limitations of the experiment, the instantaneous friction and wear boundary of the flank of the cutter teeth cannot be obtained through the milling experiment, which leads to the inability to pass the milling process. The experimental method revealed the relationship between the instantaneous flank face and the cumulative friction and wear boundary curve and the formation process of the cumulative friction and wear boundary curve.
发明内容:Invention content:
本发明为克服现有铣刀刀齿后刀面累积摩擦磨损边界的形成与演变过程,难以准确地评估和预测高效铣刀整体的使用寿命的缺陷,提供了一种高效铣刀累积摩擦磨损边界识别与验证方法,该方法利用铣刀刀齿后刀面摩擦磨损边界测量与表征方法、刀齿后刀面瞬时摩擦磨损仿真边界识别方法、铣刀刀齿后刀面累积摩擦磨损仿真边界解算方法,对铣刀刀齿后刀面累积摩擦磨损边界解算方法的正确性进行验证。The invention provides a high-efficiency milling cutter cumulative friction and wear boundary in order to overcome the defect that it is difficult to accurately evaluate and predict the overall service life of the high-efficiency milling cutter in the formation and evolution process of the existing milling cutter tooth flank cumulative friction and wear boundary. Identification and verification method. This method uses the measurement and characterization method of the friction and wear boundary of the flank of the milling cutter tooth, the identification method of the instantaneous friction and wear simulation boundary of the cutter tooth flank, and the calculation of the cumulative friction and wear simulation boundary of the milling cutter tooth flank. method to verify the correctness of the calculation method of the cumulative friction and wear boundary of the flank of the milling cutter.
本发明采用的技术方案在于:一种高效铣刀累积摩擦磨损边界识别与验证方法,具体包括以下步骤:The technical scheme adopted in the present invention is: a high-efficiency milling cutter cumulative friction and wear boundary identification and verification method, which specifically includes the following steps:
步骤A:提出高效铣刀刀齿后刀面摩擦磨损边界测量与表征方法,构建铣刀刀齿后刀面摩擦磨损边界测量坐标系,在测量坐标系中表征铣刀刀齿后刀面摩擦磨损边界的曲线形状和位置;Step A: Propose a method for measuring and characterizing the friction and wear boundary of the flank flank of a high-efficiency milling cutter, construct a boundary measurement coordinate system for the friction and wear of the flank of the milling cutter, and characterize the friction and wear of the flank of the milling cutter in the measurement coordinate system the shape and position of the curve of the boundary;
步骤B:将铣削方式、铣削振动、刀齿误差和不同层待加工表面差异性考虑为仿真边界条件,提出高效铣刀轴向分层铣削仿真模型及边界条件构建方法;Step B: Consider the milling method, milling vibration, cutter tooth error and the difference of the surface to be machined in different layers as simulation boundary conditions, and propose a high-efficiency milling cutter axial layer milling simulation model and boundary condition construction method;
步骤C:利用高效铣刀轴向分层铣削仿真模型进行高效铣刀轴向分层铣削热力耦合场仿真,并提取铣刀刀齿后刀面瞬时温度、等效应力、等效应变分布场,根据铣刀刀齿后刀面瞬时摩擦磨损边界形成过程,提出刀齿后刀面瞬时摩擦磨损仿真边界识别方法;Step C: Use the high-efficiency milling cutter axial layered milling simulation model to simulate the thermal coupling field of high-efficiency milling cutter axial layered milling, and extract the instantaneous temperature, equivalent stress, and equivalent strain distribution field of the milling cutter tooth flank, According to the formation process of the instantaneous friction and wear boundary of the flank of the milling cutter tooth, a method for identifying the boundary of the instantaneous friction and wear simulation of the flank of the cutter tooth is proposed.
步骤D:根据步骤C获取不同时刻刀齿后刀面瞬时摩擦磨损边界,分析刀齿后刀面累积摩擦磨损边界形成过程,提出铣刀刀齿后刀面累积摩擦磨损仿真边界解算方法,并计算通过仿真和实验获得的铣刀刀齿后刀面累积摩擦磨损边界曲线灰色关联度和平均相对误差,验证铣刀刀齿后刀面累积摩擦磨损边界解算方法的正确性。Step D: According to step C, the instantaneous friction and wear boundary of the flank face of the cutter tooth at different times is obtained, the formation process of the cumulative friction and wear boundary of the flank face of the cutter tooth is analyzed, and the calculation method of the cumulative friction and wear simulation boundary of the flank face of the milling cutter is proposed. The grey correlation degree and average relative error of the cumulative friction and wear boundary curve of the milling cutter tooth flank obtained through simulation and experiment are calculated to verify the correctness of the calculation method of the milling cutter tooth flank cumulative friction and wear boundary.
优选地,所述步骤A,还包括以下步骤:Preferably, the step A further comprises the following steps:
步骤A1:振动作用下高效铣刀的铣削方式:在工件坐标系中,利用铣刀坐标系的偏置和振动作用下铣刀切削状态反映铣削振动对高效铣刀切削过程的影响;Step A1: The milling method of the high-efficiency milling cutter under the action of vibration: in the workpiece coordinate system, the cutting state of the milling cutter under the action of the offset and vibration of the milling cutter coordinate system is used to reflect the impact of milling vibration on the cutting process of the high-efficiency milling cutter;
步骤A2:根据高效铣刀结构特点,确定铣刀刀齿误差及其后刀面摩擦磨损测量方法;Step A2: According to the structural characteristics of the high-efficiency milling cutter, determine the cutter tooth error of the milling cutter and the measurement method for the friction and wear of the flank;
步骤A3:利用铣削振动与刀齿误差影响下的工件坐标系、铣刀坐标系、刀齿坐标系和后刀面摩擦磨损边界测量坐标系之间的关系矩阵,反映铣削振动与刀齿误差对刀齿后刀面与工件瞬时接触关系的影响;Step A3: Use the relationship matrix between the workpiece coordinate system, the milling cutter coordinate system, the cutter coordinate system and the flank friction and wear boundary measurement coordinate system under the influence of the milling vibration and the cutter tooth error to reflect the relationship between the milling vibration and the cutter tooth error. The influence of the instantaneous contact relationship between the tooth flank and the workpiece;
步骤A4:通过在刀齿后刀面磨损边界测量坐标系中依次度量出刀齿后刀面摩擦磨损边界端点和拐点,构建铣刀刀齿后刀面摩擦磨损上下边界曲线,确定刀齿后刀面摩擦磨损边界。Step A4: Construct the upper and lower boundary curves of the friction and wear of the flank of the milling cutter by measuring the end points and inflection points of the friction and wear boundary of the flank of the cutter in turn in the measurement coordinate system of the tooth flank, and determine the flank of the cutter Surface friction and wear boundary.
优选地,所述步骤B,还包括以下步骤:Preferably, the step B further comprises the following steps:
步骤B1:将铣削方式、刀齿误差、铣削振动和不同层待加工表面差异性考虑为仿真边界条件,设计高效铣刀轴向分层铣削仿真模型整体结构,确定并解算仿真模型设计变量;Step B1: Consider the milling method, cutter tooth error, milling vibration and the difference of the surface to be machined in different layers as simulation boundary conditions, design the overall structure of the high-efficiency milling cutter axial layered milling simulation model, and determine and solve the simulation model design variables;
步骤B2:仿真模型铣削单元构建,求解振动作用下铣刀运动轨迹和瞬时切削姿态角,确定铣削单元横截面模型,考虑铣削振动和刀齿误差构建铣削单元;Step B2: Build a milling unit of simulation model, solve the motion trajectory and instantaneous cutting attitude angle of the milling cutter under the action of vibration, determine the cross-sectional model of the milling unit, and construct the milling unit considering the milling vibration and cutter tooth error;
步骤B3:在仿真过程中,根据轴向分层铣削实验中振动作用下的铣刀运动轨迹,确定铣刀运动轨迹点的离散时间间隔,对仿真过程中振动作用下的铣刀运动轨迹点进行解算。Step B3: In the simulation process, according to the movement trajectory of the milling cutter under the action of vibration in the axial layered milling experiment, the discrete time interval of the movement trajectory points of the milling cutter is determined, and the movement trajectory points of the milling cutter under the action of vibration in the simulation process are carried out. Solve.
优选地,所述步骤C,还包括以下步骤:Preferably, described step C, also comprises the following steps:
步骤C1:通过高效铣刀轴向分层铣削仿真模型进行热力耦合场仿真,提取铣刀刀齿后刀面任意瞬时切削等效应力场、等效应变场、温度场分布情况;Step C1: perform thermal coupling field simulation through the high-efficiency milling cutter axial layered milling simulation model, and extract the equivalent stress field, equivalent strain field and temperature field distribution of any instantaneous cutting of the milling cutter tooth flank;
步骤C2:提取铣刀刀齿后刀面不同位置处节点等效应力、等效应变和温度数值,构建刀齿后刀面等效应力、等效应变、温度分布曲面,平行于刀齿后刀面投影面上切削刃中点法矢量方向过切削刃上节点依次作垂直于分布曲面的截平面,并构建所有刀齿后刀面不同位置处截平面上节点等效应力、等效应变、温度分布曲线;Step C2: Extract the equivalent stress, equivalent strain and temperature values of nodes at different positions on the flank of the milling cutter tooth, and construct the equivalent stress, equivalent strain, and temperature distribution surfaces of the flank of the milling cutter, which are parallel to the flank of the tooth. The normal vector direction of the midpoint of the cutting edge on the projection surface of the cutting edge is taken as the section plane perpendicular to the distribution surface in turn through the nodes on the cutting edge, and the equivalent stress, equivalent strain, temperature of the nodes on the section plane at different positions on the flank of all teeth are constructed. distribution curve;
步骤C3:采用一元高次多项式依次拟合刀齿后刀面各截平面处等效应力、等效应变和温度分布曲线,通过求导得节点变化率分布曲线;Step C3: Use a univariate high-order polynomial to fit the equivalent stress, equivalent strain and temperature distribution curves at each section plane of the flank face of the cutter in turn, and obtain the node change rate distribution curve by derivation;
步骤C4:以刀齿材料屈服强度为判据,识别出所有截平面上等效应力大于且最接近于屈服的临界特征节点在刀齿后刀面摩擦磨损测量坐标系中,不同位置处截平面上特征点共同构成刀齿后刀面瞬时摩擦磨损仿真上边界;Step C4: Using the yield strength of the tooth material as the criterion, identify the critical feature nodes whose equivalent stress is greater than and closest to yield on all section planes. In the tool tooth flank friction and wear measurement coordinate system, the section planes at different positions The upper feature points together constitute the upper boundary of the instantaneous friction and wear simulation of the tooth flank;
步骤C5:依次识别出刀齿后刀面各截平面处温度、等效应力、等效应变变化率发生明显突变的节点,以突变节点为特征点,在刀齿后刀面摩擦磨损边界测量坐标系中,构建刀齿后刀面上节点温度、等效应力、等效应变突变特征曲线,节点温度、等效应变、等效应力突变特征曲线共同形成的最外侧轮廓构成刀齿后刀面瞬时摩擦磨损仿真下边界。Step C5: Identify the nodes where the temperature, equivalent stress, and equivalent strain rate of change at each section plane of the tooth flank have abrupt changes in turn, take the sudden change node as the feature point, and measure the coordinates on the friction and wear boundary of the tooth flank In the system, the node temperature, equivalent stress, and equivalent strain mutation characteristic curves on the flank of the tooth flank are constructed. Lower boundary of friction and wear simulation.
优选地,所述步骤D,还包括以下步骤:Preferably, the step D further comprises the following steps:
步骤D1:通过不同铣削时刻刀齿后刀面瞬时摩擦磨损边界曲线,以累积摩擦磨损下边界曲线为例,提出铣刀刀齿后刀面累积摩擦磨损边界解算方法;Step D1: According to the instantaneous friction and wear boundary curve of the flank of the cutter teeth at different milling times, and taking the lower boundary curve of the cumulative friction and wear as an example, a method for calculating the boundary of the cumulative friction and wear of the flank of the milling cutter tooth is proposed;
步骤D2:设计刀齿后刀面累积摩擦磨损上下边界仿真与实验对比方案,采用铣刀刀齿后刀面累积摩擦磨损边界曲线解算方法,构建铣刀刀齿后刀面累积摩擦磨损仿真上下边界曲线,将实验摩擦磨损上下边界曲线与仿真上下边界曲线进行对比;Step D2: Design the simulation and experimental comparison scheme of the upper and lower boundary of the cumulative friction and wear of the flank of the cutter tooth, and use the calculation method of the boundary curve of the cumulative friction and wear of the cutter tooth flank to construct the upper and lower simulation of the cumulative friction and wear of the milling cutter tooth flank. Boundary curve, compare the upper and lower boundary curves of experimental friction and wear with the upper and lower boundary curves of simulation;
步骤D3:计算仿真与实验累积摩擦磨损边界曲线之间的灰色关联度值和平均相对误差,分析仿真与实验累积摩擦磨损边界曲线之间形状相似度和位置重合度,验证累积摩擦磨损边界解算方法的正确性。Step D3: Calculate the grey correlation value and the average relative error between the simulation and the experimental cumulative friction and wear boundary curve, analyze the shape similarity and position coincidence between the simulation and the experimental cumulative friction and wear boundary curve, and verify the cumulative friction and wear boundary solution correctness of the method.
本发明的有益效果是:The beneficial effects of the present invention are:
1、本发明充分考虑了铣刀刀齿后刀面整个磨损边界的特征信息,通过建立刀齿后刀面摩擦磨损边界测量坐标系,以摩擦磨损边界上左右两侧端点和拐点为特征点,构建铣刀刀齿后刀面摩擦磨损上下边界特征点曲线,能够更加完整的表征出刀齿后刀面摩擦磨损边界曲线的形状和位置信息,通过铣削振动和刀齿误差影响下的工件坐标系、铣刀坐标系、刀齿坐标系、和后刀面摩擦磨损边界测量坐标系之间的关系矩阵,有利于揭示铣削过程中的刀工接触关系和摩擦磨损边界的形成过程。1. The present invention fully considers the characteristic information of the entire wear boundary of the flank face of the milling cutter, by establishing a coordinate system for measuring the friction and wear boundary of the flank of the cutter tooth, and taking the endpoints and inflection points on the left and right sides of the friction and wear boundary as feature points, Constructing the upper and lower boundary characteristic point curves of the friction and wear of the flank face of the milling cutter can more completely characterize the shape and position information of the boundary curve of the friction and wear of the flank of the cutter tooth. , the relationship matrix between the cutter coordinate system, the cutter tooth coordinate system, and the flank friction and wear boundary measurement coordinate system, which is beneficial to reveal the cutter-worker contact relationship and the formation process of the friction and wear boundary during the milling process.
2、本发明充分考虑了铣削方式、铣削振动、刀齿误差、不同层待加工表面差异性对仿真模型结构的影响,将铣削方式、铣削振动、刀齿误差和不同层待待加工表面考虑为仿真边界条件,提出一种新的高效铣刀轴向分层铣削仿真模型设计方法,采用该方法设计的仿真模型进行高效铣刀轴向分层铣削热力耦合场仿真,其仿真结果更加准确。2. The present invention fully considers the influence of milling method, milling vibration, cutter tooth error, and the difference of surfaces to be machined in different layers on the structure of the simulation model, and considers the milling method, milling vibration, cutter tooth error and surfaces to be machined in different layers as By simulating the boundary conditions, a new simulation model design method for high-efficiency milling cutter axial layered milling is proposed. The simulation model designed by this method is used to simulate the thermal coupling field of high-efficiency milling cutter axial layered milling, and the simulation results are more accurate.
3、本发明基于铣刀刀齿后刀面瞬时摩擦磨损仿真边界形成过程,根据刀齿后刀面瞬时切削节点温度场、等效应力、等效应变分布特性,提出了铣刀刀齿后刀面瞬时摩擦磨损仿真上下边界识别方法,铣刀刀齿后刀面累积摩擦磨损仿真边界解算方法,揭示了铣刀刀齿后刀面累积摩擦磨损边界形成过程,提出了铣刀刀齿后刀面累积摩擦磨损边界解算方法验证方法。3. Based on the simulation boundary formation process of the instantaneous friction and wear of the flank of the milling cutter, the present invention proposes a flank of the milling cutter according to the instantaneous cutting node temperature field, equivalent stress and equivalent strain distribution characteristics of the flank of the milling cutter. The upper and lower boundary identification method of instantaneous friction and wear simulation of face, and the calculation method of cumulative friction and wear simulation boundary of milling cutter tooth flank, revealing the formation process of milling cutter tooth flank cumulative friction and wear boundary, and putting forward the milling cutter tooth flank Verification method of surface cumulative friction and wear boundary solution method.
附图说明:Description of drawings:
图1为高效铣刀刀齿后刀面累积摩擦磨损边界解算方法的流程图;Fig. 1 is a flow chart of a method for calculating the boundary of cumulative friction and wear on the flank of a high-efficiency milling cutter;
图2为振动作用下高效铣刀铣削的状态图;Figure 2 is a state diagram of high-efficiency milling cutter milling under the action of vibration;
图3为振动作用下铣刀瞬时切削的状态图;Fig. 3 is the state diagram of the instantaneous cutting of the milling cutter under the action of vibration;
图4为高效铣刀结构及其刀齿误差测量方法,其中图4(a)为高效铣刀整体结构示意图、图4(b)为铣刀刀齿结构示意图;Fig. 4 is a high-efficiency milling cutter structure and a cutter tooth error measurement method, wherein Fig. 4(a) is a schematic diagram of the overall structure of a high-efficiency milling cutter, and Fig. 4(b) is a schematic diagram of the milling cutter cutter tooth structure;
图5为刀齿后刀面摩擦磨损边界测量与磨损深度表征方法示意图,其中图5(a)为刀齿后刀面摩擦磨损边界测量方法、图5(b)为刀齿后刀面磨损深度;Figure 5 is a schematic diagram of the method for measuring the friction and wear boundary of the tooth flank and characterizing the wear depth, in which Figure 5(a) is the measurement method of the friction and wear boundary of the tooth flank, and Figure 5(b) is the wear depth of the tooth flank. ;
图6为铣刀刀齿后刀面摩擦磨损上下边界的曲线图;Figure 6 is a graph of the upper and lower boundaries of friction and wear of the flank of the milling cutter tooth;
图7为分层铣削第4层有效铣削行程振动加速度时域信号示意图;Fig. 7 is a schematic diagram of the vibration acceleration time domain signal of the fourth layer of effective milling in layered milling;
图8为铣削行程5m的高效铣刀刀齿底刃后刀面磨损状态图,其中图8(a)为刀齿一后刀面磨损状态图、图8(b)为刀齿二后刀面磨损状态图、图8(c)为刀齿三后刀面磨损状态图;Figure 8 is a diagram of the wear state of the flank of the bottom edge of a high-efficiency milling cutter with a milling stroke of 5m, in which Figure 8(a) is a diagram of the wear state of the first tooth flank, and Figure 8(b) is the second tooth flank. The wear state diagram, Figure 8(c) is the wear state diagram of the third flank of the cutter teeth;
图9为高效铣刀刀齿后刀面摩擦磨损上下边界特征点分布曲线图,其中图9(a)为铣刀刀齿后刀面摩擦磨损上边界特征点分布曲线图、图9(b)为铣刀刀齿后刀面摩擦磨损下边界特征点分布曲线图;Figure 9 is the distribution curve of the upper and lower boundary feature points of the friction and wear of the flank of the high-efficiency milling cutter tooth, of which Figure 9(a) is the distribution curve of the upper boundary characteristic points of the friction and wear of the flank of the milling cutter tooth, and Figure 9(b) is the distribution curve of the lower boundary feature points of the friction and wear of the flank of the milling cutter tooth;
图10为高效铣刀轴向分层铣削仿真工件整体结构及单元横截面模型图,其中图10(a)为高效铣刀轴向分层铣削仿真工件模型整体结构、图10(b)为第i(i≥2)层铣削单元横截面模型图;Figure 10 shows the overall structure of the high-efficiency milling cutter axial layered milling simulation workpiece overall structure and unit cross-sectional model diagram, in which Figure 10(a) is the overall structure of the high-efficiency milling cutter axial layered milling simulation workpiece model, and Figure 10(b) is the first part. Model diagram of the cross-sectional model of the milling unit at layer i (i≥2);
图11为铣刀第i(i≥2)层切入切出的识别图;Fig. 11 is the identification diagram of the cutting in and out of the i-th (i≥2) layer of the milling cutter;
图12为振动作用下铣刀运动的轨迹示意图;Figure 12 is a schematic diagram of the trajectory of the milling cutter movement under the action of vibration;
图13为第i层铣刀瞬时切削姿态角的信号图,其中图13(a)为铣刀沿进给方向的偏转角、图13(b)为铣刀沿切削宽度方向的偏转角;Fig. 13 is the signal diagram of the instantaneous cutting attitude angle of the i-th layer milling cutter, wherein Fig. 13(a) is the deflection angle of the milling cutter along the feeding direction, and Fig. 13(b) is the deflection angle of the milling cutter along the cutting width direction;
图14为第i层铣削单元形成的实体图,其中图14(a)为第i层铣削单元的表面示意图,图14(b)为图14(a)中A-A处的剖面图,图14(c)为图14(a)中B-B处的剖面图;Fig. 14 is a solid view of the formation of the i-th layer milling unit, wherein Fig. 14(a) is a schematic surface view of the i-th layer milling unit, Fig. 14(b) is a cross-sectional view at A-A in Fig. 14(a), Fig. 14 ( c) is a sectional view at B-B in Figure 14(a);
图15为高效铣刀轴向分层铣削过程仿真图,其中图15(a)为铣刀切入阶段仿真图、图15(b)为铣刀切削中段仿真图、图15(c)为铣刀切出阶段仿真图;Figure 15 is the simulation diagram of the high-efficiency milling cutter axial layered milling process, in which Figure 15(a) is the simulation diagram of the milling cutter cutting stage, Figure 15(b) is the simulation diagram of the middle cutting section of the milling cutter, and Figure 15(c) is the milling cutter Cut-out stage simulation diagram;
图16为铣刀刀齿后刀面瞬时等效应力、等效应变场和温度场分布图,其中图16(a)为等效应力分布图、图16(b)为等效应变场分布图、图16(c)为温度场分布分布图;Figure 16 is the distribution diagram of instantaneous equivalent stress, equivalent strain field and temperature field on the flank of the milling cutter tooth, in which Figure 16(a) is the equivalent stress distribution diagram, and Figure 16(b) is the equivalent strain field distribution diagram , Figure 16(c) is the temperature field distribution map;
图17为平行于刀齿切削刃中点法矢量方向等间隔垂直于刀齿后刀面等效应力、等效应变、温度分布曲面的截平面,其中图17(a1)为刀齿后刀面等效应力分布图、图17(a2)为平行于刀齿后刀面投影面上切削刃中点法矢量方向过切削刃上节点依次作垂直于图17(a1)曲面的截平面、图17(b1)刀齿后刀面等效应变分布图、图17(b2)为平行于刀齿后刀面投影面上切削刃中点法矢量方向等间隔过切削刃上节点依次作垂直于图17(b1)曲面的截平面、图17(c1)为刀齿后刀面温度分布图、图17(c2)为平行于刀齿后刀面投影上切削刃中点法矢量方向过切削刃上节点依次作垂直于图17(c1)曲面的截平面;Figure 17 is a section plane parallel to the normal vector direction of the midpoint of the cutting edge of the tooth and perpendicular to the equivalent stress, equivalent strain and temperature distribution surface of the tooth flank at equal intervals, of which Figure 17(a1) is the tooth flank surface Equivalent stress distribution diagram, Fig. 17(a2) is parallel to the cutting edge normal vector direction on the projection surface of the flank face of the cutter, and the upper node of the cutting edge is in turn taken as a cross-sectional plane perpendicular to the surface of Fig. 17(a1), Fig. 17 (b1) The equivalent strain distribution diagram of the tooth flank, Figure 17 (b2) is parallel to the projection plane of the tooth flank and the normal vector direction of the midpoint of the cutting edge is equally spaced through the upper node of the cutting edge. (b1) The section plane of the curved surface, Fig. 17(c1) is the temperature distribution of the tooth flank, and Fig. 17(c2) is the normal vector direction of the midpoint of the cutting edge parallel to the projection of the tooth flank over the upper node of the cutting edge Make a section plane perpendicular to the curved surface of Figure 17(c1) in turn;
图18为刀齿后刀面不同位置处节点等效应力分布曲线图,其中图18(a)为截面A处节点应力分布曲线图、图18(b)为截面B处节点应力分布曲线图、图18(c)为截面C处节点应力分布曲线图;Figure 18 is a graph of the equivalent stress distribution of nodes at different positions on the flank of the cutter teeth, of which Figure 18(a) is a graph of the stress distribution of nodes at section A, and Figure 18(b) is a graph of the stress distribution of nodes at section B. Figure 18(c) is a graph showing the stress distribution curve of the node at section C;
图19为刀齿后刀面截面B处节点等效应变、温度分布曲线图,其中图19(a)为截面B处等效应变分布曲线图、图19(b)为截面B处温度分布曲线图;Figure 19 is a graph of the equivalent strain and temperature distribution at the node of the tooth flank section B, in which Figure 19(a) is the equivalent strain distribution curve at the section B, and Figure 19(b) is the temperature distribution curve at the section B. picture;
图20为刀齿后刀面截面B处等效应力、等效应变、温度变化率分布曲线图,其中图20(a)为等效应力变化率分布曲线图、图20(b)为等效应变变化率分布曲线图、图20(c)为温度变化率分布曲线图;Figure 20 is the distribution curve diagram of equivalent stress, equivalent strain and temperature change rate at the section B of the flank face of the cutter tooth, in which Figure 20(a) is the distribution curve diagram of the equivalent stress change rate, and Figure 20(b) is the equivalent effect Change rate distribution curve, Figure 20(c) is the temperature change rate distribution curve;
图21为铣刀刀齿后刀面瞬时摩擦磨损仿真上边界曲线图;Figure 21 is the upper boundary curve graph of the instantaneous friction and wear simulation of the flank of the milling cutter tooth;
图22为铣刀刀齿后刀面瞬时摩擦磨损仿真下边界识别曲线图;Fig. 22 is a graph showing the lower boundary identification curve of the instantaneous friction and wear simulation of the flank of the milling cutter tooth;
图23为铣刀刀齿后刀面累积摩擦磨损边界形成演变过程曲线图;Figure 23 is a graph showing the evolution process of the cumulative friction and wear boundary formation on the flank of the milling cutter tooth;
图24为铣刀刀齿一后刀面瞬时摩擦磨损上下边界曲线图,其中图24(a)为刀齿后刀面瞬时摩擦磨损上边界曲线图、图24(b)为刀齿后刀面瞬时摩擦磨损下边界曲线图;Figure 24 is the upper and lower boundary curves of the instantaneous friction and wear between the cutter tooth and the flank of the milling cutter, in which Figure 24(a) is the upper boundary curve of the instantaneous friction and wear of the cutter tooth flank, and Figure 24(b) is the cutter tooth flank surface. The lower boundary curve of instantaneous friction and wear;
图25为铣刀刀齿一仿真与实验累积摩擦磨损上下边界曲线对比曲线图。Fig. 25 is a graph comparing the upper and lower boundary curves of cumulative friction and wear of milling cutter teeth-simulation and experiment.
具体实施方式:Detailed ways:
如图1所示,本发明一种高效铣刀累积摩擦磨损边界识别与验证方法,是利用铣削振动和刀齿误差影响下的工件坐标系、铣刀坐标系、刀齿坐标系和后刀面摩擦磨损边界测量坐标系之间的关系矩阵,反映铣削振动和刀齿误差对刀齿后刀面与工件瞬时接触关系的影响,建立刀齿后刀面摩擦磨损边界测量坐标系,以摩擦磨损边界上左右两侧端点和拐点为特征点,提出铣刀刀齿后刀面摩擦磨损边界测量与表征方法,并在测量坐标系中完整表征了铣刀刀齿后刀面摩擦磨损边界曲线形状和位置;将铣削振动、刀齿误差、铣削方式、不同层待加工表面差异性考虑为仿真边界条件,提出高效铣刀轴向分层铣削仿真模型及边界条件构建方法;基于刀齿后刀面瞬时摩擦磨损边界形成过程和刀齿后刀面瞬时节点温度场、等效应力、等效应变分布特性,提出铣刀刀齿后刀面瞬时摩擦磨损上下边界仿真判据;基于瞬时摩擦磨损仿真判据,提出铣刀刀齿后刀面累积摩擦磨损边界解算方法,通过计算仿真与实验累积摩擦磨损边界曲线之间的灰色关联度值和平均相对误差,分析仿真与实验累积摩擦磨损边界曲线之间形状和位置的相似程度,验证铣刀刀齿后刀面累积摩擦磨损边界解算方法的正确性。As shown in Figure 1, a high-efficiency milling cutter cumulative friction and wear boundary identification and verification method of the present invention is to utilize the workpiece coordinate system, milling cutter coordinate system, cutter tooth coordinate system and flank under the influence of milling vibration and cutter tooth error. The relationship matrix between the friction and wear boundary measurement coordinate systems reflects the influence of milling vibration and cutter tooth error on the instantaneous contact relationship between the cutter tooth flank and the workpiece, and establishes the cutter tooth flank friction and wear boundary measurement coordinate system. The endpoints and inflection points on the upper left and right sides are the characteristic points, and a method for measuring and characterizing the friction and wear boundary of the flank of the milling cutter tooth is proposed, and the shape and position of the friction and wear boundary curve of the flank face of the milling cutter tooth are completely characterized in the measurement coordinate system. ; Considering the milling vibration, cutter tooth error, milling method, and the difference of the surface to be machined in different layers as simulation boundary conditions, a simulation model and boundary condition construction method for high-efficiency milling cutter axial layered milling are proposed; based on the instantaneous friction of the cutter tooth flank The formation process of the wear boundary and the instantaneous node temperature field, equivalent stress and equivalent strain distribution characteristics of the flank of the cutter tooth, and the simulation criteria of the upper and lower boundaries of the instantaneous friction and wear of the flank of the milling cutter are proposed. Based on the simulation criterion of instantaneous friction and wear, A method for calculating the cumulative friction and wear boundary of the flank of milling cutter teeth is proposed. By calculating the gray correlation value and the average relative error between the simulation and experimental cumulative friction and wear boundary curves, the shape between the simulation and experimental cumulative friction and wear boundary curves is analyzed. The degree of similarity with the position verifies the correctness of the calculation method of the cumulative friction and wear boundary on the flank of the milling cutter tooth.
所述一种高效铣刀累积摩擦磨损边界识别与验证方法,具体包括以下步骤:The method for identifying and verifying the cumulative friction and wear boundary of a high-efficiency milling cutter specifically includes the following steps:
步骤A:高效铣刀刀齿后刀面摩擦磨损边界测量与表征方法Step A: Measurement and characterization method of friction and wear boundary of high-efficiency milling cutter tooth flank
步骤A1振动作用下高效铣刀的铣削方式Step A1 Milling method of high-efficiency milling cutter under the action of vibration
在工件坐标系中,利用铣刀坐标系的偏置和作用下铣刀切削状态反映铣削振动对高效铣刀切削过程的影响,如图2和图3所示,图2和图3中变量参数解释如表1-1所示。In the workpiece coordinate system, the influence of milling vibration on the cutting process of the high-efficiency milling cutter is reflected by the offset of the milling cutter coordinate system and the cutting state of the milling cutter under the action, as shown in Figure 2 and Figure 3, the variable parameters in Figure 2 and Figure 3 The explanation is shown in Table 1-1.
表1-1振动作用下高效铣刀铣削方式及瞬时切削状态变量解释Table 1-1 High-efficiency milling cutter milling mode and instantaneous cutting state variable explanation under the action of vibration
步骤A2铣刀刀齿误差及其后刀面摩擦磨损测量Step A2 Milling cutter tooth error and flank friction and wear measurement
为研究铣刀刀齿误差对刀齿后刀面摩擦磨损边界的影响特特性,在铣刀坐标系中以铣刀刀齿最大外圆半径为基准度量刀齿径向误差,以铣刀刀齿最低点为基准度量刀齿轴向误差。如图4所示,图4中变量参数的解释如表1-2所示。In order to study the influence characteristics of the cutter tooth error on the friction and wear boundary of the cutter tooth flank, in the milling cutter coordinate system, the maximum outer radius of the cutter tooth is used as the benchmark to measure the cutter tooth radial error. The lowest point is the benchmark to measure the axial error of the cutter teeth. As shown in Figure 4, the explanation of the variable parameters in Figure 4 is shown in Table 1-2.
表1-2高级给铣刀结构及其刀齿误差测量方法变量解释Table 1-2 The structure of the advanced milling cutter and the variable explanation of the measurement method of the cutter tooth error
使刀齿安装基准面垂直于测量水平面,在刀齿后刀面在测量水平面的上的投影上,以过刀齿外侧刀尖点平行于切削刃中点的法矢量与沿切削刃中点的切矢量的反向交点Oj为原点,建立铣刀刀齿后刀面摩擦磨损边界测量坐标系Oj-UVS,轴U平行于投影面上刀齿切削刃中点切矢量方向,轴V平行于投影面上切削刃中点法矢量方向,轴S垂直于测量水平面。并在与刀齿切削刃中点切平面相平行的测量坐标系平行平面上度量后刀面摩擦磨损上下边界,如图5(a)所示。刀齿后刀面摩擦磨损深度表征方法如图5(b)所示。Make the installation reference plane of the cutter tooth perpendicular to the measuring horizontal plane, and on the projection of the flank face of the cutter tooth on the measuring horizontal plane, take the normal vector of the outer tip of the cutter tooth parallel to the midpoint of the cutting edge and the normal vector along the midpoint of the cutting edge. The reverse intersection point O j of the tangent vector is the origin, and the measurement coordinate system O j -UVS of the friction and wear boundary of the flank face of the milling cutter tooth is established. In the direction of the normal vector of the midpoint of the cutting edge on the projection plane, the axis S is perpendicular to the measurement horizontal plane. And measure the upper and lower boundaries of the friction and wear of the flank on the parallel plane of the measurement coordinate system parallel to the tangent plane of the midpoint of the cutting edge of the cutter teeth, as shown in Figure 5(a). Figure 5(b) shows the characterization method of the friction and wear depth on the flank of the cutter teeth.
图5中变量参数的解释如表1-3所示。The interpretation of the variable parameters in Figure 5 is shown in Table 1-3.
表1-3刀齿后刀面摩擦磨损边界测量方法变量解释Table 1-3 Explanation of variables for measuring the boundary of friction and wear on the flank of the cutter tooth
步骤A3铣削振动与刀齿误差对刀齿后刀面与工件瞬时接触关系的影响Step A3 Influence of milling vibration and cutter tooth error on the instantaneous contact relationship between cutter tooth flank and workpiece
利用铣削振动与刀齿误差影响下的工件坐标系、铣刀坐标系、刀齿坐标系和后刀面摩擦磨损测量坐标系之间的关系矩阵,反映铣削振动和刀齿误差对刀齿后刀面与工件瞬时接触关系的影响。无振动作用下铣刀坐标系与工件坐标系之间的关系矩阵如式(1-1)所示,振动与无振动作用下铣刀坐标系之间的关系矩阵如式(1-3)所示,刀齿坐标系与铣刀坐标系之间的关系矩阵如式(1-5)所示,刀齿坐标系与后刀面摩擦磨损测量坐标系之间的关系矩阵如式(1-7)。Using the relationship matrix between the workpiece coordinate system, the milling cutter coordinate system, the cutter coordinate system and the flank friction and wear measurement coordinate system under the influence of milling vibration and cutter tooth error, it reflects the impact of milling vibration and cutter tooth error on the cutter tooth flank. The influence of the instantaneous contact relationship between the surface and the workpiece. The relationship matrix between the milling cutter coordinate system and the workpiece coordinate system under the action of no vibration is shown in formula (1-1), and the relationship matrix between the milling cutter coordinate system under the action of vibration and no vibration is shown in formula (1-3) The relationship matrix between the cutter tooth coordinate system and the milling cutter coordinate system is shown in formula (1-5), and the relationship matrix between the cutter tooth coordinate system and the flank friction and wear measurement coordinate system is shown in formula (1-7) ).
[xg yg zg 1]T=M1[a b c 1]T (1-1)[x g y g z g 1] T =M 1 [abc 1] T (1-1)
[a b c 1]T=M2T1[av bv cv 1]T (1-3)[abc 1] T = M 2 T 1 [a v b v c v 1] T (1-3)
[a b c 1]T=M3T2[aq bq cq 1]T (1-5)[abc 1] T = M 3 T 2 [a q b q c q 1] T (1-5)
[aq bq cq 1]T=M4T4T3[U V S 1]T (1-7)[a q b q c q 1] T = M 4 T 4 T 3 [UVS 1] T (1-7)
其中,M1、M2、M3、M4为平移矩阵,T1、T2、T3、T4为旋转矩阵,M1中,x(t)、y(t)、z(t)为铣刀中心t时刻在工件坐标中的位置,T4中,为刀齿定位面与铣刀刀槽安装面之间的夹角,为π/2。Among them, M 1 , M 2 , M 3 , and M 4 are translation matrices; T 1 , T 2 , T 3 , and T 4 are rotation matrices; in M 1 , x(t), y(t), and z(t) is the position of the milling cutter center in workpiece coordinates at time t, in T 4 , is the angle between the positioning surface of the cutter teeth and the mounting surface of the milling cutter groove, is π/2.
步骤A4刀齿后刀面摩擦磨损边界的确定Step A4 Determination of the friction and wear boundary of the tooth flank
为完整表征出铣刀刀齿后刀面摩擦磨损上下边界的形状和位置信息,在超景深显微镜下测得的铣刀刀齿后刀面在磨损边界测量坐标系Oj-UV平面上的投影面,沿切削刃方向依次识别出刀齿后刀面摩擦磨损上下边界左右两侧端点和拐点,并在刀齿后刀面磨损边界测量坐标系中依次度量出端点和拐点沿V轴方向的坐标值,并以左右两侧端点和拐点为特征点构建铣刀刀齿后刀面摩擦磨损上下边界曲线。如图6所示,其中,V=V0(U)为铣刀刀齿后刀面原始切削刃轮廓曲线V=Vrq(U,L)(1≤q≤3)为铣削行程为L的铣刀刀齿q后刀面摩擦磨损上边界特征点曲线,V=Vhq(U,L)为铣削行程为L的铣刀刀齿q后刀面摩擦磨损下边界特征点曲线。In order to fully characterize the shape and position information of the upper and lower boundaries of the friction and wear of the milling cutter tooth flank, the projection of the milling cutter tooth flank on the wear boundary measurement coordinate system O j -UV plane measured under a super-depth microscope The end points and inflection points on the left and right sides of the upper and lower boundaries of the friction and wear of the flank flank of the cutter teeth and the inflection points are sequentially identified along the cutting edge direction, and the coordinates of the endpoints and the inflection points along the V-axis are measured in turn in the measurement coordinate system of the flank flank wear boundary of the cutter teeth. value, and the upper and lower boundary curves of friction and wear on the flank of the milling cutter tooth are constructed with the endpoints and inflection points on the left and right sides as feature points. As shown in Figure 6, where V=V 0 (U) is the original cutting edge profile curve of the flank face of the milling cutter tooth V=V rq (U, L) (1≤q≤3) is the milling stroke L The upper boundary characteristic point curve of the friction and wear of the flank face of the milling cutter tooth q, V=V hq (U, L) is the lower boundary characteristic point curve of the friction and wear of the milling cutter tooth q of the flank face with the milling stroke L.
该方法实施步骤如下,进行高效铣刀轴向分层铣削钛合金实验,获取高效铣刀刀齿后刀面摩擦磨损结果,铣削方案如表1-4所示。The implementation steps of the method are as follows. An experiment of high-efficiency milling cutter axial layer milling of titanium alloy is carried out, and the friction and wear results of the flank face of the high-efficiency milling cutter are obtained. The milling scheme is shown in Table 1-4.
表1-4高效铣刀铣削钛合金实验方案Table 1-4 Experimental scheme for milling titanium alloy with high-efficiency milling cutter
表1-4中,L为总的有效铣削行程为5m,Δcz i为刀齿i的轴向误差,Δcj i为刀齿i的径向误差。以铣刀刀齿最低点基准和铣刀刀齿刀尖点最大半径为基准计算铣刀各刀齿轴向和径向误差。In Table 1-4, L is the total effective milling stroke of 5m, Δc zi is the axial error of cutter i , and Δc ji is the radial error of cutter i . The axial and radial errors of each cutter tooth of the milling cutter are calculated based on the lowest point of the milling cutter tooth and the maximum radius of the milling cutter tooth tip point.
刀齿误差测量完成后,在三轴铣削加工中心上进行铣削钛合金实验,其中,工件长宽高为250×100×10mm,轴向分层铣削实验分20层完成,铣削实验过程中,采用加速度传感器检测靠近切削区域的工件表面处受切削力激励所产生的振动信号,并采用DH5922瞬态信号测试分析系统对所采集的振动信号进行数据分析。如图7所示。After the tool tooth error measurement is completed, the titanium alloy milling experiment is carried out on the three-axis milling machining center. The length, width and height of the workpiece are 250 × 100 × 10 mm. The acceleration sensor detects the vibration signal generated by the cutting force at the workpiece surface close to the cutting area, and uses the DH5922 transient signal test and analysis system to analyze the collected vibration signal. As shown in Figure 7.
实验结束后,使用超景深显微镜,依次对铣刀各刀齿后刀面摩擦磨损边界进行检测,如图8所示。After the experiment, the super-depth-of-field microscope was used to detect the friction and wear boundary of each tooth flank of the milling cutter in turn, as shown in Figure 8.
采用上述摩擦磨损上下边界曲线特征点识别方法依次提取图8中铣刀刀齿后刀面磨损磨损上下边界曲线特征点,如图9所示。The upper and lower boundary curve feature points identification method of friction and wear is used to extract the characteristic points of the upper and lower boundary curves of wear and wear on the flank of the milling cutter tooth in Figure 8 in turn, as shown in Figure 9.
由图9可知,该测量方法能完整表征出铣刀刀齿后刀面磨损上下边界形状和位置;受刀齿误差和铣削振动的影响,铣刀不同刀齿后刀面磨损上下边界分布呈现出明显的差异性。It can be seen from Figure 9 that this measurement method can completely characterize the shape and position of the upper and lower boundaries of the flank wear of the milling cutter teeth; affected by the cutter tooth error and milling vibration, the upper and lower boundary distributions of the flank wear of different teeth of the milling cutter show obvious differences.
步骤B:高效铣刀轴向分层铣削仿真模型及边界条件构建方法Step B: High-efficiency milling cutter axial layer milling simulation model and boundary condition construction method
步骤B1高效铣刀轴向分层铣削仿真模型整体结构设计Step B1 Overall structure design of high-efficiency milling cutter axial layered milling simulation model
将铣削方式、刀齿误差、铣削振动和不同层待加工表面差异性考虑为仿真边界条件,设计高效铣刀轴向分层铣削仿真模型整体结构,确定并解算仿真模型设计变量,如图10所示。Considering the milling method, cutter tooth error, milling vibration and the difference of the surface to be machined in different layers as simulation boundary conditions, design the overall structure of the simulation model of high-efficiency milling cutter axial layered milling, and determine and solve the design variables of the simulation model, as shown in Figure 10 shown.
图10(a)中,该仿真工件由若干个铣削单元构成,每一层铣削单元对应一种轴向分层铣削过程中的工件形态,铣削单元之间为矩形凹槽结构,不同层待加工表面轮廓存在差异性但处于同一水平高度上,Li (i+1)为第i-1层铣削单元与第i层铣削单元之间的距离,H'为矩形凹槽的高度。In Figure 10(a), the simulated workpiece is composed of several milling units. Each layer of milling units corresponds to a workpiece shape in an axial layered milling process. The milling units have a rectangular groove structure, and different layers are to be processed. The surface profiles are different but at the same level, Li (i+1) is the distance between the i-1th layer milling unit and the i-th layer milling unit, and H' is the height of the rectangular groove.
图10(b)中,Oi为铣刀刀齿底刃切削刃所在圆圆心,r0为铣刀刀齿切削刃所在圆半径,dm为铣刀轴线C1C2沿切削宽度方向到工件之间的距离,Jd为铣刀刀齿切削刃上参与切削的最低点,Jh为刀齿侧刃与待加工表面最高接触点,Δrz max和Δrj max分别为铣刀刀齿最大轴向和径向误差;Hi表示第i(i≥2)层铣削单元的高度,其中第一层铣削单元横截面与工件的初始横截面尺寸相同。In Figure 10(b), O i is the center of the circle where the cutting edge of the milling cutter tooth bottom edge is located, r 0 is the radius of the circle where the milling cutter tooth cutting edge is located, and d m is the milling cutter axis C 1 C 2 along the cutting width direction to The distance between the workpieces, J d is the lowest point on the cutting edge of the milling cutter tooth that participates in cutting, J h is the highest contact point between the side edge of the cutter tooth and the surface to be machined, Δr z max and Δr j max are the milling cutter teeth, respectively Maximum axial and radial errors; H i represents the height of the i-th (i≥2) layer of milling units, where the cross-section of the first layer of milling units is the same as the initial cross-sectional dimension of the workpiece.
该仿真工件的控制参数有:Li (i+1)、H'、Li-1 i、Hi、W、L0、dm、ap、Δrz max、Δrj max、r0、r、α、Jd、Jh。Li (i+1)的求解方法如下:The control parameters of the simulated workpiece are: L i (i+1) , H', L i-1 i , H i , W, L 0 , d m , a p , Δr z max , Δr j max , r 0 , r, a, Jd , Jh . The solution method of L i (i+1) is as follows:
其中,Ti i+1为实际加工过程中铣刀从P3i到P2(i+1)的时间。Among them, T i i+1 is the time from P 3i to P 2(i+1) of the milling cutter in the actual machining process.
Hi=H+(i-1)ap(i≥2) (2-3)H i =H+(i-1)a p (i≥2) (2-3)
步骤B2仿真模型铣削单元构建Step B2 Simulation Model Milling Unit Construction
(1)根据步骤A中图2和图3所示,求解铣削过程中t时刻铣刀沿进给速度方向和切削宽度方向的瞬时偏转角δ1(t)、δ2(t),如式(2-4)、(2-5)。(1) According to Fig. 2 and Fig. 3 in step A, solve the instantaneous deflection angles δ 1 (t) and δ 2 (t) of the milling cutter along the feed speed direction and the cutting width direction at time t during the milling process, as shown in the formula (2-4), (2-5).
求解振动作用下铣刀中心点在工件坐标系的轨迹方程fv=(xgv,ygv,zgv,t),如式(2-6)所示。Solve the trajectory equation f v = (x gv , y gv , z gv , t) of the milling cutter center point in the workpiece coordinate system under the action of vibration, as shown in formula (2-6).
式(2-6)中,Ax(t)、Ay(t)、Az(t)通过对实验过程中测得的振动加速度信号进行二次积分,再使用Matlab中的8次Sun of sin Function函数拟合获得。te为铣刀从P4i到P1(i+1)所用的时间,如式(2-7)所示,Ti为铣刀从P1(i+1)运动到P1(i+1)的时间段,如式(2-8)所示。In formula (2-6), A x (t), A y (t), and A z (t) are obtained by quadratic integration of the vibration acceleration signal measured during the experiment, and then using the 8th Sun of Matlab obtained by fitting the sin Function function. t e is the time it takes for the milling cutter to move from P 4i to P 1(i+1) , as shown in formula (2-7), T i is the milling cutter moving from P 1(i+1) to P 1(i+ 1) , as shown in formula (2-8).
Ti=[Tsi,Tei]=[(i-1)(Tc+te),iTc+(i-1)te] (2-8)T i =[T si ,T ei ]=[(i-1)(T c +t e ),iT c +(i-1)t e ] (2-8)
其中,Tsi为铣刀沿xg轴方向第i次进给的刀路起点时刻,Tei为铣刀沿xg轴方向第i次进给的抬刀时刻;Tc铣刀从P1i运动至P4i所用的时间,如式(2-9)。Among them, T si is the starting point of the tool path of the i-th feed of the milling cutter along the x g -axis direction, T ei is the tool-lifting time of the i-th feed of the milling cutter along the x g -axis direction; T c milling cutter starts from P 1i The time it takes to move to P 4i is shown in formula (2-9).
(2)识别出第i(i≥2)层铣削振动信号中的切入切出时刻,如图11所示,根据式(2-6)求解第i层切削段Tm的铣刀运动轨迹方程fv=(xgv,ygv,zgv,t),如图12所示。(2) Identify the cut-in and cut-out time in the milling vibration signal of the i-th layer (i≥2), as shown in Figure 11, and solve the milling cutter motion trajectory equation of the cutting section T m of the i-th layer according to equation (2-6). f v = (x gv , y gv , z gv , t), as shown in FIG. 12 .
(3)根据式(2-4)、式(2-5),使用matlab求解i层铣刀瞬时切削姿态角,结果如图13所示。(3) According to equations (2-4) and (2-5), use matlab to solve the instantaneous cutting attitude angle of the i-layer milling cutter, and the results are shown in Figure 13.
(4)如图14所示,根据图10(b)中铣削单元横截面模型,构建第i层铣削单元横截面,将NJh段沿铣刀轨迹曲线进行扫掠形成铣削单元待加工表面,扫掠过程中引入第i层铣刀瞬时切削姿态角,最后依次构建铣削单元其他表面并形成实体,图14中变量的参数解释如表2-1所示。(4) As shown in Figure 14, according to the cross-sectional model of the milling unit in Figure 10(b), the cross-section of the i-th layer of milling unit is constructed, and the NJ h segment is swept along the milling cutter trajectory curve to form the surface to be machined by the milling unit, In the sweeping process, the instantaneous cutting attitude angle of the i-th layer milling cutter is introduced, and finally other surfaces of the milling unit are constructed in turn to form entities. The parameter explanations of the variables in Figure 14 are shown in Table 2-1.
表2-1第i层铣削铣削单元变量解算Table 2-1 Calculation of variables for the i-th layer milling and milling unit
步骤B3振动作用下高效铣刀轴向分层铣削仿真运动轨迹解算Step B3 Simulation motion trajectory calculation of high-efficiency milling cutter axial layered milling under the action of vibration
在仿真过程中,为引入铣削过程中铣刀振动位移,根据轴向分层铣削实验中振动作用下的铣刀运动轨迹,基于不改变进给速度的前提下,确定铣刀运动轨迹点的离散时间间隔,对仿真过程中振动作用下的铣刀运动轨迹点进行解算,步骤如下:In the simulation process, in order to introduce the vibration displacement of the milling cutter during the milling process, according to the movement trajectory of the milling cutter under the action of vibration in the axial layered milling experiment, and based on the premise of not changing the feed rate, the discrete points of the movement trajectory of the milling cutter are determined. The time interval is used to calculate the movement trajectory points of the milling cutter under the action of vibration during the simulation process. The steps are as follows:
(1)确定轨迹点离散时间间隔,为了使得仿真结果更加准确,以工件最小网格单元最小单元尺寸Lmin的1/3为仿真步距,td如式(2-10)。(1) Determine the discrete time interval of the trajectory points. In order to make the simulation results more accurate, 1/3 of the minimum unit size L min of the minimum grid unit of the workpiece is used as the simulation step, and t d is as formula (2-10).
(2)对铣刀任意第i(i≥2)层切入到切出的坐标值进行重新解算,如式(2-11)所示。(2) Recalculate the coordinate values from the cut-in to the cut-out at any i-th layer (i≥2) of the milling cutter, as shown in formula (2-11).
(3)补齐任意第i层切出点到第i+1层切入点之间铣刀轨迹点坐标,求解第i层切出点到第i+1层切入点之间的时间步数如式(2-12),第i层切出点到第i+1层切入点之间相邻轨迹点步距如式(2-13),相邻两层轨迹点之间坐标求解方法如式(2-14)。(3) Fill up the coordinates of the milling cutter trajectory point between the cut-out point of any i-th layer and the cut-in point of the i+1-th layer, and solve the time steps between the cut-out point of the i-th layer and the cut-in point of the i+1-th layer as follows: Formula (2-12), the step distance of adjacent trajectory points between the cut-out point of the i-th layer and the cut-in point of the i+1-th layer is as shown in formula (2-13), and the coordinate solution method between the adjacent two-layer track points is as formula (2-14).
步骤C:刀齿后刀面瞬时摩擦磨损边界识别方法Step C: Recognition method of instantaneous friction and wear boundary of tooth flank
步骤C1高效铣刀轴向分层铣削热力耦合场仿真Step C1 High-efficiency milling cutter axial layer milling thermal coupling field simulation
为获取铣刀刀齿后刀面瞬时摩擦磨损仿真边界,利用步骤A中铣削加工实验参数,根据步骤B中设计的高效铣刀轴向分层铣削仿真模型并进行热力耦合场仿真,仿真所采用的Johnson-Cook本构模型如式(3-1)所示,仿真铣削加工参数采用表1-4中的高效铣刀铣削钛合金实验加工参数。In order to obtain the simulation boundary of instantaneous friction and wear on the flank of the milling cutter tooth, using the milling experimental parameters in step A, according to the high-efficiency milling cutter axial layered milling simulation model designed in step B, the thermal coupling field simulation was carried out. The Johnson-Cook constitutive model is shown in formula (3-1), and the simulated milling processing parameters use the high-efficiency milling cutter milling titanium alloy experimental processing parameters in Table 1-4.
式中:σ—等效流动应力;A,B,n,c,m—屈服应力强度,应变强化常数,应变强化指数,应变率强化参数,温度应变率灵敏度;T,Tr,Tm—实时温度、室温,取值25℃、材料融化温度/℃;ε,—材料等效应变,等效应变率/s-1,参考应变率取值1s-1;其中,高效铣刀基体材料为硬质合金WC,刀具涂层材料为TiN,厚度为5μm,屈服强度σs为5.9GPa,TC4钛合金材料J-C本构参数如表3-1所示。铣削仿真过程如图15所示。In the formula: σ—equivalent flow stress; A, B, n, c, m—yield stress intensity, strain hardening constant, strain hardening index, strain rate hardening parameter, temperature strain rate sensitivity; T, T r , T m — Real-time temperature, room temperature, value 25℃, material melting temperature/℃; ε, —Material equivalent strain, equivalent strain rate/s -1 , the reference strain rate is 1s -1 ; among them, the base material of the high-efficiency milling cutter is cemented carbide WC, the tool coating material is TiN, the thickness is 5 μm, and the yield strength σ s is 5.9GPa, and the JC constitutive parameters of TC4 titanium alloy material are shown in Table 3-1. The milling simulation process is shown in Figure 15.
表3-1 TC4钛合金J-C本构参数表Table 3-1 TC4 titanium alloy J-C constitutive parameter table
仿真结束后,提取铣刀刀齿后刀面任意瞬时切削等效应力场、等效应变场和温度场,如图16所示。After the simulation, the equivalent stress field, equivalent strain field and temperature field of any instantaneous cutting on the flank of the milling cutter tooth are extracted, as shown in Figure 16.
步骤C2刀齿后刀面不同位置处瞬时节点等效应力、等效应变、温度变化曲线Step C2 Instantaneous node equivalent stress, equivalent strain and temperature change curve at different positions on the flank of the tooth
提取铣刀刀齿后刀面不同位置处节点等效应力、等效应变和温度数值,构建刀齿后刀面等效应力、等效应变、温度分布曲面,平行于刀齿后刀面投影面上切削刃中点法矢量方向过切削刃上节点依次作若干垂直于分布曲面的截平面,如图17所示,构建后刀面不同位置处截平面上节点等效应力分布曲线、并识别出截平面上等效应力大于且最接近于屈服强度σs的临界节点,如图18所示,构建后刀面不同位置处截平面上节点等效应变、温度分布曲线,以截平面B为例,如图19所示。Extract the equivalent stress, equivalent strain and temperature values of nodes at different positions on the flank of the milling cutter, and construct the equivalent stress, equivalent strain, and temperature distribution surfaces of the flank, which are parallel to the projection surface of the flank. The normal vector direction of the midpoint of the upper cutting edge passes through the upper node of the cutting edge to make several section planes perpendicular to the distribution surface in turn, as shown in Figure 17, to construct the equivalent stress distribution curves of the nodes on the section plane at different positions of the flank face, and identify them. The equivalent stress on the section plane is greater than and closest to the critical node of the yield strength σ s , as shown in Figure 18, the equivalent strain and temperature distribution curves of the section plane on the section plane at different positions of the flank are constructed, taking section plane B as an example , as shown in Figure 19.
步骤C3刀齿后刀面不同位置处应力、应变、温度变化率曲线Step C3 Curves of stress, strain and temperature change rate at different positions on the flank of the tooth
在高效铣刀切削过程中,刀齿后刀面与工件过渡表面之间存在剪切与摩擦,使得在接触区域,刀齿内部形成应力集中,并向非接触区域扩散分布,其接触与未接触边界节点处,应力数值急剧下降;同时,在剪切摩擦过程中产生的热量使刀齿后刀面接触表面温度急剧升高,产生的热量在刀齿后刀面处向非接触区域传递,接触与未接触边界节点处,温度数值急剧下降,使得刀齿后刀面上产生非常大的温度梯度;在铣削过程中,刀齿后刀面与工件过渡表面之间相互作用发生挤压变形,使得刀齿后刀面与工件过渡表面接触边界节点处产生应变集中,其未接触边界处产生微小应变;基于上述分析,刀齿后刀面瞬时摩擦磨损边界节点处应力、应变、温度数值发生突变,因此,通过一元高次多项式依次拟合各截面处节点温度、等效应力、等效应变分布曲线,采用式(3-2)、式(3-3)、式(3-4)求导得节点变化率分布曲线,以截面B为例,分别识别出后刀面各截面位置处温度、应力、应变变化率发生突变的节点拐点,如图20所示。In the cutting process of high-efficiency milling cutters, there is shear and friction between the flank of the cutter teeth and the transition surface of the workpiece, so that the stress concentration is formed inside the cutter teeth in the contact area and spreads to the non-contact area. At the boundary node, the stress value drops sharply; at the same time, the heat generated in the process of shear friction causes the temperature of the contact surface of the tooth flank to rise sharply, and the generated heat is transferred to the non-contact area at the flank of the tooth, and the contact At the non-contact boundary node, the temperature value drops sharply, resulting in a very large temperature gradient on the flank of the tooth; during the milling process, the interaction between the flank of the tooth and the transition surface of the workpiece produces extrusion deformation, making the Strain concentration occurs at the contact boundary node between the tooth flank and the transition surface of the workpiece, and a small strain occurs at the non-contact boundary; Therefore, the node temperature, equivalent stress, and equivalent strain distribution curves at each section are fitted sequentially by a high-order polynomial of one variable, and the derivation is obtained by using formula (3-2), formula (3-3), and formula (3-4). The node change rate distribution curve, taking section B as an example, identifies the node inflection points where the change rate of temperature, stress and strain changes abruptly at each section position of the flank face, as shown in Figure 20.
步骤C4刀齿后刀面瞬时摩擦磨损仿真上边界识别Step C4 Recognition of upper boundary of instantaneous friction and wear simulation of tooth flank
基于材料的力学性能,当刀齿切削刃上的等效应力大于屈服强度时,切削刃上材料会发生剥落被切屑带走,导致刃形发生改变形成刀齿后刀面瞬时摩擦磨损上边界,基于上述分析,以刀齿材料屈服强度σs为铣刀刀齿后刀面瞬时摩擦磨损仿真上边界识别判剧,如式(3-5)所示,识别出所有截平面上等效应力大于且最近于屈服强度的临界特征节点,在后刀面摩擦磨损测量坐标系中,截平面上特征节点所围成的边界即为刀齿后刀面瞬时摩擦磨损上边界,如图21所示。Based on the mechanical properties of the material, when the equivalent stress on the cutting edge of the tooth is greater than the yield strength, the material on the cutting edge will peel off and be taken away by the chips, resulting in the change of the blade shape and the upper boundary of instantaneous friction and wear on the flank of the tooth. Based on the above analysis, the yield strength σ s of the cutter tooth material is used as the upper boundary of the instantaneous friction and wear simulation of the cutter tooth flank to identify and judge the drama. And the critical feature node closest to the yield strength, in the flank friction and wear measurement coordinate system, the boundary surrounded by the feature nodes on the section plane is the upper boundary of the instantaneous friction and wear of the tooth flank, as shown in Figure 21.
σ≥σs (3-5)σ≥σ s (3-5)
步骤C5刀齿后刀面瞬时摩擦磨损仿真下边界识别Step C5 Identification of the lower boundary of the instantaneous friction and wear simulation of the flank of the cutter tooth
在瞬时切削过程中,铣刀刀齿后刀面与工件过渡表面发生剪切、挤压和摩擦形成铣刀刀齿后刀面瞬摩擦磨损下边界,其刀齿后刀面与工件过渡表面接触与未接触区域之间的温度、等效应力、等效应变会存在明显差异性,基于以上分析,提出铣刀刀齿后刀面瞬时摩擦磨损仿真下边界识别方法,采用上述方法依次识别出刀齿后刀面各截平面处温度、等效应力、等效应变发生明显突变的节点,刀齿后刀面上节点温度、等效应力、等效应变突变曲线如图22所示,铣刀刀齿后刀面瞬时摩擦磨损仿真下边界曲线解算方法如式(3-6)、式(3-7)所示。During the instantaneous cutting process, the flank of the milling cutter tooth and the transition surface of the workpiece are sheared, squeezed and rubbed to form the lower boundary of instantaneous friction and wear of the flank of the milling cutter tooth, and the flank of the cutter tooth is in contact with the transition surface of the workpiece. There will be obvious differences in temperature, equivalent stress, and equivalent strain between the non-contact area and the non-contact area. Based on the above analysis, a method for identifying the lower boundary of the instantaneous friction and wear simulation of the flank of the milling cutter tooth flank is proposed. The above method is used to identify the cutter in turn. The nodes where the temperature, equivalent stress, and equivalent strain of each section plane of the tooth flank have abrupt changes are shown in Figure 22. The calculation method of the lower boundary curve of the instantaneous friction and wear simulation of the flank flank is shown in equations (3-6) and (3-7).
Vm(Ui,t)=max{VT(Ui,t),Vσ(Ui,t),Vε(Ui,t)}(0≤i≤10.7) (3-6)V m (U i ,t)=max{V T (U i ,t),V σ (U i ,t),V ε (U i ,t)}(0≤i≤10.7) (3-6)
Vm(U,t)={Vm(U1,t),Vm(U2,t)...Vm(Ui,t)} (3-7)V m (U,t)={V m (U 1 ,t),V m (U 2 ,t)...V m (U i ,t)} (3-7)
如图22中,V=V0(U)为铣刀原始轮廓边界,V=VT(U,t)、V=Vσ(U,t)、Vε(U,t)分别为t时刻铣刀刀齿后刀面温度、应力、应变突变节点曲线,式(3-6)中,VT(Ui,t)、Vσ(Ui,t)、Vε(Ui,t)分别为t时刻刀齿后刀面温度、应力、应变突变节点曲线Ui位置处V向坐标值。式(3-7)中Vm(Ui,t)为t时刻铣刀刀齿后刀面瞬时摩擦磨损下边界曲线Ui位置处V向坐标值,Vm(U,t)为t时刻铣刀刀齿后刀面瞬时摩擦磨损下边界曲线。As shown in Figure 22, V=V 0 (U) is the original contour boundary of the milling cutter, V=V T (U, t), V=V σ (U, t), V ε (U, t) are the time t respectively Temperature, stress, strain sudden change node curve of milling cutter tooth flank, in formula (3-6), V T (U i ,t), V σ (U i ,t), V ε (U i ,t) are the V-coordinate values of the tooth flank temperature, stress and strain mutation node curve U i at time t, respectively. In Equation (3-7), V m (U i ,t) is the coordinate value in the V direction at the lower boundary curve U i of the instantaneous friction and wear of the milling cutter tooth flank at time t, and V m (U, t) is the time t The lower boundary curve of instantaneous friction and wear on the flank of milling cutter teeth.
步骤D:刀齿后刀面累积摩擦磨损边界解算及验证方法Step D: Boundary solution and verification method of cumulative friction and wear on the flank of the cutter tooth
步骤D1铣刀刀齿后刀面累积摩擦磨损边界解算Step D1 Calculate the boundary of cumulative friction and wear on the flank of milling cutter teeth
在铣削过程中,受铣削振动和刀齿误差的影响,铣刀刀齿后刀面与工件过渡表面之间瞬时接触边界随时间发生动态变化,铣刀刀齿后刀面累积摩擦磨损边界由若干瞬时摩擦磨损边界累积形成,基于以上分析,以铣刀刀齿后刀面累积摩擦磨损仿真下边界为例,提取不同铣削时刻刀齿后刀面瞬时摩擦磨损下边界曲线,如图23所示,提出铣刀刀齿后刀面累积摩擦磨损边界解算方法,如式(4-1)和式(4-2)所示。During the milling process, affected by the milling vibration and cutter tooth error, the instantaneous contact boundary between the milling cutter tooth flank and the transition surface of the workpiece changes dynamically with time, and the cumulative friction and wear boundary of the milling cutter tooth flank is determined by several The instantaneous friction and wear boundary is cumulatively formed. Based on the above analysis, taking the lower boundary of the simulation of the cumulative friction and wear of the flank of the milling cutter as an example, the lower boundary curve of the instantaneous friction and wear of the flank of the cutter tooth at different milling times is extracted, as shown in Figure 23. The calculation method of the cumulative friction and wear boundary of the flank of the milling cutter tooth flank is proposed, as shown in equations (4-1) and (4-2).
Va(Ui,Δt)=max{Vm(Ui,t1),Vm(Ui,t2)...Vm(Ui,tn)}Δt=tn-t1 (4-1)V a (U i ,Δt)=max{V m (U i ,t 1 ),V m (U i ,t 2 )...V m (U i ,t n )}Δt=t n -t 1 (4-1)
Va(U,Δt)={Va(U1,Δt),Va(U2,Δt)...Va(Ui,Δt)} (4-2)V a (U,Δt)={V a (U 1 ,Δt),V a (U 2 ,Δt)...V a (U i ,Δt)} (4-2)
式(4-1)中,Va(Ui,Δt)为Δt时间段内铣刀刀齿后刀面累积摩擦磨损边界曲线在Ui位置处V向坐标值,式(4-2)中,Va(U,Δt)为Δt时间段内铣刀刀齿后刀面累积摩擦磨损边界曲线。In formula (4-1), V a (U i ,Δt) is the V coordinate value of the cumulative friction and wear boundary curve of the flank of the milling cutter tooth flank at the position U i in the Δt time period, in formula (4-2) , V a (U, Δt) is the cumulative friction and wear boundary curve of the flank of the milling cutter tooth in the Δt time period.
步骤D2刀齿后刀面累积摩擦磨损下边界仿真与实验对比方案Step D2 The lower boundary simulation and experiment comparison scheme of the cumulative friction and wear of the flank of the cutter tooth
为验证高效铣刀累积摩擦磨损边界解算方法的正确性,先采用步骤C4和步骤C5中的铣刀刀齿后刀面瞬时摩擦磨损上下边界仿真识别判据,提取步骤C1中热力耦合场仿真中0.5m、1m、1.5m、2m、2.5m、3m、3.5m、4m、4.5m、5m时铣刀各刀齿刀面瞬时摩擦磨损仿真上下边界,以刀齿一为例,如图24所示。In order to verify the correctness of the calculation method of the cumulative friction and wear boundary of the high-efficiency milling cutter, the upper and lower boundary simulation and identification criteria of the instantaneous friction and wear of the milling cutter tooth flank in steps C4 and C5 are firstly used, and the thermal coupling field simulation in step C1 is extracted. At 0.5m, 1m, 1.5m, 2m, 2.5m, 3m, 3.5m, 4m, 4.5m, and 5m, the upper and lower boundaries of instantaneous friction and wear of each cutter face of the milling cutter are simulated. Take
采用步骤D1中铣刀刀齿后刀面累积摩擦磨损边界曲线解算方法,构建铣刀刀齿后刀面累积摩擦磨损仿真上下边界曲线,与-步骤A4中铣削行程为5m的铣刀刀齿后刀面实验摩擦磨损上下边界曲线进行对比,以刀齿一为例,如图25所示。Using the calculation method of the cumulative friction and wear boundary curve of the flank of the milling cutter tooth in step D1, the upper and lower boundary curves of the cumulative friction and wear simulation of the flank face of the milling cutter tooth are constructed, and the milling cutter tooth with a milling stroke of 5m in step A4 The upper and lower boundary curves of the friction and wear of the flank face experiment are compared, taking the first tooth as an example, as shown in Figure 25.
步骤D3仿真与实验累积摩擦磨损边界形状与位置相似度分析Step D3 Simulation and Experiment Cumulative Friction and Wear Boundary Shape and Position Similarity Analysis
计算仿真与实验累积摩擦磨损边界曲线之间的灰色关联度值和平均相对误差,分析仿真与实验磨损边界曲线之间形状相似度和位置重合度,验证累积摩擦磨损边界解算方法的正确性,灰色关联分析方法步骤如下:Calculate the grey correlation value and average relative error between the simulation and experimental cumulative friction and wear boundary curves, analyze the shape similarity and position coincidence between the simulation and experimental wear boundary curves, and verify the correctness of the cumulative friction and wear boundary solution method. The steps of the grey relational analysis method are as follows:
首先分别以实验和仿真刀齿后刀面磨损边界特征点V轴坐标值作为参考序列和比较序列,如式(4-3)、(4-4)。First, the V-axis coordinate values of the experimental and simulated tool tooth flank wear boundary feature points are used as the reference sequence and the comparison sequence, such as formulas (4-3) and (4-4).
Vq=(Vq(1),Vq(2),...,Vq(m),...,Vq(n)) (4-3) Vq =( Vq (1), Vq (2),..., Vq (m),..., Vq (n)) (4-3)
Vq*=(Vq(1)*,Vq(2)*,...,Vq(m)*,...,Vq(n)*) (4-4) Vq *=( Vq (1)*, Vq (2)*,..., Vq (m)*,..., Vq (n)*) (4-4)
式(4-3)、(4-4)中Vq(m)和Vq(m)*分别为实验和仿真第q(0≤q≤3)个刀齿后刀面磨损边界曲线上沿切削刃方向第m个特征点V向坐标值Vm。In formulas (4-3) and (4-4), V q (m) and V q (m)* are the upper edge of the flank wear boundary curve of the qth (0≤q≤3) tooth in the experiment and simulation, respectively. The coordinate value V m of the m-th feature point in the cutting edge direction in the V direction.
仿真与实验刀齿后刀面磨损曲线平均相对误差计算方法如式(4-5)。The calculation method of the average relative error of the flank wear curve of the simulated and experimental tool teeth is shown in formula (4-5).
式(4-5)中,δq为仿真与实验第q个刀齿后刀面磨损边界曲线相对误差。In formula (4-5), δ q is the relative error of the flank wear boundary curve of the qth tooth in the simulation and experiment.
采用以上关联度和平均相对误差计算方法,仿真与实验刀齿后刀面磨损上下边界曲线关联度及相对误差如表4-1所示,其中γ(Vq,Vq*)(1≤q≤3)表示刀齿q仿真与实验累积摩擦磨损边界曲线之间的灰色关联度,δq表示刀齿q仿真与实验累积磨损边界曲线之间的平均相对误差。Using the above correlation degree and average relative error calculation method, the correlation degree and relative error of the upper and lower boundary curves of simulated and experimental tool tooth flank wear are shown in Table 4-1, where γ(V q ,V q *)(1≤q ≤3) represents the grey correlation degree between the tool tooth q simulation and the experimental cumulative friction and wear boundary curve, and δ q represents the average relative error between the tool tooth q simulation and the experimental cumulative wear boundary curve.
表4-1仿真与实验刀齿后刀面磨损上下边界曲线关联度及相对误差Table 4-1 Correlation degree and relative error of upper and lower boundary curves of flank wear between simulation and experiment
由表4-1,铣刀各刀齿后刀面仿真与实验累积摩擦磨损上下边界曲线关联度值均在0.75以上,平均相对误差均小于20%,其结果表明,铣刀各刀齿后刀面仿真与实验累积摩擦磨损边界曲线形状和位置具有很高的相似度,高效铣刀轴向分层铣削仿真模型设计方法的有效性,同时也验证了铣刀刀齿后刀面累积摩擦磨损边界解算方法的正确性。From Table 4-1, the correlation degree between the simulation of the flank of each tooth of the milling cutter and the upper and lower boundary curves of the cumulative friction and wear of the milling cutter are all above 0.75, and the average relative error is less than 20%. The results show that the flank of each tooth of the milling cutter The surface simulation and the experimental cumulative friction and wear boundary curve shape and position have high similarity, the effectiveness of the high-efficiency milling cutter axial layer milling simulation model design method, and also verifies the milling cutter tooth flank cumulative friction and wear boundary The correctness of the solution method.
以上所述,仅为本发明较佳的具体实施方式,这些具体实施方式都是基于本发明整体构思下的不同实现方式,而且本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。The above are only preferred specific embodiments of the present invention, and these specific embodiments are based on different implementations under the overall concept of the present invention, and the protection scope of the present invention is not limited to this. Anyone familiar with the technical field Changes or substitutions that can be easily conceived by a skilled person within the technical scope disclosed by the present invention shall be covered within the protection scope of the present invention. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011439558.1A CN112372372B (en) | 2020-12-07 | 2020-12-07 | Efficient milling cutter accumulated friction wear boundary identification and verification method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011439558.1A CN112372372B (en) | 2020-12-07 | 2020-12-07 | Efficient milling cutter accumulated friction wear boundary identification and verification method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112372372A CN112372372A (en) | 2021-02-19 |
CN112372372B true CN112372372B (en) | 2022-09-02 |
Family
ID=74590798
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011439558.1A Active CN112372372B (en) | 2020-12-07 | 2020-12-07 | Efficient milling cutter accumulated friction wear boundary identification and verification method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112372372B (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114896753B (en) * | 2021-12-09 | 2024-06-07 | 哈尔滨理工大学 | Identification method of instantaneous friction energy consumption distribution characteristics of the flank surface of energy-efficient milling cutters |
CN114580090A (en) * | 2021-12-28 | 2022-06-03 | 哈尔滨理工大学 | Dynamic characteristic resolving method for friction coefficient of rear cutter face of cutter tooth pair of square shoulder milling cutter |
CN114429064B (en) * | 2022-01-05 | 2024-10-15 | 哈尔滨理工大学 | Identification method for fractal characteristics of friction boundary of rear cutter surface of cutter tooth of high-energy-efficiency milling cutter |
CN114769692B (en) * | 2022-04-19 | 2023-01-03 | 哈尔滨理工大学 | Efficient cross-scale identification method for friction wear boundary of rear cutter face of milling cutter |
CN114952420B (en) * | 2022-05-23 | 2023-04-11 | 大连理工大学 | Wear prediction calculation method for milling CFRP (carbon fiber reinforced plastics) of ball end mill based on path dispersion |
CN116061004B (en) * | 2022-10-14 | 2023-08-08 | 哈尔滨理工大学 | Identification method for change characteristics of transient cutting energy efficiency distribution of milling cutter under vibration effect |
CN116738730B (en) * | 2023-06-17 | 2023-12-05 | 哈尔滨理工大学 | Method for solving propagation and attenuation characteristics of friction stress wave of rear cutter surface of efficient milling cutter |
CN116690313B (en) * | 2023-08-01 | 2023-12-12 | 成都飞机工业(集团)有限责任公司 | Failure monitoring method for machining cutter of web plate of aircraft structural member |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5251144A (en) * | 1991-04-18 | 1993-10-05 | Texas Instruments Incorporated | System and method utilizing a real time expert system for tool life prediction and tool wear diagnosis |
CN109434564A (en) * | 2018-12-21 | 2019-03-08 | 哈尔滨理工大学 | A kind of cutter wear state monitoring method based on deep neural network |
CN109940461A (en) * | 2019-04-10 | 2019-06-28 | 哈尔滨理工大学 | A method for detecting the wear characteristics of the flank face of a high-feed milling cutter |
CN109940460A (en) * | 2019-04-10 | 2019-06-28 | 哈尔滨理工大学 | A detection method for the geometric error distribution characteristics of milling machined surface under the action of vibration |
CN110103079A (en) * | 2019-06-17 | 2019-08-09 | 中国科学院合肥物质科学研究院 | The on-line monitoring method of tool wear in a kind of micro- milling process |
CN110704969A (en) * | 2019-09-24 | 2020-01-17 | 哈尔滨理工大学 | Analytical Model Construction and Verification Method for Revealing the Wear Variation of High-feed Milling Cutter |
CN111002104A (en) * | 2019-12-16 | 2020-04-14 | 哈尔滨理工大学 | A method for detecting and solving the friction and wear boundary of the flank face of a high-feed milling cutter |
CN111644900A (en) * | 2020-05-21 | 2020-09-11 | 西安交通大学 | Tool damage real-time monitoring method based on spindle vibration characteristic fusion |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE50001492D1 (en) * | 2000-05-15 | 2003-04-24 | Prometec Gmbh | Method and device for monitoring the state of wear of a tool |
KR100517880B1 (en) * | 2003-11-20 | 2005-09-28 | 학교법인 포항공과대학교 | Off-line feed rate scheduling for reduction of machining time and enhancement of machining accuracy in cnc machining |
WO2015079164A1 (en) * | 2013-11-29 | 2015-06-04 | Snecma | Method and device for determining the wear of a cutting tool flank |
WO2020075191A1 (en) * | 2018-10-12 | 2020-04-16 | Tata Consultancy Services Limited | Method and system for monitoring tool wear to estimate rul of tool in machining |
-
2020
- 2020-12-07 CN CN202011439558.1A patent/CN112372372B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5251144A (en) * | 1991-04-18 | 1993-10-05 | Texas Instruments Incorporated | System and method utilizing a real time expert system for tool life prediction and tool wear diagnosis |
CN109434564A (en) * | 2018-12-21 | 2019-03-08 | 哈尔滨理工大学 | A kind of cutter wear state monitoring method based on deep neural network |
CN109940461A (en) * | 2019-04-10 | 2019-06-28 | 哈尔滨理工大学 | A method for detecting the wear characteristics of the flank face of a high-feed milling cutter |
CN109940460A (en) * | 2019-04-10 | 2019-06-28 | 哈尔滨理工大学 | A detection method for the geometric error distribution characteristics of milling machined surface under the action of vibration |
CN110103079A (en) * | 2019-06-17 | 2019-08-09 | 中国科学院合肥物质科学研究院 | The on-line monitoring method of tool wear in a kind of micro- milling process |
CN110704969A (en) * | 2019-09-24 | 2020-01-17 | 哈尔滨理工大学 | Analytical Model Construction and Verification Method for Revealing the Wear Variation of High-feed Milling Cutter |
CN111002104A (en) * | 2019-12-16 | 2020-04-14 | 哈尔滨理工大学 | A method for detecting and solving the friction and wear boundary of the flank face of a high-feed milling cutter |
CN111644900A (en) * | 2020-05-21 | 2020-09-11 | 西安交通大学 | Tool damage real-time monitoring method based on spindle vibration characteristic fusion |
Non-Patent Citations (1)
Title |
---|
高进给铣刀刀齿后刀面磨损差异性研究;左林晗;《中国优秀硕士学位论文全文数据库 工程科技I辑(月刊)》;20190815(第08期);7-65 * |
Also Published As
Publication number | Publication date |
---|---|
CN112372372A (en) | 2021-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112372372B (en) | Efficient milling cutter accumulated friction wear boundary identification and verification method | |
WO2021174525A1 (en) | Parts surface roughness and cutting tool wear prediction method based on multi-task learning | |
Duan et al. | Milling force model for aviation aluminum alloy: academic insight and perspective analysis | |
CN111366123B (en) | Part surface roughness and cutter wear prediction method based on multi-task learning | |
CN111002104B (en) | A method for detecting and solving the friction and wear boundary of the flank face of a high-feed milling cutter | |
CN102622489B (en) | Cutting force prediction method for five-axis side milling based on ACIS platform | |
CN114818201A (en) | A mechanism and data-driven method for predicting surface quality of engine cylinder head milling | |
Miko et al. | Analysis and verification of surface roughness constitution model after machining process | |
CN102736558A (en) | Numerical control machine tool thermal error real-time compensation modeling method based on time series algorithm | |
CN107627152B (en) | Numerical control machining chip control method based on BP neural network | |
CN112199632B (en) | A method for predicting the surface hardness of laser-textured aluminum alloys | |
CN114429064B (en) | Identification method for fractal characteristics of friction boundary of rear cutter surface of cutter tooth of high-energy-efficiency milling cutter | |
CN102866672A (en) | Online detecting method of numerical control machining middle state of plane structural member | |
CN104794337A (en) | Modeling method for judging cutting force of orthogonal turn-milling machining end face on basis of boundary conditions | |
CN114896753A (en) | Method for identifying instantaneous friction energy consumption distribution characteristics of rear cutter face of high-energy-efficiency milling cutter | |
CN118905468B (en) | Intelligent laser cutting control method and system | |
CN102566492A (en) | Method for forecasting maximum milling force for plunge milling of metal difficult-to-cut materials | |
CN110704969A (en) | Analytical Model Construction and Verification Method for Revealing the Wear Variation of High-feed Milling Cutter | |
CN106424969A (en) | Method for precisely predicting dynamic cutting force of plunge milling by taking cutter eccentricity into consideration | |
CN105242546B (en) | A kind of High Speed Milling Force modeling method based on material property | |
CN108038333A (en) | A kind of Forecasting Methodology of the Flexible Manufacture surface error of large slot milling cutter slabbing | |
CN118797844B (en) | A method for determining the sickle edge line of a sickle front cutter end mill | |
CN115204029A (en) | Method for evaluating and optimizing process reliability of spiral bevel gear milling machine tool | |
CN106944880B (en) | Method for detecting blade shape retentivity of cutting edge of turning coarse-pitch internal thread tool | |
CN116910928A (en) | Micro milling force measurement system and calculation method based on processing parameter inversion |
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 |