Journal of Applied Mathematics

Volume 2014 (2014), Article ID 430731, 14 pages

http://dx.doi.org/10.1155/2014/430731

Research Article

Innovation in Active Vibration Control Strategy of Intelligent Structures

A. Moutsopoulou,¹ G. E. Stavroulakis,² and A. Pouliezos²

¹Department of Engineering (Civil), Technological Educational Institute of Crete, 71004 Heraklion, Crete, Greece

²Department of Production Engineering and Management, Technical University of Crete, 73100 Chania, Crete, Greece

Received 19 July 2013; Revised 13 November 2013; Accepted 14 November 2013; Published 6 March 2014

Academic Editor: Weihai Zhang

Copyright © 2014 A. Moutsopoulou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Large amplitudes and attenuating vibration periods result in fatigue, instability, and poor structural performance. In light of past approaches in this field, this paper intends to discuss some innovative approaches in vibration control of intelligent structures, particularly in the case of structures with embedded piezoelectric materials. Control strategies are presented, such as the linear quadratic control theory, as well as more advanced theories, such as robust control theory. The paper presents sufficiently a recognizable advance in knowledge of active vibration control in intelligent structures.

1. Introduction

Modeling of an ideal system consisting of a beam with piezoelectric layers was done using one-dimensional finite elements with two degrees of freedom per node. Cubic and quadratic Hermit polynomials were employed for the representation of nodal rotations and vertical displacements. System differential equations are derived from Euler Bernoulli theory [1, 2].

Setting up the problem in a two-port diagram (input-output) was not a trivial task. The classic control problem was transformed into a two-port problem. The goal of nominal design was to keep error magnitude small, despite perturbations and noise in measurements. Moreover, controller size had to be contained so as to lower energy consumption and maintain piezoelectric materials within operating limits (±500 V). By transforming transfer functions to state space equations and by using input and output equations, state space matrices have been derived; these matrices are used for finding the optimal controller according to the LQR and control criterion.

Selection of the weights involved in the controller we studied was done through optimization, while wind loading and noise in measurements were appropriately modelled for this particular problem. The obtained results were very satisfactory; beam vibration is reduced even for realistic wind measurements. Beam response results, with as well as without control, were compared for all presented control strategies.

In this paper, we address the problem of vibrations of intelligent structures. Stimuli may come from external perturbations of the system, disturbances, or excitation that may cause structural vibrations, such as wind loading or earthquakes. An intelligent structure is expected to be able to sense the vibration and counteract it in a controlled fashion, so that vibration of the system can be reduced and contained. To that end, a number of intelligent materials may be used as actuators and sensors. Piezoelectric materials, memory materials, and electrostrictive and magnetostrictive materials are such materials. In this work, we focus on the use of piezoelectric materials, given that they exhibit good sensing and actuation properties.

2. Research on Intelligent Structures

The following paragraphs give a deep insight into the research work done on the intelligent structures so far. Culshaw discussed the concept of smart structure, its benefits, and applications [3]. Rao and Sunar explained the use of piezo materials as sensors and actuators in sensing vibrations in their survey paper [4]. Hubbard and Baily have studied the application of piezoelectric materials as sensor/actuator for flexible structures [5]. Hanagud developed a Finite Element Model (FEM) for a beam with many distributed piezoceramic sensors/actuators [6].

Hwang and Park presented a new finite element (FE) modeling technique for flexible beams [7]. Continuous time and discrete time algorithms were proposed to control a thin piezoelectric structure by Bona et al. [8]. Schiehlen and Schonerstedt reported the optimal control designs for the first few vibration modes of a cantilever beam using piezoelectric sensors/actuators [9]. Choi have shown a design of position tracking sliding mode control for a smart structure [10]. Distributed controllers for flexible structures can be seen in Pourki [11].

A FEM approach was used by Benjeddou to model a sandwich beam with shear and extension piezoelectric elements [12]. The finite element model employed the displacement field of Zhang and Sun [13]. It was shown that the finite element results agree quite well with the analytical results. Raja et al. extended the finite element model of Benjeddou’s research team to include a vibration control scheme [14].

3. Mathematical Modeling

A cantilever slender beam with rectangular cross-sections is considered. Four pairs of piezoelectric patches are embedded symmetrically at the top and the bottom surfaces of the beam, as shown in Figure 1.

Figure 1: Smart beam.

