Cutting force modeling method for micro milling considering tool deformation

Abstract

The present invention provides a micro-milling cutting force modeling method that considers the influence of tool deformation. The method first calculates the real-time position of the cutting edge on the milling cutter that has undergone weak rigidity deformation, and thereby determines the actual tool tip cutting motion direction and the theoretical instantaneous undeformed chip thickness; then, according to the definition, the actual rake angle, shear angle and instantaneous undeformed chip thickness when the cutting direction is offset are calculated; then, considering the metal dead zone and material accumulation phenomenon, the cutting force at this time is calculated, and the milling process is modeled using an iterative algorithm, that is, the cutting force of the micro-milling process is obtained. The cutting force modeling method for the micro-milling process of a weak rigidity tool provided by the present invention considers the influence of the actual cutting direction on the actual value of the cutting parameters, and realizes the accurate prediction of the cutting force of the micro-milling process of a weak rigidity tool.

Description

Technical Field

[0001] The invention relates to a micro-milling cutting force modeling method, and in particular to a cutting force modeling method suitable for a micro-milling process using a weak-rigidity tool. Background Art

[0002] Micro-milling tools usually have weak rigidity due to their own slender structure, and are prone to severe deformation during the cutting process. When the tool deforms, the actual motion trajectory of the tool tip will change, causing the actual cutting direction to deviate from the theoretical cutting direction. However, by definition, the determination of parameters such as the rake angle, back angle, shear angle, and instantaneous undeformed chip thickness are all related to the cutting direction. When the cutting direction deviates, these parameters will change synchronously and affect the cutting force. Therefore, in the process of cutting force modeling for micro-milling processes using weak-rigidity tools, it is necessary to consider the impact of changes in cutting parameters caused by changes in the actual cutting direction.

[0003] The paper "TMMoges, KADesai, PVMRao, Modeling of cutting force, tool deflection, and surface error in micro-milling operation, The International Journal of Advanced Manufacturing Technology 98 (2018) 2865–2881." discloses a micro-milling cutting force modeling method that considers the tool's weak stiffness deformation. This method regards the micro-milling cutter as a cantilever beam, and after calculating the tool deformation during the machining process, the instantaneous undeformed thickness is determined based on the adjacent tooth trajectory, and then the milling process is modeled in combination with an iterative algorithm.

[0004] The typical characteristics of the above references are: in the weak rigidity milling process, only the position change of the milling cutter cutting edge is considered, while the change in cutting direction and its influence are ignored. In the micro-milling process, the deformation of the milling cutter is greater than its turning radius, and the change in cutting direction will be more drastic, which will have a serious impact on the cutting force and cutting process, and interfere with the prediction accuracy of the existing methods. Summary of the invention

[0005] Technical issues to be solved

[0006] In order to overcome the defect that the existing methods ignore the influence of actual cutting direction changes in the micro-milling modeling process of weak stiffness tools and cannot accurately model the cutting force, the present invention provides a micro-milling cutting force prediction method that takes into account the actual cutting direction changes and their influence.

[0007] Technical Solution

[0008] A micro-milling cutting force modeling method considering the influence of tool deformation is characterized by the following steps:

[0009] Step 1: The radius r of the tool used for microscopic observation is in millimeters; the radius r of the blunt circle of the cutting edge e , in millimeters; and the tool design rake angle α, in degrees;

[0010] Step 2: Consult the material manual and supplier data to obtain the shear yield strength τ of the workpiece material in MPa and the Young's modulus of the milling cutter material in MPa;

[0011] Step 3: Divide the tool into H equal units along the axial direction. For the mth unit, at time t=t i , the position change d caused by tool deformation m And the deformation direction α d ;

[0012] Step 4: Take the center of the machine tool spindle as the origin, define the feed direction as the positive direction of the x-axis, and the normal direction on the counterclockwise side of the feed direction as the positive direction of the y-axis, and establish a coordinate system; let at time t = t i , the tool rotation angle is φ(t i ), the unit is degree; the tool tip position coordinates at this moment (x(t i ),y(t i )) is calculated by the following formula:

[0013]

[0014] Where v f is the feed speed of the milling cutter;

[0015] Step 5: Let the previous moment t i-1 The tool tip position at this time is (x(t i-1 ),y(t i-1 )), calculate the time t by the following formula i Cutting direction offset φ d (t i ):

[0016]

[0017] Step 6: From the tool tip position (x(t i ),y(t i ), draw a ray along the tool rotation angle to the tool center, find the intersection of the ray and the previous tooth trajectory, calculate the distance between the tool tip position and the intersection, and use it as the current time t i The instantaneous undeformed chip thickness h 0 ;

[0018] Step 7: Calculate the nominal undeformed chip thickness h by the following formula n :

[0019] h n =h 0 +h p

[0020] Among them, h p is the material accumulation height generated at the previous moment;

[0021] Step 8: Actual Rake Angle α of the Milling Process a Calculated by the following formula:

[0022] α a =α+φ d (t i )

[0023] Step 9: Actual Shear Angle of Milling Process Calculated by the following formula:

[0024]

[0025] Where α s is the equivalent rake angle, in degrees; β s is the friction angle, in degrees;

[0026] Step 10: Cut the area thickness h s and the plowing area thickness h e Determined by the following formula:

[0027]

[0028] Step 11: Calculate the cutting force using the following formula:

[0029]

[0030] Among them, F t is the tangential cutting force, F r is the normal cutting force, in Newton; ι is the bevel cutting inclination angle, in degrees, which is consistent with the helix angle of the milling cutter during milling; B is the cutting width, in millimeters; K te , K re 、p te 、p re ,q te and q re is the plow shear force coefficient;

[0031] Step 12: For a milling cutter with K teeth, a total of J axial segments are involved in the milling process. The cutting force on the milling cutter is calculated by the following formula:

[0032]

[0033] where φ i,j (θ) is the instantaneous tooth position angle of the cutting edge on the jth axial segment of the kth tooth of the milling cutter, in degrees; F t,k,j (φ k,j (θ)) is the tangential cutting force on the cutting edge of the kth tooth of the milling cutter on the jth axial segment, F r,k,j (φ k,j (θ)) is the normal cutting force on the cutting edge of the jth axial segment of the kth tooth of the milling cutter;

[0034] Step 13: Repeat steps 3 to 12 for iteration until the cutter deformation d is m The change value is less than the given allowable error value;

[0035] Step 14: Enter the next moment t = t i+1 , let i=i+1, and repeat steps 3 to 14 until the cutting process is completed.

[0036] A computer system, characterized in that it comprises: one or more processors, and a computer-readable storage medium for storing one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement the above-mentioned method.

[0037] A computer-readable storage medium, characterized in that it stores computer-executable instructions, which are used to implement the above method when executed.

[0038] Beneficial Effects

[0039] The present invention provides a micro-milling cutting force modeling method that considers the influence of tool deformation. The method first uses a coordinate algorithm to model the instantaneous undeformed chip thickness at the cutting edge during the milling process under the condition of tool deformation, and calculates the deviation between the actual cutting direction and the theoretical cutting direction; then based on the actual cutting direction, according to the definition, the real values ​​of the instantaneous undeformed chip thickness, the milling cutter rake angle, and the shear angle during the milling process are recalculated; finally, the cutting force coefficient is updated using the real instantaneous undeformed chip thickness, the milling cutter rake angle, and the shear angle during the milling process, and the cutting force is calculated. Compared with the given literature, the present method takes the influence of the real cutting direction on the cutting force coefficient into the milling force model, can more accurately predict the influence of the weak stiffness deformation of the milling cutter on the cutting force, and improves the cutting force prediction accuracy under the condition of small tools and large deformation during micro-milling. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] The drawings are only for the purpose of illustrating particular embodiments and are not to be considered limiting of the present invention. Like reference symbols denote like components throughout the drawings.

[0041] Figure 1 Comparison chart of predicted cutting force results and measured cutting force data;

[0042] Figure 2 Comparison of the predicted milling cutter tip trajectory results and the tool mark image of the experimentally machined surface. DETAILED DESCRIPTION

[0043] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0044] The present invention provides a micro-milling cutting force prediction method that considers the actual cutting direction change and its influence. The method first calculates the real-time position of the cutting edge on the milling cutter that has undergone weak rigidity deformation, and determines the actual tool tip cutting motion direction and the theoretical instantaneous undeformed chip thickness based on this; then, according to the definition, the actual rake angle, shear angle and instantaneous undeformed chip thickness when the cutting direction is offset are calculated; then, considering the metal dead zone and material accumulation phenomenon, the cutting force at this time is calculated, and the milling process is modeled using an iterative algorithm, that is, the cutting force of the micro-milling process is obtained.

[0045] The cutting force modeling method for the micro-milling process of a weak-rigidity tool provided by the present invention takes into account the influence of the actual cutting direction on the actual value of the cutting parameter, thereby realizing accurate prediction of the cutting force in the micro-milling process of a weak-rigidity tool.

[0046] A cutting force modeling method for micro-milling process with weak-rigidity tool considering the real cutting parameter variation includes the following steps:

[0047] Step 1: Microscopic observation: The radius r of the tool used, in millimeters, and the radius r of the blunt circle of the cutting edge e , in millimeters, and the tool design rake angle α, in degrees.

[0048] Step 2: Consult the material manual and supplier data to obtain the shear yield strength τ of the workpiece material in MPa and the Young's modulus of the milling cutter material in MPa.

[0049] Step 3: Divide the tool into H equal units along the axial direction. For the mth unit, at time t=t i , the position change d caused by tool deformation mAnd the deformation direction α d Determined by referring to the method disclosed in the literature “TMMoges, KADesai, PVMRao, Modeling of cutting force, tool deflection, and surface error in micro-milling operation, The International Journal of Advanced Manufacturing Technology 98 (2018) 2865–2881.”

[0050] Step 4: Take the center of the machine tool spindle as the origin, define the feed direction as the positive direction of the x-axis, and the normal direction on the counterclockwise side of the feed direction as the positive direction of the y-axis, and establish a coordinate system. Let at time t = t i , the tool rotation angle is φ(t i ), the unit is degree. The tool tip position coordinates at this moment (x(t i ),y(t i )) is calculated by the following formula:

[0051]

[0052] Where v f is the feed speed of the milling cutter.

[0053] Step 5: Let the previous moment t i-1 The tool tip position at this time is (x(t i-1 ),y(t i-1 )), calculate the time t by the following formula i Cutting direction offset φ d (t i ):

[0054]

