Creating Nodes at Selected Locations in a Harmonically Excited Structure Using Feedback Control and Green ’ s Function

)e paper deals with designing a control force to create nodal point(s) having zero displacements and/or zero slopes at selected locations in a harmonically excited vibrating structure. It is shown that the steady-state vibrations at desired points can be eliminated using feedback control forces. )ese control forces are constructed from displacement and/or velocity measurements using sensors located either at the control force position or at some other locations. Dynamic Green’s function is exploited to derive a simple and exact closed from expression for the control force. Under a certain condition, this control force can be generated using passive elements such as springs and dampers. Numerical examples demonstrate the applicability of the method in various cases.


Introduction
In recent years, there has been a considerable interest in modeling and controlling flexible structures. is is due to the use of lightweight materials for the purposes of speed and fuel efficiency.ese structures are characterized by low internal damping leading to prolonged vibrations, resulting in problems such as human discomfort, component failure, performance degradation, noise, and many other problems.Moreover, the performance is substantially degraded in some structures that are equipped with sensitive elements due to the occurrence of vibration.erefore, it is desirable to eliminate the vibration from a certain part of the structure more than the other.For example, large flexible space structures are usually built from lightweight materials having low damping.Any excitation can cause vibrations that may propagate throughout the whole structure.In such a case, it is desirable to confine the vibrations in some chosen insensitive part while keeping other parts, for instance, an extremely sensitive antenna, relatively undisturbed.
ere are a number of ways in which the vibration in a flexible structure can be eliminated or minimized.e first approach is to reduce or eliminate the external disturbance that excites the structure.Although this is always the preferred counter measure, but in many cases, there is little that can be done to change the nature of the driving force.e second approach is to modify the system to avoid vibrations, especially, during resonance.
is may entail significant design, and it may be difficult to be applied for existing structures.A third applicable approach is to apply one of the vibration control approaches that can either absorb the vibrational energy or load transmission path for the disturbing vibration.In conventional terminology, passive, active, hybrid, and semiactive are listed in the literature as four common classifications of vibration control methods.Passive control involves some form of structural modification or redesign, often including the use of springs and dampers that leads to the mitigation of vibrations.An active control scheme augments the structure with sensors, actuators, and some form of an electronic control system to generate active force that can be controlled with different algorithms to make the system more responsive to disturbances.e hybrid approach is a combination of active and passive ones intended to reduce the amount of external power necessary to achieve system control in comparison with the potentially large control effort required in the active scheme.
e semiactive (also known as adaptive-passive) method integrates tuning control schemes with tunable passive devices by replacing the active force generators with modulated variable components such as variable rate damping and stiffness elements.
ere is a wealth of literature dealing with vibration control using the aforementioned approaches in the form of an vibration absorber, which is a device originally invented by Frahm [1].He has registered a US patent in 1909 for "Device of Damping Vibrations of Bodies." e first mathematical theory for a dynamic vibration absorber (DVA) was presented by Ormondroyd and Den Hartog [2].Later, Brock [3] derived an analytical expression for the optimal damping ratio of the DVA.e working principle of the DVA did not change since Ormondroyd and Den Hartog's study [2].However, different designs have been developed that suit different applications.ere are many interesting practical applications described by Den Hartog [4] and Hunt [5].A comprehensive study on the theory and practice of vibration absorbers is given in the book by Korenev and Reznikov [6].An excellent survey of passive, semiactive, and active DVAs was prepared by Sun et al. [7].A review of adaptive vibration absorbers, which can adjust their own parameters, can be found in the work of von Flotow et al. [8].
e first researcher that applied the principle of the DVA to a continuous system (i.e., beam) is Young [9].Snowdon [10] considered the optimization of a discrete absorber on beams with various boundary conditions when structural damping was present.Jacquot [11] has established an analogy between a beam and a SDOF system that allowed the use of the optimal parameters for the latter to determine the corresponding for the beam.Özgüven and Candir [12] have investigated the use of two DVAs to suppress the resonant vibrations associated with the first two natural frequencies.Manikanahally and Crocker [13] extended the work in reference [12] to include point masses in their analysis.Different applications of vibration absorbers can be found in references [14][15][16].
e problem of a vibrating beam with sprung masses falls in the class of the dynamics of combined systems, which consist of linear elastic structures carrying lumped attachments.e free vibrations of such systems are studied extensively in the literature [17][18][19][20].Literature review indicates that the authors of these papers have generally directed their investigations into finding the natural frequencies and the corresponding mode shapes to evaluate the free vibrations.
e most common approach used is the assumed mode method which is one of the approximate methods best suited for solving diverse problems in dynamics and structures (for example, Ginsberg [21]).Other approaches of finding the response of dynamic systems such as finite element, Laplace transform, transfer matrix, and Lagrange multiplier are also used.Cha and Ren [22], utilizing an assumed mode method, used a chain of oscillators to impose nodes (i.e., points of zero displacements) to the normal modes of a beam by properly selecting values for the oscillator parameters (i.e., masses and spring constants) and by properly choosing the attachment location for the oscillator chain.e desired nodes can either coincide with the location of the oscillator chain, or it can be located elsewhere.Using the same method, Cha and Wong [23] analyzed a freely vibrating structure with a series of attached sprung masses aiming to impose multiple nodes for any normal mode.Cha [24] used a set of sprung masses to impose a single or multiple nodes anywhere along a vibrating beam.Alsaif and Foda [25] have applied the dynamic green function to calculate the optimal values of masses and/or springs and their locations along a beam in order to confine the vibration at an arbitrary location.Similarly, Foda and Albassam [26] have exploited Green's function to develop an exact method of analysis for the steady-state response of a harmonically excited Timoshenko beam with attached springs and/or masses and having different boundary conditions with the objective of confining the vibrations along the segment of the beam.Cha and Zhou [27] have created nodes with zero displacements and zero slopes in a harmonically excited linear structure using a set of sprung masses and rotational springs.Cha and Chen [28] have used lumped masses to impose nodes in a harmonically excited flexible beam.
On the contrary, active control methods can be used to create nodal points in continuous systems.Ram [29] has considered the problem of creating a node at the end of an axially vibrating rod using displacement feedback control.He derived a formula for the feedback control gain as a function of the product of the system eigenvalues.is result is based on a relationship derived by Elhay et al. [30] relating the system eigenfunctions to a product of eigenvalues.Singh and Ram [31] have created a node in a vibrating beam using feedback control with a gain calculated from either the infinite product of eigenvalues or deflection data.e implementation of acceleration and displacement feedback control on a beam-type tunable vibration absorber has been considered by Alujevic et al. [32].ey have studied the stability of the feedback loop using the Nyquist criterion.
In the present work, a displacement and velocity feedback active control force is designed for the purpose of enforcing one node, two nodes, and one node with zero slope along any beam during harmonic excitation.e beam is modeled as an Euler-Bernoulli beam that included the effects of internal and external damping.A key feature of this work is that the problem is formulated in terms of dynamic Green's function, which is exact and straightforward.e method is chosen for its freedom from numerical accuracy when compared with the standard application of the modal superposition technique.Furthermore, the boundary conditions are embedded into Green's function of the corresponding beam, and it is not necessary to solve the free vibration problem in order to obtain the eigenvalues and the corresponding eigenfunctions, which are required for modal superposition solution.Equally important, this method exhibits appreciably greater computational efficiency when compared with other methods.e rest of the paper is organized as follows: Section 2 describes the formulation of the problem and the derivation for the feedback control gain(s) to create one node, two nodes, and one node with zero slope in a harmonically excited beam.Stability analysis and sensor noise rejection 2 Shock and Vibration are presented in Sections 3 and 4, while numerical results for the aforementioned cases are displayed in Section 5. Section 6 gives some general conclusions.

