diff --git a/figs/3d-cubic-stewart-aligned.png b/figs/3d-cubic-stewart-aligned.png new file mode 100644 index 0000000..f8d4d7b Binary files /dev/null and b/figs/3d-cubic-stewart-aligned.png differ diff --git a/figs/dvf_root_locus_ws.png b/figs/dvf_root_locus_ws.png new file mode 100644 index 0000000..8d53867 Binary files /dev/null and b/figs/dvf_root_locus_ws.png differ diff --git a/figs/iff_root_locus_ws.png b/figs/iff_root_locus_ws.png new file mode 100644 index 0000000..bdc784c Binary files /dev/null and b/figs/iff_root_locus_ws.png differ diff --git a/figs/mobility_translations_null_rotation.png b/figs/mobility_translations_null_rotation.png new file mode 100644 index 0000000..34b10ef Binary files /dev/null and b/figs/mobility_translations_null_rotation.png differ diff --git a/figs/preumont07_flexible_joints.png b/figs/preumont07_flexible_joints.png new file mode 100644 index 0000000..9cd54b2 Binary files /dev/null and b/figs/preumont07_flexible_joints.png differ diff --git a/figs/root_locus_dvf_rot_stiffness.png b/figs/root_locus_dvf_rot_stiffness.png new file mode 100644 index 0000000..2d5c6f2 Binary files /dev/null and b/figs/root_locus_dvf_rot_stiffness.png differ diff --git a/figs/root_locus_iff_rot_stiffness.png b/figs/root_locus_iff_rot_stiffness.png new file mode 100644 index 0000000..15140e8 Binary files /dev/null and b/figs/root_locus_iff_rot_stiffness.png differ diff --git a/figs/yang19_flexible_joints.png b/figs/yang19_flexible_joints.png new file mode 100644 index 0000000..5a1467f Binary files /dev/null and b/figs/yang19_flexible_joints.png differ diff --git a/index.html b/index.html index 28ead57..d9035a1 100644 --- a/index.html +++ b/index.html @@ -4,7 +4,7 @@ "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"> - + Design of the Nano-Hexapod and associated Control Architectures - Summary @@ -35,94 +35,119 @@

Table of Contents

-
-

1.1.1 Simplified Feedback Control Diagram for the NASS

+
+

1.1.1 Simplified Feedback Control Diagram for the NASS

-Let’s consider the block diagram shown in Figure 1 where the signals are: +Let’s consider the block diagram shown in Figure 1 where the signals are:

  • \(y\): the relative position of the sample with respect to the granite (the quantity to be controlled)
    • @@ -248,7 +273,7 @@ The dynamical blocks are:
-
+

classical_feedback_small.png

Figure 1: Block Diagram of a simple feedback system

@@ -270,11 +295,11 @@ In the next section, we see how the use of the feedback system permits to lower
-
-

1.1.2 How does the feedback loop is modifying the system behavior?

+
+

1.1.2 How does the feedback loop is modifying the system behavior?

-If we write down the position error signal \(\epsilon = r - y\) as a function of the reference signal \(r\), the disturbances \(d\) and the measurement noise \(n\) (using the feedback diagram in Figure 1), we obtain: +If we write down the position error signal \(\epsilon = r - y\) as a function of the reference signal \(r\), the disturbances \(d\) and the measurement noise \(n\) (using the feedback diagram in Figure 1), we obtain: \[ \epsilon = \frac{1}{1 + GK} r + \frac{GK}{1 + GK} n - \frac{G_d}{1 + GK} d \]

@@ -315,8 +340,8 @@ Ideally, we would like to design the controller \(K\) such that:
-
-

1.1.3 Trade off: Disturbance Reduction / Noise Injection

+
+

1.1.3 Trade off: Disturbance Reduction / Noise Injection

We have from the definition of \(S\) and \(T\) that: @@ -334,7 +359,7 @@ There is therefore a trade-off between the disturbance rejection and the meas

-Typical shapes of \(|S|\) and \(|T|\) as a function of frequency are shown in Figure 2. +Typical shapes of \(|S|\) and \(|T|\) as a function of frequency are shown in Figure 2. We can observe that \(|S|\) and \(|T|\) exhibit different behaviors depending on the frequency band:

@@ -360,7 +385,7 @@ We can observe that \(|S|\) and \(|T|\) exhibit different behaviors depending on -
+

h-infinity-2-blocs-constrains.png

Figure 2: Typical shapes and constrain of the Sensibility and Transmibility closed-loop transfer functions

@@ -368,11 +393,11 @@ We can observe that \(|S|\) and \(|T|\) exhibit different behaviors depending on
-
-

1.1.4 Trade off: Robustness / Performance

+
+

1.1.4 Trade off: Robustness / Performance

- +

@@ -393,11 +418,11 @@ The main issue it that for stability reasons, the system dynamics must be kno

-For mechanical systems, this generally means that the control bandwidth should take place before any appearing of flexible dynamics (right part of Figure 3). +For mechanical systems, this generally means that the control bandwidth should take place before any appearing of flexible dynamics (right part of Figure 3).

-
+

oomen18_next_gen_loop_gain.png

Figure 3: Envisaged developments in motion systems. In traditional motion systems, the control bandwidth takes place in the rigid-body region. In the next generation systemes, flexible dynamics are foreseen to occur within the control bandwidth. oomen18_advan_motion_contr_precis_mechat

@@ -427,11 +452,11 @@ This problem of robustness represent one of the main challenge for the de
-
-

1.2 Dynamic error budgeting

+
+

1.2 Dynamic error budgeting

- +

The dynamic error budgeting is a powerful tool to study the effects of multiple error sources (i.e. disturbances and measurement noise) and to predict how much these effects are reduced by a feedback system. @@ -442,19 +467,19 @@ The dynamic error budgeting uses two important mathematical functions: the Po

-After these two functions are introduced (in Sections 1.2.1 and 1.2.2), is shown how do multiple error sources are combined and modified by dynamical systems (in Section 1.2.3 and 1.2.4). +After these two functions are introduced (in Sections 1.2.1 and 1.2.2), is shown how do multiple error sources are combined and modified by dynamical systems (in Section 1.2.3 and 1.2.4).

-Finally, the dynamic noise budgeting for the NASS is derived in Section 1.2.5. +Finally, the dynamic noise budgeting for the NASS is derived in Section 1.2.5.

-
-

1.2.1 Power Spectral Density

+
+

1.2.1 Power Spectral Density

- +

@@ -483,11 +508,11 @@ One can also integrate the infinitesimal power \(S_{xx}(\omega)d\omega\) over a

-
-

1.2.2 Cumulative Power Spectrum

+
+

1.2.2 Cumulative Power Spectrum

- +

@@ -521,19 +546,19 @@ The Cumulative Power Spectrum is generally shown as a function of frequency, and

-
-

1.2.3 Modification of a signal’s PSD when going through a dynamical system

+
+

1.2.3 Modification of a signal’s PSD when going through a dynamical system

- +

-Let’s consider a signal \(u\) with a PSD \(S_{uu}\) going through a LTI system \(G(s)\) that outputs a signal \(y\) with a PSD (Figure 4). +Let’s consider a signal \(u\) with a PSD \(S_{uu}\) going through a LTI system \(G(s)\) that outputs a signal \(y\) with a PSD (Figure 4).

-
+

psd_lti_system.png

Figure 4: LTI dynamical system \(G(s)\) with input signal \(u\) and output signal \(y\)

@@ -548,15 +573,15 @@ The Power Spectral Density of the output signal \(y\) can be computed using:
-
-

1.2.4 PSD of combined signals

+
+

1.2.4 PSD of combined signals

- +

-Let’s consider a signal \(y\) that is the sum of two uncorrelated signals \(u\) and \(v\) (Figure 5). +Let’s consider a signal \(y\) that is the sum of two uncorrelated signals \(u\) and \(v\) (Figure 5).

@@ -565,7 +590,7 @@ We have that the PSD of \(y\) is equal to sum of the PSD and \(u\) and the PSD o

-
+

psd_sum.png

Figure 5: \(y\) as the sum of two signals \(u\) and \(v\)

@@ -573,15 +598,15 @@ We have that the PSD of \(y\) is equal to sum of the PSD and \(u\) and the PSD o
-
-

1.2.5 Dynamic Noise Budgeting

+
+

1.2.5 Dynamic Noise Budgeting

- +

-Let’s consider the Feedback architecture in Figure 1 where the position error \(\epsilon\) is equal to: +Let’s consider the Feedback architecture in Figure 1 where the position error \(\epsilon\) is equal to: \[ \epsilon = S r + T n - G_d S d \]

