#### Abstract

This article introduces new types of rational approximations of the inverse involute function, widely used in gear engineering, allowing the processing of this function with a very low error. This approximated function is appropriate for engineering applications, with a much reduced number of operations than previous formulae in the existing literature, and a very efficient computation. The proposed expressions avoid the use of iterative methods. The theoretical foundations of the approximation theory of rational functions, the Chebyshev and Jacobi polynomials that allow these approximations to be obtained, are presented in this work, and an adaptation of the Remez algorithm is also provided, which gets a null error at the origin. This way, approximations in ranges or degrees different from those presented here can be obtained. A rational approximation of the direct involute function is computed, which avoids the computation of the tangent function. Finally, the direct polar equation of the circle involute curve is approximated with some application examples.

#### 1. Introduction

The circle involute function, used to mathematically describe the flank of the gear tooth, is formulated as *x* = inv(*u*) = tan(*u*) − *u* (see Figure 1), where *u* is the pressure angle and *x* is the polar angle. On the other hand, the inverse involute function that computes the angle from the involute is used, for example, in the computation of the center distances between the gear axes, and cannot be expressed mathematically in a direct way. So, either tables or the best approximate formulae obtained to date should be used (see Section 2). The shape of the inverse involute function, until 70° (1.5 rad approx.), for example (see Figure 2), shows an elbow in its initial part, starting at the origin with a very high value in its derivative and then changing until almost horizontal. It is known that rational functions are especially appropriate to approximate functions with that elbow shape (see [1]). The authors of this article worked previously with the applications of the approximation theory by rational functions applied to tyre models, with elbow-shaped curves similar to the inverse involute, and obtained very efficient and accurate approximation results (see [2–6]). This article first presents information from previous works and then develops the theoretical foundation and algorithms used to obtain the approximate expressions of both the inverse involute and direct functions proposed at the end.

#### 2. Previous Works

Several authors have been publishing formulations that compute the inverse involute function. In the following sections, we will call *x* the involute of an angle *u*(rad): *x* = inv(*u*); *u* = inv − 1(*x*). To compare the errors, we will take as a reference a maximum value of the angle *u* = 45° = 0.785398 rad, which has an involute value of *x* = 0.2146. This value fulfills a major part of the engineering application requirements.

Initially, we must cite the work of Dudley [7], who proposed the following expression:which is valid for *x* ≤ 0.5 with a maximum error for *u* = 45° and *E*_{max} = 0.0012 rad = 4.12529′. This error has been calculated when comparing it with equation (3). Due to its error, this expression is only valid as a first approximation.

The iterative iteration by Laskin is far more accurate, which is as follows [8]:

This iterative formulation is valid for angles up to 65°. Its error grows with *x*, and it takes values for *u* = 45° and *E*_{max} = 0.000017 rad = 3.506″ with two iterations (computation until *u*_{2}). This error has been calculated when comparing it with equation (3). Its error is reduced when the number of iterations increases.

On the other hand, Cheng [9] used the perturbation theory, implementing the method of dominant balance in order to obtain approximations to the inverse involute function, proposing the following expression:

The previous expression in these 9 terms has a very low error, which takes values for *u* = 45° and *E*_{max} = 1.58*E* − 9 rad = 3.2589*E* − 4″ (see [10]). For its computation, it requires 38 basic floating-point operations (FLOPs) plus a cube root.

In his article [9], Cheng economized (see Section 3.1) the previous truncated series while taking into account only the first 4 terms (until *x*^{7/3}) when performing the transformation *t* = *x*^{2/3}, so the resulting series is as follows:

As the term “*t*^{1/2}” multiplies the brackets, the result will pass through the origin. Then, the 3-degree series within the brackets is economized until 1 degree using the Chebyshev polynomials, undoing the change afterwards, *x* = *t*^{3/2}, so that the following equation is obtained:

Very similar to calculation by Dudley, valid for 0 ≤ *u* ≤ 45°, with a maximum error *E*_{max} = −0.001308 rad = 4.49979′. When the angle is 45°, the value of (Inv(*x*))^{1/3} is 0.5987026.

