# A Novel Rotor Profile Error Tracing and Compensation Strategy for High Precision Machining of Screw Rotor Based on Trial Cutting of Limited Samples

Authors: ZhiHuang Shen; Bin Yao; BinQiang Chen; Wei Feng; XiangLei Zhang
Journal: Shock and Vibration
 (2015)

---

## Abstract

The machining precision plays an important role in the operation reliability and service life of screw rotors. However, the actual machined rotor profile of screw rotors is different from its theoretical profile due to the errors of the machine tool. This paper proposes a novel method of error tracing and compensation to reduce the machined errors of rotor profile on the basis of the limited sample trials, and the method is based on a matrix of error compensation. The errors of rotor profile are obtained based on limited sample trials machining of screw rotors,and fitted into piecewise smooth data. A matrix of error compensation is established to compute the errors of rotor profile to generate a compensated rotor profile. The compensated rotor profile is then used to regenerate forming tool and a new rotor product is processed on the same machine tool. And the errors of new rotor profile are smaller and can be reduced within (−0.01 mm, 0.01 mm) after compensations. Finally, the actual machining examples illustrate that the method of error compensation can not only satisfy the machining requirement of high-precision rotors, but also has the characteristics of good stability and applicability.


---

## Body

## 1. Introduction

The advantages of twin-screw compressor include high reliability, good adaptability, ease of operation and maintenance, and the ability to transport various mixed liquids. Therefore, the twin-screw compressor has been widely used in automotive, compression, refrigeration, metallurgy, buildings, and other industrial fields. The design, manufacturing, and operational maintenance of screw rotor have been a topic that attracts more and more attention. As for bearings [1], gears [2], and shafts [3–5], their operation conditions can be investigated by many advanced monitoring technologies. The screw rotor is a kind of nonstandard part whose manufacturing quality is not so easy to guarantee. And a screw rotor of inferior quality has a negative effect on the fault diagnosis, condition monitoring, and operation maintenance of screw rotors. Therefore, it is of significant importance to ensure high machining precision of screw rotors.A pair of conjugate screw rotors is the most crucial component of the twin-screw compressor, whereby the machining precision of the screw rotors can greatly affect performance of the compressor. Screw rotors are treated as helical gears with nonparallel and nonintersecting or crossed axes. And the machining of helical gears can be applied to the machining of screw rotors. Screw rotors are machined today mainly by grinding or milling. Usually two stages are needed: the first is the roughing, where the workpiece is machined to its approximate size; the second is the finishing, where the rotors are machined to their final dimensions.Due to a variety of internal and external factors of processing system of the machine tool, the machining errors seriously affected the precision and quality of rotor profile. There are two common technologies to improve the accuracy of rotor profile: direct compensation method and indirect compensation method. The machine tool is the compensation object for the indirect compensation method which improves the processing precision of rotor profiles by enhancing the precision of machine tool. Meanwhile, rotor products are the object for the direct compensation which makes direct error compensation to rotor profile. Indirect compensation method is realized by creating a new error artificially to offset the original one, which can decrease the machining error and improve the accuracy. However, this method has three shortcomings: lack of an ideal modeling, obvious technical defects in hardware error compensation, and immature software error compensation. For the direct compensation method, recent researches mainly focus on the change of rotor profile under different errors. There are many errors being discussed, including machine tool errors, parameter errors of rotor structure, and wear of forming tool. Therefore, the theory on direct compensation method of rotor profiles could be further investigated.Stosic et al. [6] presented a method of general geometrical definitions of screw machine rotors and their manufacturing tools. Wu et al. [7] proposed a design method of forming wheel edge shape for screw rotors machining when arbitrary rotor profile is given. A contact line equation between the rotor tooth surface and the forming wheel revolving surface is derived based on the meshing theory, and then the calculation method for theoretical edge shape of forming wheel is proposed, which is applicable to both analytical and discrete forms of the rotor profile. These methods are extensively applied in this study to calculate forming tool profile and normal vector of each point of rotor profile.Guo and Tang [8] presented a model-based simulation to analyze the influence of manufacturing process variations on rotor profile, this simulation serves as a virtual machine tool to machine rotors. By varying machine parameters within their tolerance bands, the simulation predicts rotor profile errors caused by various process parameter errors. Han and Meng [9] presented a new ZK worm mathematical model which involves several dressing parameters of cone grinding wheel, and then the effects of worm parameters on the normal tooth profile error along with the changing of the grinding wheel radius are analyzed. Stosic [10, 11] introduced three possible errors of forming tool set errors, namely, mismatch of the angle between the tool and rotor shafts, mismatch of the centre distance between the tool and the rotor, and mismatch in the axial tool position, to quantify the tool set errors as well as compensate the errors of rotor profile. Wu [12] also defined three error factors based on the coordinate system established, and then the effects of rotor profile error along with the changing of these three error factors are analyzed. These studies can only compensate the errors of rotor profile which distribute regularly. These methods focus on finding the changing law of rotor profile under different errors, and limited compensation case only fits the situation when errors distribute regularly, not suitable for random distribution.In the study presented here, a novel method of error compensation is proposed to reduce the errors of rotor profile. The errors of rotor profile are obtained by trial machining of rotors and then fitted into piecewise smooth data. A matrix of error compensation is established firstly, and then the errors of rotor profile are computed by this matrix to generate a compensated rotor profile. The compensated rotor profile is used to recalculate a new forming tool profile for machining rotor. The errors of new machined rotor profile are smaller than before and can be reduced within (−0.01 mm, 0.01 mm) after compensations. Finally, the actual machining examples illustrate the method of error compensation which can not only satisfy the machining requirement of high-precision rotors, but also has the characteristics of good stability and applicability.