@@ -605,25 +630,25 @@ To estimate the PSD of the position error \(\epsilon\) and thus the RMS residual
  • The Power Spectral Densities of the signals affecting the system:
      -
    • The disturbances \(S_{dd}\): this will be done in Section 3
      • +
    • The disturbances \(S_{dd}\): this will be done in Section 3
    • The sensor noise \(S_{nn}\): this can be estimated from the sensor data-sheet
    • The wanted sample’s motion \(S_{rr}\): this is a deterministic signal that we choose. For a simple tomography experiment, we can consider that it is equal to \(0\) as we only want to compensate all the sample’s vibrations
  • The dynamics of the complete system comprising the micro-station and the nano-hexapod: \(G\), \(G_d\). -To do so, we need to identify the dynamics of the micro-station (Section 2), include this dynamics in a model (Section 4) and add a model of the nano-hexapod to the model (Section 5)
    • -
  • The controller \(K\) that will be designed in Section 6
    • +To do so, we need to identify the dynamics of the micro-station (Section 2), include this dynamics in a model (Section 4) and add a model of the nano-hexapod to the model (Section 5) +
  • The controller \(K\) that will be designed in Section 6
-
-

2 Identification of the Micro-Station Dynamics

+
+

2 Identification of the Micro-Station Dynamics

- +

As explained before, it is very important to have a good estimation of the micro-station dynamics as it will be used: @@ -641,7 +666,7 @@ All the measurements performed on the micro-station are detailed in 6. +The general procedure to identify the dynamics of the micro-station is shown in Figure 6. The steps are:

    @@ -651,7 +676,7 @@ The steps are:
-
+

vibration_analysis_procedure.png

Figure 6: Vibration Analysis Procedure

@@ -663,11 +688,11 @@ Instead, the model will be tuned using both the modal model and the response mod

-
-

2.1 Experimental Setup

+
+

2.1 Experimental Setup

- +

@@ -693,13 +718,13 @@ In order to perform the modal analysis, the following devices were used: The measurement consists of:

    -
  • Exciting the structure at the same location with the instrumented hammer (Figure 7)
    • +
  • Exciting the structure at the same location with the instrumented hammer (Figure 7)
  • Fix the accelerometers on each of the stages to measure all the DOF of the structure. The position of the accelerometers are:
    • 4 on the first granite
    • 4 on the second granite
      • -
    • 4 on top of the translation stage (Figure 8)
      • +
    • 4 on top of the translation stage (Figure 8)
    • 4 on top of the tilt stage
    • 3 on top of the spindle
    • 4 on top of the hexapod
      • @@ -715,14 +740,14 @@ It was chosen to have some redundancy in the measurement to be able to verify th

      -
      +

      hammer_z.gif

      Figure 7: Example of one hammer impact

      -
      +

      accelerometers_ty_overview.jpg

      Figure 8: 3 tri axis accelerometers fixed to the translation stage

      @@ -730,11 +755,11 @@ It was chosen to have some redundancy in the measurement to be able to verify th
      -
      -

      2.2 Results

      +
      +

      2.2 Results

      - +

      @@ -743,18 +768,18 @@ From the measurements are extracted all the transfer functions from forces appli

      Modal shapes and natural frequencies are then computed. -Example of the obtained micro-station’s mode shapes are shown in Figures 9 and 10. +Example of the obtained micro-station’s mode shapes are shown in Figures 9 and 10.

      -
      +

      mode1.gif

      Figure 9: First mode that shows a suspension mode, probably due to bad leveling of one Airloc

      -
      +

      mode6.gif

      Figure 10: Sixth mode

      @@ -782,12 +807,12 @@ This thus means that a multi-body model can be used to correctly represent th

      Many Frequency Response Functions (FRF) are obtained from the measurements. -Examples of FRF are shown in Figure 11. +Examples of FRF are shown in Figure 11. These FRF will be used to compare the dynamics of the multi-body model with the micro-station dynamics.

      -
      +

      frf_all_bodies_one_direction.png

      Figure 11: Frequency Response Function from forces applied by the Hammer in the X direction to the acceleration of each solid body in the X direction

      @@ -795,8 +820,8 @@ These FRF will be used to compare the dynamics of the multi-body model with the
      -
      -

      2.3 Conclusion

      +
      +

      2.3 Conclusion

      @@ -804,7 +829,7 @@ The dynamical measurements made on the micro-station confirmed the fact that a m

      -In Section 4, the obtained Frequency Response Functions will be used to compare the model dynamics with the micro-station dynamics. +In Section 4, the obtained Frequency Response Functions will be used to compare the model dynamics with the micro-station dynamics.

      @@ -812,11 +837,11 @@ In Section 4, the obtained Frequency Response Function
      -
      -

      3 Identification of the Disturbances

      +
      +

      3 Identification of the Disturbances

      - +

      In this section, all the disturbances affecting the system are identified and quantified. @@ -830,13 +855,13 @@ Note that the low frequency disturbances such as static guiding errors and therm The main challenge is to reduce the disturbances containing high frequencies, and thus efforts are made to identify these high frequency disturbances such as:

        -
      • Ground motion (Section 3.1)
        • -
      • Vibration introduced by control systems (Section 3.2)
        • -
      • Vibration introduced by the motion of the spindle and of the translation stage (Section 3.3)
        • +
      • Ground motion (Section 3.1)
        • +
      • Vibration introduced by control systems (Section 3.2)
        • +
      • Vibration introduced by the motion of the spindle and of the translation stage (Section 3.3)

      -A noise budgeting is performed in Section 3.4, the vibrations induced by the disturbances are compared and the required control bandwidth is estimated. +A noise budgeting is performed in Section 3.4, the vibrations induced by the disturbances are compared and the required control bandwidth is estimated.

      @@ -844,11 +869,11 @@ The measurements are presented in more detail in -

      3.1 Ground Motion

      +
      +

      3.1 Ground Motion

      - +

      @@ -856,12 +881,12 @@ Ground motion can easily be estimated using an inertial sensor with sufficient s

      -To verify that the inertial sensors are sensitive enough, a Huddle test has been performed (Figure 12). +To verify that the inertial sensors are sensitive enough, a Huddle test has been performed (Figure 12). The details of the Huddle Test can be found here.

      -
      +

      geophones.jpg

      Figure 12: Huddle Test Setup

      @@ -873,7 +898,7 @@ The low frequency differences between the ground motion at ID31 and ID09 is just

      -
      +

      ground_motion_compare.png

      Figure 13: Comparison of the PSD of the ground motion measured at different location

      @@ -881,11 +906,11 @@ The low frequency differences between the ground motion at ID31 and ID09 is just
      -
      -

      3.2 Stage Vibration - Effect of Control systems

      +
      +

      3.2 Stage Vibration - Effect of Control systems

      - +

      @@ -908,11 +933,11 @@ Complete reports on these measurements are accessible -

      3.3 Stage Vibration - Effect of Motion

      +
      +

      3.3 Stage Vibration - Effect of Motion

      - +

      We consider here the vibrations induced by scans of the translation stage and rotation of the spindle. @@ -923,15 +948,15 @@ Details reports are accessible -

      Spindle and Slip-Ring

      -
      +
      +

      Spindle and Slip-Ring

      +

      -The setup for the measurement of vibrations induced by rotation of the Spindle and Slip-ring is shown in Figure 14. +The setup for the measurement of vibrations induced by rotation of the Spindle and Slip-ring is shown in Figure 14.

      -
      +

      rz_meas_errors.gif

      Figure 14: Measurement of the sample’s vertical motion when rotating at 6rpm

      @@ -947,7 +972,7 @@ A geophone is fixed at the location of the sample and the motion is measured:

    -The obtained Power Spectral Densities of the sample’s absolute velocity are shown in Figure 15. +The obtained Power Spectral Densities of the sample’s absolute velocity are shown in Figure 15.

    @@ -964,7 +989,7 @@ Its cause has not been identified yet

-
+

sr_sp_psd_sample_compare.png

Figure 15: Comparison of the ASD of the measured voltage from the Geophone at the sample location

@@ -979,19 +1004,19 @@ Some investigation should be performed to determine where does this 23Hz motion
-
-

Translation Stage

-
+
+

Translation Stage

+

The same setup is used: a geophone is located at the sample’s location and another on the granite.

-A 1Hz triangle motion with an amplitude of \(\pm 2.5mm\) is sent to the translation stage (Figure 16), and the absolute velocities of the sample and the granite are measured. +A 1Hz triangle motion with an amplitude of \(\pm 2.5mm\) is sent to the translation stage (Figure 16), and the absolute velocities of the sample and the granite are measured.

-
+

ty_position_time.png

Figure 16: Y position of the translation stage measured by the encoders

@@ -999,26 +1024,27 @@ A 1Hz triangle motion with an amplitude of \(\pm 2.5mm\) is sent to the translat

-The time domain absolute vertical velocity of the sample and granite are shown in Figure 17. +The time domain absolute vertical velocity of the sample and granite are shown in Figure 17. It is shown that quite large motion of the granite is induced by the translation stage scans. This could be a problem if this is shown to excite the metrology frame of the nano-focusing lens position stage.

-
+

ty_z_time.png

Figure 17: Vertical velocity of the sample and marble when scanning with the translation stage

-The Amplitude Spectral Densities of the measured absolute velocities are shown in Figure 18. +The Amplitude Spectral Densities of the measured absolute velocities are shown in Figure 18. The ASD contains any peaks starting from 1Hz showing the large spectral content of the motion which is probably due to the triangular reference of the translation stage.