In 2017, more recent works have been published by Liu [11], using a method similar to the one used by Cheng, but applying the integral expression of the arctan functionwhere *c* is an intermediate variable (*c* = tan(*u*)). Then, Liu expanded the tangent function in Taylor Series, obtaining finally the following two expressions: the first one, more simple but less accurate, is similar to Cheng’s (3):with an error (at *u* = 45^{o}) *E*_{max} = 2.75*E* − 5 rad = 1.73*E* − 3ʺ. The error continues growing for bigger values of *u*. The second expression, which is very accurate, but not too simple to compute, proposed in Liuʼs work iswith an error (at *u* = 45°) *E*_{max} = 1.73*E* − 11 rad = 3.57*E* − 6″ (see [11] for more details).

In this article, we start from Cheng’s expression (3) but we obtain approximations truncating expansions in a series of Jacobi polynomials and a rational series; in order to get that, we start from the Chebyshev–Padé rational approximations (see [12]) and the original formulation of the Remez algorithm, which is modified thus imposing the condition of a null error at the origin. The resulting expressions are very compact and quickly executed, with a very low error. To justify the resulting final expressions and compute some others, the concepts of expansion in a series of the Chebyshev polynomials, economization, and approximation of functions are briefly reviewed, and the series of Jacobi polynomials allowing for the selection of areas with a lower error in a flexible way are also described. The theory of rational functions, the Remez algorithm, and its adaptation proposed in this article are reviewed. Finally, the proposed formulae are presented.

#### 3. Theoretical Basis of the Approximations Used in This Article

##### 3.1. Approximation of a Function in the Chebyshev Series

The Chebyshev polynomials (see [10]) of the first kind are defined by , and they are orthogonal regarding the function in the interval . For a detailed explanation of the concepts of orthogonality, orthogonal polynomials, and orthogonality intervals, see [10]. To work in different [*a*, *b*] intervals, shifted polynomials with the following change must be used:

Successive terms of the Chebyshev polynomials are as follows:

In [10], we can find formulae to calculate products, integrals, and derivatives of the Chebyshev polynomials. The expansion of a function in the Chebyshev series has the following form:in which the prime symbol in the summation indicates that the first term must be divided by 2, beingwhere is the weight function:

As *N* grows, the series summation approaches to ; that is, it converges to . If we truncate the series in degree *N*, we get an approximation to the function, and it is more accurate with a higher *N*. Due to the properties of the Chebyshev polynomials, truncating in *N* − 1 is the best *N* − 1 degree polynomial approximation to the development at *N* degree, always working in the orthogonality interval [−1, 1]; this is called economization (see Burden and Faires [13]) and is the method used by Cheng to obtain equation (5) from equation (4). *r*_{n} is the norm of the function (*π*/2 for the Chebyshev polynomials).

In the following sections, Jacobi polynomials are introduced because they allow a flexible error choice in some areas of the curve, which will be useful. Then, we will present approximations in rational functions because they are more accurate and easier to compute.

##### 3.2. Expansion in a Series of Jacobi Polynomials

Jacobi polynomials come from the group of classic orthogonal polynomials obtained from the Sturm–Liouville differential equation, from which the Chebyshev polynomials are also obtained (see [10]). The expansion of a function in a series of Jacobi polynomials uses a Jacobi weight function this time:

The weight function in Jacobi orthogonal polynomials is as follows:

This function is controlled by two parameters, *δ* and *γ*, that allow choosing the area of a best approximation at the orthogonality interval. The minimum value of both *δ* and *γ* is −1. If *δ* = *γ* = −0.5, then we obtain the weight function of the Chebyshev polynomials (13). If both are equal, *δ* = *γ*, the error is symmetric regarding the *y*-axis. If *δ* = *γ* = −1, then we get a null error at both ends of the interval. If *δ* = −1, then the error is null at the right end. If *γ* = −1, then the error is null at the left side of the interval.

##### 3.3. Approximation in a Jacobi Series of the Inverse Involute Function

We can expand the complete *f*(*x*) polynomial obtained by Cheng (see (3)) in a Jacobi series with *δ* = *γ* = −1, as it is shown in (14), which gives a near-zero error at both ends; however, we have to previously convert equation (3) to the *y* domain by means of the transformation *y* = *x*^{1/3}, obtaining a degree 17 polynomial in *y*. Then, we convert this polynomial from the *y* domain to the domain [−1, 1] using (9). After that, integral (14) must be used to calculate the coefficients of the Jacobi series in terms of . Then, we can truncate the calculation at the desired degree. Finally, we have to undo the previous transformations to go back to the *y* domain and then to the *x* domain to obtain the following approximation:

