Next Article in Journal
Investigation of the Efficiency of Shielding Gamma and Electron Radiation Using Glasses Based on TeO2-WO3-Bi2O3-MoO3-SiO to Protect Electronic Circuits from the Negative Effects of Ionizing Radiation
Next Article in Special Issue
Meshless Chebyshev RPIM Solution for Free Vibration of Rotating Cross-Ply Laminated Combined Cylindrical-Conical Shells in Thermal Environment
Previous Article in Journal
Use of New Green Bitumen Modifier for Asphalt Mixtures Recycling
Previous Article in Special Issue
Rough Surface Characterization Parameter Set and Redundant Parameter Set for Surface Modeling and Performance Research
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Measurement Solution of Face Gears with 3D Optical Scanning

1
State Key Laboratory of High Performance Complex Manufacturing, Central South University, Changsha 410083, China
2
College of Mechanical and Electrical Engineering, Central South University, Changsha 410083, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(17), 6069; https://doi.org/10.3390/ma15176069
Submission received: 21 July 2022 / Revised: 19 August 2022 / Accepted: 30 August 2022 / Published: 1 September 2022

Abstract

:
Gears are usually measured by the contact metrology method in gear measuring centers or coordinate measuring machines. Recently, three-dimensional (3D) optical scanning, a non-contact metrology method, has been applied in the industry as an advanced measurement technology mainly due to its high efficiency. However, its applications to gears with complicated geometry, such as face gears, are still limited due to its relatively low accuracy and the void of related measurement solutions. In this work, an accurate measurement solution with 3D optical scanning is proposed for the tooth surface deviations of orthogonal face gears. First, point cloud collection is carried out by the 3D scanner. Furthermore, the measurement solution is implemented with a three-stage algorithm by aligning point clouds with the design model. Subsequently, 3D modeling is studied by numbering the points and reconstructing the real tooth surfaces. An example with a measurement experiment and loaded tooth contact analysis is given to show the validity of the proposed method.

1. Introduction

Face gears [1,2,3,4] have been applied in some transmission systems due to the advantages of smooth transmission, large coincidence, and small footprint. Compared with the bevel gear drives, face gear drives have smaller axial forces and lower sensitivity to axial installation errors. Moreover, the support apparatus of face gear drives can be designed in a relatively more compact way than spiral bevel gear drives, and it can significantly save weight. The above strengths make face gears receive attention and research. Some car companies, such as Audi, have developed the advanced technology of face gears and applied it to automobile differentials. In the aerospace industry, face gears have been successfully used in helicopter transmission systems. More types of face gears are studied nowadays. According to the macroscopic form of the tooth surface, the types of face gears can be divided into spur face gears, helical face gears, conical face gears, worm face gears, asymmetric helical face gears [5,6], and so on. The performance of face gear drives is constantly improved, and the application scenarios are wider.
In terms of processing, face gears can be machined by grinding, CNC plunge milling, hobbing, and other processing methods [7,8,9]. The spur face gears, a simple type of face gears, can reach the desired accuracy level in the aerospace industry by grinding with “grinding-measurement-correction” [10,11]. The precision grade is evaluated by the measurement, so the reliable measurement of face gears is the basis of the correction of the machine tool settings in the processing. The measurement method of face gears generally adopts the contact measuring machines [12,13,14], such as coordinate measuring machines (CMMs) and gear measuring centers (GMCs). CMMs are universal measuring tools and can measure various workpieces. GMCs are specialized machines for measuring gears, which are made by adding rotary tables to CMMs. They both utilize the touch signal emitted by a probe when the tip ball comes into contact with the surface of a workpiece and records the probe coordinates at that time. The recorded coordinates are not valid points on the surface of the workpiece due to the signal delay, the inaccurate measurement coordinate system, and the radius of the probe, so the measured data need to be compensated further by the numerical algorithm.
The real face gears are the workpieces after processing, and they usually contain machining errors, which will affect the meshing performance. During the meshing process, the tooth surfaces of face gears and pinions will contact at a point at one moment, and the contact path will generate by connecting all contact points at every moment. Tooth surface deviations will alter the original tooth surface morphology, which will change the original contact path and further has a significant impact on the meshing performance of gear drives. In the industry, the meshing performance of real face gears can usually be analyzed by roll testing experiment, but this method is costly and timely. Therefore, some scholars have conducted tooth contact analysis for the real tooth surfaces, reconstructed with the measured data using CMMs or GMCs. Wang et al. [15] used the compensated measured data to reconstruct the tooth surface of face gears and conducted the digital tooth contact analysis using the proposed robust algorithm. Lin et al. [16] presented a least rotation angle method and the improved quad-tree search algorithm to compute the contact pattern and transmission error of a single flank fitted by the data measured on GMCs. Din et al. [17] proposed a collaborative manufacturing system for hypoid gears based on data-driven programming that the data are collected using CMMs, and the machine setting is modified by the tooth geometric and meshing performance.
The above researchers all adopted the contact metrology method. Although the contact metrology method, with the aid of compensation technology, can reach high accuracy, it has two disadvantages stated as follows.
(1)
The border area of the tooth surface cannot be measured accurately since other locations would be touched. Moreover, the accuracy of the extrapolation method is insufficient when reconstructing the tooth surfaces. Hence, it is impossible to accurately determine whether edge contact occurs during the meshing of the real face gear drives. Edge contact refers to the meshing of the border areas of the tooth surface between face gears and pinions, which leads to stress concentration and meshing impact [18,19]. It will cause obvious vibration and noise during the meshing and reduces the service life of gears.
(2)
The efficiency is low when measuring multiple teeth since every point needs to be contacted.
3D optical scanning, as a non-contact metrology method, has become a new trend in measuring the macroscopic morphology of surfaces with the emergence and development of computer vision [20,21,22]. Recently, this technique has been studied to inspect and analyze gear quality. Chen et al. [23] proposed an optical inspection method for tooth surfaces using the projection moire technique and measured the deviation based on the reference gear tooth profiles measured by CMM. Gonzalez-Perez et al. [24] performed tooth contact analysis of cylindrical gears reconstructed from point clouds. Urbas et al. [25,26] studied improved methods for optical measurement accuracy of spur gears. Three-dimensional (3D) optical scanning can be a concern by scholars due to the advantages as follows.
(1)
As a non-contact metrology method, it is more efficient than the contact metrology method when measuring multiple teeth.
(2)
The collected data using this method obtain the complete geometry of a tooth, including the top land surface, the bottom land surface, and both sides of the holistic tooth surface.
However, 3D optical scanning for measurement has not been widely used in the gear industry due to the low measurement accuracy. Its accuracy typically varies from 10 μm up to 100 μm [27,28,29,30]. The measurement errors come from two aspects. One is the systematic errors, and the other is the alignment errors between the measured workpiece and the design model. The former can only be reduced by upgrading the hardware specifications, while the latter can be reduced by alignment using the numerical algorithm. Thus, it is possible to achieve an accuracy of several microns with some numerical algorithms for compensation. Moreover, the previous research objects are cylindrical gears, while geometric characteristics of face gears are different, and the macroscopic morphology of tooth surfaces is extremely complex.
In order to solve the problem of contact and non-contact metrology method and meet the increasing need for measuring holistic geometric information, a 3D optical scanning method is adopted for measuring the tooth surface of orthogonal face gears in this work, and the measurement solution is presented to ensure the measurement accuracy. In Section 2, the method of point cloud collection is introduced, and point clouds are processed by filtering and down-sampling. In Section 3, the measurement solution for tooth surfaces of face gears with 3D optical scanning is proposed to obtain the tooth surface deviations based on the alignment, and three-stage algorithm (3SA) is presented to align three important features between point clouds and the design model of face gears. In Section 4, the 3D modeling method is further investigated. An example is given to show the validity of the measurement solution with 3D optical scanning, and loaded tooth contact analysis (LTCA) for the real tooth surfaces of face gears by finite element method (FEM) is conducted to show the influence of the deviations on meshing performance in Section 5.