## 2. The Calculation of Forming Tool Profile

As is shown in Figure1, a coordinate system of the forming tool machining for screw rotor is presented. The coordinate system Sr is attached to the rotor, rotating around the zr-axis, while the coordinate system Sg is affixed to the forming tool, rotating around the zg-axis. The center distance A is defined as the minimum distance between the zr-axis and the zg-axis, and Σ is the shaft angle between two axes. The coordinate transformation from Sr to Sg can be derived as(1)xgygzg=1000cosΣ-sinΣ0sinΣcosΣxryrzr+-a00.Figure 1
Coordinate system of forming tool for machining screw rotor.There is a contact line between the forming tool contour surface and the rotor helicoids during machining process. When the contact line spirally rotates around thezr-axis, rotor helicoids are attained. And when the contact line rotates around zg-axis, the contour surface of the forming tool would be shaped. The rotor profile consisted of discrete points, and the derivative of each discrete point of rotor profile is not given. To calculate the derivative, the method of accumulated chord length parameter, which can solve the problem of the large deflection of rotor profiles, is used in this paper. Assuming the coordinate of the discrete points is P={(xi,yi),i=1,2,…,n}, the equation of rotor profile can be obtained as(2)r0=xs,ys,where s is the accumulated chord length parameter. And then the derivative of rotor profile can be attained as(3)T=∂xs∂s,∂ys∂s.The helicoids of the rotors can be expressed by the following equation:(4)rs,θ=x,y,z=xscosθ-eyssinθ,exssinθ+yscosθ,pθ,where e=1 for the right-hand rotor, e=-1 for the left-hand rotor, θ is the rotation angle of the rotor, and p=Pz/2π is helix parameter, which means that the rotor translates a distance of lead Pz when it rotates one turn. The equation of helicoids of the rotor has been obtained, and the partial derivatives of the helicoids can be calculated as follows:(5)rs=∂xs∂scosθ-e∂ys∂ssinθ,e∂xs∂ssinθ+∂ys∂scosθ,0,rθ=-xssinθ-eyscosθ,exscosθ-yssinθ,p.And then the normal vector of the helicoids of the rotor can be presented as(6)n=nx,ny,nz=rs×rθT=pe∂xs∂ssinθ+∂ys∂scosθ,p-∂xs∂scosθ+e∂ys∂ssinθ,exs∂xs∂s+ys∂ys∂s.According to envelope principle, the relative velocityv between the rotor and the forming tool can be presented as(7)v=vx,vy,vz=ycosΣ+zsinΣ,x-acosΣ,x-asinΣ.In order to determine the relationship equation between the profile parameters and the rotational angle θ, the enveloping condition [13] is that the common normal vector n at the contact point is perpendicular to the relative velocity v direction between two conjugate surfaces and can be shown as follows:(8)v·n=-pcosΣexssinθ+yscosθ×e∂xs∂ssinθ+∂ys∂scosθ-pcosΣxscosθ-eyssinθ-a×∂xs∂scosθ-e∂ys∂ssinθ+esinΣxscosθ-eyssinθ-a×xs∂xs∂s+ys∂ys∂s-p2θsinΣe∂xs∂ssinθ+∂ys∂scosθ=0.Substituting the value ofs into (6) yields the corresponding value θ=fs, for the rotational angle. Substituting the value of s and θ back into (4), the equation of contact line can be expressed as(9)rcr=xcr,ycr,zcr=xscosfs-eyssinfs,exssinfs+yscosfs,pfs.Then transforming the equation of contact linercr from rotor coordinate system Sr to forming tool coordinate system Sg,(10)rcgT=xscosfs-eyssinfs-aexssinfs+yscosfscosΣ-pfssinΣexssinfs+yscosfssinΣ+pfscosΣ.The contour surface of the forming tool can be generated by rotating the contact line around its rotating axis, and the equation of forming tool can be presented as(11)rg=Zg,Rg=zcg,xcg2+ycg2.