The beam is from graphite-epoxy T300–976 and the piezoelectric patches are PZT G1195N. The top patches act like sensors and the bottom like actuators. The resulting composite beam is modelled by means of the classical laminated technical theory of bending. Let us assume that the mechanical properties of both the piezoelectric material and the host beam are independent in time. The thermal effects are considered to be negligible as well [15].

The beam has length , width , and thickness . The sensors and the actuators have width and and thickness and , respectively. The electromechanical parameters of the beam of interest are given in Table 1.

Table 1: Parameters of the composite beam.

Parameters
Values

Beam length,
0.3 m

Beam width,
0.04 m

Beam thickness,
0.0096 m

Beam density, ρ
1600 kg/m³

Young’s modulus of the beam,
 N/m²

Piezoelectric constant,
 m/V

Electric constant,
 V m/N

Young’s modulus of the piezoelectric element
 N/m²

Width of the piezoelectric element
  m

Thickness of the piezoelectric element
   m

In order to derive the basic equations for piezoelectric sensors and actuators [1], we assume that

(i)the piezoelectric sensors actuators (S/A) are bonded perfectly on the host beam;(ii)the piezoelectric layers are much thinner than the host beam;(iii)the piezoelectric material is homogeneous, transversely isotropic, and linearly elastic;(iv)the piezoelectric S/A are transversely polarized [1, 3, 16].

3.1. Finite Element Formulation

We consider a beam element of length , which has two mechanical degrees of freedom at each node: one translational (resp. ) in direction and one rotational (resp., ). The vector of nodal displacements and rotations is defined as [17]

The beam element transverse deflection and the beam element rotation of the beam are continuous and they are interpolated within by Hermitian linear shape functions and as follows [18, 19]:

This classical finite element procedure leads to the approximate discretized variational problem. For a finite element, the discrete differential equations are obtained by substituting the discretized expressions into the first variation of the kinetic energy and strain energy [18, 20] to evaluate the kinetic and strain energies. Integrating over spatial domains and using the Hamiltons principle [20], the equation of motion for a beam element is expressed in terms of nodal variable as follows:

where is the generalized mass matrix, the viscous damping matrix, the generalized stiffness matrix, the external loading vector, and the generalized control force vector, produced by electromechanical coupling effects. The independent variable is composed of transversal deflections and rotations ; that is, [16, 18]

where is the number of nodes used in analysis. Vectors and are positive upwards. To transform to state-space control representation, let (in the usual manner)

Furthermore, to express as , we write it as , where the piezoelectric force is for a unit applied on the corresponding actuator and represents the voltages on the actuators. Furthermore, is the disturbance vector [16, 18].

Then,

The previous description of the dynamical system will be augmented with the output equation (displacements only measured) [17] as follows:

In this formulation, is (at most, but can be smaller), while d is . The units used are compatible for instance m, rad, sec, and N [21, 22].

4. Linear Quadratic Regulator: LQR Control

It is well known [23, 24] that constant input disturbances can be eliminated only if the controller has a zero at infinity (i.e., it integrates). Another useful interpretation is that an integrator is a disturbance estimator. Hence we do not expect a zero steady-state error using an LQR controller.

The structure of LQR control with reduced order observer is shown in Figure 2.

Figure 2: LQR controller with state estimator.

Here, are the disturbances, is the noise, and the controller defines,

where

minimizes the weighted performance index as follows:

and , and are design weight matrices.

The controller is given by relation

which is the solution of the Riccati equation.

The weight matrices and are used in order to:

(i)normalize the state and control vector;(ii)assess the relative influence of deflection from equilibrium position and magnitude of control on the determination of a global criterion. Matrices and are diagonal with positive diagonal inputs, so that

The value sets the maximum desirable output value . The value has similar significance for input .

Matrix sets the weight for each state, while matrix holds the weight for each actuator’s voltage. The LQR problem requires that the state be known [23].

5. Inputs

A typical wind load (Figure 3) Acting on the side of the structure. The wind load is a real-life wind speed measurement in relevance with time that took place in Estavromenos of Heraklion, Crete. We transform the wind speed in wind pressure.

Figure 3: Wind load.

Loading corresponds to the wind excitation. The function has been obtained from the wind velocity record, through the relation

where velocity, density, and .

Moreover, in all simulations, random noise has been introduced to measurements at system output locations within a probability interval of . Due to small displacements of system nodal points, noise amplitude is taken to be small, of the order of . On the other hand, the signal is introduced at each node of the beam by a different percentage, that percentage being lower at the first node due to the fact that the beam end point is clamped.

