WO2007132793A1 - 姿勢制御データ生成方法、及びそれを適用した姿勢制御装置 - Google Patents

姿勢制御データ生成方法、及びそれを適用した姿勢制御装置 Download PDF

Info

Publication number
WO2007132793A1
WO2007132793A1 PCT/JP2007/059809 JP2007059809W WO2007132793A1 WO 2007132793 A1 WO2007132793 A1 WO 2007132793A1 JP 2007059809 W JP2007059809 W JP 2007059809W WO 2007132793 A1 WO2007132793 A1 WO 2007132793A1
Authority
WO
WIPO (PCT)
Prior art keywords
attitude
control
posture
target amount
control target
Prior art date
Application number
PCT/JP2007/059809
Other languages
English (en)
French (fr)
Inventor
Toshio Kamiya
Ken Maeda
Tatsuaki Hashimoto
Shinichiro Sakai
Original Assignee
Nec Toshiba Space Systems, Ltd.
Japan Aerospace Exploration Agency
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nec Toshiba Space Systems, Ltd., Japan Aerospace Exploration Agency filed Critical Nec Toshiba Space Systems, Ltd.
Priority to EP07743244A priority Critical patent/EP2022718B1/en
Priority to CA2651814A priority patent/CA2651814C/en
Priority to US12/300,449 priority patent/US8186627B2/en
Publication of WO2007132793A1 publication Critical patent/WO2007132793A1/ja

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/28Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect
    • B64G1/283Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect using reaction wheels
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/244Spacecraft control systems
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/28Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect
    • B64G1/286Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect using control momentum gyroscopes (CMGs)
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/36Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/38Guiding or controlling apparatus, e.g. for attitude control damping of oscillations, e.g. nutation dampers

