Update Content - 2021-09-28

content/phdthesis/rankers98_machin.md

@@ -8,7 +8,7 @@ Tags
 : [Finite Element Model]({{<relref "finite_element_model.md#" >}})
 
 Reference
-: ([Rankers 1998](#org6436602))
+: ([Rankers 1998](#orgf6a233a))
 
 Author(s)
 : Rankers, A. M.
@@ -164,13 +164,13 @@ The basic questions that are addressed in this thesis are:
 
 ### Basic Control Aspects {#basic-control-aspects}
 
-A block diagram representation of a typical servo-system is shown in Figure [1](#orga7794da).
+A block diagram representation of a typical servo-system is shown in Figure [1](#orgd21c65c).
 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.
 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.
 
-<a id="orga7794da"></a>
+<a id="orgd21c65c"></a>
 
 {{< figure src="/ox-hugo/rankers98_basic_el_mech_servo.png" caption="Figure 1: Basic elements of mechanical servo system" >}}
 
@@ -181,10 +181,10 @@ The correction force \\(F\\) is defined by:
 F = k\_p \epsilon + k\_d \dot{\epsilon} + k\_i \int \epsilon dt
 \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](#org3f53580)).
+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](#orgcea9a7c)).
 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="org3f53580"></a>
+<a id="orgcea9a7c"></a>
 
 {{< figure src="/ox-hugo/rankers98_basic_elastic_struct.png" caption="Figure 2: Basic Elastic Structure" >}}
 
@@ -200,9 +200,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.
 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](#org4a1de51).
+Most modern servo-systems have not only a feedback section, but also a **feedforward** section, as indicated in Figure [3](#orga68790e).
 
-<a id="org4a1de51"></a>
+<a id="orga68790e"></a>
 
 {{< figure src="/ox-hugo/rankers98_feedforward_example.png" caption="Figure 3: Mechanical servo system with feedback and feedforward control" >}}
 
@@ -247,9 +247,9 @@ Basically, machine dynamics can have two deterioration effects in mechanical ser
 
 #### 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](#orgd374c3a)) 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](#org77ec3b1)) due to compliance between the motor and the load.
 
-<a id="orgd374c3a"></a>
+<a id="org77ec3b1"></a>
 
 {{< figure src="/ox-hugo/rankers98_actuator_flexibility.png" caption="Figure 4: Actuator Flexibility" >}}
 
@@ -259,9 +259,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).
 
 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](#orgc0e4433)).
+In the present discussion, a planar actuator with three degrees of freedom will be considered (Figure [5](#org420ac3c)).
 
-<a id="orgc0e4433"></a>
+<a id="org420ac3c"></a>
 
 {{< figure src="/ox-hugo/rankers98_guiding_flexibility_planar.png" caption="Figure 5: Planar actuator with guiding system flexibility" >}}
 
@@ -281,9 +281,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.
 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.
-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](#org22ae9b4)).
+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](#orgf8032be)).
 
-<a id="org22ae9b4"></a>
+<a id="orgf8032be"></a>
 
 {{< figure src="/ox-hugo/rankers98_limited_m_k_stationary_machine_part.png" caption="Figure 6: Limited Mass and Stiffness of Stationary Machine Part" >}}
 
@@ -296,9 +296,9 @@ The effect of frame vibrations is even worse where the quality of positioning of
 
 To understand and describe the behaviour of a mechanical system in a quantitative way, one usually sets up a model of the system.
 The mathematical description of such a model with a finite number of DoF consists of a set of ordinary differential equations.
-Although in the case of simple systems, such as illustrated in Figure [7](#orgb3fe7a3) these equations may be very understandable, in the case of complex systems, the set of differential equations itself gives only limited insight, and mainly serves as a basis for numerical simulations.
+Although in the case of simple systems, such as illustrated in Figure [7](#org36b7f2f) these equations may be very understandable, in the case of complex systems, the set of differential equations itself gives only limited insight, and mainly serves as a basis for numerical simulations.
 
-<a id="orgb3fe7a3"></a>
+<a id="org36b7f2f"></a>
 
 {{< figure src="/ox-hugo/rankers98_1dof_system.png" caption="Figure 7: Elementary dynamic system" >}}
 
@@ -456,26 +456,26 @@ The overall transfer function can be found by summation of the individual modal
 Due to the equivalence with the differential equations of a single mass spring system, equation \eqref{eq:eoq_modal_i} is often represented by a single mass spring system on which a force \\(f^\prime = \phi\_i^T f\\) acts.
 However, this representation implies an important loss of information because it neglects all information about the mode-shape vector.
 
-Consider the system in Figure [8](#org04e36cc) for which the three mode shapes are depicted in the traditional graphical representation.
+Consider the system in Figure [8](#orgdb6f949) for which the three mode shapes are depicted in the traditional graphical representation.
 In this representation, the physical DoF are located at fixed positions and the mode shapes displacement is indicated by the length of an arrow.
 
-<a id="org04e36cc"></a>
+<a id="orgdb6f949"></a>
 
 {{< figure src="/ox-hugo/rankers98_mode_trad_representation.png" caption="Figure 8: System and traditional graphical representation of modes" >}}
 
-Alternatively, considering that for each mode the mode shape vector defined a constant relation between the various physical DoF, one could also represent a mode shape by a lever (Figure [9](#org17d6c66)).
+Alternatively, considering that for each mode the mode shape vector defined a constant relation between the various physical DoF, one could also represent a mode shape by a lever (Figure [9](#org3b3e340)).
 For each individual mode \\(i\\), each physical DoF \\(x\_k\\) is indicated on the lever at a position with respect to the point of rotation that corresponds to the amplitude and sign of that DoF in the mode shape vector (\\(\phi\_{ik}\\)).
 System with no, very little, or proportional damping exhibit real mode shape vectors, and thus the various DoF each their maximum values at the same moment of the cycle.
 Consequently, the respective DoF can only be in phase or in opposite phase.
 All DoF on the same side of the rotation point have identical phases, whereas DoF on opposite sides have opposite phases.
 
-The modal DoF \\(q\_i\\) can be interpreted as the displacement at a distance "1" from the pivot point (Figure [9](#org17d6c66)).
+The modal DoF \\(q\_i\\) can be interpreted as the displacement at a distance "1" from the pivot point (Figure [9](#org3b3e340)).
 
-<a id="org17d6c66"></a>
+<a id="org3b3e340"></a>
 
 {{< figure src="/ox-hugo/rankers98_mode_new_representation.png" caption="Figure 9: System and new graphical representation of mode-shape" >}}
 
-In the case of a lumped mass model, as in the previous example, it is possible to indicate at each physical DoF on the modal lever the corresponding physical mass, as shown in Figure [10](#org0de4b98) (a).
+In the case of a lumped mass model, as in the previous example, it is possible to indicate at each physical DoF on the modal lever the corresponding physical mass, as shown in Figure [10](#org48d23d9) (a).
 The resulting moment of inertia \\(J\_i\\) of the i-th modal lever then is:
 
 \begin{equation}
@@ -490,20 +490,20 @@ m\_i = \phi\_j^T M \phi\_j = \sum\_{k=1}^n m\_k \phi\_{ik}^2
 
 As a result of this, the modal mass \\(m\_i\\) could be interpreted as the resulting mass moment of inertia of the modal lever, or alternatively as a mass located at a distance "1" from the pivot point.
 
-The transition from physical masses to modal masses is illustrated in Figure [10](#org0de4b98) for the mode 2 of the example system.
+The transition from physical masses to modal masses is illustrated in Figure [10](#org48d23d9) for the mode 2 of the example system.
 The modal stiffness \\(k\_2\\) is simply calculated via the relation between natural frequency, mass and stiffness:
 
 \begin{equation}
 k\_i = \omega\_i^2 m\_i
 \end{equation}
 
-<a id="org0de4b98"></a>
+<a id="org48d23d9"></a>
 
 {{< figure src="/ox-hugo/rankers98_mode_2_lumped_masses.png" caption="Figure 10: Graphical representation of mode 2 with (a.) lumped masses and (b.) modal mass and stiffness" >}}
 
 Let's now consider the effect of excitation forces that act on the physical DoF.
 The scalar product \\(\phi\_{ik}f\_k\\) of each force component with the corresponding element of the mode shape vector can be seen as the moment that acts on the modal level, or as an equivalent force that acts at the location of \\(q\_i\\) on the lever.
-Based on the graphical representation in Figure [11](#orgef7eb39), it is not difficult to understand the contribution of mode i to the transfer function \\(x\_l/f\_k\\):
+Based on the graphical representation in Figure [11](#org9f1d26b), it is not difficult to understand the contribution of mode i to the transfer function \\(x\_l/f\_k\\):
 
 \begin{equation}
 \left( \frac{x\_l}{f\_k} \right)\_i = \frac{\phi\_{ik}\phi\_{il}}{m\_i s^2 + k\_i}
@@ -511,7 +511,7 @@ Based on the graphical representation in Figure [11](#orgef7eb39), it is not dif
 
 Hence, the force \\(f\_k\\) must be multiplied by the distance \\(\phi\_{ik}\\) in order to find the equivalent excitation force at the location of \\(q\_i\\) on the lever, whereas the resulting modal displacement \\(q\_i\\) must be multiplied by the distance \\(\phi\_{il}\\) in order to obtain the displacement of the physical DoF \\(x\_l\\).
 
-<a id="orgef7eb39"></a>
+<a id="org9f1d26b"></a>
 
 {{< figure src="/ox-hugo/rankers98_lever_representation_with_force.png" caption="Figure 11: Graphical representation of mode \\(i\\), including the proper location of a force component \\(f\_k\\) that acts on physical DoF \\(x\_k\\)" >}}
 
@@ -522,14 +522,14 @@ This linear combination of physical DoF, which will be called "User DoF" can be
 x\_u = b\_1 x\_1 + \dots + b\_n x\_n = b^T x
 \end{equation}
 
-User DoF can be indicated on the modal lever, as illustrated in Figure [12](#orga82e374) for a user DoF \\(x\_u = x\_3 - x\_2\\).
+User DoF can be indicated on the modal lever, as illustrated in Figure [12](#orge1022d0) for a user DoF \\(x\_u = x\_3 - x\_2\\).
 The location of this user DoF \\(x\_u\\) with respect to the pivot point of modal lever \\(i\\) is defined by \\(\phi\_{iu}\\):
 
 \begin{equation}
 \phi\_{iu} = b^T \phi\_i
 \end{equation}
 
-<a id="orga82e374"></a>
+<a id="orge1022d0"></a>
 
 {{< figure src="/ox-hugo/rankers98_representation_user_dof.png" caption="Figure 12: Graphical representation of mode including user DoF \\(x\_u = x\_3 - x\_2\\)" >}}
 
@@ -539,13 +539,13 @@ Even though the dimension mode vector can be very large, only three user DoF are
 -   measured output (displacement that is measured by the position sensor)
 -   actual output (displacement that determines the accuracy of the machine)
 
-To illustrate this, a servo controlled positioning device is shown in Figure [13](#orgae9e430).
+To illustrate this, a servo controlled positioning device is shown in Figure [13](#orgb19de84).
 The task of the device is to position the payload with respect to a tool that is mounted to the machine frame.
 The actual accuracy of the machine is determined by the relative motion of these two components (actual output).
 However, direct measurement of the distance between the tool and the payload is not possible and therefore the control action is based on the measured distance between a sensor and the slide on which the payload is mounted (measured output).
 The slide is driven by a linear motor which transforms the output of the controller into a force on the slide and a reaction force on the stator (input).
 
-<a id="orgae9e430"></a>
+<a id="orgb19de84"></a>
 
 {{< figure src="/ox-hugo/rankers98_servo_system.png" caption="Figure 13: Schematic representation of a servo system" >}}
 
@@ -568,17 +568,17 @@ These effective modal parameters can be used very effectively in understanding t
 <div class="exampl">
   <div></div>
 
-The eigenvalue analysis of the two mass spring system in Figure [14](#orgbf0f71b) leads to the modal results summarized in Table [1](#table--tab:2dof-example-modal-params) and which are graphically represented in Figure [15](#org9a6d0a0).
+The eigenvalue analysis of the two mass spring system in Figure [14](#org8c59f86) leads to the modal results summarized in Table [1](#table--tab:2dof-example-modal-params) and which are graphically represented in Figure [15](#org3085b08).
 
-<a id="orgbf0f71b"></a>
+<a id="org8c59f86"></a>
 
 {{< figure src="/ox-hugo/rankers98_example_2dof.png" caption="Figure 14: Two mass spring system" >}}
 
 The modal masses can be easily found from the mode shape vectors:
 
 \begin{align}
-m\_1 &= \phi\_1^T M \phi\_1 = 50.8\,kg \\\\\\
+m\_1 &= \phi\_1^T M \phi\_1 = 50.8 kg \\\\\\
-m\_2 &= \phi\_2^T M \phi\_2 = 11.1\,kg
+m\_2 &= \phi\_2^T M \phi\_2 = 11.1 kg
 \end{align}
 
 whereas the modal stiffnesses follow from \\(k\_i = \omega\_i^2 m\_i\\).
@@ -596,7 +596,7 @@ whereas the modal stiffnesses follow from \\(k\_i = \omega\_i^2 m\_i\\).
 | Modal Mass [kg]   | \\(m\_1 = 50.8\\)             | \\(m\_2 = 11.1\\)               |
 | Modal Stiff [N/m] | \\(k\_1 = 0.46\cdot 10^7\\)   | \\(k\_2 = 1.23\cdot 10^7\\)     |
 
-<a id="org9a6d0a0"></a>
+<a id="org3085b08"></a>
 
 {{< figure src="/ox-hugo/rankers98_example_2dof_modal.png" caption="Figure 15: Graphical representation of modes and modal parameters of the two mass spring system" >}}
 
@@ -609,25 +609,25 @@ The results are summarized in Table [2](#table--tab:2dof-example-modal-params-ef
   Effective modal parameters for the two mass spring system
 </div>
 
-|                         | Mode 1                                          | Mode 2                                          |
+|                         | Mode 1                                         | Mode 2                                         |
-|-------------------------|-------------------------------------------------|-------------------------------------------------|
+|-------------------------|------------------------------------------------|------------------------------------------------|
-| Effective mass - DoF 1  | \\(m\_{\text{eff},11} = 112.1\,kg\\)            | \\(m\_{\text{eff},21} = 927.9\,kg\\)            |
+| Effective mass - DoF 1  | \\(m\_{\text{eff},11} = 112.1 kg\\)            | \\(m\_{\text{eff},21} = 927.9 kg\\)            |
-| Effective mass - DoF 2  | \\(m\_{\text{eff},12} = 92.8\,kg\\)             | \\(m\_{\text{eff},22} = 11.2\,kg\\)             |
+| Effective mass - DoF 2  | \\(m\_{\text{eff},12} = 92.8 kg\\)             | \\(m\_{\text{eff},22} = 11.2 kg\\)             |
-| Effective stiff - DoF 1 | \\(k\_{\text{eff},11} = 1.02 \cdot 10^7\,N/m\\) | \\(k\_{\text{eff},21} = 1.02 \cdot 10^9\,N/m\\) |
+| Effective stiff - DoF 1 | \\(k\_{\text{eff},11} = 1.02 \cdot 10^7 N/m\\) | \\(k\_{\text{eff},21} = 1.02 \cdot 10^9 N/m\\) |
-| Effective stiff - DoF 2 | \\(k\_{\text{eff},12} = 0.84 \cdot 10^7\,N/m\\) | \\(k\_{\text{eff},22} = 1.25 \cdot 10^7\,N/m\\) |
+| Effective stiff - DoF 2 | \\(k\_{\text{eff},12} = 0.84 \cdot 10^7 N/m\\) | \\(k\_{\text{eff},22} = 1.25 \cdot 10^7 N/m\\) |
 
-The effective modal parameters can then be used in the graphical representation of Figure [16](#orgad02a5c).
+The effective modal parameters can then be used in the graphical representation of Figure [16](#org1398e20).
-Based on this representation, it is now very easy to construct the individual modal contributions to the frequency response function \\(x\_1/F\_1\\) of the example system (Figure [17](#org5c3b059)).
+Based on this representation, it is now very easy to construct the individual modal contributions to the frequency response function \\(x\_1/F\_1\\) of the example system (Figure [17](#orgb620f7d)).
 
-<a id="orgad02a5c"></a>
+<a id="org1398e20"></a>
 
 {{< figure src="/ox-hugo/rankers98_example_2dof_effective_modal.png" caption="Figure 16: Alternative graphical representation of modes of two mass spring system based on the effective modal mass and stiffnesses in DoF \\(x\_1\\)" >}}
 
 One can observe that the low frequency part of each modal contribution corresponds to the inverse of the calculated effective modal mass stiffness at DoF \\(x\_1\\) whereas the high frequency contribution is defined by the effective modal mass.
 
-In the final Bode diagram (Figure [17](#org5c3b059), below) one can observe an interference of the two modal contributions in the frequency range of the second natural frequency, which in this example leads to a combination of an anti-resonance an a resonance.
+In the final Bode diagram (Figure [17](#orgb620f7d), below) one can observe an interference of the two modal contributions in the frequency range of the second natural frequency, which in this example leads to a combination of an anti-resonance an a resonance.
 
-<a id="org5c3b059"></a>
+<a id="orgb620f7d"></a>
 
 {{< figure src="/ox-hugo/rankers98_2dof_example_frf.png" caption="Figure 17: Frequency Response Function \\(x\_1/f\_1\\)" >}}
 
@@ -649,10 +649,10 @@ The technique furthermore gives an indication of the amount of frequency shift t
 <div class="exampl">
   <div></div>
 
-Assuming that one is asked to increase the natural frequency of the mode corresponding to Figure [18](#org78df8fe) by attaching a linear spring \\(k\\) between two of the three represented DoF.
+Assuming that one is asked to increase the natural frequency of the mode corresponding to Figure [18](#org015ebc2) by attaching a linear spring \\(k\\) between two of the three represented DoF.
 As the relative motion between \\(x\_A\\) and \\(x\_B\\) is the largest of all possible combinations, this is the choice that will maximize the natural frequency of the mode.
 
-<a id="org78df8fe"></a>
+<a id="org015ebc2"></a>
 
 {{< figure src="/ox-hugo/rankers98_example_3dof_sensitivity.png" caption="Figure 18: Graphical representation of a mod with 3 DoF" >}}
 
@@ -673,17 +673,17 @@ f\_{\text{new},i}(\Delta k) &= \frac{1}{2\pi}\sqrt{\frac{k\_{\text{eff},i} + \De
 <div class="exampl">
   <div></div>
 
-Let's use the two mass spring system in Figure [14](#orgbf0f71b) as an example.
+Let's use the two mass spring system in Figure [14](#org8c59f86) as an example.
 
 In order to analyze the effect of an extra mass at \\(x\_2\\), the effective modal mass at that DoF needs to be known for both modes (see Table [2](#table--tab:2dof-example-modal-params-eff)).
-Then using equation \eqref{eq:sensitivity_add_m}, one can estimate the effect of an extra mass \\(\Delta m = 1\,kg\\) added to \\(m\_2\\).
+Then using equation \eqref{eq:sensitivity_add_m}, one can estimate the effect of an extra mass \\(\Delta m = 1 kg\\) added to \\(m\_2\\).
 
 To estimate the influence of extra stiffness between the two DoF, one needs to calculate the effective modal stiffness that corresponds to the relative motion between \\(x\_2\\) and \\(x\_1\\).
-This can be graphically done as shown in Figure [19](#org4893595):
+This can be graphically done as shown in Figure [19](#orgf24a43b):
 
 \begin{align}
-k\_{\text{eff},1,(2-1)} &= 0.46 \cdot 10^7 / 0.07^2 = 93.9 \cdot 10^7 \, N/m \\\\\\
+k\_{\text{eff},1,(2-1)} &= 0.46 \cdot 10^7 / 0.07^2 = 93.9 \cdot 10^7   N/m \\\\\\
-k\_{\text{eff},2,(2-1)} &= 1.23 \cdot 10^7 / 1.1^2 = 1.0 \cdot 10^7 \, N/m
+k\_{\text{eff},2,(2-1)} &= 1.23 \cdot 10^7 / 1.1^2 = 1.0 \cdot 10^7   N/m
 \end{align}
 
 And using equation \eqref{eq:sensitivity_add_m}, the effect of additional stiffness on the frequency of the two modes can be computed.
@@ -696,13 +696,13 @@ The results are summarized in Table [3](#table--tab:example-sensitivity-2dof-res
   Sensitivity analysis results
 </div>
 
-|                                                                     | f1 [Hz] | f2 [Hz] |
+|                                                                    | f1 [Hz] | f2 [Hz] |
-|---------------------------------------------------------------------|---------|---------|
+|--------------------------------------------------------------------|---------|---------|
-| Original                                                            | 47.8    | 167.7   |
+| Original                                                           | 47.8    | 167.7   |
-| \\(\Delta m = 1\,kg\\) added to \\(m\_2\\)                          | 47.5    | 160.7   |
+| \\(\Delta m = 1 kg\\) added to \\(m\_2\\)                          | 47.5    | 160.7   |
-| \\(\Delta k = 10^7\, N/m\\) added between \\(x\_2\\) and \\(x\_1\\) | 48.1    | 237.2   |
+| \\(\Delta k = 10^7  N/m\\) added between \\(x\_2\\) and \\(x\_1\\) | 48.1    | 237.2   |
 
-<a id="org4893595"></a>
+<a id="orgf24a43b"></a>
 
 {{< figure src="/ox-hugo/rankers98_example_sensitivity_2dof.png" caption="Figure 19: Graphical representation of modes and modal parameters of two mass spring system" >}}
 
@@ -712,9 +712,9 @@ The results are summarized in Table [3](#table--tab:example-sensitivity-2dof-res
 ### Modal Superposition {#modal-superposition}
 
 Previously, the lever representation was used only to represent the individual mode shapes.
-In the mechanism shown in Figure [20](#org3df11f4), the motion of the output \\(y\\) is equals to the sum of the motion of the two inputs \\(x\_1\\) and \\(x\_2\\).
+In the mechanism shown in Figure [20](#orgc27801a), the motion of the output \\(y\\) is equals to the sum of the motion of the two inputs \\(x\_1\\) and \\(x\_2\\).
 
-<a id="org3df11f4"></a>
+<a id="orgc27801a"></a>
 
 {{< figure src="/ox-hugo/rankers98_addition_of_motion.png" caption="Figure 20: Addition of motion" >}}
 
@@ -724,9 +724,9 @@ This approach can be applied to the concept of modal superposition, which expres
 x\_k(t) = \sum\_{i=1}^n \phi\_{ik} q\_i(t) = \sum\_{i=1}^n x\_{ki}(t)
 \end{equation}
 
-Combining the concept of summation of modal contribution with the lever representation of mode shapes leads to Figure [21](#orgb8f94d2), which is a visualization of the transformation between the modal and the physical domains.
+Combining the concept of summation of modal contribution with the lever representation of mode shapes leads to Figure [21](#org67908e9), which is a visualization of the transformation between the modal and the physical domains.
 
-<a id="orgb8f94d2"></a>
+<a id="org67908e9"></a>
 
 {{< figure src="/ox-hugo/rankers98_conversion_modal_to_physical.png" caption="Figure 21: Conversion between modal DoF to physical DoF" >}}
 
@@ -737,9 +737,9 @@ The "rigid body modes" usually refer to the lower natural frequencies of a machi
 This is misleading at it suggests that the structure exhibits no internal deformation.
 A better term for such a mode would be **suspension mode**.
 
-To illustrate the important of the internal deformation, a very simplified physical model of a precision machine is considered (Figure [22](#orgc5ca877)).
+To illustrate the important of the internal deformation, a very simplified physical model of a precision machine is considered (Figure [22](#org44924f2)).
 
-<a id="orgc5ca877"></a>
+<a id="org44924f2"></a>
 
 {{< figure src="/ox-hugo/rankers98_suspension_mode_machine.png" caption="Figure 22: Simplified physical model of a precision machine" >}}
 
@@ -762,11 +762,11 @@ This equation shows that the internal deformation associated with the suspension
 <div class="exampl">
   <div></div>
 
-As an example of a situation in which the internal deformation associated with the suspension mode is of significant importance, one could consider a high precision machine that is excited due to floor vibrations such that it vibrates on its suspension with an amplitude of \\(100\, \mu m\\) and a frequency of 3 Hz.
+As an example of a situation in which the internal deformation associated with the suspension mode is of significant importance, one could consider a high precision machine that is excited due to floor vibrations such that it vibrates on its suspension with an amplitude of \\(100  \mu m\\) and a frequency of 3 Hz.
 Assuming that the internal frequency of the manipulator is equal to 150 Hz, the internal deformation of the machine is:
 
 \begin{equation}
-\epsilon = \frac{3^2}{150^2} 100 \, \mu m = 40\, nm
+\epsilon = \frac{3^2}{150^2} 100   \mu m = 40  nm
 \end{equation}
 
 which can be a lot for high precision machines.
@@ -776,14 +776,16 @@ which can be a lot for high precision machines.
 
 ## Modes and Servo Stability {#modes-and-servo-stability}
 
-One of the two limiting effects
+The effect of machine dynamics on the servo control loop stability is discussed in this chapter.
+
+The interaction between the desired (rigid body) motion and the dynamics of one additional mode and its effect on the freuency response function \\(x\_{\text{servo}}/F\_{\text{servo}}\\) is the basis of this chapter.
 
 
 ### Basic Characteristics of Mechanical FRF {#basic-characteristics-of-mechanical-frf}
 
-Consider the position control loop of Figure [23](#org18af44d).
+Consider the position control loop of Figure [23](#orgcd512cf).
 
-<a id="org18af44d"></a>
+<a id="orgcd512cf"></a>
 
 {{< figure src="/ox-hugo/rankers98_mechanical_servo_system.png" caption="Figure 23: Mechanical position servo-system" >}}
 
@@ -793,7 +795,9 @@ In the ideal situation the mechanical system behaves as one rigid body with mass
 \frac{x\_{\text{servo}}}{F\_{\text{servo}}} = \frac{1}{m s^2}
 \end{equation}
 
-<a id="orgac4be89"></a>
+The corresponding Bode and Nyquist plots and shown in Figure [24](#orgb3ba1d2).
+
+<a id="orgb3ba1d2"></a>
 
 {{< figure src="/ox-hugo/rankers98_ideal_bode_nyquist.png" caption="Figure 24: FRF of an ideal system with no resonances" >}}
 
@@ -813,45 +817,184 @@ Let's introduce a variable \\(\alpha\\), which relates the high-frequency contri
 
 which simplifies equation \eqref{eq:effect_one_mode} to:
 
-\begin{equation}
+\begin{equation} \label{eq:effect\_one\_mode\_simplified}
 \frac{x\_{\text{servo}}}{F\_{\text{servo}}} = \frac{1}{ms^2} + \frac{\alpha}{m s^2 + m \omega\_i^2}
 \end{equation}
 
-<a id="orgda4ef86"></a>
+Equation \eqref{eq:effect_one_mode_simplified} will be the basis for the discussion of the various patterns that can be observe in the frequency response functions and the effect of resonances on servo stability.
+
+Three different types of intersection pattern can be found in the amplitude plot as shown in Figure [25](#orgf995437).
+Depending on the absolute value of \\(\alpha\\) one can observe:
+
+-   \\(|\alpha| < 1\\): two intersections
+-   \\(|\alpha| = 1\\): one intersection and asymptotic approach at high frequencies
+-   \\(|\alpha| > 1\\): one intersection
+
+The interaction between the rigid body motion and the additional mode will not only depend on \\(|\alpha|\\) but also on the sign of \\(\alpha\\), which determined the phase relation between the two contributions.
+
+<a id="orgf995437"></a>
 
 {{< figure src="/ox-hugo/rankers98_frf_effect_alpha.png" caption="Figure 25: Contribution of rigid-body motion and modal dynamics to the amplitude and phase of FRF for various values of \\(\alpha\\)" >}}
 
-<a id="org6c0aa47"></a>
+The general shape of the overall FRF can be constructed for all cases (Figure [26](#org3659f0b)).
+Interesting points are the interaction of the two parts at the frequency that corresponds to an intersection in the amplitude plot.
+At this frequency the magnitudes are equal, so it depends on the phase of the two contributions whether they cancel each other, thus leading to a zero, or just add up.
+
+<a id="org3659f0b"></a>
 
 {{< figure src="/ox-hugo/rankers98_final_frf_alpha.png" caption="Figure 26: Bode diagram of final FRF (\\(x\_{\text{servo}}/F\_{\text{servo}}\\)) for six values of \\(\alpha\\)" >}}
 
+When analyzing the plots of Figure [26](#org3659f0b), four different types of FRF can be found:
 
-#### FRF with "-2 slope/zero/pole/-2 slope" behavior {#frf-with-2-slope-zero-pole-2-slope-behavior}
+-   -2 slope / zero / pole / -2 slope (\\(\alpha > 0\\))
+-   -2 slope / pole / zero / -2 slope (\\(-1 < \alpha < 0\\))
+-   -2 slope / pole / -4 slope (\\(\alpha = -1\\))
+-   -2 slope / pole / -2 slope (\\(\alpha < -1\\))
 
+All cases are shown in Figure [27](#org5f42df0).
 
-#### FRF with "-2 slope/pole/zero/-2 slope" behavior {#frf-with-2-slope-pole-zero-2-slope-behavior}
+<a id="org5f42df0"></a>
 
-
+{{< figure src="/ox-hugo/rankers98_interaction_shapes.png" caption="Figure 27: Bode plot of the different types of FRF" >}}
-#### 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}
-
-
-#### Summary {#summary}
-
-|           | << wb | ~ wb | >> wb |
-|-----------|-------|------|-------|
-| Pole Zero | OK    | NOK  | OK    |
-| Zero Pole | NOK   | OK   | NOK   |
-| Pole      | NOK   | NOK  | OK    |
 
 
 ### Destabilising Effect of Modes {#destabilising-effect-of-modes}
 
+In this section the effect of each basic mechanical FRF on the stability of the control loop will be discussed.
+A PID controller with additional second order low pass filter will be the basis for the discussion of stability.
+
+Typical crossover frequencies for a PID controller with second order low pass filtering are:
+
+\begin{align\*}
+f\_i &= f\_b/10 \\\\\\
+f\_d &= f\_b/3 \\\\\\
+f\_{lp} &= 4 \cdot f\_b
+\end{align\*}
+
+with \\(f\_b\\) the bandwidth frequency.
+The asymptotic amplitude plot is shown in Figure [28](#orgc6fa8a4).
+
+<a id="orgc6fa8a4"></a>
+
+{{< figure src="/ox-hugo/rankers98_pid_amplitude.png" caption="Figure 28: Typical crossover frequencies of a PID controller with 2nd order low pass filtering" >}}
+
+With these settings, the open loop response of the position loop (controller + mechanics) looks like Figure [29](#orgdf31690).
+
+<a id="orgdf31690"></a>
+
+{{< figure src="/ox-hugo/rankers98_ideal_frf_pid.png" caption="Figure 29: Ideal open loop FRF of a position servo without mechanical resonances (\\(f\_b = 30\text{ Hz}\\))" >}}
+
+Conclusions are:
+
+-   A "-2 slope / zero / pole / -2 slope" characteristic leads to a phase lead, and is therefore potentially destabilizing in the low frequency (Figure [30](#org9480ff1)) and high frequency (Figure [32](#org9b86fef)) regions.
+    In the medium frequency region (Figure [31](#orga4bae4e)), it adds an extra phase lead to the already existing margin, which does not harm the stability.
+-   A "-2 slope / pole / zero / -2 slope" combination has the reverse effect.
+    It is potentially destabilizing in the medium frequency range (Figure [34](#orga71718d)) and is harmless in the low (Figure [33](#org9b0e1e4)) and high frequency (Figure [35](#orgc0f181f)) ranges.
+-   The "-2 slope / pole / -4 slope" behavior always has a devastating effect on the stability of the loop if located in the low of medium frequency ranges.
+
+These conclusions may differ for different mass ratio \\(\alpha\\).
+
+<a id="org9480ff1"></a>
+
+{{< figure src="/ox-hugo/rankers98_zero_pole_low_freq.png" caption="Figure 30: Open Loop FRF of type \"-2 slope / zero / pole / -2 slope\" with low frequency resonance" >}}
+
+<a id="orga4bae4e"></a>
+
+{{< figure src="/ox-hugo/rankers98_zero_pole_medium_freq.png" caption="Figure 31: Open Loop FRF of type \"-2 slope / zero / pole / -2 slope\" with medium frequency resonance" >}}
+
+<a id="org9b86fef"></a>
+
+{{< figure src="/ox-hugo/rankers98_zero_pole_high_freq.png" caption="Figure 32: Open Loop FRF of type \"-2 slope / zero / pole / -2 slope\" with high frequency resonance" >}}
+
+<a id="org9b0e1e4"></a>
+
+{{< figure src="/ox-hugo/rankers98_pole_zero_low_freq.png" caption="Figure 33: Open Loop FRF of type \"-2 slope / pole / zero / -2 slope\" with low frequency resonance" >}}
+
+<a id="orga71718d"></a>
+
+{{< figure src="/ox-hugo/rankers98_pole_zero_medium_freq.png" caption="Figure 34: Open Loop FRF of type \"-2 slope / pole / zero / -2 slope\" with medium frequency resonance" >}}
+
+<a id="orgc0f181f"></a>
+
+{{< figure src="/ox-hugo/rankers98_pole_zero_high_freq.png" caption="Figure 35: Open Loop FRF of type \"-2 slope / pole / zero / -2 slope\" with high frequency resonance" >}}
+
 
 ### Design for Stability {#design-for-stability}
 
+
+#### Actuator Flexibility {#actuator-flexibility}
+
+Figure [36](#org5c9d040) shows the schematic representation of a system with a certain compliance between the motor and the load.
+
+<a id="org5c9d040"></a>
+
+{{< figure src="/ox-hugo/rankers98_2dof_actuator_flexibility.png" caption="Figure 36: Servo system with actuator flexibility - Schematic representation" >}}
+
+The corresponding modes are shown in Figure [37](#orgc92a13c).
+
+<a id="orgc92a13c"></a>
+
+{{< figure src="/ox-hugo/rankers98_2dof_modes_act_flex.png" caption="Figure 37: Servo System with Actuator Flexibility - Modes" >}}
+
+Assuming first that the **servo position is measured at the motor**.
+The following transfer function must be considered:
+
+\begin{align}
+\frac{x\_1}{F\_{\text{servo}}} &= \frac{1}{m\_{\text{eff},11} s^2} + \frac{1}{m\_{\text{eff,21}}s^2 + \omega\_2^2 m\_{\text{eff},21}} \\\\\\
+&= \frac{1}{m\_1 + m\_2} \left( \frac{1}{s^2} + \frac{\alpha}{s^2 + \omega\_2^2} \right)
+\end{align}
+
+with \\(\alpha = m\_2/m\_1\\) (mass ratio) relates the "mass" of the additional modal contribution to the mass of the rigid body motion.
+The resulting FRF exhibit a "-2 slope / zero / pole / -2 slope" (Figure [38](#orgb7822b6)).
+
+<a id="orgb7822b6"></a>
+
+{{< figure src="/ox-hugo/rankers98_2dof_act_flex_frf.png" caption="Figure 38: Mechanical FRF of a system with actuator flexibility and position measurement at motor" >}}
+
+The asymptotes at low and high frequencies are:
+
+\begin{align}
+\left( \frac{x\_1}{F\_{\text{servo}}} \right)\_{s \to 0} &= \frac{1}{(m\_1 + m\_2) s^2} \\\\\\
+\left( \frac{x\_1}{F\_{\text{servo}}} \right)\_{s \to 0} &= \frac{1}{(m\_1 + m\_2) s^2} + \frac{1}{m\_1/m\_2(m\_1 + m\_2) s^2} = \frac{1}{m\_1 s^2}
+\end{align}
+
+which corresponds to the engineering feeling that at very low frequencies the two masses move as one single mass, whereas at very high frequencies the mass \\(m\_2\\) of the load is completely decoupled such that the servo system only "feels and sees" the motion of the motor mass \\(m\_1\\).
+
+<div class="important">
+  <div></div>
+
+Guideline in presence of actuator flexibility with measurement at the motor position:
+
+-   The motor inertia should be one to three times to inertial of the load
+-   The resonance frequency should either be near the bandwidth frequency and much above
+
+</div>
+
+Now assume that the **servo position is measured at the load**.
+Now we are interested by the following transfer function:
+
+\begin{equation}
+\frac{x\_2}{F\_{\text{servo}}} = = \frac{1}{m\_1 + m\_2} \left( \frac{1}{s^2} - \frac{1}{s^2 + \omega\_2^2} \right)
+\end{equation}
+
+The mass ratio \\(\alpha\\) equal -1, and thus the FRF will be of type "-2 slope / pole / -4 slope" (Figure [39](#orgbe8ffaf)).
+
+<a id="orgbe8ffaf"></a>
+
+{{< figure src="/ox-hugo/rankers98_2dof_act_flex_meas_load_frf.png" caption="Figure 39: FRF \\(k\_p \cdot (x\_{\text{servo}}/F\_{\text{servo}})\\) of a system with actuator flexibility and position measurement at the load" >}}
+
+<div class="important">
+  <div></div>
+
+Guideline in presence of actuator flexibility with measurement at the load position:
+
+-   Resonance frequency larger than 5 to 10 times the wanted bandwidth
+
+</div>
+
+
+#### Guiding System Flexibility {#guiding-system-flexibility}
+
 <div class="important">
   <div></div>
 
@@ -865,6 +1008,12 @@ Guiding system flexibility:
 </div>
 
 
+#### Limited Mass and Stiffness of Stationary Machine Part {#limited-mass-and-stiffness-of-stationary-machine-part}
+
+
+#### General Guidelines {#general-guidelines}
+
+
 ## Predictive Modelling {#predictive-modelling}
 
 
@@ -907,4 +1056,4 @@ Through the enormous performance drive in mechatronics systems, much has been le
 
 ## Bibliography {#bibliography}
 
-<a id="org6436602"></a>Rankers, Adrian Mathias. 1998. “Machine Dynamics in Mechatronic Systems: An Engineering Approach.” University of Twente.
+<a id="orgf6a233a"></a>Rankers, Adrian Mathias. 1998. “Machine Dynamics in Mechatronic Systems: An Engineering Approach.” University of Twente.