High-precision laser de-weight method for de-weight dynamic balance of rotor

Abstract

A high-precision laser de-weight method for de-weight dynamic balance of rotors comprises the following steps: firstly, establishing a standard de-weight blind hole parameter library for different rotors, and removing any unbalance vector in a specified range of the corresponding rotors through a plurality of groups of standard de-weight blind hole combinations; determining various laser processing parameters for processing corresponding de-weight blind holes according to laser simulation and experiments; selecting corresponding de-weighting blind hole parameters by utilizing a standard de-weighting parameter library according to the numerical value of the rotor; performing trial de-weighting on the de-weighting blind holes by using the selected laser parameters on a test piece made of the same material as the to-be-balanced rotor, and correcting a real de-weighting value; according to a real de-weight value obtained after the de-weight blind hole is corrected, the actual de-weight position is corrected through an integration method, and de-weight is carried out on the batch rotors; after deduplication, detecting a deduplication effect, and if the deduplication effect is not qualified, measuring an unbalanced vector and carrying out deduplication again; according to the method, the influence of nonlinear factors on the de-weight quality in the laser de-weight process is reduced, and the precision and efficiency of de-weight dynamic balance of batch micro rotors are improved.

Description

Technical Field

[0001] The invention relates to the technical field of laser deweighting, and in particular to a high-precision laser deweighting method for rotor deweighting dynamic balancing. Background Art

[0002] As a high-speed rotating precision component, the balance of the micro gyro rotor is crucial to ensure the normal operation of the equipment, improve work efficiency and extend service life. The traditional method uses mechanical tools to remove weight: after calculating the rotor imbalance, a hole of a specified depth is drilled at a certain location to remove the unbalanced mass. However, this method has problems such as low deweighting accuracy and high technical requirements for operators.

[0003] Laser deweighting (Wang Qingkai, Zhang Yu, Chen Dongsuo, et al. Research on precision control of ultra-high-speed lightweight rotor dynamic balancing technology [J]. Daily Electrical Appliances, 2022 (2): 67-69.) Compared with traditional contact methods, laser deweighting has the advantages of non-contact deweighting and high deweighting accuracy, and is more suitable for use in ultra-precision deweighting scenarios of micro rotors; however, due to the different properties of the rotor surface oxide layer and the metal in the laser processing deweighting process, the complex parameters lead to nonlinearity between the laser processing time and the removal quality, resulting in low laser deweighting accuracy. The field of laser deweighting requires a method for precise removal of tiny masses. Summary of the invention

[0004] In order to overcome the above-mentioned shortcomings of the prior art, the purpose of the present invention is to provide a high-precision laser deweighting method for rotor deweighting and dynamic balancing, so as to achieve the processing requirements of efficient, high-precision and high-quality rotor dynamic balancing.

[0005] In order to achieve the above object, the present invention is implemented by the following technical solutions:

[0006] A high-precision laser deweighting method for rotor deweighting dynamic balancing comprises the following steps:

[0007] Step 1: Establish a standard de-duplication blind hole parameter library for different rotor sizes, and remove any imbalance vector within the specified range of the corresponding rotor through multiple groups of standard de-duplication blind hole combinations. The de-duplication selection principle of the standard de-duplication blind hole parameter library is shown in formula (1):

[0008] U i =M i ·(X i -s i )~M i ·(X i +s i ) (1)

[0009] U i M is the range of imbalance vectors that can be corrected by single blind hole deweighting.i is the deduplication quality of laser single blind hole, which are M1, M2, M3...M n , are all selected fixed values, X i The weight M of the single punching i The distance between the center of the corresponding deweighting range and the center of the deweighting workpiece rotor, s i For X i Corresponding to the range, duplicate blind holes are removed by selecting the appropriate position Xs within the range i To remove the appropriate deduplication vector;

[0010] For any deweighting vector U within the specified range, by selecting different blind hole qualities M i Composition of different U i , making U=U1+U2+...+U n , thereby achieving the removal of the specified mass;

[0011] According to X i The deduplication range determines the diameter R of the corresponding deduplication blind hole i , ensure that R i +R i-1 <(X i -s i )-(X i-1 +s i-1 ), so that different deduplication blind holes will not affect each other. The relationship between the deduplication weight and the radius is: Where ρ1 is the density of the rotor substrate, ρ2 is the density of the rotor surface layer; h1 is the de-weighting depth of the blind hole on the rotor substrate, and h2 is the de-weighting depth of the blind hole on the rotor surface layer;

