Abstract
Hysteresis is highly undesired for the vibration control of piezoelectric beams especially in high-precision applications. Current-controlled piezoelectric beams cope with hysteresis substantially in comparison to the voltage-controlled counterparts. However, the existing low fidelity current-controlled beam models are finite dimensional, and they are either heuristic or mathematically over-simplified differential equations. In this paper, novel infinite-dimensional models, by a thorough variational approach, are introduced to describe vibrations on a piezoelectric beam. Electro-magnetic effects due to Maxwell’s equations factor in the models via the electric and magnetic potentials. Both models are written in the standard state-space formulation (A, B, C), and are shown to be well-posed in the energy space by fixing the so-called Coulomb and Lorenz gauges. Different from the voltage-actuated counterparts, the control operator B is compact in the energy space, i.e. the exponential stabilizability is not possible. Considering the compact \(C=B^*-\)type state feedback controller (induced voltage), both models fail to be asymptotically stable if the material parameters satisfy certain conditions. To achieve at least asymptotic stability, we propose an additional controller. Finally, the stabilizability of infinite-dimensional electrostatic/quasi-static model (no magnetic effects) is analyzed for comparison. The biggest contrast is that the asymptotic stability is achieved by an electro-mechanical state feedback controller for all material parameters. Our findings are simulated by the filtered semi-discrete Finite Difference Method.
Similar content being viewed by others
1 Introduction
Piezoelectric materials have become more and more promising in aeronautic, civil, defense, biomedical, and space structures [8, 47] due to their small size and high power density. These materials can be used as both sensors and actuators since they have the ability to interchange mechanical energy, electrical energy, and magnetic energy during the motion. The amount of magnetic energy produced/stored can be minor or major in certain applications [55, Chap. 8]. Recent cutting-edge applications include cardiac pacemakers [7], course-changing bullets [10], structural health monitoring [17], nano-positioners [18], ultrasound imagers [45], ultrasonic welders and cleaners [51], energy harvesting [14] due to the excellent advantages of the fast response time, large mechanical force, and extremely fine resolution [18, 42].
There are mainly three ways to electrically actuate piezoelectric materials: supplying voltage, current or charge to its electrodes, see Fig. 1. Piezoelectric materials have been traditionally activated by a voltage source [4, 13, 41, 47, 48, 50]. In fact, it is stated in [1] that “From the very beginning of applying piezoelectric materials as actuators, they are voltage-controlled, and this is still the standard way of electrical steering.” However, the biggest disadvantage of voltage-controlled piezoelectric materials is that there is huge hysteresis between voltage source and mechanical strain [9]. An empirical approach to avoid hysteresis is simply applying only low voltages, but this prevents these structures from being used at their maximum potential. Another approach, which is gaining popularity [3], is current or charge-control designs since these types of actuation lead to considerably less hysteresis [5, 15, 24,25,26]. In fact, charge or current control enhances the dynamics, reduces the costs, and performs better at low frequencies [39]. This motivates the consideration of current-actuated piezoelectric beams.
Accurate modeling of piezoelectric beams are critical to achieve control objectives for high-precision motion [4]. However, the existing current-controlled beam models in the literature are generally based on ordinary differential equations, and these equations model only the first several vibrational models of devices. Moreover, majority of existing models considers only electrostatic or quasi-static assumptions due to Maxwell’s equations). Involving electrical effects in the modeling is done by coupling a circuit equation to the mechanical equations. This way, the current input is accounted for, see [19].
Even though it is known that the electrostatic theory is sufficient for many applications in piezoelectric acoustic wave devices, there are situations in which the full electro-magnetic coupling needs to be considered [21, 55, 56]. To our knowledge, there is no fully-dynamic infinite-dimensional model for “current-controlled” piezoelectric beams in the literature accounting for the full electro-magnetic effects. By incorporating mathematically-thorough Lagrangian formalism involving all electro-magneto-mechanical interactions, and using the variational approach [23], “partial” differential equations (PDE) can be used to accurately model vibrational dynamics of these devices. In this approach, the states of the electro-magnetic interaction due to the Maxwell’s equations are chosen in terms of the scalar electric potential and vector magnetic potential. However, it is well-known that the electric and magnetic potentials are not uniquely determined for this formulation (for instance see [11, p. 80]). To obtain a unique solution to the initial-boundary value problem, an appropriate gauge condition is chosen. Once a gauge condition is fixed, the electro-magnetic coupling in the charge equation is eliminated. There are many gauges used for Maxwell’s equations [12]. One of the most commonly-used one is called Coulomb gauge, and it assumes that electrical potential has infinite velocity. The other commonly-used one is called Lorenz gauge which assumes that the electrical field equation is a wave-type equation so electrical field propagates with the speed of light. In many applications, Coulomb gauge is preferred since proving uniqueness of the solutions is less challenging [29].
The controllability of piezoelectric beams is mostly studied for voltage control in finite or infinite dimensional setting [4, 9]. Letting \(\mathbf{y}, \mathcal A, B, u\) be the system state, system operator, control operator, and electric controllers, the PDE models can be transformed into the state-space formulation
The control operator B is unbounded in the infinite dimensional energy space if the piezoelectric material is controlled by a voltage or a charge source with [27, 28] or without [2, 20, 33, 47] magnetic effects. The controllability of electrostatic, quasi-static and fully dynamic models for voltage control is thoroughly studied in [27]. It is shown that only electrostatic/quasi-static model is exponentially stable by a \(B^*-\)type feedback controller. The fully dynamic model is not exponentially stabilizable by any state feedback, and not even asymptotically stable for many combinations of material parameters by the \(B^*-\)feedback controller [27, 30]. On the other hand, the control operator B is compact in the energy space if the piezoelectric material is controlled by a current source. This is first observed in [28, 37] where a piezoelectric beam model is considered in the standard state-space formulation with electromagnetic potentials and the Coulomb gauge. The current-driven piezoelectric beams have different controllability characteristics though since B is compact. Exponential stability is not possible. A lack of asymptotic stability is shown with the \(B^*-\)feedback controller provided that a number-theoretical condition is satisfied for certain subclass of material parameters. In this work, a thorough analysis for the description of the spectrum is not possible due to the complexity of coupling and the gauge condition. For the finite dimensional circuit-model setting, there are several results available in the literature [9, 19, 24, 25].
Following the modeling and control methodology of [27, 34, 35, 37], modeling and stabilization current and charge-controlled smart composites are also studied with the Coulomb gauge formalism in [32]. For the inclusion of electro-magnetic effects in piezoelectric layers, all three approaches: electrostatic, quasi-static, and fully dynamic, are considered. Even though a dramatic reflection of magnetic effects on the stabilizability characteristics of the composite with the \(B^*-\)type feedback is shown, the models provided alternate electro-magnetic feedback controllers. For the electrostatic and quasi-static approaches, exponential stability (charge control) and asymptotic stability (current control) results are proved.
There are also recent Port-Hamiltonian models of piezoelectric beams where the electro-magnetic modeling of the piezoelectric modeling is based on the models for the transmission lines. These models are automatically passive and with the extra dissipation these model can be made strictly passive. When these models are obtained, purely mechanical part can be interconnected to the piezoelectric part. In [54], it is shown that quasi-static Port-Hamiltonian modeling does not fulfill the necessary conditions of stabilization. In [52], fully dynamic model is obtained by blending the full set of Maxwell’s equations to the P-H structure for mechanical equations but without using the potential theory or gauge formalism. In [53], a nonlinear piezoelectric beam model is obtained with and without magnetic effects. This model is shown to fulfill the stability conditions with full mechanical controllers. Agreeing with the results of [27, 32], it is shown that the fully dynamic models of piezo-patch system with electro-magnetic and mechanical controllers, beams achieves stability [52].
To the best of our knowledge, there are no other results reported in the literature for the infinite dimensional models including all types of electro-magnetic effects. The fundamental goals of this paper are
-
(i)
to introduce fully-dynamic and electrostatic/quasi-static current-controlled PDE models obtained by a consistent variational approach. Novelty of the models is due to the inclusion of electro-magnetic effects by scalar electric and vector magnetic potentials.
-
(ii)
to prove that Lorenz gauge and Coulomb gauge models have well-posed state-space formulations in the corresponding energy spaces with compatibility conditions. In particular, the proposed Lorenz gauge formulation is quite new yet requires a special class of initial conditions for the electric potential. Analogous to [27, 37], both models provide alternative measurable electro-magnetic \(B^*-\)feedback controllers.
-
(iii)
to show that the closed-loop system for each model with the corresponding \(B^*-\)feedback controller lacks asymptotic stability for certain combinations of material parameters as in [32, 37]. Those combinations correspond to high-frequency vibrational modes. In this paper, we propose a remedial mechanical boundary controller to achieve at least asymptotic stability.
-
(iv)
to introduce the current-actuated electrostatic/quasi-static model which can never be exponentially stable. However, a feedback controller, requiring measurements of both mechanical and electrical states, can be designed to achieve polynomial stability. This type of control law has much to do with energy harvesting applications [14, 46].
-
(v)
to simulate vibrations with proposed control laws by the semi-discrete Finite Difference method with indirect filtering by adding relevant viscosity terms to each equation. This is a discretization technique where the so-called spill-over issue and the effects of spurious vibrational modes are discarded [49]. It is also safe to say that, to our knowledge, the numerical schemes using the filtering technique for a piezoelectric beam are preliminary and have never been proposed for the current-actuated piezoelectric beam models with the gauge formalism. However, further analysis is necessary.
Note that the compatibility condition arises for both gauges. The one found in the Lorenz gauge is quite new and it results from the semigroup formulation. Note that compatibility conditions are very similar to the “incompressibility condition” in hydrodynamics. This also lines up with the \(H(curl, \varOmega )\) and \(H(div, \varOmega )\) spaces in [11, 12] for which the global existence and uniqueness theorems are proved for Maxwell’s equations.
Note also that fully-dynamic modeling and stabilizability results for Coulomb gauge and free-free boundary conditions are first presented in [28, 37]. However, in this paper, we include the discussion of well-posed semigroup formulation with the Lorenz gauge, the model with the electrostatic/quasi-static assumptions, and the related numerical investigations for clamped-free boundary conditions. The numerical investigations are all based on the filtering technique developed in [49].
2 Modeling Assumptions and Gauge Conditions
Consider a piezoelectric beam of dimensions \(\varOmega =[0,L]\times [-r,r]\times [-\frac{h}{2}, \frac{h}{2}]\) where \(L> > h.\) Let x, y be the longitudinal directions, z the transverse direction, and \(\partial \varOmega \) the boundary. Since \(L>>h\), beam theory is appropriate and Euler-Bernoulli small displacement assumptions are good enough to model mechanical effects. To model the electro-magnetic properties, there are mainly three approaches [48]: electrostatic, quasi-static, and dynamic. Here, a fully dynamic approach is used. The Gauss’s law of magnetism \(\nabla \cdot B=0\) implies that there exists a magnetic potential vector A such that \(B=\nabla \times A\) by Poincáre’s lemma (p. 214, [12]). Substituting B into the Faraday’s law \(\dot{B}= -\nabla \times E\) yields the existence of a scalar electric potential \(\phi \) such that
To include the induced effects of the electromagnetism, we assume a quadratic-through thickness scalar electrical potential distribution and a linear-through thickness magnetic potential vector distribution:
Therefore, (1) is rewritten as the following:
We also assume that there are no mechanical external forces acting on the beam, and polarization and isotropy are taken in the transverse direction. Denoting v the stretching of the centerline of the beam in \(x-\)direction (longitudinal vibrations), the application of the Hamilton’s principle to the Lagrangian (kinetic energy-enthalpy+work due to \(i_s(t)\)) with respect to all possible admissible displacements yields two decoupled systems of equations; one set of equations for stretching another set of equations for bending. For details about modeling, refer to [28, 32, 37]. The applied current \(i_s(t)\) at the electrodes affects only the stretching motion. Therefore, only stretching equations are considered in this paper:
where \(\phi :=\phi ^1, \theta :=A_1^1,\) \(\eta :=A_3^0,\) \(\xi :=\frac{\varepsilon _{1}h^2}{12{\varepsilon _{3}}},\) and \(\rho , \alpha , \gamma , \varepsilon _1, {\varepsilon _{3}}, \mu \) denote the mass density per unit volume, elastic stiffness, piezoelectric coupling coefficient, permittivity in x and z directions, and magnetic permeability, respectively. It will also be assumed the surface electrical continuity condition [27, 50] is satisfied
since there is no external body charge or current applied to the beam. We consider the beam clamped on the left end, i.e. \(v(0)=0,\) and free on the right end. By the Hamilton’s principle, the natural boundary conditions are
The PDE system (5)–(7) does not yield a unique solution since (i) the magnetic potential vector components \(\theta ,\eta \) and the electric potential \(\phi \) are not uniquely defined due to (1) [11, 27], (ii) the Lagrangian is invariant under certain transformations [27, 37]. The ambiguity arising from this invariance must be fixed. To obtain a unique solution, particular gauge conditions are presented in electro-magnetic theory to completely decouple the electromagnetic equations in (5). Two of the most widely-used gauges are Lorenz and Coulomb gauges [12, 27, 29, 37, 44]. For the piezoelectric beam model, these gauges are
with the boundary conditions \(\theta (0)=\theta (L)= 0.\) The difference between these two gauges arise from assumption whether the electric potential \(\phi \) is static or dynamic.
Notice that the Coulomb gauge (8) makes the term \( -\xi {\dot{\theta }}_x+ {\dot{\eta }}\) in the \(\phi -\)equation in (5) zero. Therefore, \(\phi -\)equation in (5) is transformed into an elliptic equation:
Defining the operator \(P_\xi \) by
It can be shown that \({P_\xi }\) is a compact and non-negative operator on \({\mathcal {L}}_{2} (0,L)\). Therefore, the equation (9) has the solution \(\phi =\frac{\gamma }{{\varepsilon _{3}}}P_\xi v_x.\)
In the case of Lorenz gauge (8), the term \( -\xi {\dot{\theta }}_x+ {\dot{\eta }}\) in the \(\phi -\)equation in (5) is transformed into \(\frac{\xi {\varepsilon _{3}}}{\mu } \ddot{\phi }.\) As well, the terms \( {\dot{\phi }}_x -\frac{\mu }{\xi {\varepsilon _{3}}} \left( \eta _x- \theta \right) \) and \( {\dot{\phi }}-\frac{\mu }{{\varepsilon _{3}}} \left( \eta _{xx}-\theta _x\right) \) in \(\theta \) and \(\eta \) equations in (5) are transformed into \( -\frac{\mu }{\xi {\varepsilon _{3}}} \left( \xi \theta _{xx}- \theta \right) \) and \(-\frac{\mu }{\xi {\varepsilon _{3}}} \left( \xi \eta _{xx}- \eta \right) \), respectively. This transformation not only converts the \(\phi -\)equation to a wave equation but also the \(\theta \) and \(\eta \) equations. Therefore, both electric and magnetic equations are wave equations.
Taking into the argument above, the equations of motion (5)–(10) respectively reduce to
Note that distributed mechanical damping such as viscous damping \(\delta \dot{v}\) or Kelvin-Voigt damping \(-\delta \dot{v}_{xx}\) with \(\delta > 0,\) can also be considered in both models. However, the goal of the paper is to see the effect of the active electrical controller \(i_s(t)\) on the stabilization of the piezoelectric beam model.
3 Well-Posedness of the Lorenz Gauge Model (11)
The natural energy E(t) associated with (11) is the sum of kinetic, potential, magnetic, and electrical energies, respectively:
where \(\theta -\eta _x,\) \({\dot{\theta }}+\phi _x,\) and \({\dot{\eta }}+\phi \) refer to magnetic field component in \(y-\)direction, electrical field component in \(x-\)direction, and electric field component in \(z-\)direction. Therefore,
For any scalar \(C^1-\)function \(\chi =\chi (x, t),\) consider the transformation
Since
E(t) is invariant under the transformation (14). Once the Lorenz gauge (8) is used for \(\{{\tilde{\theta }},{\tilde{\eta }},{\tilde{\phi }}\}:\)
By (14), the boundary conditions for \(\chi \) is chosen in such a way that they align with the boundary conditions of \(\theta \) and \(\eta \) in (11)and the gauge condition (8). Therefore, \(\chi \) in (16) is considered together with the Neumann boundary conditions:
By (14), \({\tilde{\theta }} (x,0)= \theta (x,0)+ \chi _x (x,0),\quad {\tilde{\eta }}(x,0)= \eta (x,0)+ \chi (x,0),\quad {\tilde{\phi }} (x,0)= \phi (x,0)-{\dot{\chi }} (x,0).\) To obtain a zero solution for (16)–(17), we choose the zero initial conditions \(\phi \):
This implies that \(\chi (x,0)={\dot{\chi }}(x,0)=0.\) Therefore, the initial conditions for \(\{\theta ,\eta \}\) satisfy the following condition:
By (13), \(E(t)\equiv 0\) implies that
Finally, we obtain two sets of equations for \(\theta \) and \(\phi :\)
The solution of (19) is \(\phi \equiv 0.\) Feeding \(\phi \equiv 0\) into (20) yields \(\theta \equiv 0.\) Finally by (8), \(\eta =0.\)
Note that, in the case of a Coulomb gauge in Sect. 4, we obtain in Sect. 4 that \(\chi \) satisfies an elliptic equation, not a wave equation, with the same boundary conditions. The extra condition is not needed (18).
Considering the components of E(t) in (13), now define the states
with \( \mathbf{y}(x,0)=\mathbf{y}_0=[\theta (x,0)-\eta _x(x,0), {\dot{\theta }}(x,0), {\dot{\eta }}(x,0), v_x(x,0), \dot{v}(x,0) ]^{\mathrm{T}}.\)
By these choices of the states, note that the following compatibility condition arises
due to the second equation in (11) and the gauge condition (8); taking the time derivative of (8) and substituting it in the second equation in (11) and using (21) yields
Let \(H^1_L(0,L)=\{f\in H^1(0,L)~:~ f(0)=0 \}.\) Define the linear space
and the bilinear form on \(\mathrm H\times \mathrm H:\)
The term \({\varepsilon _{3}} y_3{\bar{z}}_3\) in (23) can be replaced by the compatibility condition in (23) so that
The bilinear form (24) is therefore symmetric, continuous and coercive on \(\mathrm H\times \mathrm H\). Continuity of \(a(\cdot ,\cdot )\) follows from the Poincare’s inequlaity for the \(y_2\)-term. The coercivity of \(a(\cdot ,\cdot )\) follows from the generalized Young’s inequality \(ab\le \frac{a^2}{\kappa } + \kappa b^2\) as the following
where \(\kappa \) is chosen as \( \frac{\gamma \xi }{\alpha +\frac{\gamma ^2}{{\varepsilon _{3}}}}<\kappa < \frac{\xi {\varepsilon _{3}}}{2\gamma }\) so that the coefficients \(\frac{\xi ^2{\varepsilon _{3}}}{2}-\gamma \xi \kappa , {\alpha }+\frac{\gamma ^2}{2{\varepsilon _{3}}}-\frac{\gamma \xi }{\kappa }>0,\) and \(C=\mathrm{min}(\mu , \xi {\varepsilon _{3}}, \left( \frac{\xi ^2{\varepsilon _{3}}}{2}-\gamma \xi \kappa \right) , \frac{{\varepsilon _{3}}}{2}, ({\alpha }+\frac{\gamma ^2}{2{\varepsilon _{3}}}-\frac{\gamma \xi }{\kappa }), \rho ).\)
Therefore, the bilinear form (24) defines an inner product on the linear space \(\mathrm H\). The following result is immediate:
Theorem 1
The energy E in (13) is the norm induced by this inner product, and \(\mathrm H\) is a Hilbert space with this norm.
Proof
We only prove the completeness of \(\mathrm H\) here. Let \(\mathbf{y}^k\) be a Cauchy sequence in \(\mathrm H\) so that \(y_1^k, y_2^k, (y_2^k)_x, y_3^k, y_4^k, y_5^k\) are Cauchy sequences in \({\mathcal {L}}_{2} (0,L)\) with \(y_2^k(0)=y_2^k(L)=0.\) By the completeness of \({\mathcal {L}}_{2} (0,L)\) and \(H^1_0(0,L)\) there exists \(z_1, z_2, z_{22}, z_3, z_4, z_5\in {\mathcal {L}}_{2} (0,L)\) such that \(y_1^k\rightarrow z_1,\) \(y_2^k\rightarrow z_2,~(y_2^k)_x\rightarrow z_{22}, ~y_3^k\rightarrow z_3,~ y_4^k\rightarrow z_4,~ y_5^k\rightarrow z_5\) as \(k\rightarrow \infty .\) Moreover, \(z_2(0)=z_2(L)=0.\) Therefore, for \(\varphi \in C^\infty _0(0,L),\)
Therefore, \(z_2\in H^1(0,L)\) with \((z_2)_x=z_{22}\) as distributions, and by construction, \(z_{22}\in {\mathcal {L}}_{2} (0,L).\) Therefore, \(z_2, (z_2)_x \in {\mathcal {L}}_{2} (0,L).\) Similarly, for all \(\varphi \in C^\infty _0(0,L),\)
as \(k\rightarrow \infty .\) Hence, \(\xi (z_2)_x -z_3+\frac{\gamma }{{\varepsilon _{3}}}z_4=0\) as distributions. \(\square \)
Define the operator
with
and the B and \(B^*\) operators with the new state are
and \(\mathcal B \mathcal B^*\in \mathcal L( \mathrm H, \mathrm H).\)
Lemma 1
The operator \(\mathcal A : \mathrm{Dom} ({\mathcal A}) \rightarrow \mathrm H.\)
Proof
Let \(\mathbf{y}\in \mathrm{Dom}(\mathcal A).\) Then
For given \(\mathbf{y}\in \mathrm{Dom}(\mathcal A),\) it is obvious that \(y_2-(y_3)_x\in {\mathcal {L}}_{2} (0,L), -\frac{\mu }{\xi {\varepsilon _{3}}}y_1\in H^1_0(0,L), ~ -\frac{\mu }{{\varepsilon _{3}}}(y_1)_x+\frac{\gamma }{{\varepsilon _{3}}}(y_5)_x,~ (y_5)_x,~ (y_3)_x+\frac{\alpha }{\rho }(y_4)_x \in {\mathcal {L}}_{2} (0,L).\) Next we check if \(\mathcal A \mathbf{y}\) satisfies the compatibility condition (22) for \(\mathcal A \mathbf{y}:\)
This completes the proof of \(\mathcal A \mathbf{y}\in \mathrm H.\) \(\square \)
Theorem 2
For any \(\mathbf{g}\in \mathrm {H}\) there is \(\mathbf{y}\in \mathrm{Dom}(\mathcal {A} )\) so that \(\mathcal {A} \mathbf{y}=\mathbf{g}. \) That is, \(0\in \rho (\mathcal A).\)
Proof
This is equivalent to solving the system of equations for \( \mathbf{y}\in \text {Dom}(\mathcal {{\tilde{A}}} ).\) Simplifying the equations leads to
Since \(\mathbf{g}\in \mathrm H,\) the gauge condition \(\xi (g_2)_x -g_3+\frac{\gamma }{{\varepsilon _{3}}}g_4=0\) is satisfied. By using the Green’s function theory, we obtain \( \mathbf{y}\in \text {Dom}(\mathcal {{\tilde{A}}} )\) where
where \(\mathbf{y}\) clearly satisfies the compatibility condition \(\xi (y_2)_x -y_3+\frac{\gamma }{{\varepsilon _{3}}}y_4=0.\) \(\square \)
Theorem 3
The operator \(\mathcal A\) satisfies \(\mathcal {A}^*=-\mathcal {A}\) on \(\mathrm {H},\) and
\({\mathcal {A}}: \mathrm{Dom} (\mathcal A) \subset \mathrm H \rightarrow \mathrm H\) defined by (26) is the generator of a unitary semigroup \(\{e^{{\mathcal {A}}t}\}_{t\ge 0}.\)
Proof
Let \(\mathbf{y}, \mathbf{z}\in \mathrm{Dom}({\mathcal {A}}).\) Then, by integration by parts and rearranging the terms, we have
with \(\mathrm{Dom}(\mathcal A^*)=\mathrm{Dom}(\mathcal A).\) By Theorem 2, is the generator of a unitary semigroup.
\(\square \)
The system (11) can be written as
where \(\mathbf{y}_0\) is defined by (18) and (21), the control operator B is defined by
where I is the identity operator.
Theorem 4
Let \(T>0,\) and \(i_s(t)\in L^2(0,T).\) For any \(\mathbf{y}_0 \in \mathrm { H},\) \(\mathbf{y}\in C[[0,T]; \mathrm H],\) and there exists a positive constants c(T) such that (29) satisfies
Proof
The operator \(\mathcal A: {\mathrm{Dom}}(\mathcal A) \rightarrow \mathcal H\) is a unitary semigroup by Lemma 3. Therefore, it is an infinitesimal generator of \(C_0-\)semigroup of contractions by Lümer-Phillips theorem [38]. Since with \(\frac{di_s}{dx}=0\) by the surface continuity condition (6), B is admissible for \(\{e^{\mathcal A t}\}_{t\ge 0}:\)
Hence the conclusion (30) follows. \(\square \)
3.1 Lack of Stabilization of the Lorenz Gauge Model (11)
Along the trajectories of (11), the energy (13) satisfies
The adjoint operator \(B^*\) of B is found by
where \(B^*\mathbf{y}= \frac{1}{h}\int _0^L {\bar{y}}_2(x) dx.\) We investigate the asymptotic stability of the closed-loop system with the \(B^*-\)feedback.
This leads to the feedback control
where \(K_1 >0\) is arbitrary. This can be regarded as induced total charge (voltage) accumulated at the electrodes of the beam.
The uncontrolled system (11) has an infinite number of eigenvalues on the imaginary axis. Moreover, the feedback control operator is compact and rank 1. Thus, it is not possible to exponentially stabilize the piezo-electric beam; see for instance [6].
The abstract state-space description (29) with (32) becomes
The operator \(\mathcal {{\tilde{A}}}: \mathrm{Dom}(\mathcal {{\tilde{A}}})\subset \mathrm H \rightarrow \mathrm H \) defined by (33) with domain \(\mathrm{Dom}(\mathcal {{\tilde{A}}})=\mathrm{Dom}(\mathcal { A})\) is densely defined in \( \mathrm H.\)
Lemma 2
The infinitesimal generator \(\mathcal {{\tilde{A}}}\) satisfies \(\mathcal {{\tilde{A}}}^*=-\mathcal {{\tilde{A}}}(-K_1)\) on \(\mathrm {H}.\) Moreover it is dissipative and it satisfies \(\mathrm{Re}\left\langle \mathcal {{\tilde{A}}} \mathbf{y}, \mathbf{y}\right\rangle _{{\mathrm H}}\le 0.\)
Theorem 5
\(\mathcal {{\tilde{A}}}: \mathrm{Dom}( {\mathcal A}) \rightarrow {\mathrm H}\) is the infinitesimal generator of a \(C_0-\)semigroup of contractions. Therefore for every \(T\ge 0,\) and \(\mathbf{y}^{\mathbf{0}}\in \mathrm{Dom}( {\mathcal A})\) solves (33) and \( \mathbf{y}\in C([0,T]; \mathrm{Dom}( {\mathcal A}))\cap C^1\left( [0,T]; {\mathrm H}\right) .\) Moreover, the spectrum \(\sigma ({\mathcal {{\tilde{A}}}})\) of \(\mathcal {{\tilde{A}}}\) has all isolated eigenvalues.
Proof
(outline): The domain \(\mathrm{Dom}({\mathcal {{\tilde{A}}}})\) is densely defined and is compact in \({\mathrm H}\) due to the compact embedding of \(H^1(0,L)\subseteq L^2(0,L),\) and \(H^2(0,L)\subseteq H^1(0,L).\) Moreover, \(0\in \rho ({ \mathcal A} )\) by Theorem 2. Therefore, \((\lambda I -{\mathcal {{\tilde{A}}}})^{-1}\) is compact at \(\lambda =0, \) (and thus compact for all \(\lambda \in \rho (\tilde{\mathcal A}) \)). Hence the spectrum of \({\mathcal {{\tilde{A}}}}\) has all isolated eigenvalues. \(\square \)
Theorem 6
Let \(a_1=\frac{ 2n\pi }{L}, a_2=\frac{ 2m \pi }{L},\) and \(m^2+ n^2>\frac{\mu \rho L^2 }{16 \alpha \xi {\varepsilon _{3}} \pi ^2}\) for \(m,n\in {\mathbb {N}}\) where \(a_1\) and \(a_2\) are defined by (41). For \(\mathbf{y}^{\mathbf{0}}\in \mathrm H ,\) the semigroup \(\{e^{\mathcal {{\tilde{A}}}t}\}_{t\ge 0}\) corresponding to (33) is not asymptotically stable in \(\mathrm {H}\), i.e. \(\Vert e^{\mathcal {{\tilde{A}}}t}\mathbf{y}^{\mathbf{0}}\Vert _{\mathrm {H}}\nrightarrow 0, ~~~ t\rightarrow \infty .\)
Proof
Since \(0\in \rho (\mathcal {{\tilde{A}}}),\) we show that there are eigenvalues on the imaginary axis, or in other words, the set
has nontrivial solutions \(\mathbf{y}\ne 0.\) It will then follow that the system (33) is NOT strongly stable.
Let \(\lambda =\imath \tau .\) The eigenvalue problem is
together with the boundary conditions
and the compatibility condition (22). By using (22), we can eliminate the variables \(y_1, y_3, y_4,\) and the system (35) can be reduced to
Now we let \(Z=[y_2~ (y_2)_x~ y_5 ~(y_5)_x]^{\mathrm{T}}.\) We write the system (37)
with the boundary conditions
The solution to (38) is \(Z=e^{\mathcal D x} C\) where \(C=[c_1 ~ c_2 ~ c_3 ~ c_4]^{\mathrm{T}}\) is a vector of arbitrary coefficients. The characteristic equation for \(\mathcal D\) is
Let
and
Note that in the description of \(B^2-4AC,\) the term \({\varepsilon _{3}}\left( \alpha -\frac{\gamma ^2}{{\varepsilon _{3}}}\right) >0\) since \(\frac{\gamma ^2}{{\varepsilon _{3}}}\) is due to the quadratic-through thickness electrical field assumption. Indeed, the piezoelectric beam becomes relatively less-stiff even though \(\frac{\gamma ^2}{{\varepsilon _{3}}} \ll \alpha .\) We consider the case
and therefore, \(B^2-4AC>0.\) This generates four complex conjugate roots of the characteristic equation:
The solution of (38) is written as \(Z=P e^{Jx}P^{-1} C\) where
the constants \(b_1,b_2\) are defined as in
where \(b_1-b_2\ne 0.\) By using the boundary conditions \(z_1(0)=z_3(0)=0\) in (39), we get \(c_1=c_3=0.\) Therefore, the solution to (38) is
Choosing \(a_1\) and \(a_2\) as in the statement of the theorem, we immediately get \(c_2=\mathrm{free}:\)
For simplicity take \(c_2=1.\) Therefore, a family of nontrivial solutions of the overdetermined eigenvalue problem (together with the feedback condition \(\int _0^L z_1 (x)dx=0\)) is
and therefore, the solution of the system (35) is
Note that for the choices of \(a_1=\frac{2n \pi }{L}\) and \(a_2=\frac{ 2m \pi }{L},\) \(\tau _{nm}\) can be determined as the following
The condition (42) is equivalent to
Let the initial condition of the control-free problem (29) be
where \(y_{i,nm}\)’s are defined by (47). Then the solution of the initial value problem (29) is
where n, m can be chosen to satisfy (48). This type of solution results in \(B^*\mathbf{y}_{nm} \equiv 0.\)
\(\square \)
Remark 1
Note that the wave speeds of mechanical and electro-magnetic equations in (11) are \(\left\{ \sqrt{\frac{\alpha }{\rho }},\sqrt{\frac{\mu }{{\varepsilon _{3}}}}\right\} ,\) respectively, where \(\sqrt{\frac{\mu }{{\varepsilon _{3}}}}\sim 10^8\) is the speed of light with the material parameters in (70), and \(\frac{\mu }{{\varepsilon _{3}}}>> \frac{\alpha }{\rho }.\) Since \(\xi \) in (70) is around \(10^{-6}\) the condition \(\tau>\frac{\mu }{\xi {\varepsilon _{3}}}>>\frac{\mu }{{\varepsilon _{3}}}\) implies the high-frequency eigenvalues \(\lambda =\imath \tau \) in (35). Thus, Theorem (5) proves the lack of stability for the high-frequency vibration modes. As \(\tau <\frac{\mu }{\xi {\varepsilon _{3}}},\) the form of the eigenvalues don’t look as good. For that reason, the lack of stability for these modes are shown through simulations in Sect. 6.
3.2 Additional Remedial Controller
In addition to the current controller \(i_s(t)\), an additional mechanical controller g(t), as in [22], can be considered at the tip \(x=L\) so that the boundary condition \({\alpha } v_{x}(L) + \gamma \phi (L) +\gamma {\dot{\eta }}(L)=0\) in (11) is replaced by
Then the system (29) is replaced by
where \( B= \left[ \begin{array}{cc} 0 &{} 0 \\ \frac{1}{ \xi {\varepsilon _{3}} h} &{} 0 \\ 0 &{} 0 \\ 0 &{} 0 \\ 0 &{} \frac{1}{h}\delta (x-L) \end{array} \right] ,\) with \(B^*\mathbf{y}= \left[ \begin{array}{ccccc} 0 &{} \frac{1}{ \xi {\varepsilon _{3}} h} &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0 &{} \frac{1}{h}\delta (x-L) \end{array} \right] \mathbf{y}= \left[ \begin{array}{c} \frac{1}{\xi {\varepsilon _{3}} h}\int _0^L y_2(x) dx\\ \frac{1}{h}y_5(L) \end{array} \right] ,\) and \(\delta (x-L)\) is the Dirac-Delta distribution at \(x=L.\) Then, the energy of the closed-loop system
where \(K= \left[ \begin{array}{ccccc} 0 &{} 0 &{}0&{}0 &{}0\\ -\xi ^2 {\varepsilon _{3}}^3 h^2 K_1&{} 0 &{} 0 &{} 0 &{}0 \\ 0 &{} 0 &{} 0&{} 0 &{}0\\ 0 &{} 0 &{} 0&{} 0 &{}0\\ 0 &{} 0&{} 0&{} 0&{} -h^2 K_2\end{array} \right] \) with \(K_1,K_2>0\) satisfies \(\frac{d{{\mathrm {E}}}(t)}{dt}=-K_1 \left( \int _0^L \psi _2(x) dx\right) ^2 -K_2 |y_5(L)|^2.\) By using the similar methodology in Sect. 3.1, it can be shown that there are no eigenvalues on the imaginary axis, or in other words, the set
has only the trivial solution \(\mathbf{y}=0.\) Hence, the following theorem is immediate:
Theorem 7
For \(a_1\ne \frac{n \pi }{L}\) and \(a_2\ne \frac{ m \pi }{L}\) for \(m,n\in {\mathbb {N}},\) and \(\mathbf{y}^{\mathbf{0}}\in \mathrm H ,\) the semigroup \(\{e^{(\mathcal A + K B B^*)t}\}_{t\ge 0}\) corresponding to (33) is asymptotically stable in \(\mathrm {H}\), i.e. \(\Vert e^{(\mathcal A + K B B^*)t}\mathbf{y}^{\mathbf{0}}\Vert _{\mathrm {H}}\rightarrow 0, ~~~ t\rightarrow \infty .\)
Remark 2
Following up with Remark 1, the additional remedial controller help asymptotically stabilize both mechanical and electro-magnetic states. This includes both high and low-frequency vibration modes. See Sect. 6 for the power of the remedial controller.
Remark 3
An additional feedback controller can be considered in the distributed (stronger) sense, i.e., viscous damping [40] or fractional damping [16]. This type of feedback controller makes the voltage-controlled beam equations exponentially stable. It is an open problem whether the exponential stability is retained for the current-controlled beam equations (33).
4 Well-Posedness of the Coulomb Gauge Model (12)
Note that a similar analysis of this section is carried out for (12) with free-free boundary conditions [27, 37] and different choice of states. Our paper deals with clamped-free boundary conditions. For \(t\in {\mathbb {R}}, \) the natural energy associated with (12) is
Define the linear space
with the bilinear form
Since \(P_\xi \) is a nonnegative operator, the following result is straightforward as in Theorem 1:
Lemma 3
The form (24) defines an inner product on the linear space \(\mathrm H\). Moreover, E is the norm induced by this inner product and \(\mathrm H\) is complete.
Let
Then (12) can be written as the following
where \( A= \left( {\begin{array}{*{20}c} 0 &{} I &{} -D_x &{} 0 &{} 0 \\ \frac{-\mu }{\xi {\varepsilon _{3}}} I &{} 0 &{} 0 &{} 0 &{} \frac{-\gamma }{{\varepsilon _{3}}}D_xP_\xi D_x \\ \frac{-\mu }{{\varepsilon _{3}}}D_x &{} 0 &{} 0 &{} 0&{} \frac{\gamma }{{\varepsilon _{3}}} (D_x-D_x P_\xi ) \\ 0&{} 0 &{} 0 &{} 0 &{} D_x \\ 0&{} 0 &{} \frac{\gamma }{\rho } D_x &{} \frac{\alpha }{\rho } D_x +\frac{\gamma ^2}{\rho {\varepsilon _{3}}} P_\xi D_x &{} 0 \\ \end{array}} \right) \quad \mathrm{and}\)
and the B and \(B^*\) operators with the new state are
and \(\mathcal B \mathcal B^*\in \mathcal L( \mathrm H, \mathrm H).\) We recall the following result from [27]:
Lemma 4
Let \(\mathrm{Dom} (D_x^2)=\{w\in H^2(0,L) ~:~ w_x(0)=w_x(L)=0\}.\) The operator \(\frac{1}{\xi }(P_\xi -I)\) is continuous, self-adjoint and non-positive on \({\mathcal {L}}_{2} (0,L).\) Moreover, for all \(w\in \mathrm{Dom} (P_\xi ), \) \( J=D_x^2 ~P_\xi = D_x^2 (I-\xi D_x^2)^{-1}w.\)
Lemma 5
The operator A maps \(\mathrm{Dom}(A)\subset \mathrm H\) to H, and is densely defined in \( \mathrm H.\)
Proof
Let \(\mathbf{y}\in \mathrm{Dom}(A).\) Then
Since \(\mathbf{y}\in W_1,\) \(y_1\in H^1(0,L)\) and \(y_4\in H^1_0(0,L),\) \((y_1)_x\in {\mathcal {L}}_{2} (0,L),\) and by the definition of the operator in (10) \(P_\xi (y_1)_x\in H^2(0,L).\) Therefore \(({P_\xi } (y_1)_x)_x + y_4\in H^1_0(0,L).\) Hence, \(A \mathbf{y}\in W_1.\) \(\square \)
Let’s also verify that the gauge condition (52) for \(A\mathbf{y}\) is satisfied:
where we use Lemma 4.
The following result is immediate, and it can obtained in a similar fashion as in Sect. 3:
Theorem 8
The operator \({A}: \mathrm{Dom} (A)\subset \mathrm H \rightarrow \mathrm H\) satisfies \({A}^*=-{A}\) on \(\mathrm H,\) and A
defined by (26) is the generator of a unitary semigroup \(\{e^{{A}t}\}_{t\ge 0}\) on \(\mathrm H.\) Letting \(T>0\) and \(i_s(t)\in L^2(0,T),\) for any \(\mathbf{y}_0 \in \mathrm { H},\) \(\mathbf{y}\in C[[0,T]; \mathrm H],\) and there exists a positive constants c(T) such that (55) satisfies
4.1 Lack of Stabilization of the Coulomb Gauge Model (12)
Similar to the Lorenz gauge case, along the trajectories of (12), the energy (51) satisfies \(\frac{d{{\mathrm {E}}}(t)}{dt}=i_s(t) \left( \int _0^L \psi _2(x) dx\right) .\)
We investigate the asymptotic stability of (12) for the same \(B^*\)-feedback (31). This leads to the feedback control \(i_s^1(t)=-K_1 h^2 \xi ^2{\varepsilon _{3}}^2 \int _0^L ~ {\dot{\theta }}(z)dz \) where \(K_1 >0\) is arbitrary.
The abstract state-space description (29) becomes
The operator \({{\tilde{A}}}: \mathrm{Dom}( {{\tilde{A}}})\subset \mathrm H \rightarrow \mathrm H \) defined by (57) with domain \(\mathrm{Dom}({{\tilde{A}}})=\mathrm{Dom}({ A})\) is densely defined in \( \mathrm H.\)
Lemma 6
The infinitesimal generator \( {{\tilde{A}}}\) satisfies \( {{\tilde{A}}}^*=- {{\tilde{A}}}(-K_1)\) on \(\tilde{\mathrm H}.\) Moreover it is dissipative and it satisfies \({\mathrm{Re}}\left\langle {{\tilde{A}}} \mathbf{y}, \mathbf{y}\right\rangle _{\tilde{\mathrm {H}}}\le 0.\)
Theorem 9
\(\tilde{{\mathcal {A}}}: \mathrm{Dom}( {{\tilde{A}}}) \rightarrow \tilde{\mathrm H}\) is the infinitesimal generator of a \(C_0-\)semigroup of contractions. Therefore for every \(T\ge 0,\) and \(\mathbf{y}^{\mathbf{0}}\in \mathrm{Dom}( {{\tilde{A}}})\) we have \( \mathbf{y}\in C\left( [0,T]; \mathrm{Dom}({{\tilde{A}}})\right) \cap C^1\left( [0,T]; {\mathrm H}\right) .\) Moreover, the spectrum \(\sigma (\tilde{{A}})\) of \(\tilde{{\mathcal {A}}}\) has all isolated eigenvalues.
Proof
By Lemma 6, \(\tilde{{A}}: \mathrm{Dom}( {{\tilde{A}}}) \rightarrow {\mathrm H} \) is the infinitesimal generator of a \(C_0-\)semigroup of contraction by Lümer-Phillips theorem [38]. Therefore for every \(T\ge 0,\) and \(\mathbf{y}^{\mathbf{0}}\in \mathrm{Dom}(\tilde{{A}})\) solves (55), and we have \( \mathbf{y}\in C\left( [0,T]; \mathrm{Dom}(\tilde{{A}})\right) \cap C^1\left( [0,T]; \tilde{{\mathrm H}}\right) .\)
Finally, since \(\mathrm{Dom}(\tilde{ A})\) is densely defined and compact in \({\mathrm H},\) \(0\in \rho ({\tilde{A}}), \) we have \((\lambda I -{\tilde{ A}})^{-1}\) is compact at \(\lambda =0,\) thus compact for all \(\lambda \in \rho ({\tilde{ A}}).\) Hence the spectrum of \({\tilde{A}}\) has all isolated eigenvalues. \(\square \)
We have the following result that lists some choices of material parameters causing instabilities.
Theorem 10
Let \(a_1=\frac{2n\pi }{L}, a_2=\frac{2m\pi }{L}\) in (41), and \(m^2+ n^2>\frac{\mu \rho L^2 }{16 \alpha \xi {\varepsilon _{3}} \pi ^2}\) for \(m,n\in {\mathbb {N}}\) where \(a_1\) and \(a_2\) are defined by (41). For \(\mathbf{y}^{\mathbf{0}}\in { {\mathrm H}},\) the semigroup \(\{e^{ {{\tilde{A}}}t}\}_{t\ge 0}\) corresponding to (33) is not asymptotically stable in \({\mathrm {H}}\), i.e. \(\Vert e^{{{\tilde{A}}}t}\mathbf{y}^{\mathbf{0}}\Vert _{\tilde{\mathrm {H}}}\nrightarrow 0, ~~~ t\rightarrow \infty .\)
Proof
Since \(0\in \rho ( {{\tilde{A}}}),\) we show that there are eigenvalues on the imaginary axis, or in other words, the set
has nontrivial solutions. It will then follow that the system (33) is not strongly stable. This is equivalent to saying that the current-controlled problem is not approximately controllable on \(\mathrm H.\) \(\square \)
Let \(\lambda =\imath \tau .\) The the eigenvalue problem \({\tilde{A}} \mathbf{y}=\lambda \mathbf{y}=\imath \tau \mathbf{y}\) is
with the boundary conditions
As well, (58) implies that \(\int _0^L y_2 (x) dx=0.\) Let \(r:=P_\xi (y_5)_x.\) By (10), \(r-\xi r_{xx}=(y_5)_x. \) By using the gauge condition, i.e. \(y_3=\xi (y_2)_x\), (59) reduces to
with the boundary conditions
and \(\int _0^L y_2(x)dx=0.\)
Solving (60) for \(r_{xx}, (y_1)_{xx}\) and \((y_3)_{xx}\) yields
Now we let \(Z=[y_5~ (y_5)_x~ y_2 ~(y_2)_x~ r ~r_x]^{\mathrm{T}}.\) We write the system (61)
The solution to (62) is \(Z=e^{{\tilde{\mathcal D}} x} C\) where \(C=[c_1 ~ c_2 ~ c_3 ~ c_4~c_5 ~c_6]^{\mathrm{T}}\) is a vector of arbitrary coefficients. The characteristic equation is
For the same A, B, C in (40), and real \(\tau \), define
For
the first four eigenvalues turn to be similar to those obtained in Sect. 3.1:
where \(a=\frac{1}{\sqrt{\xi }},\) and \(a_1,a_2\) are defined in (41). Since \(a_3 \ne a_4\), in (63), \(b_3-b_4\ne 0, b_{33}-b_{44}\ne 0,\) the solution of (62) is written as \(Z=P e^{Jx}P^{-1} K\) where \(P,e^{Jx}\) are defined by
where \(b=\frac{-\imath \tau {\varepsilon _{3}}a}{\gamma }.\) Using the boundary conditions at \(x=0\) yields \(c_1=c_3=c_6=0.\) By using the other three boundary conditions at \(x=L,\) the feedback condition \(\int _0^L z_3 (x)dx=0,\) and \(a_1 =\frac{2n\pi }{L}\) and \(a_2=\frac{2m\pi }{L}\), we obtain an overdetermined system of equations for the unknowns \(\{c_2,c_4,c_5\}.\) After row reducing it, we obtain that \(c_5\) is a free variable. The solution of the eigenvalue problem
where \(C(a_k)\ne 0,\)
This implies that the eigenvalue problem has nontrivial solutions. It is the same conclusion in Sect. 3.1 that the integral-type observation \(\int _0^L z_3(x)dx\) does not help to show approximate controllability. Since \(c_3\) can be chosen to be a free variable, the solution of the eigenvalue problem can be chosen.
For \(a_1\ne \frac{n \pi }{L}\) and \(a_2\ne \frac{ m \pi }{L},\) the remedial controller (50) discussed in Sect. 3.2 can be also considered here to obtain exactly the same asymptotic stability result, see Theorem 7.
5 Electrostatic/Guasi-Static Model
To see the difference between the electrostatic/quasi-staic and fully magnetic models, the charge or voltage-controlled model [27], obtained by a thorough variational approach, is chosen to start with. The magnetic effects are totally freed. The mechanical equation is coupled to the dynamic charge equation:
Here \(\sigma _s(t)\) is the applied charge to the electrodes on the top and bottom surfaces of piezoelectric beam. The model (66) is a well-studied wave-equation with a dynamic boundary controller. This type of model is first proposed and analyzed by [43].
The natural energy associated with (66) is
Define the linear space \(H^1_L(0,L)\times L^2(0,L)\times {\mathbb {C}}\) and the bilinear form
Let \(\mathbf{y}=[y,\dot{y}, \sigma _s]^{\mathrm{T}}\) be the states, \(\mathcal A\mathbf{y}= [y_2, \frac{\alpha +\frac{\gamma ^2}{{\varepsilon _{3}}}}{\rho }(y_1)_{xx}, 0],\) and \(B=[0 ~0 ~ 1]\) where
Then (66) can be written as
The model (69) is well studied in the literature. This model is never exponentially stabilizable [57]. However, by the designing the feedback controller \(i_s(t)=\frac{\gamma }{{\varepsilon _{3}} h}\left( \dot{v}(L,t)-K_3\sigma _s(t)\right) ,~~K_3>0,\) we have the following optimal polynomial decay result
Theorem 11
(Thm. 2.1, [57]) Let \(T>0,\) and \(i_s(t)=\left( \dot{v}(L,t)-K_3\sigma _s(t)\right) ,\) \(K_3\in {\mathbb {R}}^+.\) For any \(\mathbf{y}_0 \in \mathrm{Dom}(\mathcal A),\) there exists a constant \( C(\mathbf{y}_0,K_3)>0\) such that (67) satisfies
Note that the feedback \(i_s(t)=-K_3B^*\mathbf{y}=-K_3\sigma _s\) does not make the energy dissipative at all.
6 Simulations
Now, we illustrate a numerical experiment by a sample example for the material parameters
Our goal is to simulate the low-frequency vibrational states. First, to carry out numerical calculations, we non-dimensionalize the systems of Eqs. (11)–(12) as in Table 1. We also non-dimensionalize the model (66) in Table 2.
Consider the discretization of the interval [0, L] as the following
The following are the second order finite differences where \(i=,0,1,\ldots , N.\) Henceforth, to simplify the notation, we use \(z(x_{i})=z_i.\) The approximations for different order derivatives are
We consider simulations for \(N=40,\) and \(T_{\mathrm{final}}=300.\) We adopt the methodology to derive a filtered semi-discrete scheme in finite differences as it is applied successfully for a boundary-controlled wave equation [49]. The high frequency spurious solutions, causing artificial instability (unobservability), are filtered by adding a viscosity term to the wave equation. After this type of filtering, the discrete control-free problem becomes uniformly exactly observable with respect to the mesh parameter, and therefore, the controlled problem mimics the same exponential stability behavior of the original PDE problem. By adopting the same strategy, additional viscosity terms
are added to the \(\{v,\phi , \theta ,\eta \}-\)equations with appropriate coefficients in Tables 1 and 2, respectively.
The simulations for the models in Tables 1 and 2are performed for four different control strategies as shown in Table 3. The effect of numerical viscosity terms can be seen in Tables 4, 5, 6, 7, 8 and 9. The effect of these terms diminishes as the mesh size N increases. These effects are getting bigger in the Coulomb and Lorenz gauge models since there are viscosity terms as many as the equations.
One general observation in Tables 4, 5, 6, 7 and 8 that the integral controller \(i_s^1(t)=-K_1\int _0^L y_2(x,t) dx,~ K_1>0\) itself is not even able to stabilize any state with low-frequency vibration modes for the Coulomb and Lorenz gauge models. The remedial mechanical feedback controller \(g(t)=-K_2 \dot{v}(L,t), ~K_2>0\) acting on the boundary in each case is able to stabilize the vibrations very fast by itself. In some cases, the current controller has a minor effect to improve the controllers’ performance for the electro-magnetic states \(y_1, y_2, y_3.\) The polynomial stability of the electrostatic/quasi-static model reflects in the simulations, see Tables 7 and 8.
The mechanical states are stabilized faster for the electrostatic/quasi-static equations regardless of optimal choice of the feedback gain \(K_3>0\) in Table 2.
7 Conclusion
In this paper, we consider two well-posed models for current-controlled piezoelectric beams including full electro-magnetic effects; one uses the Coulomb gauge and the other one uses Lorenz gauge. Both models are shown to be not asymptotically stable with only one feedback controller if the material parameters satisfy certain conditions. These conditions correspond to the high-frequency vibration modes. This is due to the addition of the magnetic effects back into to the model; these effects are minor in comparison to the mechanical and electrical effects. The lack of stability for the low-frequency modes is almost impossible to show. Therefore, low frequency vibrational modes are chosen as initial conditions for simulations. It is observed that the choice of a gauge condition does not play an important role in terms of the stabilization, which is good since one can go with either model depending on the physical assumptions for the electromagnetic equations (hyperbolic or elliptic).
A parallel discussion is carried out in [52,53,54] where both quasi-static and fully dynamic cases are handled for voltage or mechanical actuation in the port-Hamiltonian framework. It is shown that asymptotic stability is at stake for the quasi-static model whereas the asymptotic stability can be achievable by retaining the magnetic effects in the model and with additional feedback controllers. Relevant problems have been studied in [31, 32, 36] in details where the quasi-static model with enough number of controllers provide an asymptotic stability result.
References
Adriaens, H., de Koning, W.L., Banning, R.: Modeling piezoelectric actuators. IEEE/ASME Trans. Mechtron. 5–4, 331–341 (2000)
Banks, H.T., Smith, R.C., Wang, Y.: Smart Material Structures: Modelling Estimation and Control. Mason, Paris (1996)
Bazghaleh, M., Grainger, S., Mohammadzaheri, M.: A review of charge methods for driving piezoelectric actuators. J. Intell. Mater. Syst. 15, 29–110 (2017)
Cao, Y., Chen, X.B.: A survey of modeling and control issues for piezo-electric actuators. J. Dyn. Syst. Meas. Control 137–1, 014001 (2014)
Comstock, R.: Charge control of piezoelectric actuators to reduce hysteresis effects, (United States Patent # 4,263,527, Assignee: The Charles Stark Draper Labrotary)
Curtain, R.F., Zwart, H.J.: An Introduction to Infinite-Dimensional Linear Systems Theory. Springer, Berlin (1995)
Dagdeviren, C., et al.: Conformal piezoelectric energy harvesting and storage from motions of the heart, lung, and diaphragm. Proc. Natl. Acad. Sci. US A 111, 1927–1932 (2014)
Dagdeviren, C., et al.: Recent progress in flexible and stretchable piezoelectric devices for mechanical energy harvesting. Sens. Actuation Extreme Mech. Lett. 9(1), 269–281 (2016)
Devasia, S., Eleftheriou, E., Reza Moheimani, S.O.: A Survey of control issues in nanopositioning. IEEE Trans. Control Syst. Technol. 15, 802–823 (2007)
EXACTO Guided Bullet Demonstrates Repeatable Performance against Moving Targets, DARPA, 27 April 2015. Last accessed: (5 January 2018)
Dautray, R., Lions, J.L.: Spectral Theory and Applications. Mathematical Analysis and Numerical Methods for Science and Technology, vol. 1. Springer, Berlin (1988)
Dautray, R., Lions, J.L.: Spectral Theory and Applications. Mathematical Analysis and Numerical Methods for Science and Technology, vol. 1. Springer, Berlin (1990)
Destuynder, Ph, Legrain, I., Castel, L., Richard, N.: Theoretical, numerical and experimental discussion of the use of piezoelectric devices for control-structure interaction. Eur J. Mech. A Solids 11, 181–213 (1992)
Erturk, A., Inman, D.: A distrbiuted parameter model for cantilever model for piezoelectric energy harvesting from base excitations. J. Vib. Acoust. 130, 041002 (2008)
Fleming, A.J., Moheimani, S.O.R.: Precision current and charge amplifiers for driving highly capacitive piezoelectric loads. Electron. Lett. 39–3, 282–284 (2003)
Freitas, M.M., Ramos, A.J.A., Almeida Jr., D.S., Özer, A.Ö.: Long time dynamics for fractional piezoelectric system with magnetic effects and Fourier Law, submitted
Giurgiutiu, V.: Structural Health Monitoring/with Piezoelectric Wafer Active Sensors. Elsevier Academic Press, Amsterdam (2008)
Gu, G.Y., Zhu, L.M., Su, C.Y., Ding, H., Fatikow, S.: Modeling and control of piezo-actuated nanopositioning stages: a survey. IEEE Trans. Autom. Sci. Eng. 13–1, 313–332 (2016)
Hagood, N.W., Chung, W.H., Von Flotow, A.: Modelling of piezoelectric actuator dynamics for active structural control. J. Intell. Mater. Syst. Struct. 1, 327–354 (1990)
Hansen, S. W.: Analysis of a plate with a localized piezoelectric patch, Conference on Decision & Control, Tampa, Florida, 2952–2957 (1998)
Jiang, S.N., Jiang, Q., Li, X.F., Guo, S.H., Zhou, H.G., Yang, J.S.: Piezoelectro-magnetic waves in a ceramic plate between two ceramic half-spaces. Int. J. Solids Struct. 43, 5799–5810 (2006)
Lagnese, J.E., Lions, J.-L.: Modeling Analysis and Control of Thin Plates. Masson, Paris (1988)
Lee, P.C.Y.: A variational principle for the equations of piezoelectromagnetism in elastic dielectric crystals. J. Appl. Phys. 69–11, 7470–7473 (1991)
Main, J.A., Garcia, E.: Design impact of piezoelectric actuator nonlinearities. J. Guid. Control Dyn. 20–2, 327–332 (1997)
Main, J.A., Garcia, E., Newton, D.V.: Precision position control of piezoelectric actuators using charge feedback. J. Guid. Control Dyn. 18–5, 1068–73 (1995)
Newcomb, C., Flinn, I.: Improving the linearity of piezoelectric ceramic actuators. Electron. Lett. 18–11, 442–443 (1982)
Morris, K.A., Özer, A.Ö.: Modeling and stabilizability of voltage-actuated piezoelectric beams with magnetic effects. SIAM J. Control Optim. 52–4, 2371–2398 (2014)
Morris, K.A., Özer, A.Ö.: Comparison of stabilization of current-actuated and voltage-actuated piezoelectric beams, the \(53^{{\rm rd}}\) Proceedings of the IEEE Conference on Decision & Control, Los Angeles, CA, 571–576 (2014)
Oh, S.J., Tataru, D.: Local well-posedness of the \((4 + 1)-\)dimensional Maxwell–Klein–Gordon equation at energy regularity. Ann. PDE 2, 2 (2016)
Özer, A.Ö.: Further stabilization and exact observability results for voltage-actuated piezoelectric beams with magnetic effects. Math. Control Signals Syst. 27, 219–244 (2015)
Özer, A.Ö.: Modeling and control results for an active constrained layered (ACL) beam actuated by two voltage sources with/without magnetic effects. IEEE Trans. Autom. Control 62–12, 6445–6450 (2017)
Özer, A.Ö.: Potential formulation for charge or current-controlled piezoelectric smart composites and stabilization results: Electrostatic vs. Quasi-static vs. Fully-dynamic approaches. IEEE Trans. Autom. Control 64(3), 989–1002 (2018)
Özer, A.Ö.: Exponential stabilization of the smart piezoelectric composite beam with only one boundary controller. Proc. IFAC Conf. 51–3, 80–85 (2018)
Özer, A.Ö., Hansen, S.W.: Uniform stabilization of a multi-layer Rao–Nakra sandwich beam. Evol. Equ. Control Theory 2–4, 195–210 (2013)
Özer, A.Ö., Hansen, S.W.: Exact boundary controllability results for a multilayer Rao–Nakra sandwich beam. SIAM J. Control Optim. 52–2, 1314–1337201 (2014)
Özer, A.Ö., Khenner, M.: An alternate numerical treatment for the nonlinear PDE models of piezoelectric laminates, SPIE Proceeding of Active and Passive Smart Structures and Integrated Systems XIII, 109671R (2019)
Özer, A.Ö., Morris, K.A.: Modeling and stabilization of current-controlled piezoelectric beams with dynamic electro-magnetic field. ESAIM Control Optim. Calc. Var. 24–8, 1–24 (2019)
Pazy, A.: Semigroups of Linear Operators and Applications to Partial Differential Equations. Springer, New York (1983)
PI Inc., Piezolectric Actuators: components, technologies, operation, http://www.piezo.ws/pdf/Piezo.pdf, p. 63
Ramos, A.J.A., Goncalves, C.S.L., Neto, S.S.C.: Exponential stability and numerical treatment for piezoelectric beams with magnetic effect. ESAIM 52(1), 255–284 (2018)
Rogacheva, N.: The Theory of Piezoelectric Shells and Plates. CRC Press, Boca Raton (1994)
Ru, C., Liu, X., Sun, Y.: Nanopositioning Technologies: Fundamentals and Applications. Springer International, Berlin (2016)
Russell, D.L.: A General framework for the study of indirect damping mechanisms in elastic systems. J. Math. Anal. Appl. 173–2, 339–358 (1993)
Selberg, S., Tesfahun, A.: Finite-energy global well-posedness of the maxwell-klein-gordon system in lorenz gauge. Commun. Partial Differ. Equ. 35–6, 1029–1057 (2010)
Ultrasound imaging of the brain and liver, Science Daily Magazine, Source: Acoustical Society of America, 26 June 2017. Accessed: 5 January 2018
Shubov, M.: Spectral analysis of a non-selfadjoint operator generated by an energy harvesting model and application to an exact controllability problem. Asymptot. Anal. 102(3–4), 119–156 (2017)
Smith, R.C.: Smart Material Systems. Society for Industrial and Applied Mathematics, Philadelphia (2005)
Tiersten, H.F.: Linear Piezoelectric Plate Vibrations. Plenum Press, New York (1969)
Tebou, L.T., Zuazua, E.: Uniform boundary stabilization of the finite difference space discretization of the 1-d wave equation. Adv. Comput. Math. 26, 337–365 (2007)
Tzou, H.S.: Piezoelectric shells, Solid Mechanics and Its Applications. Kluwer Academic, Dordrecht (1993)
Ultrasonic Welding: Basic Knowledge, Beijing Ultrasonic [PDF document], uploaded on 07 Jan 2013.
Voss, T., Scherpen, J.M.A.: Stabilization and shape control of a 1-D piezoelectric timoshenko beam. Automatica 47, 2780–2785 (2011)
Voss, T., Scherpen, J.M.A.: Port-hamiltonian modeling of a nonlinear timoshenko beam with piezo actuation. SIAM J. Control Optim. 52–1, 493–519 (2014)
Voss, T., Sherpen, J.M.A.: Stability analysis of piezoelectric beams. In 50th IEEE Conference on Decision and Control and European Control Conference, 3758–3763 (2011)
Yang, J.: Special Topics in the Theory of Piezoelectricity. Springer, New York (2009)
Yang, G., Du, J., Wang, J., Yang, J.: Frequency dependence of electromagnetic radiation from a finite vibrating piezoelectric body. Mech. Res. Commun. 93, 163–168 (2018)
Wehbe, A.: Rational energy decay rate in a wave equation with dynamical control. Appl. Math. Lett. 16–3, 357–364 (2003)
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The financial support of KY NSF EPSCoR Grant (\(\#1514712-1\)) and Western Kentucky University RCAP grant (\(\#20-8038\)) for this research is gratefully acknowledged.
Rights and permissions
About this article
Cite this article
Özer, A.Ö. Stabilization Results for Well-Posed Potential Formulations of a Current-Controlled Piezoelectric Beam and Their Approximations. Appl Math Optim 84, 877–914 (2021). https://doi.org/10.1007/s00245-020-09665-4
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00245-020-09665-4