It requires 4 products and 2 sums and has an error distributed along the interval with a maximum value *E*_{max} = 0.000222 rad = 0.7631′ for 45°, much smaller than expression (5). Figure 3 shows error curves of both equations (5) and (16).

##### 3.4. Approximation Using Rational Functions (RA)

Approximations in rational functions (see [13]) with polynomials of degree *n* at the numerator and degree *m* at the denominator (we will denote them as [*n*, *m*]) are more efficient when the function varies in a quick way in some areas but not in some others [1], as happens in the inverse involute function. As an example, in [14], we can find a nice application of rational functions to the vibrational frequency of a one-dimensional bar. In [13], we can see the theory of Pade and Chebyshev–Pade rational approximations. In MAPLE, the Chebpade function implements the rational approximations and it allows the numerical conversion between normal polynomials and Chebyshev polynomials too. The Chebyshev–Padé functions obtain good approximations, but not those of the minimum value of the maximum error (known as minimax). To find the latter, the iterative Remez algorithm is used [10, 12, 15–17], which in turn stems from the Chebyshev–Padé approximation; it fine tunes the result by numeric iterations and converges to an improved minimax approximation of both rational and polynomial functions.

According to Tawfik [17], supposing an interval [*a*, *b*], Chebyshev’s criterion proposes that if *P*_{n}(*x*) is the minimax polynomial of an *n* degree, there must be at least (*n* + 2) points in this interval where the error function *E*(*x*) obtains the maximum absolute value *E*_{max} with an alternating sign. Chebyshev demonstrated the existence and unicity of that minimax polynomial *P*_{n}(*x*) (or rational function *R*_{n}(*x*_{i})). This criterion is mathematically defined with the following equations:

In this manner, we can obtain a system of (*n* + 2) equations with (*n* + 2) unknown (the coefficients of *P*_{n}(*x*), or *R*(*x*), and *E*_{max}). The system will be nonlinear in the case of a rational approximant *R*_{n}(*x*). Solving it, we can get a new set of points *x*_{i} with alternating the maximum error *E*_{max}, which will be more reduced in this situation. Repeating this iteration several times, we can reach the Chebyshev minimax condition so that the value of the error function in the final set of (*n* + 2) points (*E*_{max}) represents the absolute maximum value of the approximate error (see [17] for more details).

Sometimes imposed constraints are needed for the approximating function, for example, when the original function takes a null value at the origin and it seems conceptually unacceptable to take a nonzero value at the origin in the approximated function. In this paper, the previously developed Remez algorithm has been adapted to obtain rational approximations with null error at the origin just adding that mathematical constraint to the general algorithm presented in [17]:

To do this, the iterations can start with the solution achieved through Jacobi polynomials in the case *δ* = −0.5 and *γ* = −1, which gives a null error at the origin; however, minimax is not taken as the local maxima of the error function, which take different values. In summary,(1)We started from the Cheng 17-degree polynomial (see equation (3)).(2)We convert it into a rational function [1, 2] using the Chebyshev–Padé method (see [13]). Previously, the result has been transformed to the *y* domain.(3)Then, we apply the adapted Remez algorithm, starting the iteration with the resulting roots in the previous step, imposing the constraint of null error at the origin.

We will apply this 3-step method to the inverse involute function in Section 4 to obtain Apsol8.

To either manipulate the polynomials or approximate them to rational functions and vice versa, the direct and inverse Chebyshev–Padé method can be used (see [13]).

Applications of expansions in a Chebyshev and Jacobi series of rational functions to special functions, heating and air conditioning, vehicle dynamics, steering mechanisms, and tyre nonlinear behavior can be found in [18].

#### 4. Rational Approximation to the Inverse Involute Function by the Adapted Remez Algorithm

Once we have the required theoretical background, we apply the method described in Section 3.4 to obtain Remez minimax approximations passing through the origin for the inverse involute function from the Cheng polynomial (3) and obtain the following:

We could see graphically regarding the 17-degree polynomial that this solution has a maximum error of 8.34*E* − 6 rad (1.72″) and the error at the origin is not null.

Applying the Remez algorithm with the null error at the origin constraint, the following can be obtained:

In Figures 4 and 5, we can see the errors of the approximations Apsol 1 and Apsol2. Apsol 1 is obtained using the Chebyshev–Padé algorithm, it does not pass at the origin as shown more clearly in Figure 5. This process has a maximum error of *E*_{max} = 7.5*E* − 6 rad (1.54″) and null error at the origin, and it requires two operations less than the previous process (Apsol1) to include, 5 FLOPs in total plus the cube root.