Problem Formulation
We consider the problem of transverse vibration of a uniform elastic beam of finite length L originally at rest, as shown in Figure 1. e beam can have different classical or unconventional boundary conditions at x � 0 and x � L. e beam is excited by a sinusoidal external force f ext (x, t) that causes it to vibrate.A control force u(x, t) is designed and applied to the beam in order to reduce or completely eliminate the induced vibrations.e lateral vibration of the beam is governed by the following equations [33]: where EI, ρ, A, and w(x, t) denote, respectively, the flexural rigidity of the beam, the density, the cross-sectional area, and the transverse deflection of the beam at point x and time t. e second term represents the external damping effect which is assumed of the viscous type due to friction with the surrounding and has a coefficient c 1 .e third term represents the internal damping effect which is assumed to be viscoelastic of the Kelvin-Voigt type with coefficient c 2 .
Since damping exists, we will use complex forms for the input force, control force, and the response of the beam.e external excitation force, f ext (x, t), is assumed to be given by where F 0 is the amplitude of the external force, ω is the excitation frequency, i � �� � −1 √ , and x f is the location of the excitation force.
e symbol δ() denotes the Dirac delta function.It is understood that the excitation force is given by the real part of the right-hand side of equation ( 2). e control force, u(x, t), which is applied at x � x a is given by where α 1 , α 2 , and x s are the displacement control gain, velocity control gain, and displacement and velocity sensor location, respectively.At this stage, one can recognize two types of controllers: the collocated type is when the control actuator location (x a ) coincides with the sensor location (x s ) and the noncollocated type is when the actuator and the sensor have different locations along the beam.
Assuming an applied external force with monochromatic frequency, given by equation ( 2), the solution of equation (1), using the separation of variables, is sought to have the following form: where W(x) is the spatial beam deflection at point x.Substitution of equation (4) into equation (1) reduces it to Once again, it is assumed that we are only interested in the real part of the spatial beam deflection W(x).
To make the analysis more general and the problem solution is applicable to a large combinations of system parameters, it is preferred to work with quantities that are nondimensional.erefore, the following nondimensional variables and coefficients are defined: where Substituting the variables and coefficients, in equations ( 6)-( 8), into equation ( 5) results in the nondimensional form given by It is to be noted that the prime denotes the derivative with respect to x and  δ � (1/L)δ.For the sake of convenience, the hat symbol on the dimensionless quantities can be omitted proceeding to find the solution of equation ( 9) using dynamic Green's function.Denoting G(x, u) as the dynamic Green's function for the stated problem, which is not yet known, the solution of equation ( 5) can be given by Shock and Vibration 3 Performing the integration using the properties of the Dirac delta function to obtain where G(x, u) is the solution of the differential equation where the coefficient q 4 is defined by e solution of equation ( 12) is reported in references [25,34] for conventional and nonconventional boundary conditions and will not be repeated here.For example, for a simply supported beam, Green's function takes the following simple form: where g(x, u) � sinh(qL)sin(qx)sin(qL − qu) − sin(qL)sinh(qx)sinh(qL − qu), (15) and g(u, x) is obtained by replacing x with u in g(x, u).Dynamic Green's function provides afundamental solution for a differential equation.is function is very useful when dealing with linear systems, for example, books by Roach [35,36].Green's function G(x, u) is sometimes called source function, impulse response function, or receptance.It is considered a two-point function of position and is independent of the forcing term.Green's function G(x, u) for a system is defined as the displacement of the system at location x which is due to a unit force applied at location u.It is a symmetric function, in which, G(x, u) is equal to G(u, x), where the latter represents the displacement at point u due to a unit force applied at point x.erefore, Green's function satisfies the Maxwell-Rayleigh reciprocity law.
To this end, one should note that equation ( 11) is valid at x � x s to give By solving equation (16) for W(x s ), one gets Substituting W(x s ), from equation (17), into equation (11), the steady-state deflection at any point x is, therefore, given by e expression for W(x), in equation (18), gives the exact closed form solution for the steady-state displacement field of a beam deflection.