6. Results of Application of LQR Control

The and that were used are

Since and (11), matrix is the design matrix. Its eigenvalues are chosen in such a way that the observer subsystem can be about two times faster than the observed system. The selected values for our simulation are

These values have been obtained by trial and error, given the poor numerical properties of the system. To find these values, we have used a robust pole computation algorithm included in MATLAB [25, 26].

The controller [Klqr] is given by relation (11) which is the solution of the Riccati equation, where and are respective state and control matrices of the system and and are weight matrices of the performance criterion (regulator) (14) and (15), respectively.

For the simulation, beam nodal displacements and rotations with and without control are displayed in Figures 4 and 5, while Figure 6 presents actuator voltage values for control of all beam nodes.

(a)

(b)

(c)

(d)

Figure 4: Displacement at all beam nodal points, with and without LQR control.

(a)

(b)

(c)

(d)

Figure 5: Angle of rotation at all beam nodal points, with and without LQR control.

Figure 6: Control produced voltage at all beam nodal points with LQR control.

6.1. Discussion of the Results of the Linear Quadratic Regulator (LQR)

Using the linear quadratic controller criterion LQR, beam vibration reduction is observed at all nodal points, for both constant and sinusoidal mechanical input, as well as for realistic wind loading. LQR control achieves reduction of vibration but at the same time requires the entire system state time history as well as an extensive sensor distribution.

We encounter the following difficulties:

(i)system disturbances are unknown and unpredictable;(ii)the state vector is not measurable in its entirety, which in turn necessitates the use of an observer. This setup is problematic, as the observer has no information on the disturbance, which results in erroneous estimates. A way to circumvent this problem is the use of an unknown input observer. Unfortunately, this approach is not feasible, as one of the prerequisite conditions is not met. This situation complicates the problem, making the application of classic controllers such as LQR difficult, since its performance is directly related to the availability of the state vector, or at best of a reliable estimator of the state vector.

For the reasons mentioned above, we will continue with a discussion of more advanced control techniques for this particular problem, such as the control.

7. Control

To relate the structures used in classical and control, let us look at Figure 7, in the frequency domain [16, 23, 24].

Figure 7: control bloc diagram in the frequency domain.

In this diagram, all inputs and outputs of interest are included, along with their respective weighs , where , , , and are the weighs for the disturbances, control, noise, and outputs, respectively. The exogenous inputs are the noise and the disturbances . is the controller, , , , , define at the relation (6, 7, 8), and is the transfer function of our system.

To find the necessary transfer functions consider the following:

Combining all these gives

Note that the plant transfer function matrix, , is deduced from the suitably reformulated plant equations as follows:

where .

Hence,

The equivalent two-port diagram in the state space form is shown in Figure 8 for the close loop, and with more details in Figure 9, with

where are the output variables to be controlled and the exogenous inputs.

Figure 8: Two-port diagram.

Figure 9: Details of structure.

Given that has two inputs and two outputs, it is, as usual, naturally partitioned as

Also,

Using (18) the transfer function for which is

while the closed loop transfer function for is

or

Equation (25) is the well known lower LFT for .

To express in state space form, the natural partitioning

is used (where the packed form has been used), while the corresponding form for the controller is [27–29]

Equation (27) defines the following equations:

To find the matrices involved, we break the feedback loop and use the relevant equations.

Therefore the equations relating the inputs, outputs, states, and input/output to the controller are

From (30), we use ,    and take our initial state space equation in the form of (6, 7, 8), as follows:

Therefore the matrices are

7.1. Results with

Figure 10 presents the Bode diagrams of diagonal elements of the above weight matrices. These matrices have been obtained through a number of tests, to ensure the feasibility of finding a controller .

(a)

(b)

(c)

(d)

Figure 10: Bode diagrams of diagonal elements of weight matrices.

The controller obtained by applying control is 36 order. For this controller, . A plot of the maximum singular value of the weighted closed loop system (beam plus controller) is given in Figure 11, where we can clearly note that the value remains below at all frequencies.

Figure 11: Maximum singular value of the unweighted closed loop system.

Figures 12, 13, and 14 further show the maximum singular values of transfer functions of the unweighted closed loop system (i.e., the initial one) that are of interest.

Figure 12: Maximum singular value disturbance to error.

Figure 13: Maximum singular value noise to error.