A smoother motion for the translation stage (such as a sinus motion, of a filtered triangular signal) could help reducing much of the vibrations. +The goal is to inject no motion outside the control bandwidth.

@@ -1029,7 +1055,7 @@ Thus, if the detector is only used in between the triangular peaks, the vibratio

-
+

asd_z_direction.png

Figure 18: Amplitude spectral density of the measure velocity corresponding to the geophone in the vertical direction located on the granite and at the sample location when the translation stage is scanning at 1Hz

@@ -1038,11 +1064,11 @@ Thus, if the detector is only used in between the triangular peaks, the vibratio
-
-

3.4 Open Loop noise budgeting

+
+

3.4 Open Loop noise budgeting

- +

@@ -1050,7 +1076,7 @@ We can now compare the effect of all the disturbance sources on the position err

-The Power Spectral Density of the motion error due to the ground motion, translation stage scans and spindle rotation are shown in Figure 19. +The Power Spectral Density of the motion error due to the ground motion, translation stage scans and spindle rotation are shown in Figure 19.

@@ -1058,26 +1084,26 @@ We can see that the ground motion is quite small compare to the translation stag

-
+

dist_effect_relative_motion.png

Figure 19: Amplitude Spectral Density fo the motion error due to disturbances

-The Cumulative Amplitude Spectrum is shown in Figure 20. +The Cumulative Amplitude Spectrum is shown in Figure 20. It is shown that the motion induced by translation stage scans and spindle rotation are in the micro-meter range for frequencies above 1Hz.

-
+

dist_effect_relative_motion_cas.png

Figure 20: Cumulative Amplitude Spectrum of the motion error due to disturbances

-From Figure 20, required bandwidth can be estimated by seeing that \(10\ nm [rms]\) motion is induced by the perturbations above 100Hz. +From Figure 20, required bandwidth can be estimated by seeing that \(10\ nm [rms]\) motion is induced by the perturbations above 100Hz.

@@ -1090,8 +1116,8 @@ From that, it can be concluded that control bandwidth will have to be around 100

-
-

3.5 Better estimation of the disturbances

+
+

3.5 Better estimation of the disturbances

All the disturbance measurements were made with inertial sensors, and to obtain the relative motion sample/granite, two inertial sensors were used and the signals were subtracted. @@ -1111,8 +1137,8 @@ The detector requirement would need to have a sample frequency above \(400Hz\) a

-
-

3.6 Conclusion

+
+

3.6 Conclusion

@@ -1137,14 +1163,14 @@ This should however not change the conclusion of this study nor significantly ch

-
-

4 Multi Body Model

+
+

4 Multi Body Model

- +

-As was shown during the modal analysis (Section 2), the micro-station behaves as multiple rigid bodies (granite, translation stage, tilt stage, spindle, hexapod) connected with some discrete flexibility (stiffnesses and dampers). +As was shown during the modal analysis (Section 2), the micro-station behaves as multiple rigid bodies (granite, translation stage, tilt stage, spindle, hexapod) connected with some discrete flexibility (stiffnesses and dampers).

@@ -1157,11 +1183,11 @@ A small summary of the multi-body Simscape is available -

4.1 Multi-Body model

+
+

4.1 Multi-Body model

- +

@@ -1185,11 +1211,11 @@ Then, the values of the stiffnesses and damping properties of each joint is manu

-The 3D representation of the simscape model is shown in Figure 21. +The 3D representation of the simscape model is shown in Figure 21.

-
+

simscape_picture.png

Figure 21: 3D representation of the simscape model

@@ -1197,11 +1223,11 @@ The 3D representation of the simscape model is shown in Figure -

4.2 Validity of the model’s dynamics

+
+

4.2 Validity of the model’s dynamics

- +

@@ -1209,7 +1235,7 @@ Tuning the dynamics of such model is very difficult as there are more than 50 pa

-The comparison of three of the Frequency Response Functions are shown in Figure 22. +The comparison of three of the Frequency Response Functions are shown in Figure 22.

@@ -1221,7 +1247,7 @@ We believe that the model is representing the micro-station dynamics with suffic

-
+

identification_comp_top_stages.png

Figure 22: Frequency Response function from Hammer force in the X,Y and Z directions to the X,Y and Z displacements of the micro-hexapod’s top platform. The measurements are shown in blue and the Model in red.

@@ -1254,11 +1280,11 @@ Then, using the model, we can
-
-

4.3 Wanted position of the sample and position error

+
+

4.3 Wanted position of the sample and position error

- +

