# (UPDATED) How do I find the intersection of this involute and trochoid, given their parametric equations?

I have two parametric curves. The first, an involute of a circle with parameter theta: $$x_1(θ) = Z × \cos(α) × (\cos(θ − \tan(α) + α) + θ × \sin(θ − \tan(α) + α)),$$ $$y_1(θ) = Z × \cos(α) × (\sin(θ − \tan(α) + α) − θ × \cos(θ − \tan(α) + α)),$$ $$0 ≤ θ_{min} ≤ θ ≤ θ_{max}$$ (note that theta must always be greater than or equal to zero)

And the second, a trochoid with parameter gamma: $$x_2(γ) = (Z − 2 × (C − X)) × \cos(γ) + (2 × \tan(α) × (C − X) + Z × γ) × \sin(γ),$$ $$y_2(γ) = (Z − 2 × (C − X)) × \sin(γ) − (2 × \tan(α) × (C − X) + Z × γ) × \cos(γ),$$ $$γ_{min} ≤ γ ≤ γ_{max}$$ With four independent variables:

- $Z$ is a positive integer,
- $X$ is a real number in the interval $[−1, 1]$,
- $C$ is a real number in the interval $[1, 1.5]$, and
- $α$ is an angle in the interval $\left[0, \arctan\left(\frac{π}{4 × C}\right)\right]$.

- $θ_{max} = \frac{\sec(α)}{Z} × \sqrt{(2 × X + Z + 2)² − Z² × \cos(α)²}$,
- $γ_{max} = \frac{−2}{Z} × (C − X) × \tan(α)$, and
- $θ_{min}$ and $γ_{min}$ are the parameter values at the intersection that I want to find.

I know that the curves have a

*touching*intersection (where their tangents are the same) at the parameter values $θ = \tan(α) − \frac{4}{Z} × (C − X) × \csc(2 × α)$ and $γ = \frac{−4}{Z} × (C − X) × \csc(2 × α)$. See the red and blue curves in this picture: https://i.stack.imgur.com/XFcNp.png

However, when the values of the four variables are such that $α < \arcsin\left(\sqrt{\frac{2}{Z} × (C − X)}\right)$, two things happen. First, the expression for theta in the previous paragraph goes negative, which makes it invalid for my purposes. Second, a