Definitions

  • the present invention relates to a posture control data generation method for generating posture control data during posture maneuver for a movable object such as a flying object, and a posture control device to which the posture control data is applied. Itoda
  • attitude change control of a flying object such as an artificial satellite having a flexible object
  • it is required to change to a desired attitude faster and with higher accuracy.
  • vibration of a flexible object is regarded as the most problematic. If the flexible object is shaken by the posture change, the vibration is transmitted to the main body of the flying object and an error occurs in the pointing direction of the main body. Also, once vibration occurs in a flexible object, it takes a considerable amount of time for the vibration to attenuate.
  • an important issue for attitude change control of a flying object having a flexible object is how to suppress the vibration of the flexible object.
  • Japanese Patent No. 3 6 2 3 7 4 7 discloses, for example, an attitude change control method for a triaxial satellite as a technique related to attitude change (attitude control) with respect to an artificial satellite or the like. Yes.
  • this attitude change control method frequency separation is not performed for the frequency components included in the target quantity in the process of generating the attitude control target quantity. For this reason, it cannot be avoided that a flexible object (artificial satellite flexible variable body) mounted on an artificial satellite or the like is excited, and a steady vibration of the flexible object remains even after the maneuver ends. In other words, it is necessary to wait for the attitude to settle (vibration convergence) until the observation mission can be started after the end of the maneuver. Therefore, the observation mission cannot be started because the directional stability is not ensured until the vibration converges. As a result, the total time required for the attitude maneuver (time from the start of attitude maneuver to the start of observation mission) Takes a long time.
  • Attitude maneuver control method of flexible satellite by gBang control (Reference 3: “L iu, Q., Wie, B .: R obust Ti me— Op tima 1 C ontrolof Uncertain F lexible S pacecraft, Journal of Guidance , Contro 1, and Dynamics, 15 (1992), pp. 597-604 ”). Disclosure of the invention
  • the Bang-Bang control in Reference 2 or Reference 3 described above is intended to use a thruster as a control actor, and for controlling the attitude of a control actor other than a thruster, for example, an artificial satellite. It is difficult to apply to Reaction Wheel (RW) or Control Moment Gy ro (CM G) to move the attitude of the satellite.
  • RW Reaction Wheel
  • CM G Control Moment Gy ro
  • the object of the present invention is to provide all the maneuver required for the posture of a movable object having a flexible object.
  • An object of the present invention is to provide a posture control data generation method that enables shortening of the body time.
  • Another object of the present invention is to provide a posture control data generation method that enables effective suppression of vibration in a posture maneuver even when a flexible object of a movable object has a plurality of flexible modes.
  • Another object of the present invention is to provide an attitude control device capable of reducing the overall time required for attitude maneuver of a movable object having a flexible object.
  • Another object of the present invention is to provide an attitude control device capable of effectively suppressing vibrations caused by an attitude maneuver even when a flexible object of a movable object has a plurality of flexible modes.
  • a posture control data generation method for generating posture control data in a posture maneuver of a movable object having a soft material, the vibration control unit generating vibrations during the posture maneuver. Based on this, a sampling function acquisition stage for obtaining a sampling function having no frequency component above a specific frequency, and a control target quantity generation stage for generating a control target quantity as a pre-frequency shaping type feed-forward control law using the sampling function. And a control data generation step of generating the attitude control data based on the control target amount.
  • control target amount generation step the control target amount may be generated using a s inc function capable of generating a posture profile waveform in which frequency components are completely separated as the sampling function.
  • acceleration / deceleration data may be generated as the control target amount by combining two positive and negative waveforms in the sampling function by the s inc function while shifting by one cycle.
  • control target amount generation step the control target amount may be generated by multiplying the sync function by a window function in order to minimize the vibration of the flexible object after the posture maneuver.
  • an attitude control device that controls the attitude of a movable object having a flexible object by an attitude maneuver, and uses a sampling function that does not have a frequency component greater than a specific frequency.
  • Pre-frequency shaping die feed forward quantity An apparatus characterized by having a control target amount generation means for generating as a control law is obtained.
  • the control target amount generation means may obtain the sampling function based on vibration of the flexible object during the posture maneuver.
  • it may have posture changing means for changing the posture of the movable object using posture control data generated based on the control target amount.
  • it may include a posture rotating means installed on the movable object for rotating the movable object.
  • the control target amount generation means may generate the control target amount using a sinc function that can generate a waveform of a posture profile in which frequency components are completely separated as the sampling function.
  • the control target amount generation means may generate acceleration / deceleration data as the control target amount by combining two positive and negative waveforms in the sampling function based on the sinc function while shifting by one cycle.
  • the control target amount generation means may generate the control target amount by multiplying the s inc function by a window function.
  • the movable object may be an artificial satellite, and the flexible object may be a flexible variable body part mounted on the artificial satellite.
  • the flexible variable body section may include a reaction wheel for rotating the attitude of the artificial satellite and a control moment gyro for moving the attitude of the artificial satellite.
  • the attitude control data generation method makes it possible to reduce the overall time required for the attitude maneuver of a movable object having a flexible object, and the flexible object of the movable object has a plurality of flexible modes. Even in this case, it enables effective suppression of vibration in the posture maneuver.
  • the attitude control device can reduce the overall time required for the attitude maneuver of a movable object having a flexible object, and even when the flexible object of the movable object has a plurality of flexible modes. Vibration due to the posture maneuver can be effectively suppressed.
  • FIG. 1 is a time chart for explaining a sampling function used in an attitude control data generation method according to an embodiment of the present invention.
  • (A) is two samplings of positive and negative polarities shifted by one period. The characteristics of the function are shown, and (b) shows the combined characteristics of the two sampling functions shown in (a).
  • FIG. 2 is a time chart for explaining the processing of the sampling function used in the attitude control data generation method according to the embodiment of the present invention.
  • A shows the composite characteristic shown in FIG. 1 (b).
  • B shows the window function characteristics for the peak part in the composite characteristics shown in (a)
  • c shows the characteristics obtained by multiplying the composite characteristics shown in (a) by the window function characteristics in (b). Show.
  • FIG. 3 is a time chart showing the peak portion in the window function characteristics shown in FIG. 2 (c) in relation to the amplitude with respect to the frequency in comparison with the peak portions in the component characteristics.
  • FIG. 4 is a time chart illustrating various characteristics of a maneuver profile that is a control target amount generated based on a sampling function in the attitude control data generation method according to an embodiment of the present invention.
  • FIG. 5 is a time chart illustrating the characteristics of the flexible mode displacement during attitude maneuver executed on a flexible object (satellite flexible variable section) by the attitude control device according to an embodiment of the present invention. a) shows the fundamental mode characteristic expressed by the relationship of the flexible mode displacement to time, and (b) shows the differential mode characteristic expressed by the relationship of the differential flexible mode displacement to time.
  • FIG. 6 is a functional block diagram showing the basic configuration of the attitude control device for artificial satellite to which the attitude control data generation method according to one embodiment of the present invention is applied.
  • This attitude control data generation method is a method for generating attitude control data in the attitude maneuver of a movable object having a flexible object, and has a frequency component greater than a specific frequency based on vibration of the flexible object during the attitude maneuver.
  • a sampling function acquisition stage that obtains a sampling function that does not exist, and a control target quantity generation stage that generates a control target quantity as a pre-frequency-shaped feed-forward control law using the sampling function.
  • attitude control data is generated at the control data generation stage based on this control target amount.
  • Figure 1 exemplifies the characteristics of the sampling function used in the control target quantity generation stage [expressed in relation to the amplitude with respect to time (s)].
  • the figure (a) shows the positive / negative / reverse polarity shifted by one period.
  • (B) is a composite characteristic (s 1 f 1 1 + 1) that combines and combines the characteristics of (a). 2).
  • FIG. 1 (a) the positive and negative polarities are shifted by one period that can generate the posture profile waveform as the characteristic s-f 1 1 and s-f 1 2 of the sampling function in the control target quantity generation stage described above. Use the sinc function.
  • Figure 1 (b) shows how the sampling function based on these sinc functions is added and combined to generate acceleration / deceleration data that is the control target amount as the combined characteristic (s—f 1 1 + 1 2). Yes.
  • Figure 2 shows a comparison of various characteristics [expressed by the relationship of amplitude to time (s)] when the sampling function used in the control target quantity generation stage is applied.
  • Figure 3 shows the peak part in the component characteristic (s—f 1 1 + 1 2) as the peak part in the characteristic (s—f 1 1 + 1 2) X (w ⁇ f) shown in FIG. 2 (c). In contrast to this, it is shown in relation to the amplitude with respect to the frequency (Hz).
  • the frequency components are completely separated even if the above-mentioned additive synthesis characteristic (s-f 1 1 + 1 2) is multiplied by the window function characteristic (w-f), that is, the generation of the control target amount
  • the frequency separation can be performed completely by using the sinc function as the sampling function in the generation of the control target amount at the stage.
  • the steady residual vibration of a flexible object after posture maneuver can be brought close to zero, and the observation missions performed through various posture maneuvers can be effectively implemented by the shortest time posture maneuver. Suggests that you can.
  • simply superimposing the sampling functions as shown in Fig. 1 (b) (s-fll + 1 2) will cause vibration to remain after the posture maneuver.
  • the characteristics (w-f) of the window function as shown in Fig. 3
  • the control target amount can be smoothly attenuated from the attitude maneuver to the end of the attitude maneuver. Yes.
  • the sampling function (sine function) is a function that is generally known not to include any frequency component above a specific frequency band ⁇ s (Hz).
  • y (t) si ⁇ ( ⁇ s t) Z w s t preparative Ru.
  • the first consideration is to use an acceleration / deceleration profile that does not cause steady vibration of the flexible object.
  • it is necessary to perform a posture maneuver that maximizes the torque capacity that can be generated by the actuator (torque generator).
  • the first consideration is the angular acceleration during the attitude maneuver of the satellite using the waveform of the sampling function (sinc function) as shown in Fig. 1 (a). This can be achieved by generating a degree profile (control target amount).
  • the angular acceleration profile generated using the sampling function is viewed in the frequency distribution, the distribution shown in Fig. 3 is obtained, and the waveform does not contain frequency components above a specific frequency determined by the sampling function. I understand that.
  • the peaks of the two sampling functions shown in Fig. 1 (a) are shifted by one period, and positive waves are used for acceleration and negative waves are used for deceleration.
  • This can be realized by generating an angular acceleration profile of acceleration / deceleration of the attitude of the satellite as a characteristic as shown in Fig. 2 (c) by superimposing the two waves.
  • the actuator is utilized while taking advantage of the first consideration. Acceleration / deceleration can be performed with the maximum capacity of
  • the attitude maneuver rule executed in the present invention is to perform control so that the Euler axis is obtained from the attitude before the attitude change and the target attitude after the attitude change, and the attitude of the artificial satellite is rotated around the Euler axis.
  • Figure 4 exemplifies the characteristics of the maneuver profile, which is the control target amount after the Euler axis determination that is generated in the control target amount generation stage.
  • observation mission period t Assuming that bs and posture maneuver period t mn v are alternately repeated, an offset angle of 3 deg (3 degrees) as shown in the lower row.
  • an offset angle of 3 deg (3 degrees) As shown in the lower row.
  • the attitude maneuver that rotates the attitude of the artificial satellite around the Euler axis is shown based on the angular acceleration characteristics shown in the upper row.
  • the posture control device to which the posture control data generation method is applied sends a control signal to the posture rotation means for rotating the posture installed on the flexible object and the posture change means for moving the posture.
  • It includes a basic configuration that performs posture maneuver from one stationary state to another.
  • this attitude control device is included as a control signal control element in order to suppress vibration generated by a flexible object during attitude maneuver, and also uses an attitude function that does not have a frequency component above a specific frequency. It includes posture maneuver control target amount generation means for generating a control target amount serving as posture control data during maneuver as a pre-frequency shaping type feedforward control law.
  • the posture maneuver control target amount generation means has a function of generating a control target amount by using a sinc function capable of generating a posture profile waveform that completely separates frequency components as a sampling function, and a sinc function. It has a function to generate acceleration / deceleration data as a control target by combining two positive and negative waveforms in the sampling function by 1 cycle, and further suppresses vibration of flexible objects after the posture maneuver ends To multiply the sinc function by a window function It is preferable to have a function of generating a control target amount by combining them.
  • a flexible object is effective when it is a flexible variable body part for an artificial satellite mounted on an artificial satellite. It is preferable that the flexible variable body portion for the artificial satellite includes a reaction wheel for rotating the attitude of the artificial satellite, a control moment moment jar for moving the attitude of the artificial satellite, and the like.
  • Figure 5 illustrates the characteristics of the flexible mode displacement during the attitude maneuver performed by the attitude control device for the flexible satellite body, and ( a ) shows the flexible mode displacement with respect to time (s). (B) relates to the fractional mode characteristic indicated by the relationship between the differential flexible mode displacement d / dt with respect to time (a).
  • Figs. 5 (a) and 5 (b) show the time history of the mode displacement for the satellite flexible variable part during attitude maneuver.
  • zero characteristic SI is obtained in the first mode (1 st mode), and negative between 10 and 20 seconds in the second mode (2 st mode).
  • Has two large peaks facing positive and has a characteristic S 2 with two large peaks facing positive and negative between 40 seconds and 50 seconds.
  • 10 seconds to The characteristic S 2 has two small peaks facing negative and positive in 20 seconds, and two small peaks facing positive and negative in 40 to 50 seconds.
  • the third mode (3 st mo de), the fifth mode (5 st mo de) to the eighth mode (8 st mo de)] show almost zero.
  • the force is zero in the first mode (1 st mode), and the force is S 1 in the second mode (2 st mode), 10 seconds. It has negative, positive, negative, medium, large, and medium peaks for 20 seconds, and positive, negative, and positive medium, large, and medium peaks for 40 to 50 seconds.
  • the fourth mode (4 st mode) there are three very small peaks that are negative, positive, and negative in the fourth mode (4 st mode).
  • attitude maneuver is performed using the control target amount generated based on the characteristic (s—f 1 1 +12) X (w ⁇ f) shown in FIG.
  • the attitude maneuver can be completed in the shortest required time, and the observation missions performed through various attitude maneuvers can be effectively performed by the shortest time maneuver.
  • Attitude Control System (ACS) 40 may be referred to as a disturbance torque (DT) signal (a suppression torque signal) from an external electric device (not shown).
  • DT disturbance torque
  • a suppression torque signal a suppression torque signal
  • DT disturbance torque
  • CMG control moment gyro
  • the reaction wheel (RW) and control moment gyroscope (CMG) can generate its reaction torque by rotating the flywheel. This is an actuator that can rotate the attitude of an artificial satellite.
  • the control moment moment gyro uses a gimbal that rotates at a certain rotational speed to rotate it using a gimbal that has a rotation axis perpendicular to the rotation axis of the flywheel. This is an actuator that can output the torque of the gyro effect and thereby control the rotation of the satellite.
  • the hardware part of the attitude control device (ACS) 40 is basically composed of a gyroscope (Gy ro S copes: GS) 1 as a gyro sensor that works as an attitude information sensor in the flexible variable body 50 for artificial satellites, and an artificial satellite.
  • a gyroscope Gy ro S copes: GS 1 as a gyro sensor that works as an attitude information sensor in the flexible variable body 50 for artificial satellites, and an artificial satellite.
  • the gyroscope (GS) 1 outputs an angular velocity signal co GS containing information obtained by detecting the angular velocity ⁇ rotated around the three orthogonal axes of the satellite relative to the inertial space (which can be considered as space here)
  • a sensor or an angle detection signal 0c S including information indicating that the angle 0 is detected can be output).
  • the optical sensor group (O S) 2 is a set of sensors that output an angle detection signal e os including information obtained by detecting angles with a star, the earth, and the sun.
  • the reaction wheel (RW), control moment gyro (CMG), gyroscope (GS) 1 and optical sensor group (OS) 2 described above are widely used in the field of satellite attitude control systems. Therefore, the detailed configuration and functions are not described here.
  • the gimbal angle / angular velocity sensor 3 outputs an angle Z angular velocity ratio signal 0 ⁇ / ⁇ ⁇ indicating the ratio of the angle Z angular velocity detected from the control moment gyro (CMG) gimbal.
  • the attitude control circuit (ACE) 30 is an angular velocity required for the attitude control of the satellite based on the angular velocity signal co GS from the gyroscope (GS) 1 and the angle detection signal 0 os from the optical sensor group (OS) 2.
  • attitude estimator for outputting a position estimation signal indicating the estimated value of the angle 6 beta: and (At titude E sti ma taor AE ) 4, attitude estimator (AE) from 4
  • the control target amount required for the attitude maneuver is indicated according to the Euler axis obtained from the initial attitude and the changed attitude before executing the attitude maneuver.
  • the target moment profile is an angle target signal indicating the target value of the angle e M , the output of the angular velocity target signal indicating the target value of the angular velocity ⁇ M , and the control moment based on the detection result of the attitude flow (At titude D rift)
  • CMG gimbal
  • Angle 0 angle target signal indicating a target value of e
  • generating an output of the angular velocity target signal indicating a target value of the angular velocity co G and the reaction product output wheel FF Te control signal Fidofu Owado control law for (RW) (RW) of Attitude Guidace Law (AG L) 20 and the angle target signal from the attitude guidance unit (AGL) 20
  • a subtractor that subtracts each estimated signal and the subtracted value from this subtractor are input as a signal indicating the posture flow (AD), and a control signal ⁇ FB (
  • the reaction wheel control unit (RW Controller) 5 that generates and outputs by PID control, and the control signal ⁇ FF ( R3 ⁇ 4 and control signal ⁇ ⁇ (RW from the reaction wheel control unit 5) ) and adding the And adder inputs the sum of this adder, a reaction wheel distributor (RW D istributor) 10 which control signal by pseudo-inverse Te produces outputs RW, angle from the posture guiding portion (AGL) 20 Target signal, angular speed Gimbal angle from target signal Angular speed ratio signal from angular speed sensor 3 Subtractor for subtracting ⁇ G / co G , and the subtraction value from this subtractor are input and controlled by PD control It includes a control moment gyro gimbal controller (CMG—G imb a 1 Controller) 1 1 that generates and outputs CMG as a control signal for the moment gyro (CMG) gimbal.
  • CMG—G imb a 1 Controller control moment gyro gimbal controller
  • the posture guidance unit (AGL) 20 performs an initial posture in kinematics before executing the posture maneuver based on the estimated value of the angular velocity ⁇ ⁇ ⁇ and the estimated value of the angle ⁇ ⁇ ⁇ from the posture estimator (AE) 4.
  • the angle as a target posture profile that determines the Euler axis from the changed posture and the control target amount required for the posture maneuver according to the Euler axis
  • attitude maneuver time control target amount profile generator as a control target amount generating means for outputting an angular velocity target signal indicating a target value of the angle target signal and the angular velocity omega Micromax indicating a target value of Micromax and (NME P rofiler) 6,
  • Control signal ⁇ FF ( c MG) for feedforward control law for control moment jar (CMG) by inputting angle target signal and angular velocity target signal from attitude maneuver control target amount profile generator 6
  • Torque profile generator / distributor 7 that generates and outputs a control signal ⁇ ' FF (RW) for the feedforward control law for the reaction wheel (RW) by reverse fluctuation, and the reaction wheel control unit
  • torque profile And additional torque profile generator for generating and outputting a signal (To rque P rofiler) 9, and a control signal ⁇ FF (RW) from the torque
  • each part of the attitude control circuit (ACE) 30 can be configured with only a hardwired logic circuit, and has a configuration in which a computer such as a CPU is built in and control logic is executed in software.
  • a computer such as a CPU
  • control logic is executed in software.
  • the attitude guidance circuit (AGL) 20 in the attitude control circuit unit (ACE) 30 the attitude maneuver control target amount profile generator 6, the torque profile generator / distributor 7, Jinnol guidance part 8 for control moment gai mouth, attitude guidance part (AGL) 20 Outside the attitude estimator (AE) 4, reaction wheel control part 5, reaction wheel distributor 10 and control moment gyro
  • Each part of the gimbal control part 1 1 is assumed to be realized by software.
  • the output of the gyroscope (GS) 1 and the output of the optical sensor group (OS) 2 are used for artificial satellites. It has a function to estimate the current posture in flexible variable body part 50.
  • Attitude maneuver control target amount profile generation unit 6 is a part that executes a pre-frequency shaping type feed-forward control law for posture maneuver, and generates a control target amount (target posture) profile for that purpose.
  • Torque profile generator Z distributor 7 functions to calculate the control torque required to achieve the target attitude profile, and to avoid the amount of inter-axis interference due to the angular momentum of the satellite. The required decoupling amount is calculated at the same time, and the distribution function has a function to distribute the calculated control torque to the control moment gyro (CMG) and the reaction wheel (RW).
  • the control moment gyro ginnol guide 8 calculates the gimbal angle profile and gimbal angular velocity profile of the control moment gyro (CMG) required to output the control torque distributed for the control moment gyro (CMG). To do.
  • the gimbal control unit 1 1 is based on the gimbal angle error and the gimbal angular velocity error obtained by comparing the current gimbal state obtained from the gimbal angle angular velocity sensor 3 with the target gimbal angle state. Performs minor loop control to perform the target movement.
  • the reaction wheel control unit 5 calculates the feed pack control torque for the posture to perform the target motion based on the posture angle error and the posture angular velocity error.
  • the reaction wheel distributor 10 contributes by distributing the feed pack control torque necessary for the posture to perform the desired movement to each reaction wheel.
  • the components of the above parts are the attitude maneuver control target amount profile generation unit for calculating the optimal feedforward control law that does not excite the natural frequency of the flexible object. Except for the six components, the well-known configuration frequently used for attitude control of artificial satellites is introduced.
  • the operation of the posture target manipulator profile profile generator 6 will be described in detail.
  • the Euler axis obtained from the posture before the posture change and the target posture after the posture change is shown in FIG.
  • the attitude maneuver that rotates the attitude of the artificial satellite is executed. That is, the attitude maneuver control target amount profile generation unit 6 obtains the Euler axis from the initial attitude and the changed attitude before executing the attitude maneuver, and shows the Euler axis in the upper part of FIG. 4 according to the Euler axis.
  • the characteristics obtained by calculating the change in the angular acceleration of the satellite around the Euler axis that can realize attitude maneuver with the shortest path like this are obtained, and the attitude control of the satellite with such characteristics one by one during attitude maneuver
  • the target amount is generated and given to the attitude guidance unit (AGL) 20 that executes the attitude guidance law by the control moment gyro (CMG) and the V action wheel control unit 5 that forms the PID control loop. If the satellite rotates according to this target amount, the vibration of the reaction wheel (RW) and control moment gyroscope (CMG) in the flexible satellite body 50 for artificial satellites, which is the flexible object for attitude change control of the shortest path, is minimized. It will be possible to achieve it while keeping it down.
  • the satellite can excite only the minimum vibrations in the flexible object through the shortest path that is theoretically required, and this makes it possible to perform highly accurate attitude maneuvers. it can.
  • reaction wheel (RW), control moment jay port (CMG)] in the flexible variable body part 50 for artificial satellite performing attitude maneuver by the attitude control device 40 shown in FIG. 6 is reacted. It is also possible to replace the wheel (RW) with a thruster and replace the gimbal guide 8 for the control mouth moment jay mouth with a thruster for the reaction wheel-(RW). A similar effect can be expected. However, in this case, the reaction wheel (RW) generally generates a small control torque more precisely, although the maximum output torque is smaller than the control moment jay port (CMG). Therefore, such a configuration is suitable when it is required to change the posture at a small angle with higher accuracy.
  • attitude control device 40 can be used without any limitation on the type of the actuator simply by making a slight change, and can perform the attitude maneuver that effectively utilizes the performance of each of the actuators. .
  • the feedforward control law generally used for the high-precision attitude maneuver of an artificial satellite is not less than a specific frequency as shown in Fig. 1 (b) and Fig. 3.
  • the characteristic of the control target quantity (angular acceleration profile) as a feed-forward control law determined in advance as shown in the upper part of Fig. 4 is generated. Therefore, it is characterized by the fact that a posture maneuver that suppresses the vibration of a flexible object as much as possible is possible, as shown in the mode characteristics during posture maneuver shown in Figs. 5 (a) and (b).
  • FIG. 1 (a) has two positive and negative polarities that are shifted by one period and superimposed, and frequency separation as shown in Fig. 1 (b) and Fig. 3 is performed.
  • the posture can be accelerated and decelerated.
  • the window function as shown in Fig. 2 (b) is multiplied by the sampling function to make both ends of the waveform from zero to the end point (deceleration end) from the start point (start acceleration) to the end point (deceleration end).
  • the forced solution (residual vibration) after the end of the attitude maneuver can be made zero, as shown in the mode characteristics for the attitude maneuver.
  • the frequency component above a specific frequency is Using a sampling function that does not have a sinc function, the control target amount that is the attitude control data during attitude maneuver is generated as a pre-frequency shaping type feed-forward control law. Even if posture maneuver is performed on the target object, vibration can be effectively suppressed (this is an effective excitation avoidance technique for frequencies above a certain frequency), and the total time required for posture maneuver is unprecedented. It becomes possible to shorten it significantly. In the attitude control device to which such a technique is applied, it is possible to quickly and accurately perform attitude maneuvers that use reaction wheels, control moment jacks, etc., with a simple configuration, and especially for large flexible objects. This is very effective when high-precision attitude maneuvers are performed on satellites that have this. Industrial applicability
  • the attitude control data generation method of the present invention and the attitude control apparatus to which the attitude control data generation method is applied can be applied to the attitude maneuver of an artificial satellite that is a movable object on which a flexible object is mounted, and controls a robot manipulator that is a movable object. Application in usage is also possible.

Abstract

 柔軟物(50)を有する可動物体の姿勢を姿勢マヌーバにより制御するため、まず、例えば姿勢マヌーバ時の柔軟物の振動に基き、特定の周波数以上の周波数成分を持たないサンプリング関数を得る。このサンプリング関数を用いて制御目標量を事前周波数成形型フィードフォワード制御則として生成する。この制御目標量に基いて姿勢制御データを生成する。この姿勢制御データは可動物体の姿勢マヌーバに利用することができる。

Description

姿勢制御データ生成方法、 及びそれを適用した姿勢制御装置 技術分野
本発明は、 飛翔体等の可動物体を対象とする姿勢マヌーバ時における姿勢制御 データを生成する姿勢制御データ生成方法、 及びそれを適用した姿勢制御装置に 関する。 糸田
背景技術
一般に、 柔軟物を有する人工衛星等の飛翔体の姿勢変更制御においては、 多く の場合、 より早く、 より高精度に所望の姿勢へと変更を行うことが要求される。 この際、 最も問題視されるのが柔軟物の振動である。 柔軟物が姿勢変更によって 揺さぶられると、 飛翔体の本体にその振動が伝わって本体の指向方向に誤差が生 .じてしまう。 また、 一旦、 柔軟物に振動が発生すると、 その振動が減衰するまで には相当な時間がかかってしまう。 このように、 柔軟物を有する飛翔体の姿勢変 更制御については、 柔軟物の振動を如何に抑制できるかが重要な課題となってい る。
例えば日本国特許第 3 6 2 3 7 4 7号公報 (文献 1 ) は、 人工衛星等に対する 姿勢変更 (姿勢制御) に関連する技術として、 例えば三軸衛星の姿勢変更制御方 式を開示している。 この姿勢変更制御方式では、 姿勢制御目標量の生成過程で、 目標量に含まれる周波数成分については、周波数分離を行っていない。このため、 人工衛星等に搭載された柔軟物 (人工衛星柔軟変動体部) が励振することを回避 できず、 マヌーバ終了後にも柔軟物の定常的な振動が残留する。 即ち、 マヌーバ 終了後、 観測ミッションを開始できる状態に至るまでに姿勢静定 (振動収束) の 待機時間が必要である。 したがって、 振動が収束するまでは指向安定度が確保さ れないため、 観測ミッションを開始することができない。 結果として、 姿勢マヌ ーバに要する全体の時間(姿勢マヌーパ開始から観測ミッション開始までの時間) が長く掛かってしまう。
そこで、 姿勢変更制御に際して柔軟物の振動を抑制するための対策も提案され ている。 例えば、 N. C. S i n g e r等により提案された B a n g— B a n g 制御による柔軟衛星の姿勢マヌーバ制御手法(文献 2:「S i n g e r, N. C., S e e r i n g , W. P. : P r e s h a p i n g Co mm a n d I n p u t s t o R e d u c e S y s t em V i b r a t i o n, J o u r n a l o f Dyn am i c S y s t ems a n d Me a s u r eme n t C o n t r o l , 1 12 (1 990), p p. 76— 82」 参照)、 及ぴ B . W i e等により提案された B a n g-B a n g制御による柔軟衛星の姿勢マヌーバ 制御手法 (文献 3 :「L i u, Q., W i e, B. : R o b u s t T i me— Op t i m a 1 C o n t r o l o f Un c e r t a i n F l e x i b l e S p a c e c r a f t, J o u r n a l o f Gu i d a n c e, C o n t r o 1 , a n d Dyn am i c s, 15 (1 992), p p. 597— 604」 参 照) 等がある。 発明の開示
上述した文献 2または文献 3の B a n g— B a n g制御は、 何れも制御ァクチ ユエータとしてスラスタを使用対象としたものであり、 スラスタ以外の制御ァク チユエータ、例えば人工衛星の姿勢を回転させるためのリアクションホイール(R e a c t i o n Wh e e l : RW) または人工衛星の姿勢を動かすためのコン トローノレモーメントジャィ口 (C o n t r o l Mome n t Gy r o : CM G) 等へは適用し難い。
それらの B a n g-B a n g制御は、 柔軟物が 1つの柔軟モードのみを有する 場合には有効であると言える。しかし、実際の柔軟物は分布定数系であると共に、 複数の支配的な柔軟モードを有する振動系である場合が殆どであるため、 或る特 定の柔軟モードの励振を回避できるとしても、 その他の大部分の大きな柔軟モー ドを回避することはできない。
それ故に本発明の目的は、 柔軟物を有する可動物体の姿勢マヌーバに要する全 体の時間の短縮を可能にする姿勢制御データ生成方法を提供することにある。 本発明の他の目的は、 可動物体の柔軟物が複数の柔軟モードを有する場合でも 姿勢マヌーバにおける振動の効果的な抑制を可能にする姿勢制御データ生成方法 を提供することにある。
本発明の別の目的は、 柔軟物を有する可動物体の姿勢マヌーバに要する全体の 時間を短縮できる姿勢制御装置を提供することにある。
本発明の他の目的は、 可動物体の柔軟物が複数の柔軟モードを有する場合でも 姿勢マヌーパによる振動を効果的に抑制することができる姿勢制御装置を提供す ることにある。
本発明の第 1の態様によれば、 柔軟物を有する可動物体の姿勢マヌーバにおけ る姿勢制御データを生成する姿勢制御データ生成方法であって、 前記姿勢マヌー バ時の前記柔軟物の振動に基き、 特定の周波数以上の周波数成分を持たないサン プリング関数を得るサンプリング関数取得段階と、 前記サンプリング関数を用い て制御目標量を事前周波数成形型フィードフォヮ一ド制御則として生成する制御 目標量生成段階と、 前記制御目標量に基いて前記姿勢制御データを生成する制御 データ生成段階とを有することを特徴とする方法が得られる。
前記制御目標量生成段階では、 前記サンプリング関数として、 周波数成分を完 全分離した姿勢プロファイルの波形を生成可能な s i n c関数を用いて前記制御 目標量を生成してもよい。
前記制御目標量生成段階では、 前記 s i n c関数による前記サンプリング関数 における正負の 2つの波形を 1周期分ずらして組み合わせることにより前記制御 目標量として加速減速用データを生成してもよい。
前記制御目標量生成段階では、 前記姿勢マヌーバ後の前記柔軟物の振動を最大 限に抑制するために、 前記 s i n c関数に対して窓関数を掛け合わせて前記制御 目標量を生成してもよい。
本発明の第 2の態様によれば、 柔軟物を有する可動物体の姿勢を姿勢マヌーパ により制御する姿勢制御装置であって、 特定の周波数以上の周波数成分を持たな いサンプリング関数を用いて制御目標量を事前周波数成形型フィードフォヮ一ド 制御則として生成する制御目標量生成手段を有することを特徴とする装置が得ら れる。
前記制御目標量生成手段は、 前記姿勢マヌーバ時の前記柔軟物の振動に基いて 前記サンプリング関数を得るものであってもよい。
さらに、 前記制御目標量に基いて生成した姿勢制御データを用いて前記可動物 体の姿勢を変更する姿勢変更手段を有してもよい。
さらに、 前記可動物体に設置された、 前記可動物体を回転させるための姿勢回 転手段を含むものでもよい。
前記制御目標量生成手段は、 前記サンプリング関数として、 周波数成分を完全 分離した姿勢プロファイルの波形を生成可能な s i n c関数を用いて前記制御目 標量を生成するものであってもよい。
前記制御目標量生成手段は、 前記 s i n c関数による前記サンプリング関数に おける正負の 2つの波形を 1周期分ずらして組み合わせることにより前記制御目 標量として加速減速用データを生成するものであってもよい。
前記制御目標量生成手段は、 前記 s i n c関数に対して窓関数を掛け合わせて 前記制御目標量を生成するものであってもよい。
前記可動物体は人工衛星であり、 前記柔軟物は前記人工衛星に搭載された柔軟 変動体部であってもよい。
前記柔軟変動体部は、 前記人工衛星の姿勢を回転させるためのリアクションホ ィールと、 前記人工衛星の姿勢を動かすためのコントロールモーメントジャイロ とを含むものでもよい。
本発明の第 1の態様による姿勢制御データ生成方法は、 柔軟物を有する可動物 体の姿勢マヌーバに要する全体の時間の短縮を可能にし、 また、 可動物体の柔軟 物が複数の柔軟モードを有する場合でも姿勢マヌーバにおける振動の効果的な抑 制を可能にする。
本発明の第 2の態様による姿勢制御装置は、 柔軟物を有する可動物体の姿勢マ ヌーバに要する全体の時間を短縮でき、 また、 可動物体の柔軟物が複数の柔軟モ ードを有する場合でも姿勢マヌーバによる振動を効果的に抑制することができる。 図面の簡単な説明
図 1は、 本発明の一実施例に係る姿勢制御データ生成方法で用いられるサンプ リング関数を説明するためのタイムチャートであり、 (a ) は 1周期分ずらした正 負逆極性の 2つのサンプリング関数の特性を示し、 ( b ) は ( a ) に示す 2つのサ ンプリング関数の特性を加算合成した合成特性を示す。
図 2は、 本発明の一実施例に係る姿勢制御データ生成方法で用いられるサンプ リング関数の処理を説明するためのタイムチャートであり、 (a ) は図 1 ( b ) に 示す合成特性を示し, (b ) は (a ) に示す合成特性におけるピーク部分に対する 窓関数の特性を示し, (c ) は (a ) に示す合成特性に (b ) の窓関数の特性を掛 け合わせた特性を示す。
図 3は、 図 2 ( c ) に示す窓関数の特性におけるピーク部分をその成分特性に おけるピーク部分と対比して周波数に対する振幅の関係で示すタイムチャートで ある。
図 4は、 本発明の一実施例に係る姿勢制御データ生成方法でサンプリング関数 に基いて生成される制御目標量であるマヌーバプロフアイルの諸特性を例示した タイムチャートである。
図 5は、 本発明の一実施例に係る姿勢制御装置により柔軟物 (人工衛星柔軟変 動体部) を対象として実行する姿勢マヌーバ時の柔軟モード変位の特性を例示し たタイムチャートであり、 (a )は時間に対する柔軟モード変位の関係で示される 基本モード特性を示し, (b ) は時間に対する微分柔軟モード変位の関係で示され る微分モード特性を示す。
図 6は、 本発明の一実施例に係る姿勢制御データ生成方法を適用した人工衛星 用姿勢制御装置の基本構成を示した機能ブロック図である。 発明を実施するための最良の形態
先ず、 図 1一 5を参照して、 本発明の実施の形態に係る姿勢制御データ生成方 法とそれを適用した姿勢制御装置の概要とについて説明する。 この姿勢制御データ生成方法は、 柔軟物を有する可動物体の姿勢マヌーバ時に おける姿勢制御データを生成する方法であって、 姿勢マヌーパ時の柔軟物の振動 に基き、 特定の周波数以上の周波数成分を持たないサンプリング関数を得るサン プリング関数取得段階と、 サンプリング関数を用いて制御目標量を事前周波数成 形型フイードフォヮ一ド制御則として生成する制御目標量生成段階とを有してい る。 さらに、 この制御目標量に基いて制御データ生成段階で姿勢制御データを生 成する。
図 1は、制御目標量生成段階で用いられるサンプリング関数の特性 [時間 (s) に対する振幅の関係で表わされる] を例示したものであり、 同図 (a) は 1周期 分ずらした正負逆極性の 2つのサンプリング関数の特性 ( s— f 1 1 , s - f 1 2) に関するもの, 同図 (b) は同図 (a) の特性を加算合成した合成特性 (s 一 f 1 1 + 1 2) に関するものである。
図 1 (a) では、 上述した制御目標量生成段階において、 サンプリング関数の 特性 s— f 1 1, s— f 1 2として、 姿勢プロファイルの波形を生成可能な 1周 期分ずらした正負逆極性の s i n c関数を用いる。 図 1 (b) では、 これらの s i n c関数によるサンプリング関数を加算合成することにより、 合成特性 (s— f 1 1 + 1 2) として制御目標量となる加速減速用データを生成した様子を示し ている。
図 2は、 制御目標量生成段階で用いられるサンプリング関数を応用した場合の 諸特性 [時間 (s) に対する振幅の関係で表わされる] を対比して例示し、 同図
(a) は図 1 (b) に示す合成特性 (s - f 1 1 + 1 2) に関するもの, 同図 (b) は同図 (a) の特性におけるピーク部分に対する窓関数の特性 (w— f ) に関す るもの, 同図 (c) は同図 (a) の特性に同図 (b) の特性を掛け合わせた特性
(s - f 1 1 + 1 2) X (w- f ) に関するものである。
ここでは、 更に制御目標量生成段階において、 姿勢マヌーバ終了後の柔軟物の 振動を最大限に抑制するために、 図 2 (a) に示されるような s i n c関数に基 づいて合成して得た合成特性 (s— f 1 1 + 1 2) に対して図 2 (b) に示され るようなハミング窓関数の特性 (w— f ) を掛け合わせた結果として得られる図 2 (c) に示されるような特性 (s— f 1 1 + 1 2) X (w- f ) として制御目 標量を生成することを示している。
図 3は、 図 2 (c) に示す特性 ( s— f 1 1 + 1 2) X (w- f ) におけるピ ーク部分をその成分特性 (s— f 1 1 + 1 2) におけるピーク部分と対比して周 波数 (Hz) に対する振幅の関係で示したものである。
ここでは、 上述した加算合成の特性 (s - f 1 1 + 1 2) に窓関数の特性 (w - f ) を掛け合わせても、 周波数成分が完全分離された様子、 即ち、 制御目標量 生成段階での制御目標量の生成において、 サンプリング関数として s i n c関数 を用いることにより周波数分離を完全に行うことができることを示している。 こ のことは、 姿勢マヌーバ後の柔軟物の定常的な残留振動を殆ど零に近付けること ができ、 各種の姿勢マヌーバを経て行われる観測ミッシヨンを最短時間の姿勢マ ヌーバにより効果的に実施することができることを示唆している。特に、図 1 ( b ) に示した合成特性 (s— f l l + 1 2) のようにサンプリング関数を重ね合わせ るだけでは、 姿勢マヌーバ後に振動が残ってしまうことになるため、 図 2 (b) に示したような窓関数の特性 (w- f ) を掛け合わせることによって、 制御目標 量を姿勢マヌーバから姿勢マヌーバ終了に及んで滑らかに減衰させることができ るので、 この点が重要になっている。
因みに、サンプリング関数(s i n e関数) とは、特定の周波数帯域 ω s (Hz) 以上の周波数成分を全く含まないことで一般に知られる関数である。 但し、 ここ で特定の周波数帯域 c s (H z) のサンプリング関数を数式で表わすと、 y (t) = s i η (ω s t ) Z w s tと る。
柔軟物を有する可動物体である人工衛星を最短時間で姿勢マヌーバを完了する ためには、 第 1の留意事項として、 定常的な柔軟物の振動を引き起こさない加減 速プロファイルを用いること、 第 2の留意事項として、 ァクチユエータ (トルク 発生装置) が発生できるトルク能力を最大限使うような姿勢マヌーバを行うこと が必要であり、 これらを実現できれば良いことになる。
このうち、 第 1の留意事項については、 図 1 (a) に示すようなサンプリング 関数 (s i n c関数) の波形を用いて人工衛星の姿勢マヌーバに際しての角加速 度プロファイル (制御目標量) の生成を行えば実現できる。 サンプリング関数を 利用して生成した角加速度プロフアイルを周波数分布で見ると、 図 3に示したよ うな分布となり、 サンプリング関数によつて決まる或る特定の周波数以上の周波 数成分を含まない波形であることが分かる。
又、 第 2の留意事項については、 図 1 ( a ) に示される 2つのサンプリング関 数のピークの山同士を 1周期分ずらして正の波を加速用, 負の波を減速用の角加 速度波形と考え、 2つの波を重ね合わせることで人工衛星の姿勢の加速減速の角 加速度プロファイルを図 2 ( c ) に示されるような特性として生成すれば実現で きる。 ここで、 最大限にァクチユエータの能力を使うために、 角加速度のピーク をァクチユエータ能力の最大トルクから求められる姿勢の角加速度に設定するこ とで、 第 1の留意事項の利点を活かしたままァクチユエータの能力を最大限利用 した加速減速を行うことができる。
尚、 サンプリング関数 (s i n e関数) を 1周期ずらして重ね合わせることに より加速減速のプロファイルを生成する際、 例えば図 2 ( a ) に示されるような 2つのサンプリング関数の重ね合わせた特性の場合、 そのままでは姿勢マヌーバ 終了後に僅かに波が残ってしまう。 このため、 姿勢マヌーバ終了後の姿勢の振動 を減らすために、 更に図 2 ( b ) に示されるような窓関数を掛け合わせることで 得られる図 2 ( c ) に示されるような波形を人工衛星の姿勢の加速減速の角加速 度プロファイルとして利用することが重要となる。 こうした窓関数を掛け合わせ ることにより、 姿勢マヌーバ終了後の姿勢に振動が残らず、 より精度の高い指向 性能が確保できることになる。
ところで、 実際に人工衛星等に搭載される柔軟物を対象として姿勢マヌーバを 行う場合、 或る二つの姿勢を最短経路で変更 (姿勢マヌーバ) するためには、 一 般的に二つの姿勢から唯一求められるオイラー軸と呼ばれる軸回りの姿勢回転を 行うことが最短経路であることが知られている。 そこで、 本発明で実行する姿勢 マヌーバ則は、 姿勢変更前の姿勢と姿勢変更後の目標姿勢とからオイラー軸を求 め、 そのオイラー軸回りに人工衛星の姿勢を回転させるような制御を行うもので ある。 図 4は、 制御目標量生成段階で生成するオイラー軸決定後の制御目標量である マヌーバプロファイルの諸特性を例示したものである。 但し、 図 4中の下段は時 間 (s) に対する角度 0 (d e g) の関係で示される角度特性に関するもの, 中 段は時間 (s) に対する角速度 d eZd t (d e g/s) の関係で示される角速 度特性に関するもの, 上段は時間 (s) に対する角加速度 d20Zd t 2 (d e g
/ S 2) の関係で示される角加速度特性に関するものである。
ここでは、観測ミツション周期 t。 b sと姿勢マヌーバ周期 t mn vとが交互に繰り 返されることを前提として、 下段に示されるような 3 d e g (3度) のオフセッ ト角 0。f f setを持つ波形の角度特性を時刻 (t) で微分して得られる中段に示さ れるような角速度特性について、 更に時刻 (t) で微分 [即ち、 角度特性を時刻 ( t )で 2階微分]して得られる上段に示されるような角加速度特性に基づいて、 オイラー軸回りに人工衛星の姿勢を回転させるような姿勢マヌーバを実行するこ とを示している。
一方、 姿勢制御データ生成方法を適用した姿勢制御装置は、 柔軟物に設置され た姿勢を回転させるための姿勢回転手段と姿勢を動かすための姿勢変更手段とへ 制御信号を送出して柔軟物を或る静止状態から別の静止状態へと姿勢マヌーバを 行う基本構成を含んでいる。 さらに、 この姿勢制御装置は、 姿勢マヌーバ時に柔 軟物により発生する振動を抑制するために制御信号の制御要素として含まれると 共に、 特定の周波数以上の周波数成分を持たないサンプリング関数を用いて姿勢 マヌーバ時の姿勢制御データとなる制御目標量を事前周波数成形型フィードフォ ヮード制御則として生成する姿勢マヌーバ時制御目標量生成手段を備えたもので ある。
ここで、 姿勢マヌーバ時制御目標量生成手段は、 サンプリング関数として、 周 波数成分を完全分離した姿勢プロファイルの波形を生成可能な s i n c関数を用 いて制御目標量を生成する機能を持つこと、 s i n c関数によるサンプリング関 数における正負の 2つの波形を 1周期分ずらして組み合わせることにより制御目 標量として加速減速用データを生成する機能を持つこと、 更に姿勢マヌーバ終了 後の柔軟物の振動を最大限に抑制するために s i n c関数に対して窓関数を掛け 合わせることで制御目標量を生成する機能を持つことがそれぞれ好ましい。 何れ にしても、 柔軟物については、 人工衛星に搭載された人工衛星用柔軟変動体部で ある場合には有効となる。 人工衛星用柔軟変動体部としては、 人工衛星の姿勢を 回転させるためのリアクションホイールや、 人工衛星の姿勢を動かすためのコン ト口ールモーメントジャィ口等が含まれることが好ましい。
この姿勢制御装置の場合においても、 人工衛星用柔軟変動体部を対象として或 る静止状態から別の静止状態への姿勢マヌーバを行う際、 姿勢マヌーバ時の制御 目標量の生成則として、 上述した図 1 (a), (b) 乃至図 4で説明したような特 定の周波数以上の周波数成分を持たないサンプリング関数 (s i n c関数) を利 用すれば、 姿勢マヌーバにより発生する柔軟物の振動を大きく低減することが可 肯 gとなる。
図 5は、 係る姿勢制御装置で人工衛星柔軟変動体部を対象として実行する姿勢 マヌーバ時の柔軟モード変位の特性を例示したもので、 同図 (a) は時間 (s) に対する柔軟モード変位 の関係で示される基本モード特性に関するもの, 同図 ( b ) は時間 ( a ) に対する微分柔軟モード変位 d / d tの関係で示される微 分モード特性に関するものである。
即ち、 図 5 (a), (b) は、 姿勢マヌーバ時の人工衛星柔軟変動体部について のモード変位の時間履歴を示している。 図 5 (a) に示される基本モード特性の 場合、 1回目モード (1 s t mo d e) では零の特性 S Iとなり、 2回目モー ド (2 s t mo d e) では 10秒〜 20秒の間で負, 正を向く大きな二つのピ ークを持つと共に、 40秒〜 50秒の間でも正, 負を向く大きな二つのピークを 持つた特性 S 2となり、 4回目モード ( 4 s t m o d e ) では 10秒〜 20秒 の間で負, 正を向く二つの小さなピークを持つと共に、 40秒〜 50秒の間でも 正, 負を向く小さな二つのピークを持った特性 S 2となり、 その他の回数のモー ド [3回目モード (3 s t mo d e), 5回目モード (5 s t mo d e) 〜8 回目モード (8 s t mo d e)] では殆ど零となっている様子を示している。 図 5 (b) に示される微分モード特性の場合、 1回目モード (1 s t mo d e) では零の特性 S 1となる力 2回目モード (2 s t mo d e) では 10秒 〜 20秒の間で負, 正, 負を向く中, 大, 中の三つのピークを持つと共に、 40 秒〜 50秒の間でも正, 負, 正を向く中, 大, 中の三つのピークを持った特性 S 2となり、 4回目モード( 4 s t mo d e)では 10秒〜 20秒の間で負, 正, 負を向く三つの非常に小さなピーク (正を向く中央のピークは他のものよりも大 きい) を持つと共に、 40秒〜 50秒の間でも正, 負, 正を向く非常に小さな三 つのピーク (負を向く中央のピークは他のものよりも大きい) を持った特性 S 2 となり、 その他の回数のモード [3回目モード (3 s t mo d e), 5回目モー ド (5 s t mo d e) 〜8回目モード (8 s t m o d e )] では殆ど零となつ ている様子を示している。
このように、 上述した図 2 (c) に示す特性 (s— f 1 1 + 12) X (w- f ) に基づいて生成した制御目標量を用いて姿勢マヌーバを行うと、 或る静止状態か ら別の静止状態へと姿勢マヌーバした後、 定常的な振動が殆ど残留せず、 直ちに 姿勢の振動状態が収束することになり、 姿勢マヌーバ直後に指向安定度が確保で きるようになる。 このため、 最短の所要時間で姿勢マヌーバを完了することがで き、 各種の姿勢マヌーバを経て行われる観測ミッシヨンを最短時間の姿勢マヌー バにより効果的に実施することが可能となる。
次に、 図 6をも参照して、 上述した姿勢制御データ生成方法を適用した人工衛 星用姿勢制御装置について説明する。
姿勢制御装置 (At t i t u d e C o n t r o l S y s t ern : ACS) 40は、 外部の他の図示されない電気系装置から妨害トルク (D i s t u r b a n c e To r qu e : DT) 信号 (抑制トルク信号と呼ばれても良い) が入力 されると共に、 人工衛星の姿勢を回転させるための略図するリアクションホイ一 ル(RW)、並びに人工衛星の姿勢を動かすための略図するコントロールモーメン トジャイロ (CMG) が設置された人工衛星用柔軟変動体部 50におけるリアク シヨンホイール (RW), コントロールモーメントジャイロ (CMG) へそれぞれ 制御信号て RW, て CMGを送出することにより図示されない人工衛星の本体を姿勢 制御するものである。 因みに、 リアクションホイール (RW) は、 フライホイ一 ルを回転させることにより、 その反力トルクを発生させることができ、 これによ つて人工衛星の姿勢を回転させることができるァクチユエータである。 コント口 ールモーメントジャイロ (CMG) は、 或る一定の回転スピードで回転させたフ ライホイールをフライホイールの回転軸に対して直交方向に回転軸を持つジンバ ルを用いて回転させることにより、 ジャイロ効果のトルクを出力でき、 これによ つて人工衛星の姿勢を回転制御することができるァクチユエータである。
姿勢制御装置 (ACS) 40におけるハードウェア部分は、 基本構成として、 人工衛星用柔軟変動体部 50における姿勢情報センサとして働くジャイロセンサ としてのジャイロスコープ (Gy r o S c o p e s : GS) 1と、 人工衛星用 柔軟変動体部 50における恒星,地球,太陽の角度検出機能を持つ光センサ群(O p t i c a l S e n s o r s : OS) 2と、 人工衛星柔軟変動部 50における コントロールモーメントジャイロ (CMG) のジンバノレについての角度 Z角速度 の比を検出するジンバル角度 Z角速度センサ (G imb a l An g 1 e/R a t e S e n s o r) 3と、 各センサからの情報を用いてリアクションホイール (RW), コントロールモーメントジャイロ (CMG) へ制御信号 T rw, τ CMG¾r 送出する姿勢制御回路部 (A t t i t u d e C o n t r o l E l e c t r o n i c s : ACE) 30とを含んでいる。
ジャイロスコープ (GS) 1は、 慣性空間 (ここでは宇宙空間と考えて良い) に対して人工衛星の直交三軸回りに回転した角速度 ωを検出した情報を含む角速 度信号 coGSを出力するセンサ (或いは角度 0を検出した情報を含む角度検出信号 0cSを出力することも可能) である。 光センサ群 (〇S) 2は、 恒星, 地球, 太 陽との角度を検出した情報を含む角度検出信号 eosを出力するセンサの集合体 である。 尚、 上述したリアクションホイール (RW), コントロールモーメントジ ャイロ (CMG), ジャイロスコープ (G S) 1, 光センサ群 (OS) 2について は、 人工衛星の姿勢制御系の分野では汎用的であり、 ここでの技術的要旨とは直 接関係しないため、 その詳細な構成や機能についての説明は省略する。
ジンバル角度/角速度センサ 3は、 コントロールモーメントジャイロ (CMG) のジンバルから検出した角度 Z角速度の比を示す角度 Z角速度比信号 0 ^/ω^ を出力するセンサである。 姿勢制御回路部 (ACE) 30は、 ジャイロスコープ (GS) 1からの角速度 信号 coGS, 光センサ群 (OS) 2からの角度検出信号 0 osに基づいて人工衛星の 本体姿勢制御に必要な角速度 ωΒの推定値を示す角速度推定信号, 角度 6 Βの推定 値を示す角度推定信号を出力する姿勢推定器 (At t i t u d e E s t i ma t a o r : AE) 4と、 姿勢推定器 (AE) 4からの角速度 ω Bの推定値及び角度 0Bの推定値に基づいて姿勢マヌーバを実行する前に初期の姿勢と変更後の姿勢 とから求めたオイラー軸に応じて姿勢マヌーバに必要な制御目標量を示す目標姿 勢プロファイルとしての角度 eMの目標値を示す角度目標信号, 角速度 ω Mの目標 値を示す角速度目標信号の生成出力、 姿勢の流れ (At t i t u d e D r i f t) の検出結果に基づいてコントロールモーメントジャイロ (CMG) のジンバ ル向けの角度 0eの目標値を示す角度目標信号, 角速度 coGの目標値を示す角速度 目標信号の生成出力、 並びにリアクションホイール (RW) に対するフィードフ ォヮード制御則の制御信号て FF (RW)の生成出力を担う姿勢誘導部(A t t i t u d e Gu i d a c e Law : AG L) 20と、 姿勢誘導部 (AGL) 20か らの角度目標信号, 角速度目標信号から姿勢推定器 (AE) 4からの角速度推定 信号, 角度推定信号をそれぞれ減算する減算器と、 この減算器での減算値を姿勢 の流れ (AD) を示す信号として入力し、 リアクションホイール (RW) に対す るフィードバック量の制御信号 τ FB (RW)を P I D制御により生成出力するリアク シヨンホイール制御部 (RW Co n t r o l l e r) 5と、 姿勢誘導部 (AG L) 20からの制御信号 τ FF ( とリアクションホイール制御部 5からの制御信 号 τΡΒ (RW) とを加算する加算器と、 この加算器での加算値を入力し、擬似反転さ せて制御信号て RWを生成出力するリアクションホイール分配器 (RW D i s t r i b u t o r) 10と、 姿勢誘導部 (AGL) 20からの角度目標信号, 角速 度目標信号からジンバル角度 Z角速度センサ 3からの角度ノ角速度比信号 Θ G/ coGを減算する減算器と、 この減算器での減算値を入力し、 PD制御によりコント ロールモーメントジャイロ (CMG) のジンバル向けの制御信号て CMGを生成出 力するコントロールモーメントジャイロ向けジンパル制御部 (CMG— G i mb a 1 Co n t r o l l e r) 1 1とを含んでいる。 更に、 姿勢誘導部 (AGL) 20は、 姿勢推定器 (AE) 4からの角速度 ωΒ の推定値及び角度 θ Βの推定値に基づいて姿勢マヌーバを実行する前に運動学上 で初期の姿勢と変更後の姿勢とからオイラー軸を求めると共に、 オイラー軸に応 じて姿勢マヌーバに必要な制御目標量を示す目標姿勢プロファイルとしての角度
0 Μの目標値を示す角度目標信号と角速度 ω Μの目標値を示す角速度目標信号と を出力する制御目標量生成手段としての姿勢マヌーバ時制御目標量プロファイル 生成部 (NME P r o f i l e r) 6と、 姿勢マヌーバ時制御目標量プロファ ィル生成部 6からの角度目標信号, 角速度目標信号を入力してコント口ールモー メントジャィ口 (CMG) に対するフィードフォヮ一ド制御則の制御信号 τ FF (c MG), 並びに反転変動によりリアクションホイール (RW) に対するフィードフォ ワード制御則の制御信号 τ ' FF(RW)を生成出力するトルクプロファイル生成/分 配器 (To r q u e P r o f i l e r /D i s t r i b u t o r) 7と、 リア クションホイール制御部 5の入力側からの姿勢の流れを示す信号 (点線で示す) を入力し、 その信号のバイアス流れを減少させることにより トルクプロフアイノレ 信号を生成出力するトルクプロファイル生成器 (To r q u e P r o f i l e r ) 9と、 トルクプロファイル生成器 9からのトルクプロファイル信号と トルク プロフアイル生成 Ζ分配器 7からの制御信号 τ F F (RW)とを加算してリアクショ ンホイール(RW) に対するフィードフォヮ一ド制御則の制御信号 τ FF (RW)を出 力する加算器と、 トルクプロファイル生成/分配器 7からの制御信号 T FF (CMG) に基づいてコントロールモーメントジャイロ (CMG) のジンバル向けの角度 0G の目標値を示す角度目標信号,角速度 ω Gの目標値を示す角速度目標信号を生成出 力するコントロールモーメントジャイロ向けジンバル誘導部 (CMG-G i mb a 1 S t e e r i n g L aw) 8とを含んでいる。
尚、 姿勢制御回路部 (ACE) 30の各部については、 ハードワイア一ドロジ ック回路のみで構成することも可能である他、 CPU等のコンピュータを内蔵し て制御ロジックをソフトゥエアで実行する構成とすることもできる。 ここでは、 姿勢制御回路部 (ACE) 30における姿勢誘導部 (AGL) 20内では姿勢マ ヌーバ時制御目標量プロフアイル生成部 6, トルクプロフアイル生成/分配器 7, コントロールモーメントジャィ口向けジンノ ル誘導部 8の各部、 姿勢誘導部 (A G L) 2 0外では姿勢推定器 (A E ) 4, リアクションホイール制御部 5 , リア クシヨンホイール分配器 1 0, コントロールモーメントジャイロ向けジンバル制 御部 1 1の各部がソフトウエアにより実現されているものとする。
ソフトウェアにより具現される各部の技術的機能を簡単に説明すれば、 姿勢推 定器 (A E ) 4の場合、 ジャイロスコープ (G S ) 1の出力, 光センサ群 (O S ) 2の出力から人工衛星用柔軟変動体部 5 0における現在姿勢を推定する機能を持 つ。 姿勢マヌーバ時制御目標量プロファイル生成部 6は、 姿勢マヌーバ用の事前 周波数整形型フイードフォヮ一ド制御則を実行する部分であり、 そのための制御 目標量 (目標姿勢) のプロファイルを生成する。 トルクプロファイル生成 Z分配 器 7は、 目標姿勢のプロファイルを実現するために必要な制御トルクを計算する 機能を担う他、 衛星が角運動量を持っていることによる軸間干渉量を回避するた めに必要なデカップリング量も同時に計算すると共に、 分配機能では計算で得ら れた制御トルクをコントロールモーメントジャイロ ( CMG) 向けとリアクショ ンホイール (RW) 向けとに分配する機能を持つ。
コントローノレモーメントジャイロ向けジンノ ル誘導部 8は、 コントロールモー メントジャイロ ( CMG) 向けとして分配された制御トルクを出力するために必 要なコントロールモーメントジャイロ (CMG) のジンバル角プロファイル及び ジンバル角速度プロファイルを計算する。 コントロールモーメントジャイロ向け ジンパル制御部 1 1は、 ジンバル角度ノ角速度センサ 3から得られた現在のジン バル状態と目標のジンバノレ状態との対比で得られるジンバル角誤差及ぴジンバル 角速度誤差に基づいて、 ジンバルが目標の運動を行うためのマイナーループ制御 を行う。 リアクションホイール制御部 5は、 姿勢角誤差及ぴ姿勢角速度誤差に基 づいて姿勢が目標の運動を行うためのフィードパック制御トルクの計算を行う。 リアクションホイール分配器 1 0は、 姿勢が目標の運動を行うために必要なフィ 一ドパック制御トルクを各リアクションホイールに分配して寄与する。
以上の各部の構成部分は、 柔軟物の固有振動数を励起しない最適なフィードフ ォヮード制御則を算出するための姿勢マヌーバ時制御目標量プロファイル生成部 6の構成部分を除けば、 人工衛星の姿勢制御に頻繁に用いられている周知構成を 導入したものである。
そこで、 姿勢マヌーパ時制御目檩量プロフアイル生成部 6の動作について詳細 に説明すれば、 ここでは姿勢変更前の姿勢と姿勢変更後の目標姿勢とから求めら れたオイラー軸の回りに図 4に示した角加速度プロファイル (上段の特性) を用 いるものであり、 これによつて人工衛星の姿勢を回転させるような姿勢マヌーバ を実行する。 即ち、 姿勢マヌーバ時制御目標量プロファイル生成部 6は、 姿勢マ ヌーバを実行する前に初期の姿勢と変更後の姿勢とからオイラー軸を求め、 この オイラー軸に応じて図 4の上段に示したような最短経路で姿勢マヌーバを実現で きるオイラー軸回りの人工衛星の角加速度の変化を計算して得た特性を得ており、 姿勢マヌーバ中は逐一、 このような特性の人工衛星の姿勢制御目標量を生成して コントロールモーメントジャイロ ( CMG) による姿勢誘導則を実行する姿勢誘 導部 (A G L ) 2 0、 及び P I D制御ループを成す Vアクションホイール制御部 5へ与える。 この目標量通りに人工衛星が回転すれば、 最短経路の姿勢変更制御 が柔軟物である人工衛星用柔軟変動体部 5 0におけるリアクションホイール (R W) , コントロールモーメントジャイロ (CMG) の振動を最低限に抑えながら実 現できるようになる。
以上の諸機能の組み合わせにより、 人工衛星は理論的に求められる最短な経路 を通って柔軟物にあっての最低限の振動しか励起させず、 これによつて高精度な 姿勢マヌーバを行うことができる。
尚、 図 6に示した姿勢制御装置 4 0により姿勢マヌーバを行う人工衛星用柔軟 変動体部 5 0におけるァクチユエータ部分 [リアクションホイ一ル (RW) , コン- トロールモーメントジャィ口 ( CMG) ] をリアクションホイール (RW) ゃスラ スタに置き換え、 コント口ールモーメントジャィ口向けジンバル誘導部 8をリア クションホイール- (RW) 向けゃスラスタ向けのものに代えてシステム構成する ことも可能であり、こうした場合にも同様な効果を期待できる。但し、この場合、 一般にリアクションホイール (RW) はコントロールモーメントジャィ口 ( CM G) よりも最大出力トルクは小さいものの、 小さい制御トルクをより精密に発生 させることができるため、 係る構成は小さい角度の姿勢変更をより高精度に行う ことが要求される場合に適している。 又、 一般にスラスタはリアクションホイ一 ル (RW) よりも大きく、 コントロールモーメントジャイロ (CMG) よりも小 さい制御トルクを発生するが、 リアクションホイール (RW) やコントロールモ —メントジャイロ (CMG) のように滑らかな出力トルク特性でないため、 係る 構成は一般に姿勢精度が若干粗くなる。 このように、 姿勢制御装置 40では、 若 干の変更を行うだけでァクチユエータの種類を限定することなく使用することが でき、 それぞれのァクチユエータの性能を有効に活用した姿勢マヌーバを行うこ とができる。
要するに、 上述した姿勢制御装置 40の場合、 一般に人工衛星の高精度な姿勢 マヌーバのために用いられるフィードフォワード制御則について、 図 1 (b) 及 ぴ図 3に示されるような特定の周波数以上の周波数を持たないサンプリング関数 (s i n c関数) の波形成形したものを用いて図 4の上段に示されるような事前 に決めるフィードフォヮ一ド制御則としての制御目標量の特性 (角加速度プロフ アイル) を生成しているため、 図 5 (a), (b) に示される姿勢マヌーバ時のモ 一ド特性のように柔軟物の振動を極力抑制した姿勢マヌーバが可能となることを 特色としている。 その他、 図 1 (a) に示されるサンプリング関数 (s i n c関 数) の正負の 2つの極性の波を 1周期分ずらして重ね合わせ、 図 1 (b) 及ぴ図 3に示されるような周波数分離効果を活かしたままの角加速度プロフアイルを利 用することにより、姿勢の加速と減速とが可能になる。又、サンプリング関数 (s i n c関数) のみを用いた場合には定常的な柔軟物の振動は励起しないものの、 強制力により強制解として発生する振動成分が姿勢マヌーバ終了後にも残ってし まうので、 これについては姿勢マヌーバの始点 (加速開始) から終点 (減速終了) に及んで波形の両端を零にするために図 2 (b) に示されるような窓関数をサン プリング関数に掛け合わせることにより、 図 5 (a), (b) に示される姿勢マヌ ーバ時のモード特性のように姿勢マヌーバ終了後の強制解 (残留振動) を零にす ることが可能になる。
上述した姿勢制御データ生成方法によると、 特定の周波数以上の周波数成分を 持たないサンプリング関数 (s i n c関数) を用いて姿勢マヌーバ時の姿勢制御 データとなる制御目標量を事前周波数成形型フィードフォヮ一ド制御則として生 成するため、 多数の柔軟モードを有する一般的な柔軟物を使用対象にして姿勢マ ヌーパを行っても有効に振動を抑制することができ (特定の周波数以上の周波数 に有効な励振回避手法となる)、従来に無く姿勢マヌーパに要する全体の時間を顕 著に短縮できるようになる。 係る手法を適用した姿勢制御装置では、 簡素な構成 でリアクションホイールやコントロールモーメントジャィ口等を使用対象にした 姿勢マヌーバを迅速にして、 精度良く適確に行うことが可能となり、 特に大きな 柔軟物を有する人工衛星を対象として高精度な姿勢マヌーパを行う際に非常に有 効となる。 産業上の利用可能性
本発明の姿勢制御データ生成方法及ぴそれを適用した姿勢制御装置は、 柔軟物 が搭載された可動物体である人工衛星の姿勢マヌーバへの適用できる他、 可動物 体であるロボットマニピュレータを制御する用途での適用も可能である。

Claims

請 求 の 範 囲
1 . 柔軟物を有する可動物体の姿勢マヌーバにおける姿勢制御データを生成す る姿勢制御データ生成方法であって、 前記姿勢マヌーバ時の前記柔軟物の振動に 基き、 特定の周波数以上の周波数成分を持たないサンプリング関数を得るサンプ リング関数取得段階と、 前記サンプリング関数を用いて制御目標量を事前周波数 成形型フィードフォヮ一ド制御則として生成する制御目標量生成段階と、 前記制 御目標量に基いて前記姿勢制御データを生成する制御データ生成段階とを有する ことを特徴とする方法。
2 . 請求の範囲 1に記載の姿勢制御データ生成方法であって、 前記制御目標量 生成段階では、 前記サンプリング関数として、 周波数成分を完全分離した姿勢プ 口ファイルの波形を生成可能な s i n c関数を用いて前記制御目標量を生成する ことを特徴とする方法。
3 . 請求の範囲 2に記載の姿勢制御データ生成方法であって、 前記制御目標量 生成段階では、 前記 s i n c関数による前記サンプリング関数における正負の 2 つの波形を 1周期分ずらして組み合わせることにより前記制御目標量として加速 減速用データを生成することを特徴とする方法。
4 . 請求の範囲 2又は 3に記載の姿勢制御データ生成方法において、 前記制御 目標量生成段階では、 前記姿勢マヌーバ後の前記柔軟物の振動を最大限に抑制す るために、 前記 s i n c関数に対して窓関数を掛け合わせて前記制御目標量を生 成することを特徴とする方法。
5 . 柔軟物を有する可動物体の姿勢を姿勢マヌーバにより制御する姿勢制御装 置であって、 特定の周波数以上の周波数成分を持たないサンプリング関数を用い て制御目標量を事前周波数成形型フイードフォヮ一ド制御則として生成する制御 目標量生成手段を有することを特徴とする装置。
6 .請求の範囲 5に記載の姿勢制御装置であって、前記制御目標量生成手段は、 前記姿勢マヌーバ時の前記柔軟物の振動に基いて前記サンプリング関数を得るも のであることを特徴とする装置。
7 . 請求の範囲 5又は 6に記載の姿勢制御装置であって、 さらに、 前記制御目 標量に基いて生成した姿勢制御データを用いて前記可動物体の姿勢を変更する姿 勢変更手段を有することを特徴とする装置。
8 . 請求の範囲 7に記載の姿勢制御装置であって、 さらに、 前記可動物体に設 置された、 前記可動物体を回転させるための姿勢回転手段を含むことを特徴とす る装置。
9 . 請求の範囲 5— 8のいずれかに記載の姿勢制御装置であって、 前記制御目 標量生成手段は、 前記サンプリング関数として、 周波数成分を完全分離した姿勢 プロファイルの波形を生成可能な s i n c関数を用いて前記制御目標量を生成す ることを特徴とする装置。
1 0 . 請求の範囲 9に記載の姿勢制御装置であって、 前記制御目標量生成手段 は、 前記 s i n c関数による前記サンプリング関数における正負の 2つの波形を 1周期分ずらして組み合わせることにより前記制御目標量として加速減速用デー タを生成することを特徴とする装置。
1 1 . 請求の範囲 9又は 1 0に記載の姿勢制御装置であって、 前記制御目標量 生成手段は、 前記 s i n c関数に対して窓関数を掛け合わせて前記制御目標量を 生成することを特徴とする装置。
1 2 . 請求の範囲 5— 1 1のいずれかに記載の姿勢制御装置であって、 前記可 動物体は人工衛星であり、 前記柔軟物は前記人工衛星に搭載された柔軟変動体部 であることを特徴とする装置。
1 3 . 請求の範囲 1 2に記載の姿勢制御装置であって、 前記柔軟変動体部は、 前記人工衛星の姿勢を回転させるためのリアクションホイールと、 前記人工衛星 の姿勢を動かすためのコントロールモーメントジャイロとを含むことを特徴とす る装置。
PCT/JP2007/059809 2006-05-12 2007-05-02 姿勢制御データ生成方法、及びそれを適用した姿勢制御装置 WO2007132793A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP07743244A EP2022718B1 (en) 2006-05-12 2007-05-02 Attitude control data creating method, and attitude control device applying the method
CA2651814A CA2651814C (en) 2006-05-12 2007-05-02 Attitude control data creating method, and attitude control device applying the method
US12/300,449 US8186627B2 (en) 2006-05-12 2007-05-02 Attitude control data creating method, and attitude control system applying the method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2006-133696 2006-05-12
JP2006133696A JP4189604B2 (ja) 2006-05-12 2006-05-12 柔軟構造物用姿勢マヌーバ時における姿勢制御データ生成方法、及びそれを適用した姿勢制御装置

Publications (1)

Publication Number Publication Date
WO2007132793A1 true WO2007132793A1 (ja) 2007-11-22

Family

ID=38693884

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2007/059809 WO2007132793A1 (ja) 2006-05-12 2007-05-02 姿勢制御データ生成方法、及びそれを適用した姿勢制御装置

Country Status (5)

Country Link
US (1) US8186627B2 (ja)
EP (1) EP2022718B1 (ja)
JP (1) JP4189604B2 (ja)
CA (1) CA2651814C (ja)
WO (1) WO2007132793A1 (ja)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101008176B1 (ko) * 2008-12-30 2011-01-13 한국항공우주연구원 반작용휠과 추력기 기반 자세제어기를 동시에 이용한 자세기동 및 가제어성 향상 방법
US8209065B2 (en) * 2009-01-21 2012-06-26 The Boeing Company Steering logic for spacecraft slew maneuvers
KR101319714B1 (ko) 2012-01-19 2013-10-17 한양대학교 산학협력단 리액션휠을 이용한 진동제거장치
JP5955037B2 (ja) * 2012-03-07 2016-07-20 三菱電機株式会社 アクチュエータ駆動装置
CN103808323B (zh) * 2012-11-07 2016-09-07 上海航天控制工程研究所 一种卫星姿态跟踪机动的余弦过渡角加速度路径方法
CN104249816B (zh) * 2013-06-27 2017-12-19 上海新跃仪表厂 非合作目标绕飞悬停的姿轨协同控制方法
CN104267732B (zh) * 2014-09-29 2017-07-28 哈尔滨工业大学 基于频域分析的挠性卫星高稳定度姿态控制方法
US9764858B2 (en) * 2015-01-07 2017-09-19 Mitsubishi Electric Research Laboratories, Inc. Model predictive control of spacecraft
CN105045270B (zh) * 2015-06-17 2017-09-29 北京控制工程研究所 一种基于振动补偿与状态反馈的刚柔系统姿态控制方法
CN105371851B (zh) * 2015-11-20 2018-08-07 国家测绘地理信息局卫星测绘应用中心 一种基于频域分析的卫星姿态模型构建方法
CN105786036B (zh) * 2016-04-05 2018-08-31 北京控制工程研究所 一种抑制转子动不平衡扰动的控制力矩陀螺框架控制系统及方法
JP6556371B2 (ja) * 2016-09-09 2019-08-07 三菱電機株式会社 姿勢制御装置、姿勢制御システム、地上局、人工衛星、姿勢制御方法及びプログラム
JP7199071B2 (ja) 2017-07-06 2023-01-05 国立研究開発法人宇宙航空研究開発機構 移動型撮像装置
CN107515612B (zh) * 2017-10-20 2021-03-02 湖北航天技术研究院总体设计所 基于侧喷流控制的弹性振动抑制方法
CN107632612B (zh) * 2017-10-30 2024-04-26 华中科技大学 一种内外联合式姿态控制水下航行器
CN109655218B (zh) * 2019-01-08 2020-10-09 上海卫星工程研究所 用卫星陀螺数据辨识整星挠性振动模态频率的方法及系统
US20220402633A1 (en) * 2020-02-26 2022-12-22 Mitsubishi Electric Corporation Orbital attitude control device, satellite, orbital attitude control method, and recording medium
CN113391638A (zh) * 2021-05-18 2021-09-14 北京航空航天大学 一种用于绳系拖曳离轨的系统组合体姿态协同控制方法
CN114229039B (zh) * 2021-12-14 2023-09-15 内蒙古工业大学 一种充液柔性航天器自适应无角速度复合控制方法
CN114578713B (zh) * 2022-03-17 2022-09-13 山东拓新电气有限公司 一种顶管机姿态控制方法及其装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11129997A (ja) * 1997-11-04 1999-05-18 Mitsubishi Electric Corp 人工衛星の姿勢制御装置
JP2001063699A (ja) * 1999-08-30 2001-03-13 Mitsubishi Electric Corp 人工衛星の姿勢制御装置
JP2001097291A (ja) * 1999-07-29 2001-04-10 Space Syst Loral Inc 適応性質量放出姿勢制御装置及び方法
JP2002255098A (ja) * 2001-03-01 2002-09-11 Nec Corp 人工衛星のマヌーバ制御装置
JP3623747B2 (ja) 2001-03-19 2005-02-23 日本電気航空宇宙システム株式会社 三軸衛星の姿勢変更制御方式
JP2006008132A (ja) * 2005-09-22 2006-01-12 Mitsubishi Electric Corp 人工衛星の姿勢制御装置

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0046151B1 (de) * 1980-08-19 1985-11-21 Messerschmitt-Bölkow-Blohm Gesellschaft mit beschränkter Haftung Vorrichtung zur Lagestabilisierung von elastischen Fahrzeugen
US4687161A (en) * 1985-09-30 1987-08-18 Ford Aerospace & Communications Corporation Pointing compensation system for spacecraft instruments
US4892273A (en) * 1988-12-30 1990-01-09 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Active damping of spacecraft structural appendage vibrations
US5311435A (en) * 1991-11-27 1994-05-10 Hughes Aircraft Company Method for station keeping control of flexible spacecraft using onboard gain scheduling scheme
US5257802A (en) * 1992-06-25 1993-11-02 The United States Of America As Represented By The Secretary Of The Air Force Zero-residual-energy minimum-time slew of a flexible space structure with damping
US5806804A (en) * 1996-11-12 1998-09-15 Lockheed Martin Corp. Adaptive harmonic disturbance compensation system
US5957411A (en) * 1997-01-31 1999-09-28 Space Systems/Loral, Inc. Method using double thruster firings to deadbeat flexible solar array structural oscillations
US6125310A (en) * 1997-07-18 2000-09-26 Space Systems/Loral, Inc. Thruster on time signaling with flexure compensation avoidance
JP3185738B2 (ja) * 1997-12-25 2001-07-11 日本電気株式会社 移動物体の状態制御装置及びその状態制御方法
US6347262B1 (en) * 2000-01-05 2002-02-12 Hughes Electronics Corporation Minimum fuel attitude and nutation controller for spinning spacecraft
US6845951B2 (en) * 2002-04-08 2005-01-25 The Aerospace Corporation Method of controlling pointing for a flexible structure
US7917256B2 (en) * 2007-05-02 2011-03-29 Lockheed Martin Corporation Active vibration damping (AVD) system for precision pointing spacecraft

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11129997A (ja) * 1997-11-04 1999-05-18 Mitsubishi Electric Corp 人工衛星の姿勢制御装置
JP2001097291A (ja) * 1999-07-29 2001-04-10 Space Syst Loral Inc 適応性質量放出姿勢制御装置及び方法
JP2001063699A (ja) * 1999-08-30 2001-03-13 Mitsubishi Electric Corp 人工衛星の姿勢制御装置
JP2002255098A (ja) * 2001-03-01 2002-09-11 Nec Corp 人工衛星のマヌーバ制御装置
JP3623747B2 (ja) 2001-03-19 2005-02-23 日本電気航空宇宙システム株式会社 三軸衛星の姿勢変更制御方式
JP2006008132A (ja) * 2005-09-22 2006-01-12 Mitsubishi Electric Corp 人工衛星の姿勢制御装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LIU, Q.; WIE, B.: "Robust Time-Optimal Control of Uncertain Flexible Spacecraft", JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS, vol. 15, 1992, pages 597 - 604, XP000287994
SINGER, N.C.; SEERING, W.P.: "Preshaping Command Inputs to Reduce System Vibration", JOURNAL OF DYNAMIC SYSTEMS AND MEASUREMENT CONTROL, vol. 112, 1990, pages 76 - 82
TOSHIAKI YAMASHITA ET AL.: "Improved Satellite Attitude Control Using A Disturbance Compensator", vol. 55, 2004, ACTA ASTRONAUTICA, pages: 15 - 25

Also Published As

Publication number Publication date
CA2651814C (en) 2013-02-19
US8186627B2 (en) 2012-05-29
EP2022718A1 (en) 2009-02-11
EP2022718A4 (en) 2010-07-07
JP2007302142A (ja) 2007-11-22
US20090218449A1 (en) 2009-09-03
JP4189604B2 (ja) 2008-12-03
EP2022718B1 (en) 2012-01-11
CA2651814A1 (en) 2007-11-22

Similar Documents

Publication Publication Date Title
WO2007132793A1 (ja) 姿勢制御データ生成方法、及びそれを適用した姿勢制御装置
JP5042906B2 (ja) 船舶用自動操舵装置
JP4679439B2 (ja) 人工衛星の姿勢制御装置
JP5993125B2 (ja) 船舶用自動操舵装置
JP3859454B2 (ja) 人工衛星のマヌーバ制御装置
US9505074B2 (en) Weaving control device for articulated robot
JP2004090796A (ja) 飛翔体の姿勢変更制御装置及び姿勢変更制御方法
US11669055B2 (en) Vibration suppression device, method and computer-readable medium using estimated vibration torque
Tomić Evaluation of acceleration-based disturbance observation for multicopter control
KR20140094313A (ko) 마찰력 보상기 및 이를 이용한 마찰력 보상방법
Goodarzi et al. Extended Kalman filter on SE (3) for geometric control of a quadrotor UAV
TW442425B (en) State control device of moving body and its state control method
Lee et al. Decentralized guidance and control for spacecraft formation flying using virtual leader configuration
US6860451B1 (en) Spacecraft spin axis reorientation method
JP5240365B2 (ja) 移動方向制御装置及びコンピュータプログラム
JP5955037B2 (ja) アクチュエータ駆動装置
Ireland et al. A comparison of inverse simulation-based fault detection in a simple robotic rover with a traditional model-based method
JP5126107B2 (ja) 人工衛星の姿勢制御装置
JP3905259B2 (ja) 人工衛星の姿勢制御装置
JPH11129997A (ja) 人工衛星の姿勢制御装置
JP2006008132A (ja) 人工衛星の姿勢制御装置
JP2001153680A (ja) 慣性装置
JP5963582B2 (ja) 船舶用自動操舵装置
JPS6138076B2 (ja)
JP2013056658A (ja) 船舶用自動操舵装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 07743244

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2651814

Country of ref document: CA

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2007743244

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 12300449

Country of ref document: US