## 3. The Error Compensation Method of Rotor Profile

### 3.1. Analysis of Error Source

The errors of rotor profile are mainly caused by the comprehensive errors of machine tool which can be divided into two major categories of external and internal errors according to its own attributes, as is shown in Figure2. The internal error is composed of geometric errors, force errors, thermal errors, wear of forming tool, and its own vibration. The external error mainly considers the influence of external environment including temperature, humidity, and adjacent vibration. Ignoring the external error, according to “machine tool, fixture, rotor, and forming tool” machining system, the internal error can be divided into more detailed categories, and the distribution of each error factor is roughly shown in Table 1.Table 1
The proportion of internal error factors.
Factors
Proportion
Geometric errors
35%
Thermal errors
18%
Force errors
12%
Wear of forming tool
10%
Deformation error of workpiece
8%
Setting error
3%
Fixture error
5%
Others
9%Figure 2
Classification of errors.


### 3.2. Method of Error Compensation

To improve the machining precision of rotor profile, this study presents a novel compensation method for rotor profile errors. In order to compare with the errors between different rotor profiles, the overall error of rotor profileη is introduced as an evaluation index:(12)η=∑i=1nδi21/2n,where δi is the error of each point of rotor profile.As is shown in Figure3, the compensation procedure can be divided into the following steps.(1)
According to the theoretical rotor profile, the profile of forming tool is first calculated. On the basis of this profile, a forming tool is dressed and then used to process a rotor workpiece to obtain a rotor product. Errors of rotor profile are obtained by measuring the machined rotor, and the overall errorη1 can be calculated.(2)
The measured error of rotor profile needs to be processed on account of its large fluctuation feature. According to the trend of error curveab, bc, cd, de, ef, fg, gh, hi, and ij, errors are divided into 9 sections. If the radius of probe is larger than that of curvature at the root of the rotor, the points are unable to be measured; see ef section. A linear transition can be used to fit the section ef. So the fitting equation of ab, bc, cd, de, ef, fg, gh, hi, and ij can be expressed as follows:(13)δi=0.00166i+0.00387,i∈1,10,-0.0016i+0.03639,i∈10,25,-0.000065i+0.00706,i∈25,106,0.00296i-0.31445,i∈106,117,-0.001i+0.148,i∈117,147,0,i∈147,169,0.01i-1.69,i∈169,174,-0.00286i+0.56302,i∈174,212,0.00115i-0.2879,i∈212,325,-0.01091i+3.63172,i∈325,335,0.000366i-0.13513,i∈335,340.(3)
The errors data has been reconstructed, and then the matrix of error compensation is established to generate a compensated rotor profile:(14)
r
s
=
x
s
1
y
s
1
x
s
2
y
s
2
⋮
⋮
x
s
i
y
s
i
⋮
⋮
x
s
n
-
1
y
s
n
-
1
x
s
n
y
s
n
1
0
0
1
=
δ
1
0
0
0
0
0
0
x
1
y
1
0
δ
2
0
0
0
0
0
x
2
y
2
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
0
0
0
δ
i
0
0
0
x
i
y
i
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
0
0
0
0
0
δ
n
-
1
0
x
n
-
1
y
n
-
1
0
0
0
0
0
0
δ
n
x
n
y
n
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
0
1
×
x
1
′
y
1
′
x
2
′
y
2
′
⋮
⋮
x
i
′
y
i
′
⋮
⋮
x
n
-
1
′
y
n
-
1
′
x
n
′
y
n
′
1
0
0
1
,wherers=(xsi,ysi) is the compensated rotor profile and r0=(xi,yi) is the theoretical rotor profile. n0=xi′,yi′ is the normal vector of theoretical rotor profile and can be calculated as follows: (15)n0=xi′,yi′=∂yi/∂s∂xi/∂s2+∂yi/∂s2,-∂xi/∂s∂xi/∂s2+∂yi/∂s2.(4)
Finally, the rotor profilers can be used to recalculate the forming tool profile by (2)–(11). After a new rotor is machined, the overall error of new rotor profile η2 is calculated. By comparing the values of η1 and η2, the errors after compensation are smaller than the errors before compensation, and that can be reduced within (−0.01 mm, 0.01 mm). Thus, the error compensation of rotor profile is achieved.Figure 3
Procedure of error compensation.