@@ -1266,7 +1292,7 @@ For the control of the nano-hexapod, the sample position error (the motion to be

-To do so, several computations are performed (summarized in Figure 23): +To do so, several computations are performed (summarized in Figure 23):

  • First, the wanted pose (3 translations and 3 rotations) of the sample with respect to the granite is computed. @@ -1280,7 +1306,7 @@ Both computation are performed
-
+

control-schematic-nass.png

Figure 23: Figure caption

@@ -1292,11 +1318,11 @@ More details about these computations are accessible -

4.4 Simulation of Experiments

+
+

4.4 Simulation of Experiments

- +

@@ -1306,20 +1332,21 @@ Now that the dynamics of the model is tuned and the disturbances included in the

A first simulation is done with the nano-hexapod modeled as a rigid-body. This does represent the system without the NASS and permits to estimate the sample’s vibrations using the micro-station alone. -The results of this simulation will be compared to simulations using the NASS in Section 6.4. +The results of this simulation will be compared to simulations using the NASS in Section 6.4.

-An 3D animation of the simulation is shown in Figure 24. +An 3D animation of the simulation is shown in Figure 24.

-A zoom in the micro-meter ranger on the sample’s location is shown in Figure 25 with two frames: +A zoom in the micro-meter ranger on the sample’s location is shown in Figure 25 with two frames:

    -
  • a non-rotating frame corresponding to the wanted position of the sample. -Note that this frame is moving with the granite.
    • +
  • a non-rotating frame corresponding to the focusing point of the X-ray. +It does in that case correspond to the wanted position of the sample. +Note that this frame is moving with the granite as the nano-focusing optics are fixed to the granite.
  • a rotating frame that corresponds to the actual pose of the sample
@@ -1328,7 +1355,7 @@ The motion of the sample follows the wanted motion but with vibrations in the mi

-
+

open_loop_sim.gif

Figure 24: Tomography Experiment using the Simscape Model

@@ -1336,14 +1363,14 @@ The motion of the sample follows the wanted motion but with vibrations in the mi -
+

open_loop_sim_zoom.gif

Figure 25: Tomography Experiment using the Simscape Model - Zoom on the sample’s position (the full vertical scale is \(\approx 10 \mu m\))

-The position error of the sample with respect to the granite are shown in Figure 26. +The position error of the sample with respect to the granite are shown in Figure 26. It is confirmed that the X-Y-Z position errors are in the micro-meter range.

@@ -1361,7 +1388,7 @@ The vertical rotation error is meaningless for two reasons: -
+

exp_scans_rz_dist.png

Figure 26: Position error of the Sample with respect to the granite during a Tomography Experiment with included disturbances

@@ -1369,8 +1396,8 @@ The vertical rotation error is meaningless for two reasons:
-
-

4.5 Conclusion

+
+

4.5 Conclusion

@@ -1395,14 +1422,14 @@ In the next sections, it will allows to optimally design the nano-hexapod, to de

-
-

5 Optimal Nano-Hexapod Design

+
+

5 Optimal Nano-Hexapod Design

- +

-As explain before, the nano-hexapod properties (mass, stiffness, architecture, …) will influence: +As explain before, the nano-hexapod properties (mass, stiffness, legs’ orientation, …) will influence:

  • the effect of disturbances \(G_d\) (important for the rejection of disturbances)
    • @@ -1413,24 +1440,29 @@ As explain before, the nano-hexapod properties (mass, stiffness, architecture, & Thus, we here wish to find the optimal nano-hexapod properties such that:

      -
    • the effect of disturbances is minimized (Section 5.1)
      • -
    • the plant uncertainty due to a change of payload mass and experimental conditions is minimized (Section 5.2)
      • +
    • the effect of disturbances is minimized (Section 5.1)
      • +
    • the plant uncertainty due to a change of payload mass and experimental conditions is minimized (Section 5.2)
      • +
    • the plant has nice dynamical properties for control (Section 5.3)

    -The study presented here only consider changes in the nano-hexapod stiffness for two reasons: +In this study, the effect of the nano-hexapod’s mass characteristics is not taken into account because: +

    +
      +
    1. it cannot be changed a lot
      2. +
    2. it is quite negligible compare the to metrology reflector and the payload’s masses that is fixed to nano-hexapod’s top platform
      3. +
    + +

    +Also, the effect of the nano-hexapod’s damping properties will be studied when applying active damping techniques.

    -
      -
    • the nano-hexapod mass cannot be change too much, and will anyway be negligible compare to the metrology reflector and the payload masses
      • -
    • the choice of the nano-hexapod architecture (e.g. orientations of the actuators and implementation of sensors) will be further studied in accord with the control architecture
      • -
-
-

5.1 Optimal Stiffness to reduce the effect of disturbances

+
+

5.1 Optimal Stiffness to reduce the effect of disturbances

- +

As will be seen, the nano-hexapod stiffness have a large influence on the sensibility to disturbance (the norm of \(G_d\)). @@ -1442,11 +1474,11 @@ A complete study of the optimal nano-hexapod stiffness for the minimization of d

-
-

Sensibility to stage vibrations

-
+
+

Sensibility to stage vibrations

+

-The sensibility to the spindle’s vibration for all the considered nano-hexapod stiffnesses (from \(10^3\,[N/m]\) to \(10^9\,[N/m]\)) is shown in Figure 27. +The sensibility to the spindle’s vibration for all the considered nano-hexapod stiffnesses (from \(10^3\,[N/m]\) to \(10^9\,[N/m]\)) is shown in Figure 27. It is shown that a softer nano-hexapod it better to filter out vertical vibrations of the spindle. More precisely, is start to filters the vibration at the first suspension mode of the payload on top of the nano-hexapod.

@@ -1456,7 +1488,7 @@ The same conclusion is made for vibrations of the translation stage.

-
+

opt_stiff_sensitivity_Frz.png

Figure 27: Sensitivity to Spindle vertical motion error to the vertical error position of the sample

@@ -1464,21 +1496,21 @@ The same conclusion is made for vibrations of the translation stage.
-
-

Sensibility to ground motion

-
+
+

Sensibility to ground motion

+

-The sensibilities to ground motion in the Y and Z directions are shown in Figure 28. +The sensibilities to ground motion in the Y and Z directions are shown in Figure 28. We can see that above the suspension mode of the nano-hexapod, the norm of the transmissibility is close to one until the suspension mode of the granite. Thus, a stiff nano-hexapod is better for reducing the effect of ground motion at low frequency.

-It will be suggested in Section 7.4 that using soft mounts for the granite can greatly lower the sensibility to ground motion. +It will be suggested in Section 7.6 that using soft mounts for the granite can greatly lower the sensibility to ground motion.

-
+

opt_stiff_sensitivity_Dw.png

Figure 28: Sensitivity to Ground motion to the position error of the sample

@@ -1486,13 +1518,13 @@ It will be suggested in Section 7.4 that using soft mo
-
-

Dynamic Noise Budgeting considering all the disturbances

-
+
+

Dynamic Noise Budgeting considering all the disturbances

+

However, lowering the sensibility to some disturbance at a frequency where its effect is already small compare to the other disturbances sources is not really interesting. What is more important than comparing the sensitivity to disturbances, is thus to compare the obtain power spectral density of the sample’s position error. -From the Power Spectral Density of all the sources of disturbances identified in Section 3, we compute what would be the Power Spectral Density of the vertical motion error for all the considered nano-hexapod stiffnesses (Figure 29). +From the Power Spectral Density of all the sources of disturbances identified in Section 3, we compute what would be the Power Spectral Density of the vertical motion error for all the considered nano-hexapod stiffnesses (Figure 29).

@@ -1500,7 +1532,7 @@ We can see that the most important change is in the frequency range 30Hz to 300H

-
+

opt_stiff_psd_dz_tot.png

Figure 29: Amplitude Spectral Density of the Sample vertical position error due to Vertical vibration of the Spindle for multiple nano-hexapod stiffnesses

@@ -1508,13 +1540,13 @@ We can see that the most important change is in the frequency range 30Hz to 300H

-If we look at the Cumulative amplitude spectrum of the vertical error motion in Figure 30, we can observe that a soft hexapod (\(k < 10^5 - 10^6\,[N/m]\)) helps reducing the high frequency disturbances, and thus a smaller control bandwidth will suffice to obtain the wanted performance. +If we look at the Cumulative amplitude spectrum of the vertical error motion in Figure 30, we can observe that a soft hexapod (\(k < 10^5 - 10^6\,[N/m]\)) helps reducing the high frequency disturbances, and thus a smaller control bandwidth will suffice to obtain the wanted performance.

-
+

opt_stiff_cas_dz_tot.png

Figure 30: Cumulative Amplitude Spectrum of the Sample vertical position error due to all considered perturbations for multiple nano-hexapod stiffnesses

@@ -1523,11 +1555,11 @@ If we look at the Cumulative amplitude spectrum of the vertical error motion in
-
-

5.2 Optimal Stiffness to reduce the plant uncertainty

+
+

5.2 Optimal Stiffness to reduce the plant uncertainty

- +

One of the most important design goal is to obtain a system that is robust to all changes in the system. @@ -1559,15 +1591,15 @@ However, the dynamics from forces to sensors located in the nano-hexapod legs, s

-
-

Effect of Payload

-
+
+

Effect of Payload

+

The most obvious change in the system is the change of payload.

-In Figure 31 the dynamics is shown for payloads having a first resonance mode at 100Hz and a mass equal to 1kg, 20kg and 50kg. +In Figure 31 the dynamics is shown for payloads having a first resonance mode at 100Hz and a mass equal to 1kg, 20kg and 50kg. On the left side, the change of dynamics is computed for a very soft nano-hexapod, while on the right side, it is computed for a very stiff nano-hexapod.

@@ -1585,14 +1617,14 @@ For the stiff-nano-hexapod, the change of payload mass has very little effect (t

-
+

opt_stiffness_payload_mass_fz_dz.png

Figure 31: Dynamics from \(\mathcal{F}_z\) to \(\mathcal{X}_z\) for varying payload mass, both for a soft nano-hexapod (left) and a stiff nano-hexapod (right)

-In Figure 32 is shown the effect of a change of payload dynamics. +In Figure 32 is shown the effect of a change of payload dynamics. The mass of the payload is fixed and its resonance frequency is changing from 50Hz to 500Hz.

@@ -1601,7 +1633,7 @@ We can see (more easily for the soft nano-hexapod), that resonance of the payloa

-
+

opt_stiffness_payload_freq_fz_dz.png

Figure 32: Dynamics from \(\mathcal{F}_z\) to \(\mathcal{X}_z\) for varying payload resonance frequency, both for a soft nano-hexapod and a stiff nano-hexapod

@@ -1609,7 +1641,7 @@ We can see (more easily for the soft nano-hexapod), that resonance of the payloa

-The dynamics for all the payloads (mass from 1kg to 50kg and first resonance from 50Hz to 500Hz) and all the considered nano-hexapod stiffnesses are display in Figure 33. +The dynamics for all the payloads (mass from 1kg to 50kg and first resonance from 50Hz to 500Hz) and all the considered nano-hexapod stiffnesses are display in Figure 33.

@@ -1630,7 +1662,7 @@ For nano-hexapod stiffnesses above \(10^7\,[N/m]\): -

+

opt_stiffness_payload_impedance_all_fz_dz.png

Figure 33: Dynamics from \(\mathcal{F}_z\) to \(\mathcal{X}_z\) for varying payload dynamics, both for a soft nano-hexapod and a stiff nano-hexapod

@@ -1659,11 +1691,11 @@ Heavy samples with low first resonance mode will be very problematic.
-
-

Effect of Micro-Station Compliance

-
+
+

Effect of Micro-Station Compliance

+

-The micro-station dynamics is quite complex as was shown in Section 2, moreover, its dynamics can change due to: +The micro-station dynamics is quite complex as was shown in Section 2, moreover, its dynamics can change due to:

  • a change in some mechanical elements
    • @@ -1684,7 +1716,7 @@ This as several other advantages:

    -To identify the effect of the micro-station compliance on the system dynamics, for each nano-hexapod stiffness, we identify the plant dynamics in two different case (Figure 34): +To identify the effect of the micro-station compliance on the system dynamics, for each nano-hexapod stiffness, we identify the plant dynamics in two different case (Figure 34):

    • without the micro-station (solid curves)
      • @@ -1700,7 +1732,7 @@ For nano-hexapod stiffnesses above \(10^7\,[N/m]\), the micro-station compliance

      -
      +

      opt_stiffness_micro_station_fx_dx.png

      Figure 34: Change of dynamics from force \(\mathcal{F}_x\) to displacement \(\mathcal{X}_x\) due to the micro-station compliance

      @@ -1720,15 +1752,15 @@ If a stiff nano-hexapod is used, the control bandwidth should probably be limite
      -
      -

      Effect of Spindle Rotating Speed

      -
      +
      +

      Effect of Spindle Rotating Speed

      +

      Let’s now consider the rotation of the Spindle.

      -The plant dynamics for spindle rotation speed from 0rpm up to 60rpm are shown in Figure 35. +The plant dynamics for spindle rotation speed from 0rpm up to 60rpm are shown in Figure 35.

      @@ -1740,7 +1772,7 @@ For very soft nano-hexapods, the main resonance is split into two resonances and

      -
      +

      opt_stiffness_wz_fx_dx.png

      Figure 35: Change of dynamics from force \(\mathcal{F}_x\) to displacement \(\mathcal{X}_x\) for a spindle rotation speed from 0rpm to 60rpm

      @@ -1759,16 +1791,16 @@ A very soft (\(k < 10^4\,[N/m]\)) nano-hexapod should not be used due to the eff
      -
      -

      Total Plant Uncertainty

      -
      +
      +

      Total Plant Uncertainty

      +

      -Finally, let’s combined all the uncertainties and display the plant dynamics “spread” for all the nano-hexapod stiffnesses (Figure 36). +Finally, let’s combined all the uncertainties and display the plant dynamics “spread” for all the nano-hexapod stiffnesses (Figure 36). This show how the dynamics evolves with the stiffness and how different effects enters the plant dynamics.

      -
      +

      opt_stiffness_plant_dynamics_task_space.gif

      Figure 36: Variability of the dynamics from \(\bm{\mathcal{F}}_x\) to \(\bm{\mathcal{X}}_x\) with varying nano-hexapod stiffness

      @@ -1801,22 +1833,213 @@ In such case, the main limitation will be heavy samples with small stiffnesses.
      -
      -

      5.3 Conclusion

      +
      +

      5.3 Nano-Hexapod Architecture

      +

      + +

      +
      + +
      +

      Kinematic Analysis - Jacobian Matrix

      +
      +

      +The kinematic analysis of a parallel manipulator is well described in taghirad13_paral: +

      +
      +

      +Kinematic analysis refers to the study of the geometry of motion of a robot, without considering the forces an torques that cause the motion. +In this analysis, the relation between the geometrical parameters of the manipulator with the final motion of the moving platform is derived and analyzed. +

      +
      + + +

      +From taghirad13_paral: +

      +
      +

      +The Jacobian matrix not only reveals the relation between the joint variable velocities of a parallel manipulator to the moving platform linear and angular velocities, it also constructs the transformation needed to find the actuator forces from the forces and moments acting on the moving platform. +

      +
      + +

      +The Jacobian matrix \(\bm{\mathcal{J}}\) can be computed form the orientation of the legs and the position of the flexible joints. +

      + +

      +If we note: +

      +
        +
      • \(\delta\bm{\mathcal{L}} = [ \delta l_1, \delta l_2, \delta l_3, \delta l_4, \delta l_5, \delta l_6 ]^T\): the vector of small legs’ displacements
        • +
      • \(\delta \bm{\mathcal{X}} = [\delta x, \delta y, \delta z, \delta \theta_x, \delta \theta_y, \delta \theta_z ]^T\): the vector of small mobile platform displacements
        • +
      + +

      +The Jacobian matrix links the two vectors: +

      +\begin{align*} + \delta\bm{\mathcal{L}} &= \bm{J} \delta\bm{\mathcal{X}} \\ + \delta\bm{\mathcal{X}} &= \bm{J}^{-1} \delta\bm{\mathcal{L}} +\end{align*} + + +

      +If we note: +

      +
        +
      • \(\bm{\tau} = [\tau_1, \tau_2, \cdots, \tau_6]^T\): vector of actuator forces applied in each strut
        • +
      • \(\bm{\mathcal{F}} = [\bm{f}, \bm{n}]^T\): external force/torque action on the mobile platform
        • +
      + +\begin{equation*} + \bm{\mathcal{F}} = \bm{J}^T \bm{\tau} +\end{equation*} + + + +\begin{equation*} + \bm{\mathcal{F}} = \bm{K} \delta \bm{\mathcal{X}} +\end{equation*} + +\begin{equation*} + \bm{K} = \bm{J}^T \mathcal{K} \bm{J} +\end{equation*} +\begin{equation*} + \bm{C} = \bm{K}^{-1} = (\bm{J}^T \mathcal{K} \bm{J})^{-1} +\end{equation*} + + + +

      +Kinematic Study https://tdehaeze.github.io/stewart-simscape/kinematic-study.html +

      + + +

      +Mobility can be estimated from the architecture of the Stewart platform and the leg’s stroke. +

      + + +

      +Stiffness properties is estimated from the architecture and leg’s stiffness +

      +
      +
      + +
      +

      Kinematic Analysis - Mobility

      +
      + +
      +

      mobility_translations_null_rotation.png +

      +

      Figure 37: Figure caption

      +
      +
      +
      + +
      +

      Kinematic Study

      +
      +
      +
      + +
      +

      Flexible Joints

      +
      +

      +Active Damping Study https://tdehaeze.github.io/stewart-simscape/control-active-damping.html +Flexible Joint stiffness => not problematic for the chosen active damping technique +

      + + + + +++ ++ + + + + + + + + + + + + + + + + + + + + + +
      Table 1: Stiffness of flexible joints
      \(k_{\theta u},\ k_{\psi u}\)\(72 Nm/rad\)
      \(k_{\theta s}\)\(51 Nm/rad\)
      \(k_{\psi s}\)\(62 Nm/rad\)
      \(k_{\gamma s}\)\(64 Nm/rad\)
      + + +
      +

      preumont07_flexible_joints.png +

      +

      Figure 38: Figure caption preumont07_six_axis_singl_stage_activ

      +
      + + + +
      +

      yang19_flexible_joints.png +

      +

      Figure 39: Figure caption

      +
      +
      +
      + + +
      +

      Cubic Architecture

      +
      +

      +Study of cubic architecture https://tdehaeze.github.io/stewart-simscape/cubic-configuration.html +Has some advantages such as uniform stiffness and uniform mobility. +It can have very nice properties in specific conditions that will not be the case for this application. +The cubic configuration also puts much restriction on the position and orientation of each leg. +This configuration is such not recommended. +

      + + +
      +

      3d-cubic-stewart-aligned.png +

      +

      Figure 40: Figure caption

      +
      +
      +
      +
      + +
      +

      5.4 Conclusion

      +

      -In Section 5.1, it has been concluded that a nano-hexapod stiffness below \(10^5-10^6\,[N/m]\) helps reducing the high frequency vibrations induced by all sources of disturbances considered. +In Section 5.1, it has been concluded that a nano-hexapod stiffness below \(10^5-10^6\,[N/m]\) helps reducing the high frequency vibrations induced by all sources of disturbances considered. As the high frequency vibrations are the most difficult to compensate for when using feedback control, a soft hexapod will most certainly helps improving the performances.

      -In Section 5.2, we concluded that a nano-hexapod leg stiffness in the range \(10^5 - 10^6\,[N/m]\) is a good compromise between the uncertainty induced by the micro-station dynamics and by the rotating speed. +In Section 5.2, we concluded that a nano-hexapod leg stiffness in the range \(10^5 - 10^6\,[N/m]\) is a good compromise between the uncertainty induced by the micro-station dynamics and by the rotating speed. Provided that the samples used have a first mode that is sufficiently high in frequency, the total plant dynamic uncertainty should be manageable.

      -Thus, a stiffness of \(10^5\,[N/m]\) will be used in Section 6 to develop the robust control architecture and to perform simulations. +Thus, a stiffness of \(10^5\,[N/m]\) will be used in Section 6 to develop the robust control architecture and to perform simulations.

      @@ -1828,11 +2051,11 @@ A more detailed study of the determination of the optimal stiffness based on all

      -
      -

      6 Robust Control Architecture

      +
      +

      6 Robust Control Architecture

      - +

      Before designing the control system, let’s summarize what has been done: @@ -1859,8 +2082,8 @@ This would require to measure the mass/inertia of each used payload and manually

      -
      -

      6.1 High Authority Control / Low Authority Control Architecture

      +
      +

      6.1 High Authority Control / Low Authority Control Architecture

      Many control architecture could be used for the control of the nano-hexapod. @@ -1876,7 +2099,7 @@ Some properties of the HAC-LAC architecture are explained below (taken from

      -The HAC/LAC approach consist of combining the two approached in a dual-loop control as shown in Figure 37. +The HAC/LAC approach consist of combining the two approached in a dual-loop control as shown in Figure 41. 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:

      @@ -1888,42 +2111,36 @@ This approach has the following advantages:
      -
      +

      control_architecture_hac_lac_one_input.png

      -

      Figure 37: HAC-LAC Architecture with a system having only one input

      +

      Figure 41: HAC-LAC Architecture with a system having only one input

      The HAC-LAC architecture thus consisted of two cascade controllers:

        -
      • a Low Authority Controller that is used to damp the system (Section 6.2)
        • -
      • a High Authority Controller used to suppress the sample’s vibration in a wide frequency range (Section 6.3)
        • +
      • a Low Authority Controller that is used to damp the system (Section 6.2)
        • +
      • a High Authority Controller used to suppress the sample’s vibration in a wide frequency range (Section 6.3)
      -
      -

      6.2 Active Damping and Sensors to be included in the nano-hexapod

      +
      +

      6.2 Active Damping and Sensors to be included in the nano-hexapod

      - +

      - -

      -Active Damping can help: -

      -
        -
      • by reducing the effect of disturbances close to the resonance of the system
        • -
      • by making the plant dynamics simpler to control for the High Authority Controller
        • -
      - -

      Depending on the chosen active damping technique, either force sensors, relative motion sensors or inertial sensors should be included in each of the nano-hexapod’s legs.

      +

      +Because of the rotation of the hexapod, +

      +

      A separate study (accessible here) about the use of all three sensors types have been done, the conclusions are:

      @@ -1933,14 +2150,64 @@ A separate study (accessible +

      Effect of the Spindle’s Rotation

      +
      + +
      +

      dvf_root_locus_ws.png +

      +

      Figure 42: Figure caption

      +
      + + +
      +

      iff_root_locus_ws.png +

      +

      Figure 43: Figure caption

      +
      +
      +
      + +
      +

      Effect of Flexible Joint

      +
      + +
      +

      root_locus_iff_rot_stiffness.png +

      +

      Figure 44: Figure caption

      +
      + + +
      +

      root_locus_dvf_rot_stiffness.png +

      +

      Figure 45: Figure caption

      +
      +
      +
      + +
      +

      Relative Direct Velocity Feedback Architecture

      +
      +

      +Active Damping can help: +

      +
        +
      • by reducing the effect of disturbances close to the resonance of the system
        • +
      • by making the plant dynamics simpler to control for the High Authority Controller
        • +

      -Relative motion sensors are then included in each of the nano-hexapod’s leg and a decentralized direct velocity feedback control architecture is applied (Figure 38). +Relative motion sensors are then included in each of the nano-hexapod’s leg and a decentralized direct velocity feedback control architecture is applied (Figure 46).

      -The signals shown in Figure 38 are: +The signals shown in Figure 46 are:

      • \(\bm{\tau}\): Actuator forces applied in each leg
        • @@ -1956,33 +2223,42 @@ The force applied in each leg being proportional to the relative velocity of the

        -
        +

        control_architecture_dvf.png

        -

        Figure 38: Low Authority Control: Decentralized Direct Velocity Feedback

        +

        Figure 46: Low Authority Control: Decentralized Direct Velocity Feedback

        The DVF gain is here chosen in such a way that the suspension modes of the nano-hexapod are critically damped whatever the sample mass. This may not be the optimal choice as will be further explained.

        +
        +
      - +
      +

      Effect of Active Damping on the Primary Plant Dynamics

      +

      -The plant dynamics before (solid curves) and after (dashed curves) the Low Authority Control implementation are compared in Figure 39. +The plant dynamics before (solid curves) and after (dashed curves) the Low Authority Control implementation are compared in Figure 47. It is clear that the use of the DVF reduces the dynamical spread of the plant dynamics between 5Hz up too 100Hz. This will make the primary controller more robust and easier to develop.

      -
      +

      opt_stiff_primary_plant_damped_L.png

      -

      Figure 39: Primary plant in the space of the legs with (dashed) and without (solid) Direct Velocity Feedback

      +

      Figure 47: Primary plant in the space of the legs with (dashed) and without (solid) Direct Velocity Feedback

      +
      +
      +
      +

      Effect of Active Damping on the Sensibility to Disturbances

      +

      -The change of sensibility to disturbances with the use of DVF is shown in Figure 40. +The change of sensibility to disturbances with the use of DVF is shown in Figure 48. It is shown that the DVF control lowers the sensibility to disturbances in the vicinity of the nano-hexapod resonance but increases the sensibility at higher frequencies.

      @@ -1991,28 +2267,29 @@ This is probably not the optimal gain that could have been used, and further ana

      -
      +

      opt_stiff_sensibility_dist_dvf.png

      -

      Figure 40: Norm of the transfer function from vertical disturbances to vertical position error with (dashed) and without (solid) Direct Velocity Feedback applied

      +

      Figure 48: Norm of the transfer function from vertical disturbances to vertical position error with (dashed) and without (solid) Direct Velocity Feedback applied

      +
      -
      -

      6.3 High Authority Control

      +
      +

      6.3 High Authority Control

      - +

      -The complete HAC-LAC architecture is shown in Figure 41 where an outer loop is added to the decentralized direct velocity feedback loop. +The complete HAC-LAC architecture is shown in Figure 49 where an outer loop is added to the decentralized direct velocity feedback loop.

      The block Compute Position Error is used to compute the position error \(\bm{\epsilon}_{\mathcal{X}_n}\) of the sample with respect to the nano-hexapod’s base platform from the actual measurement of the sample’s pose \(\bm{\mathcal{X}}\) and the wanted pose \(\bm{r}_\mathcal{X}\). -The computation done in such block was briefly explained in Section 4.3. +The computation done in such block was briefly explained in Section 4.3.

      @@ -2025,10 +2302,10 @@ Then, a diagonal controller \(\bm{K}_\mathcal{L}\) generates the required force

      -
      +

      control_architecture_hac_dvf_pos_L.png

      -

      Figure 41: Cascade Control Architecture. The inner loop consist of a decentralized Direct Velocity Feedback. The outer loop consist of position control in the leg’s space

      +

      Figure 49: Cascade Control Architecture. The inner loop consist of a decentralized Direct Velocity Feedback. The outer loop consist of position control in the leg’s space

      @@ -2036,32 +2313,32 @@ Some alternative to this control architecture have been studied, but this is the

      -The plant dynamics for each of the six legs and for the three payload’s masses is shown in Figure 42. +The plant dynamics for each of the six legs and for the three payload’s masses is shown in Figure 50. The dynamical spread is kept reasonably small thanks to both the optimal nano-hexapod design and the Low Authority Controller.

      -
      +

      opt_stiff_primary_plant_L.png

      -

      Figure 42: Diagonal elements of the transfer function matrix from \(\bm{\tau}^\prime\) to \(\bm{\epsilon}_{\mathcal{X}_n}\) for the three considered masses

      +

      Figure 50: Diagonal elements of the transfer function matrix from \(\bm{\tau}^\prime\) to \(\bm{\epsilon}_{\mathcal{X}_n}\) for the three considered masses

      The diagonal controller \(\bm{K}_\mathcal{L}\) is then tuned in such a way that the control bandwidth is around 100Hz and such that enough stability margins are obtained for all the payload’s masses used. -The obtained loop gain is shown in Figure 43. +The obtained loop gain is shown in Figure 51.

      -
      +

      opt_stiff_primary_loop_gain_L.png

      -

      Figure 43: Loop gain for the primary plant

      +

      Figure 51: Loop gain for the primary plant

      -The sensibility to disturbance after the use of HAC-LAC control is shown in Figure 44. +The sensibility to disturbance after the use of HAC-LAC control is shown in Figure 52. The change of sensibility is very typical for feedback system:

        @@ -2079,30 +2356,30 @@ This should gives slightly better performance and robustness, but should not cha

        -
        +

        opt_stiff_primary_control_L_senbility_dist.png

        -

        Figure 44: Sensibility to disturbances when the HAC-LAC control is applied (dashed) and when it is not (solid)

        +

        Figure 52: Sensibility to disturbances when the HAC-LAC control is applied (dashed) and when it is not (solid)

      -
      -

      6.4 Simulation of Tomography Experiments

      +
      +

      6.4 Simulation of Tomography Experiments

      - +

      A new simulation of a tomography is performed with the optimal nano-hexapod and the HAC-LAC architecture implemented in the model. -The results of this simulation will be compare to the simulation performed in Section 4.4 without the nano-hexapod. +The results of this simulation will be compare to the simulation performed in Section 4.4 without the nano-hexapod. All the disturbances are included such as ground motion, spindle and translation stage vibrations.

      -The Power Spectral Density of the sample’s position error is plotted in Figure 45 and the Cumulative Amplitude Spectrum is shown in Figure 46. +The Power Spectral Density of the sample’s position error is plotted in Figure 53 and the Cumulative Amplitude Spectrum is shown in Figure 54. The top three plots corresponds to the X, Y and Z translations and the bottom three plots corresponds to the X,Y and Z rotations.

      @@ -2131,49 +2408,78 @@ This increase in rotation is still very small and is not foreseen to be a proble
    -
    +

    opt_stiff_hac_dvf_L_psd_disp_error.png

    -

    Figure 45: Amplitude Spectral Density of the position error in Open Loop and with the HAC-LAC controller

    +

    Figure 53: Amplitude Spectral Density of the position error in Open Loop and with the HAC-LAC controller

    -
    +

    opt_stiff_hac_dvf_L_cas_disp_error.png

    -

    Figure 46: Cumulative Amplitude Spectrum of the position error in Open Loop and with the HAC-LAC controller

    +

    Figure 54: Cumulative Amplitude Spectrum of the position error in Open Loop and with the HAC-LAC controller

    -The time domain sample’s vibrations are shown in Figure 47. +The time domain sample’s vibrations are shown in Figure 55. The use of the nano-hexapod combined with the HAC-LAC architecture is shown to considerably reduce the sample’s vibrations.

    -An animation of the experiment is shown in Figure 48 and we can see that the actual sample’s position is more closely following the ideal position compared to the simulation of the micro-station alone in Figure 25 (same scale was used for both animations). +An animation of the experiment is shown in Figure 56 and we can see that the actual sample’s position is more closely following the ideal position compared to the simulation of the micro-station alone in Figure 25 (same scale was used for both animations).

    -
    +

    opt_stiff_hac_dvf_L_pos_error.png

    -

    Figure 47: Position Error of the sample during a tomography experiment when no control is applied and with the HAC-DVF control architecture

    +

    Figure 55: Position Error of the sample during a tomography experiment when no control is applied and with the HAC-DVF control architecture

    -
    +

    closed_loop_sim_zoom.gif

    -

    Figure 48: Tomography Experiment using the Simscape Model in Closed Loop with the HAC-LAC Control - Zoom on the sample’s position (the full vertical scale is \(\approx 10 \mu m\))

    +

    Figure 56: Tomography Experiment using the Simscape Model in Closed Loop with the HAC-LAC Control - Zoom on the sample’s position (the full vertical scale is \(\approx 10 \mu m\))

    -
    -

    6.5 Conclusion

    +
    +

    6.5 Simulation of More Complex Experiments

    +
    +
    +

    Micro-Hexapod offset

    +
    + +
    +

    tomography_dh_offset.gif +

    +

    Figure 57: Top View of a tomography experiment with a 10mm offset imposed by the micro-hexapod

    +
    +
    +
    + +
    +

    Simultaneous Translation Scans

    +
    + +
    +

    ty_scans.gif +

    +

    Figure 58: Top View of a tomography experiment combined with translation scans

    +
    +
    +
    +
    + +
    +

    6.6 Conclusion

    +

    The High Authority Control / Low Authority Control architecture has implemented on the multi-body model of the system. @@ -2207,22 +2513,26 @@ A more complete study of the control of the NASS is performed -

    7 General Conclusion and Further notes

    +
    +

    7 General Conclusion and Further notes

    - +

    -
    -

    7.1 General Conclusion

    +
    +

    7.1 Nano-Hexapod Specifications

    +
    + +
    +

    7.2 General Conclusion

    -
    -

    7.2 Sensor Noise introduced by the Metrology

    -
    +
    +

    7.3 Sensor Noise introduced by the Metrology

    +

    Say that is will introduce noise inside the bandwidth (100Hz) This should not be significant. @@ -2231,23 +2541,27 @@ This should not be significant.

    -
    -

    7.3 Further Work

    +
    +

    7.4 Further Work

    +
    + +
    +

    7.5 Cable Forces

    -
    -

    7.4 Using soft mounts for the

    -
    +
    +

    7.6 Using soft mounts for the Granite

    +

    - +

    -
    +

    opt_stiff_soft_granite_Dw.png

    -

    Figure 49: Change of sensibility to Ground motion when using a stiff Granite (solid curves) and a soft Granite (dashed curves)

    +

    Figure 59: Change of sensibility to Ground motion when using a stiff Granite (solid curves) and a soft Granite (dashed curves)

    @@ -2260,9 +2574,9 @@ Sensible to detector motion?

    -
    -

    7.5 Others

    -
    +
    +

    7.7 Others

    +

    Common metrology frame for the nano-focusing optics and the measurement of the sample position?

    @@ -2283,6 +2597,8 @@ Slip-Ring noise?

    Bibliography

    @@ -2290,7 +2606,7 @@ Slip-Ring noise?

    Date: 05-2020

    Author: Thomas Dehaeze

    -

    Created: 2020-04-29 mer. 19:02

    +

    Created: 2020-04-30 jeu. 14:42

    diff --git a/index.org b/index.org index 8c0a5fe..b9a4711 100644 --- a/index.org +++ b/index.org @@ -604,6 +604,7 @@ The ASD contains any peaks starting from 1Hz showing the large spectral content #+begin_important A smoother motion for the translation stage (such as a sinus motion, of a filtered triangular signal) could help reducing much of the vibrations. + The goal is to inject no motion outside the control bandwidth. We should also note that away from the rapid change of velocity, the sample's vibrations are much reduced. Thus, if the detector is only used in between the triangular peaks, the vibrations are expected to be much lower than those estimated. @@ -758,8 +759,9 @@ The results of this simulation will be compared to simulations using the NASS in An 3D animation of the simulation is shown in Figure [[fig:open_loop_sim]]. A zoom in the micro-meter ranger on the sample's location is shown in Figure [[fig:open_loop_sim_zoom]] with two frames: -- a non-rotating frame corresponding to the wanted position of the sample. - Note that this frame is moving with the granite. +- a non-rotating frame corresponding to the focusing point of the X-ray. + It does in that case correspond to the wanted position of the sample. + Note that this frame is moving with the granite as the nano-focusing optics are fixed to the granite. - a rotating frame that corresponds to the actual pose of the sample The motion of the sample follows the wanted motion but with vibrations in the micro-meter range as was expected. @@ -803,17 +805,20 @@ In the next sections, it will allows to optimally design the nano-hexapod, to de <> ** Introduction :ignore: -As explain before, the nano-hexapod properties (mass, stiffness, architecture, ...) will influence: +As explain before, the nano-hexapod properties (mass, stiffness, legs' orientation, ...) will influence: - the effect of disturbances $G_d$ (important for the rejection of disturbances) - the plant dynamics $G$ (important for the control robustness properties) Thus, we here wish to find the optimal nano-hexapod properties such that: - the effect of disturbances is minimized (Section [[sec:optimal_stiff_dist]]) - the plant uncertainty due to a change of payload mass and experimental conditions is minimized (Section [[sec:optimal_stiff_plant]]) +- the plant has nice dynamical properties for control (Section [[sec:nano_hexapod_architecture]]) -The study presented here only consider changes in the nano-hexapod *stiffness* for two reasons: -- the nano-hexapod mass cannot be change too much, and will anyway be negligible compare to the metrology reflector and the payload masses -- the choice of the nano-hexapod architecture (e.g. orientations of the actuators and implementation of sensors) will be further studied in accord with the control architecture +In this study, the effect of the nano-hexapod's mass characteristics is not taken into account because: +1. it cannot be changed a lot +2. it is quite negligible compare the to metrology reflector and the payload's masses that is fixed to nano-hexapod's top platform + +Also, the effect of the nano-hexapod's damping properties will be studied when applying active damping techniques. ** Optimal Stiffness to reduce the effect of disturbances <> @@ -1043,6 +1048,127 @@ This show how the dynamics evolves with the stiffness and how different effects In such case, the main limitation will be heavy samples with small stiffnesses. #+end_important +** Nano-Hexapod Architecture +<> + +*** Introduction :ignore: + +*** Kinematic Analysis - Jacobian Matrix +:PROPERTIES: +:UNNUMBERED: t +:END: + +The kinematic analysis of a parallel manipulator is well described in cite:taghirad13_paral: +#+begin_quote +Kinematic analysis refers to the study of the geometry of motion of a robot, without considering the forces an torques that cause the motion. +In this analysis, the relation between the geometrical parameters of the manipulator with the final motion of the moving platform is derived and analyzed. +#+end_quote + + +From cite:taghirad13_paral: +#+begin_quote +The Jacobian matrix not only reveals the *relation between the joint variable velocities of a parallel manipulator to the moving platform linear and angular velocities*, it also constructs the transformation needed to find the *actuator forces from the forces and moments acting on the moving platform*. +#+end_quote + +The Jacobian matrix $\bm{\mathcal{J}}$ can be computed form the orientation of the legs and the position of the flexible joints. + +If we note: +- $\delta\bm{\mathcal{L}} = [ \delta l_1, \delta l_2, \delta l_3, \delta l_4, \delta l_5, \delta l_6 ]^T$: the vector of small legs' displacements +- $\delta \bm{\mathcal{X}} = [\delta x, \delta y, \delta z, \delta \theta_x, \delta \theta_y, \delta \theta_z ]^T$: the vector of small mobile platform displacements + +The Jacobian matrix links the two vectors: +\begin{align*} + \delta\bm{\mathcal{L}} &= \bm{J} \delta\bm{\mathcal{X}} \\ + \delta\bm{\mathcal{X}} &= \bm{J}^{-1} \delta\bm{\mathcal{L}} +\end{align*} + + +If we note: +- $\bm{\tau} = [\tau_1, \tau_2, \cdots, \tau_6]^T$: vector of actuator forces applied in each strut +- $\bm{\mathcal{F}} = [\bm{f}, \bm{n}]^T$: external force/torque action on the mobile platform + +\begin{equation*} + \bm{\mathcal{F}} = \bm{J}^T \bm{\tau} +\end{equation*} + + + +\begin{equation*} + \bm{\mathcal{F}} = \bm{K} \delta \bm{\mathcal{X}} +\end{equation*} + +\begin{equation*} + \bm{K} = \bm{J}^T \mathcal{K} \bm{J} +\end{equation*} +\begin{equation*} + \bm{C} = \bm{K}^{-1} = (\bm{J}^T \mathcal{K} \bm{J})^{-1} +\end{equation*} + + + +Kinematic Study https://tdehaeze.github.io/stewart-simscape/kinematic-study.html + + +Mobility can be estimated from the architecture of the Stewart platform and the leg's stroke. + + +Stiffness properties is estimated from the architecture and leg's stiffness + +*** Kinematic Analysis - Mobility +:PROPERTIES: +:UNNUMBERED: t +:END: + + +#+name: fig:mobility_translations_null_rotation +#+caption: Figure caption +[[file:figs/mobility_translations_null_rotation.png]] + +*** Kinematic Study +:PROPERTIES: +:UNNUMBERED: t +:END: + +*** Flexible Joints +:PROPERTIES: +:UNNUMBERED: t +:END: + +Active Damping Study https://tdehaeze.github.io/stewart-simscape/control-active-damping.html +Flexible Joint stiffness => not problematic for the chosen active damping technique + +#+name: tab:yang19_stiffness_flexible_joints +#+caption: Stiffness of flexible joints +| $k_{\theta u},\ k_{\psi u}$ | $72 Nm/rad$ | +| $k_{\theta s}$ | $51 Nm/rad$ | +| $k_{\psi s}$ | $62 Nm/rad$ | +| $k_{\gamma s}$ | $64 Nm/rad$ | + +#+name: fig:preumont07_flexible_joints +#+caption: Figure caption cite:preumont07_six_axis_singl_stage_activ +[[file:figs/preumont07_flexible_joints.png]] + + +#+name: fig:yang19_flexible_joints +#+caption: Figure caption +[[file:figs/yang19_flexible_joints.png]] + + +*** Cubic Architecture +:PROPERTIES: +:UNNUMBERED: t +:END: + +Study of cubic architecture https://tdehaeze.github.io/stewart-simscape/cubic-configuration.html +Has some advantages such as uniform stiffness and uniform mobility. +It can have very nice properties in specific conditions that will not be the case for this application. +The cubic configuration also puts much restriction on the position and orientation of each leg. +This configuration is such not recommended. + +#+name: fig:3d-cubic-stewart-aligned +#+caption: Figure caption +[[file:figs/3d-cubic-stewart-aligned.png]] + ** Conclusion #+begin_important In Section [[sec:optimal_stiff_dist]], it has been concluded that a nano-hexapod stiffness below $10^5-10^6\,[N/m]$ helps reducing the high frequency vibrations induced by all sources of disturbances considered. @@ -1102,19 +1228,39 @@ The HAC-LAC architecture thus consisted of two cascade controllers: ** Active Damping and Sensors to be included in the nano-hexapod <> -Active Damping can help: -- by reducing the effect of disturbances close to the resonance of the system -- by making the plant dynamics simpler to control for the High Authority Controller - - +*** Introduction :ignore: Depending on the chosen active damping technique, either force sensors, relative motion sensors or inertial sensors should be included in each of the nano-hexapod's legs. +Because of the rotation of the hexapod, + A separate study (accessible [[https://tdehaeze.github.io/rotating-frame/index.html][here]]) about the use of all three sensors types have been done, the conclusions are: - the use of force sensors is to be avoided as it could introduce instability in the system due to the nano-hexapod's rotation - the use of inertial sensor should not be used as it would tends to decouple the motion of the sample from the motion of the granite (which is not wanted). It would also be difficult to apply in a robust way due to the non-collocation with the actuators - relative motion sensors can be used to damped the nano-hexapod's modes in a robust way but may increase the sensibility to stages vibrations +*** Effect of the Spindle's Rotation +:PROPERTIES: +:UNNUMBERED: t +:END: + +#+name: fig:dvf_root_locus_ws +#+caption: Figure caption +[[file:figs/dvf_root_locus_ws.png]] + +#+name: fig:iff_root_locus_ws +#+caption: Figure caption +[[file:figs/iff_root_locus_ws.png]] + +*** Relative Direct Velocity Feedback Architecture +:PROPERTIES: +:UNNUMBERED: t +:END: + +Active Damping can help: +- by reducing the effect of disturbances close to the resonance of the system +- by making the plant dynamics simpler to control for the High Authority Controller + *Relative motion sensors* are then included in each of the nano-hexapod's leg and a decentralized direct velocity feedback control architecture is applied (Figure [[fig:control_architecture_dvf]]). @@ -1135,6 +1281,10 @@ The force applied in each leg being proportional to the relative velocity of the The DVF gain is here chosen in such a way that the suspension modes of the nano-hexapod are critically damped whatever the sample mass. This may not be the optimal choice as will be further explained. +*** Effect of Active Damping on the Primary Plant Dynamics +:PROPERTIES: +:UNNUMBERED: t +:END: The plant dynamics before (solid curves) and after (dashed curves) the Low Authority Control implementation are compared in Figure [[fig:opt_stiff_primary_plant_damped_L]]. It is clear that the use of the DVF reduces the dynamical spread of the plant dynamics between 5Hz up too 100Hz. @@ -1144,6 +1294,11 @@ This will make the primary controller more robust and easier to develop. #+caption: Primary plant in the space of the legs with (dashed) and without (solid) Direct Velocity Feedback [[file:figs/opt_stiff_primary_plant_damped_L.png]] +*** Effect of Active Damping on the Sensibility to Disturbances +:PROPERTIES: +:UNNUMBERED: t +:END: + The change of sensibility to disturbances with the use of DVF is shown in Figure [[fig:opt_stiff_sensibility_dist_dvf]]. It is shown that the DVF control lowers the sensibility to disturbances in the vicinity of the nano-hexapod resonance but increases the sensibility at higher frequencies. @@ -1250,14 +1405,24 @@ An animation of the experiment is shown in Figure [[fig:closed_loop_sim_zoom]] a [[file:figs/closed_loop_sim_zoom.gif]] ** Simulation of More Complex Experiments +*** Introduction :ignore: + +*** Micro-Hexapod offset +:PROPERTIES: +:UNNUMBERED: t +:END: #+name: fig:tomography_dh_offset -#+caption: Figure caption +#+caption: Top View of a tomography experiment with a 10mm offset imposed by the micro-hexapod [[file:figs/tomography_dh_offset.gif]] +*** Simultaneous Translation Scans +:PROPERTIES: +:UNNUMBERED: t +:END: #+name: fig:ty_scans -#+caption: Figure caption +#+caption: Top View of a tomography experiment combined with translation scans [[file:figs/ty_scans.gif]] ** Conclusion @@ -1281,6 +1446,8 @@ A more complete study of the control of the NASS is performed [[https://tdehaeze * General Conclusion and Further notes <> +** Nano-Hexapod Specifications + ** General Conclusion @@ -1291,8 +1458,10 @@ This should not be significant. ** Further Work +** Cable Forces -** Using soft mounts for the + +** Using soft mounts for the Granite <> #+name: fig:opt_stiff_soft_granite_Dw diff --git a/ref.bib b/ref.bib index 95e2173..4c9bc14 100644 --- a/ref.bib +++ b/ref.bib @@ -3,7 +3,6 @@ title = {The Design of High Performance Mechatronics - 2nd Revised Edition}, year = {2014}, publisher = {Ios Press}, - tags = {favorite}, } @article{oomen18_advan_motion_contr_precis_mechat, @@ -16,7 +15,6 @@ year = {2018}, doi = {10.1541/ieejjia.7.127}, url = {https://doi.org/10.1541/ieejjia.7.127}, - tags = {favorite}, } @book{preumont18_vibrat_contr_activ_struc_fourt_edition, @@ -28,5 +26,40 @@ doi = {10.1007/978-3-319-72296-2}, pages = {nil}, series = {Solid Mechanics and Its Applications}, +} + +@book{taghirad13_paral, + author = {Taghirad, Hamid}, + title = {Parallel robots : mechanics and control}, + year = {2013}, + publisher = {CRC Press}, + address = {Boca Raton, FL}, + isbn = {9781466555778}, tags = {favorite, parallel robot}, } + +@article{yang19_dynam_model_decoup_contr_flexib, + author = {Yang, XiaoLong and Wu, HongTao and Chen, Bai and Kang, ShengZheng and Cheng, ShiLi}, + title = {Dynamic Modeling and Decoupled Control of a Flexible Stewart Platform for Vibration Isolation}, + journal = {Journal of Sound and Vibration}, + volume = {439}, + pages = {398-412}, + year = {2019}, + doi = {10.1016/j.jsv.2018.10.007}, + url = {https://doi.org/10.1016/j.jsv.2018.10.007}, + issn = {0022-460X}, + month = {Jan}, + publisher = {Elsevier BV}, +} + +@article{preumont07_six_axis_singl_stage_activ, + author = {A. Preumont and M. Horodinca and I. Romanescu and B. de Marneffe and M. Avraam and A. Deraemaeker and F. Bossens and A. Abu Hanieh}, + title = {A Six-Axis Single-Stage Active Vibration Isolator Based on Stewart Platform}, + journal = {Journal of Sound and Vibration}, + volume = {300}, + number = {3-5}, + pages = {644-661}, + year = {2007}, + doi = {10.1016/j.jsv.2006.07.050}, + url = {https://doi.org/10.1016/j.jsv.2006.07.050}, +}