Figure 14: Maximum singular value plots for control.

These figures show that the performance of the computed controller is satisfactory [30] since:

(i)as shown in Figure 12, there is a significant improvement in the effect of disturbance on error up to the frequency of 1000 Hz;(ii)as shown in Figure 13, there seems to be little effect of noise on error for frequencies beyond 1000 Hz;(iii)Figure 14 shows a satisfactory effect of the disturbance on the size of the control scheme (the design could be improved, if it were possible to reduce noise effect for frequencies of 1000 Hz).

To validate the above findings, system response time histories for the three input cases mentioned in this section are presented below.

Using the mechanical input, we get the following result.

Figure 15 shows the displacement time history at all nodal points of the beam, with and without control, while Figure 16 displays the angle of rotation time history at all beam nodal points, with and without control. By employing the control, vibration reduction is achieved; we observe vibration reduction of 90%. Figure 17 presents the time evolution of the produced actuator voltage, which turns out to be lower than the piezoelectric voltage limit value of 500 V.

(a)

(b)

(c)

(d)

Figure 15: Displacement at all beam nodal points, with and without control.

(a)

(b)

(c)

(d)

Figure 16: Angle of rotation at all beam nodal points, with and without control.

Figure 17: Control produced voltage at all beam nodal points with control.

7.2. Order Reduction of Controller

The controller found is of order 36. The fact that controller order, which is equal to the order of the system, is relatively higher than the order of classical controllers such as PI and LQR has led a number of researchers to develop order reduction algorithms. The most widely used such algorithm, known as HIFOO, has been implemented in a Matlab environment, and is the one used in the following procedure [31].

The general problem is to compute a controller of reduced rank/order while retaining the performance of the criterion as well as the behaviour of a full order controller for the given system [32, 33] as follows:

The state space equations for the controller are

Let be the spectral abscissa of a matrix , that is, the maximum real part of its eigenvalues. Then, we require not only that , where is the closed loop system matrix, but that as well. The feasible set of , that is the set of stable matrices, is not a convex set and has a boundary that is not smooth [34, 35].

The HIFOO procedure has two phases: stability and performance optimization [31, 36]. In the stability phase, HIFOO attempts to minimize

where is a positive parameter that will be described shortly, until a controller is found for which this quantity is negative; that is, the controller is stable and makes the closed loop system stable. In case it is unable to find such a controller, HIFOO terminates unsuccessfully.

In the performance optimization phase, HIFOO searches for a local minimizer of

where

The introduction of is motivated by the fact that the main design objective is to attain closed loop system stability and to minimize , by demonstrating that should be relatively small; the term , however, prevents the controller norm from becoming too large, in which case the stability constraint by itself would not exist. Given that it is preceded by the stability phase, the performance optimization phase is initialized with a finite value of . Consequently, when it reaches a value of for which , that value is rejected, since an objective reduction is sought at each iteration [31, 36].

7.3. Results Using Controller HIFOO

As mentioned before, the HIFOO controller is implemented in Matlab by way of appropriate routines. It is called in the following manner:

where plant is the system description in the form of (33) and is the controller order.

The resulting controller is described in state space in similar manner as ; that is,

The controller state space equation is given by (39), where controller matrices are equal to

For the purpose of comparison of HIFOO controller performance to that of , the beam free end response is examined, for the mechanical input.

For the input in Figure 18, the beam free end response is shown, initially with and then without the HIFOO controller, while Figure 19 presents produced actuator voltage using the HIFOO controller.

Figure 18: Beam free end displacement, with and without HIFOO control.

Figure 19: Stress at beam nodal points, using HIFOO.

Using the HIFOO controller for an actual wind loading, beam position control is effected with node displacements of order of , with lower produced voltage. We therefore maintain criterion performance with a lower order controller. The maximum produced voltage for the HIFOO controller is 7 V; the respective value is 45 V for the controller. In other words, beam adjustment to its equilibrium position is achieved with a lower order controller that requires lower voltage; see Figure 19.

8. Results

In the present work, the use of active control technology in intelligent structures has been presented. The goal of control is vibration reduction, while sustaining low steady state error, short recovery time, and small maximum uplift; at the same time, control energy must remain within operating limits.

The beam that was used was discretized using 1-dimensional finite elements with two degrees of freedom per node. Piezoelectric actuators were embedded in it with the objective of reducing vibrations under deterministic as well as stochastic loading conditions.