One Node.
In this section, it is desired to calculate the control gains that result in creating a node at location x n on a beam.First, we find the deflection at a specific point x n along the beam by substituting x � x n in equation ( 18) to get If, in an application, it is required to have a node at a certain location on the beam, for example, x � x n , then one would set W(x n ) equal to zero, to give At this juncture, we utilize the fact that a complex number vanishes if its real and imaginary parts vanish independently.Equation ( 20) can be arranged and simplified to give x  Shock and Vibration To calculate the gains, α 1 and α 2 , we equate the real and imaginary parts of both sides of equation ( 21), to give where Re[] and Im[] denote the real and imaginary parts of [•], respectively.Equations ( 22) and ( 23) represent the explicit formulas for the control gains which ensure that the vibrational amplitude at x � x n is equal to zero. is result is better than the one reported by Singh and Ram [31] for an undamped beam because it is exact and does not require infinite product of the eigenvalues of the system.
Equations ( 22) and ( 23) represent the general case, in which, the control sensor and actuator are at different locations, i.e., x s ≠ x a , which refers to the noncollocated control case.On the contrary, the control gains for the collocated control case can be calculated by substituting x s � x a in equations ( 22) and (23).
Examining equations ( 22) and ( 23), one can observe that if the sensor location coincides with the node location, i.e., x s � x n , then the denominator of both equations become zero and the control gains approach infinity.is case is expected since the sensor is located at a stationary point where both displacement and velocity are equal to zero.erefore, it is always required not to place the sensor at the node location.Consequently, if it is desired to create a node at the control actuator location, i.e., x n � x a , then one should avoid placing the sensor at the actuator location and use the noncollocated control case.In addition, the control gains become infinity if the actuator location coincides with the excitation force location, i.e., x a � x f .Moreover, if the system is undamped (c 1 � c 2 � 0), then Green's function becomes a real valued function where α 2 , from equation (23), is equal to zero, which means that a node cannot be imposed by a control force that is proportional to both beam displacement and velocity at the sensor position.Instead, a node can be enforced at any location by using a control force that is proportional to the beam displacement only.A comparable situation is encountered in the theory of vibration absorber, where one cannot introduce antiresonance in a system consisting of an undamped primary system and a damped DVA as an auxiliary attachment.
It is to be noted that if the sensor location coincides with the actuator location and the control gains are negative, then the control action can be implemented using passive elements.In other words, the control can be realized by a spring of stiffness equals −α 1 and a damper with damping constant equals −α 2 , both are connected in parallel between the ground and the beam at the actuator location (x a ).

Two Nodes.
In this section, we demonstrate a numerical method to impose two nodes in a vibrating beam subjected to a harmonic force by applying two control forces at specific locations on the beam using Green's function.Our goal here is to determine the values of the control forces required.
e net control force u(x, t) formed by the summation of the two control forces applied at locations x a 1 and x a 2 is given by where x s 1 and x s 2 are the locations of the sensors for measuring the displacement and velocity amplitudes for the two control forces, respectively.Similar to the one node case, the nondimensional governing equation of the system takes the following form: Similar to equation (11), the solution of the aforementioned equation takes the form For notational simplicity and convenience, we drop the hat symbol in equation ( 26) and follow the same steps as in the one node case to calculate the control gains α 1 , α 2 , α 3 , and α 4 , which are as follows: (1) Substitute x � x s 1 and x � x s 2 in equation (26)  Shock and Vibration It is to be noted that equations (27a) and (27b) generate the two desired nodes on the beam.
After applying the above steps, one can obtain the following results: where W(x s 1 ) and W(x s 2 ) are the steady-state amplitude of the beam at the sensor locations x s 1 and x s 2 , respectively, and A, B, and C are given by and Solving equations (27a) and (27b) for B and C to get From equations (29) for B and C, the control gains are given by It is to be noted that this case represents a noncollocated case, in which the locations of the control sensors are different from the locations of the control actuators.e control force gains for the collocated case can be easily obtained by substituting x s 1 � x a 1 and x s 2 � x a 2 in equations (30a) and (30b).

