Open access peer-reviewed chapter

Investigation of Singularities in a 3/4 CMG Configuration with Mixed Skew Angles

By Jonathan W. Lang

Submitted: January 4th 2019Reviewed: June 16th 2019Published: January 15th 2020

DOI: 10.5772/intechopen.88087

Downloaded: 133

Abstract

This paper discusses the singularities that exist within a 3/4 CMG configuration when the CMGs are placed at mixed skew angles. CMGs are typically mounted with the same skew angles and are fixed throughout the spacecraft’s lifetime. Changing these skew angles can bring about unique attributes for the spacecraft such as an increased pitch, roll or yaw capability. Mapping out these singularities when each CMG is mounted differently can show an engineer how to maximize these capabilities and enhance a spacecraft’s mission completion ability. Using singularity penetration logic, the spacecraft’s attitude controls system can pass through these singularities. These singularities would best be avoided to provide optimal control. Finding these limited singularity penetration regions is the focus of this paper. Different mixed skew configurations appear to be more ideal than others for spacecraft that focus on maneuvers about only one axis of rotation.

Keywords

  • mixed skew angles
  • 3/4 CMG configuration
  • singularity penetration
  • CMG configuration singularity regions

1. Introduction

CMGs have become a staple in the space community as a means to accomplish pointing, tracking, and acquiring. CMGs create a torque by rotating or “gimbaling” the CMGs angular momentum vector. The change in angular momentum is how the CMG produces torque and creates movement for the spacecraft.

The Torque axis (τ) rotation for a CMG is the axis around which the spacecraft maneuver is accomplished. The Torque Axis and the orthogonal relationship between the CMG’s gimbal axis (δ) and angular momentum axis (H) can be seen in Figure 1. These relationships demonstrate how one may discover the Torque axis and provide the intuition behind singularity generation in this configuration [1].

Figure 1.

CMG angular momentum-torque-gimbal axis representation.

2. CMG 3/4 configuration model

The 3/4 configuration will have the CMGs positioned in the same configuration around the spacecraft throughout the analysis in this paper and can be seen in Figure 2.

Figure 2.

3/4 configuration of CMGs on a spacecraft.

The 3/4 configuration [2] is a simple form of spacecraft CMG placement which provides an environment where singularities exist.

All modeling in this project will be accomplished through MATLAB—Simulink modeling software. A step size of 0.001 s was found to be the most accurate while using ODE4 Runge-Kutta Integration Solver.

3. CMG angular momentum projected on the body reference axes

The three CMG’s each have an angular momentum vector (Hi) which can be seen in Figure 2 and is projected onto the x, y, z reference frame axes in Eqs. (1)(3). These three vectors rotate about each of the CMG’s gimbal axis (δi).

Hx=H1cosθ1+H2cosβ2sinθ2+H3cosθ3E1
Hy=H1cosβ1sinθ1H2cosθ2+H3cosβ3sinθ3E2
Hz=sinβ1H1sinθ1+sinβ2H2sinθ2+sinβ3H3sinθ3E3

A new matrix [A] made up of the spatial gradient of Eqs. (1)(3) can be generated [2]. This new matrix is how singularities can be discovered numerically.

Hθ=HxθiHyθiHzθi=A
A=sinθ1cosβ2cosθ2sinθ3cosβ1cosθ1sinθ2cosβ3cosθ3sinβ1cosθ1sinβ2cosθ2sinβ3cosθ3E4

Torque is generated by the time rate of change in the angular momentum, which may be expressed in a chain rule of derivatives and solved for the gimbal rates.

T=Ḣ
Ḣ=Hθdt
Ḣ=Aθ̇
θ̇=A1ḢE5

Now there is a relationship between the time rate of change in the angular momentum (Ḣ) and the time rate of change of gimbal axis rotation (θ̇)found in Eq. (5). The required gimbal axis rotation for a commanded torque may be applied for the appropriate spacecraft maneuver. The one possible danger associated with this is the inversion of the matrix A.

4. Accuracy of inverting the Amatrix

Inverting the [A] matrix becomes an integral part of determining the creation of the gimbal axis rotation (θ̇)as seen in Eq. (5). Analysis was done on eight different inversion models, the principal one on which all others are compared was the Moore-Penrose matrix inversion [3]. The other seven cases are found in Table 1.

Case[A] Inversion modelMatrix-Norm difference
1A10.872826646563208
2invA0.872826646563208
3pinvA0
4A\eye(sizeA)0.872826646560334
5LUDecomposition0.872826646560693
6Analytical Matrix6.67204727298e+04
7Analytical Formulas6.67204727298e+04