## 3.1. Analysis of Error Source

The errors of rotor profile are mainly caused by the comprehensive errors of machine tool which can be divided into two major categories of external and internal errors according to its own attributes, as is shown in Figure2. The internal error is composed of geometric errors, force errors, thermal errors, wear of forming tool, and its own vibration. The external error mainly considers the influence of external environment including temperature, humidity, and adjacent vibration. Ignoring the external error, according to “machine tool, fixture, rotor, and forming tool” machining system, the internal error can be divided into more detailed categories, and the distribution of each error factor is roughly shown in Table 1.Table 1
The proportion of internal error factors.
Factors
Proportion
Geometric errors
35%
Thermal errors
18%
Force errors
12%
Wear of forming tool
10%
Deformation error of workpiece
8%
Setting error
3%
Fixture error
5%
Others
9%Figure 2
Classification of errors.


## 3.2. Method of Error Compensation

To improve the machining precision of rotor profile, this study presents a novel compensation method for rotor profile errors. In order to compare with the errors between different rotor profiles, the overall error of rotor profileη is introduced as an evaluation index:(12)η=∑i=1nδi21/2n,where δi is the error of each point of rotor profile.As is shown in Figure3, the compensation procedure can be divided into the following steps.(1)
According to the theoretical rotor profile, the profile of forming tool is first calculated. On the basis of this profile, a forming tool is dressed and then used to process a rotor workpiece to obtain a rotor product. Errors of rotor profile are obtained by measuring the machined rotor, and the overall errorη1 can be calculated.(2)
The measured error of rotor profile needs to be processed on account of its large fluctuation feature. According to the trend of error curveab, bc, cd, de, ef, fg, gh, hi, and ij, errors are divided into 9 sections. If the radius of probe is larger than that of curvature at the root of the rotor, the points are unable to be measured; see ef section. A linear transition can be used to fit the section ef. So the fitting equation of ab, bc, cd, de, ef, fg, gh, hi, and ij can be expressed as follows:(13)δi=0.00166i+0.00387,i∈1,10,-0.0016i+0.03639,i∈10,25,-0.000065i+0.00706,i∈25,106,0.00296i-0.31445,i∈106,117,-0.001i+0.148,i∈117,147,0,i∈147,169,0.01i-1.69,i∈169,174,-0.00286i+0.56302,i∈174,212,0.00115i-0.2879,i∈212,325,-0.01091i+3.63172,i∈325,335,0.000366i-0.13513,i∈335,340.(3)
The errors data has been reconstructed, and then the matrix of error compensation is established to generate a compensated rotor profile:(14)
r
s
=
x
s
1
y
s
1
x
s
2
y
s
2
⋮
⋮
x
s
i
y
s
i
⋮
⋮
x
s
n
-
1
y
s
n
-
1
x
s
n
y
s
n
1
0
0
1
=
δ
1
0
0
0
0
0
0
x
1
y
1
0
δ
2
0
0
0
0
0
x
2
y
2
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
0
0
0
δ
i
0
0
0
x
i
y
i
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
0
0
0
0
0
δ
n
-
1
0
x
n
-
1
y
n
-
1
0
0
0
0
0
0
δ
n
x
n
y
n
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
0
1
×
x
1
′
y
1
′
x
2
′
y
2
′
⋮
⋮
x
i
′
y
i
′
⋮
⋮
x
n
-
1
′
y
n
-
1
′
x
n
′
y
n
′
1
0
0
1
,wherers=(xsi,ysi) is the compensated rotor profile and r0=(xi,yi) is the theoretical rotor profile. n0=xi′,yi′ is the normal vector of theoretical rotor profile and can be calculated as follows: (15)n0=xi′,yi′=∂yi/∂s∂xi/∂s2+∂yi/∂s2,-∂xi/∂s∂xi/∂s2+∂yi/∂s2.(4)
Finally, the rotor profilers can be used to recalculate the forming tool profile by (2)–(11). After a new rotor is machined, the overall error of new rotor profile η2 is calculated. By comparing the values of η1 and η2, the errors after compensation are smaller than the errors before compensation, and that can be reduced within (−0.01 mm, 0.01 mm). Thus, the error compensation of rotor profile is achieved.Figure 3
Procedure of error compensation.