One Node with Zero Slope.
In this section, we demonstrate a numerical method to impose one node with zero displacements and zero slopes in a vibrating beam subjected to a harmonic force by applying two control forces at specific locations on the beam using Green's function.e advantage of having the slope of the displacement at the node to be zero is the improved robustness by increasing the interval of zero displacements near the node.Our goal here is to determine the values of the control forces required.
e equations involving u(x, t), the equation of motion, and its solution are given in equations ( 24)- (26) given in Section 2.2.
e following steps are applied in order to calculate the control gains α 1 , α 2 , α 3 , and α 4 : (1) Substitute x � x s 1 and x � x s 2 in equation (26) to form two equations (2) Solve for W(x s 1 ) and W(x s 2 ) from the two equations in Step 1 and substitute them back into W(x) in equation ( 26 Step 2 (4) Solve for the control gains α 1 , α 2 , α 3 , and α 4 from the following two equations: Equations (32a) and (32b) are responsible for generating a node at a desired location (x n ) on the beam having both zero displacements and zero slopes.Shock and Vibration After applying the above steps, one can obtain the equations for W(x s 1 ) and W(x s 2 ) which are the same as that given by equations (28a) and (28b).e solutions for B and C are given by where () ′ denotes differentiation with respect to x. e control gains, α 1 , α 2 , α 3 , and α 4 , can be calculated using equations (31a), (31b), and (29).
is case represents the noncollocated case, in which the control sensors are located in different locations from the control actuators.On the contrary, the gains for the collocated control case can be calculated by substituting x s 1 � x a 1 and x s 2 � x a 2 in equations (33a) and (33b).

Stability Analysis
It has been shown that the control gains to impose a node in a beam structure can be either positive or negative.If they are negative, then the control can be implemented by using passive elements.is is done by attaching to the ground a spring of constant k � −α 1 and a damper with damping constant c � −α 2 at the control force actuator location.
If the system is undamped, then it has been shown that a node can be imposed by using only displacement feedback with a control gain α.In this case, Singh and Ram [31] have shown that the controlled system is stable if and only if where d is the static deflection of the beam due to a unit control force applied at the control actuator location (x a ).Since d is always positive, the inequality (34) holds whenever the value of the control gain α is negative.In this case, it is concluded that the controlled system is stable.When the beam structure is damped, imposing a node requires velocity feedback signal with gain α 2 in addition to the displacement feedback signal with gain α 1 .It is assumed that α 1 satisfies the stability condition (34).If the value of α 2 is negative, then this is equivalent to attaching a damper with damping coefficient c � −α 2 at the location of the control actuator (x a ).Since the damping coefficient c is positive, the damper will dissipate beam energy, and, therefore, the controlled beam is stable.On the contrary, if the velocity feedback control gain α 2 is positive, then the stability can be examined by discretizing the equation of motion, equation (1), after substituting equations ( 2) and (3).
Using the assumed modes method, the transverse deflection of the beam w(x, t) can be approximated using where ϕ i (x) denotes the i th assumed mode shape, η i (t) denotes the i th generalized coordinate, and n denotes the number of terms retained in the approximation.e kinetic energy of the beam structure can be written as and the potential energy is given by e virtual work done by the external forces (nonconservative forces) can be expressed as where p(x, t) represents the external applied force per unit length of the beam as well as the damping effects and the control terms and is given by where all the terms are defined in the previous sections.By substituting equation (35) into equations ( 36)-(38) and performing all the integrations and applying Lagrange's equations, e following discretized equations of motion, in the matrix form, can be derived: where Shock and Vibration where M, C, K, C a , and K a are n × n matrices, Q is an n × 1 vector, () ′ denotes the differentiation with respect to x, () T denotes the transpose, and the n × 1 vectors ϕ and η are defined by Equation ( 41) can be transformed to a state-space form by defining the following states: to give the state equation where z is the 2n × 1 state vector, A is the 2n × 2n system matrix, B is the 2n × 1 input vector, and u * is a scalar input, which are given by e following comparison function ϕ i (x) is chosen as the assumed modes shape: which satisfies both the geometric and natural boundary conditions of a cantilever beam.e stability of the controlled beam structure can be evaluated by calculating the eigenvalues of the system matrix A. For a stable system, all the eigenvalues of the system matrix must have a negative real part and marginally stable if any of the eigenvalues have zero real parts.

Sensor Noise Rejection
In this section, the effect of sensor noise on the system response is analyzed by using the state-space model in equation ( 46). e states represent the generalized modal coordinates which are related to the beam deflection (w(x, t)) through equation (35).Since a sensor noise affects the measurements of the beam deflection, mathematically it will affect the values of the generalized modal coordinates.Consequently, to study the effect of sensor noise on the response and then try to eliminate its effect, a noisy input is added to the sinusoidal input of the state-space model to generate a noisy signal at the output.
e state-space model for a controlled cantilever beam is given by the state equations, equation (45), along with the output equation given by where y is an r × 1 output vector, C * is an r × 2n output matrix, and D is an r × 1 feedforward matrix.
To obtain the deflection of the beam at any location, the system represented by equations ( 45) and ( 48) is solved for the state-space vector (z(t)) due to a sinusoidal input.If the desired output of the state-space model is the deflection of the beam at the sensor location (w(x s , t)), then the 1 × 2n output matrix (C * ) and the 1 × 1 feedforward matrix (D) are given by To verify the accuracy of the state-space model, the beam deflection at the sensor location calculated using the statespace model is compared with that obtained from Green's function method.
A simple low-pass filter design is performed using Simulink.A Simulink block diagram is created that is comprised of a noisy input signal, the state-space model, a low-pass filter, and an output port.e input to the system is a sinusoidal input with an added random signal with Gaussian distribution.e low-pass filter is given by the first-order transfer function: where a is a constant selected based on the parameters of the input and the beam system and s represents the complex Laplace variable.is filter has a cutoff frequency of a rad/s.

Simulations and Discussion
In this section, different numerical experiments are carried out to demonstrate the applicability of the proposed approach for creating single node, two nodes, and a single node with zero slopes during harmonic excitations of beams.
Cases where active and passive means for vibration control, in addition to collocated and noncollocated control cases, are investigated, and system stability analysis is performed.Furthermore, the numerical technique is applied to beams with different end supports.e last section demonstrates a simple filter design using Simulink to eliminate the noisy signal produced by the sensor that can affect the performance of the controller.Nondimensional quantities, defined in equations ( 6)-( 8), are used in the numerical examples which can be 8 Shock and Vibration transformed to their dimensional counterpart once all the parameters of the beam structure are defined.

One Node. Consider a uniform cantilever beam of length
L � 1, excited by a concentrated harmonic force of amplitude F 0 � 1, an excitation frequency ω � 30, and excitation location x f � 1. e external damping coefficient c 1 � 0.1, while the internal damping coefficient c 2 � 0.001.For a given application, a node is desired at the location of the force, i.e., x n � x f .e steady-state lateral deflection of the beam is shown in Figure 2. e solid curve corresponds to the deformed shape of the beam with a control force acting at position x a � 0.55, which coincides with the location of the control sensor (x s ) (i.e., collocated).e dashed curve corresponds to the deformed shape of the beam with no control force, and the zero horizontal line represents the configuration of the undeformed beam.e control gains are α 1 � −317.159 and α 2 � 0.201, and the corresponding amplitude of the control force is |u| � 2.710.As seen in Figure 2, the point x n � x f remains stationary if the control force with the proper chosen gains is applied.e casual observer may doubt the results shown in Figure 2 arguing that since a node exists at the location of the excitation force, the deflection along the whole beam span should be zero because the work of the excitation force is zero.A similar situation occurs in the undamped DVA theory that can remove this suspicion.e excitation force that acts on the primary mass results in the vibration of both masses, namely, the primary mass and the absorber mass.However, when the tuning is satisfied, the primary mass becomes motionless, with no work being injected into the system, but the absorber mass keeps vibrating.e book by Ginsberg [21] has a good discussion on this subject.
e stability of the beam system in this example can be examined using the analysis in Section 3. e system matrix and the input vector, in equation (45), are computed using the dimensional parameters of the beam structure, EI � 5 N•m 2 , ρA � 3 kg/m, L � 1 m, and F 0 � 1 N. e number of generalized coordinates used in the simulation is 100.
e eigenvalues of the system matrix are computed using Matlab.All the calculated eigenvalues have negative real parts proving the controlled beam to be stable.e time simulation for all the generalized coordinates (η i (t)) due to a unit impulse input applied at the excitation force location (x f ) is shown in Figure 3.
Figure 4 shows the design plots of the required control gains and the amplitude of the control force versus the location of the control actuator (x a ) for the same frequency, damping, excitation force location, and node location associated with Figure 2. Considerable information related to stability, optimality, and type of control can be drawn from this figure.Here the optimality is based on the minimum control force requirement.e control gains near x a � 0 are very large because the displacement amplitude W(x a ) is almost equal to zero.Moreover, the control force at x a � 0 is equal to zero since the displacement amplitude at the fixed end is equal to zero.On the contrary, the control gains which are needed to achieve no steady-state motion at the free end are infinite because the displacement is required to be zero, but the control force amplitude has to be equal to −1.Furthermore, it is noted that the feedback control gain, α 1 , experiences local maximum and minimum, whereas the feedback control gains, α 2 , experience a local maximum at the two control force locations, x a � 0.34 and x a � 0.82. is is due to the coincidence of the control force location with a node created in the controlled beam deflection.
It is to be noted that the control gains α 1 and α 2 can have positive or negative values depending on the system and the control parameters.Consequently, the control force for which the control gains are negative can be implemented by using passive elements, i.e., springs and dampers.
e issue of robustness is more critical than the issue of optimality of the control force.is is clearly obvious because at the free end, the control force requirement is a minimum but the control gain is infinity.e range for the control gains, along the beam, having a flat shape is a sign of robustness.
Figure 5 shows the required control gains and the amplitude of the control force in terms of the excitation frequency for the same system parameters given in Figure 2. It is noted that the control gains, α 1 , experience local maximum and minimum, whereas the control gains, α 2 , experience a local maximum at the two excitation force frequencies, ω � 32.83 and ω � 144. is is due to the coincidence of the control force location with a node created in the controlled beam deflection at these two frequencies.Furthermore, the amplitude of the control force is almost zero when the excitation force frequency is ω � 15.42 because the controlled and uncontrolled beam deflection coincide with each other at this frequency and a node is created at the free end without the need for a control force.Similar behavior occurs at the excitation force frequency ω � 50.In conclusion, and for a given excitation force frequency, node location, and damping coefficients, one should place the control force actuator that results in minimizing the control force magnitude.
e steady-state deformed shape of a uniform undamped cantilever beam subjected to a localized harmonic force applied at x f � 0.5 with a forcing frequency ω � 61.68, which is near the beam third nondimensional natural frequency, is shown in Figure 6.It is desired to have a node at x n � 0.8 with a control actuator located at x a � 0.3.e control gains are α 1 � 20.347 and α 2 � 0, and the corresponding amplitude of the control force |u| � 0.0267.e beam with control force has a node at exactly x n � 0.8, and it experiences substantially less vibration compared to the uncontrolled beam.Since the value of the control gain is positive, the stability of the controlled beam needs to be examined using condition (34).e control gain (α 1 ) and the inverse of the beam static deflection due to a unit control force applied at the control actuator location (x a ) versus the control actuator location (x a ) is shown in Figure 7 and its zoomed counterpart in Figure 8.It is apparent from Figure 7 that the actively controlled system is stable when the position of control force is between 0 ≪ x a ≪ 0.47, 0.50 ≪ x a ≪ 0.72 and 0.80 ≪ x a ≪ 0.94.Furthermore, it is clear from Figure 8 that the value of the control gain (α 1 ), at the control actuator Shock and Vibration location (x a � 0.3), is less than the value of the inverse of the static deflection (1/d), and therefore, the controlled beam system is stable.
In certain applications, it is desired to place the actuator at the required node location, i.e., x a � x n .In this case, it is compulsory to place the sensor at a location x s that does not coincide with x a .
is situation is shown in Figure 9, in which a node is desired to be created at the control actuator location, namely, x n � x a � 0.4, in an undamped simply supported beam.e sensor position is at x s � 0.75 and the excitation force frequency ω � 42.
e required control gains are α 1 � −269.566 and α 2 � 0, and the corresponding amplitude of the control force is |u| � 0.8533.
Figure 10 displays the steady-state deformed shape of an undamped simply supported beam subjected to a concentrated harmonic excitation force with a frequency ω � 246.73 which nearly coincides with nondimensional fifth natural frequency of the uncontrolled beam given by 25π 2 .e external force acts at x f � 0.77, and the node is required to be imposed at the excitation force location, i.e., x n � 0.77.A control force which is proportional to the displacement is applied at x a � 0.9 to achieve the aforementioned requirement.e required control gains are α 1 � −1507.459and α 2 � 0, and the corresponding amplitude of the control force is |u| � 0.45355.
is figure indicates that the controlled beam is practically motionless.e underlying reason for this behavior is that if the system is driven at resonant frequency, only one mode that represents a flexural wave contributes to the beam deformation.
erefore, if the control force is applied, it creates an antiwave so that when it is added to the original flexural wave, the result is no vibration.
In comparison with other numerical methods to impose a node with zero displacements, Example (1) of Singh and Ram [31] is considered.
ey considered an undamped Euler-Bernoulli cantilever beam of length L � 2, EI � 5, ρA � 3, F 0 � 2, ω � 10, and x f � 2 all with appropriate units.e control force location x a � 0.75.eir objective is to determine the feedback control force gain needed to eliminate the motion at the free end, i.e., x n � 2. e  Shock and Vibration measurement location (x s ) is at 0.75 (collocated case).ey used a numerical method based on deformation that divides the beam into two parts and solves two differential equations for the beam deflection amplitudes to obtain the control gain (α).eir other method of solution is to express the control gain in terms of an infinite product of eigenvalues.e calculated control gain is α � −1781.5.Using Green's function, the above parameters are transformed into a nondimensional form using equations ( 6) and (7)

Two Nodes.
We now turn to demonstrating the effectiveness of the feedback control design method using Green's function to create two nodes in a vibrating beam using two control forces receiving displacement and velocity measurements from sensors that are either collocated or noncollocated with the control actuator.Figure 12 shows the steady-state deflection of a damped cantilever beam with a concentrated harmonic force applied at x f � 1 with frequency ω � 22. e external damping coefficient c 1 � 0.1, while the internal damping coefficient c 2 � 0.001.It is desired to introduce two nodes at locations x n 1 � 0.3 and x n 2 � 0.8.To accomplish this task, two control forces with collocated sensors are applied at x a 1 � 0.6 and x a 2 � 0.9 with calculated displacement and velocity feedback gains α 1 � −1577.876and α 2 � −1.515 for the first control force and α 3 � −2374.132and α 4 � −2.192 for the second control force.e corresponding amplitudes of the control forces are |u 1 | � 1.463 and |u 2 | � 0.396.If, in certain applications, one wishes to place the actuators at the required nodes locations, i.e., x a 1 � x n 1 and x a 2 � x n 2 , then the sensors should be placed at locations x s 1 and x s 2 different from the actuators locations x a 1 or x a 2 .is situation is illustrated in Figure 13, in which a simply supported beam is excited by a concentrated harmonic force that acts at x f � 0.5 with the excitation frequency ω � 10, which is very close to the beam nondimensional first natural frequency (π 2 ).It is required to create two nodes at x n 1 � 0.6 and x n 2 � 0.8 using two control forces with actuators located at the same locations as the nodes, i.e., x a 1 � 0.6 and x a 2 � 0.8.e two sensors, feeding back displacement signals to the two control force actuators,

One
Node with Zero Slope.In this section, numerical experiments are carried out to illustrate the applicability of the approach in Section 2.3 to impose a single node with both zero deflections and slopes using two control forces.While the zero deflections force a point on the beam to be stationary, the zero slopes create a flat region near the node having zero deflections, which makes the control design more robust.
Figure 14 displays the steady-state displacement of a controlled (solid line) and uncontrolled (dashed line) damped cantilever beam excited by a harmonic concentrated force applied at x f � 1 with a frequency ω � 20.It is desired to impose a node, having both zero deflections and slopes, at the beam location x n � 0.3, using two collocated control forces applied at x a 1 � 0.95 and x a 2 � 0.9.e external and internal damping coefficients for the beam are c 1 � 0.  Figure 12: e steady-state deformed shape of a collocated controlled (solid line) and uncontrolled (dashed line) damped cantilever beam when c 1 � 0.1, c 2 � 0.001, ω � 22, x f � 1, x n 1 � 0.3, x n 2 � 0.8, x a 1 � 0.6, and x a 2 � 0.9.e inset shows the detailed plot of the controlled shape.
14 Shock and Vibration case, it is assumed that a node is required to be enforced at the control force actuator location.In this case, the sensor for that actuator should be placed at any other different location from that of the actuator.is situation is illustrated in Figure 15, in which an undamped simply supported beam is excited by a concentrated harmonic force that acts at x f � 0.5 with excitation frequency ω � 10, which is very close to the beam nondimensional first natural frequency (π 2 ).It is required to create a node at the excitation force location, x n � 0.5, having both zero deflections and zero slopes using two control forces.e solid curve corresponds to the deflected shape of the beam with two control forces acting at locations x a 1 � 0. To show the computational efficiency of the proposed numerical method, an example reported in Cha and Zhou [27] is solved using the proposed numerical method.In Section 3.2 of their paper, they consider an undamped simply supported beam excited by a concentrated harmonic force with frequency ω � 57 and an excitation location x f � 0.77.eir goal is to impose a node with zero slopes at x n � 0.31 using tuned translational and rotational oscillators attached to the beam at x a � 0.65.All the previous parameters are nondimensional.ey have shown the simulation of the steady-state lateral displacement of the beam in Figure 7 of their paper.e simulation contains three curves: the solid line corresponds to the deformed beam with sprung mass and rotational oscillator, the dashed dotted line Figure 14: e steady-state deformed shape of a collocated controlled (solid line) and uncontrolled (dashed line) damped cantilever beam when c 1 � 0.1, c 2 � 0.001, ω � 20, x f � 1, x n � 0.3, x a 1 � 0.95, and x a 2 � 0.9.e inset shows the detailed plot of the controlled shape.

Shock and Vibration
corresponds to the deformed beam with sprung mass only, and the dashed line corresponds to the deformed beam without any oscillator.To simulate this problem, they have used 1000 terms in the assumed modes method and have to solve highly nonlinear simultaneous equations using the Matlab routine "fsolve" to search for a solution.e same example is solved using the exact analytical solution outlined in Section 2.3 to impose a node with zero displacements and zero slopes using collocated feedback control case, with x a 1 � 0.65 and x a 2 � 0.55.First, equations (33a) and (33b) are solved for B and C, and then the feedback control gains are calculated using equations (31a) and (31b).
e computations are performed using Matlab, and the execution time is about 2 seconds.e calculated displacement control gains for the two control forces that are responsible for imposing a node with zero slopes are −1.658× 10 4 and −7.621 × 10 14 , whereas that for imposing a node is −1.899 × 10 3 .e steady-state responses for this case are shown in Figure 16.For comparison, the figure has three curves: the solid line corresponds to the deformed beam having a node with zero slopes, the dashed dotted line corresponds to the deformed beam having one node, and the dashed line corresponds to the deformed beam without feedback control.e two figures are comparable.
e formulation presented in this paper is applicable to any beam structure under conventional and nonconventional boundary conditions.is is due to the fact that Green's function derivation for a specific beam includes its boundary conditions, and, therefore, they are embedded in the derived Green's function for a specific beam.erefore, the proposed feedback control design method is demonstrated only for cantilever and simply supported Figure 16: e steady-state deformed shape of a collocated controlled beam deflection having a node with zero slopes (solid line), a node (dashed dotted line), and uncontrolled beam deflection (dashed line) for the undamped simply supported beam when ω � 57, x f � 0.77, x n � 0.31, x a 1 � 0.65, and x a 2 � 0.55.16 Shock and Vibration beams in the previous examples, and other types of beams are not considered for the sake of brevity.

Sensor Noise Simulation.
e simulation of the procedure outlined in Section 4 is demonstrated on a damped cantilever beam system having a mass per unit length of 1 kg/ m, flexural rigidity of 5 N•m 2 , and length of 1 m. e beam is excited by a sinusoidal input of sin(ωt) with frequency equals to 5 rad/s.It is desired to impose a node at location of 0.9 m using a control force with a collocated actuator and sensor at 0.95 m. e coefficients of external and internal damping are 0.1 N•s/m 2 and 0.001 N•s•m 2 , respectively.e calculated displacement and velocity control gains are −9216.167N/m and −2.229 N•s/m, respectively.
e state-space model outlined in Sections 3 and 4 is developed in Matlab using n � 30 states that represent the generalized coordinates η i (t), i � 1, 2, . . ., 30. e calculated state matrix and input matrix have sizes of 60 × 60 and 60 × 1, respectively.e output variable (y(t)), considered in this example, is the beam deflection at the sensor location (w(x s , t)) so that the 1 × 60 output matrix and the feedforward matrices are given in equation (49).
To verify the accuracy of the state-space model, the beam deflection at the sensor location is calculated using Green's function and compared with that obtained using the statespace model and is shown in Figure 17. is figure shows an exact agreement for the beam deflection obtained by the two methods, which ensure the accuracy of the state-space model.
e sensor noise effect on the feedback control design using the state-space model is simulated using the Simulink block diagram shown in Figure 18. e sinusoidal excitation force with amplitude one and frequency 5 rad/s is generated by the sine wave block.To simulate a sensor noise, a random noise signal with the Gaussian distribution having a mean of zero and variance of 0.1 is created by the random source block.
Both signals are added to form the input u * to the state-space model.e state-space block contains the state, input, output, and feedforward matrices derived based on the aforementioned parameters having 30 states.e sensor noise filter transfer function, in equation ( 50), is represented by the lowpass filter block having a cutoff frequency equal to 10 rad/s.e noisy excitation force input to the beam system represented by the state-space model is shown in Figure 19, while the unfiltered controlled beam deflection at the sensor location represented by the output of the state-space model is shown in Figure 20.e filtered sensor output represented by the output of the low-pass filter is shown in Figure 21.e lowpass filter is capable of eliminating the effect of noise on the sensor output as evidenced by comparing Figures 20 and 21.

Conclusions
In this paper, a numerical method is developed to eliminate the steady-state vibrations at a prescribed point or points in a harmonically vibrating structure.ese points having zero displacements and/or zero slopes are referred to as nodes.
e paper makes use of the dynamic Green's function to analyze the problem.
e approach retains the infinitedimensional model for the structure to derive the exact solution for the steady-state deflection in terms of algebraic equations.
erefore, the method is free from numerical inaccuracy when compared with modal superposition techniques.e main advantages of this approach are its generality and its computational simplicity.e nodes are created by applying feedback control forces.ese control forces are constructed from displacement and/or velocity measurements using sensors located either at the control force position or at some other points.It is shown that if the system is damped, then the control force is proportional to the displacement and velocity, whereas it is proportional to the displacement only if the system is undamped.A closed form expression for the control gain(s) for three cases is Shock and Vibration  18 Shock and Vibration presented by utilizing dynamic Green's function.ese cases are creating one node, two nodes, and one node with zero slope at desired location(s) in a beam.e condition to create the node(s) using active or passive control is established.Several numerical examples are presented that show the effectiveness of the proposed control method.In the examples, it was shown that when the beam is excited near one of its natural frequency, then the feedback control force to create one node can lead to a motionless controlled beam.
is is due to the fact that only one mode dominates the flexural wave contributing to the beam deformation and the control force creates antiwave that cancels the original flexural wave.Furthermore, when the locations of the two nodes are closely spaced, a motionless region on the beam can be created.
e stability of the feedback control system can be guaranteed provided that the gain of the control force satisfies a stability condition (34) for the undamped system.e stability of the damped beam system has been examined using the eigenvalues of the system matrix.Furthermore, a sensor noise rejection model is implemented using a first-order filter and its effect on the response of the feedback control system is demonstrated using a Simulink block diagram.

Figure 1 :
Figure 1: Configuration of a harmonically excited cantilever beam with vibration control force.

Figure 4 :Figure 5 : 2 .Figure 6 :
Figure 4: e design plot of the required control gains and amplitude versus control actuator location.

1 and c 2
� 0.001, respectively.e calculated control gains are α 1 � −3.444 × 10 15 and α 2 � −7.066 × 10 13 for the first control force and α 3 � −96029.714and α 4 � −95.993 for the second control force.e corresponding amplitudes of the two control forces are |u 1 | � 1.022 and |u 2 | � 1.884.As seen in the figure, the beam with control forces has a node exactly at x n � 0.3 in addition to creating flat region having zero deflections over the whole beam span.Moreover, the controlled beam experiences vibrational amplitudes far less than the uncontrolled one compared to its uncontrolled counterpart.To demonstrate the noncollocated feedback control
3 and x a 2 � 0.4.e two sensors, feeding back displacement signals to the two control force actuators, are placed at the locations x s 1 � 0.6 and x s 2 � 0.7.e required control gains are α 1 � −5.630 × 10 15 and α 2 � 0 for the first control force and α 3 � −7.587 × 10 15 and α 4 � 0 for the second control force, and the corresponding control force amplitudes are |u 1 | � 4.831 × 10 12 and |u 2 | � 1.264 × 10 12 .e dashed curve corresponds to the deflected shape of the beam with no control forces.is figure indicates that the controlled beam is practically motionless.

Figure 17 :
Figure17: Beam deflection at sensor location (w(x s , t)) simulated using Green's function and state-space models.
Solve for W(x s 1 ) and W(x s 2 ) from the two equations (26)ep 1 and substitute them back into W(x) in equation(26)(3) Substitute x � x n 1 and x � x n 2 in the equation of W(x) formed in Step 2 to form the following two equations to solve for the control gains α 1 , α 2 , α 3 , and α 4