Initially, we examined the linear quadratic control criterion using a reduced rank observer, which makes the simulation more realistic. To find the observer, we employed a robust pole location algorithm. By selecting appropriate weights, beam vibration reduction was achieved for stochastic loading cases. In all simulations, random noise has been introduced in measurements, so that the system better approximate reality, given that displacement measurement by means of piezoelectric sensors is not reliable. Next we applied more advanced control techniques, such as the criterion. The controller found is of order 36.

In order to reduce computational requirements of the model, controller rank was reduced by means of non-parametric and nonconvex optimization, using the HIFOO controller. The controller exhibited good performance even for a significantly smaller system degree.

A natural consequence of the proposed research innovations is the acknowledgement of new scientific problems that can be used as the basis for further research beyond the scope of this work. The advantage of active control is the fact that it allows taking into account in the computation the worst case result of disturbances with uncertainty and system noise. Moreover, the active control can effectively cope with stronger input, permitting the design for a large frequency bandwidth. Results are noteworthy; vibration reduction is observed even for realistic wind loading, with piezoelectric component voltage kept within tolerance.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank the Wind Energy and Synthesis of Energy Systems Laboratory, of the Technological Educational Institute of Crete, for its offer of real life wind speed measurements in relevance with time that took place in Estavromenos of Heraclion, Crete.

References

B. Bandyopadhyay, T. C. Manjunath, and M. Umapathy, Modeling, Control and Implementation of Smart Structures, vol. 350 of Lecture Notes in Control and Information Sciences, Springer, 2007.

D. Halim, Vibration analysis and control of smart structures [Ph.D. thesis], 2002.

B. Culshaw, “Smart structures—a concept or a reality?” Journal of Systems and Control Engineering, vol. 206, no. 1, pp. 1–8, 1992.

S. S. Rao and M. Sunar, “Piezoelectricity and its use in disturbance sensing and control of flexible structures: a survey,” Applied Mechanics Reviews, vol. 47, no. 4, pp. 113–123, 1994.

T. Bailey and J. E. Ubbard Jr., “Distributed piezoelectric-polymer active vibration control of a cantilever beam,” Journal of Guidance, Control, and Dynamics, vol. 8, no. 5, pp. 605–611, 1985.

S. Hanagud, M. W. Obal, and A. J. Calise, “Optimal vibration control by the use of piezoceramic sensors and actuators,” Journal of Guidance, Control, and Dynamics, vol. 15, no. 5, pp. 1199–1206, 1992.

W.-S. Hwang and H. C. Park, “Finite element modeling of piezoelectric sensors and actuators,” AIAA Journal, vol. 31, no. 5, pp. 930–937, 1993.

B. Bona, M. Indri, and A. Tornambé, “Flexible piezoelectric structures-approximate motion equations and control algorithms,” IEEE Transactions on Automatic Control, vol. 42, no. 1, pp. 94–101, 1997.

W. Schiehlen and H. Schonerstedt, “Controller design for the active vibration damping of beam structure,” in Proceedings of the Smart Mechanical Systems Adaptronics SAE International USA, pp. 1137–1146, 1998.

S.-B. Choi, C.-C. Cheong, and C.-H. Lee, “Position tracking control of a smart flexible structure featuring a piezofilm actuator,” Journal of Guidance, Control, and Dynamics, vol. 19, no. 6, pp. 1364–1369, 1996.