*transversal*intersection (where the curves' tangents are distinct) appears at values of gamma and theta closer to zero than those for the touching intersection, and specifically with a positive value of theta. See the red and blue curves in this picture: https://i.stack.imgur.com/5fM6k.png

**THE PROBLEM:**

**I want to find the parameter values $θ_{min}$ and $γ_{min}$ for this**When that inequality is

*transversal intersection*in terms of $Z$, $α$, $X$, and $C$, in all cases where $α < \arcsin\left(\sqrt{\frac{2}{Z} × (C − X)}\right)$, given the stated domains of the independent variables $Z$, $α$, $X$, and $C$. Additionally, when that inequality is true, $θ_{min}$ should always be strictly greater than zero.*not*true, the value of any new expression for $θ_{min}$ can be anything, because I already have expressions that hold when the opposite of that inequality, $α ≥ \arcsin\left(\sqrt{\frac{2}{Z} × (C − X)}\right)$, is true.

I do know that on the first curve, the radius of a point for a given value of theta is $r(θ) = Z × \cos(α) × \sqrt{θ² + 1}$, and conversely the value of theta for a given radius is $θ(r) = \sqrt{\frac{r²} {Z² × \cos(α)²} − 1}$. This means that if I can find the radius of the transversal intersection point by any process, I can easily convert it into the value for $θ_{min}$, and vice versa. I also know that the angle made with the $x$-axis by a point on the curve at a given value of theta is $β(θ) = θ − \arctan(θ) + α − \tan(α)$, which is transcendental and has no closed-form inverse, so I

*cannot*use a known angle to find the value of theta.

Similarly, I know that for the second curve, the radius of a point for a given value of gamma is $r(γ) = \sqrt{(2 × \tan(α) × (C − X) + Z × γ)² + (Z −2 ×(C − X))²}$, and conversely the value of gamma for a given radius is $γ(r) = \frac{±\sqrt{r² − (Z − 2 ×(C − X))²} − 2 ×(C − X) ×\tan(α)}{Z}$. This means that if I can find the radius of the transversal intersection point by any process, I can convert it into the value for $γ_{min}$, and vice versa. I also know that the angle made with the $x$-axis by a point on the curve at a given value of gamma is $β(γ) = γ − \arctan\left(\frac{Z (γ + \frac{2}{3} ×(C − X) ×\tan(α))}{Z −2 ×(C − X)}\right)$, which is transcendental and has no closed-form inverse, so I

*cannot*use a known angle to find the value of gamma.

Thus, if I get either one of $θ_{min}$ or $γ_{min}$, I can use that value to find the other. If I can find the radius of the intersection separately, I can use it to find both values.

Given these expressions, I know that the radius of the touching intersection is $r = \sqrt{4 ×\cot(α)² ×(C − X)² + (Z − 2 ×(C − X))²}$, or equivalently $r = Z ×\cos(α) ×\sqrt{(\tan(α) − \frac{4}{Z} ×(C − X) × \csc(2α))² + 1}$. I've attempted to find similar expressions for the transversal intersection by working backwards from numerically-calculated values, without any success yet.

For example, by plotting the two curves in graphing software and numerically calculating the parameter values of their intersections to ten decimal places, I created this plot, which shows the $γ_{min}$ value against the pressure angle $α$ (including some technically-invalid negative values of $α$ to get a broader sample size) for four different values of $Z$, all with $X = 0$ and $C = 1$: https://i.stack.imgur.com/9zOb1.png

The green lines show the known expression $γ_{min}(α) = \frac{4}{Z} × (C − X) × csc(2 × α)$, while the red points are samples of the unknown expression for which I am searching. I've been trying to fit a curve to the red points, unfortunately without any success yet.

The points

*appear*to trace out a sine wave of the form $γ_{min}(α) = A × \sin(ω × α + φ)$, where

- $A$ is the amplitude,
- $ω$ is the angular frequency, and
- $φ$ is the phase.

Additionally, my attempts to fit a sine curve to the numerically-calculated points suggest that the sampled function isn't a perfect sine wave and may be modified with some other term(s), the effect of which is particularly noticeable for small values of $Z$. At the very least, the red points for $Z = 8$ have proved much more difficult to get a sine wave to conform to than the points for higher $Z$ values. My best attempts so far are:

- $0.794 × \sin(1.20 × α + 1.70)$ for $Z = 8$,
- $0.664 × \sin(1.50 × α + 1.77)$ for $Z = 12$,
- $0.484 × \sin(2.18 × α + 1.83)$ for $Z = 24$, and
- $0.248 × \sin(4.30 × α + 1.89)$ for $Z = 96$.

**AN APPEALING BUT FLAWED APPROACH:**

Another approach that has been tried to solve this is to create the new variables

- $A = Z × \cos(α)$,
- $B = Z − 2 × (C − X)$,
- $ψ(θ) = θ − \tan(α) + α$, and
- $ξ(γ) = 2 × \tan(α) × (C − X) + γ × Z$.

Because the two curves intersect when $x_1(θ) = x_2(γ)$ and $y_1(θ) = y_2(γ)$, with these re-written equations we can see that an intersection exists when $A = B$, $A × θ = ξ(γ)$, and $γ = ψ(θ)$. This gives the following equations:

- $Z × \cos(α) = Z − 2 × (C − X)$
- $Z × \cos(α) × θ = 2 × \tan(α) × (C − X) + Z × γ$
- $γ = θ − \tan(α) + α$

- $α = \arccos\left(\frac{Z − 2 × (C − X)} {Z}\right)$
- $θ = \frac{α − \sin(α)} {\cos(α) − 1}$
- $γ = θ − \tan(α) + α$

*does*produce an intersection of the two curves, but it fails for my purposes on several accounts:

- The first of these three equations fixes the value of alpha in relation to the other three variables $Z$, $X$, and $C$. In reality, alpha is an independent variable whose value in this problem is constrained by, but not dependent on, the values of the other variables. Moreover, for most combinations of values of the other three variables, the value produced by this equation is outside the domain of alpha that I specify at the beginning of this question.
- My problem has the specific constraint that a solution must be valid when $α < \arcsin\left(\sqrt{\frac{2}{Z} × (C − X)}\right)$, and the value of alpha from that first equation is always greater than or equal to the right side of this inequality, so it is always invalid.
- The second equation produces a value of theta that, using the value of alpha from the first equation, is always negative or zero. For my problem the value of theta must be strictly greater than or equal to zero, so this part is also almost always invalid.

**THE CONTEXT:**

The curve with parameter theta is the involute face curve of a tooth on an involute gear, while the curve with parameter gamma is the trochoid root curve of the same tooth. These curves are naturally generated in real life by the gear-shaping process called

*hobbing*, without needing any fancy math. Representing them in a computer, which I want to do, is more difficult. The shapes of these curves are defined by four variables:

- $α$, the angle of the contact force between meshed gear teeth, called the
*pressure angle*; - $Z$, the number of teeth on the gear;
- $X$, the
*profile shift coefficient*, specifying how far in or out the cutting tool is moved compared to cutting a standard gear profile; and - $C$, the
*clearance factor*, specifying how much clearance there is between the tooth roots on one gear and the tooth tips on a meshing gear as a multiple of the overall tooth height.

*module*or

*pitch*, which describes the overall size of the gear. Because this variable is a uniform scaling factor, it has no effect on the angles involved or on the values of theta and gamma, so I have left it out of the equations for the sake of simplicity.

When $α ≥ \arcsin\left(\sqrt{\frac{2}{Z} × (C − X)}\right)$, the involute face curve transitions smoothly into the trochoid root curve (with a touching intersection). However, when $α < \arcsin\left(\sqrt{\frac{2}{Z} × (C − X)}\right)$, the root curve cuts off some of the face curve (with a transversal intersection). This is called

*undercutting*and is in general undesirable, as it reduces the strength of the gear. However, small amounts of undercutting are tolerated in many situations. I want to find the point on each curve where this undercutting occurs so I can accurately draw an undercut gear in software.

Lawton

12

Join Matchmaticians Affiliate Marketing
Program to earn up to 50% commission on every question your affiliated users ask or answer.

- closed
- 531 views
- $50.00

### Related Questions

- I am looking for a formula to calculate the radii of an ellipsoid from coordinates of non-coplanar locii on its surface.
- Find the area of the shaded region
- Optimization problem
- Can enough pizza dough be made to cover the surface of the earth?
- Find $n$ such that $\lim _{x \rightarrow \infty} \frac{1}{x} \ln (\frac{e^{x}+e^{2x}+\dots e^{nx}}{n})=9$
- Multivariable Calc: Vector equations, parametric equations, points of intersection
- Stokes' theorem $\int_S \nabla \times F \cdot dS= \int_C F\cdot dr$ verification
- How to recalculate 2D polygon side lengths when tilt is applied in 3D space?

I guess your offer should be higher for such open ended questions. The required computations get ugly quickly.

@SavionF What do you think would be a fair amount?

Think about how long it would take someone to figure this out and write it down, and how much the time of such qualified professional with a graduate degree is worth.