If we want to reduce the number of decimals, admitting bigger errors, the approximation can be

Graphically, we could see that a slightly bigger maximum error (*E*_{max} = 8.847*E* − 6 rad = 1.82″) occurs but also passes through the origin and requires 5 FLOPs plus the cube root.

As the value of the standard pressure angle is 20°, we will frequently need approximations in its environment; therefore, if we consider the interval between 0° and 30°, approximations will be either more accurate or simple. If we repeat the previous process, we can obtainwhich passes through the origin and has a maximum error *E*_{max} = 1.124*E* − 7 rad = 0.0231″ in the interval (0°–30°).

In the previous Apsol4, the term *x*^{1/3} is negative and is almost null, but for the range 0°–45^{o} (Apsol2), the term is positive, indicating that it achieves a value very close to 0 at some intermediate point between 30° and 45°. This error is observed to be almost zero in the interval 0°–35°, so it can be rounded with a very low error, obtaining a very compact and very accurate final expression:

This also passes through the origin and has a very low maximum error *E*_{max} = 6.6*E* − 7 rad = 0.137″ in the interval (0°–35°), requiring only 3 FLOPs plus the cube root.

The limit of 35° reaches the head radius with *z* ≤ 14. The limit of 45° reaches the head radius with *z* ≤ 7, in which *z* is the number of teeth.

Table 1 shows a summary of the approximations to the inverse involute function presented in this article.

The three new formulae have been designed to work within the intervals shown at column 2: 0°–45° for the first two formulae and 0°–35° for the latter. The values of maximum errors can be found within the intervals. If they are used outside these intervals, the errors grow very fast as we can see in Figures 3 and 4, so they cannot be used at all.

The interval 0°–35° is wide enough for many engineering applications, if we work in that interval, the third expression is very compact, easy to use, and more accurate than the others.

Liu’s rational inverse trigonometric function (8) and Liu’s 6 terms (7) are more accurate, but are far more complicated. Dudley’s, Laskin’s, and Cheng’s polynomials (2 terms) are less accurate. For these reasons, we recommend the third new formula Apsol5.

#### 5. Rational Approximation of the Direct Involute Function

The approximation of the inverse involute using polynomial (3) and the rational approximations obtained from equation (3) have the drawback of computing of the cube root of *x*, which is not very computationally efficient. However, the development in series of that same inverse involute function in terms of *x*, which should be more efficient, requires polynomials of a higher degree or in the case of rational functions higher degrees at the numerator and the denominator.

The exact computation of the involute function can be directly performed, inv(*u*) = tan(*u*) − *u*, which is why its polynomial or rational approximation is less interesting. Nevertheless, as the computation of the function tan(*u*) is not very efficient and in order to complete this work, this approximation is presented as another example of the application of the methods described here. Rational expressions in terms of *u* can be obtained with a reasonable approximation for the direct involute function, with lower degree polynomials than in the inverse involute, as shown below with a [2, 2] degree example:

It is graphically shown that the maximum error up to 45° is as follows: *E*_{max} = 0.0003209, 0.15% of the full-scale value (0.2146). It requires 8 FLOPs but avoids the calculation of the tangent. Obviously, by inverting Apsol5, we could obtain another approximation to the direct involute as follows: , which is valid up to 35° with an ultralow maximum error in this interval: *E*_{max} = 8.37*E* − 11, 9.36*E* − 8% of the full-scale value (0.089342).

If we use Apsol3, we can obtain the following:which is valid up to 45° with a maximum error *E*_{max} = 6.35*E* − 8, 2.95*E* − 5% of the full-scale value (0.2146).

#### 6. Applications

The proposed expressions can simplify many calculations in gear engineering. We show here some examples:

##### 6.1. The Direct Polar Equation of the Involute Curve

As we know from Figure 1, the polar equation of the circle involute curve at a given radius *r*_{x} is *x* = Inv(*u*) = tan(*u*) − *u*; *u* = Inv^{−1}(*x*). Assuming that we use the standard pressure angle of 20°, if *r*_{x} is the pitch radius (*r* = *r*_{x}), then *u* = 20° and

We can use the previous techniques to obtain a direct relation between the radius *r*_{x} and the polar angle *x*. In this case, the approximate function does not take the value 0 at the origin and a zero error does not need to be forced at the origin. We can find the following minimax rational [2, 2] expression, assuming a standard pressure angle of 20°:

