tdehaeze/digital-brain
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Compare commits

base: tdehaeze:54103360a2bc5f23ff05334afa6f4632747611c5
Branches Tags
tdehaeze:master
..
compare: tdehaeze:38343a7d7519c84c3273392f1c8249c40cd2d87e
Branches Tags
tdehaeze:master

2 Commits
54103360a2 ... 38343a7d75

Author SHA1 Message Date
Thomas Dehaeze
38343a7d75 Update Content - 2021-04-15 2021-04-15 16:00:04 +02:00
Thomas Dehaeze
337dc99605 Update Content - 2021-04-14 2021-04-14 17:48:36 +02:00
6 changed files with 145 additions and 60 deletions
Show Stats Expand all files Collapse all files

90
content/book/preumont18_vibrat_contr_activ_struc_fourt_edition.md
View File

@ -8,7 +8,7 @@ Tags
: [Vibration Isolation]({{< relref "vibration_isolation" >}}), [Reference Books]({{< relref "reference_books" >}}), [Stewart Platforms]({{< relref "stewart_platforms" >}}), [HAC-HAC]({{< relref "hac_hac" >}}) : [Vibration Isolation]({{< relref "vibration_isolation" >}}), [Reference Books]({{< relref "reference_books" >}}), [Stewart Platforms]({{< relref "stewart_platforms" >}}), [HAC-HAC]({{< relref "hac_hac" >}})
Reference Reference
: ([Preumont 2018](#orgaa0487d)) : ([Preumont 2018](#orgd83c544))
Author(s) Author(s)
: Preumont, A. : Preumont, A.
@ -61,11 +61,11 @@ There are two radically different approached to disturbance rejection: feedback
#### Feedback {#feedback} #### Feedback {#feedback}
<a id="org5636ea9"></a> <a id="orgda21dda"></a>
{{< figure src="/ox-hugo/preumont18_classical_feedback_small.png" caption="Figure 1: Principle of feedback control" >}} {{< figure src="/ox-hugo/preumont18_classical_feedback_small.png" caption="Figure 1: Principle of feedback control" >}}
The principle of feedback is represented on figure [1](#org5636ea9). The output \\(y\\) of the system is compared to the reference signal \\(r\\), and the error signal \\(\epsilon = r-y\\) is passed into a compensator \\(K(s)\\) and applied to the system \\(G(s)\\), \\(d\\) is the disturbance. The principle of feedback is represented on figure [1](#orgda21dda). The output \\(y\\) of the system is compared to the reference signal \\(r\\), and the error signal \\(\epsilon = r-y\\) is passed into a compensator \\(K(s)\\) and applied to the system \\(G(s)\\), \\(d\\) is the disturbance.
The design problem consists of finding the appropriate compensator \\(K(s)\\) such that the closed-loop system is stable and behaves in the appropriate manner. The design problem consists of finding the appropriate compensator \\(K(s)\\) such that the closed-loop system is stable and behaves in the appropriate manner.
In the control of lightly damped structures, feedback control is used for two distinct and complementary purposes: **active damping** and **model-based feedback**. In the control of lightly damped structures, feedback control is used for two distinct and complementary purposes: **active damping** and **model-based feedback**.
@ -87,12 +87,12 @@ The objective is to control a variable \\(y\\) to a desired value \\(r\\) in spi
#### Feedforward {#feedforward} #### Feedforward {#feedforward}
<a id="org88ce537"></a> <a id="orgf75c047"></a>
{{< figure src="/ox-hugo/preumont18_feedforward_adaptative.png" caption="Figure 2: Principle of feedforward control" >}} {{< figure src="/ox-hugo/preumont18_feedforward_adaptative.png" caption="Figure 2: Principle of feedforward control" >}}
The method relies on the availability of a **reference signal correlated to the primary disturbance**. The method relies on the availability of a **reference signal correlated to the primary disturbance**.
The idea is to produce a second disturbance such that is cancels the effect of the primary disturbance at the location of the sensor error. Its principle is explained in figure [2](#org88ce537). The idea is to produce a second disturbance such that is cancels the effect of the primary disturbance at the location of the sensor error. Its principle is explained in figure [2](#orgf75c047).
The filter coefficients are adapted in such a way that the error signal at one or several critical points is minimized. The filter coefficients are adapted in such a way that the error signal at one or several critical points is minimized.
@ -123,11 +123,11 @@ The table [1](#table--tab:adv-dis-type-control) summarizes the main features of
### The Various Steps of the Design {#the-various-steps-of-the-design} ### The Various Steps of the Design {#the-various-steps-of-the-design}
<a id="org0685157"></a> <a id="org1939c0d"></a>
{{< figure src="/ox-hugo/preumont18_design_steps.png" caption="Figure 3: The various steps of the design" >}} {{< figure src="/ox-hugo/preumont18_design_steps.png" caption="Figure 3: The various steps of the design" >}}
The various steps of the design of a controlled structure are shown in figure [3](#org0685157). The various steps of the design of a controlled structure are shown in figure [3](#org1939c0d).
The **starting point** is: The **starting point** is:
@ -154,14 +154,14 @@ If the dynamics of the sensors and actuators may significantly affect the behavi
### Plant Description, Error and Control Budget {#plant-description-error-and-control-budget} ### Plant Description, Error and Control Budget {#plant-description-error-and-control-budget}
From the block diagram of the control system (figure [4](#org23c9634)): From the block diagram of the control system (figure [4](#orgaf01f6c)):
\begin{align\*} \begin{align\*}
y &= (I - G\_{yu}H)^{-1} G\_{yw} w\\\\\\ y &= (I - G\_{yu}H)^{-1} G\_{yw} w\\\\\\
z &= T\_{zw} w = [G\_{zw} + G\_{zu}H(I - G\_{yu}H)^{-1} G\_{yw}] w z &= T\_{zw} w = [G\_{zw} + G\_{zu}H(I - G\_{yu}H)^{-1} G\_{yw}] w
\end{align\*} \end{align\*}
<a id="org23c9634"></a> <a id="orgaf01f6c"></a>
{{< figure src="/ox-hugo/preumont18_general_plant.png" caption="Figure 4: Block diagram of the control System" >}} {{< figure src="/ox-hugo/preumont18_general_plant.png" caption="Figure 4: Block diagram of the control System" >}}
@ -186,12 +186,12 @@ Even more interesting for the design is the **Cumulative Mean Square** response
It is a monotonously decreasing function of frequency and describes the contribution of all frequencies above \\(\omega\\) to the mean-square value of \\(z\\). It is a monotonously decreasing function of frequency and describes the contribution of all frequencies above \\(\omega\\) to the mean-square value of \\(z\\).
\\(\sigma\_z(0)\\) is then the global RMS response. \\(\sigma\_z(0)\\) is then the global RMS response.
A typical plot of \\(\sigma\_z(\omega)\\) is shown figure [5](#orgc0a0d3d). A typical plot of \\(\sigma\_z(\omega)\\) is shown figure [5](#org7ddcf2a).
It is useful to **identify the critical modes** in a design, at which the effort should be targeted. It is useful to **identify the critical modes** in a design, at which the effort should be targeted.
The diagram can also be used to **assess the control laws** and compare different actuator and sensor configuration. The diagram can also be used to **assess the control laws** and compare different actuator and sensor configuration.
<a id="orgc0a0d3d"></a> <a id="org7ddcf2a"></a>
{{< figure src="/ox-hugo/preumont18_cas_plot.png" caption="Figure 5: Error budget distribution in OL and CL for increasing gains" >}} {{< figure src="/ox-hugo/preumont18_cas_plot.png" caption="Figure 5: Error budget distribution in OL and CL for increasing gains" >}}
@ -304,7 +304,7 @@ The mode shapes are orthogonal with respect to the stiffness and mass matrices:
With \\(\mu\_i\\) the **modal mass** (also called the generalized mass) of mode \\(i\\). With \\(\mu\_i\\) the **modal mass** (also called the generalized mass) of mode \\(i\\).
### Modal Decomposition {#modal-decomposition} ### [Modal Decomposition]({{< relref "modal_decomposition" >}}) {#modal-decomposition--modal-decomposition-dot-md}
#### Structure Without Rigid Body Modes {#structure-without-rigid-body-modes} #### Structure Without Rigid Body Modes {#structure-without-rigid-body-modes}
@ -398,11 +398,11 @@ With:
D\_i(\omega) = \frac{1}{1 - \omega^2/\omega\_i^2 + 2 j \xi\_i \omega/\omega\_i} D\_i(\omega) = \frac{1}{1 - \omega^2/\omega\_i^2 + 2 j \xi\_i \omega/\omega\_i}
\end{equation} \end{equation}
<a id="org960ee21"></a> <a id="orga618336"></a>
{{< figure src="/ox-hugo/preumont18_neglected_modes.png" caption="Figure 6: Fourier spectrum of the excitation \\(F\\) and dynamic amplitification \\(D\_i\\) of mode \\(i\\) and \\(k\\) such that \\(\omega\_i < \omega\_b\\) and \\(\omega\_k \gg \omega\_b\\)" >}} {{< figure src="/ox-hugo/preumont18_neglected_modes.png" caption="Figure 6: Fourier spectrum of the excitation \\(F\\) and dynamic amplitification \\(D\_i\\) of mode \\(i\\) and \\(k\\) such that \\(\omega\_i < \omega\_b\\) and \\(\omega\_k \gg \omega\_b\\)" >}}
If the excitation has a limited bandwidth \\(\omega\_b\\), the contribution of the high frequency modes \\(\omega\_k \gg \omega\_b\\) can be evaluated by assuming \\(D\_k(\omega) \approx 1\\) (as shown on figure [6](#org960ee21)). If the excitation has a limited bandwidth \\(\omega\_b\\), the contribution of the high frequency modes \\(\omega\_k \gg \omega\_b\\) can be evaluated by assuming \\(D\_k(\omega) \approx 1\\) (as shown on figure [6](#orga618336)).
And \\(G(\omega)\\) can be rewritten on terms of the **low frequency modes only**: And \\(G(\omega)\\) can be rewritten on terms of the **low frequency modes only**:
\\[ G(\omega) \approx \sum\_{i=1}^m \frac{\phi\_i \phi\_i^T}{\mu\_i \omega\_i^2} D\_i(\omega) + R \\] \\[ G(\omega) \approx \sum\_{i=1}^m \frac{\phi\_i \phi\_i^T}{\mu\_i \omega\_i^2} D\_i(\omega) + R \\]
@ -441,9 +441,9 @@ The open-loop FRF of a collocated system corresponds to a diagonal component of
If we assumes that the collocated system is undamped and is attached to the DoF \\(k\\), the open-loop FRF is purely real: If we assumes that the collocated system is undamped and is attached to the DoF \\(k\\), the open-loop FRF is purely real:
\\[ G\_{kk}(\omega) = \sum\_{i=1}^m \frac{\phi\_i^2(k)}{\mu\_i (\omega\_i^2 - \omega^2)} + R\_{kk} \\] \\[ G\_{kk}(\omega) = \sum\_{i=1}^m \frac{\phi\_i^2(k)}{\mu\_i (\omega\_i^2 - \omega^2)} + R\_{kk} \\]
\\(G\_{kk}\\) is a monotonously increasing function of \\(\omega\\) (figure [7](#orgcc3baba)). \\(G\_{kk}\\) is a monotonously increasing function of \\(\omega\\) (figure [7](#orgecdb253)).
<a id="orgcc3baba"></a> <a id="orgecdb253"></a>
{{< figure src="/ox-hugo/preumont18_collocated_control_frf.png" caption="Figure 7: Open-Loop FRF of an undamped structure with collocated actuator/sensor pair" >}} {{< figure src="/ox-hugo/preumont18_collocated_control_frf.png" caption="Figure 7: Open-Loop FRF of an undamped structure with collocated actuator/sensor pair" >}}
@ -457,9 +457,9 @@ For lightly damped structure, the poles and zeros are just moved a little bit in
</div> </div>
If the undamped structure is excited harmonically by the actuator at the frequency of the transmission zero \\(z\_i\\), the amplitude of the response of the collocated sensor vanishes. That means that the structure oscillates at the frequency \\(z\_i\\) according to the mode shape shown in dotted line figure [8](#org2ea7272). If the undamped structure is excited harmonically by the actuator at the frequency of the transmission zero \\(z\_i\\), the amplitude of the response of the collocated sensor vanishes. That means that the structure oscillates at the frequency \\(z\_i\\) according to the mode shape shown in dotted line figure [8](#org2e6ee6b).
<a id="org2ea7272"></a> <a id="org2e6ee6b"></a>
{{< figure src="/ox-hugo/preumont18_collocated_zero.png" caption="Figure 8: Structure with collocated actuator and sensor" >}} {{< figure src="/ox-hugo/preumont18_collocated_zero.png" caption="Figure 8: Structure with collocated actuator and sensor" >}}
@ -474,9 +474,9 @@ The open-loop poles are independant of the actuator and sensor configuration whi
</div> </div>
By looking at figure [7](#orgcc3baba), we see that neglecting the residual mode in the modelling amounts to translating the FRF diagram vertically. That produces a shift in the location of the transmission zeros to the right. By looking at figure [7](#orgecdb253), we see that neglecting the residual mode in the modelling amounts to translating the FRF diagram vertically. That produces a shift in the location of the transmission zeros to the right.
<a id="orga2c2292"></a> <a id="org8e5acfb"></a>
{{< figure src="/ox-hugo/preumont18_alternating_p_z.png" caption="Figure 9: Bode plot of a lighly damped structure with collocated actuator and sensor" >}} {{< figure src="/ox-hugo/preumont18_alternating_p_z.png" caption="Figure 9: Bode plot of a lighly damped structure with collocated actuator and sensor" >}}
@ -486,7 +486,7 @@ The open-loop transfer function of a lighly damped structure with a collocated a
G(s) = G\_0 \frac{\Pi\_i(s^2/z\_i^2 + 2 \xi\_i s/z\_i + 1)}{\Pi\_j(s^2/\omega\_j^2 + 2 \xi\_j s /\omega\_j + 1)} G(s) = G\_0 \frac{\Pi\_i(s^2/z\_i^2 + 2 \xi\_i s/z\_i + 1)}{\Pi\_j(s^2/\omega\_j^2 + 2 \xi\_j s /\omega\_j + 1)}
\end{equation} \end{equation}
The corresponding Bode plot is represented in figure [9](#orga2c2292). Every imaginary pole at \\(\pm j\omega\_i\\) introduces a \\(\SI{180}{\degree}\\) phase lag and every imaginary zero at \\(\pm jz\_i\\) introduces a phase lead of \\(\SI{180}{\degree}\\). The corresponding Bode plot is represented in figure [9](#org8e5acfb). Every imaginary pole at \\(\pm j\omega\_i\\) introduces a \\(\SI{180}{\degree}\\) phase lag and every imaginary zero at \\(\pm jz\_i\\) introduces a phase lead of \\(\SI{180}{\degree}\\).
In this way, the phase diagram is always contained between \\(\SI{0}{\degree}\\) and \\(\SI{-180}{\degree}\\) as a consequence of the interlacing property. In this way, the phase diagram is always contained between \\(\SI{0}{\degree}\\) and \\(\SI{-180}{\degree}\\) as a consequence of the interlacing property.
@ -508,12 +508,12 @@ Two broad categories of actuators can be distinguish:
A voice coil transducer is an energy transformer which converts electrical power into mechanical power and vice versa. A voice coil transducer is an energy transformer which converts electrical power into mechanical power and vice versa.
The system consists of (see figure [10](#orgd4ab71d)): The system consists of (see figure [10](#org5b9842b)):
- A permanent magnet which produces a uniform flux density \\(B\\) normal to the gap - A permanent magnet which produces a uniform flux density \\(B\\) normal to the gap
- A coil which is free to move axially - A coil which is free to move axially
<a id="orgd4ab71d"></a> <a id="org5b9842b"></a>
{{< figure src="/ox-hugo/preumont18_voice_coil_schematic.png" caption="Figure 10: Physical principle of a voice coil transducer" >}} {{< figure src="/ox-hugo/preumont18_voice_coil_schematic.png" caption="Figure 10: Physical principle of a voice coil transducer" >}}
@ -551,9 +551,9 @@ Thus, at any time, there is an equilibrium between the electrical power absorbed
#### Proof-Mass Actuator {#proof-mass-actuator} #### Proof-Mass Actuator {#proof-mass-actuator}
A reaction mass \\(m\\) is conected to the support structure by a spring \\(k\\) , and damper \\(c\\) and a force actuator \\(f = T i\\) (figure [11](#org89e3371)). A reaction mass \\(m\\) is conected to the support structure by a spring \\(k\\) , and damper \\(c\\) and a force actuator \\(f = T i\\) (figure [11](#org608f53f)).
<a id="org89e3371"></a> <a id="org608f53f"></a>
{{< figure src="/ox-hugo/preumont18_proof_mass_actuator.png" caption="Figure 11: Proof-mass actuator" >}} {{< figure src="/ox-hugo/preumont18_proof_mass_actuator.png" caption="Figure 11: Proof-mass actuator" >}}
@ -583,9 +583,9 @@ with:
</div> </div>
Above some critical frequency \\(\omega\_c \approx 2\omega\_p\\), **the proof-mass actuator can be regarded as an ideal force generator** (figure [12](#org1b03971)). Above some critical frequency \\(\omega\_c \approx 2\omega\_p\\), **the proof-mass actuator can be regarded as an ideal force generator** (figure [12](#org21ce10b)).
<a id="org1b03971"></a> <a id="org21ce10b"></a>
{{< figure src="/ox-hugo/preumont18_proof_mass_tf.png" caption="Figure 12: Bode plot \\(F/i\\) of the proof-mass actuator" >}} {{< figure src="/ox-hugo/preumont18_proof_mass_tf.png" caption="Figure 12: Bode plot \\(F/i\\) of the proof-mass actuator" >}}
@ -610,7 +610,7 @@ By using the two equations, we obtain:
Above the corner frequency, the gain of the geophone is equal to the transducer constant \\(T\\). Above the corner frequency, the gain of the geophone is equal to the transducer constant \\(T\\).
<a id="org8e3e3a9"></a> <a id="orgb548c88"></a>
{{< figure src="/ox-hugo/preumont18_geophone.png" caption="Figure 13: Model of a geophone based on a voice coil transducer" >}} {{< figure src="/ox-hugo/preumont18_geophone.png" caption="Figure 13: Model of a geophone based on a voice coil transducer" >}}
@ -619,9 +619,9 @@ Designing geophones with very low corner frequency is in general difficult. Acti
### General Electromechanical Transducer {#general-electromechanical-transducer} ### General Electromechanical Transducer {#general-electromechanical-transducer}
The consitutive behavior of a wide class of electromechanical transducers can be modelled as in figure [14](#org8d49672). The consitutive behavior of a wide class of electromechanical transducers can be modelled as in figure [14](#org98492c9).
<a id="org8d49672"></a> <a id="org98492c9"></a>
{{< figure src="/ox-hugo/preumont18_electro_mechanical_transducer.png" caption="Figure 14: Electrical analog representation of an electromechanical transducer" >}} {{< figure src="/ox-hugo/preumont18_electro_mechanical_transducer.png" caption="Figure 14: Electrical analog representation of an electromechanical transducer" >}}
@ -646,7 +646,7 @@ With:
Equation \eqref{eq:gen_trans_e} shows that the voltage across the electrical terminals of any electromechanical transducer is the sum of a contribution proportional to the current applied and a contribution proportional to the velocity of the mechanical terminals. Equation \eqref{eq:gen_trans_e} shows that the voltage across the electrical terminals of any electromechanical transducer is the sum of a contribution proportional to the current applied and a contribution proportional to the velocity of the mechanical terminals.
Thus, if \\(Z\_ei\\) can be measured and substracted from \\(e\\), a signal proportional to the velocity is obtained. Thus, if \\(Z\_ei\\) can be measured and substracted from \\(e\\), a signal proportional to the velocity is obtained.
To do so, the bridge circuit as shown on figure [15](#org9077cf9) can be used. To do so, the bridge circuit as shown on figure [15](#org3b85763) can be used.
We can show that We can show that
@ -656,7 +656,7 @@ We can show that
which is indeed a linear function of the velocity \\(v\\) at the mechanical terminals. which is indeed a linear function of the velocity \\(v\\) at the mechanical terminals.
<a id="org9077cf9"></a> <a id="org3b85763"></a>
{{< figure src="/ox-hugo/preumont18_bridge_circuit.png" caption="Figure 15: Bridge circuit for self-sensing actuation" >}} {{< figure src="/ox-hugo/preumont18_bridge_circuit.png" caption="Figure 15: Bridge circuit for self-sensing actuation" >}}
@ -664,9 +664,9 @@ which is indeed a linear function of the velocity \\(v\\) at the mechanical term
### Smart Materials {#smart-materials} ### Smart Materials {#smart-materials}
Smart materials have the ability to respond significantly to stimuli of different physical nature. Smart materials have the ability to respond significantly to stimuli of different physical nature.
Figure [16](#orga08bcd9) lists various effects that are observed in materials in response to various inputs. Figure [16](#org6279c77) lists various effects that are observed in materials in response to various inputs.
<a id="orga08bcd9"></a> <a id="org6279c77"></a>
{{< figure src="/ox-hugo/preumont18_smart_materials.png" caption="Figure 16: Stimulus response relations indicating various effects in materials. The smart materials corresponds to the non-diagonal cells" >}} {{< figure src="/ox-hugo/preumont18_smart_materials.png" caption="Figure 16: Stimulus response relations indicating various effects in materials. The smart materials corresponds to the non-diagonal cells" >}}
@ -761,7 +761,7 @@ It measures the efficiency of the conversion of the mechanical energy into elect
</div> </div>
If one assumes that all the electrical and mechanical quantities are uniformly distributed in a linear transducer formed by a **stack** (see figure [17](#orge7aeb11)) of \\(n\\) disks of thickness \\(t\\) and cross section \\(A\\), the global constitutive equations of the transducer are obtained by integrating \eqref{eq:piezo_eq_matrix_bis} over the volume of the transducer: If one assumes that all the electrical and mechanical quantities are uniformly distributed in a linear transducer formed by a **stack** (see figure [17](#org8006b4a)) of \\(n\\) disks of thickness \\(t\\) and cross section \\(A\\), the global constitutive equations of the transducer are obtained by integrating \eqref{eq:piezo_eq_matrix_bis} over the volume of the transducer:
\begin{equation} \begin{equation}
\begin{bmatrix}Q\\\Delta\end{bmatrix} \begin{bmatrix}Q\\\Delta\end{bmatrix}
@ -782,7 +782,7 @@ where
- \\(C = \epsilon^T A n^2/l\\) is the capacitance of the transducer with no external load (\\(f = 0\\)) - \\(C = \epsilon^T A n^2/l\\) is the capacitance of the transducer with no external load (\\(f = 0\\))
- \\(K\_a = A/s^El\\) is the stiffness with short-circuited electrodes (\\(V = 0\\)) - \\(K\_a = A/s^El\\) is the stiffness with short-circuited electrodes (\\(V = 0\\))
<a id="orge7aeb11"></a> <a id="org8006b4a"></a>
{{< figure src="/ox-hugo/preumont18_piezo_stack.png" caption="Figure 17: Piezoelectric linear transducer" >}} {{< figure src="/ox-hugo/preumont18_piezo_stack.png" caption="Figure 17: Piezoelectric linear transducer" >}}
@ -802,7 +802,7 @@ Equation \eqref{eq:piezo_stack_eq} can be inverted to obtain
#### Energy Stored in the Piezoelectric Transducer {#energy-stored-in-the-piezoelectric-transducer} #### Energy Stored in the Piezoelectric Transducer {#energy-stored-in-the-piezoelectric-transducer}
Let us write the total stored electromechanical energy of a discrete piezoelectric transducer as shown on figure [18](#org62300cf). Let us write the total stored electromechanical energy of a discrete piezoelectric transducer as shown on figure [18](#org7c30411).
The total power delivered to the transducer is the sum of electric power \\(V i\\) and the mechanical power \\(f \dot{\Delta}\\). The net work of the transducer is The total power delivered to the transducer is the sum of electric power \\(V i\\) and the mechanical power \\(f \dot{\Delta}\\). The net work of the transducer is
@ -810,7 +810,7 @@ The total power delivered to the transducer is the sum of electric power \\(V i\
dW = V i dt + f \dot{\Delta} dt = V dQ + f d\Delta dW = V i dt + f \dot{\Delta} dt = V dQ + f d\Delta
\end{equation} \end{equation}
<a id="org62300cf"></a> <a id="org7c30411"></a>
{{< figure src="/ox-hugo/preumont18_piezo_discrete.png" caption="Figure 18: Discrete Piezoelectric Transducer" >}} {{< figure src="/ox-hugo/preumont18_piezo_discrete.png" caption="Figure 18: Discrete Piezoelectric Transducer" >}}
@ -844,10 +844,10 @@ The ratio between the remaining stored energy and the initial stored energy is
#### Admittance of the Piezoelectric Transducer {#admittance-of-the-piezoelectric-transducer} #### Admittance of the Piezoelectric Transducer {#admittance-of-the-piezoelectric-transducer}
Consider the system of figure [19](#orga98ecb7), where the piezoelectric transducer is assumed massless and is connected to a mass \\(M\\). Consider the system of figure [19](#org5060008), where the piezoelectric transducer is assumed massless and is connected to a mass \\(M\\).
The force acting on the mass is negative of that acting on the transducer, \\(f = -M \ddot{x}\\). The force acting on the mass is negative of that acting on the transducer, \\(f = -M \ddot{x}\\).
<a id="orga98ecb7"></a> <a id="org5060008"></a>
{{< figure src="/ox-hugo/preumont18_piezo_stack_admittance.png" caption="Figure 19: Elementary dynamical model of the piezoelectric transducer" >}} {{< figure src="/ox-hugo/preumont18_piezo_stack_admittance.png" caption="Figure 19: Elementary dynamical model of the piezoelectric transducer" >}}
@ -866,9 +866,9 @@ And one can see that
\frac{z^2 - p^2}{z^2} = k^2 \frac{z^2 - p^2}{z^2} = k^2
\end{equation} \end{equation}
Equation \eqref{eq:distance_p_z} constitutes a practical way to determine the electromechanical coupling factor from the poles and zeros of the admittance measurement (figure [20](#orge87e33b)). Equation \eqref{eq:distance_p_z} constitutes a practical way to determine the electromechanical coupling factor from the poles and zeros of the admittance measurement (figure [20](#org7f3b3bf)).
<a id="orge87e33b"></a> <a id="org7f3b3bf"></a>
{{< figure src="/ox-hugo/preumont18_piezo_admittance_curve.png" caption="Figure 20: Typical admittance FRF of the transducer" >}} {{< figure src="/ox-hugo/preumont18_piezo_admittance_curve.png" caption="Figure 20: Typical admittance FRF of the transducer" >}}
@ -1566,7 +1566,7 @@ Their design requires a model of the structure, and there is usually a trade-off
When collocated actuator/sensor pairs can be used, stability can be achieved using positivity concepts, but in many situations, collocated pairs are not feasible for HAC. When collocated actuator/sensor pairs can be used, stability can be achieved using positivity concepts, but in many situations, collocated pairs are not feasible for HAC.
The HAC/LAC approach consist of combining the two approached in a dual-loop control as shown in Figure [21](#orgeb43e36). The HAC/LAC approach consist of combining the two approached in a dual-loop control as shown in Figure [21](#org62e1395).
The inner loop uses a set of collocated actuator/sensor pairs for decentralized active damping with guaranteed stability ; the outer loop consists of a non-collocated HAC based on a model of the actively damped structure. The inner loop uses a set of collocated actuator/sensor pairs for decentralized active damping with guaranteed stability ; the outer loop consists of a non-collocated HAC based on a model of the actively damped structure.
This approach has the following advantages: This approach has the following advantages:
@ -1574,7 +1574,7 @@ This approach has the following advantages:
- The active damping makes it easier to gain-stabilize the modes outside the bandwidth of the output loop (improved gain margin) - The active damping makes it easier to gain-stabilize the modes outside the bandwidth of the output loop (improved gain margin)
- The larger damping of the modes within the controller bandwidth makes them more robust to the parmetric uncertainty (improved phase margin) - The larger damping of the modes within the controller bandwidth makes them more robust to the parmetric uncertainty (improved phase margin)
<a id="orgeb43e36"></a> <a id="org62e1395"></a>
{{< figure src="/ox-hugo/preumont18_hac_lac_control.png" caption="Figure 21: Principle of the dual-loop HAC/LAC control" >}} {{< figure src="/ox-hugo/preumont18_hac_lac_control.png" caption="Figure 21: Principle of the dual-loop HAC/LAC control" >}}
@ -1818,4 +1818,4 @@ This approach has the following advantages:
## Bibliography {#bibliography} ## Bibliography {#bibliography}
<a id="orgaa0487d"></a>Preumont, Andre. 2018. _Vibration Control of Active Structures - Fourth Edition_. Solid Mechanics and Its Applications. Springer International Publishing. <https://doi.org/10.1007/978-3-319-72296-2>. <a id="orgd83c544"></a>Preumont, Andre. 2018. _Vibration Control of Active Structures - Fourth Edition_. Solid Mechanics and Its Applications. Springer International Publishing. <https://doi.org/10.1007/978-3-319-72296-2>.

115
content/phdthesis/rankers98_machin.md
View File

@ -8,7 +8,7 @@ Tags
: [Finite Element Model]({{< relref "finite_element_model" >}}) : [Finite Element Model]({{< relref "finite_element_model" >}})
Reference Reference
: ([Rankers 1998](#org1e7e43d)) : ([Rankers 1998](#org98ff031))
Author(s) Author(s)
: Rankers, A. M. : Rankers, A. M.
@ -163,13 +163,13 @@ The basic questions that are addressed in this thesis are:
### Basic Control Aspects {#basic-control-aspects} ### Basic Control Aspects {#basic-control-aspects}
A block diagram representation of a typical servo-system is shown in Figure [1](#org8d51e97). A block diagram representation of a typical servo-system is shown in Figure [1](#orgfda2012).
The main task of the system is achieve a desired positional relation between two or more components of the system. The main task of the system is achieve a desired positional relation between two or more components of the system.
Therefore, a sensor measures the position which is then compared to the desired value, and the resulting error is used to generate correcting forces. Therefore, a sensor measures the position which is then compared to the desired value, and the resulting error is used to generate correcting forces.
In most systems, the "actual output" (e.g. position of end-effector) cannot be measured directly, and the feedback will therefore be based on a "measured output" (e.g. encoder signal at the motor). In most systems, the "actual output" (e.g. position of end-effector) cannot be measured directly, and the feedback will therefore be based on a "measured output" (e.g. encoder signal at the motor).
It is important to realize that these two outputs can differ, first due to resilience in the mechanical system, and second because of geometrical imperfections in the mechanical transmission between motor and end-effector. It is important to realize that these two outputs can differ, first due to resilience in the mechanical system, and second because of geometrical imperfections in the mechanical transmission between motor and end-effector.
<a id="org8d51e97"></a> <a id="orgfda2012"></a>
{{< figure src="/ox-hugo/rankers98_basic_el_mech_servo.png" caption="Figure 1: Basic elements of mechanical servo system" >}} {{< figure src="/ox-hugo/rankers98_basic_el_mech_servo.png" caption="Figure 1: Basic elements of mechanical servo system" >}}
@ -180,10 +180,10 @@ The correction force \\(F\\) is defined by:
F = k\_p \epsilon + k\_d \dot{\epsilon} + k\_i \int \epsilon dt F = k\_p \epsilon + k\_d \dot{\epsilon} + k\_i \int \epsilon dt
\end{equation} \end{equation}
It is illustrative to see that basically the proportional and derivative part of such a position control loop is very similar to a mechanical spring and damper that connect two points (Figure [2](#org7b90722)). It is illustrative to see that basically the proportional and derivative part of such a position control loop is very similar to a mechanical spring and damper that connect two points (Figure [2](#org2d694a7)).
If \\(c\\) and \\(d\\) represent the constant mechanical stiffness and damping between points \\(A\\) and \\(B\\), and a reference position profile \\(h(t)\\) is applied at \\(A\\), then an opposing force \\(F\\) is generated as soon as the position \\(x\\) and speed \\(\dot{x}\\) of point \\(B\\) does not correspond to \\(h(t)\\) and \\(\dot{h}(t)\\). If \\(c\\) and \\(d\\) represent the constant mechanical stiffness and damping between points \\(A\\) and \\(B\\), and a reference position profile \\(h(t)\\) is applied at \\(A\\), then an opposing force \\(F\\) is generated as soon as the position \\(x\\) and speed \\(\dot{x}\\) of point \\(B\\) does not correspond to \\(h(t)\\) and \\(\dot{h}(t)\\).
<a id="org7b90722"></a> <a id="org2d694a7"></a>
{{< figure src="/ox-hugo/rankers98_basic_elastic_struct.png" caption="Figure 2: Basic Elastic Structure" >}} {{< figure src="/ox-hugo/rankers98_basic_elastic_struct.png" caption="Figure 2: Basic Elastic Structure" >}}
@ -199,9 +199,9 @@ These properties are very essential since they introduce the issue of **servo st
An important aspect of a feedback controller is the fact that control forces can only result from an error signal. An important aspect of a feedback controller is the fact that control forces can only result from an error signal.
Thus any desired set-point profile first leads to a position error before the corresponding driving forces are generated. Thus any desired set-point profile first leads to a position error before the corresponding driving forces are generated.
Most modern servo-systems have not only a feedback section, but also a feedforward section, as indicated in Figure [3](#org4adf108). Most modern servo-systems have not only a feedback section, but also a feedforward section, as indicated in Figure [3](#org0c46a44).
<a id="org4adf108"></a> <a id="org0c46a44"></a>
{{< figure src="/ox-hugo/rankers98_feedforward_example.png" caption="Figure 3: Mechanical servo system with feedback and feedforward control" >}} {{< figure src="/ox-hugo/rankers98_feedforward_example.png" caption="Figure 3: Mechanical servo system with feedback and feedforward control" >}}
@ -246,9 +246,9 @@ Basically, machine dynamics can have two deterioration effects in mechanical ser
#### Actuator Flexibility {#actuator-flexibility} #### Actuator Flexibility {#actuator-flexibility}
The basic characteristics of what is called "actuator flexibility" is the fact that in the frequency range of interest (usually \\(0-10\times \text{bandwidth}\\)) the driven system no longer behaves as one rigid body (Figure [4](#org72d5adf)) due to compliance between the motor and the load. The basic characteristics of what is called "actuator flexibility" is the fact that in the frequency range of interest (usually \\(0-10\times \text{bandwidth}\\)) the driven system no longer behaves as one rigid body (Figure [4](#org765d5f4)) due to compliance between the motor and the load.
<a id="org72d5adf"></a> <a id="org765d5f4"></a>
{{< figure src="/ox-hugo/rankers98_actuator_flexibility.png" caption="Figure 4: Actuator Flexibility" >}} {{< figure src="/ox-hugo/rankers98_actuator_flexibility.png" caption="Figure 4: Actuator Flexibility" >}}
@ -258,9 +258,9 @@ The basic characteristics of what is called "actuator flexibility" is the fact t
The second category of dynamic phenomena results from the limited stiffness of the guiding system in combination with the fact the the device is driven in such a way that it has to rely on the guiding system to suppress motion in an undesired direction (in case of a linear direct drive system this occurs if the driving force is not applied at the center of gravity). The second category of dynamic phenomena results from the limited stiffness of the guiding system in combination with the fact the the device is driven in such a way that it has to rely on the guiding system to suppress motion in an undesired direction (in case of a linear direct drive system this occurs if the driving force is not applied at the center of gravity).
In general, a rigid actuator possesses six degrees of freedom, five of which need to be suppressed by the guiding system in order to leave one mobile degree of freedom. In general, a rigid actuator possesses six degrees of freedom, five of which need to be suppressed by the guiding system in order to leave one mobile degree of freedom.
In the present discussion, a planar actuator with three degrees of freedom will be considered (Figure [5](#orgc8208b1)). In the present discussion, a planar actuator with three degrees of freedom will be considered (Figure [5](#org387cdc7)).
<a id="orgc8208b1"></a> <a id="org387cdc7"></a>
{{< figure src="/ox-hugo/rankers98_guiding_flexibility_planar.png" caption="Figure 5: Planar actuator with guiding system flexibility" >}} {{< figure src="/ox-hugo/rankers98_guiding_flexibility_planar.png" caption="Figure 5: Planar actuator with guiding system flexibility" >}}
@ -280,9 +280,9 @@ The last category of dynamic phenomena results from the limited mass and stiffne
In contrast to many textbooks on mechanics and machine dynamics, it is good practice always to look at the combination of driving force on the moving part, and reaction force on the stationary part, of a positioning device. In contrast to many textbooks on mechanics and machine dynamics, it is good practice always to look at the combination of driving force on the moving part, and reaction force on the stationary part, of a positioning device.
When doing so, one has to consider what the effect of the reaction force on the systems performance will be. When doing so, one has to consider what the effect of the reaction force on the systems performance will be.
In the discussion of the previous two dynamic phenomena, the stationary part of the machine was assumed to be infinitely stiff and heavy, and therefore the effect of the reaction force was negligible. In the discussion of the previous two dynamic phenomena, the stationary part of the machine was assumed to be infinitely stiff and heavy, and therefore the effect of the reaction force was negligible.
However, in general the stationary part is neither infinitely heavy, nor is it connected to its environment with infinite stiffness, so the stationary part will exhibit a resonance that is excited by the reaction forces (Figure [6](#org1252ea7)). However, in general the stationary part is neither infinitely heavy, nor is it connected to its environment with infinite stiffness, so the stationary part will exhibit a resonance that is excited by the reaction forces (Figure [6](#orgb3f73d2)).
<a id="org1252ea7"></a> <a id="orgb3f73d2"></a>
{{< figure src="/ox-hugo/rankers98_limited_m_k_stationary_machine_part.png" caption="Figure 6: Limited Mass and Stiffness of Stationary Machine Part" >}} {{< figure src="/ox-hugo/rankers98_limited_m_k_stationary_machine_part.png" caption="Figure 6: Limited Mass and Stiffness of Stationary Machine Part" >}}
@ -393,7 +393,7 @@ q\_1 \\ q\_2 \\ \vdots \\ q\_n
For the i-th modal coordinate \\(q\_i\\) the equation of motion is: For the i-th modal coordinate \\(q\_i\\) the equation of motion is:
\begin{equation} \begin{equation} \label{eq:eoq\_modal\_i}
m\_i \ddot{q\_i}(t) + k\_i q\_i(t) = \phi\_i^T f(t) m\_i \ddot{q\_i}(t) + k\_i q\_i(t) = \phi\_i^T f(t)
\end{equation} \end{equation}
@ -401,6 +401,35 @@ which is a simple second order differential equation similar to that of a single
Using basic formulae that are derived for a simple mass spring system, one is now able to analyse the time and frequency response of all individual modes. Using basic formulae that are derived for a simple mass spring system, one is now able to analyse the time and frequency response of all individual modes.
Having done that, the total motion of the system can simply be obtained by summing the contributions of all modes. Having done that, the total motion of the system can simply be obtained by summing the contributions of all modes.
Characterisation of the dynamics of a mechanical system in terms of frequency response behavior plays a major role in the stability analysis of the control loop of a mechatronic device.
In such an analysis one is typically interested in the transfer function between a measured displacement \\(x\_l\\) and a force \\(f\_k\\), which acts at the physical DoF \\(x\_k\\).
Applying the principle of modal decomposition, any transfer function can be derived by first calculating the behavior of the individual modes, and then summing all modal contributions.
The contribution of one single mode \\(i\\) to the transfer function \\(x\_l/f\_k\\) can be derived by first considering the response of the modal DoF \\(q\_i\\) to a force vector \\(f\\) with only one non-zero component \\(f\_k\\).
In that case, equation \eqref{eq:eoq_modal_i} is reduced to:
\begin{equation}
m\_i \ddot{q}\_i(t) + k\_i q\_i(t) = \phi\_{ik} f\_k(t)
\end{equation}
After a Laplace transformation and some rearrangement:
\begin{equation}
q\_i(s) = f\_k(s) \frac{\phi\_{ik}}{m\_i s^2 + k\_i}
\end{equation}
Once the modal response \\(q\_i\\) is known, the response of the physical DoF \\(x\_l\\) is found by a simple premultiplication with \\(\phi\_{il}\\), which finally leads to the following expression for the contribution of mode \\(i\\) to the transfer function:
\begin{equation}
\left( \frac{x\_l}{f\_k} \right)\_i = \frac{\phi\_{ik}\phi\_{il}}{m\_i s^2 + k\_i}
\end{equation}
The overall transfer function can be found by summation of the individual modal contributions, which all have the same structure:
\begin{equation}
\left( \frac{x\_l}{f\_k} \right) = \sum\_{i = 1}^n \left( \frac{x\_l}{f\_k} \right)\_i = \sum\_{i = 1}^n \frac{\phi\_{ik} \phi\_{il}}{m\_i s^2 + k\_i}
\end{equation}
### Graphical Representation {#graphical-representation} ### Graphical Representation {#graphical-representation}
@ -422,6 +451,62 @@ Having done that, the total motion of the system can simply be obtained by summi
### Basic Characteristics of Mechanical FRF {#basic-characteristics-of-mechanical-frf} ### Basic Characteristics of Mechanical FRF {#basic-characteristics-of-mechanical-frf}
Consider the position control loop of Figure [7](#orgeee8a5d).
<a id="orgeee8a5d"></a>
{{< figure src="/ox-hugo/rankers98_mechanical_servo_system.png" caption="Figure 7: Mechanical position servo-system" >}}
In the ideal situation the mechanical system behaves as one rigid body with mass \\(m\\), so the mechanical transfer function can be written as:
\begin{equation}
\frac{x\_{\text{servo}}}{F\_{\text{servo}}} = \frac{1}{m s^2}
\end{equation}
<a id="org2ffb1b1"></a>
{{< figure src="/ox-hugo/rankers98_ideal_bode_nyquist.png" caption="Figure 8: FRF of an ideal system with no resonances" >}}
In the case of one extra modal contribution, the equation for the mechanical transfer function needs to be extended with one extra term:
\begin{equation} \label{eq:effect\_one\_mode}
\frac{x\_{\text{servo}}}{F\_{\text{servo}}} = \frac{1}{m s^2} + \frac{\phi\_{i,\text{servo}} \phi\_{i,\text{force}}}{m\_i s^2 + k\_i} = \frac{1}{m s^2} + \frac{\phi\_{i,\text{servo}} \phi\_{i,\text{force}}}{m\_i s^2 + m\_i \omega\_i^2}
\end{equation}
The final transfer function and the exact interaction between the two parts depends on the values of the various parameters.
Let's introduce a variable \\(\alpha\\), which relates the high-frequency contribution of the mode to that of the rigid-body motion:
\begin{equation} \label{eq:alpha}
\alpha = \frac{\frac{\phi\_{i,\text{servo}} \phi\_{i,\text{force}}}{m\_i}}{\frac{1}{m}}
\end{equation}
which simplifies equation \eqref{eq:effect_one_mode} to:
\begin{equation}
\frac{x\_{\text{servo}}}{F\_{\text{servo}}} = \frac{1}{ms^2} + \frac{\alpha}{m s^2 + m \omega\_i^2}
\end{equation}
<a id="orgc9ecb0c"></a>
{{< figure src="/ox-hugo/rankers98_frf_effect_alpha.png" caption="Figure 9: Contribution of rigid-body motion and modal dynamics to the amplitude and phase of FRF for various values of \\(\alpha\\)" >}}
<a id="org5180ee3"></a>
{{< figure src="/ox-hugo/rankers98_final_frf_alpha.png" caption="Figure 10: Bode diagram of final FRF (\\(x\_{\text{servo}}/F\_{\text{servo}}\\)) for six values of \\(\alpha\\)" >}}
#### FRF with "-2 slope/zero/pole/-2 slope" behavior {#frf-with-2-slope-zero-pole-2-slope-behavior}
#### FRF with "-2 slope/pole/zero/-2 slope" behavior {#frf-with-2-slope-pole-zero-2-slope-behavior}
#### FRF with "-2 slope/pole/-4 slope" behavior {#frf-with-2-slope-pole-4-slope-behavior}
#### FRF with "-2 slope/pole/-2 slope" behavior {#frf-with-2-slope-pole-2-slope-behavior}
### Destabilising Effect of Modes {#destabilising-effect-of-modes} ### Destabilising Effect of Modes {#destabilising-effect-of-modes}
@ -468,4 +553,4 @@ Through the enormous performance drive in mechatronics systems, much has been le
## Bibliography {#bibliography} ## Bibliography {#bibliography}
<a id="org1e7e43d"></a>Rankers, Adrian Mathias. 1998. “Machine Dynamics in Mechatronic Systems: An Engineering Approach.” University of Twente. <a id="org98ff031"></a>Rankers, Adrian Mathias. 1998. “Machine Dynamics in Mechatronic Systems: An Engineering Approach.” University of Twente.

BIN
static/ox-hugo/rankers98_final_frf_alpha.png Normal file
View File

Binary file not shown.
Side by Side

After

Width:  |  Height:  |  Size: 565 KiB

BIN
static/ox-hugo/rankers98_frf_effect_alpha.png Normal file
View File

Binary file not shown.
Side by Side

After

Width:  |  Height:  |  Size: 319 KiB

BIN
static/ox-hugo/rankers98_ideal_bode_nyquist.png Normal file
View File

Binary file not shown.
Side by Side

After

Width:  |  Height:  |  Size: 571 KiB

BIN
static/ox-hugo/rankers98_mechanical_servo_system.png Normal file
View File

Binary file not shown.
Side by Side

After

Width:  |  Height:  |  Size: 52 KiB