2. Point Cloud Collection

Point clouds containing the holistic geometry information of face gears can be collected at a rate of millions of points per second by 3D optical scanning. A way to implement 3D optical scanning is by binocular fringe projection method [31,32,33,34], as shown in Figure 1. The main principle is to project the fringes generated by the system onto the surface of objects in the corresponding order. Binocular cameras are applied to take images, and the original fringes will be found to be deformed due to the height change of the surface of objects. Then, the phase information of the surface of objects can be solved by the phase-shifting algorithm with the collected data of the deformed fringes. According to the phase information, the coordinate value of each point in the space can be calculated by triangulation, and then the point cloud data of the surface can be obtained. The binocular fringe projection method has been widely studied, and the detailed implementation methods will not be further described in the paper.
Currently, the available 3D scanning technology has been integrated into commercial equipment by some companies. Choosing the proper scanner equipment is a requirement for obtaining high-quality point clouds. Generally, the higher the spatial resolution and measurement accuracy of the equipment, the higher the quality of the collected point clouds of face gears. This work takes ATOS 3D Scanner to collect point clouds, as shown in Figure 2a. ATOS 3D Scanner is an industrial optical system with a charge-coupled device camera at both ends. The measured data by it can be transmitted to the computer, and the coordinate information of point clouds can be shown in the special software. The main system configuration is shown in Table 1.
When performing the operation, first ensure that the experiments are in a stable working environment, such as a stable light and warm-up environment, in which the measurement results of the 3D scanner are not affected by external factors. Then, the scanner is set up on the stationary support, and the settings of equipment are adjusted, such as lens position and aperture size. Subsequently, the scanner needs to be calibrated with the calibration plates. Through calibration, the system can accurately calculate the position of the scanning equipment relative to the scanning objects. Therefore, accurate calibration is a prerequisite for measurement accuracy. Moreover, the white developers need to be sprayed onto the face gear tooth surfaces to make them diffuse reflection for a better quality of point clouds. Since the bottom surface of face gears is the processing reference with high accuracy, face gears can be placed flat on the workbench when collecting point clouds, as shown in Figure 2b. Finally, the scanner system can be operated to take pictures of the workpieces, and the system would automatically obtain the coordinate data of the workpiece so that it obtains the point cloud data. In Figure 2c, a workpiece of face gear is being scanned by the 3D scanner, and the point clouds can be seen in the scanner system.
Subsequently, point clouds would require data processing. As shown in Figure 3a, the unprocessed point clouds will contain unnecessary point data. Only the points on the tooth surface are needed, and the other points should be removed. Here take the conditional filters to remove the outliers. Since the table is horizontal, the high area is the point cloud of tooth surfaces, and the low area is the point cloud of the conical surface and upper surface. When setting the appropriate z-value, the two different point clouds could be distinguished. The needed point cloud of tooth surfaces is shown in Figure 3b. If multiple point clouds need to be stitched together, the noise points will be generated, and further noise point reduction is needed with filter algorithms, such as the Gaussian filter [35]. Subsequently, the data volume of the point cloud can be reduced by the down-sampling method [36], as shown in Figure 3c. The different grid boxes in the space are divided, and a specific density of points in the grid is taken.

3. The Measurement Solution