[0012] Step 2: Standard deduplication blind hole processing: first, evaporate and remove the oxide layer on the rotor surface, and then evaporate and remove the rotor substrate layer. The parameter table establishes the functional relationship between laser parameters and deduplication amount based on laser simulation and experiment. Based on this functional relationship, various laser processing parameters corresponding to the deduplication blind hole oxide layer and substrate layer are established, including laser power, spot radius, defocus, scanning speed, spot overlap rate, and number of scans.

[0013] Step 3: Select corresponding deduplication blind hole parameters using the standard deduplication blind hole parameter library according to the shape, mass, density and other values ​​of the rotor;

[0014] Step 4: Use the selected laser parameters to test several deweighting blind holes on the test piece with the same material as the rotor to be balanced, and correct the actual deweighting value. Before the trial deweighting, test several sets of selected values ​​on the same material sample, measure one for more than 5 times and take the average value to calibrate the parameters;

[0015] Step 5: According to the real de-weighting value obtained after the de-weighting blind hole correction, the actual de-weighting position is corrected by the integral method to de-weight the rotor;

[0016] Step 6: After deduplication, check the deduplication effect. If it is unqualified, the unbalanced vector is measured and returned to step 5 to perform deduplication again until the deduplication effect is qualified.

[0017] The standard deduplication blind hole parameter library described in step 1 must contain deduplication blind hole parameters and blind hole position ranges of various diameters and depths corresponding to various rotor materials and sizes, with a minimum deduplication weight of 0.1 mg and a minimum imbalance balance accuracy of 0.5 mg.mm.

[0018] Step 2 The simulation process of laser processing parameters includes the following steps:

[0019] 2.1) First, select an 8-node hexahedral unit, establish a 3D model of the deweighted material, and perform laser processing deweighting simulation to simulate the temperature field and thermal stress of the rotor. When the unit temperature reaches the vaporization temperature or the thermal stress between units of different materials is greater than the limit adhesion, the unit is considered to be removed;

[0020] 2.2) According to the physical constants and basic equations, the initial temperature parameters and boundary conditions of the geometric body are defined for simulation. The output of the previous time step is used as the input of the next time step. The pulsed laser is input in the form of heat flux. The heat flux distribution function is shown in formula (2):

[0021]

[0022] Where A is the absorption rate, P is the average laser power, r is the spot radius, t is the thermal conduction time, v is the scanning speed, τ is the time relative to the full width of the pulse to reach half the maximum value, τ m is the time to reach the maximum power;

[0023] 2.3) The laser heat source input is simplified to a point heat source on the spatial section and a pulse input on the time section; the load matrix list is used to provide the corresponding heat flux loads at different positions and times. The heat is radiated to the surface of the oxide layer and diffused to the bottom of the oxide layer and the substrate through heat conduction, following the heat balance equation in Fourier's law, as shown in formula (3):

[0024]

[0025] Where, the temperature field distribution is T(x,y,z,t), the heat flux distribution is q(x,y,z,t), the energy input is g(x,y,z,t), κ is the thermal conductivity, and α=κ / ρc is the thermal diffusion coefficient;

[0026] 2.4) During the laser cleaning process, there are differences in the thermal expansion coefficient γ and Young's modulus Y between the oxide film layer inside the material and the base material, which leads to the generation of thermal stress, as shown in formula (4):

[0027]

[0028] Where u is the displacement, γ is the thermal expansion coefficient, and Y is Young's modulus;

[0029] The adhesion force f between the substrate layer and the oxide layer is Z is the distance between the two layers, A 12 is the Hamaker coefficient between the oxide film and the substrate. When the thermal stress is greater than the adhesion force, the oxide film is removed;

[0030] 2.5) Based on the above model, simulation analysis and experiments are combined to determine the laser processing parameters.

[0031] The maximum deweighting mass selected in step 3 shall not exceed 1% of the rotor mass, and the position ranges of the selected standard deweighting blind holes shall not overlap.

[0032] The parameters for processing the corresponding blind holes in step 4 are based on simulation formulas. If the difference between the actual weight removal value and the simulation value is too large to meet the actual weight removal range requirements, other blind hole values ​​in the parameter table are selected for assistance.

[0033] The integration method in step five includes using a numerical integration algorithm to calculate a correction value of the actual deduplication position.

[0034] Compared with the prior art, the present invention has the following beneficial effects:

[0035] (1) The present invention divides the laser dynamic balancing deweighting process into a surface oxide layer removal mode and a substrate removal mode, and achieves deweighting of different layers of the workpiece by changing the laser processing parameters, so its deweighting effect is more accurate.

[0036] (2) The present invention determines the laser deduplication parameters of blind holes of different sizes through finite element simulation, and can more quickly establish corresponding parameters for different materials.

[0037] (3) The present invention establishes a standard deduplication blind hole parameter library and reasonably selects parameter combinations, thereby reducing the influence of various nonlinear factors on deduplication accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] Figure 1 A flowchart of an embodiment.

[0039] Figure 2 It is a schematic diagram of a model for simulating blind hole combination deduplication in an embodiment.

[0040] Figure 3 Schematic diagram of blind hole removal in an embodiment.

[0041] Figure 4 It is a schematic diagram of the equipment used for actual deduplication in the embodiment. DETAILED DESCRIPTION

[0042] The present invention will be further described below in conjunction with the accompanying drawings and embodiments.

[0043] During the laser processing deweighting process, there are reasons such as the different properties of the rotor surface oxide layer and the metal, and the complex parameters leading to nonlinearity between the laser processing time and the removal quality, which result in low laser deweighting accuracy. Therefore, it is necessary to use computer simulation methods to establish a standard deweighting blind hole parameter library and corresponding laser processing parameters. By selecting blind holes and laser processing parameters corresponding to the material size of the rotor, accurate dynamic balancing and deweighting of the micro rotor can be achieved.

[0044] like Figure 1 As shown, a high-precision laser deweighting method for rotor deweighting dynamic balancing is used. A standard deweighting blind hole parameter library for different rotor sizes and materials is established. Various laser processing parameters for processing corresponding deweighting blind holes are determined based on laser simulation and experiments. Corresponding deweighting blind hole parameters are selected according to the rotor. The selected laser parameters are used on a test piece with the same material as the rotor to be balanced to test deweight several deweighting blind holes. The actual deweighting position is corrected by the integral method based on the actual deweighting value obtained after the correction of the deweighting blind hole, and the batch rotors are deweighted to improve the deweighting accuracy. The specific steps are as follows:

[0045] Step 1: Establish a standard deduplication blind hole parameter library for different rotor sizes, and remove any imbalance vector within the specified range of the corresponding rotor through multiple sets of standard deduplication blind hole combinations, such as Figure 2 As shown;

[0046] The deduplication selection principle of the standard deduplication blind hole parameter library is shown in formula (1):

[0047] U i =M i ·(X i -s i )~M i ·(X i +s i ) (1)

[0048] U i M is the range of imbalance vectors that can be corrected by single blind hole deweighting. i is the deduplication quality of laser single blind hole, which are M1, M2, M3...M n , are all selected fixed values, X i The weight M of the single punching iThe distance between the center of the corresponding deweighting range and the center of the deweighting workpiece rotor, s i For X i Corresponding to the range, duplicate blind holes are removed by selecting the appropriate position Xs within the range i To remove the appropriate deduplication vector;

[0049] For any deweighting vector U within the specified range, by selecting different blind hole qualities M i Composition of different U i , making U=U1+U2+...+U n , thereby achieving the removal of the specified mass;

[0050] According to X i The deduplication range determines the diameter R of the corresponding deduplication blind hole i , ensure that R i +R i-1 <(X i -s i )-(X i-1 +s i-1 ), so that different deduplication blind holes will not affect each other. The relationship between the deduplication weight and the radius is: Where ρ1 is the density of the rotor substrate, ρ2 is the density of the rotor surface layer; h1 is the de-weighting depth of the blind hole on the rotor substrate, and h2 is the de-weighting depth of the blind hole on the rotor surface layer;

[0051] The standard de-weighting blind hole parameter library must contain de-weighting blind hole parameters and blind hole position intervals of various diameters and depths corresponding to various rotor materials and sizes. The minimum de-weighting weight is 0.1 mg, and the minimum unbalanced balance accuracy is 0.5 mg.mm.

[0052] The standard blind hole deduplication parameter library of this embodiment is shown in Table 1;