Table 1.

Matrix-norm difference between each case and the Moore-Penrose matrix inversion.

The pinvAmatrix inversion model turns out to be the best suited model as seen in Table 1. This seems to be intuitive since the pinvAis a reiteration of the Moore-Penrose inversion.

5. Applying CMG torque to the spacecraft

Now that the CMG torque has been developed it can now be applied to the spacecraft.

τ=Jω̇+ω×E6

Using Euler’s dynamical Eq. (6), the angular momentum and cross-coupled disturbances can be used to determine the resulting three-axis rotation [4].

6. Singularity penetration logic

Singularity existence has been discussed at length in [2].

Singularities exist when the determinant of the Amatrix from Eq. (4) approaches zero [2]. Penetration Logic takes advantage of the fact that as the inverse condition number approaches zero the singularity would be met. At a user defined threshold of 1 e-6 penetration logic tells the commanded Gimbal motor rates to repeat the last iteration (before the inverse condition number crossed the user defined threshold and reaches singularity), see Figure 3.

Figure 3.

Singularity penetration with unit delay logic.

7. PDI controller tuning

The PDI controller will be used to generate a control signal from the commanded input signal and the spacecraft feedback signal. The topography for the PDI controller can be seen in Figure 4.

Figure 4.

PDI controller topography.

The PDI controller accepts both the commanded angle and angular velocity as well as the feedback angle and angular velocity and avoids using the derivative function.

Tuning the gains become the next step in building the controller section of the spacecraft. Three different tuning techniques will be covered:

7.1 Ziegler-Nichols tuning

This tuning technique requires the gain margin (Ku) and bandwidth frequency (ωb) for the system to be found. Once these two values are discovered the equations in (7) may be used to find the required gains.

Kp=0.6Ku
Ki=2Kpωn
Kd=Kpωn8E7

7.2 Manual tuning

Several design criteria must be considered when utilizing the manual tuning technique. Rise time (tr), damping ratio (ζ), and settling time (ts) are the design characteristics required and may be used in the equation in (8) to find the new gains.

T=10ζωn
Kp=Jωn2+2ζωnT
Ki=Jωn2T
Kd=J2ζωn+1TE8

7.3 Tuning using the linear quadratic regulator function

The LQR function in MATLAB requires the state space form of the control system. Using the form ẋ=Ax+Buand y=Cx+Duthe state space form of the torque equation in the inertial frame can be constructed. The state space from in Eq. (9) was developed and used to derive the gains from the LQR function.

A=0100
B=0J1
C=10
D=0E9

The three tuning techniques were completed, and the gain values were calculated. These gains have been compiled and entered in Table 2 [5].

CaseTuning techniqueKpKiKd
1Ziegler-Nichols16.200.7884.65
2Manual36.761.4155.67
3LQR function1.000.1011.45

Table 2.

Gain solutions for different tuning techniques.

Using the gains found in Table 2, three differing step functions can be found and analyzed. Ziegler-Nichols was found to be the optimal tuning method for the controllers. The LQR function is the worst way of tuning the controller. The overshoot is largest of the three tuning methods and the settling time takes longer than 100 s [5].

8. Plotting singularity regions with mixed skew angles

A series of mixed skew profiles were selected to perform initial analysis on. Upon investigation, certain profiles were further analyzed to see how they would be advantageous for a spacecraft with a specific maneuver requirement such as maneuver in pitch, roll, or yaw. Figure 2 demonstrates how these skew angles (β) can change for each of the CMGs in the 3/4 configuration. When these skew angles are altered, new singularity regions are developed and create a hazard for spacecraft control. Five series of singularity plots have been generated to attain an initial picture of how the skew angles affect those singularity regions (Tables 37).

ProfileCMG 1-βCMG 2-βCMG 3-β
1153060
2301560
3306015
4156030
5601530
6603015

Table 3.

Series 1 CMG βangle configurations.

Series 1. Using skew angles (β)-15, 30, 60°.

ProfileCMG 1-βCMG 2-βCMG 3-β
1204080
2208040
3402080
4408020
5802040
6804020

Table 4.

Series 2 CMG βangle configurations.

Series 2. Using skew angles (β)-20, 40, 80°.

ProfileCMG 1-βCMG 2-βCMG 3-β
1306090
2309060
3603090
4609030
5903060
6906030

Table 5.

Series 3 CMG βangle configurations.

