This is the analysis that I promised.
In the drawing below are shown two views of a rigid-body curve (red), first in plan view (x, z), and second in a right-side view (y, z). The (x, y, z) axes of the coordinate system are shown. The circular arc is generated by a parameter, ϴ. ϴ is termed the generator of the curve. For compatibility with MSTS, the circular arc in plan view must be rotated to represent a track gradient. An arbitrary point on the curve is shown as P, P(x, z) in the plan view and P(y, z) in the right-side view. The curve is seen in the right-side view as a straight black line.
Quadrant.gif (5.72K)
Number of downloads: 3
To develop formulas for the shape of the curve developed by MSTS, we refer to the two views in the drawing above. First, we develop expressions for x’, y’, and z’, which are coordinates in the unrotated plan view (x, z):
x’ = -R (1 - cos ϴ)
y’ = 0 (before the rotation to achieve a gradient, y’ is 0)
z’ = R sin ϴ
Then we rotate the plane of the curve by angle ϵ to achieve our grade expressions:
x = x’
y = z’ sin ϵ
z = z’ cos ϵ
Substituting expressions for x’ and z’,
x = -R (1 - cos ϴ)
y = R sin ϴ sin ϵ
z = R sin ϴ cos ϵ
The red curve in the plot below is the formula for y.
Plywood.gif (4.32K)
Number of downloads: 0
Let’s spot-check the ends as a “sanity” check. First, for ϴ = 0:
x = y = z = 0 CHECK
Then, at ϴ = π/2 (90°) and for, say, ϵ = 3°:
x = -R = -100
y= R sin ϵ = 5.233595624 COMPUTED WITH DOUBLE-PRECISION
z = R cos ϵ = 99.86295348 COMPUTED WITH DOUBLE PRECISION
Note that neither y nor z is an expected value. The value y is way below the expected grade for a full quadrant of arc (shown by the green line), and z is 137 mm short, a consequence of the rigid body rotation.
With those unexpected results, let’s take a look at the slope of the curve.
Slope in a curve in three-space (3D space) is defined as the rate of change of vertical distance (y) with respect to the distance, call it s, travelled along the curve. Fortunately for us, the curve is a circular arc, and s = R ϴ.
Deriving an expression for the slope requires us of differential calculus. The slope is defined as:
dy/ds = (dy/dϴ) (dϴ/ds)
But y = R sin ϴ sin ϵ and s = R ϴ.
dy/ds = (R sin ϵ cos ϴ) (1 / R) = sin ϵ cos ϴ
And here’s the plot:
Slope.gif (3.75K)
Number of downloads: 0
The initial slope is sin 3° = 0.05233595624.
Understand that this consequence is true for all MSTS track – not just dynamic track. For static sections, the same thing is true once arc lengths grow long. The effect is there; the only question is whether it’s noticeable.
Somewhere in this topic in recent weeks, I said something that was all wet. I said that grades are measured by rise over run. That’s fine – and correct, but which is which? Rise is rise, the change in elevation. I claimed that run is the distance along the grade – the hypotenuse. Wrong! Trains magazine corrected me:
"In North America, gradient is expressed in terms of the number of feet of rise per 100 feet of horizontal distance. Two examples: if a track rises 1 foot over a distance of 100 feet, the gradient is said to be "1 percent;" a rise of 2 and-a-half feet would be a grade of "2.5 percent." In other parts of the world, particularly Britain and places with heavy British influence, gradients are expressed in terms of the horizontal distance required to achieve a 1-foot rise. This system would term the above examples "1 in 100" and "1 in 40," respectively."
For a 3° grade, the grade slope ratio would be tan 3° = 0.05240777928.
Please let me know of any errors or misstatements in this. In the coming week, I’ll try to crank out a formula for what the space curve should be to produce a constant specified slope on a circular arc curve.