This result is valid up to *u* ≤ 35°. It has a low maximum error. In a function of the pitch radius, we will have approximately *E*_{max} = (*r*/25)*E* − 4 (obtained by calculating the error for different values of *r*), both the pitch radius and the error being expressed in mm.

Of course, we could directly use the Apsol5, which would be much more accurate, assuming the following computation of cosine functions:

The result is valid up to *u* *≤* 35° with a maximum error *E*_{max} = (*r*/25)*E* − 6, in which both the pitch radius and the error are expressed in mm. In the range up to *u* *≤* 45°, we can use Apsol3 to obtain a maximum error *E*_{max} = (*r*/10)*E* − 5, with the pitch radius and the error expressed in mm.

Similar expressions could be obtained for different standard pressure angles (25°, for example).

##### 6.2. Center Distance When Using Shifted Gears with Zero Circular Tolerance

In this case, it is known that the functioning pitch radii arein which is the operating pressure angle that gives 0 circular tolerance and good transmission, when using toothed wheels with *z*_{1} and *z*_{2} teeth and *s*_{1} and *s*_{2} as the profile shifts. From the calculation of , we can use Apsol8 to obtain the pitch radii and the center distance.

The polar equation of the circle involute curve can be used to simplify calculations in many applications, as mentioned in previous research [19], where the authors calculate the load distribution along the line of contact for involute external gears and the elastic potential energy of a spur tooth, or the approximate method shown in another study [20].

#### 7. Conclusions

In this article, new approximations to the inverse involute function have been presented. They are very accurate for a low degree polynomial, requiring a much reduced number of operations. In general, the balance between the computation efficiency and accuracy of the resulting rational approximations is very favorable for most engineering applications, including the following example, which is valid for *u* ≤ 35°:

Although some specific approximations with a limited value of *u* = 45° and rational functions of [1, 2] and [2, 2] degrees have been presented here, we have also described the theoretical principles and the algorithms for applying them, which allow other approximations to result in larger or smaller degrees or on larger or smaller angle values, even in intervals that do not start from zero, which would generate much lower errors when reducing the approximation interval. This has been completed in this article with constraint imposition applied to the Remez algorithm. Additionally, it provides a rational approximation to the direct involute function as a new example of application, which avoids computation of the tangent function or the cube root in an eventual expansion in a series. Finally, the polar equation of the involute curve has been approximated with a very low error. It will allow a simplification of many other calculations in gear engineering. Some examples of its application have been presented.

#### Nomenclature

U: | Pressure angle in the involute curve. As we want to obtain approximations of u, we will consider it a function of the independent variable x; u = inv^{−1}(x) |

: | Operating pressure angle used at point 6 |

c: | c = tan(u) is the intermediate variable used in (6) |

u_{1}, u_{2}, u_{3}: | Intermediate values of the pressure angle u in the iterations of the Laskin formulation |

x: | Polar angle, as the involute of the pressure angle u. We use it too as an independent variable working in the interval [a, b] |

t: | Intermediate variable resulting of the transformation t = x^{2/3} used in (4) |

y: | Intermediate variable resulting of the transformation y = x^{1/3} used at Section 3.3 |

: | Variable resulting of transformation (9) working in the interval [−1, 1] |

a, b: | Limits of the interval of values of x in (9) |

T_{i}(x): | Chebyshev polynomial of degree i for the variable x |

J_{i}(x): | Jacobi polynomial of degree i for the variable x |

a_{n}: | Coefficients of a Chebyshev series |

r_{n}: | Norm of the function a_{n} in (12) |

: | Weight function in the calculation of the coefficients in orthogonal polynomial series |

f(x): | A general function in terms of the independent variable x |

δ, γ: | Parameters in the weight function of Jacobi polynomials |

E_{max}: | Maximum error of an approximation |

z_{i}: | Number of teeth in the gear i |

r: | Pitch radius of a gear |

: | Functioning pitch radius of a gear |

s_{i}: | Profile shift of the gear i. |

#### Data Availability

The data used to support the findings of this study have been deposited in the library repository of the Universidad Francisco de Vitoria at http://ddfv.ufv.es/handle/10641/1435.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Supplementary Materials

The supplementary material is an Excel file including all the calculations of the different mathematical expressions of the inverse involute present in the paper.* (Supplementary Materials)*