[0053] Table 1

[0054] <![CDATA[Deduplication quality (M i )]]> <![CDATA[M1]]> <![CDATA[M2]]> <![CDATA[M3]]> …… <![CDATA[M n ]]> <![CDATA[Blind hole radius (R i )]]> <![CDATA[R1]]> <![CDATA[R2]]> <![CDATA[R3]]> …… <![CDATA[R n ]]> <![CDATA[Blind hole depth (H i )]]> <![CDATA[H1]]> <![CDATA[H2]]> <![CDATA[H3]]> …… <![CDATA[H n ]]> <![CDATA[Deduplication range (X i )]]> <![CDATA[X1(±s1)]]> <![CDATA[X2(±s2)]]> <![CDATA[X3(±s3)]]> …… <![CDATA[X n (±s n )]]> <![CDATA[Deduplicated vector range (U i )]]> <![CDATA[U1]]> <![CDATA[U2]]> <![CDATA[U3]]> …… <![CDATA[U n ]]>

[0055] Step 2: Standard deduplication blind hole processing, refer to Figure 3 , firstly, the oxide layer on the surface of the rotor is evaporated and removed, and then the rotor substrate layer is evaporated and removed, where h1 is the de-weighting depth of the blind hole on the rotor substrate, h2 is the de-weighting depth of the blind hole on the rotor surface layer, and D is the de-weighting radius of the blind hole; the parameter table establishes the functional relationship between the laser parameters and the de-weighting amount according to laser simulation and experiments, and based on this functional relationship, establishes various laser processing parameters for processing the corresponding de-weighting blind hole oxide layer and substrate layer, including laser power, spot radius, defocus, scanning speed, spot overlap rate, and scanning times;

[0056] The simulation process of laser processing parameters includes the following steps:

[0057] 2.1) First, select an 8-node hexahedral unit, establish a 3D model of the deweighted material, and perform laser processing deweighting simulation to simulate the temperature field and thermal stress of the rotor. When the unit temperature reaches the vaporization temperature or the thermal stress between units of different materials is greater than the limit adhesion, the unit is considered to be removed;

[0058] 2.2) According to the physical constants and basic equations, the initial temperature parameters and boundary conditions of the geometric body are defined for simulation. The output of the previous time step is used as the input of the next time step. The pulsed laser is input in the form of heat flux. The heat flux distribution function is shown in formula (2):

[0059]

[0060] Where A is the absorption rate, P is the average laser power, r is the spot radius, t is the thermal conduction time, v is the scanning speed, τ is the time relative to the full width of the pulse to reach half the maximum value, τ m is the time to reach the maximum power;

[0061] 2.3) The laser heat source input is simplified to a point heat source on the spatial section and a pulse input on the time section; the load matrix list is used to provide the corresponding heat flux loads at different positions and times. The heat is radiated to the surface of the oxide layer and diffused to the bottom of the oxide layer and the substrate through heat conduction, following the heat balance equation in Fourier's law, as shown in formula (3):

[0062]

[0063] Where, the temperature field distribution is T(x,y,z,t), the heat flux distribution is q(x,y,z,t), the energy input is g(x,y,z,t), κ is the thermal conductivity, and α=κ / ρc is the thermal diffusion coefficient;

[0064] 2.4) During the laser cleaning process, there are differences in the thermal expansion coefficient γ and Young's modulus Y between the oxide film layer inside the material and the base material, which leads to the generation of thermal stress, as shown in formula (4):

[0065]

[0066] Where u is the displacement, γ is the thermal expansion coefficient, and Y is Young's modulus;

[0067] The adhesion force f between the substrate layer and the oxide layer is Z is the distance between the two layers, A 12 is the Hamaker coefficient between the oxide film and the substrate. When the thermal stress is greater than the adhesion force, the oxide film is removed;

[0068] 2.5) Based on the above model, simulation analysis and experiments are conducted to determine the laser processing parameters;

[0069] In this embodiment, the laser power is 50W, the spot radius is 40.6μm, the defocus is -0.6mm, the scanning speed is 1000mm / s, the spot overlap rate is 49.71%, and the number of scans is 10 times;

[0070] Step 3: Select corresponding de-duplication blind hole parameters using the standard de-duplication blind hole parameter library according to the shape, mass, density and other values ​​of the rotor, and the position ranges of the selected standard de-duplication blind holes cannot overlap; the selected maximum de-duplication mass shall not exceed 1% of the rotor mass;