F. Pourki, “Distributed controllers for flexible structures using piezoelectric actuators/sensors,” in Proceedings of the 32nd IEEE Conference on Decision and Control (CDC '93), vol. 2, pp. 1367–1370, San Antonio, Tex, USA, December 1993.

A. Benjeddou, M. A. Trindade, and R. Ohayon, “New shear actuated smart structure beam finite element,” AIAA Journal, vol. 37, no. 3, pp. 378–383, 1999.

X. D. Zhang and C. T. Sun, “Formulation of an adaptive sandwich beam,” Smart Materials and Structures, vol. 5, no. 6, pp. 814–823, 1996.

S. Raja, G. Prathap, and P. K. Sinha, “Active vibration control of composite sandwich beams with piezoelectric extension-bending and shear actuators,” Smart Materials and Structures, vol. 11, no. 1, pp. 63–71, 2002.

G. Zames and B. A. Francis, “Feedback, minimax sensitivity, and optimal robustness,” IEEE Transactions on Automatic Control, vol. 28, no. 5, pp. 585–601, 1983.

D. Halim and S. O. R. Moheimani, “Spatial resonant control of flexible structures—application to a piezoelectric laminate beam,” IEEE Transactions on Control Systems Technology, vol. 9, no. 1, pp. 37–53, 2001.

Z. Friedman and J. B. Kosmatka, “An improved two-node Timoshenko beam finite element,” Computers and Structures, vol. 47, no. 3, pp. 473–481, 1993.

G. E. Stavroulakis, G. Foutsitzi, E. Hadjigeorgiou, D. Marinova, and C. C. Baniotopoulos, “Design and robust optimal control of smart beams with application on vibrations suppression,” Advances in Engineering Software, vol. 36, no. 11-12, pp. 806–813, 2005.

B. Miara, G. Stavroulakis, and V. Valente, Topics on Mathematics for Smart Systems: Proceedings of the European Conference, Rome, Italy, 26–28 October 2006, World Scientific, Singapore, 2007.

A. Moutsopoulou, A. Pouliezos, and G. E. Stavroulakis, “Modelling with uncertainty and robust control of smart beams,” in Proceedings of the 9th International Conference on Computational Structures Technology (CST '09), B. H. V. Topping and M. Papadrakakis, Eds., paper 35, Civil Comp Press, Athens, Greece, September 2008.

H. F. Tiersten, Linear Piezoelectric Plate Vibrations, Plenum Press, New York, NY, USA, 1st edition, 1969.

S. M. Yang and Y. J. Lee, “Optimization of noncollocated sensor/actuator location and feedback gain in control systems,” Smart Materials and Structures, vol. 2, no. 2, pp. 96–102, 1993.

O. H. Bosgra and H. Kwakernaak, Design Methods for Control Systems, Course Notes, Dutch Institute for Systems and Control, 2001.

H. Kwakernaak, “Robust control and -optimization—tutorial paper,” Automatica, vol. 29, no. 2, pp. 255–273, 1993.

R. C. Ward, “Balancing the generalized eigenvalue problem,” SIAM Journal on Scientific and Statistical Computing, vol. 2, no. 2, pp. 141–152, 1981.

A. L. Tits and Y. Yang, “Globally convergent algorithms for robust pole assignment by state feedback,” IEEE Transactions on Automatic Control, vol. 41, no. 10, pp. 1432–1452, 1996.

H. Kimura, “Robust stabilizability for a class of transfer functions,” IEEE Transactions on Automatic Control, vol. 29, no. 9, pp. 788–793, 1984.

B. A. Francis, A Course on H∞ Control Theory, Lecture Notes in Control and Information Sciences, Springer, 1987.

J. C. Doyle, K. Glover, P. P. Khargonekar, and B. A. Francis, “State-space solutions to standard and control problems,” IEEE Transactions on Automatic Control, vol. 34, no. 8, pp. 831–847, 1989.

A. Packard, J. Doyle, and G. Balas, “Linear, multivariable robust control with a μ perspective,” Journal of Dynamic Systems, Measurement and Control, vol. 115, no. 2, pp. 426–438, 1993.

J. V. Burke, D. Henrion, A. S. Lewis, and M. L. Overton, “HIFOO—a MATLAB package for fixed-order controller design and optimization,” in Proceedings of the 5th IFAC Symposium on Robust Control Design (ROCOND '06), pp. 339–344, Toulouse, France, July 2006.

J. V. Burke, D. Henrion, A. S. Lewis, and M. L. Overton, “Stabilization via nonsmooth, nonconvex optimization,” IEEE Transactions on Automatic Control, vol. 51, no. 11, pp. 1760–1769, 2006.

J. V. Burke, A. S. Lewis, and M. L. Overton, “A robust gradient sampling algorithm for nonsmooth, nonconvex optimization,” SIAM Journal on Optimization, vol. 15, no. 3, pp. 751–779, 2005.

J. V. Burke and M. L. Overton, “Variational analysis of non-Lipschitz spectral functions,” Mathematical Programming B, vol. 90, no. 2, pp. 317–351, 2001.

D. Hinrichsen and A. J. Pritchard, Mathematical Systems Theory I: Modelling, State Space Analysis, Stability and Robustness, vol. 48 of Texts in Applied Mathematics, Springer, 2006.

M. Millston, HIFOO 1.5: structured control of linear systems with a non trivial feedthrough [M.S. thesis], Courant Institute of Mathematical Sciences, New York University, New York, NY, USA, 2006.