Series 3. Using skew angles (β)-30, 60, 90°.

ProfileCMG 1-βCMG 2-βCMG 3-β
104590
209045
345090
445900
590045
690450

Table 6.

Series 4 CMG βangle configurations.

Series 4. Using skew angles (β)-0, 45, 90°.

ProfileCMG 1-βCMG 2-βCMG 3-β
103060
206030
330060
430600
560330
660300

Table 7.

Series 5 CMG βangle configurations.

Series 5. Using skew angles (β)-0, 30, 60°.

The singularity plots associated with these series are found in Appendix A of this report. Each of these plots was analyzed to see breadth and depth of singularity surfaces internal to the saturation region of the 3/4 CMG configuration depicted in Figure 2.

9. Determining singularity free regions

The five series of singularity plots were observed to determine internal singularity free regions. The objective was to find an internal region in which the spacecraft could maneuver without running into a singularity. Also, if there was a singularity, would the spacecraft be capable of passing through a small amount of singularities in order to continue maneuvering to a commanded rotation. Singularity penetration is a feasible concept however the goal of this paper was to use it sparingly.

The singularity regions found in Figure 5 demonstrate how with varied skew angles new singularity free regions may be discovered and utilized to a spacecraft designer’s advantage. Although the angular momentum plotted in Figure 5 first interacts with a singularity at a value of 0.73, it only hits the singularity for an instance. Using singularity penetration logic discussed in part VI may allow the spacecraft to operate in an angularity momentum regime at H = 1.0 or greater.

Figure 5.

Singularity regions in the X-, Y- and, Z rotation for 0, 45, and 90° skew angles.

The maximum angular momentum achieved by the mixed skew configurations in series 1–5 has been plotted and the data was consolidated into Table 8.

SeriesCMG 1-βCMG 2-βCMG 3-βMax H*
11530600.54
3015600.51
3060150.16
1560300.16
6015300.51
6030150.54
22040800.34
2080400.17
4020800.33
4080200.17
8020400.33
8040200.53
33060900.12
3090600.41
6030900.51
6090300.41
9030600.51
9060300.25
4045900.73
090450.22
450900.41
459000.22
900450.41
904500.73
5030600.52
060300.39
300600.52
306000.39
600300.52
603000.53

Table 8.

Series 5 CMG βangle configurations.

Maximum Ang. Momentum (H) in singularity free region.


10. Conclusion

Mixed skew angles bring a new variety and flexibility in spacecraft design. These new CMG configurations enable engineers to now explore new singularity free regions and push spacecraft to possibly operate with higher levels of momentum.

As seen in Table 8, when the opposing CMG configurations have a dramatic change in configuration, the Angular Momentum is typically higher. This is not always the case, as seen in series 3, configurations 1 and 6 (the most dramatic change between CMG 1 and 3 in the series) demonstrate the worst angular momentum possibility. Several of these configurations cater to different requirements of the spacecraft. For example, Series 3 Configuration 5 (ref. Appendix A) may be more suitable for a spacecraft that requires movement about the roll and yaw axes. Series 3 Configuration 2 may be more suitable for pitch and roll spacecraft. Appendix A is meant to be used by spacecraft designers to design a spacecraft suitable to the requirements needed.

Appendix A

Singularity Plots

Series 1:

Configuration 1:

Configuration 2:

Configuration 3:

Configuration 4:

Configuration 5:

Configuration 6:

Series 2:

Configuration 1:

Configuration 2:

Configuration 3:

Configuration 4:

Configuration 5:

Configuration 6:

Series 3:

Configuration 1:

Configuration 2:

Configuration 3:

Configuration 4:

Configuration 5:

Configuration 6:

Series 4:

Configuration 1:

Configuration 2:

Configuration 3:

Configuration 4:

Configuration 5:

Configuration 6:

Series 5:

Configuration 1:

Configuration 2:

Configuration 3:

Configuration 4:

Configuration 5:

Configuration 6:

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Jonathan W. Lang (January 15th 2020). Investigation of Singularities in a 3/4 CMG Configuration with Mixed Skew Angles, Advances in Spacecraft Attitude Control, Timothy Sands, IntechOpen, DOI: 10.5772/intechopen.88087. Available from:

chapter statistics

133total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Control Moment Gyroscope Skew Angle Variation and Singularity Penetration

By Kyle A. Baker

Related Book

First chapter

Stochastic Artificial Intelligence: Review Article

By T.D. Raheni and P. Thirumoorthi

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us