[0055] Step 6: From the tool tip position (x(t i ),y(t i ), draw a ray along the tool rotation angle to the tool center, find the intersection of the ray and the previous tooth trajectory, calculate the distance between the tool tip position and the intersection, and use it as the current time t i The instantaneous undeformed chip thickness h 0 .

[0056] Step 7: Let the material accumulation height generated at the previous moment be h p, in millimeters, is determined according to the literature "D.-Y.Wen, M.Wan, Y.-Y.Ren, W.-H.Zhang, Y.Yang, Material piling up and spreading effects in the cutting processes with small feed rates, Mechanical Systems and Signal Processing 171 (2022) 108839.". The nominal undeformed chip thickness h is calculated by the following formula n :

[0057] h n =h 0 +h p

[0058] Step 8. Actual rake angle α during milling a Calculated by the following formula:

[0059] α a =α+φ d (t i )

[0060] Step 9. Actual shear angle during milling Calculated by the following formula:

[0061]

[0062] Where α s is the equivalent rake angle, in degrees, and the determination method refers to the literature "G. Bissacco, HN Hansen, J. Slunsky, Modelling the cutting edge radius size effect for force prediction in micro milling, CIRP Annals-Manufacturing Technology 57 (1) (2008) 113–116.". β sis the friction angle, in degrees, determined by referring to the method disclosed in the document “M. Wan, D.-Y. Wen, Y.-C. Ma, W.-H. Zhang, On materials separation and cutting force prediction in micro milling through involving the effect of dead metal zone, International Journal of Machine Tools and Manufacture 146 (2019) 103452.”

[0063] Step 10: Shearing area thickness h s and the plowing area thickness h e Determined by the following formula

[0064]

[0065] Step 11. Calculate the cutting force using the following formula:

[0066]

[0067] Among them, F t is the tangential cutting force, F r is the normal cutting force, in Newton. ι is the bevel cutting inclination angle, in degrees, which is consistent with the helix angle of the milling cutter during milling. B is the cutting width, in millimeters. K te , K re 、p te 、p re ,q te and q re is the plowing force coefficient, which is determined by referring to the method disclosed in the document “M. Wan, D.-Y. Wen, Y.-C. Ma, W.-H. Zhang, On material separation and cutting force prediction in micro milling through involving the effect of dead metal zone, International Journal of Machine Tools and Manufacture 146 (2019) 103452.”

[0068] Step 12: For a milling cutter with K teeth, a total of J axial segments are involved in the milling process. The cutting force on the milling cutter is calculated by the following formula:

[0069]

[0070] where φ i,j (θ) is the instantaneous tooth position angle of the cutting edge on the jth axial segment of the kth tooth of the milling cutter, in degrees. t,k,j (φ k,j (θ)) is the tangential cutting force on the cutting edge of the kth tooth of the milling cutter on the jth axial segment, F r,k,j (φ k,j (θ)) is the normal cutting force on the cutting edge of the j-th axial segment of the k-th tooth of the milling cutter.

[0071] Step 13: Repeat steps 3 to 12 for iteration until the cutter deformation d is equal to the value of the cutter deformation d between two adjacent iterations. m The change value is less than the given allowable error value.

[0072] Step 14: Enter the next moment t=t i+1 , let i=i+1, and repeat steps 3 to 14 until the cutting process is completed.

[0073] In order to enable those skilled in the art to better understand the present invention, the present invention is described in detail below in conjunction with specific embodiments.

[0074] Embodiment 1:

[0075] In the experiment, a flat-bottomed carbide end mill with the number of teeth K = 2, radius r = 0.5 mm, helix angle β = 30 degrees, and a designed rake angle α = 5 degrees was selected to carry out milling tests on aluminum alloy thin-walled parts using the down-milling cutting method on a three-coordinate vertical machining center, with a radial cutting depth of 0.1 mm and a feed per tooth of 0.01 mm / tooth.

[0076] Step 1: Obtain the radius r of the blunt circle of the cutting edge through microscope observation e =0.012 mm.

[0077] Step 2: By consulting the material manual, determine that the shear yield strength of aluminum alloy Al 7050-T7451 is τ = 500 MPa. By consulting the supplier data, confirm that the Young's modulus of the milling cutter tool material is 211000 MPa.

[0078] Step 3: Divide the tool into 200 equal units along the axial direction. For the mth unit, at time t=t i , the position change d caused by tool deformation m And the deformation direction α dDetermined by referring to the method disclosed in the literature “TMMoges, KADesai, PVMRao, Modeling of cutting force, tool deflection, and surface error in micro-milling operation, The International Journal of Advanced Manufacturing Technology 98 (2018) 2865–2881.”

[0079] Step 4: Take the center of the machine tool spindle as the origin, define the feed direction as the positive direction of the x-axis, and the normal direction on the counterclockwise side of the feed direction as the positive direction of the y-axis, and establish a coordinate system. Let at time t = t i , the tool rotation angle is φ(t i ), the unit is degree. The tool tip position coordinates at this moment (x(t i ),y(t i )) is calculated by the following formula:

[0080]

[0081] Where v f is the feed speed of the milling cutter.

[0082] Step 5: Let the previous moment t i-1 The tool tip position at this time is (x(t i-1 ),y(t i-1 )), the cutting direction offset φ at time t is calculated by the following formula d (t i ):

[0083]

[0084] Step 6: From the tool tip position (x(t i ),y(t i ), draw a ray toward the tool center along the tool rotation angle, find the intersection of the ray and the previous tooth trajectory, calculate the distance between the tool tip position and the intersection, and use it as the instantaneous undeformed chip thickness h at the current time t 0 .

[0085] Step 7: Let the material accumulation height generated at the previous moment be h p, in millimeters, is determined according to the literature "D.-Y.Wen, M.Wan, Y.-Y.Ren, W.-H.Zhang, Y.Yang, Material piling up and spreading effects in the cutting processes with small feed rates, Mechanical Systems and Signal Processing 171 (2022) 108839.". The nominal undeformed chip thickness h is calculated by the following formula n :

[0086] h n =h 0 +h p

[0087] Step 8. Actual rake angle α during milling a Calculated by the following formula:

[0088] α a =α+φ d (t i )

[0089] Step 9. Actual shear angle during milling Calculated by the following formula:

[0090]

[0091] Where α s is the equivalent rake angle, in degrees, and the determination method refers to the literature "G. Bissacco, HN Hansen, J. Slunsky, Modelling the cutting edge radius size effect for force prediction in micro milling, CIRP Annals-Manufacturing Technology 57 (1) (2008) 113–116.". β sis the friction angle, in degrees, determined by referring to the method disclosed in the document “M. Wan, D.-Y. Wen, Y.-C. Ma, W.-H. Zhang, On materials separation and cutting force prediction in micro milling through involving the effect of dead metal zone, International Journal of Machine Tools and Manufacture 146 (2019) 103452.”

[0092] Step 10: Shearing area thickness h s and the plowing area thickness h e Determined by the following formula

[0093]

[0094] Step 11. Calculate the cutting force using the following formula:

[0095]

[0096] Among them, F t is the tangential cutting force, F r is the normal cutting force, in Newton. ι is the bevel cutting inclination angle, in degrees, which is equal to the milling cutter helix angle, which is 30 degrees. B is the cutting width, in millimeters. K te , K re 、p te 、p re ,q te and q re is the plowing force coefficient, which is determined by referring to the method disclosed in the literature “M. Wan, D.-Y. Wen, Y.-C. Ma, W.-H. Zhang, On material separation and cutting force prediction in micro milling through involving the effect of dead metal zone, International Journal of Machine Tools and Manufacture 146 (2019) 103452.”: K te =1182, K re =1584, p te =0.8532, p re =0.6071, q te =0.01174, qre =0.2832.

[0097] Step 12: For the current example, the milling cutter has 2 teeth and 2 axial segments involved in the cutting process. The cutting force on the milling cutter is calculated using the following formula:

[0098]

[0099] where φ k,j (θ) is the instantaneous tooth position angle of the cutting edge on the jth axial segment of the kth tooth of the milling cutter, in degrees. t,k,j (φ k,j (θ)) is the tangential cutting force on the cutting edge of the kth tooth of the milling cutter on the jth axial segment, F r,k,j (φ k,j (θ)) is the normal cutting force on the cutting edge of the j-th axial segment of the k-th tooth of the milling cutter.

[0100] Step 13: Repeat steps 3 to 12 for iteration until the cutter deformation d is equal to the value of the cutter deformation d between two adjacent iterations. m The change value is less than the given allowable error of 0.00005mm.

[0101] Step 14: Enter the next moment t=t i+1 , let i=i+1, and repeat steps 3 to 14 until the cutting process is completed.

[0102] After the calculation is completed, draw the x-direction cutting force F of the milling cutter x (θ) and y-direction cutting force F y (θ) changes with the milling cutter rotation angle θ, and the attached Figure 1 . Plot the coordinates x(t) and y(t) of the milling cutter tip position during the milling process and obtain the attached Figure 2 .

[0103] Through the above steps, the theoretical prediction results of cutting force for cutting aluminum alloy Al7050-T7451 in a down milling mode can be obtained using a two-tooth micro-milling cutter with a diameter of 1mm, an overhang length of 10mm, a design rake angle of 5 degrees, a helix angle of 30 degrees, and a cutting edge tip bluntness of 0.012mm at a spindle speed of 5000 rpm, an axial cutting depth of 0.5mm, a radial cutting depth of 0.1mm, and a feed per tooth of 0.01mm. Figure 1 , and tool tip trajectory prediction results, refer to the attached Figure 2 . By attaching Figure 1 It can be seen that the cutting force results predicted by the method of the present invention are consistent with the experimental results. Figure 2It can be seen that the trajectory of the milling cutter tip predicted by the method of the present invention is consistent with the morphology of the tool mark on the machined surface, which means that the trajectory of the milling cutter tip predicted by the method of the present invention is consistent with the actual trajectory of the tool tip.

[0104] Attached Figure 1 and attached Figure 2 The effectiveness of the proposed cutting force modeling method for micro-milling process with weak-rigidity tools considering the real cutting parameters variation is demonstrated.

[0105] Figure 1 The predicted cutting force results and the measured cutting force data are compared when a two-tooth micro-milling cutter with a diameter of 1 mm, an overhang length of 10 mm, a designed rake angle of 5 degrees, a helix angle of 30 degrees, and a cutting edge tip bluntness of 0.012 mm is used to cut aluminum alloy Al 7050-T7451 in a down milling manner at a spindle speed of 5000 rpm, an axial cutting depth of 0.5 mm, a radial cutting depth of 0.1 mm, and a feed per tooth of 0.01 mm.

[0106] Figure 2 The predicted milling cutter tip trajectory results of the two-tooth micro-milling cutter with a diameter of 1 mm, an overhang length of 10 mm, a designed rake angle of 5 degrees, a helix angle of 30 degrees, and a cutting edge tip bluntness of 0.012 mm in the embodiment of the present invention are compared with the experimental machining surface tool mark image when cutting aluminum alloy Al 7050-T7451 in a down milling manner at a spindle speed of 5000 rpm, an axial cutting depth of 0.5 mm, a radial cutting depth of 0.1 mm, and a feed per tooth of 0.01 mm.

[0107] The above description is only a specific implementation mode of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of various equivalent modifications or substitutions within the technical scope disclosed in the present invention, and these modifications or substitutions should be included in the protection scope of the present invention.

Claims

1. A micro-milling cutting force modeling method considering the influence of tool deformation, characterized in that Here are the steps: Step 1: The radius r of the tool used for microscopic observation is in millimeters; the radius r of the blunt circle of the cutting edge e , in millimeters; and the tool design rake angle α, in degrees; Step 2: Consult the material manual and supplier data to obtain the shear yield strength τ of the workpiece material in MPa; and Young’s modulus of the milling cutter material in MPa; Step 3: Divide the tool into H equal units along the axial direction. For the mth unit, at time t=t i , the position change d caused by tool deformation m And the deformation direction α d ; Step 4: Take the center of the machine tool spindle as the origin, define the feed direction as the positive direction of the x-axis, and the normal direction on the counterclockwise side of the feed direction as the positive direction of the y-axis, and establish a coordinate system; let at time t = t i , the tool rotation angle is φ(t i ), the unit is degree; At this moment, the tool tip position coordinates (x(t i ),y(t i )) is calculated by the following formula: Where v f is the feed speed of the milling cutter; Step 5: Let the previous moment t i-1 The tool tip position at this time is (x(t i-1 ),y(t i-1 )), the time t is calculated by the following formula i Cutting direction offset φ d (t i ): Step 6: From the tool tip position (x(t i ),y(t i ), draw a ray along the tool rotation angle to the tool center, find the intersection of the ray and the previous tooth trajectory, calculate the distance between the tool tip position and the intersection, and use it as the current time t i The instantaneous undeformed chip thickness h0; Step 7: Calculate the nominal undeformed chip thickness h by the following formula n : h n =h0+h p Among them, h p is the material accumulation height generated at the previous moment; Step 8: Actual Rake Angle α of the Milling Process a Calculated by the following formula: a a =a+f d (t i ) Step 9: Actual Shear Angle of Milling Process Calculated by the following formula: Where α s is the equivalent rake angle, in degrees; β s is the friction angle, in degrees; Step 10: Cut the area thickness h s and the plowing area thickness h e Determined by the following formula: Step 11: Calculate the cutting force using the following formula: Among them, F t is the tangential cutting force, F r is the normal cutting force, in Newton; ι is the bevel cutting inclination angle, in degrees, which is consistent with the helix angle of the milling cutter during milling; B is the cutting width, in millimeters; K te , K re 、p te 、p re ,q te and q re is the plow shear force coefficient; Step 12: For a milling cutter with K teeth, a total of J axial segments are involved in the milling process. The cutting force on the milling cutter is calculated by the following formula: where φ i,j (θ) is the instantaneous tooth position angle of the cutting edge on the jth axial segment of the kth tooth of the milling cutter, in degrees; F t,k,j (φ k,j (θ)) is the tangential cutting force on the cutting edge of the kth tooth of the milling cutter on the jth axial segment, F r,k,j (φ k,j (θ)) is the normal cutting force on the cutting edge of the jth axial segment of the kth tooth of the milling cutter; Step 13: Repeat steps 3 to 12 for iteration until the cutter deformation d is m The change value is less than the given allowable error value; Step 14: Enter the next time t = t i+1 , let i=i+1, and repeat steps 3 to 14 until the cutting process is completed.

2. A computer system, characterized in that include: One or more processors, and a computer-readable storage medium for storing one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors are enabled to implement the method of claim 1.

3. A computer-readable storage medium, characterized in that Computer executable instructions are stored, and when the instructions are executed, they are used to implement the method of claim 1.

Images