The point cloud collected in Section 2 is in the measurement coordinate system (MCS), which is different from the design coordinate system (DCS). DCS is a fixed coordinate system determined in the design process of face gears, while MCS is the coordinate system set randomly in the measurement process. In order to obtain the tooth surface deviations of face gears, the point cloud containing the geometry information must be compared to the design model in the same coordinate system. The tooth surface deviation of one point can be represented as
d e v i = ( p m i p t i ) n i
where p t i and n i are the coordinates and the unit normal of one point on the theoretical tooth surface in DCS, respectively. p m i is the measured point which is in the direction of n i . Before computing the tooth surface deviations, alignment is required, and the geometric features between the point cloud and the design model are aligned, which is mathematically equivalent to transforming the coordinate system. The alignment is taken into three stages: the alignment for the top land surface, the center of circles, and the tooth surface, as shown in Figure 4, and 3SA is presented as follows.
(1)
The first stage
The top land surface of the designed orthogonal face gears is a plane, which is parallel to the xoy plane in DCS, and the points on it are the same z-value. The collected points on the top land surface will not be the same z-value in MCS for the measurement and machining errors. Hence, the plane of the measured top land surface is fitted by the following equation as
z ls = a x ls + b y ls + c
The subscript ls represents the top land surface of the workpiece, and a, b, c are the coefficients. The solution can be obtained by solving the overdetermined systems of equations as
A X = b X = ( A A T ) 1 A T b
where
X = ( a b c ) , A = [ x 1 y 1 1 x 2 y 2 1 x 3 y 3 1 x n y n 1 ] , b = ( z 1 z 2 z 3 z n )
( x 1 , y 1 , z 1 ) ,   ( x 2 , y 2 , z 2 ) ,   ( x 3 , y 3 , z 3 ) ,   ,   ( x n , y n , z n ) are the coordinates of the points on the top land surface of the workpiece. Then, the point cloud is translated along the z-axis as
Δ z = z ls z t
where z t is the value of the theoretical top land surface along the z-axis in DCS.
(2)
The second stage
The center of circles Ot of the designed face gears is zero point, while the center of circles Op of the point cloud is set roughly. Hence, they are aligned by the pattern search technique [37,38]. The maximum and the minimum radius in the point cloud on the xoy plane are respectively Rrmax and Rrmin, and their coordinates are respectively (xrmax, yrmax, zrmax) and (xrmin, yrmin, zrmin). The maximum and the minimum radius of the designed face gears are Rtmax and Rtmin, respectively. The zero point in MCS is set on the inner side of the point cloud, so the initial value for the translation vector of the point cloud is set by the following equation as
I n i   =   ( ( R r max R tmax ) R r max x r max , ( R r max R tmax ) R r max y r max , 0 )
The following optimization objective is established as
min f = norm ( R r max R tmax , R r min R tmin )
Then the optimal solution is searched as follows.
Step 1. Set the initial translation vector q 0 = I n i , the initial step δ 0 = 1 , the contraction factor α = 0.6 , and the tolerance error ε = 10 6 .
Step 2. Let p = q k   and   j = 1 .
Step 3. If f ( p + δ k e j )   <   f ( p ) , p = p + δ k e j and turn to Step 5. Otherwise, turn to Step 4.
Step 4. If   f ( p     δ k e j )   <   f ( p ) , p = p     δ k e j and turn to Step 5. Otherwise, turn to Step 5.
Step 5. If j   <   n ,   j = j + 1 and turn to Step 3. Otherwise, q k + 1 = p and turn to Step 6.
Step 6. If   f (   q k + 1 )   <   f (   q k ) , let p = q k + 1 + ( q k + 1     q k ) , δ k + 1 = α δ k , k = k + 1 , j = 1 and turn to Step 3. Otherwise, turn to Step 7.
Step 7. If f (   q k )   <   ε , stop computation, and q k is the optimal solution. Otherwise, turn to Step 8.
Step 8. If q k + 1 = q k , let δ k + 1 = α δ k , k = k + 1 and turn to Step 2. Otherwise, let q k + 1 = q k , δ k + 1 = δ k , k = k + 1 and turn to Step 2.
The optimal solution is obtained through the above steps, as
q * = ( q 1 * , q 2 * , 0 )
where the superscript * represents the optimal solution. Then, the point cloud is translated according to the optimal solution as follows.
{ Δ x = q 1 * Δ y = q 2 *
(3)
The third stage
In the third stage, the points after being processed in the second stage will be aligned with the theoretical tooth surface. The 5 × 9 grid on each side of the measured tooth surface of a tooth is selected randomly. The grids are translated and rotated around the z-axis. The transformed grids are represented as
[ p tr 1 ] = M t r [ p mg 1 ]
where
p t r = [ x tr y tr z tr ] , p mg = [ x mg y mg z mg ] , M tr = [ cos α sin α 0 d x sin α cos α 0 d y 0 0 1 d z 0 0 0 1 ]
p mg represents the points on the grids before the transformation, Mtr represents the transformation matrix, and p tr represents the point on the transformed grids. α is the rotation angle around the z-axis, and dx, dy, and dz are the translation values along the x-, y-, and z-axis, respectively. They are all chosen as the optimization variables. The reason for the second translation here is that the data processed in the first two stages will contain systematic errors from the top land surface and center of circles, so the tooth surfaces need to be further aligned to reduce the systematic errors.
Then, the rotation radius r and height h of each point of grids can be obtained after rotating the grids around the z-axis, and the points are mapped to the mapping plane. Every point on the theoretical tooth surface with the same rotation radius and height is calculated according to the tooth surface equation, as follows.
{ r ( g , i , j ) = x tr ( g , i , j ) ( v a r ) 2 + y tr ( g , i , j ) ( v a r ) 2 h ( g , i , j ) = z tr ( g , i , j ) ( v a r ) , g = 1 , 2 ; i = 1 , 2 , , 5 ; j = 1 , 2 , , 9
{ p x ( g , i , j ) ( u v 1 , u v 2 ) 2 + p y ( g , i , j ) ( u v 1 , u v 2 ) 2 = r ( g , i , j ) p z ( g , i , j ) ( u v 1 , u v 2 ) = h ( g , i , j ) , g = 1 , 2 ; i = 1 , 2 , , 5 ; j = 1 , 2 , , 9
where the different values of the subscript g represent the grid on the different tooth surfaces of a tooth. The subscript i and j represent the point number on the grid. ( p x ( g , i , j ) , p y ( g , i , j ) , p z ( g , i , j ) ) represents the coordinate of the point on the theoretical tooth surface. var represents the unknown variables (α,dx,dy,dz). (uv1,uv2) represents two unknown variables of the tooth surface equation. The tooth surface equation of face gears has been derived in our previous literature [39], and here it is not stated in detail. The vector between the point p tr on the transformed grid and the corresponding theoretical point p tp can be yielded as
ε p ( g , i , j ) = p tr ( g , i , j ) ( v a r ) p tp ( g , i , j ) ( v a r )
In order to achieve the best alignment between the transformed grid and the theoretical tooth surface, a nonlinear least square minimization optimization is established, and it is represented as
m i n f ε   =   1 2 g , i , j ε p ( g , i , j ) ( v a r ) T ε p ( g , i , j ) ( v a r )
The above over-determined nonlinear equation is expanded by second-order Taylor to transform it into a trust region problem, as
min O ε k ( d ) = f ε k + f ε k T d + 1 2 d T 2 f ε k d
where d is the optimal step of the iteration. f ε k is the value of f ε in the kth iteration. f ε k and 2 f ε k are the first and second gradient, respectively. For convenience, f ε k is marked as gk. Since the calculation for 2 f ε k is too complex, the approximate Hessian matrix B k is used to replace, as
B k = J k T J k
where J k is the Jacobi matrix. Subsequently, the trust region dogleg method is used to solve the optimization function. The optimization problem is represented as
d = arg min f ε k + g k T d + 1 2 d T B k d s . t . d   Δ k
where Δ k is the trust region radius. The solution to Equation (18) can be given as
d * = { τ C , 0 τ 1 C + ( τ 1 ) ( G N + C ) , 1 < τ 2
where
C = g k T g k g k T B k g k g k
G N = ( B k ) 1 g k
C is the Cauchy step, and GN is the Gauss-Newton step. τ can be solved by the following scalar quadratic equation as
C + ( τ 1 ) ( G N + C ) 2 = Δ k
The ratio ρ k is used to measure the approximation between Equations (15) and (18), as follows.
ρ k = f ε ( v a r k ) f ε ( v a r k + d ) O ε k ( 0 ) O ε k ( d )
where v a r k is the value of optimization variables in the kth iteration. When 0.9   <   ρ k     1 , the approximation is well and Δ k is increased. When 0 . 1   <   ρ k     0.9 ,   Δ k remains the same. When ρ k     0.1 ,   Δ k is reduced.
When the optimal solution is solved, the point cloud can be transformed with
M tr * = [ cos α * sin α * 0 d x * sin α * cos α * 0 d y * 0 0 1 d z * 0 0 0 1 ]

4. 3D Modelling

To better visualize the tooth surface deviations using 3D software and further carry out contact simulation by FEM [40,41,42,43], the 3D model of real five-pairs-of-teeth face gears needs to be modeled. For complex surface reconstruction, non-uniform rational b-spline (NURBS) [44,45,46] is the common modeling method. However, the disorderly scattered points make it difficult to reconstruct the face gear tooth surface accurately. Therefore, it is necessary to determine the regularly arranged points in the range of the tooth surface and label the points with serial numbers.
As shown in Figure 5a, the theoretical points arranged equidistantly within the range of tooth surface are calculated first according to the equation of face gear tooth surface [39], and (a,b) is used to record the sequence number of points. The nearest points cp from the point cloud are found by calculating the minim Euler distance E from the theoretical points tp(a,b), and cp will be numbered as same as tp(a,b), as follows.
min E = c p ( a , b ) t p ( a , b ) 2
The point cloud with numbered points is obtained, as shown in Figure 5b. Subsequently, the tooth surface is represented as a bicubic NURBS surface, as
S ( u , v ) = i = 0 m j = 0 n N i , 3 ( u ) N j , 3 ( v ) ω i , j d i , j i = 0 m j = 0 n N i , 3 ( u ) N j , 3 ( v ) ω i , j
where
{ N i , 0 ( u ) = { 1 , u i u < u i + 1 0 , otherwise N i , p ( u ) = u u i u i + p u i N i , p 1 ( u ) + u i + p + 1 u u i + p + 1 u i + 1 N i + 1 , p 1 ( u ) , p = 1 , 2 , 3
{ N j , 0 ( v ) = { 1 , v j v < v j + 1 0 , otherwise N j , p ( v ) = v v j v j + p v j N j , p 1 ( v ) + v j + p + 1 v v j + p + 1 v j + 1 N j + 1 , p 1 ( v ) , p = 1 , 2 , 3
u is along the radius of the face gear, and v is along the depth of the tooth. m and n are the number of control vertex of u and v, respectively. Set m = 21 and n = 15. ωi,j is the weight factor. Ni,3 and Nj,3 are the basis functions of u and v, respectively. di,j is the control vertex.
According to Equation (26), a tooth surface can be reconstructed, and a tooth space surface is composed of two tooth surfaces. Another tooth surface can be reconstructed by the same method. Then, the tooth space model can be obtained by Boolean operation between the tooth space and the face gear blank. Further, the 3D model of real five-pairs-of-teeth face gears can be obtained by segmenting the blank with the other tooth space surfaces, as shown in Figure 6.

5. Example and Discussion

Experiment with point cloud collection on the workpiece of face gear machined by grinding is conducted as stated in Section 2, and the parameters of a face gear are shown in Table 2. After being downsampled, the point cloud with 40,377 points is obtained. Subsequently, the point cloud is aligned using the measurement solution proposed in Section 3. In Figure 7a, the curve of objective value versus iteration number shows that the objective value of Equation (7) can reach the set accuracy at the 775th iteration. Then, the approximation of the solution at each iteration to the solution at the 775th iteration is defined by the 2-norm distance between them in the parametric space. The curve in Figure 7b shows that the solution tends to converge. Figure 8a shows that the change of the objective value of Equation (15) is extremely little after the 1000th iteration, reaching below 10−4 mm2. Figure 8b shows that the approximation between the solution at each iteration and the solution at the 1000th iteration is convergent. The convergence accuracy achieved at the 1000th iteration has almost no effect on the measurement results requiring micron-level accuracy. The above shows that 3SA has convergence characteristics.
As shown in Table 3, the objective value after the second stage of 3SA reaches 9.537 × 10−7 mm, which is a great improvement compared with 0.121 mm at the beginning, and the objective value of the tooth surface in the third stage increases from 2.23 mm2 to 0.0457 mm2.
The tooth surface reconstructed from the point cloud is compared with the theoretical tooth surface in the CATIA, as shown in Figure 9a. When viewing from top to bottom, the left side of the tooth is defined as the left tooth surface, and the right side is the right tooth surface. The result is stated as follows.
(1) Maximum deviation, not more than 0.030 mm, is on the border areas. The deviations above 0.022 mm are extremely few, and most are below 0.022 mm.
(2) The tooth surface deviations are inclined. The deviations of the inner side on the left tooth surface are negative, the deviations of the middle area trend to 0, and the deviations of the outer side are positive. It indicates that the inner side is overcut, and the outer side is undercut. While the deviations of the inner side on the right tooth surface are positive, and the deviations of the outer side are negative, which indicates that the inner side is undercut, and the outer side is overcut.
(3) The deviation trends of five teeth are almost the same.
The probable reason for the above result of deviation is that there is an angle error in the installation of the tool. Compared the measurement result of the proposed method with that of GMC shown in Figure 9b, the trend of deviation distribution is consistent. Although the maximum deviation is larger, it is distributed in the border areas that cannot be measured by GMC. The error between tooth surface deviations obtained by GMC and the proposed method is less than 7 μm, which meets the measurement requirements for most face gears. It demonstrates that the measurement solution with 3D optical scanning is valid.
Moreover, the measurement solution proposed in the paper can show deviations in the border areas and more details of deviations. For example, the deviations increase significantly in some local areas shown in Figure 9a, which may be caused by the profile error of the worm wheel.
In order to further study the effect of the deviations on the contact performance of the face gear drives in this example, the 3D model of real five-pairs-of-teeth face gear is obtained according to Section 4, and LTCA is conducted further by FEM. It should be noted that the model here contains the tooth surface deviations, and the deviations of different teeth are a little different, so the previous method cannot be used to generate the mesh. In the previous mesh generation method, the finite element meshes of a tooth model are generated first, and the finite element meshes of other teeth models are obtained by rotation. While the five teeth model contains tooth surface deviations, the finite element meshes of every tooth needs to be generated, and all finite element meshes are combined, as shown in Figure 10a. Due to the proven processing technology of the pinions, the designed model of the pinion without deviation in this example is applied, and the face gear drive is assembled, as shown in Figure 10b. A torque load of 160,000 N·mm is applied to the pinion to conduct the simulated meshing process of the real face gear drive.
The result of combining all the simulation frames is shown in Figure 11a,b and Figure 12a,b. The first and fifth teeth are the incoming and outgoing teeth of the model, respectively, so they are excluded for comparison. Figure 11 shows the contact pattern of real face gear tends to the area with the maximum tooth surface deviation and the outer side, which will cause edge contact. The maximum contact stress of real face gear is about 200 Mpa, larger than that of designed face gear. Moreover, the maximum contact stress of real face gear occurs in the border areas, which is the area with the maximum deviation meanwhile. The reason is that the “bulge” formed by the positive deviations of the tooth surface in this example will make contact earlier than the original contact points during the meshing of face gear drives. Further, it destroys the original contact path, and the “bulge” bears more contact load so that the contact stress is more concentrated in the border areas.
From Figure 12, the root bending stress of real face gear is 330.9 Mpa, while that of designed face gear is about 200 Mpa. The root bending stress increased by about 0.65 times, which will greatly affect the service life of face gear. The deviation difference between the different teeth has almost no influence on the meshing of face gear drive in this example.
The result of LTCA shows that the deviations in the border areas cause the change of contact pattern, contact stress, and root bending stress. The specific influence in this example is stated as follows.
(1) Edge contact occurs due to the deviations. The real face gear drive cannot make smooth contact as the designed model, and the tooth surfaces of the face gear and pinion collide with each other, which will cause vibration and noise. Hence, the stability of the transmission of the face gear drive will be reduced.
(2) Edge contact causes the contact stress concentrated in the border areas, and contact stress is increased by approximately 0.4, resulting in more easily occurring tooth surface pitting failure. The root bending stress is increased by approximately 0.65, resulting in more easily occurring tooth root breakage failure. The increase in contact stress and root bending stress will reduce the service life of the face gear.
Therefore, further processing should be conducted for the areas with large positive deviations, especially the border areas of the tooth surface. A simple method is to use an end mill cutter to mill the border areas of the tooth surface. The more flexible and practical method is to correct the settings of the grinding machine tool to reduce the machining errors or to modify the profile of the grinding tool to process the higher-performance tooth surface according to the measurement result of the proposed method.

6. Conclusions

In this paper, the measurement solution with 3D optical scanning is proposed to measure the tooth surfaces of orthogonal face gears. Aiming at the problem that the contact metrology method cannot measure the border areas of tooth surfaces and inefficiently acquires the huge number of points, this paper adopts the 3D optical scanning method to acquire point clouds, and the measurement solution is proposed to ensure measurement accuracy for the tooth surfaces. A 3SA is presented to align the features between point clouds and the designed face gears, including top land surfaces, centers of circles, and tooth surfaces. Moreover, the 3D modeling method is studied by reconstructing the real tooth surface with the points after being numbered. An example is given to show that the measurement solution with 3D optical scanning is valid. With the comparison, the proposed measurement solution can show the deviations in the border areas and more details that the previous measurement method could not reveal. LTCA is further conducted to study the meshing performance of the real tooth surface. The results of LTAC show that the deviations in the border areas cause edge contact, which will significantly reduce the contact strength and bending strength of real face gear. Summarily, the proposed measurement solution can provide a more reliable base for correction processing.
The proposed measurement solution can meet the measurement required for the tooth surfaces of face gears applied in automobile and other industries. In the future, the numerical algorithm of precise multi-point cloud splicing will be studied to measure face gears with all teeth, and the measurements of more precision face gears that can be applied to the aerospace industry will be conducted.

Author Contributions

Conceptualization, X.L. and Y.Z.; methodology, X.L.; software, X.L. and Z.C.; validation, X.Z. and B.H.; formal analysis, X.L.; investigation, X.L., X.Z. and B.H.; resources, X.Z., Y.Z. and J.T.; data curation, X.Z., B.H. and Y.Z.; writing—original draft preparation, X.L.; writing—review and editing, X.Z. and B.H.; visualization, X.L. and X.Z.; supervision, Y.Z. and J.T.; project administration, Y.Z. and J.T.; funding acquisition, Y.Z. and J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science and Technology Major Project of China through Grant No. 2019-VII-0017-0158, the National Natural Science Foundation of China through Grant No. 51705541, Industrial Technology Development Program No. JCKY2020213B006, Hunan Provincial Natural Science Foundation of China through Grant No. 2021JJ20071, and State Key Laboratory of High Performance Complex Manufacturing through Grant No. Kfkt2019-08 and No. ZZYJKT2020-06.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors gratefully acknowledge the support from the National Science and Technology Major Project of China, National Natural Science Foundation of China, Industrial Technology Development Program, Hunan Provincial Natural Science Foundation of China, State Key Laboratory of High Performance Complex Manufacturing, and the assistance of Zhimin Guo in the revised version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Litvin, F.L.; Fuentes, A.; Howkins, M. Design, generation and TCA of new type of asymmetric face-gear drive with modified geometry. Comput. Methods Appl. Mech. Eng. 2001, 190, 5837–5865. [Google Scholar] [CrossRef]
  2. Litvin, F.L.; Wang, J.-C.; Bossler, R.B., Jr.; Chen, Y.-J.D.; Heath, G.; Lewicki, D.G. Application of face-gear drives in helicopter transmissions. J. Mech. Des. 1994, 116, 672–676. [Google Scholar] [CrossRef]
  3. Lu, X.; Zhou, Y.; He, D.; Zheng, F.; Tang, K.; Tang, J. A novel two-variable optimization algorithm of TCA for the design of face gear drives. Mech. Mach. Theory 2022, 175, 104960. [Google Scholar] [CrossRef]
  4. Mo, S.; Song, W.; Zhu, S.; Feng, Z.; Tang, W.; Gao, H. Complex geometric modeling and tooth contact analysis of a helical face gear pair with arc-tooth. J. Cent. South Univ. 2022, 29, 1213–1225. [Google Scholar] [CrossRef]
  5. Mo, S.; Zhang, Y.; Luo, B.; Bao, H.; Cen, G.; Huang, Y. The global behavior evolution of non-orthogonal face gear-bearing transmission system. Mech. Mach. Theory 2022, 175, 104969. [Google Scholar] [CrossRef]
  6. Mo, S.; Luo, B.; Song, W.; Zhang, Y.; Cen, G.; Bao, H. Geometry design and tooth contact analysis of non-orthogonal asymmetric helical face gear drives. Mech. Mach. Theory 2022, 173, 104831. [Google Scholar] [CrossRef]
  7. Tang, J.; Cui, W.; Zhou, H.; Yin, F. Integrity of grinding face-gear with worm wheel. J. Cent. South Univ. 2016, 23, 77–85. [Google Scholar] [CrossRef]
  8. Wang, Y.; Chu, X.; Zhao, W.; Wang, Z.; Su, G.; Huang, Y. A precision generating hobbing method for face gear with assembly spherical hob. J. Cent. South Univ. 2019, 26, 2704–2716. [Google Scholar] [CrossRef]
  9. Zschippang, H.A.; Weikert, S.; Wegener, K. Face-gear drive: Simulation of shaping as manufacturing process of face-gears. Mech. Mach. Theory 2022, 172, 104791. [Google Scholar] [CrossRef]
  10. Zhou, Y.; Tang, Z.; Shi, X.; Tang, J.; Li, Z. Efficient and accurate worm grinding of spur face gears according to an advanced geometrical analysis and a closed-loop manufacturing process. J. Cent. South Univ. 2022, 29, 1–13. [Google Scholar] [CrossRef]
  11. Tang, Z.; Zhou, Y.; Wang, S.; Zhu, J.; Tang, J. An innovative geometric error compensation of the multi-axis CNC machine tools with non-rotary cutters to the accurate worm grinding of spur face gears. Mech. Mach. Theory 2022, 169, 104664. [Google Scholar] [CrossRef]
  12. Xu, H.; Zhou, Y.; He, Y.; Tang, J. An efficient approach to the five-axis flank milling of non-ferrous spiral bevel gears. Materials 2021, 14, 4848. [Google Scholar] [CrossRef] [PubMed]
  13. Litvin, F.L.; Fuentes, A. Gear Geometry and Applied Theory; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  14. Zhou, Z.; Zhang, Y.; Tang, K. Sweep scan path planning for efficient freeform surface inspection on five-axis CMM. Comput.-Aided Des. 2016, 77, 1–17. [Google Scholar] [CrossRef]
  15. Wang, S.; Zhou, Y.; Tang, J.; Tang, K.; Li, Z. Digital tooth contact analysis of face gear drives with an accurate measurement model of face gear tooth surface inspected by CMMs. Mech. Mach. Theory 2022, 167, 104498. [Google Scholar] [CrossRef]
  16. Lin, C.-H.; Fong, Z.-H. Numerical tooth contact analysis of a bevel gear set by using measured tooth geometry data. Mech. Mach. Theory 2015, 84, 1–24. [Google Scholar] [CrossRef]
  17. Ding, H.; Wan, Z.; Zhou, Y.; Tang, J. A data-driven programming of the human-computer interactions for modeling a collaborative manufacturing system of hypoid gears by considering both geometric and physical performances. Robot. Comput.-Integr. Manuf. 2018, 51, 121–138. [Google Scholar] [CrossRef]
  18. Stadtfeld, H.J. Handbook of Bevel and Hypoid Gears: Calculation, Manufacturing, Optimization; Rochester Institute of Technology: Rochester, NY, USA, 1993. [Google Scholar]
  19. Rong, K.; Ding, H.; Chen, S.; Wang, Y.; Zhou, Z.; Lu, R. Top-Rem grinding tool modification considering loaded edge contact for spiral bevel gears. Adv. Eng. Inform. 2022, 53, 101697. [Google Scholar] [CrossRef]
  20. Kuş, A. Implementation of 3D optical scanning technology for automotive applications. Sensors 2009, 9, 1967–1979. [Google Scholar] [CrossRef]
  21. Hawryluk, M.; Ziemba, J. Application of the 3D reverse scanning method in the analysis of tool wear and forging defects. Measurement 2018, 128, 204–213. [Google Scholar] [CrossRef]
  22. Valigi, M.C.; Logozzo, S.; Affatato, S. New challenges in tribology: Wear assessment using 3D optical scanners. Materials 2017, 10, 548. [Google Scholar] [CrossRef]
  23. Chen, Y.-C.; Chen, J.-Y. Optical inspection system for gear tooth surfaces using a projection moiré method. Sensors 2019, 19, 1450. [Google Scholar] [CrossRef]
  24. Gonzalez-Perez, I.; Fuentes-Aznar, A. Tooth contact analysis of cylindrical gears reconstructed from point clouds. In New Approaches to Gear Design and Production; Goldfarb, V., Trubachev, E., Barmina, N., Eds.; Springer: Cham, Switzerland, 2020; pp. 219–237. [Google Scholar]
  25. Urbas, U.; Zorko, D.; Černe, B.; Tavčar, J.; Vukašinović, N. A method for enhanced polymer spur gear inspection based on 3D optical metrology. Measurement 2021, 169, 108584. [Google Scholar] [CrossRef]
  26. Urbas, U.; Hrga, T.; Povh, J.; Vukašinović, N. Novel alignment method for optical 3D gear metrology of spur gears with a plain borehole. Measurement 2022, 192, 110839. [Google Scholar] [CrossRef]
  27. Brajlih, T.; Valentan, B.; Drstevenser, I.; Balic, J.; Pogacar, V. Testing the accuracy of ATOS[TM] 3d optical scanner measuring volumes. Ann. DAAAM Proc. 2007, 2007, 111–113. [Google Scholar]
  28. Hwang, J.; Yu, S.; Kim, J.; Lee, S. 3D face modeling using the multi-deformable method. Sensors 2012, 12, 12870–12889. [Google Scholar] [CrossRef]
  29. Vagovský, J.; Buranský, I.; Görög, A. Evaluation of measuring capability of the optical 3D scanner. Procedia Eng. 2015, 100, 1198–1206. [Google Scholar] [CrossRef]
  30. Hawryluk, M.; Ziemba, J.; Zwierzchowski, M.; Janik, M. Analysis of a forging die wear by 3D reverse scanning combined with SEM and hardness tests. Wear 2021, 476, 203749. [Google Scholar] [CrossRef]
  31. Zhong, K.; Li, Z.; Li, R.; Shi, Y.; Wang, C. Pre-calibration-free 3D shape measurement method based on fringe projection. Opt. Express 2016, 24, 14196–14207. [Google Scholar] [CrossRef]
  32. Huang, J.; Wu, Q. A new reconstruction method based on fringe projection of three-dimensional measuring system. Opt. Lasers Eng. 2014, 52, 115–122. [Google Scholar] [CrossRef]
  33. He, W.; Zhong, K.; Li, Z.; Meng, X.; Cheng, X.; Liu, X.; Shi, Y. Accurate calibration method for blade 3D shape metrology system integrated by fringe projection profilometry and conoscopic holography. Opt. Lasers Eng. 2018, 110, 253–261. [Google Scholar] [CrossRef]
  34. Valagiannopoulos, C.A. On smoothening the singular field developed in the vicinity of metallic edges. Int. J. Appl. Electromagn. Mech. 2009, 31, 67–77. [Google Scholar] [CrossRef]
  35. Young, I.T.; van Vliet, L.J. Recursive implementation of the Gaussian filter. Signal Process. 1995, 44, 139–151. [Google Scholar] [CrossRef]
  36. Al-Rawabdeh, A.; He, F.; Habib, A. Automated feature-based down-sampling approaches for fine registration of irregular point clouds. Remote Sens. 2020, 12, 1224. [Google Scholar] [CrossRef]
  37. Rizk-Allah, R.M.; Hassanien, A.E. A hybrid equilibrium algorithm and pattern search technique for wind farm layout optimization problem. ISA Trans. 2022, in press. [CrossRef]
  38. Melo, I.D.; Antunes, M.P. Bad data correction in harmonic state estimation for power distribution systems: An approach based on generalised pattern search algorithm. Electr. Power Syst. Res. 2022, 204, 107684. [Google Scholar] [CrossRef]
  39. Zhou, Y.; Wu, Y.; Wang, L.; Tang, J.; Ouyang, H. A new closed-form calculation of envelope surface for modeling face gears. Mech. Mach. Theory 2019, 137, 211–226. [Google Scholar] [CrossRef]
  40. Barone, S.; Borgianni, L.; Forte, P. Evaluation of the effect of misalignment and profile modification in face gear drive by a finite element meshing simulation. J. Mech. Des. 2004, 126, 916–924. [Google Scholar] [CrossRef]
  41. Pietrusiak, D.; Wróbel, J.; Czechowski, M.; Fiebig, W. Dynamic NVH numerical analysis of power steering in the presence of lubricant in the system. Materials 2022, 15, 2406. [Google Scholar] [CrossRef]
  42. Fu, X.; Wang, B.; Wang, T. Experimental and numerical investigations on the deformation defects by gear rolling process with local induction heating. Int. J. Adv. Manuf. Technol. 2021, 117, 835–844. [Google Scholar] [CrossRef]
  43. Ge, Y.; Liu, Z.; Zhou, K.; Yang, J.; Liu, D. Investigation of the electromagnetic characteristics and operating performance of a bidirectional PM excited machine. Sci. Rep. 2022, 12, 8408. [Google Scholar] [CrossRef]
  44. Piegl, L.; Tiller, W. The NURBS Book; Springer: Berlin/Heidelberg, Germany, 1996. [Google Scholar]
  45. Oesterle, B.; Geiger, F.; Forster, D.; Fröhlich, M.; Bischoff, M. A study on the approximation power of NURBS and the significance of exact geometry in isogeometric pre-buckling analyses of shells. Comput. Methods Appl. Mech. Eng. 2022, 397, 115144. [Google Scholar] [CrossRef]
  46. Pham-Tan, H.; Thai, C.H.; Phung-Van, P. NURBS-based refined plate theory for metal foam plates with porosities. Thin-Walled Struct. 2022, 175, 109246. [Google Scholar] [CrossRef]
Figure 1. Fringe projection.
Figure 1. Fringe projection.
Materials 15 06069 g001
Figure 2. (a) ATOS 3D Scanner, (b) a face gear, and (c) the point cloud collection.
Figure 2. (a) ATOS 3D Scanner, (b) a face gear, and (c) the point cloud collection.
Materials 15 06069 g002
Figure 3. (a) Unprocessed point cloud, (b) filtering, and (c) down-sampling.
Figure 3. (a) Unprocessed point cloud, (b) filtering, and (c) down-sampling.
Materials 15 06069 g003
Figure 4. The three-stage alignment for top land surface, center of circles, and tooth surface.
Figure 4. The three-stage alignment for top land surface, center of circles, and tooth surface.
Materials 15 06069 g004
Figure 5. (a) The regularly discrete points on the theoretical tooth surface, and (b) the point cloud with numbered points.
Figure 5. (a) The regularly discrete points on the theoretical tooth surface, and (b) the point cloud with numbered points.
Materials 15 06069 g005
Figure 6. A 3D model of real five-pairs-of-teeth face gear.
Figure 6. A 3D model of real five-pairs-of-teeth face gear.
Materials 15 06069 g006
Figure 7. Convergence curves of (a) objective value and (b) iterative solution in the second stage of the algorithm.
Figure 7. Convergence curves of (a) objective value and (b) iterative solution in the second stage of the algorithm.
Materials 15 06069 g007
Figure 8. Convergence curves of (a) objective value and (b) iterative solution in the third stage of the algorithm.
Figure 8. Convergence curves of (a) objective value and (b) iterative solution in the third stage of the algorithm.
Materials 15 06069 g008
Figure 9. (a) Holistic deviation of the tooth surfaces obtained by the proposed method and (b) the measurement report of GMC.
Figure 9. (a) Holistic deviation of the tooth surfaces obtained by the proposed method and (b) the measurement report of GMC.
Materials 15 06069 g009
Figure 10. FEM: (a) Mesh generation, (b) the finite element model of real five-pairs-of-teeth face gear drive.
Figure 10. FEM: (a) Mesh generation, (b) the finite element model of real five-pairs-of-teeth face gear drive.
Materials 15 06069 g010
Figure 11. The simulation results for the contact pattern and contact stress of (a) the real face gear drive and (b) the designed face gear drive.
Figure 11. The simulation results for the contact pattern and contact stress of (a) the real face gear drive and (b) the designed face gear drive.
Materials 15 06069 g011aMaterials 15 06069 g011b
Figure 12. The simulation results for the bending stress of (a) the real face gear drive and (b) the designed face gear drive.
Figure 12. The simulation results for the bending stress of (a) the real face gear drive and (b) the designed face gear drive.
Materials 15 06069 g012
Table 1. The configuration of ATOS 3D Scanner.
Table 1. The configuration of ATOS 3D Scanner.
SystemATOS 3D Scanner
Camera pixels 2   ×   3296   ×   2472 pixels
Measuring volumes 38   ×   29   ×   15     2000   ×   1500   ×   1500   mm 3
Point spacing0.01–0.61 mm
Working distance490–2000 mm
Projected light sourceStructured blue light
Operating Temperature5–40 °C
Table 2. The parameters of a face gear.
Table 2. The parameters of a face gear.
ParametersValues
Shaper teeth number Ns22
Pinion teeth number N119
Face gear teeth number N2142
Module m1.95 mm
Pressure angle of the rack cutter α25°
Shaft angle γm90°
Inner radius L1128 mm
Outer radius L2152 mm
Table 3. The comparison of objective value before and after computing in 3SA.
Table 3. The comparison of objective value before and after computing in 3SA.
Second StageThird Stage
Objective valueBefore0.868 mm2.23 mm2
After9.537 × 10−7 mm0.0043 mm2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lu, X.; Zhao, X.; Hu, B.; Zhou, Y.; Cao, Z.; Tang, J. A Measurement Solution of Face Gears with 3D Optical Scanning. Materials 2022, 15, 6069. https://doi.org/10.3390/ma15176069

AMA Style

Lu X, Zhao X, Hu B, Zhou Y, Cao Z, Tang J. A Measurement Solution of Face Gears with 3D Optical Scanning. Materials. 2022; 15(17):6069. https://doi.org/10.3390/ma15176069

Chicago/Turabian Style

Lu, Xinxin, Xing Zhao, Bo Hu, Yuansheng Zhou, Zhezhen Cao, and Jinyuan Tang. 2022. "A Measurement Solution of Face Gears with 3D Optical Scanning" Materials 15, no. 17: 6069. https://doi.org/10.3390/ma15176069

APA Style

Lu, X., Zhao, X., Hu, B., Zhou, Y., Cao, Z., & Tang, J. (2022). A Measurement Solution of Face Gears with 3D Optical Scanning. Materials, 15(17), 6069. https://doi.org/10.3390/ma15176069

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop