Publication number | USRE39906 E1 |
Publication type | Grant |
Application number | US 09/888,291 |
Publication date | Nov 6, 2007 |
Filing date | Jun 21, 2001 |
Priority date | Oct 26, 1995 |
Fee status | Paid |
Also published as | USRE44396 |
Publication number | 09888291, 888291, US RE39906 E1, US RE39906E1, US-E1-RE39906, USRE39906 E1, USRE39906E1 |
Inventors | Gerald P. Roston, Charles J. Jacobus |
Original Assignee | Immersion Corporation |
Export Citation | BiBTeX, EndNote, RefMan |
Patent Citations (62), Non-Patent Citations (41), Referenced by (9), Classifications (28), Legal Events (4) | |
External Links: USPTO, USPTO Assignment, Espacenet | |
This is a Continuation application of U.S. Pat. No. RE37,374, filed on Nov. 30, 1999, in the name of Roston, et al., which is a reissue of U.S. Pat. No. 5,754,023 filed on Oct. 22, 1996, which claims priority to Provisional Application No. 60/005,861, filed Oct. 26, 1995.
The present invention relates generally to force feedback and, more particularly, to the use of gyroscopic stabilization to provide an inertial frame against which a force-reflecting device react.
Force-feedback technology and related devices may be divided into four broad application areas: medical, entertainment, teleoperations, and virtual reality. Teleoperations, the research of which provided the foundation for the development of force-feedback devices, is the process of locally controlling a remote device. The primary difference between virtual reality and teleoperations is in the objects which they control. With teleoperations, actual physical robots are manipulated in the real world, whereas virtual reality involves simulated devices in synthetic worlds. Force-feedback for telerobotics has evolved large and bulky mechanical arms to more joystick-like designs. In general, these devices are designed for six degree-of-freedom (6DOF) force feedback, and have the capability to provide high levels of force. More recently, finger-operated devices have also been introduced for use in teleoperations applications.
The use of force feedback in medical training, simulation, and teleoperations is also increasing, with the primary application being minimally invasive surgical techniques which use laparscopic tools to perform intricate tasks when inserted into body cavities through small incisions. To realistically simulate laparoscopic tool forces, special-purpose force-feedback devices are currently under development.
The entertainment field is very difficult to address with force-feedback technology, since the applications demand both higher performance and lower costs. There are three primary markets for force feedback devices in entertainment: location-based entertainment (LBE), arcades, and home entertainment. LBE demands the highest performance while home entertainment demands the lowest cost. Despite the conflicting demands, progress is being made in each of these fields.
It may be argued that each of the application domains just described has its roots in virtual reality, which is becoming dominant in all immersive applications. As a consequence, on-going research in immersive applications is often termed “virtual reality,” whereas, when the research is completed, the application is given a specific name, such as a surgical simulator. Overall, virtual reality is becoming increasingly popular as a preferred means of interacting with many scientific and engineering applications. To cite two of many examples, molecular modeling and automobile design are moving from standard graphics, carried out on conventional graphics terminals, to more interactive environments utilizing 3-D stereo graphics, head-mounted displays and force feedback.
As visualization is a very important aspect of these applications, interesting and useful technologies are being developed, including graphical object representations and large working volumes (CAVES). Concurrently, haptic interfaces are being perfected, which enable manual interactions with virtual environments or teleoperated remote systems. The haptic system is a unique sensory system in that it can both sense the environment and allow a user to react accordingly. As a result, haptic devices not only stimulate the user with realistic sensor input (forces, tactile sensations, heat, slip, etc.), but also sense the user's actions so that realistic sensory inputs can be generated. Haptic devices are divided into two classes, depending upon the type of sensory information being simulated. The first, tactile, refers to the sense of contact with the object. The second, kinesthetic, refers to the sense of position and motion of a user's limbs along with associated forces.
Broadly, these approaches point toward the same goal: to immerse a person in a seemingly visual reality, complete with haptic feedback. However, a major deficiency with all existing force-generating devices is the requirement that they be connected to a fixed frame, thus forcing immobility on the user. State-of-the-art force-feedback devices, for example, are table mounted, requiring the device to be mounted to an immobile object in order to generate a fixed point of leverage for forces and/or torques. Consequently, no existing force feedback device allows for easy mobility and force generation. This problem is fundamental, since many virtual reality applications require large working volumes and the ability to move freely within these volumes, to provide realistic visual and audio feedback during walk-through scenarios, for example.
In summary, large, immersive environments such as CAVES currently lack haptic feedback, primarily because the existing technology will not support unrestricted motion. This leads to one conclusion that force-feedback devices must migrate as visual technologies have, that is, from the desktop to large-volume, immersive environments. However, the design of a hand-held, spatially unrestricted force-feedback device is fundamentally different from existing devices, which typically use primarily electromechanical or pneumatic actuators operating against fixed supports to achieve active force feedback. Nor is the realization of such a device intuitively obvious. To construct an n-axis joystick, requiring 1, 2, 3 to n+3 motors, presents significant challenges, for example, since the additional motors may significantly increase the cost and/or weight of the device.
The present invention addresses the need for force feedback in large, immersive environments by providing a device that uses a gyro-stabilization to generate a fixed point of leverage for the requisite forces and/or torques. In one embodiment, one or more orthogonally oriented rotating gyroscopes are used to provide a stable body or platform to which a force-reflecting device can be mounted, thereby coupling reaction forces to the user without the need for connection to a fixed frame. In one embodiment, a user-interactable member is physically coupled to a stabilized body, with the control structure used for stabilization and that used to mitigate force-feedback being substantially independent of one another, enabling different stabilization mechanisms as described herein to be used with existing force-feedback capabilities. In alternative embodiments, inventive apparatus and methods are used which take into account both the movements associated with the gyroscopic stabilization, a user's movements, and the application of torques and forces to realize a spatially unrestricted force-feedback device requiring fewer motors and structural elements. Specifically, an inventive control scheme is used in these cases to accelerate and decelerate the motor(s) associated with providing the gyroscopic stabilization such that only the desired tactile feedback is experienced by the user. All of the various approaches are applicable to single and multiple degrees of freedom.
A three-axis implementation includes a set of three, mutually perpendicular momentum wheels which form the gyro-stabilized platform, an attitude measuring device, and a control system. The attitude measuring device is employed to detect disturbances to the gyro-stabilized platform, including reaction torques due to a user's interactions with the device. The control system varies the speed the momentum wheels in order to maintain the gyro-stabilized platform in a fixed position. In an alternative embodiment, a reaction sphere is used to produce the requisite inertial stabilization. Since the sphere is capable of providing controlled torques about three arbitrary, linearly independence axes, it can be used in place of three reaction wheels to provide three-axis stabilization for a variety of space-based and terrestrial applications.
According to the invention, programmed amounts of rotary force are used for motion compensating and/or the stabilization of free-flying platforms, or to provide force/torque outputs from platforms to attached frames. Specific embodiments are disclosed with respect to spacecraft stabilization, as well as to the application of forces and/or torques to hand-held force generating devices, including joysticks, steering wheels, and implements of arbitrary shape for specific applications, such as sports simulations.
By way of introduction, reaction wheels use the inertia of one or more (typically up to three) rotating flywheels to generate torques. These wheels are typically accelerated using electric motors which can be controlled to increase or decrease rotary speed, thus changing rotational momentum. When the wheel on a particular axis is accelerated through increased motor torque, an equal and opposite reaction torque is generated and applied to the base upon which the wheel is mounted.
Reaction wheels are the most precise type of attitude control mechanism. However, when called upon to provide non-cyclic torques, they must be periodically unloaded by other means (i.e. when the motors have accelerated to maximum RPM in any direction, no additional acceleration can be realized in that direction unless the motors are slowed, generating torques in the opposite direction). Moreover, to provide arbitrary torques, three wheel axes must be provided.
This application describes how reaction wheels as currently only applied only to spacecraft can be extended into several other related terrestrial applications, including gyro-stabilized bodies and tethered, force-generating/reflective input devices. Preliminarily, the following description will demonstrate and how three axes of reaction wheel can be reduced into a single, reaction sphere, useful either in the space-based or terrestrial applications. A description of reaction wheels and spheres will first be presented, followed by a discussion of the extensions to such technology made possible by the invention.
The singe plane model for a torque consists of a spinning wheel attached to a frame. There is a reference frame, B, embedded in the frame and a fixed reference frame, A, in the world. Reference frame B is aligned with the axis of the spinning wheel. The system is shown in FIG. 1.
The reference frame B has two degrees of freedom with respect to reference frame A. These degrees of freedom are described by generalized coordinates q_{1}, q_{4}, where q_{1 }represents the angular degree of freedom about unit vector a_{1 }and q_{4 }represents the linear degree of freedom along unit vector a_{1}. With this model, the unit vectors in frames A and B are related by
a_{1}=b_{1}. (1)
The center of mass of the frame is located at the origin of the B reference frame. The frame is assumed to be a cube with a mass of M and height of R.
The location of the center of mass of the wheel, d, is given by the vector p^{d},
^{B}p^{d}=lb_{1 } (2)
where l is the offset from the center of the frame (in meters) and is a run-time parameter. The notation ^{A}v^{B }denotes some vector v in reference frame B with respect to reference frame A. The mass of the (assuming without loss of generality a solid cylindrical disk) wheel is given by
m=πρr^{2}h (3)
where ρ is the density of the material of the wheel, r is the radius of the disk and h is the height of the wheel; and have the units of kg/m^{3}, meters and meters respectively.
The central inertia dyadic of the wheel is given by
I^{d/d+}=I_{2}b_{1}b_{1 } (4)
where
The orientation of the wheel with respect to the frame is given by the generalized coordinates q_{3}. This generalized coordinate is about the b_{1 }axis.
The central inertia dyadic of the frame is given by
I^{f/f+}=I_{6}b_{1}b_{1 } (6)
where (assuming without loss of generality that the frame is cubic)
Since there are two rigid bodies in this model (the frame and the wheel), the angular velocities and accelerations for both must be developed.
The reference frame B is said to have a simple angular velocity in the reference frame A because there exists for all time a unit vector whose orientation in both the reference frame B and reference frame A is independent of time. This allows writing the angular velocity of reference frame B as the magnitude of its angular velocity times the fixed unit vector
^{A}ω^{B}={dot over (q)}_{1}b_{1}. (8)
To make the equations of motion concise, a generalized velocity will be defined as
u_{1}={dot over (q)}_{1}. (9)
Using the definition in Equation (8), ^{A}ω^{B }can be rewritten as
^{A}ω^{B}=u_{1}b_{1 } (10)
The wheel is said to have a simple angular velocity in the reference frame B because there exists for all time a unit vector whose orientation in both the wheel reference frame and reference frame B is independent of time. This allows writing the angular velocity of the wheel as the magnitude of its angular velocity times the fixed unit vector
^{B}ω^{d}={dot over (q)}_{7}b_{1}. (11)
The angular velocity of the wheel in reference frame A is given by
^{A}ω^{d}=(u_{1}+{dot over (q)}_{7})b_{1}. (12)
Defining another generalized velocity,
u_{7}={dot over (q)}_{7 } (13)
allows expressing Equation (12) in terms of generalized velocities only as
^{A}ω^{d}=(u_{1}+u_{7})b_{1}. (14)
The angular acceleration of reference frame B is found to be
^{A}α^{B}={dot over (u)}_{1}b_{1}. (15)
and the angular acceleration of the wheel can be written as
^{A}α^{d}=({dot over (u)}_{1}+{dot over (u)}_{7})b_{1}. (16)
The location of reference frame B is given by
^{A}p^{B}=q_{4}a_{1}. (17)
The velocity ^{A}v^{B }and acceleration ^{A}a^{B }of this frame are found to be
^{A}v^{B}={dot over (q)}_{4}a_{1 } ^{A}a^{B}={dot over (q)}_{4}a_{1 } (18)
since the unit vectors a_{i }are fixed in reference frame A. Defining a generalized velocity
u_{4}={dot over (q)}_{4 } (19)
allows rewriting Equation (18) as
^{A}v^{B}=u_{4}a_{1 } ^{A}a^{B}={dot over (u)}_{4}a_{1 } (10)
By defining the disturbance forces acting at the origin of reference frame B as
F=f_{x}a_{1}, (21)
the disturbance torque acting on the frame as
T_{d}=τ_{x}a_{1}, and (22)
and the motor torque, applied to the wheel, as
T_{d}=τ_{a}b_{1 } (23)
the equations of motion are found to be:
τ_{x}−((I_{2}+I_{6}){dot over (u)}_{1}+I_{2}{dot over (u)}_{7})=0
f_{m}−(M+m){dot over (u)}_{4}=0
τ_{d}−(I_{2}{dot over (u)}_{2}+I_{2}{dot over (u)}_{7})=0 (24)
Rewriting Equation (24) in matrix form yields
To control this system, an expression for τ_{d }that allows the system to move from any value of {q_{1},u_{1}} to any other value of {q_{1},u_{1}} in the presence of disturbance torques τ_{x }must be developed (see Section below). To gain a understanding of the system, first set τ_{x}=0. Equation (25) can now be written as
τ_{d}+I_{2}{dot over (u)}_{1}=0. (26)
From control theory, it is known that this equation is not stable since the poles lie on the imaginary axis. Thus, the form of τ_{d }required to satisfy stability criteria must meet the following two criteria:
1. It must move the poles of Equation (26) into the left-half plane.
2. It should utilize values of {q_{1}, u_{1}, q_{7}, u_{7}} to control the system as these state variables can be measured.
If the disturbance torque is not set equal to zero, then Equation (26) is rewritten as
τ_{d}+I_{2}{dot over (u)}_{1}=τ_{x } (27)
and a third requirement for the control torque is added:
3. It must be robust for a specified set of disturbance torque values and functional forms.
Some simple relationships are also developed to suggest appropriate motor parameter values and sizes for the momentum wheels. For real world application, it is important to be able to specify certain aspects of the problem, such as force produced, the period of time for which it is produced and the mass of the device. Also, to stay firmly rooted in reality, it is important to specify the power output of the motor.
Equation (28) shows the basic equations
I=mr^{2 }τ=Iα
ω=αt P=τω (28)
where I is the moment of inertial of the momentum wheel (assuming that it is a thin hoop; for a solid disk, I=mr^{2}/2 and, in reality, the actual value will fall some place in between), m is the mass of the momentum wheel, r is the radius of the momentum wheel, τ is the torque applied to the operator (which is the same as the torque produce by the motor), (τ is the angular acceleration of the momentum wheel, ω is the angular velocity of the momentum wheel, t is the period of time for which the torque is felt and P is the power output of the motor.
To feel a torque produced by a motor that is not attached to some fixed structure, the motor rotor must be accelerating. The rotor will continue to accelerate until the motor reaches its maximum angular velocity, a value that is determined by motor parameters (but the calculation of which is not important for this analysis). To increase the amount of time during which the torque can be felt, it is necessary to slow down the angular acceleration of the motor by increasing the moment of inertia of the rotor.
Equation (28) has four equations and eight parameters. Of these parameters, an equation is formed that relates m,r,t,τ and P because these are the parameters that can be controlled during the design of the device. One such form of this equation is
To generate arbitrary torques, a 3D platform consisting of three spinning wheels attached to three non-coplanar axes of a frame is required. (For simplicity, and without loss of generality, this work assumes that the axes are mutually perpendicular.) There is a reference frame, B, embedded in the frame and a fixed reference frame, A, in the world. Reference frame B is aligned with the axes of the three spinning wheels, thus defining a set of mutually perpendicular unit vectors. The system is shown in FIG. 2.
The reference frame B has six degrees of freedom with respect to reference frame A. These degrees of freedom are described by generalized coordinates q_{1}, . . . , q_{6}, where q_{1}, . . . , q_{8 }represent the angular degrees of freedom about unit vectors a_{1},a_{2},a_{3 }respectively and q_{4}, . . . , q_{6 }represent the linear degrees of freedom along unit vectors a_{1},a_{2},a_{3 }respectively. The orientation of reference frame B with respect to reference frame A is described using a Body 3: 1-2-3 representation. Table 1 shows the relationship between the unit vectors a_{1}, a_{2}, a_{3 }and b_{1}, b_{2}, b_{3}.
TABLE 1 | |||||
Direction cosines | |||||
b_{1} | b_{2} | b_{3} | |||
a_{1} | c_{2}c_{3} | −c_{2}s_{3} | s_{2} | ||
a_{2} | s_{1}s_{2}c_{3 }+ s_{3}c_{1} | −s_{1}s_{2}s_{3 }+ c_{3}c_{1} | −s_{1}c_{2} | ||
a_{3} | −c_{1}s_{2}c_{3 }+ s_{3}s_{1} | c_{1}s_{2}s_{3 }+ c_{3}s_{1} | c_{1}c_{2} | ||
The terms c_{i},s_{i }are defined as cos (q_{i}) and sin (q_{i}) respectively.
To simplify some expressions, the following terms are defined:
Since the equations of motion will be developed using the unit vectors in reference from B, the unit vectors in reference frame A are explicitly presented using the terms Z_{t }defined in Equation (30).
a_{1}=Z_{1}b_{1}+Z_{4}b_{2}+Z_{7}b_{3 }
a_{2}=Z_{2}b_{1}+Z_{5}b_{2}+Z_{8}b_{3 }
a_{3}=Z_{3}b_{1}+Z_{6}b_{2}+Z_{9}b_{3 } (31)
For simplicity, and without loss of generality, the center of mass of the frame is located at the origin of the B reference frame and the frame is assumed to be cubical with a mass of M and height of R.
The central inertia dyadics of the frame is given by
I^{ƒ/ƒ+}=I_{6}b_{1}b_{1}+I_{6}b_{2}b_{2}+I_{6}b_{3}b_{3 } (32)
where
The locations of the center of mass of the wheels, d_{i}, are given by the vectors p^{d} ^{ 1 }, where
^{B}p^{d} ^{ 1 }=lb_{1 } ^{B}p^{d} ^{ 2 }=lb_{2 } ^{B}p^{d} ^{ d }=lb_{3 } (34)
where 1 is the offset from the center of the frame (in meters). The mass of each wheel (assuming without loss of generality that each wheel is a solid cylinder) is given by
m=πρr^{2}h (35)
where ρ is the density of the material of the wheel, r is the radius of the fisk and h is the height of the wheel and have the units of kg/m^{3}, meters and meters respectively.
The central inertia dyadics of the wheels are given by
I^{d} ^{ 1 } ^{/d+} ^{ 1 }=I_{2}b_{1}b_{1}+I_{1}b_{2}b_{2}+I_{1}b_{3}b_{3 }
I^{d} ^{ 2 } ^{/d+} ^{ 2 }=I_{2}b_{1}b_{1}+I_{2}b_{2}b_{2}+I_{1}b_{3}b_{3 }
I^{d} ^{ 3 } ^{/d+} ^{ 3 }=I_{1}b_{1}b_{1}+I_{1}b_{2}b_{2}+I_{2}b_{3}b_{3 } (36)
where
The orientation of the wheels with respect to the frame are given by the generalized coordinates q_{7}, . . . , q_{9}. These generalized coordinates are about the b_{1}, b_{2}, b_{3 }axes respectively.
Since there are four rigid bodies in this model (the frame and the three wheels), the angular velocities and accelerations for all four must be developed.
The angular velocity of the frame, ^{A}ω^{B }is found to be
^{A}ω^{B}=(c_{2}c_{3}{dot over (q)}_{1}+s_{3}{dot over (q)}_{2})/b_{2}+(−c_{2}s_{3}{dot over (q)}_{1}+c_{3}{dot over (q)}_{2})b_{2}+(s_{2}{dot over (q)}_{1}+{dot over (q)}_{3})b_{3}. (38)
To make the equations of motion concise, three generalized velocities will be defined as
u_{1}=c_{2}c_{3}{dot over (q)}_{1}+s_{3}{dot over (q)}_{2 }u_{2}=−c_{2}s_{3}{dot over (q)}_{1}+c_{3}{dot over (q)}_{2 }u_{3}=s_{2}{dot over (q)}_{1}+{dot over (q)}_{3}. (39)
Using the definition in Equation (39) ^{A}ω^{B }can be rewritten as
^{A}ω^{B}=u_{1}b_{1}+u_{2}b_{2}+u_{3}b_{3 } (40)
The wheels are said to have a simple angular velocity in the reference frame B because there exists for all time a unit vector whose orientation in both the wheel reference frames and reference frame B is independent of time. This allows writing the angular velocities of the wheels as the magnitude of their angular velocity times the fixed unit vector
^{B}ω^{d} ^{ 1 }={dot over (q)}_{7}b_{1 } ^{B}ω^{d} ^{ 2 }={dot over (q)}_{8}b_{2 } ^{B}ω^{d} ^{ 2 }={dot over (q)}_{9}b_{3}. (41)
The angular velocities of the wheels in reference frame A are given by
^{A}ω^{d} ^{ 1 }=(u_{1}+{dot over (q)}_{7})b_{1}+u_{2}b_{2}+u_{3}b_{3 }
^{A}ω^{d} ^{ 2 }=u_{1}b_{1}+(u_{2}+{dot over (q)}_{8})b_{2}+u_{3}b_{3}.
^{A}ω^{d} ^{ 3 }=u_{1}b_{1}+u_{2}b_{2}+(u_{3}+{dot over (q)}_{9})b_{3 } (42)
Defining three more generalized velocities,
u_{7}={dot over (q)}_{7 }u_{8}={dot over (q)}_{8 }u_{9}={dot over (q)}_{9 } (43)
allows expressing Equation (42) in terms of generalized velocities only as
^{A}ω^{d} ^{ 1 }=(u_{1}+u_{7})b_{1}+u_{2}b_{2}+u_{3}b_{3 }
^{A}ω^{d} ^{ 2 }=u_{1}b_{1}+(u_{2}+u_{8})b_{2}+u_{3}b_{3}.
^{A}ω^{d} ^{ 3 }=u_{1}b_{1}+u_{2}b_{2}+(u_{3}+u_{9})b_{3 } (44)
The angular acceleration of reference frame B is found to be
^{A}α^{B}={dot over (u)}_{1}b_{1}+{dot over (u)}_{2}b_{2}+{dot over (u)}_{3}b_{3}. (45)
The angular accelerations of the wheel can be written as
^{A}α^{d} ^{ 1 }=({dot over (u)}_{1}+{dot over (u)}_{7})b_{1}+({dot over (u)}_{2}+u_{3}u_{7})b_{2}+({dot over (u)}_{3}−u_{2}u_{7})b_{3 }
^{A}α^{d} ^{ 2 }=({dot over (u)}_{1}−u_{3}u_{8})b_{1}+({dot over (u)}_{2}+{dot over (u)}_{8})b_{2}+({dot over (u)}_{3}+u_{1}u_{8})b_{3 }
^{A}α^{d} ^{ 3 }=({dot over (u)}_{1}+u_{2}u_{9})b_{1}+({dot over (u)}_{2}−u_{1}u_{9})b_{2}+({dot over (u)}_{3}+{dot over (u)}_{9})b_{3 } (46)
The following terms are defined to simplify the equations
thus allowing Equation (46) to be rewritten as
^{A}α^{d} ^{ 1 }=Z_{10}b_{1}+Z_{13}b_{2}+Z_{16}b_{3 }
^{A}α^{d} ^{ 2 }=Z_{11}b_{1}+Z_{14}b_{2}+Z_{17}b_{3 }
^{A}α^{d} ^{ 3 }=Z_{12}b_{1}+Z_{15}b_{2}+Z_{18}b_{3 } (48)
There are four points of interest in this problem: the location of reference frame B and the locations of the centers of mass for each of the wheels. The location of reference frame B is given by
^{A}p^{B}=q_{4}a_{1}+q_{5}a_{2}+q_{6}a_{3}. (49)
and its viscosity ^{A}v^{B }and acceleration ^{A}a^{B }are found to be
^{A}v^{B}={dot over (q)}_{4}a_{1}+{dot over (q)}_{5}a_{2}+{dot over (q)}_{6}a_{3 } ^{A}α^{B}=
since the unit vectors a_{i }are fixed in reference frame A.
Defining generalized velocities
u_{4}={dot over (q)}_{4 }u_{5}={dot over (q)}_{5 }n_{6}={dot over (q)}_{6 } (51)
allows rewriting Equation (50) as
^{A}v^{B}u_{4}a_{1}+u_{5}a_{2}+u_{6}a_{3 } ^{A}a^{B}{dot over (u)}_{4}a_{1}+{dot over (u)}_{5}a_{2}+{dot over (u)}_{6}a_{3 } (52)
The velocities of the centers of the masses of the wheels are found to be
^{A}v^{d} ^{ 1 }=^{A}v^{B}+^{A}ω^{B}×lb_{1}=u_{4}a_{1}+u_{5}a_{2}+u_{6}a_{3}+lu_{3}b_{2}−lu_{2}b_{3 }
^{A}v^{d} ^{ 2 }=^{A}v^{B}+^{A}ω^{B}×lb_{2}=u_{4}a_{1}+u_{5}a_{2}+u_{6}a_{3}−lu_{3}b_{1}+lu_{1}b_{3}.
^{A}v^{d} ^{ 3 }=^{A}v^{B}+^{A}ω^{B}×lb_{3}=u_{4}a_{1}+u_{5}a_{2}+u_{6}a_{3}+lu_{2}b_{1}−lu_{1}b_{2 } (53)
and the accelerations of the centers of the masses of the wheels are found to be
^{A}a^{d} ^{ 1 }={dot over (u)}_{4}a_{1}+{dot over (u)}_{5}a_{2}+{dot over (u)}_{6}a_{3}−l(u_{2} ^{2}+u_{3} ^{2})b_{2}+l(u_{1}u_{2}+{dot over (u)}_{3})b_{2}+l(u_{1}u_{3}−{dot over (u)}_{2})b_{3 }
^{A}a^{d} ^{ 2 }={dot over (u)}_{4}a_{1}+{dot over (u)}_{5}a_{2}+{dot over (u)}_{6}a_{3}+l(u_{2}u_{2}−{dot over (u)}_{2})b_{2}+l(u_{2}u_{3}−{dot over (u)}_{1})b_{2}−l(u_{1} ^{2}+u_{2} ^{2})b_{3 } (54)
There are three sets of forces acting on this system: the applied disturbance forces and torques applied to reference frame B that represent, the gravity forces acting on the wheel and frame masses and the motor torques applied to the wheels.
The disturbance forces acting at the origin of reference frame B is defined as
F=f_{x}a_{1}+f_{y}a_{2}+f_{z}a_{3 } (55)
the gravity force on the frame is defined as
G_{ƒ}=−Mga_{3 } (56)
or equivalency as
G_{ƒ}=−Mg(Z_{3}b_{1}+Z_{6}b_{2}+Z_{9}b_{3}); (57)
and three gravity forces, which act at the center of the wheels, are defined as
G_{d} _{ 1 }=−mga_{3 }G_{d} _{ 2 }=−mga_{3 }G_{d} _{ 3 }=−mga_{3 } (58)
or equivalently as
G_{d} _{ 1 }=−mg(Z_{3}b_{1}+Z_{6}b_{2}+Z_{9}b_{3})
G_{d} _{ 2 }=−mg(Z_{3}b_{1}+Z_{6}b_{2}+Z_{9}b_{3})
G_{d} _{ 3 }=−mg(Z_{3}b_{1}+Z_{6}b_{2}+Z_{9}b_{3}) (59)
The disturbance torque acting on the frame is defined as
T_{D}=τ_{x}a_{1}+τ_{y}a_{2}+τ_{z}a_{3 } (60)
and the three motor torques are defined as
T_{d} _{ 1 }=τ_{d} _{ 1 }b_{1 }T_{d} _{ 2 }=τ_{d} _{ 2 }b_{2 }T_{d} _{ 3 }=τ_{d} _{ 3 }b_{3 } (61)
with the positive sense of the torque being applied to the wheel. However, Newton's second law demands that there be an equal and opposite torque applied to the frame, body B. Thus, the resultant acting on body B is given by
T=τ_{x}a_{1}+τ_{y}a_{2}+τ_{z}a_{3}−τ_{d} _{ 1 }b_{1}−τ_{d} _{ 2 }b_{2}−τ_{d} _{ 3 }b_{3}. (62)
The definitions of the generalized inertia forces was facilitated by defining the following terms:
The equations of motion are found to be:
f_{y}−m(3{dot over (u)}_{5}+(−Z_{19}+Z_{20}+Z_{21})Z_{2}+(Z_{22}−Z_{23}+Z_{24})Z_{5}+(Z_{25}+Z_{26}−Z_{27})Z_{8})−M{dot over (u)}_{5}=0
f_{z}−m(3g+3{dot over (u)}_{6}+(−Z_{19}+Z_{20}Z_{21})Z_{3}+(Z_{22}−Z_{23}+Z_{24})Z_{6}+(Z_{25}+Z_{26}−Z_{27})Z_{9})−M{dot over (u)}_{6}=0
τ_{d} _{ 1 }−I_{2}Z_{10}=0
τ_{d} _{ 2 }−I_{2}Z_{14}=0
τ_{d} _{ 3 }−I_{2}Z_{18}=0
Since there are nine generalized coordinates, there are 18 equations of motion, nine kinematics and nine dynamic. To solve these equations numerically, they must be written in the form
{dot over (y)}=f(y) (65)
where the state vector y has the form y={q_{1}, . . . , q_{n}, u_{1}, . . . , u_{n}}. This necessitates rewriting Equation (39) to solve for the q_{i }in terms of the u_{i}. The nine kinematic equation of motion can now be written as
To write the dynamical equations of motion in the same manner, Equation (64) must be solved for the {dot over (u)}_{i}. Since several of the Z_{i }include {dot over (u)}_{t}, these terms will need to be expanded. As a first step, the following terms are defined
Z_{28}=a_{1}u_{2 }Z_{29}=u_{1}u_{3 }Z_{30}=u_{2}u_{3 } (67)
K_{1}=I_{2}(u_{2}u_{9}−u_{3}u_{8})−mlg(Z_{6}−Z_{9})−τ_{x}Z_{1}−τ_{y}Z_{2}−τ_{z}Z_{3 }
K_{2}=I_{2}(u_{3}u_{7}−u_{1}u_{9})−mlg(Z_{9}−Z_{3})−τ_{x}Z_{4}−τ_{y}Z_{5}−τ_{z}Z_{6 }
K_{3}=I_{2}(u_{1}u_{8}−u_{2}u_{7})−mlg(Z_{3}−Z_{6})−τ_{x}Z_{7}−τ_{y}Z_{8}−τ_{z}Z_{9 } (68)
K_{4}=m((l(Z_{28}+Z_{29})−Z_{19})Z_{1}+(l(Z_{28}+Z_{30})−Z_{23})Z_{4}+(l(Z_{29}+Z_{30})−Z_{27})Z_{7})−f_{x }
K_{5}=m((l(Z_{28}+Z_{29})−Z_{19})Z_{2}+(l(Z_{28}+Z_{30})−Z_{25})Z_{5}+(l(Z_{29}+Z_{30})−Z_{27})Z_{8})−f_{y }
K_{6}=m(3g+(l(Z_{28}+Z_{29})−Z_{19})Z_{3}+(l(Z_{28}+Z_{30})−Z_{23})Z_{6}+(l(Z_{29}+Z_{30})−Z_{27})Z_{9})−f_{z }
K_{7}=−τ_{d} _{ 1 }K_{8}=−τ_{d} _{ 2 }K_{9}=−τ_{d} _{ 3 }
Next, Equation (64) is rewritten in the form A{dot over (u)}_{i}=K_{i}, thus providing a means for solving for the {dot over (u)}_{i}. K_{i }is given in Equation (68) and A is defined as
where I_{4}=2I_{1}+I_{2}+I_{6}.
A discussion of the control system is presented in below. To simplify the equations of motion to facilitate control development, those terms and equations that deal with the linear position/force are eliminated because a gyro-stabilized platform can only counteract torques, not forces. Rewriting Equation (69) as indicated yields
This equation can also be rewritten to explicitly express {dot over (u)}_{i }
where I_{5}=(I_{4}+2ml^{2})−I_{2}.
These equations of motion for the sphere, see
τ_{x}Z_{1}+τ_{y}Z_{2}+τ_{z}Z_{3}−[(I_{1}+I_{6}){dot over (u)}_{2}+I_{1}{dot over (u)}_{7}+I_{1}(u_{2}u_{9}−u_{3}u_{8})]=0
τ_{x}Z_{4}+τ_{y}Z_{5}+τ_{z}Z_{6}−[(I_{1}+I_{6}){dot over (u)}_{2}+I_{1}{dot over (u)}_{8}+I_{1}(u_{3}u_{7}−u_{1}u_{9})]=0
τ_{x}Z_{7}+τ_{y}Z_{8}+τ_{z}Z_{9}−[(I_{1}+I_{6}){dot over (u)}_{3}+I_{1}{dot over (u)}_{9}+I_{1}(u_{1}u_{8}−u_{2}u_{7})]=0
f_{x}−(m+M){dot over (u)}_{4}=0
f_{y}−(m+M){dot over (u)}_{5}=0
f_{z}−(m+M)3{dot over (u)}_{6}−mg=0
τ_{d} _{ 1 }−I_{2}Z_{11}=0
τ_{d} _{ 2 }−I_{2}Z_{14}=0
τ_{d} _{ 3 }−I_{2}Z_{18}=0
Because the crossing-coupling through the linear velocity terms does not exist for this device, controlling a system that employs this device for stabilization is easier than controlling a system that employs three reactions wheels for stabilization.
Control theory is defined as a division of engineering mathematics that attempts, through modeling, to analyze and to command a system in a desired manner. Of particular interest are closed-loop systems. In a closed-loop system, the forcing signals of the system (calling inputs) are determined (at least partially) by the responses (or output) of the system. In this manner, the inputs and outputs are interrelated. In
Control theory can be classified in two categories: classical and modern. Classical control theory is generally a trial-and-error system in which various type of analyses are used iteratively to force a electromechanical system to behave in an acceptable manner. In classical control design, the performance of a system is measured by such elements as setting time, overshoot and bandwidth. However, for highly complex, multi-input/multi-output (MIMO) systems entirely different methods of control system design should be implemented to meet the demands of modern technology. Modern control has seen wide-spread usage within the last fifteen years or so. Advancements in technology, such as faster computers, cheaper and more reliable sensors and the integration of control considerations in product design, have made it possible to extend the practical applications of automatic control to systems that were impossible to deal with in the past using classical approaches. Modern control theories are capable of dealing with issues such as performance and robustness. The spatially-unrestricted force-feedback system makes use of two modern control design methods: disturbance rejection and optimal control.
In the design of electromechanical systems, one can consider that the system will be exposed to disturbances. A disturbance may be defined as an unwanted input. In
Disturbance rejection design can be used to create a compensator which is able to ignore the disturbance and cause the desired plant output. In this section, the basic method of disturbance rejection design is presented during a MIMO model. For this model, notation must be established to designate the various elements of the control device; let:
[A, B, C, D] be a state-space representation of the plant (with state x), assuming (A,B) is completely controllable,
x(t)ε
The lumped MIMO linear, time-invariant (LTI) system, may be expressed as:
{dot over (x)}=Ax+Bu+Ew (73)
y=Cx+Du+Fw (74)
e=r−y (75)
The model for the input (Equations (76)-(77)) and the noise (Equations (78)-(79)) are:
{dot over (x)}_{r}=A_{r}x_{r } (76)
r=C_{r}x_{r } (77)
{dot over (x)}_{w}=A_{w}x_{w } (78)
w=C_{w}x_{w } (79)
The objectives in the design of the feedback system in
If this is true, then for all initial states of the system, e(t)→0ε
^{n} ^{ 0 }as t→∞.Given the system [A,B,C,D], suppose it is minimal. Let the compensator be given by
{dot over (x)}_{o}=A_{c}x+B_{c}e (80)
where
A_{c}=diag[Γ_{1},Γ_{2}, . . . , Γ_{n} _{ o }]ε
Since A_{w }and A_{r }are known, [83] can be derived from the equation
Φ(A_{w};A_{r})=s^{Q}+α_{1}s^{Q−1}+. . .+α_{q−1}+s+α_{1 } (85)
which is the least-common multiple of the characteristic equations of A_{w }and A_{r}.
Under these conditions, if
(which generates that the system is still completely controllable with the addition of the compensator) then
The discussion contained here is establishes a mathematical basis for the invention. Control of a gyro-stabilized force feedback device is based on its ability to respond robustly to a control signal and to respond correctly despite system noise. For the single-input/single-output (SISO) case, this theorem reduces to the classical control case where an integrator is required for robust performance. This result is used in the design of the 1D experiment which is similar to the classical satellite control problem.
Optimal control theory can be used to design compensators which are able to take into account the cost of performing a particular action. A classical example of optimal control is the use of fuel to maneuver a satellite in orbit above the earth. Two extreme scenarios are possible: movement taking minimum time or movement taking minimum fuel. In the following section, discussion will focus on the fundamental principles of optimal-control design.
The optimal control problem is to find a control u*(t) which causes the system {dot over (x)}(t)=a[x(t), u(t), t] to follow a desired trajectory x* that minimizes the performance measure
Other names for J include cost function, penalty function, and performance index. Assume that the admissible state and control regions are not bounded. (This removes all mechanical constraints; these can be included in later development) Let the initial states, x(t_{o})=x_{o}, of the system and initial time, t_{o}, be known. Also, let xε
^{n }and uε ^{m}. The goal now is to establish tote necessary conditions for optimality:Assuming that h is differentiable and that initial conditions are fixed and do not affect minimization, [b 87] can be expressed as
For generality, apply the chain rule and include differential equation constraints to form an augmented cost function:
using Lagrange multipliers p_{1}(t), . . . , p_{n}(t). To simplify the notation, rewrite [] as follows:
The necessary conditions for optimal control can be derived using calculus of variations. Specifically, take the variations of the functional J_{a}(u) by δx, δ{dot over (x)}, δu, δp and δt_{f}. (Increment of the functional J is defined as: ΔJ(x, δx)=δJ(x, δx)+g(x, δx)·∥δx∥; δJ is linear with respect to δx; δx is called the variation of the function x.) From this, the necessary conditions may be derived:
for all tε[t_{o}, t_{f}], and
where
N(x(t),u(t),p(t),t)≡g(x(t),u(t),p(t),t)+p^{2}[a(x(t),u(t),p(t),t] (95)
The principles of calculus of variations are applied to the design of a linear regulator. The linear regulator is used in the control of the motors used to spin the inertial masses to change the attitude of the satellite system. The regulator design is particularly useful in controlling unstable systems through optimal pole placement. First, recall the state equation of a linear, time-varying plant:
{dot over (x)}(t)=A(t)x(t)+B(t)u(t) (96)
The cost function to be used is
where t_{f }is fixed, H and Q are real, positive-semi-definite matrices, and R is a real, positive-definite matrix. The purpose of the regulator is to maintain the state of the system as close to a desired set of parameters as possible without excessive control effort. The necessary conditions for optimality to be used are:
{dot over (x)}*(t)=A(t)x*(f)+B(t)u*(t) (98)
{dot over (p)}*(t)=−Q(t)x*(t)−A^{r}(t)p*(t) (99)
O=R(t)u*(t)+B^{T}(t)p*(t) (100)
where the Hamiltonian is defined as
Equation (100) is easily solved for the optimal input for the regulation, yielding
u*(t)=−R^{−1}(t)B^{T}(t)p*(t) (102)
It is now possible to form an augmented, closed-loop state-space equation of the regulated system:
These 2n differential equations have a solution of the form:
Note: φ(t_{f},t) is called the transition matrix, define by d/dtφ(t_{f},t)=A(t)φ(t_{f},t) with the initial condition of φ(t_{o},t_{o})=1 and is solved through numerical integration. By partitioning the state transition matrix, φ(t_{f},t), the following solution for p*(t) can be reached:
p*(t)=K(t)x*(t). (105)
Therefore, the optimal control law is
u*(t)=−R^{−1}(t)B^{T}(t)K(t)x*(t)=F(t)(x*(t)); (106)
the next step is to define a method of solving for K. This is achieved using a Riccati-type differential equation:
{dot over (K)}(t)=−K(t)A(t)−A^{T}(t)K(t)−Q(t)+K(t)B(t)R^{−1}(t)B^{T}(t)K(t) (107)
which involves solving n(n+1)/2 first-order differential equations. Fortunately, the motor system involved in the hand controller control system can be considered time invariant. This simplifies the previous equations, which can be summarized as:
{dot over (x)}(t)=Ax(t)+Bu(t) (108)
and the optimal control law is
u*(t)=−R^{−1}B^{T}Kx*(t)=Fx*(t). (111)
As long as Q is positive definite, the closed-loop system is guaranteed to be stable and the controller may be used for pole placement design of the system, as shown in FIG. 6.
The design of the controller system for the 1D model is now presented. The first segment of the design is a optimal pole-placement. This is needed because a the 1D model of the spatially unrestricted force feedback device (which is a simplified version of the actual 3D version), which can be considered a second-order system, is inherently unstable. Definitions of “stable” vary; here, “stable” is considered any plant which has only poles and zeros to the left of the imagery axis in the complex plane (i.e., left-hand poles and zeros). Using previously established results, the poles of the system are placed optimally based on the inertia of a second-order linear model. Lastly, disturbance rejection is augmented to the control system for robustness.
The plant for a single DOF hand controller,
where α(t) is the angular acceleration, I_{6 }is the inertial mass, and τ(t) is the torque.
Since the stability of this system is (at best) marginal, a pole placement is performed. Further, optimal methods are employed for placing these poles at the best locations. The new plant will follow the model in FIG. 8. The optimal design will give the “best” values to use for K_{1 }and K_{2}.
The first step is to choose the cost function to minimize, set initial conditions, and select the necessary conditions and boundary conditions which apply to this problem. Let the initial states of the satellite be zero: x(0)=0; {dot over (x)}(0)=0. The cost function for minimal control effect is
such that the amount of acceleration of the system, whether it is positive or negative, for all time is minimal. This is frequently used for satellites because the amount of acceleration is the magnitude of the control input, or for satellite, the amount of fuel, which is a limited resource. For this system, the following parameters are known:
with a state defined by
and choose
and R=1.
For the LTI Ricatti equation, [110], K has four solutions, but the only positive-definite solution is
which results in a regulator, F, of the following form:
The two terms of the F vector are the position feedback and velocity feedback required for optimal tracking, as in FIG. 8.
The final step is to include an integrator which provides the SISO case with robustness. The final controller design is shown in FIG. 9.
There are some control issues that are specific to the momentum wheel concept. These issues are those that deal with determining the state of the sphere, which must be known to calculate the sphere's angular momentum vector. Since the nature of a spherical object allows it to be at any orientation relative to it's cavity, a method that can detect the sphere's exact orientation relative to the three fixed orthogonal axis of the sphere housing is used. This is illustrated in FIG. 10.
Each of the three sphere housing axis is outfitted with a band of optical infra-red emitters to detect the relative position of the sphere. Each emitter will be placed between two (or more) infra-red detectors as shown in FIG. 11. This technique will enable fine position sensing and simultaneously minimize power requirements since a single emitter will service two (or more) detectors.
The sphere is equipped with a single great circle band of reflective material as illustrated in FIG. 12. As shown in
The IR emitter/detector sensors are located directly on the cavity face to simplify construction of the sphere housing. Each emitter and detector is directly interfaced to the housing cavity by a fiber optic cable that ends at a lens mounted on the cavity face as shown in FIG. 14. Using a lens permits the use of lower power infra-red emitters.
As shown in
There are two inter-related branches of mechanics that are used to spacecraft control: celestial mechanics and attitude mechanics. The former deals with the position and velocity of the center of mass of the spacecraft as it travels through space, whereas the latter deals with the motion of the spacecraft about its center of mass, see FIG. 16.
Attitude mechanics is divided into three components: determination, prediction and control. Attitude determination is the process of computing the current orientation of the spacecraft with respect to some specified inertial frame. Attitude prediction is the process of computing the future attitude of the spacecraft based on its current state and motion. Attitude control is the process of applying torques to the spacecraft to reorient it into some desired future state. The devices mentioned in this patent deal primarily with the control aspect of attitude mechanics.
For most modern spacecraft applications, three-axis control is required. This method of control allows mission planners to specify the orientation of the spacecraft at all times during the course of a mission. Missions that employ this type of control include all communications satellites, the space shuttle and earth-orbiting scientific satellites.
To function properly, three axis stabilized spacecraft employ sensing devices that identify the spacecraft's attitude by determining two mutually perpendicular orientation vectors. Some typical examples include two-axis sun sensors and magnetic field sensors. Once the spacecraft's attitude is determined, the mission profile determines the control requirements. Certain scientific satellites require extremely precise attitude control (arc-seconds) for the purpose of data collection. Others, such as C-band television satellite, require less precise control (arc-minutes). Since all satellites are subject to disturbances, some method of maintaining proper orientation is required.
There are three primary means for controlling a satellites attitude: gas jets, electromagnets and reaction wheels. Reaction jets operate by expelling gas through an orifice to impart a moment on the spacecraft. These devices can produce large (but imprecise) torques, but since they expend fuel, there on-station operating time is limited. Electromagnets operate by creating magnet fields that interact with the magnetic field of a nearby body to produce a torque on the satellite. Although these systems do no expend fuel, they only function near bodies with large magnetic fields. Reaction wheels operate by way of Newton's third law by accelerating a wheel to absorb torque that is applied to the satellite. If the applied disturbances are cyclic, these systems can operate indefinitely since there is not net gain/loss of energy. For real-world systems, reaction wheels typically operate in conjunction with gas jets, which are used to bleed off excess momentum as the wheels approach their operating condition boundaries. Reaction wheels provide a very fine degree of attitude control.
What differentiates space-based applications from other applications is not the lack of gravity but rather the fact that gravity is the same in all directions. Similar situations can occur on the Earth: system with neutral buoyancy in a liquid and system that are fixed in the direction of gravity operate under similar principles as space-based systems, see FIG. 17.
For example, consider the case where a sensor platform is to collect data from a lake over a period of time. If this platform is required to maintain a particular attitude, a gyroscopic system can be used for stabilization. Similarly, a sensor platform mounted on a research balloon may be required to maintain two-axis attitude control for the duration of the mission. Again, a gyroscopic system can be used to stabilize the two rotational degrees of freedom of this system.
Two sets of experiments were carried out with the single degree of freedom device. The first experiment was intended to validate Equation (119). A second experiment was intended to demonstrate a control system for a three DOF system.
To carry out these tests, a test stand was developed, as shown in FIG. 18.
This test setup consists of the following components:
The position, velocity, and/or acceleration provided on a user-interactable member is sensed and transmitted as a command to a computer model or simulation which implements a virtual reality force field. In turn, the force field value for the given position, velocity, and/or acceleration is sent back to the member, which generates a force command, thereby providing the user with direct kinesthetic feedback from the virtual environment traversed. Although applicable to controlling a virtual or simulated environment, the technology is also well suited to the control of a remote or physical device. Further, the present invention is suited for application to any number of axes.
The operation of the IMC system and PC interface will be best understood by referring to commonly assigned U.S. Pat. Nos. 5,389,865 and 5,459,382, and pending applications Ser. Nos. 08/513,488 and 08/543,606, the contents of each of which are incorporated herein in their entirety by reference. These patents and co-pending applications describe systems and methods for presenting forces to an operator of a remote device or to a user interacting with a virtual environment in multiple axes simultaneously mediated through a computer controlled interface system which provides a position, velocity, and/or acceleration (to be referred to generally as “force”) to a user interface which, in turn, generates an electrical signal for each of a plurality of degrees of freedom. These electrical signals are fed to a virtual reality force field generator which calculates force field values for a selected force field. These force field values are fed to the force signal generator which generates a force signal for each of the plurality of degrees of freedom of the user as a function of the generated force field. These motion commands are fed back to actuators of the user interface which provide force to the user interface and, thus, to the user in contact with the interface device.
Before discussing these applications in further detail, a background will be provided with respect to inertial stabilization as it relates to reaction wheels and space-based applications, as certain principles of spacecraft platform stabilization have, for the first time according to this invention, been applied to spatially unrestricted terrestrial control.
Since a known momentum wheel was used, the form of Equation (29) is not quite right for this experiment. Instead, this equation is rewritten as
where the factor of two is used because a solid disk, not a thin hoop, was used. The mass of the momentum wheel is 0.277 kg (measured) and the inertia of the motor rotor is ignored.
Using the motor electrical parameters and the electrical characteristics of the IMC chassis, the maximum torque that can be applied by the motor is known to be 0.18 Nm. Inserting these values into Equation (119) yields a time of 0.09 seconds.
To measure the torque, a spring with a spring constant 110 N/m was attached to the adaptor block by way of a bolt, at a distance of 0.050 m from the center of rotation. Since F=kx and τ=Fd, these terms can be related in the following manner
Of course, the equations used are very primitive and do not account for many of the real-world affects. The affects, which are primarily frictive in nature, should tend to make the displacement less than predicted and the time greater than predicted. The results of these experiments are shown in FIG. 19 and FIG. 20. The picture on the left shows the system just before current is applied to the motor. The picture on the right shows the system at maximum spring extension, which occurred five video frames, at 30 frames per second, later. The results show a displacement of 1.2 inches (0.030 m) and a time of 0.16 seconds. Given the experimental setup, these results are well within the range of experimental error, thus giving credence to the model.
Experiments were also performed to control the position of the turntable, in the face of disturbances, by controlling the speed of the momentum wheel. The equations and methods used to develop this control scheme was discussed previously. For this experiment, the same setup was used as for the previous experiment with several small modifications:
Using MATLAB, which is a PC based mathematical tool designed to aid engineers in the development of complex mathematical systems, the controller and plant were simulated. Since the amount of control input is not a particular concern, optimal control parameters were selected to produce a system that responds quickly. In the following experiments, the values q_{1}=q_{2}=10 and R=1 were selected. To select an appropriate value for the disturbance rejection gain, a root locus plot of the system,
The control parameters determined using the optimal control techniques and the root-locus method were applied to the system shown in
In this particular device, there was a great deal of Coulomb friction in the base bearing. The components mounted to the base would not complete a single rotation before coming to a halt after an initial spin. This has the effect of adding instability to the system. In particular, what tends to happen is that the system will stay at some point for some period of time while the integrator error (the disturbance rejection) adds up. At some point, there is sufficient energy to overcome the static friction, which is less than the dynamic friction. Once moving, the system will tend to overshoot the desired point and try to compensate, but the same sequence of events occurs.
An experiment was also performed to determine if the forces generated were noticeable by a human. To perform this experiment, three motors with momentum wheels were mounted onto the adapter block used in the previous experiments.
The motors were spun up to a speed of 5000 RPM. Individuals were asked to handle the device and to make subjective evaluations of the torques felt as the device was moved about. In all cases, the subjects reported feeling appreciable forces that were deemed to be sufficient for carrying out meaningful tasks. A picture of the device is shown in FIG. 24.
The torques felt were generated because the control system had been commanded to maintain the momentum wheels at a constant angular velocity. By moving the device about, the angular momentum vectors were changed, thus causing a torque. The control system compensated for these motions by adjusting the output to the motors. Since the motors were already spinning at high speed, the period of time for which a torque could be applied was far more limited than for the case where the motor is initially at rest.
Having demonstrated that forces can be generated in any direction, the final task is to control the motors in an appropriate manner so as to provide haptic feedback to the user. This task requires a sophisticated control algorithm for two reasons: first, the platform will be grossly displaced from its nominal operating orientation, and second, for any motion of the platform (for simplicity consider just rotations about the world coordinate axes with which the device is initially oriented), some subset of the motors will produce torques (due to changes in the orientation of the angular momentum vectors) that are undesired. To counteract these undesired torques, some subset of the motors will need to be accelerated to produce counter torques. The control system must model the full, non-linear dynamics of the system, have a high speed attitude sensor and possibly a control to smoothly generate the prescribed forces. A block diagram of the system is shown in FIG. 25.
As discussed above, one family of applications for the devices described above utilizes inputs received from a virtual environment. For this type of application, the virtual environment models some set of objects, and hand controller or other force-reflection device produces forces that are representative of some activity within the virtual environment. Since it is not required that the forces produced correlate to any specific activity, the only restriction placed on the commands sent to the gyro-stabilized device is that the output forces be within the range of forces that the device can produce. An alternative family of applications for these devices produces forces in accordance with inputs received from a (possibly remote) physical device. For this type of application, the forces produced are typically a scaled representation of the actual forces produced at some point on the actual physical device. To provide the widest range of haptic input, the scaling is typically designed such that the maximum force that can be applied to the physical device is mapped into the maximum force that the haptic device can produce.
To the first order, the devices described are marginally stable at best. To control these devices to produce desire torque outputs in the face of input disturbances, a two step controller is preferably utilized. The first step stabilizes the controller by doing a pole placement. The location of the poles can be determined using any applicable method although optimal control is preferred. The second step creates a robust controller by canceling out disturbance inputs. Robust control theory is applied for this task.
With specific regard to platform stabilization, the desired input is typically a zero input, i.e., that the system should not change state. For these applications, sensor are employed to determine when the system changes state due to disturbances and the controller acts to return the system to the zero state.
The human operator who controls the haptic device is, from the perspective of the momentum device, equivalent to group. Although any amount of angular momentum can be removed from the device when it is coupled to ground, since this is a haptic device, the strategy is to slowly and continually remove angular momentum so as to have as minimal affect on the user as possible. In particular, the momentum sphere has a maximum speed at which it can operate due to the materials and construction techniques employed. When the sphere approaches this maximum velocity, momentum must be unloaded from the sphere for it to continue to function. To do this requires the application of an external torque that will cause the angular momentum vector to be diminished. This can be accomplished in three ways: reaction jets, magnetic field torques and/or spacecraft reorientation. The first two methods work by applying a torque to the spacecraft that diminished the angular momentum of the sphere. The third method works if the following two conditions are met: the disturbances to the spacecraft are primarily applied in the same direction and the spacecraft can continue to operate at different attitudes. If these conditions are met, the spacecraft can be reoriented such that the disturbance torque act to cancel the sphere's angular momentum. It may also be feasible to rigidly couple the platform to ground for a brief period of time. While coupled to ground, any amount of angular momentum can be removed from the stabilized platform.
Cited Patent | Filing date | Publication date | Applicant | Title |
---|---|---|---|---|
US3157853 | Dec 6, 1957 | Nov 17, 1964 | Joseph Hirsch | Tactile communication system |
US3220121 | May 20, 1963 | Nov 30, 1965 | Communications Patents Ltd | Ground-based flight training or simulating apparatus |
US3497668 | Aug 25, 1966 | Feb 24, 1970 | Joseph Hirsch | Tactile control system |
US3517446 | Apr 19, 1967 | Jun 30, 1970 | Singer General Precision | Vehicle trainer controls and control loading |
US3623064 | Oct 11, 1968 | Nov 23, 1971 | Bell & Howell Co | Paging receiver having cycling eccentric mass |
US3902687 | Jun 25, 1973 | Sep 2, 1975 | Robert E Hightower | Aircraft indicator system |
US3903614 | Mar 27, 1970 | Sep 9, 1975 | Singer Co | Apparatus for simulating aircraft control loading |
US3911416 | Aug 5, 1974 | Oct 7, 1975 | Motorola Inc | Silent call pager |
US4160508 | Aug 19, 1977 | Jul 10, 1979 | Nasa | Controller arm for a remotely related slave arm |
US4236325 | Dec 26, 1978 | Dec 2, 1980 | The Singer Company | Simulator control loading inertia compensator |
US4262549 | May 10, 1978 | Apr 21, 1981 | Schwellenbach Donald D | Variable mechanical vibrator |
US4333070 | Feb 6, 1981 | Jun 1, 1982 | Barnes Robert W | Motor vehicle fuel-waste indicator |
US4443952 * | Nov 9, 1981 | Apr 24, 1984 | The Bendix Corporation | Gyroscopic apparatus |
US4464117 | Aug 26, 1981 | Aug 7, 1984 | Dr. Ing. Reiner Foerst Gmbh | Driving simulator apparatus |
US4484191 | Jun 14, 1982 | Nov 20, 1984 | Vavra George S | Tactile signaling systems for aircraft |
US4513235 | Jan 24, 1983 | Apr 23, 1985 | British Aerospace Public Limited Company | Control apparatus |
US4581491 | May 4, 1984 | Apr 8, 1986 | Research Corporation | Wearable tactile sensory aid providing information on voice pitch and intonation patterns |
US4599070 | Jul 29, 1981 | Jul 8, 1986 | Control Interface Company Limited | Aircraft simulator and simulated control system therefor |
US4708656 | Feb 4, 1986 | Nov 24, 1987 | Fokker B.V. | Simulator of mechanical properties of a steering system |
US4713007 | Oct 11, 1985 | Dec 15, 1987 | Alban Eugene P | Aircraft controls simulator |
US4794392 | Feb 20, 1987 | Dec 27, 1988 | Motorola, Inc. | Vibrator alert device for a communication receiver |
US4891764 | Dec 11, 1987 | Jan 2, 1990 | Tensor Development Inc. | Program controlled force measurement and control system |
US4930770 | Dec 1, 1988 | Jun 5, 1990 | Baker Norman A | Eccentrically loaded computerized positive/negative exercise machine |
US4934694 | Mar 9, 1988 | Jun 19, 1990 | Mcintosh James L | Computer controlled exercise system |
US5019761 | Feb 21, 1989 | May 28, 1991 | Kraft Brett W | Force feedback control for backhoe |
US5022384 | May 14, 1990 | Jun 11, 1991 | Capitol Systems | Vibrating/massage chair |
US5022407 | Jan 24, 1990 | Jun 11, 1991 | Topical Testing, Inc. | Apparatus for automated tactile testing |
US5035242 | Apr 16, 1990 | Jul 30, 1991 | David Franklin | Method and apparatus for sound responsive tactile stimulation of deaf individuals |
US5038089 | Oct 28, 1988 | Aug 6, 1991 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Synchronized computational architecture for generalized bilateral control of robot arms |
US5044956 * | Jan 12, 1989 | Sep 3, 1991 | Atari Games Corporation | Control device such as a steering wheel for video vehicle simulator with realistic feedback forces |
US5078152 | Dec 25, 1988 | Jan 7, 1992 | Loredan Biomedical, Inc. | Method for diagnosis and/or training of proprioceptor feedback capabilities in a muscle and joint system of a human patient |
US5165897 | Aug 10, 1990 | Nov 24, 1992 | Tini Alloy Company | Programmable tactile stimulator array system and method of operation |
US5175459 | Aug 19, 1991 | Dec 29, 1992 | Motorola, Inc. | Low profile vibratory alerting device |
US5186695 | Oct 26, 1990 | Feb 16, 1993 | Loredan Biomedical, Inc. | Apparatus for controlled exercise and diagnosis of human performance |
US5212473 | Feb 21, 1991 | May 18, 1993 | Typeright Keyboard Corp. | Membrane keyboard and method of using same |
US5240417 | Mar 14, 1991 | Aug 31, 1993 | Atari Games Corporation | System and method for bicycle riding simulation |
US5271290 | Apr 14, 1993 | Dec 21, 1993 | United Kingdom Atomic Energy Authority | Actuator assembly |
US5275174 | Jul 16, 1992 | Jan 4, 1994 | Cook Jonathan A | Repetitive strain injury assessment |
US5283970 | Sep 25, 1992 | Feb 8, 1994 | Strombecker Corporation | Toy guns |
US5299810 | Jun 23, 1992 | Apr 5, 1994 | Atari Games Corporation | Vehicle simulator including cross-network feedback |
US5309140 | Nov 26, 1991 | May 3, 1994 | The United States Of America As Represented By The Secretary Of The Navy | Feedback system for remotely operated vehicles |
US5334027 | Feb 25, 1991 | Aug 2, 1994 | Terry Wherlock | Big game fish training and exercise device and method |
US5389865 * | Dec 2, 1992 | Feb 14, 1995 | Cybernet Systems Corporation | Method and system for providing a tactile virtual reality and manipulator defining an interface device therefor |
US5436622 | Jul 6, 1993 | Jul 25, 1995 | Motorola, Inc. | Variable frequency vibratory alert method and structure |
US5466213 | Jan 6, 1994 | Nov 14, 1995 | Massachusetts Institute Of Technology | Interactive robotic therapist |
US5481914 * | Mar 28, 1994 | Jan 9, 1996 | The Charles Stark Draper Laboratory, Inc. | Electronics for coriolis force and other sensors |
US5547382 | Apr 10, 1995 | Aug 20, 1996 | Honda Giken Kogyo Kabushiki Kaisha | Riding simulation system for motorcycles |
US5575761 | Jul 27, 1994 | Nov 19, 1996 | Hajianpour; Mohammed-Ali | Massage device applying variable-frequency vibration in a variable pulse sequence |
US5577981 * | Aug 4, 1995 | Nov 26, 1996 | Jarvik; Robert | Virtual reality exercise machine and computer controlled video system |
US5709219 * | Aug 1, 1996 | Jan 20, 1998 | Microsoft Corporation | Method and apparatus to create a complex tactile sensation |
US5766016 | Nov 14, 1994 | Jun 16, 1998 | Georgia Tech Research Corporation | Surgical simulator and method for simulating surgical procedure |
US5785630 | Nov 6, 1996 | Jul 28, 1998 | Tectrix Fitness Equipment, Inc. | Interactive exercise apparatus |
US5872438 * | Jul 1, 1997 | Feb 16, 1999 | Cybernet Systems Corporation | Whole-body kinesthetic display |
US6004134 * | May 19, 1994 | Dec 21, 1999 | Exos, Inc. | Interactive simulation including force feedback |
US6111577 | Apr 4, 1996 | Aug 29, 2000 | Massachusetts Institute Of Technology | Method and apparatus for determining forces to be applied to a user through a haptic interface |
US6160489 | Jun 23, 1994 | Dec 12, 2000 | Motorola, Inc. | Wireless communication device adapted to generate a plurality of distinctive tactile alert patterns |
US6219034 | Feb 23, 1998 | Apr 17, 2001 | Kristofer E. Elbing | Tactile computer interface |
US6422941 | Sep 23, 1997 | Jul 23, 2002 | Craig Thorner | Universal tactile feedback system for computer video games and simulations |
JPH047371A | Title not available | |||
JPH02109714A | Title not available | |||
JPH05193862A | Title not available | |||
JPS643664A | Title not available |
Reference | ||
---|---|---|
1 | "Coaxial Control Shaker Part No. C-25502," Safe Flight Instrument Corporation, 26 pages, Jul. 1, 1967; Revised Jan. 28, 2002. | |
2 | "Taking a Joystick Ride", Computer Currents, Tim Scannell, Nov. 1994, Boston Edition, vol. 9 No. 11. | |
3 | Adelstein, "A Virtual Environment System For The Study of Human Arm Tremor," Ph.D. Dissertaton, Dept. of Mechanical Engineering, MIT, Jun. 1989. | |
4 | Adelstein, "Design and Implementation of a Force Reflecting Manipulandum for Manual Control research," DSC-vol. 42, Advances in Robotics, Edited by H. Kazerooni, pp. 1-12, 1992. | |
5 | Baigrie, "Electric Control Loading-A Low Cost, High Performance Alternative," Proceedings, pp. 247-254, Nov. 6-8, 1990. | |
6 | Bejczy et al., "A Laboratory Breadboard System For Dual-Arm Teleoperation," SOAR '89 Workshop, JSC, Houston, TX, Jul. 25-27, 1989. | |
7 | Bejczy et al., "Kinesthetic Coupling Between Operator and Remote Manipulator," International Computer Technology Conference, The American Society of Mechanical Engineers, San Francisco, CA, Aug. 12-15, 1980. | |
8 | Bejczy, "Generalization of Bilateral Force-Reflecting Control of Manipulators," Proceedings Of Fourth CISM-IFToMM, Sep. 8-12, 1981. | |
9 | Bejczy, "Sensors, Controls, and Man-Machine Interface for Advanced Teleoperation," Science, vol. 208, No. 4450, pp. 1327-1335, 1980. | |
10 | Bejczy, et al., "Universal Computer Control System (UCCS) For Space Telerobots," CH2413-3/87/0000/0318501.00 1987 IEEE, 1987. | |
11 | Bliss, "Optical-to-Tactile Image Conversion for the Blind," IEEE Transactions on Man-Machine Systems, vol. MMS-11, No. 1, Mar 1970. | |
12 | Brooks et al., "Hand Controllers for Teleoperation-A State-of-the Art Technology Survey and Evaluation," JPL Publication 85-11; NASA-CR-175890; N85-28559, pp. 1-84, Mar. 1, 1985. | |
13 | Burdea et al., "Distributed Virtual Force Feedback, Lecture Notes for Workshop on Force Display in Virtual Environments and its Application to Robotic Teleoperation," 1993 IEEE International Conference on Robotics and Automation, pp. 25-44, May 2, 1993. | |
14 | Cadler, "Design of A Force-Feedback Touch-Introducing Actuator For Teleoperator Robot Control," Bachelor of Science Thesis, MIT, Jun. 23, 1983. | |
15 | Caldwell et al., "Enhanced Tactile Feedback (Tele-Taction) Using a Multi-Functional Sensory System," 1050-4729/93, pp. 955-960, 1993. | |
16 | Eberhardt et al., "Omar-A Haptic display for speech perception by deaf and deaf-bind individuals," IEEE Virtual Reality Annual International Symposium, Seattle, WA, Sep. 18-22, 1993. | |
17 | Gotow et al., "Controlled Impedance Test Apparatus for Studying Human Interpretation of Kinesthetic Feedback," WA11-11:00, pp. 332-337. | |
18 | Howe, "A Force-Reflecting Teleoperated Hand System for the Study of Tactile Sensing in Precision Manipulation," Proceedings of the 1992 IEEE International Conference on Robotics and Automation, Nice, France, May 1992. | |
19 | IBM Technical Disclosre Bullein, "Mouse Ball-Actuating Device With Force and Tactile Feedback," vol. 32, No. 9B, Feb. 1990. | |
20 | Iwata, "Pen-based Haptic Virtual Environment," 0-7803-1363-1/93 IEEE, pp. 287-292, 1993. | |
21 | Jacobsen et al., "High Performance, Dextrous Telerobotic Manipulator With Force Reflection," Intervention/ROV '91 Conference & Exposition, Hollywood, Florida, May 21-23, 1991. | |
22 | Johnson, "Shape-Memory Alloy Tactile Feedback Actuator," Armstrong Aeorspace Medical Research Laboratory, AAMRL-TR-90-039, Aug. 1990. | |
23 | Jones et al., "A perceptual analysis of stiffness," ISSN 0014-4819 Springer International (Springer-Verlag); Experimental Brain Research, vol. 79, No. 1, pp. 150-156, 1990. | |
24 | Kontarinis et al., "Display of High-Frequency Tactile Information to Teleoperators," Telemanipulator Technology and Space Telerobotics, Won S. Kim, Editor, Proc. SPIE vol. 2057, pp. 40-50, Sep. 7-9, 1993. | |
25 | Kontarinis et al., "Tactile Display of Vibratory Information in Teleoperation and Virtual Environments," Presence, 4(4):387-402, 1995. | |
26 | Marcus, "Touch Feedback in Surgery," Proceedings of Virtual Reality and Medicine The Cutting Edge, Sep. 8-11, 1994. | |
27 | McAffee, "Teleoperator Subsystem/Telerobot Demonsdtrator: Force Reflecting Hand Controller Equipment Manual," JPL D-5172, pp. 1-50, A1-A36, B1-B5, C1-C36, Jan. 1988. | |
28 | Minsky, "Computational Haptics: The Sandpaper System for Synthesizing Texture for a Force-Feedback Display," Ph.D. Dissertation, Mit, Jun. 1995. | |
29 | Ouh-Young, "A Low-Cost Force Feedback Joystick and Its Use in PC Video Games," IEEE Transactions on Consumer Electronics, vol. 41, No. 3, Aug. 1995. | |
30 | Ouh-Young, "Force Display in Molecular Docking," Order No. 9034744, p. 1-369, 1990. | |
31 | Patrick et al., "Design and Testing of A Non-reactive, Fingertip, Tactile Display for Interaction with Remote Environments," Cooperative Intelligent Robotics in Space, Rui J. deFigueiredo et al., Editor, Proc. SPIE vol. 1387, pp. 215-222, 1990. | |
32 | Patrick, "Design, Construction, and Testing of a Fingertip Tactile Display for Interaction with Virtual and Remote Environments," Master of Science Thesis, MIT, Nov. 8, 1990. | |
33 | Rabinowitz et al., "Multidimensional tactile displays: Identification of vibratory intensity, frequency, and contactor area," Journal of The Acoustical Society of America, vol. 82, No. 4, Oct. 1987. | |
34 | Russo, "Controlling Dissipative Magnetic Particle Brakes in Force Reflective Devices," DSC-vol. 42, Advances in Robotics, pp. 63-70, ASME 1992. | |
35 | Russo, "The Design and Implementation of a Three Degree of Freedom Force Output Joystick," MIT Libraries Archives Aug. 14, 1990, pp. 1-131, May 1990. | |
36 | Shimoga "Finger Force and Touch Feedback Issues in Dexterous Telemanipulation," Proceedings of Fourth Annual Conference on Intelligent Robotic Systems for Space Expploration, Rensselaer Polytechnic Institute, Sep. 30-Oct. 1, 1992. | |
37 | Snow et al., "Model-X Force-Reflecting-Hand-Controller," NT Control No. MPO-17851; JPL Case No. 5348, pp. 1-4, Jun. 15, 1989. | |
38 | Stanley et al., "Computer Simulation of Interacting Dynamic Mechanical Systems Using Distributed Memory Parallel Processors," DSC-vol. 42, Advances in Robotics, pp. 55-61, ASME 1992. | |
39 | Tadros, Circuit System Design for a Three Degree of Freedom Virtual Environment Simulator Using Motor/Brake Pair Actuators, MIT Archive (C) Massachusetts Institute of Technology, pp. 1-88, Feb. 1990. | |
40 | Terry et al., "Tactile Feedback In A Computer Mouse," Proceedings of Fouteenth Annual Northeast Bioengineering Conference, University of New Hampshire, Mar. 10-11, 1988. | |
41 | Wiker, "Teletouch Display Development Phase 1 Report," Technical Report 1230, Naval Ocean Systems Center, San Diego, Apr. 17, 1989. |
Citing Patent | Filing date | Publication date | Applicant | Title |
---|---|---|---|---|
US7867141 * | Mar 7, 2005 | Jan 11, 2011 | Panasonic Electric Works Co., Ltd. | Physical activity measuring system |
US8088042 * | Dec 7, 2004 | Jan 3, 2012 | Elisa Oyj | Method, system, measurement device and receiving device for providing feedback |
US8160731 * | Aug 18, 2009 | Apr 17, 2012 | Empire Technology Development Llc | Stabilized platform system |
US8210942 | Mar 30, 2007 | Jul 3, 2012 | Wms Gaming Inc. | Portable wagering game with vibrational cues and feedback mechanism |
US8500534 | Nov 7, 2007 | Aug 6, 2013 | Wms Gaming Inc. | Gaming machine having display with sensory feedback |
US8882575 | Jul 3, 2013 | Nov 11, 2014 | Wms Gaming Inc. | Gaming machine having display with sensory feedback |
US9058714 | May 22, 2012 | Jun 16, 2015 | Wms Gaming Inc. | Wagering game systems, wagering gaming machines, and wagering gaming chairs having haptic and thermal feedback |
US9142083 | Jun 4, 2012 | Sep 22, 2015 | Bally Gaming, Inc. | Convertible gaming chairs and wagering game systems and machines with a convertible gaming chair |
US20060020174 * | Mar 7, 2005 | Jan 26, 2006 | Yoshihiro Matsumura | Physical activity measuring system |
U.S. Classification | 318/561, 414/5, 318/568.11, 318/649 |
International Classification | G05B5/01, G06F3/01, B25J9/00, G05B13/02, F16F15/00, G05B13/04, G05G5/03, G06F3/00 |
Cooperative Classification | G05B5/01, G05B2219/37164, G05B2219/40122, G06F2203/013, G05B2219/37174, G05G5/03, F16F15/00, G06F3/016, G05B13/042, G06F2203/015, G05B2219/41274 |
European Classification | F16F15/00, G06F3/01F, G05B13/04B, G05G5/03, G05B5/01 |
Date | Code | Event | Description |
---|---|---|---|
Aug 19, 2008 | CC | Certificate of correction | |
Nov 19, 2009 | FPAY | Fee payment | Year of fee payment: 12 |
May 18, 2010 | AS | Assignment | Owner name: IMMERSION CORPORATION,CALIFORNIA Free format text: MERGER;ASSIGNOR:CYBERNET HAPTIC SYSTEMS CORPORATION;REEL/FRAME:024402/0612 Effective date: 20021002 Owner name: CYBERNET HAPTIC SYSTEMS CORPORATION,CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:ROSTON, GERALD P.;REEL/FRAME:024402/0537 Effective date: 19991104 Owner name: IMMERSION CORPORATION DELAWARE D/B/A IMMERSION COR Free format text: MERGER;ASSIGNOR:IMMERSION CORPORATION;REEL/FRAME:024402/0497 Effective date: 19991102 Owner name: IMMERSION CORPORATION, CALIFORNIA Free format text: CHANGE OF NAME;ASSIGNOR:IMMERSION HUMAN INTERFACE CORPORATION;REEL/FRAME:024402/0432 Effective date: 19980406 Owner name: IMMERSION CORPORATION,CALIFORNIA Free format text: CHANGE OF NAME;ASSIGNOR:IMMERSION HUMAN INTERFACE CORPORATION;REEL/FRAME:024402/0432 Effective date: 19980406 |
Apr 23, 2013 | CC | Certificate of correction |