Title of the UAIDE Paper (1964: Los Angeles, August 12-14) was THE ACTUAL MOTION AND ORIENTATION OF AN ORBITING SPACECRAFT VIEWED IN THE FORM OF A MOTION PICTURE.
NASA Technical Note received August 28, 1964)
The motion and orientation of an orbiting spacecraft are normally represented by vectors and angles in the various coordinate systems (celestial inertial coordinates, spacecraft coordinates, geodetic coordinates, etc.); and it is difficult to interpret spacecraft behavior in terms of actual orbital position and attitude. Motion and orientation information for a particular satellite (OGO) have been analyzed by the 1107 computer (Univac) and the 4020 microfilm plotter (Stromberg-Carlson), and numerical data on satellite behavior were made available. The data were then used to create a motion picture film illustrating satellite attitude in orbit.
OGO I (1964 54A) was successfully launched on 5 September 1964 from Cape Kennedy.
The National Aeronautics and Space Administration will soon orbit the first of a series of satellites called Orbiting Geophysical Observatories (OGO), which fall into the category of second-generation satellites. Satellites launched in the past have carried few scientific experiments, and the spacecraft was designed around these. Most of the early satellites were characterized by low weight, inflexible data systems, and spin-stabilization.
The era of an observatory in space is now here. The design of a spacecraft observatory may incorporate fifty or more diversified experiments, have versatile data-handling equipment to accommodate any experiment requirements, and employ sophisticated attitude-control systems. Close control of spacecraft attitude is one of the most important requirements of an orbiting observatory. The OGO (Figure 1 and Figure 2) employs a system allowing orientation of any particular experiment with the earth, the sun, or the orbital plane.
Experiments may be mounted on any one of three sections of the spacecraft, each of which maintains a unique orientation. The main body of the satellite is oriented towards the earth; the solar paddles which contain solar-oriented experiment packages (SOEPs) are oriented toward the sun; and the orbital-plane experiment package (OPEP) is oriented into the orbital plane. The OPEP maintains its front surface into the orbital plane, using a caged gyro as a sensing element (Figure 3).
The sophisticated attitude-control system of OGO, while one of its most important characteristics, is also one of the most critical. The success of most experiments depends upon accurate spacecraft orientation, and in addition, some experiments may be rendered useless by a complete, or even a partial failure of the attitude-control system.
The Orbiting Geophysical Observatory has three independent sets of spacecraft coordinates, which are maintained with respect to the body, the SOEPs and the OPEP, respectively. The body is oriented towards the earth; the SOEPs are oriented toward the sun; and the OPEP is oriented in the direction of motion of the spacecraft. The three axes of each set of coordinates are mutually perpendicular, and are defined as follows:
We must define two angles between these coordinate systems:
These nine coordinate axes and two angles can be used to determine the attitude of the three spacecraft subsystems, but their utility is limited since they are coordinate axes for the spacecraft only.
Since the attitude of the spacecraft must be presented with respect to the earth and the sun; the use of spacecraft coordinates is of limited value. This use would correspond to keeping the spacecraft fixed and having the earth and sun revolve about it. The solution is to rotate the spacecraft coordinates into celestial inertial coordinates, having the earth at the center of the system (Figure 4). The Xc axis lies in the plane of the equator and is directed towards the first point of the star Aries. The Zc axis is perpendicular to the plane of the equator, pointing towards Polaris, the north star. The Yc axis lies in the plane of the equator such that it is mutually perpendicular to the Xc and Zc axes.
This system, using celestial inertial coordinates, has meaning to any observer since the satellite motion and orientation are, in fact, related to the earth rather than any other reference. The system is complicated, however, by the necessity to use 27 vectors in replacing the 9 of the spacecraft coordinate system, since each of the 9 spacecraft vectors has an X, a Y, and a z value in the new system. Tables 1 and 2 contain sample data which might be telemetered from a satellite.
These data might also be computed from satellite tracking data.
Table 1 contains information on the satellite orbital position and Table 2 contains information on the conversion of orbital position data to satellite attitude data.
Time of Data 183 120000. Local Time 2 hrs 48 min 1.00 sec Right Ascension 41.35395 Declination -15.74462° Position Vector PX PY PZ Magnitude 4877. km 4293.km -1832. km Velocity Vector VX VY VZ Magnitude -4.3564 km/sec 8.7097 km/sec 4.2609 km/sec Solar Vector SX SY SZ Magnitude -24,285,378.km 137,822,250. km 59,765,857. km Longitude 121. 77864° Height 373.81 km Latitude -15.74462° True Anomaly 13.95780° Sun Earth Satellite Angle 69.26325°
IDEAL MAIN BODY | |||||||||
---|---|---|---|---|---|---|---|---|---|
Coord Axis | XbXi | YbXi | ZbXi | XbYi | YBYi | ZbYi | XbZi | YbZi | ZbZi |
Dir Cosine | .5298 | -.2571 | .8082 | -.4442 | .7277 | .5227 | -.7225 | -.6359 | .2714 |
IDEAL PADDLE | |||||||||
Coord Axis | XpXi | YpXi | ZpXi | XpYi | YBYi | ZpYi | XpZi | YpZi | ZpZi |
Dir Cosine | .5298 | -.2571 | .8082 | -.1596 | .9057 | .3928 | -.8329 | -.3371 | .4388 |
IDEAL OPEP | |||||||||
Coord Axis | XeXi | YeXi | ZeXi | XeYi | YBYi | ZeYi | XeZi | YeZi | ZeZi |
Dir Cosine | -.4971 | .7506 | .4353 | -.4805 | -1796 | -.8584 | -.7225 | -.6359 | .2714 |
ACTUAL MAIN BODY | |||||||||
Coord Axis | XbX | XbY | XbZ | YbX | YbY | YbX | ZbX | ZbY | ZbZ |
Dir Cosine | .5298 | -.2571 | .8082 | -.4442 | .7277 | .5227 | -.7225 | -.6359 | .2714 |
ACTUAL PADDLE | |||||||||
Coord Axis | XpX | XpY | XpZ | YpX | YpY | YpZ | ZpX | ZpY | ZpZ |
Dir Cosine | .5298 | -.2571 | .8082 | -.1596 | .9057 | .3928 | -.8329 | -.3371 | .4388 |
ACTUAL OPEP | |||||||||
Coord Axis | XeX | XeY | XeZ | YeX | YeY | YeZ | ZeX | ZeY | ZeZ |
Dir Cosine | -.4971 | .7506 | .4353 | -.4805 | .1796 | -.8584 | -.7225 | -.6359 | .2714 |
After defining spacecraft coordinates, we may calculate their values in the celestial inertial coordinate system. Data from the expected orbit of OGO (after it is launched) will be used in the following analysis. From this data we may obtain all the information necessary to calculate the attitude of the spacecraft. Three vectors in the celestial inertial coordinate system are of primary concern:
The body Zb axis of the observatory is oriented towards the center of the earth's disc. Therefore, the unit vector Zb, is
Equation (1):
Zb = - P / |P|
The SOEP Yp axis points toward the center of the sun's disc, and the position of the sun S (as measured from the earth) is known. Thus the solar position vector from the spacecraft becomes
(S - P) / |S - P|, a unit vector, and the unit vector Yp, is equal to the solar position vector
Equation (2):
Yp = (S - P) / |S - P|
The OPEP Xe. axis is always oriented in the direction of the spacecraft velocity vector; so the unit vector Xe, is
Equation (3):
Xe = V / |V|
Xb can be calculated from the known values of Zb, and Yp Since
Equation (4):
cos φp = Zb • Yp / |Zb × Yp|
it follows that
Equation (5):
Xb = (Yp × Zb ) / sin φp = (Yp × Zb ) / sqrt(1-cos2φp)
Now the quantities Zb, Yp, Xe, and Xb have been determined, and Yb is the only unknown body coordinate remaining. It is seen to be:
Equation (6):
Yb = Zb × Xb
For the SOEPs, Yp has already been calculated, and
Equation (7):
Xp = Xb,
so it follows that
Equation (8):
Zp = Xp × Yp
For the OPEP, Xe has been calculated, and
Equation (9):
Ze = Zb
and from this one sees that
Equation (10):
Ye = Ze × Xe
Thus the actual orientation of the spacecraft can be determined from telemetered data on position vectors and the two angles between the coordinate axes, for any instant in time. The data contained in Tables 1 and 2, while accurately representing spacecraft motion and orientation, are difficult to interpret, and it is suggested that the data could be computer-reduced and presented for examination and analysis in visual form by employing a motion picture film. This technique, using an 1107 computer (Univac) and a 4020 microfilm plotter (Stromberg-Carlson) to aid in preparing the film, has been experimentally employed for the OGO, and could be extended to visually analyze the motion and attitude of a manned spacecraft such as Gemini.
The spacecraft is dimensioned in terms of the celestial inertial coordinate system, where the BODY, SOEP, and OPEP X, Y, and z axes are collinear with the celestial coordinate system.
Every intersection (In) of 2 or 3 lines is positioned in terms, then, of the celestial coordinate system. These intersections or corners are assigned a position in an array, so that point In (having position X= X0, Y = Y0, Z= Z0) would be assigned:
Equation 11:
Px = X0;
Py = Y0;
Pz = Z0.
The lines between points are defined by means of a double-subscripted array L(J, I): where J is the line number, and K has the value of 1 or 2 (indicating the start and end, respectively, of a line).
Thus, L(J, 1) equals the point (Im) where the line J starts and L(J, 2) equals the point (In) where the line J ends.
In order to describe the spacecraft, then, it is merely necessary to progress down the list of points L(J, K), starting with J = 1, and draw lines from point I1=L(J,1) whose dimensions are:
Equation 12:
Px (I1) = X1
Py (I1) = Y1
Pz (I1) = Z1
to point I2=L(J,2) whose dimensions are:
Equation 13:
Px (I2) = X2
Py (I2) = Y2
Pz (I2) = Z2
Now that the technique of dimensioning and defining the spacecraft has been explained, it is possible to show how the spacecraft model is rotated and translated so as to agree with its actual orientation.
This is accomplished by means of the actual BODY, SOEP, and OPEP X, Y, and z vectors as supplied by the OGO attitude-determination program in the computer.
Spacecraft motion and attitude may be easily visualized if the spacecraft coordinates (rectangular) are converted into celestial inertial coordinates {also rectangular) by relatively simple transformation techniques for each of the three sets of axes BODY, SOEPs and OPEP). In order to convert from one rectangular coordinate system to another, direction cosines are employed; and they relate the two coordinate systems by a series of cosine functions (Figure 5).
Since the body zb axis is always oriented to the earth, and is coincident with the spacecraft position vector; we shall represent the transformation of this axis from spacecraft coordinates into celestial inertial coordinates by Equation 16 which is in turn defined by Equations 14 and 15. The transformation of all other coordinate axes may be performed in a similar manner.
Let the position vector P (which is expressed in spacecraft coordinates) be
Equation 14:
P = δ(Xc, Yc, Zc)
and let the direction cosines of each component be defined by
Equation 15:
B1 3 = cos φ
B2 3 = cos θ
B3 3 = cos ψ
Therefore the position vector P is
Equation 16:
P = Xc cos φ + Yc cos θ + Zc cos ψ
or
P = B1 3Xc + B2 3Yc + B2 3Zc
This procedure, when applied to the other spacecraft coordinate axes, yields similar results.
The actual body Xb, Yb, Zb coordinates are defined in terms of the celestial coordinates system, where:
the Xb axis has direction cosines B1 1, B1 2, B1 3;
the Yb axis has direction cosines B2 1, B2 2, B2 3;
the Zb axis has direction cosines B3 1, B3 2, B3 3;
Figure 6 shows how the OGO spacecraft would appear if plotted as it is defined in the computer memory using the technique just described. The points, Ib, located on the body are rotated into the celestial coordinate system, using the following transformation:
Equation 17:
PRx (Ib) = B1 1Xb + B2 1Yb + B3 1Zb
PRy (Ib) = B1 2Xb + B2 2Yb + B3 2Zb
PRz (Ib) = B1 3Xb + B2 3Yb + B3 3Zb
When this transformation has been completed, the OGO spacecraft as defined in the computer memory, if plotted, would appear as shown in Figure 7. The actual coordinates xp, Yp, zp of the OPEP are also defined in terms of the celestial coordinate system where:
the Xρ axis has direction cosines Ρ1 1, Ρ1 2, Ρ1 3;
the Yρ axis has direction cosines Ρ2 1, Ρ2 2, Ρ2 3;
the Xρ axis has direction cosines Ρ3 1, Ρ3 3, Ρ3 3.
The points Iρ located on the SOEP are rotated into the celestial coordinate system using the following transformation:
Equation 18:
PRx(Iρ) = Ρ1 1 Xρ + Ρ2 1 Yρ + Ρ3 1 Zρ
PRy(Iρ) = Ρ1 2 Xρ + Ρ2 2 Yρ + Ρ3 2 Zρ
PRz(Iρ) = Ρ1 3 Xρ + Ρ2 3 Yρ + Ρ3 3 Zρ
When this transformation has been completed, the OGO spacecraft as defined at this point in the computer memory, when plotted in the Y, Z plane, would appear as shown in Figure 8.
The actual coordinates Xe, Ye, Ze of the OPEP are also defined in terms of the celestial coordinate system where:
the Xe axis has direction cosines E11, E12, E13;
the Ye axis has direction cosines E21, E22, E23;
the Ze axis has direction cosines E31, E32, E33;
The OPEP system as initially defined is displaced from its actual position by the following amounts:
Equation 19:
X = 0
Y = a
Z = 0
The points, Ie, located on the OPEP are rotated and translated into the celestial coordinate system using the following transformation:
Equation 20:
PRx (Ie) = E1 1Xe + E2 1Ye + E 3 1Ze + 0
PRy (Ie) = E1 2Xe + E2 2Ye + E 3 2Ze + a
PRz (Ie) = E1 3Xe + E2 3Ye + E 3 3Ze + 0
This last transformation completes the task of properly orienting the satellite in space. A completed plot of this object, as defined in the computer memory and plotted in the Y, Z plane, would appear as shown in Figure 9.
Once these three transformations (Equations 18, 19, and 21) have been performed, the spacecraft can be plotted from any vantage point by means of one further transformation.
Let the vantage point be XL, the vertical coordinate be ZL, and YL defined to make a right-handed system. (In other words, the YL and ZL coordinates define the plane of the plot as viewed from XL.)
Then
XL has direction cosines L1 1, L1 2, L1 3;
YL has direction cosines L2 1, L2 2, L2 3;
ZL has direction cosines L3 1, L3 2, L3 3.
The transformation for the entire spacecraft is:
Equation 21:
PTx(IL) = L1 1 PRx(IL) + L2 1 PRy(IL) + L3 1 PRz(IL)
PTy(IL) = L1 2 PRx(IL) + L2 2 PRy(IL) + L3 2 PRz(IL)
PTz(IL) = L1 3 PRx(IL) + L2 3 PRy(IL) + L3 3 PRz(IL)
In the present film, the observatory is shown maintaining its orientation throughout its first few orbits. Results were based upon the satellite's achieving its nominal orbit. No attempt was made to put into the film any other features, except for a true anomaly clock which shows the satellite position in orbit, with respect to the perigee. Figure 10 is a sample frame from the film.
The next step would be to incorporate the earth into the picture. Following from this, inclusion of the position of the sun would aid significantly in interpreting observatory behavior; and finally, the spacecraft should be made to orbit the earth, using the proper aspect ratio of the orbiting satellite and the earth.
The implementation of these additions has been calculated and, after the OGO has been launched, will be incorporated along with the satellite's telemetered attitude information. One use of this technique of data display is to analyze large quantities of orbit and attitude data and to determine certain effects upon the spacecraft not easily evaluated, such as: (1) the effect of darkness on the ability of the solar arrays to maintain orientation, (2) the orientation ability of the OPEP gyro at apogee (80,000 nautical miles) when the satellite velocity is at a minimum, (3) the reaction of the satellite during a yaw maneuver, and (4) the reaction of the horizon sensors, the SOEPs and the OPEP during perigee.
Since the orbital period of the OGO is more than 60 hours, the behavior of the satellite for this time can be evaluated in the 2 minutes required to display the film. In addition, this system is not limited to one spacecraft. One interesting application would be to display the attitude and motion of manned spacecraft such as Gemini or Apollo, where the capsule behavior is displayed in real time. Minor system modifications would allow pertinent system data to be simultaneously displayed, as shown by a hypothetical spacecraft display in Figure 11.
Furthermore, for a rendezvous of two orbiting spacecraft, maneuvers could be initiated prior to visual contact by means of a television link such that each pilot-astronaut could view the orientation and motion of both spacecraft simultaneously, and that the command station on the ground could view the television display.