The subject matter disclosed herein relates generally to active vibration control systems and methods. More particularly, the subject matter disclosed herein relates to implementation and control schemes for an active vibration control system, such as is used to control vibration in a helicopter.
It is sometimes desired to place multiple actuators, such as linear actuators or circular force generators (CFGs), close together at particular locations to increase controllability of certain modes of vibration. When this is done, however, adaptive algorithms that are commonly used to control such modes of vibration (e.g., filtered least mean squares) can have difficulty finding the optimal solution. These difficulties can generally arise either because the algorithm takes a significantly longer path to find the minimal solution (i.e., slow convergence) or because it can have a difficult time finding a unique solution, and it will thus oscillate back and forth looking for the minimum (i.e., poor performance).
As a result, it would be advantageous for systems and methods for controlling multiple actuators to quickly and accurately identify an optimal solution to generate the desired force output from the combined operation of the multiple actuators.
In accordance with this disclosure, systems, methods, and computer program products for directional force weighting of an active vibration control system are provided. In one aspect, an active vibration control system includes a plurality of force generators arranged in an array, with each of the plurality of force generators being configured to generate individual component force outputs. An even number of the plurality of force generators are arranged in pairs that are placed in close proximity to one another for multi-directional force generation. A controller is configured to individually control each of the plurality of force generators to achieve a combination of the individual component force outputs that produces a desired total force vector.
In another aspect, a method for directional force weighting of an active vibration control system is provided. The method involves arranging a plurality of force generators in an array, identifying individual component forces corresponding to force outputs of each of the plurality of force generators, determining a combination of the individual component forces that will produce a desired total force vector, and adjusting the outputs of each of the plurality of force generators such that the combination of the individual component forces are at least substantially similar to the desired force vector.
Although some of the aspects of the subject matter disclosed herein have been stated hereinabove, and which are achieved in whole or in part by the presently disclosed subject matter, other aspects will become evident as the description proceeds when taken in connection with the accompanying drawings as best described hereinbelow.
Numerous objects and advantages of the subject matter will become apparent as the following detailed description of the preferred embodiments is read in conjunction with the drawings, which illustrate such embodiments.
Actuator mapping is able to transform redundant and/or poorly conditioned degrees of freedom into simpler primary Degrees of Freedom (DOFs) in a very simple way. As shown in
In one non-limiting example, a system of linear actuators is controllable to achieve an aggregate force output that includes both linear and rotational modes of vibration. As shown in
To achieve these complex modes of vibration, a transformation matrix is applied to the inputs from each of the linear actuators to achieve a desired output. In general, an active vibration control system operating at a single frequency is described as
e=Cf+d
where e is a [n×1] complex vector of vibration signals at the frequency of interest and measured by the vibration sensors, f is a [m×1] complex vector of input force commands at the frequency of interest, C is the [n×m] complex transfer function matrix between f and e, and d is the n×1 complex vector of vibration signals measured when there is no control. The control system is functional to adapt f such that the product of C and f looks as close to −d as possible such that e is minimal (in a least squares sense).
In this regard, a force transform vector is produced:
M1=M2M3
where an output force vector M1 represents a complex vector with elements for which control weighting is desired, and a transformation matrix M2 maps the natural modes of vibration generated by an input force vector M3 to achieve the desired control weighting.
Referring again to the actuator configuration shown in
Similarly, the pairing of proximal circular force generators (CFGs) enables bidirectional force generation. There may be situations where a systems engineer will want to create a single direction force using two CFGs. For example, multiple circular forces can be mapped to independent linear forces (and vice versa). To this end, the vibration control algorithm implicitly will converge to an elliptical resultant force profile for pairs of CFGs such that a weighted sensor set is minimized. The following provides a manner for doing so by penalizing or applying control weighting to various rectilinear directions while maintaining independent CFG control.
For example, as shown in
Again, identifying the proper control weighting for each of the five CFGs in this exemplary configuration is achieved by transforming the input forces generated individually into an aggregate output force vector having the desired modes of vibration. In particular, for example,
Using this form of force transform, a cost function is defined as follows:
where Q is a sensor weighting matrix and R is a control weighting matrix. With respect to the configuration discussed above with respect to the arrangement shown in
R=diag{0,ray,rby,0,r5}
where ray, rby, are adjustable to ensure unidirectionality of first CFG pair 21a and second CFG pair 21b, respectively, and r5 provides control weighting on fifth CFG 11e.
The adaptation algorithm has the following form:
fk+1=(I−{circumflex over (R)})fk−μC*Qek
In addition, in the configuration shown in
In yet a further configuration, actuators are mounted near the transmission of a helicopter to suppress the primary DOFs: X, Y, Z, pitch, and roll. Specifically, as shown in
An exemplary mapping matrix for such a configuration is designed as shown in
In any configuration, if the control authority of a particular DOF is significantly larger or smaller than the others, it can also cause poor transient performance. A simple way to improve this is to normalize the actuator response in the plant model (C-model):
for n=1: nact
end
The present subject matter can be embodied in other forms without departure from the spirit and essential characteristics thereof. The embodiments described therefore are to be considered in all respects as illustrative and not restrictive. Although the present subject matter has been described in terms of certain preferred embodiments, other embodiments that are apparent to those of ordinary skill in the art are also within the scope of the present subject matter.
The present application claims the benefit of U.S. Provisional Patent Application Ser. No. 61/754,234, filed Jan. 18, 2013, the disclosure of which is incorporated herein by reference in its entirety.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US2014/011956 | 1/17/2014 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2014/113619 | 7/24/2014 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
5310137 | Yoerkie, Jr. | May 1994 | A |
5750897 | Kato | May 1998 | A |
6467723 | Rossetti | Oct 2002 | B1 |
6536735 | Haga | Mar 2003 | B2 |
6644590 | Terpay | Nov 2003 | B2 |
7017857 | Hill | Mar 2006 | B2 |
7822509 | Heiland | Oct 2010 | B2 |
8267652 | Jolly | Sep 2012 | B2 |
9174739 | Jamshidiat | Nov 2015 | B2 |
9284048 | Harrigan | Mar 2016 | B2 |
9488247 | Rijnveld | Nov 2016 | B2 |
Number | Date | Country |
---|---|---|
1 400 885 | Mar 2004 | EP |
0039480 | Jul 2000 | WO |
2008079431 | Jul 2008 | WO |
2012021202 | Feb 2012 | WO |
2012125924 | Sep 2012 | WO |
Number | Date | Country | |
---|---|---|---|
20150353191 A1 | Dec 2015 | US |
Number | Date | Country | |
---|---|---|---|
61754234 | Jan 2013 | US |