[0071] The maximum radius of the batch rotor to be deweighted in this embodiment is 22 mm, the mass is 25.1 g, the material is 6061 aluminum alloy, the surface anodized layer is 10 μm, and the blind holes are selected as 0.1 mg, 1 mg, 10 mg, 50 mg, 100 mg, and 250 mg respectively;

[0072] Step 4: Use the selected laser parameters to test several deweighting blind holes on the test piece with the same material as the rotor to be balanced, and correct the actual deweighting value. Before the test, test the selected values ​​of several blind holes on the same material sample; m1, m2, m3, m4, m5; measure one 5 times and take the average value to calibrate the parameters;

[0073] This embodiment takes a 100mg standard blind hole as an example, and uses the selected laser parameters for processing. The average value of the processing is 91.3mg after 5 times of processing. The calibration of other blind holes is carried out according to this method. The actual weight removal is 0.1mg, 0.9mg, 9.7mg, 47.1mg, 91.3mg, and 221.9mg. The actual weight removal selection and calculation are based on this result.

[0074] Step 5: According to the real de-weighting value obtained after the de-weighting blind hole correction, the actual de-weighting position is corrected by the integral method to de-weight the rotor;

[0075] In this embodiment, the unbalance vector of the rotor is measured to be 1471.3 mg.mm. The blind holes and deweighting positions are selected according to the real deweighting values ​​obtained after the deweighting blind holes are corrected, which are: 0.1 mg*5.00 mm, 0.9 mg*7.52 mm, 9.7 mg*10,22 mm, 91.3 mg*14.95 mm, respectively. The actual deweighting position is corrected by the integral method to deweight the rotor.

[0076] Step 6: After deduplication, check the deduplication effect. If it is unqualified, the unbalanced vector is measured and returned to step 5 to perform deduplication again until the deduplication effect is qualified.

[0077] Reference Figure 4In this embodiment, the device for actual deweighting controls the processing laser height through the z-direction lead screw motor, and the z-direction lead screw motor adjusts the horizontal position through the x-direction lead screw motor. The processing laser is irradiated on the rotor to be balanced, and the rotor to be balanced is installed on the dynamic balancing spindle.

[0078] The beneficial effects of this embodiment are as follows: the laser dynamic balancing deweighting process of this embodiment is divided into a surface oxide layer removal mode and a substrate removal mode, and deweighting of different layers of the workpiece is achieved by changing the laser processing parameters, so its deweighting effect is more accurate. This embodiment determines the laser deweighting parameters of blind holes of different sizes through finite element simulation, and can more quickly establish corresponding parameters for different materials. This embodiment establishes a standard deweighting blind hole parameter library and reasonably selects parameter combinations, thereby reducing the influence of various nonlinear factors on deweighting accuracy.

Claims

1. A high-precision laser deweighting method for rotor deweighting dynamic balancing, characterized in that: The steps include: Step 1: Establish a standard de-duplication blind hole parameter library for different rotor sizes, and remove any imbalance vector within the specified range of the corresponding rotor through multiple groups of standard de-duplication blind hole combinations. The de-duplication selection principle of the standard de-duplication blind hole parameter library is shown in formula (1): U i =M i ·(X i -s i )~M i ·(X i +s i )(1) U i M is the range of imbalance vectors that can be corrected by single blind hole de-weighting; i is the deduplication quality of laser single blind hole, which are M1, M2, M3...M n , are all selected fixed values, X i The weight M of the single punching i The distance between the center of the corresponding deweighting range and the center of the deweighting workpiece rotor, s i For X i Corresponding to the range, duplicate blind holes are removed by selecting the appropriate position Xs within the range i To remove the appropriate deduplication vector; For any deweighting vector U within the specified range, by selecting different blind hole qualities M i Composition of different U i , making U=U1+U2+...+U n , thereby achieving the removal of the specified mass; According to X i The deduplication range determines the diameter R of the corresponding deduplication blind hole i , ensure that R i +R i-1 <(X i -s i )-(X i-1 +s i-1 ), so that different deduplication blind holes will not affect each other. The relationship between the deduplication weight and the radius is: Where ρ1 is the density of the rotor substrate, ρ2 is the density of the rotor surface layer; h1 is the de-weighting depth of the blind hole on the rotor substrate, and h2 is the de-weighting depth of the blind hole on the rotor surface layer; Step 2: Standard deduplication blind hole processing: first, evaporate and remove the oxide layer on the rotor surface, and then evaporate and remove the rotor substrate layer. The parameter table establishes the functional relationship between laser parameters and deduplication amount based on laser simulation and experiment. Based on this functional relationship, various laser processing parameters corresponding to the deduplication blind hole oxide layer and substrate layer are established, including laser power, spot radius, defocus, scanning speed, spot overlap rate, and scanning times. Step 3: Select corresponding deduplication blind hole parameters using the standard deduplication blind hole parameter library according to the shape, mass and density values ​​of the rotor; Step 4: Use the selected laser parameters to test several deweighting blind holes on the test piece with the same material as the rotor to be balanced, and correct the actual deweighting value. Before the trial deweighting, test several sets of selected values ​​on the same material sample, measure one for more than 5 times and take the average value to calibrate the parameters; Step 5: According to the real de-weighting value obtained after the de-weighting blind hole correction, the actual de-weighting position is corrected by the integral method to de-weight the rotor; Step 6: After deduplication, check the deduplication effect. If it is unqualified, return to step 5 to perform deduplication again until the deduplication effect is qualified.