## 4. Example

The method above has been successfully applied to many screw compressor manufacturers. Taking a machining example of the male rotor and female rotor which were supplied from a screw compressor manufacturer, the parameters of the rotor are given in Table2. Figure 4 shows the cross-section contour of the male rotor and female rotor; they are composed of discrete points. The forming tool profiles of male and female rotors are calculated by (2)–(11) and shown in Figure 5.Table 2
The parameters of rotor.
Type
Revolving direction
Number of teeth
Outer diameter (mm)
Root diameter (mm)
Lead (mm)
Helical angle (deg)
Material
Male
Right-hand
4
68.25
39.7635
128
45
42CrMo
Female
Left-hand
5
54.653
27.248
160Figure 4
Cross-section contour of rotors.
(a)
Male rotor
(b)
Female rotorFigure 5
Forming tool profile of rotors.
(a)
Male rotor
(b)
Female rotorFigure6 shows the process of machining rotors, and the machined rotor products are shown in Figure 7. To assess the precision of rotor profile, the end face of rotor is measured by 3D CMM, as is shown in Figure 8.Figure 6
The machining of rotors.
(a)
(b)Figure 7
Rotor products.Figure 8
Measurement of rotors.The measurement results of male and female rotors profile are shown in Figures9 and 10, in which the errors are magnified by 50 times, and the tolerance of rotor profile is (−0.01 mm, 0.01 mm). Figure 9 shows the rotor profiles after trial machining, of which the overall errors are ηmale_1 = 0.06675 mm and ηfemale_1 = 0.1839 mm under the effect of comprehensive errors of machine tool. The method of error compensation is used to reduce the errors as is shown in Figure 10; the overall errors of male rotor and female rotor are ηmale_2 = 0.00785 mm and ηfemale_2 = 0.0127 mm. By comparing the value of ηmale_1 and ηmale_2 and ηfemale_1 and ηfemale_2, the errors after compensation are far less than the errors before compensation. Thus, the method of error compensation can satisfy the machining requirement of high-precision rotors.Figure 9
Rotor profile before compensation.
(a)
Male rotor
(b)
Female rotorFigure 10
Rotor profile after compensation.
(a)
Male rotor
(b)
Female rotorIn the mass production of screw rotors, the precision of rotor profiles is tracked. Figure11 shows 20 consecutive machined products of the same rotor profile, and the overall errors distribute within (0.005 mm, −0.0094 mm). Figure 12 shows the comparison of profile errors of 5 different rotor products before and after compensation. The result shows that, after compensation, all errors of 5 different rotor products are reduced. Figures 11 and 12 illustrate that the method of error compensation is stable and applicable.Figure 11
Overall errors of 20 consecutive machining rotor products.Figure 12
Comparison of profile errors of 5 different rotor products.


## 5. Conclusions

(1)
This paper investigates the envelope condition to derive an equation of the contact line, and then the contact line is used to generate the forming tool profile.
(2)
After analyzing the error sources of machine tool, a method of error compensation is presented to improve the accuracy of rotor profile. The errors of rotor profile are obtained by small sample machining trials of rotors and then fitted into piecewise smooth data. A matrix of error compensation is established firstly, and then the errors of rotor profile are computed by this matrix to generate a compensated rotor profile. The compensated rotor profile is then used to regenerate forming tool and a new rotor product is processed on the same machine tool. And the errors of new rotor profile are smaller and can be reduced within (−0.01 mm, 0.01 mm) after compensations.
(3)
Finally, the actual five machining examples in the high precision, stability, and applicability illustrate that the method of error compensation can not only satisfy the machining requirements of high-precision rotors, but also has the characteristics of good stability and applicability.



---

*Source: 978325-2015-09-27.xml*