2. The method according to claim 1, characterized in that: The standard deduplication blind hole parameter library described in step 1 must contain deduplication blind hole parameters and blind hole position ranges of various diameters and depths corresponding to various rotor materials and sizes, with a minimum deduplication weight of 0.1 mg and a minimum imbalance balance accuracy of 0.5 mg.mm.

3. The method according to claim 1, characterized in that: Step 2 The simulation process of laser processing parameters includes the following steps: 2.1) First, select an 8-node hexahedral unit, establish a 3D model of the deweighted material, and perform laser processing deweighting simulation to simulate the temperature field and thermal stress of the rotor. When the unit temperature reaches the vaporization temperature or the thermal stress between units of different materials is greater than the limit adhesion, the unit is considered to be removed; 2.2) According to the physical constants and basic equations, the initial temperature parameters and boundary conditions of the geometric body are defined for simulation. The output of the previous time step is used as the input of the next time step. The pulsed laser is input in the form of heat flux. The heat flux distribution function is shown in formula (2): Where A is the absorption rate, P is the average laser power, r is the spot radius, t is the thermal conduction time, v is the scanning speed, τ is the time relative to the full width of the pulse to reach half the maximum value, τ m is the time to reach the maximum power; 2.3) The laser heat source input is simplified to a point heat source on the spatial section and a pulse input on the time section; the load matrix list is used to provide the corresponding heat flux loads at different positions and times. The heat is radiated to the surface of the oxide layer and diffused to the bottom of the oxide layer and the substrate through heat conduction, following the heat balance equation in Fourier's law, as shown in formula (3): Where, the temperature field distribution is T(x,y,z,t), the heat flux distribution is q(x,y,z,t), the energy input is g(x,y,z,t), κ is the thermal conductivity, and α=κ / ρc is the thermal diffusion coefficient; 2.4) During the laser cleaning process, there are differences in the thermal expansion coefficient γ and Young's modulus Y between the oxide film layer inside the material and the base material, which leads to the generation of thermal stress, as shown in formula (4): Where u is the displacement, γ is the thermal expansion coefficient, and Y is Young's modulus; The adhesion force f between the substrate layer and the oxide layer is Z is the distance between the two layers, A 12 is the Hamaker coefficient between the oxide film and the substrate. When the thermal stress is greater than the adhesion force, the oxide film is removed; 2.5) Based on the above model, simulation analysis and experiments are combined to determine the laser processing parameters.

4. The method according to claim 1, characterized in that: The maximum deweighting mass selected in step 3 shall not exceed 1% of the rotor mass, and the position ranges of the selected standard deweighting blind holes shall not overlap.

5. The method according to claim 1, characterized in that: The parameters for processing the corresponding blind holes in step 4 are based on simulation formulas. If the difference between the actual weight removal value and the simulation value is too large to meet the actual weight removal range requirements, other blind hole values ​​in the parameter table are selected for assistance.

6. The method according to claim 1, characterized in that: The integration method in step five includes using a numerical integration algorithm to calculate a correction value of the actual deduplication position.

Images