% Created 2024-11-05 Tue 13:26
% Intended LaTeX compiler: pdflatex
\documentclass[a4paper, 10pt, DIV=12, parskip=full, bibliography=totoc]{scrreprt}

\input{preamble.tex}
\newacronym{haclac}{HAC-LAC}{High Authority Control - Low Authority Control}
\newacronym{hac}{HAC}{High Authority Control}
\newacronym{lac}{LAC}{Low Authority Control}
\newacronym{nass}{NASS}{Nano Active Stabilization System}
\newacronym{asd}{ASD}{Amplitude Spectral Density}
\newacronym{psd}{PSD}{Power Spectral Density}
\newacronym{cps}{CPS}{Cumulative Power Spectrum}
\newacronym{cas}{CAS}{Cumulative Amplitude Spectrum}
\newacronym{frf}{FRF}{Frequency Response Function}
\newacronym{iff}{IFF}{Integral Force Feedback}
\newacronym{rdc}{RDC}{Relative Damping Control}
\newglossaryentry{ms}{name=\ensuremath{m_s},description={{Mass of the sample}}}
\newglossaryentry{mn}{name=\ensuremath{m_n},description={{Mass of the nano-hexapod}}}
\newglossaryentry{mh}{name=\ensuremath{m_h},description={{Mass of the micro-hexapod}}}
\newglossaryentry{mt}{name=\ensuremath{m_t},description={{Mass of the micro-station stages}}}
\newglossaryentry{mg}{name=\ensuremath{m_g},description={{Mass of the granite}}}
\newglossaryentry{xf}{name=\ensuremath{x_f},description={{Floor motion}}}
\newglossaryentry{ft}{name=\ensuremath{f_t},description={{Disturbance force of the micro-station}}}
\newglossaryentry{fs}{name=\ensuremath{f_s},description={{Direct forces applied on the sample}}}
\newglossaryentry{d}{name=\ensuremath{d},description={{Measured motion between the nano-hexapod and the granite}}}
\newglossaryentry{fn}{name=\ensuremath{f_n},description={{Force sensor on the nano-hexapod}}}
\newglossaryentry{psdx}{name=\ensuremath{\Phi_{x}},description={{Power spectral density of signal $x$}}}
\newglossaryentry{asdx}{name=\ensuremath{\Gamma_{x}},description={{Amplitude spectral density of signal $x$}}}
\newglossaryentry{cpsx}{name=\ensuremath{\Phi_{x}},description={{Cumulative Power Spectrum of signal $x$}}}
\newglossaryentry{casx}{name=\ensuremath{\Gamma_{x}},description={{Cumulative Amplitude Spectrum of signal $x$}}}
\input{preamble_extra.tex}
\bibliography{nass-uniaxial-model.bib}
\author{Dehaeze Thomas}
\date{\today}
\title{Nano Active Stabilization System - Uniaxial Model}
\hypersetup{
 pdfauthor={Dehaeze Thomas},
 pdftitle={Nano Active Stabilization System - Uniaxial Model},
 pdfkeywords={},
 pdfsubject={},
 pdfcreator={Emacs 29.4 (Org mode 9.6)}, 
 pdflang={English}}
\usepackage{biblatex}

\begin{document}

\maketitle
\tableofcontents

\clearpage

In this report, a uniaxial model of the \acrfull{nass} is developed and used to obtain a first idea of the challenges involved in this complex system.
Note that in this study, only the vertical direction is considered (which is the most stiff), but other directions were considered as well, yielding to similar conclusions.
The model is schematically shown in Figure \ref{fig:uniaxial_overview_model_sections} where the colors represent the parts studied in different sections.

To have a relevant model, the micro-station dynamics is first identified and its model is tuned to match the measurements (Section \ref{sec:micro_station_model}).
Then, a model of the nano-hexapod is added on top of the micro-station.
With the added sample and sensors, this gives a uniaxial dynamical model of the \acrshort{nass} that will be used for further analysis (Section \ref{sec:nano_station_model}).

The disturbances affecting position stability are identified experimentally (Section \ref{sec:uniaxial_disturbances}) and included in the model for dynamical noise budgeting (Section \ref{sec:uniaxial_noise_budgeting}).
In all the following analysis, three nano-hexapod stiffnesses are considered to better understand the trade-offs and to find the most adequate nano-hexapod design.
Three sample masses are also considered to verify the robustness of the applied control strategies with respect to a change of sample.

To improve the position stability of the sample, an \acrfull{haclac} strategy is applied.
It consists of first actively damping the plant (the \acrshort{lac} part), and then applying a position control on the damped plant (the \acrshort{hac} part).

Three active damping techniques are studied (Section \ref{sec:uniaxial_active_damping}) which are used to both reduce the effect of disturbances and make the system easier to control afterwards.
Once the system is well damped, a feedback position controller is applied and the obtained performance is analyzed (Section \ref{sec:uniaxial_position_control}).

Two key effects that may limit that positioning performances are then considered: the limited micro-station compliance (Section \ref{sec:uniaxial_support_compliance}) and the presence of dynamics between the nano-hexapod and the sample's point of interest (Section \ref{sec:uniaxial_payload_dynamics}).

\begin{figure}[htbp]
\centering
\includegraphics[scale=1]{figs/uniaxial_overview_model_sections.png}
\caption{\label{fig:uniaxial_overview_model_sections}Uniaxial Micro-Station model in blue (Section \ref{sec:micro_station_model}), Nano-Hexapod models in red (Section \ref{sec:nano_station_model}), Disturbances in yellow (Section \ref{sec:uniaxial_disturbances}), Active Damping in green (Section \ref{sec:uniaxial_active_damping}), Position control in purple (Section \ref{sec:uniaxial_position_control}) and Sample dynamics in cyan (Section \ref{sec:uniaxial_payload_dynamics})}
\end{figure}

\chapter{Micro Station Model}
\label{sec:micro_station_model}
In this section, a uniaxial model of the micro-station is tuned to match measurements made on the micro-station.
The measurement setup is shown in Figure \ref{fig:uniaxial_ustation_first_meas_dynamics} where several geophones\footnote{Mark Product L4-C geophones are used with a sensitivity of \(171\,\frac{V}{m/s}\) and a natural frequency of \(\approx 1\,\text{Hz}\)} are fixed to the micro-station and an instrumented hammer is used to inject forces on different stages of the micro-station.

From the measured frequency response functions (FRF), the model can be tuned to approximate the uniaxial dynamics of the micro-station.

\begin{figure}[htbp]
\centering
\includegraphics[scale=1,width=\linewidth]{figs/uniaxial_ustation_first_meas_dynamics.jpg}
\caption{\label{fig:uniaxial_ustation_first_meas_dynamics}Experimental setup used for the first dynamical measurements on the Micro-Station. Geophones are fixed to different stages of the micro-station.}
\end{figure}
\section{Measured dynamics}

The measurement setup is schematically shown in Figure \ref{fig:uniaxial_ustation_meas_dynamics_schematic} where two vertical hammer hits are performed, one on the Granite (force \(F_{g}\)) and the other on the micro-hexapod's top platform (force \(F_{h}\)).
The vertical inertial motion of the granite \(x_{g}\) and the top platform of the micro-hexapod \(x_{h}\) are measured using geophones.
Three frequency response functions were computed: one from \(F_{h}\) to \(x_{h}\) (i.e., the compliance of the micro-station), one from \(F_{g}\) to \(x_{h}\) (or from \(F_{h}\) to \(x_{g}\)) and one from \(F_{g}\) to \(x_{g}\).

Due to the poor coherence at low frequencies, these frequency response functions will only be shown between 20 and 200Hz (solid lines in Figure \ref{fig:uniaxial_comp_frf_meas_model}).

\begin{figure}[htbp]
\begin{subfigure}{0.69\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_ustation_meas_dynamics_schematic.png}
\end{center}
\subcaption{\label{fig:uniaxial_ustation_meas_dynamics_schematic}Measurement setup - Schematic}
\end{subfigure}
\begin{subfigure}{0.29\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_model_micro_station.png}
\end{center}
\subcaption{\label{fig:uniaxial_model_micro_station}Uniaxial model of the micro-station}
\end{subfigure}
\caption{\label{fig:micro_station_uniaxial_model}Schematic of the Micro-Station measurement setup and uniaxial model.}
\end{figure}

\section{Uniaxial Model}
The uniaxial model of the micro-station is shown in Figure \ref{fig:uniaxial_model_micro_station}.
It consists of a mass spring damper system with three degrees of freedom.
A mass-spring-damper system represents the granite (with mass \(m_g\), stiffness \(k_g\) and damping \(c_g\)).
Another mass-spring-damper system represents the different micro-station stages (the \(T_y\) stage, the \(R_y\) stage and the \(R_z\) stage) with mass \(m_t\), damping \(c_t\) and stiffness \(k_t\).
Finally, a third mass-spring-damper system represents the micro-hexapod with mass \(m_h\), damping \(c_h\) and stiffness \(k_h\).

The masses of the different stages are estimated from the 3D model, while the stiffnesses are from the data-sheet of the manufacturers.
The damping coefficients were tuned to match the damping identified from the measurements.
The parameters obtained are summarized in Table \ref{tab:uniaxial_ustation_parameters}.

\begin{table}[htbp]
\centering
\begin{tabularx}{0.9\linewidth}{lXXX}
\toprule
\textbf{Stage} & \textbf{Mass} & \textbf{Stiffness} & \textbf{Damping}\\
\midrule
Micro-Hexapod & \(m_h = 15\,\text{kg}\) & \(k_h = 61\,N/\mu m\) & \(c_h = 3\,\frac{kN}{m/s}\)\\
\(T_y\), \(R_y\), \(R_z\) & \(m_t = 1200\,\text{kg}\) & \(k_t = 520\,N/\mu m\) & \(c_t = 80\,\frac{kN}{m/s}\)\\
Granite & \(m_g = 2500\,\text{kg}\) & \(k_g = 950\,N/\mu m\) & \(c_g = 250\,\frac{kN}{m/s}\)\\
\bottomrule
\end{tabularx}
\caption{\label{tab:uniaxial_ustation_parameters}Physical parameters used for the micro-station uniaxial model}

\end{table}

Two disturbances are considered which are shown in red: the floor motion \(x_f\) and the stage vibrations represented by \(f_t\).
The hammer impacts \(F_{h}, F_{g}\) are shown in blue, whereas the measured inertial motions \(x_{h}, x_{g}\) are shown in black.

\section{Comparison of model and measurements}
The transfer functions from the forces injected by the hammers to the measured inertial motion of the micro-hexapod and granite are extracted from the uniaxial model and compared to the measurements in Figure \ref{fig:uniaxial_comp_frf_meas_model}.

Because the uniaxial model has three degrees of freedom, only three modes with frequencies at \(70\,\text{Hz}\), \(140\,\text{Hz}\) and \(320\,\text{Hz}\) are modeled.
Many more modes can be observed in the measurements (see Figure \ref{fig:uniaxial_comp_frf_meas_model}).
However, the goal is not to have a perfect match with the measurement (this would require a much more complex model), but to have a first approximation.
More accurate models will be used later on.

\begin{figure}[htbp]
\centering
\includegraphics[scale=1]{figs/uniaxial_comp_frf_meas_model.png}
\caption{\label{fig:uniaxial_comp_frf_meas_model}Comparison of the measured FRF and identified ones from the uniaxial model}
\end{figure}

\chapter{Nano-Hexapod Model}
\label{sec:nano_station_model}
A model of the nano-hexapod and sample is now added on top of the uniaxial model of the micro-station (Figure \ref{fig:uniaxial_model_micro_station_nass}).
Disturbances (shown in red) are \gls{fs} the direct forces applied to the sample (for example cable forces), \gls{ft} representing the vibrations induced when scanning the different stages and \gls{xf} the floor motion.
The control signal is the force applied by the nano-hexapod \(f\) and the measurement is the relative motion between the sample and the granite \(d\).
The sample is here considered as a rigid body and rigidly fixed to the nano-hexapod.
The effect of resonances between the sample's point of interest and the nano-hexapod actuator will be considered in Section \ref{sec:uniaxial_payload_dynamics}.

\begin{figure}[htbp]
\begin{subfigure}{0.39\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_model_micro_station_nass.png}
\end{center}
\subcaption{\label{fig:uniaxial_model_micro_station_nass}Uniaxial mass-spring-damper model of the NASS}
\end{subfigure}
\begin{subfigure}{0.59\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_plant_first_params.png}
\end{center}
\subcaption{\label{fig:uniaxial_plant_first_params}Bode Plot of the transfer function from actuator forces $f$ to measured displacement $d$ by the metrology}
\end{subfigure}
\caption{\label{fig:uniaxial_model_micro_station_nass_with_tf}Uniaxial model of the NASS (\subref{fig:uniaxial_model_micro_station_nass}) with the micro-station shown in black, the nano-hexapod represented in blue and the sample represented in green. Disturbances are shown in red. Extracted transfer function from \(f\) to \(d\) (\subref{fig:uniaxial_plant_first_params}).}
\end{figure}
\section{Nano-Hexapod Parameters}
The nano-hexapod is represented by a mass spring damper system (shown in blue in Figure \ref{fig:uniaxial_model_micro_station_nass}).
Its mass \gls{mn} is set to \(15\,\text{kg}\) while its stiffness \(k_n\) can vary depending on the chosen architecture/technology.
The sample is represented by a mass \gls{ms} that can vary from \(1\,\text{kg}\) up to \(50\,\text{kg}\).

As a first example, the nano-hexapod stiffness of is set at \(k_n = 10\,N/\mu m\) and the sample mass is chosen at \(m_s = 10\,\text{kg}\).

\section{Obtained Dynamic Response}
The sensitivity to disturbances (i.e., the transfer functions from \(x_f,f_t,f_s\) to \(d\)) can be extracted from the uniaxial model of Figure \ref{fig:uniaxial_model_micro_station_nass} and are shown in Figure \ref{fig:uniaxial_sensitivity_dist_first_params}.
The \emph{plant} (i.e., the transfer function from actuator force \(f\) to measured displacement \(d\)) is shown in Figure \ref{fig:uniaxial_plant_first_params}.

For further analysis, 9 ``configurations'' of the uniaxial NASS model of Figure \ref{fig:uniaxial_model_micro_station_nass} will be considered: three nano-hexapod stiffnesses (\(k_n = 0.01\,N/\mu m\), \(k_n = 1\,N/\mu m\) and \(k_n = 100\,N/\mu m\)) combined with three sample's masses (\(m_s = 1\,kg\), \(m_s = 25\,kg\) and \(m_s = 50\,kg\)).

\begin{figure}[htbp]
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_dist_first_params_fs.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_dist_first_params_fs}Direct forces}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_dist_first_params_ft.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_dist_first_params_ft}$\mu\text{-station}$ disturbances}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_dist_first_params_xf.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_dist_first_params_xf}Floor motion}
\end{subfigure}
\caption{\label{fig:uniaxial_sensitivity_dist_first_params}Sensitivity of the relative motion \(d\) to disturbances: \(f_s\) the direct forces applied on the sample (\subref{fig:uniaxial_sensitivity_dist_first_params_fs}), \(f_t\) disturbances from the micro-station stages (\subref{fig:uniaxial_sensitivity_dist_first_params_ft}) and \(x_f\) the floor motion (\subref{fig:uniaxial_sensitivity_dist_first_params_fs})}
\end{figure}
\chapter{Disturbance Identification}
\label{sec:uniaxial_disturbances}
To quantify disturbances (red signals in Figure \ref{fig:uniaxial_model_micro_station_nass}), three geophones\footnote{Mark Product L-22D geophones are used with a sensitivity of \(88\,\frac{V}{m/s}\) and a natural frequency of \(\approx 2\,\text{Hz}\)} are used.
One is located on the floor, another one on the granite, and the last one on the micro-hexapod's top platform (see Figure \ref{fig:uniaxial_ustation_meas_disturbances}).
The geophone located on the floor was used to measure the floor motion \(x_f\) while the other two geophones were used to measure vibrations introduced by scanning of the \(T_y\) stage and \(R_z\) stage (see Figure \ref{fig:uniaxial_ustation_dynamical_id_setup}).

\begin{figure}[htbp]
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_ustation_meas_disturbances.png}
\end{center}
\subcaption{\label{fig:uniaxial_ustation_meas_disturbances}Disturbance measurement setup - Schematic}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_ustation_dynamical_id_setup.jpg}
\end{center}
\subcaption{\label{fig:uniaxial_ustation_dynamical_id_setup}Two geophones are used to measure vibrations induced by $T_y$ and $R_z$ scans}
\end{subfigure}
\caption{\label{fig:uniaxial_ustation_meas_disturbances_setup}Identification of the disturbances coming from the micro-station. The measurement schematic is shown in (\subref{fig:uniaxial_ustation_meas_disturbances}). A picture of the setup is shown in (\subref{fig:uniaxial_ustation_dynamical_id_setup})}
\end{figure}
\section{Ground Motion}
To acquire the geophone signals, the measurement setup shown in Figure \ref{fig:uniaxial_geophone_meas_chain} is used.
The voltage generated by the geophone is amplified using a low noise voltage amplifier\footnote{DLPVA-100-B from Femto with a voltage input noise is \(2.4\,nV/\sqrt{\text{Hz}}\)} with a gain of 60dB before going to the ADC.
This is done to improve the signal-to-noise ratio.

To reconstruct the displacement \(x_f\) from the measured voltage \(\hat{V}_{x_f}\), the transfer function of the measurement chain from \(x_f\) to \(\hat{V}_{x_f}\) needs to be estimated.
First, the transfer function \(G_{geo}\) from the floor motion \(x_{f}\) to the generated geophone voltage \(V_{x_f}\) is shown in \eqref{eq:uniaxial_geophone_tf}, with \(T_g = 88\,\frac{V}{m/s}\) the sensitivity of the geophone, \(f_0 = \frac{\omega_0}{2\pi} = 2\,\text{Hz}\) its resonance frequency and \(\xi = 0.7\) its damping ratio.
This model of the geophone was taken from \cite{collette12_review}.
The gain of the voltage amplifier is \(V^{\prime}_{x_f}/V_{x_f} = g_0 = 1000\).

\begin{equation}\label{eq:uniaxial_geophone_tf}
G_{geo}(s) = \frac{V_{x_f}}{x_f}(s) = T_{g} \cdot s \cdot \frac{s^2}{s^2 + 2 \xi \omega_0 s + \omega_0^2} \quad \left[ V/m \right]
\end{equation}

\begin{figure}[htbp]
\centering
\includegraphics[scale=1]{figs/uniaxial_geophone_meas_chain.png}
\caption{\label{fig:uniaxial_geophone_meas_chain}Measurement setup for one geophone. The inertial displacement \(x\) is converted to a voltage \(V\) by the geophone. This voltage is amplified by a factor \(g_0 = 60\,dB\) using a low-noise voltage amplifier. It is then converted to a digital value \(\hat{V}_x\) using a 16bit ADC.}
\end{figure}

The amplitude spectral density of the floor motion \(\Gamma_{x_f}\) can be computed from the amplitude spectral density of measured voltage \(\Gamma_{\hat{V}_{x_f}}\) using \eqref{eq:uniaxial_asd_floor_motion}.
The estimated amplitude spectral density \(\Gamma_{x_f}\) of the floor motion \(x_f\) is shown in Figure \ref{fig:uniaxial_asd_floor_motion_id31}.

\begin{equation}\label{eq:uniaxial_asd_floor_motion}
\Gamma_{x_f}(\omega) = \frac{\Gamma_{\hat{V}_{x_f}}(\omega)}{|G_{geo}(j\omega)| \cdot g_0} \quad \left[ m/\sqrt{\text{Hz}} \right]
\end{equation}

\begin{figure}[htbp]
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_asd_floor_motion_id31.png}
\end{center}
\subcaption{\label{fig:uniaxial_asd_floor_motion_id31}Estimated ASD of $x_f$}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_asd_disturbance_force.png}
\end{center}
\subcaption{\label{fig:uniaxial_asd_disturbance_force}Estimated ASD of $f_t$}
\end{subfigure}
\caption{\label{fig:uniaxial_asd_disturbance}Estimated amplitude spectral density of the floor motion \(x_f\) (\subref{fig:uniaxial_asd_floor_motion_id31}) and of the stage disturbances \(f_t\) (\subref{fig:uniaxial_asd_disturbance_force})}
\end{figure}

\section{Stage Vibration}
To estimate the vibrations induced by scanning the micro-station stages, two geophones are used, as shown in Figure \ref{fig:uniaxial_ustation_dynamical_id_setup}.
The vertical relative velocity between the top platform of the micro hexapod and the granite is estimated in two cases: without moving the micro-station stages, and then during a Spindle rotation at 6rpm.
The vibrations induced by the \(T_y\) stage are not considered here because they have less amplitude than the vibrations induced by the \(R_z\) stage and because the \(T_y\) stage can be scanned at lower velocities if the induced vibrations are found to be an issue.

The amplitude spectral density of the relative motion with and without the Spindle rotation are compared in Figure \ref{fig:uniaxial_asd_vibration_spindle_rotation}.
It is shown that the spindle rotation increases the vibrations above \(20\,\text{Hz}\).
The sharp peak observed at \(24\,\text{Hz}\) is believed to be induced by electromagnetic interference between the currents in the spindle motor phases and the geophone cable because this peak is not observed when rotating the spindle ``by hand''.

\begin{figure}[htbp]
\centering
\includegraphics[scale=1]{figs/uniaxial_asd_vibration_spindle_rotation.png}
\caption{\label{fig:uniaxial_asd_vibration_spindle_rotation}Amplitude Spectral Density \(\Gamma_{R_z}\) of the relative motion measured between the granite and the micro-hexapod's top platform during Spindle rotating}
\end{figure}

To compute the equivalent disturbance force \(f_t\) (Figure \ref{fig:uniaxial_model_micro_station}) that induces such motion, the transfer function \(G_{f_t}(s)\) from \(f_t\) to the relative motion between the micro-hexapod's top platform and the granite \((x_{h} - x_{g})\) is extracted from the model.
The amplitude spectral density \(\Gamma_{f_{t}}\) of the disturbance force is them computed from \eqref{eq:uniaxial_ft_asd} and is shown in Figure \ref{fig:uniaxial_asd_disturbance_force}.

\begin{equation}\label{eq:uniaxial_ft_asd}
\Gamma_{f_{t}}(\omega) = \frac{\Gamma_{R_{z}}(\omega)}{|G_{f_t}(j\omega)|}
\end{equation}

\chapter{Open-Loop Dynamic Noise Budgeting}
\label{sec:uniaxial_noise_budgeting}
Now that a model of the \acrshort{nass} has been obtained (see section \ref{sec:nano_station_model}) and that the disturbances have been estimated (see section \ref{sec:uniaxial_disturbances}), it is possible to perform an \emph{open-loop dynamic noise budgeting}.

To perform such noise budgeting, the disturbances need to be modeled by their spectral densities (done in section \ref{sec:uniaxial_disturbances}).
Then, the transfer functions from disturbances to the performance metric (here the distance \(d\)) are computed (Section \ref{ssec:uniaxial_noise_budget_sensitivity}).
Finally, these two types of information are combined to estimate the corresponding spectral density of the performance metric.
This is very useful to identify what is limiting the performance of the system, or the compare the achievable performance with different system parameters (Section \ref{ssec:uniaxial_noise_budget_result}).
\section{Sensitivity to disturbances}
\label{ssec:uniaxial_noise_budget_sensitivity}
From the uniaxial model of the \acrshort{nass} (Figure \ref{fig:uniaxial_model_micro_station_nass}), the transfer function from the disturbances (\(f_s\), \(x_f\) and \(f_t\)) to the displacement \(d\) are computed.

This is done for two extreme sample masses \(m_s = 1\,\text{kg}\) and \(m_s = 50\,\text{kg}\) and three nano-hexapod stiffnesses:
\begin{itemize}
\item \(k_n = 0.01\,N/\mu m\) that represents a voice coil actuator with soft flexible guiding
\item \(k_n = 1\,N/\mu m\) that represents a voice coil actuator with a stiff flexible guiding or a mechanically amplified piezoelectric actuator
\item \(k_n = 100\,N/\mu m\) that represents a stiff piezoelectric stack actuator
\end{itemize}

The obtained sensitivity to disturbances for the three nano-hexapod stiffnesses are shown in Figure \ref{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses} for the sample mass \(m_s = 1\,\text{kg}\) (the same conclusions can be drawn with \(m_s = 50\,\text{kg}\)):
\begin{itemize}
\item The soft nano-hexapod is more sensitive to forces applied on the sample (cable forces for instance), which is expected due to its lower stiffness (Figure \ref{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_fs})
\item Between the suspension mode of the nano-hexapod (here at 5Hz for the soft nano-hexapod) and the first mode of the micro-station (here at 70Hz), the disturbances induced by the stage vibrations are filtered out (Figure \ref{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_ft})
\item Above the suspension mode of the nano-hexapod, the sample's inertial motion is unaffected by the floor motion; therefore, the sensitivity to floor motion is close to \(1\) (Figure \ref{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_xf})
\end{itemize}

\begin{figure}[htbp]
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_fs.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_fs}Direct forces}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_ft.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_ft}$\mu\text{-station}$ disturbances}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_xf.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_xf}Floor motion}
\end{subfigure}
\caption{\label{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses}Sensitivity of \(d\) to disturbances for three different nano-hexpod stiffnesses. \(f_s\) the direct forces applied on the sample (\subref{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_fs}), \(f_t\) disturbances from the micro-station stages (\subref{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_ft}) and \(x_f\) the floor motion (\subref{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_fs})}
\end{figure}

\section{Open-Loop Dynamic Noise Budgeting}
\label{ssec:uniaxial_noise_budget_result}
Now, the amplitude spectral densities of the disturbances are considered to estimate the residual motion \(d\) for each nano-hexapod and sample configuration.
The Cumulative Amplitude Spectrum of the relative motion \(d\) due to both floor motion \(x_f\) and stage vibrations \(f_t\) are shown in Figure \ref{fig:uniaxial_cas_d_disturbances_stiffnesses} for the three nano-hexapod stiffnesses.
It is shown that the effect of floor motion is much less than that of stage vibrations, except for the soft nano-hexapod below \(5\,\text{Hz}\).

The total cumulative amplitude spectrum of \(d\) for the three nano-hexapod stiffnesses and for the two samples masses are shown in Figure \ref{fig:uniaxial_cas_d_disturbances_payload_masses}.
The conclusion is that the sample mass has little effect on the cumulative amplitude spectrum of the relative motion \(d\).

\begin{figure}[htbp]
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_cas_d_disturbances_stiffnesses.png}
\end{center}
\subcaption{\label{fig:uniaxial_cas_d_disturbances_stiffnesses}Effect of floor motion $x_f$ and stage disturbances $f_t$}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_cas_d_disturbances_payload_masses.png}
\end{center}
\subcaption{\label{fig:uniaxial_cas_d_disturbances_payload_masses}Effect of nano-hexapod stiffness $k_n$ and payload mass $m_s$}
\end{subfigure}
\caption{\label{fig:uniaxial_cas_d_disturbances}Cumulative Amplitude Spectrum of the relative motion \(d\). The effect of \(x_f\) and \(f_t\) are shown in (\subref{fig:uniaxial_cas_d_disturbances_stiffnesses}). The effect of sample mass for the three hexapod stiffnesses is shown in (\subref{fig:uniaxial_cas_d_disturbances_payload_masses}). The control objective of having a residual error of 20 nm RMS is shown by the horizontal black dashed line.}
\end{figure}

\section*{Conclusion}
The open-loop residual vibrations of \(d\) can be estimated from the low-frequency value of the cumulative amplitude spectrum in Figure \ref{fig:uniaxial_cas_d_disturbances_payload_masses}.
This residual vibration of \(d\) is found to be in the order of \(100\,nm\,\text{RMS}\) for the stiff nano-hexapod (\(k_n = 100\,N/\mu m\)), \(200\,nm\,\text{RMS}\) for the relatively stiff nano-hexapod (\(k_n = 1\,N/\mu m\)) and  \(1\,\mu m\,\text{RMS}\) for the soft nano-hexapod (\(k_n = 0.01\,N/\mu m\)).
From this analysis, it may be concluded that the stiffer the nano-hexapod the better.

However, what is more important is the \emph{closed-loop} residual vibration of \(d\) (i.e., while the feedback controller is used).
The goal is to obtain a closed-loop residual vibration \(\epsilon_d \approx 20\,nm\,\text{RMS}\) (represented by an horizontal dashed black line in Figure \ref{fig:uniaxial_cas_d_disturbances_payload_masses}).
The bandwidth of the feedback controller leading to a closed-loop residual vibration of \(20\,nm\,\text{RMS}\) can be estimated as the frequency at which the cumulative amplitude spectrum crosses the black dashed line in Figure \ref{fig:uniaxial_cas_d_disturbances_payload_masses}.
A closed loop bandwidth of \(\approx 10\,\text{Hz}\) is found for the soft nano-hexapod (\(k_n = 0.01\,N/\mu m\)), \(\approx 50\,\text{Hz}\) for the relatively stiff nano-hexapod (\(k_n = 1\,N/\mu m\)), and \(\approx 100\,\text{Hz}\) for the stiff nano-hexapod (\(k_n = 100\,N/\mu m\)).
Therefore, while the \emph{open-loop} vibration is the lowest for the stiff nano-hexapod, it requires the largest feedback bandwidth to meet the specifications.

The advantage of the soft nano-hexapod can be explained by its natural isolation from the micro-station vibration above its suspension mode, as shown in Figure \ref{fig:uniaxial_sensitivity_disturbances_nano_hexapod_stiffnesses_ft}.

\chapter{Active Damping}
\label{sec:uniaxial_active_damping}
In this section, three active damping techniques are applied to the nano-hexapod (see Figure \ref{fig:uniaxial_active_damping_strategies}): Integral Force Feedback (IFF) \cite{preumont91_activ}, Relative Damping Control (RDC) \cite[Chapter 7.2]{preumont18_vibrat_contr_activ_struc_fourt_edition} and Direct Velocity Feedback (DVF) \cite{karnopp74_vibrat_contr_using_semi_activ_force_gener,serrand00_multic_feedb_contr_isolat_base_excit_vibrat,preumont02_force_feedb_versus_accel_feedb}.

These damping strategies are first described (Section \ref{ssec:uniaxial_active_damping_strategies}) and are then compared in terms of achievable damping of the nano-hexapod mode (Section \ref{ssec:uniaxial_active_damping_achievable_damping}), reduction of the effect of disturbances (i.e., \(x_f\), \(f_t\) and \(f_s\)) on the displacement \(d\) (Sections \ref{ssec:uniaxial_active_damping_sensitivity_disturbances}).

\begin{figure}[htbp]
\begin{subfigure}{0.37\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_active_damping_strategies_iff.png}
\end{center}
\subcaption{\label{fig:uniaxial_active_damping_strategies_iff}IFF}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_active_damping_strategies_rdc.png}
\end{center}
\subcaption{\label{fig:uniaxial_active_damping_strategies_rdc}RDC}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_active_damping_strategies_dvf.png}
\end{center}
\subcaption{\label{fig:uniaxial_active_damping_strategies_dvf}DVF}
\end{subfigure}
\caption{\label{fig:uniaxial_active_damping_strategies}Three active damping strategies. Integral Force Feedback (\subref{fig:uniaxial_active_damping_strategies_iff}) using a force sensor, Relative Damping Control (\subref{fig:uniaxial_active_damping_strategies_rdc}) using a relative displacement sensor, and Direct Velocity Feedback (\subref{fig:uniaxial_active_damping_strategies_dvf}) using a geophone}
\end{figure}
\section{Active Damping Strategies}
\label{ssec:uniaxial_active_damping_strategies}
\paragraph{Integral Force Feedback (IFF)}
The Integral Force Feedback strategy consists of using a force sensor in series with the actuator (see Figure \ref{fig:uniaxial_active_damping_iff_schematic}) and applying an ``integral'' feedback controller \eqref{eq:uniaxial_iff_controller}.

\begin{equation}\label{eq:uniaxial_iff_controller}
\boxed{K_{\text{IFF}}(s) = \frac{g}{s}}
\end{equation}

The mechanical equivalent of this IFF strategy is a dashpot in series with the actuator stiffness with a damping coefficient equal to the stiffness of the actuator divided by the controller gain \(k/g\) (see Figure \ref{fig:uniaxial_active_damping_iff_equiv}).

\begin{figure}[htbp]
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_active_damping_iff_schematic.png}
\end{center}
\subcaption{\label{fig:uniaxial_active_damping_iff_schematic}Integral Force Feedback}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_active_damping_iff_equiv.png}
\end{center}
\subcaption{\label{fig:uniaxial_active_damping_iff_equiv}Equivalent mechanical representation}
\end{subfigure}
\caption{\label{fig:uniaxial_active_damping_iff}Integral Force Feedback (\subref{fig:uniaxial_active_damping_iff_schematic}) is equivalent to a damper in series with the actuator stiffness (\subref{fig:uniaxial_active_damping_iff_equiv})}
\end{figure}

\paragraph{Relative Damping Control (RDC)}
For the Relative Damping Control strategy, a relative motion sensor that measures the motion of the actuator is used (see Figure \ref{fig:uniaxial_active_damping_rdc_schematic}) and a ``derivative'' feedback controller is used \eqref{eq:uniaxial_rdc_controller}.

\begin{equation}\label{eq:uniaxial_rdc_controller}
\boxed{K_{\text{RDC}}(s) = - g \cdot s}
\end{equation}

The mechanical equivalent of \acrshort{rdc} is a dashpot in parallel with the actuator with a damping coefficient equal to the controller gain \(g\) (see Figure \ref{fig:uniaxial_active_damping_rdc_equiv}).

\begin{figure}[htbp]
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_active_damping_rdc_schematic.png}
\end{center}
\subcaption{\label{fig:uniaxial_active_damping_rdc_schematic}Relative motion control}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_active_damping_rdc_equiv.png}
\end{center}
\subcaption{\label{fig:uniaxial_active_damping_rdc_equiv}Equivalent mechanical representation}
\end{subfigure}
\caption{\label{fig:uniaxial_active_damping_rdc}Relative Damping Control (\subref{fig:uniaxial_active_damping_rdc_schematic}) is equivalent to a damper in parallel with the actuator (\subref{fig:uniaxial_active_damping_rdc_equiv})}
\end{figure}

\paragraph{Direct Velocity Feedback (DVF)}
Finally, the direct velocity feedback strategy consists of using an inertial sensor (usually a geophone) that measures the ``absolute'' velocity of the body fixed on top of the actuator (see Figure \ref{fig:uniaxial_active_damping_dvf_schematic}).
This velocity is fed back to the actuator with a ``proportional'' controller \eqref{eq:uniaxial_dvf_controller}.

\begin{equation}\label{eq:uniaxial_dvf_controller}
\boxed{K_{\text{DVF}}(s) = - g}
\end{equation}

This is equivalent to a dashpot (with a damping coefficient equal to the controller gain \(g\)) between the body (on which the inertial sensor is fixed) and an inertial reference frame (see Figure \ref{fig:uniaxial_active_damping_dvf_equiv}).
This is usually referred to as ``\emph{sky hook damper}''.

\begin{figure}[htbp]
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_active_damping_dvf_schematic.png}
\end{center}
\subcaption{\label{fig:uniaxial_active_damping_dvf_schematic}Direct velocity feedback}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_active_damping_dvf_equiv.png}
\end{center}
\subcaption{\label{fig:uniaxial_active_damping_dvf_equiv}Equivalent mechanical representation}
\end{subfigure}
\caption{\label{fig:uniaxial_active_damping_dvf}Direct velocity Feedback (\subref{fig:uniaxial_active_damping_dvf_schematic}) is equivalent to a ``sky hook damper'' (\subref{fig:uniaxial_active_damping_dvf_equiv})}
\end{figure}

\section{Plant Dynamics for Active Damping}
\label{ssec:uniaxial_active_damping_plants}
The plant dynamics for all three active damping techniques are shown in Figure \ref{fig:uniaxial_plant_active_damping_techniques}.
All have \emph{alternating poles and zeros} meaning that the phase does not vary by more than \(180\,\text{deg}\) which makes the design of a \emph{robust} damping controller very easy.

This alternating poles and zeros property is guaranteed for the IFF and RDC cases because the sensors are collocated with the actuator \cite[Chapter 7]{preumont18_vibrat_contr_activ_struc_fourt_edition}.
For the DVF controller, this property is not guaranteed, and may be lost if some flexibility between the nano-hexapod and the sample is considered \cite[Chapter 8.4]{preumont18_vibrat_contr_activ_struc_fourt_edition}.

When the nano-hexapod's suspension modes are at frequencies lower than the resonances of the micro-station (blue and red curves in Figure \ref{fig:uniaxial_plant_active_damping_techniques}), the resonances of the micro-stations have little impact on the IFF and DVF transfer functions.
For the stiff nano-hexapod (yellow curves), the micro-station dynamics can be seen on the transfer functions in Figure \ref{fig:uniaxial_plant_active_damping_techniques}.
Therefore, it is expected that the micro-station dynamics might impact the achievable damping if a stiff nano-hexapod is used.

\begin{figure}[htbp]
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.99\linewidth]{figs/uniaxial_plant_active_damping_techniques_iff.png}
\end{center}
\subcaption{\label{fig:uniaxial_plant_active_damping_techniques_iff}IFF}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.99\linewidth]{figs/uniaxial_plant_active_damping_techniques_rdc.png}
\end{center}
\subcaption{\label{fig:uniaxial_plant_active_damping_techniques_rdc}RDC}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.99\linewidth]{figs/uniaxial_plant_active_damping_techniques_dvf.png}
\end{center}
\subcaption{\label{fig:uniaxial_plant_active_damping_techniques_dvf}DVF}
\end{subfigure}
\caption{\label{fig:uniaxial_plant_active_damping_techniques}Plant dynamics for the three active damping techniques (IFF: \subref{fig:uniaxial_plant_active_damping_techniques_iff}, RDC: \subref{fig:uniaxial_plant_active_damping_techniques_rdc}, DVF: \subref{fig:uniaxial_plant_active_damping_techniques_dvf}), for three nano-hexapod stiffnesses (\(k_n = 0.01\,N/\mu m\) in blue, \(k_n = 1\,N/\mu m\) in red and \(k_n = 100\,N/\mu m\) in yellow) and three sample's masses (\(m_s = 1\,kg\): solid curves, \(m_s = 25\,kg\): dot-dashed curves, and \(m_s = 50\,kg\): dashed curves).}
\end{figure}

\section{Achievable Damping and Damped Plants}
\label{ssec:uniaxial_active_damping_achievable_damping}

To compare the added damping using the three considered active damping strategies, the root locus plot is used.
Indeed, the damping ratio \(\xi\) of a pole in the complex plane can be estimated from the angle \(\phi\) it makes with the imaginary axis \eqref{eq:uniaxial_damping_ratio_angle}.
Increasing the angle with the imaginary axis therefore means that more damping is added to the considered resonance.
This is illustrated in Figure \ref{fig:uniaxial_root_locus_damping_techniques_micro_station_mode} by the dashed black line indicating the maximum achievable damping.

\begin{equation}\label{eq:uniaxial_damping_ratio_angle}
\xi = \sin(\phi)
\end{equation}

The Root Locus for the three nano-hexapod stiffnesses and the three active damping techniques are shown in Figure \ref{fig:uniaxial_root_locus_damping_techniques}.
All three active damping approaches can lead to \emph{critical damping} of the nano-hexapod suspension mode (angle \(\phi\) can be increased up to 90 degrees).
There is even some damping authority on micro-station modes in the following cases:
\begin{description}
\item[{IFF with a stiff nano-hexapod (Figure \ref{fig:uniaxial_root_locus_damping_techniques_stiff})}] This can be understood from the mechanical equivalent of IFF shown in Figure \ref{fig:uniaxial_active_damping_iff_equiv} considering an high stiffness \(k\).
The micro-station top platform is connected to an inertial mass (the nano-hexapod) through a damper, which dampens the micro-station suspension suspension mode.
\item[{DVF with a stiff nano-hexapod (Figure \ref{fig:uniaxial_root_locus_damping_techniques_stiff})}] In that case, the ``sky hook damper'' (see mechanical equivalent of DVF in Figure \ref{fig:uniaxial_active_damping_dvf_equiv}) is connected to the micro-station top platform through the stiff nano-hexapod.
\item[{RDC with a soft nano-hexapod (Figure \ref{fig:uniaxial_root_locus_damping_techniques_micro_station_mode})}] At the frequency of the micro-station mode, the nano-hexapod top mass behaves as an inertial reference because the suspension mode of the soft nano-hexapod is at much lower frequency.
The micro-station and the nano-hexapod masses are connected through a large damper induced by RDC (see mechanical equivalent in Figure \ref{fig:uniaxial_active_damping_rdc_equiv}) which allows some damping of the micro-station.
\end{description}

\begin{figure}[htbp]
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_root_locus_damping_techniques_soft.png}
\end{center}
\subcaption{\label{fig:uniaxial_root_locus_damping_techniques_soft}$k_n = 0.01\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_root_locus_damping_techniques_mid.png}
\end{center}
\subcaption{\label{fig:uniaxial_root_locus_damping_techniques_mid}$k_n = 1\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_root_locus_damping_techniques_stiff.png}
\end{center}
\subcaption{\label{fig:uniaxial_root_locus_damping_techniques_stiff}$k_n = 100\,N/\mu m$}
\end{subfigure}
\caption{\label{fig:uniaxial_root_locus_damping_techniques}Root Loci for the three active damping techniques (IFF in blue, RDC in red and DVF in yellow). This is shown for the three nano-hexapod stiffnesses. The Root Loci are zoomed in the suspension mode of the nano-hexapod.}
\end{figure}

\begin{figure}[htbp]
\centering
\includegraphics[scale=1]{figs/uniaxial_root_locus_damping_techniques_micro_station_mode.png}
\caption{\label{fig:uniaxial_root_locus_damping_techniques_micro_station_mode}Root Locus for the three damping techniques applied with the soft nano-hexapod. It is shown that the RDC active damping technique has some authority on one mode of the micro-station. This mode corresponds to the suspension mode of the micro-hexapod.}
\end{figure}

The transfer functions from the plant input \(f\) to the relative displacement \(d\) while active damping is implemented are shown in Figure \ref{fig:uniaxial_damped_plant_three_active_damping_techniques}.
All three active damping techniques yielded similar damped plants.

\begin{figure}[htbp]
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_damped_plant_three_active_damping_techniques_vc.png}
\end{center}
\subcaption{\label{fig:uniaxial_damped_plant_three_active_damping_techniques_vc}$k_n = 0.01\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_damped_plant_three_active_damping_techniques_md.png}
\end{center}
\subcaption{\label{fig:uniaxial_damped_plant_three_active_damping_techniques_md}$k_n = 1\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_damped_plant_three_active_damping_techniques_pz.png}
\end{center}
\subcaption{\label{fig:uniaxial_damped_plant_three_active_damping_techniques_pz}$k_n = 100\,N/\mu m$}
\end{subfigure}
\caption{\label{fig:uniaxial_damped_plant_three_active_damping_techniques}Obtained damped transfer function from \(f\) to \(d\) for the three damping techniques.}
\end{figure}

\section{Sensitivity to disturbances and Noise Budgeting}
\label{ssec:uniaxial_active_damping_sensitivity_disturbances}

Reasonable gains are chosen for the three active damping strategies such that the nano-hexapod suspension mode is well damped.
The sensitivity to disturbances (direct forces \(f_s\), stage vibrations \(f_t\) and floor motion \(x_f\)) for all three active damping techniques are compared in Figure \ref{fig:uniaxial_sensitivity_dist_active_damping}.
The comparison is done with the nano-hexapod having a stiffness \(k_n = 1\,N/\mu m\).

Several conclusions can be drawn by comparing the obtained sensitivity transfer functions:
\begin{itemize}
\item IFF degrades the sensitivity to direct forces on the sample (i.e., the compliance) below the resonance of the nano-hexapod (Figure \ref{fig:uniaxial_sensitivity_dist_active_damping_fs}).
This is a well-known effect of using IFF for vibration isolation \cite{collette15_sensor_fusion_method_high_perfor}.
\item RDC degrades the sensitivity to stage vibrations around the nano-hexapod's resonance as compared to the other two methods (Figure \ref{fig:uniaxial_sensitivity_dist_active_damping_ft}).
This is because the equivalent damper in parallel with the actuator (see Figure \ref{fig:uniaxial_active_damping_rdc_equiv}) increases the transmission of the micro-station vibration to the sample which is not the same for the other two active damping strategies.
\item both IFF and DVF degrade the sensitivity to floor motion below the resonance of the nano-hexapod (Figure \ref{fig:uniaxial_sensitivity_dist_active_damping_xf}).
\end{itemize}

\begin{figure}[htbp]
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_dist_active_damping_fs.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_dist_active_damping_fs}Direct forces}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_dist_active_damping_ft.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_dist_active_damping_ft}$\mu\text{-station}$ disturbances}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_dist_active_damping_xf.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_dist_active_damping_xf}Floor motion}
\end{subfigure}
\caption{\label{fig:uniaxial_sensitivity_dist_active_damping}Change of sensitivity to disturbance with all three active damping strategies. \(f_s\) the direct forces applied on the sample (\subref{fig:uniaxial_sensitivity_dist_active_damping_fs}), \(f_t\) disturbances from the micro-station stages (\subref{fig:uniaxial_sensitivity_dist_active_damping_ft}) and \(x_f\) the floor motion (\subref{fig:uniaxial_sensitivity_dist_active_damping_fs})}
\end{figure}

From the amplitude spectral density of the disturbances (computed in Section \ref{sec:uniaxial_disturbances}) and the sensitivity to disturbances estimated using the three active damping strategies, a noise budget can be calculated.
The cumulative amplitude spectrum of the distance \(d\) with all three active damping techniques is shown in Figure \ref{fig:uniaxial_cas_active_damping} and compared with the open-loop case.
All three active damping methods give similar results.

\begin{figure}[htbp]
\begin{subfigure}{0.37\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_cas_active_damping_soft.png}
\end{center}
\subcaption{\label{fig:uniaxial_cas_active_damping_soft}$k_n = 0.01\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_cas_active_damping_mid.png}
\end{center}
\subcaption{\label{fig:uniaxial_cas_active_damping_mid}$k_n = 1\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_cas_active_damping_stiff.png}
\end{center}
\subcaption{\label{fig:uniaxial_cas_active_damping_stiff}$k_n = 100\,N/\mu m$}
\end{subfigure}
\caption{\label{fig:uniaxial_cas_active_damping}Comparison of the cumulative amplitude spectrum (CAS) of the distance \(d\) for all three active damping techniques (OL in black, IFF in blue, RDC in red and DVF in yellow).}
\end{figure}

\section*{Conclusion}
Three active damping strategies have been studied for the \acrfull{nass}.
Equivalent mechanical representations were derived in Section \ref{ssec:uniaxial_active_damping_strategies} which are helpful for understanding the specific effects of each strategy.
The plant dynamics were then compared in Section \ref{ssec:uniaxial_active_damping_plants} and were found to all have alternating poles and zeros, which helps in the design of the active damping controller.
However, this property is not guaranteed for DVF.
The achievable damping of the nano-hexapod suspension mode can be made as large as possible for all three active damping techniques (Section \ref{ssec:uniaxial_active_damping_achievable_damping}).
Even some damping can be applied to some micro-station modes in specific cases.
The obtained damped plants were found to be similar.
The damping strategies were then compared in terms of disturbance reduction in Section \ref{ssec:uniaxial_active_damping_sensitivity_disturbances}.

The comparison between the three active damping strategies is summarized in Table \ref{tab:comp_active_damping}.
It is difficult to conclude on the best active damping strategy for the \acrfull{nass} yet.
Which one will be used will be determined by the use of more accurate models and will depend on which is the easiest to implement in practice

\begin{table}[htbp]
\centering
\scriptsize
\begin{tabularx}{0.9\linewidth}{Xccc}
\toprule
 & \textbf{IFF} & \textbf{RDC} & \textbf{DVF}\\
\midrule
\textbf{Sensor} & Force sensor & Relative motion sensor & Inertial sensor\\
\midrule
\textbf{Damping} & Up to critical & Up to critical & Up to Critical\\
\midrule
\textbf{Robustness} & Requires collocation & Requires collocation & Impacted by geophone resonances\\
\midrule
\(f_s\) \textbf{Disturbance} & \(\nearrow\) at low frequency & \(\searrow\) near resonance & \(\searrow\) near resonance\\
\(f_t\) \textbf{Disturbance} & \(\searrow\) near resonance & \(\nearrow\) near resonance & \(\searrow\) near resonance\\
\(x_f\) \textbf{Disturbance} & \(\nearrow\) at low frequency & \(\searrow\) near resonance & \(\nearrow\) at low frequency\\
\bottomrule
\end{tabularx}
\caption{\label{tab:comp_active_damping}Comparison of active damping strategies}

\end{table}

\chapter{Position Feedback Controller}
\label{sec:uniaxial_position_control}
The \gls{haclac} architecture is shown in Figure \ref{fig:uniaxial_hac_lac_architecture}.
This corresponds to a \emph{two step} control strategy:
\begin{itemize}
\item First, an active damping controller \(\bm{K}_{\textsc{LAC}}\) is implemented (see Section \ref{sec:uniaxial_active_damping}).
It allows the vibration level to be reduced, and it also makes the damped plant (transfer function from \(u^{\prime}\) to \(y\)) easier to control than the undamped plant (transfer function from \(u\) to \(y\)).
This is called \emph{low authority} control as it only slightly affects the system poles \cite[Chapter 14.6]{preumont18_vibrat_contr_activ_struc_fourt_edition}.
\item Then, a position controller \(\bm{K}_{\textsc{HAC}}\) is implemented and is used to control the position \(d\).
This is called \emph{high authority} control as it usually relocates the system's poles.
\end{itemize}

In this section, Integral Force Feedback is used as the Low Authority Controller (the other two damping strategies would lead to the same conclusions here).
This control architecture applied to the uniaxial model is shown in Figure \ref{fig:uniaxial_hac_lac_model}.

\begin{figure}[htbp]
\begin{subfigure}{0.54\textwidth}
\begin{center}
\includegraphics[scale=1,width=1.0\linewidth]{figs/uniaxial_hac_lac_architecture.png}
\end{center}
\subcaption{\label{fig:uniaxial_hac_lac_architecture}Typical HAC-LAC Architecture}
\end{subfigure}
\begin{subfigure}{0.45\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_hac_lac_model.png}
\end{center}
\subcaption{\label{fig:uniaxial_hac_lac_model}Uniaxial model with HAC-IFF strategy}
\end{subfigure}
\caption{\label{fig:uniaxial_hac_lac}\acrfull{haclac}}
\end{figure}
\section{Damped Plant Dynamics}
\label{ssec:uniaxial_position_control_damped_dynamics}
The damped plants obtained for the three nano-hexapod stiffnesses are shown in Figure \ref{fig:uniaxial_hac_iff_damped_plants_masses}.
For \(k_n = 0.01\,N/\mu m\) and \(k_n = 1\,N/\mu m\), the dynamics are quite simple and can be well approximated by a second-order plant (Figures \ref{fig:uniaxial_hac_iff_damped_plants_masses_soft} and \ref{fig:uniaxial_hac_iff_damped_plants_masses_mid}).
However, this is not the case for the stiff nano-hexapod (\(k_n = 100\,N/\mu m\)) where two modes can be seen (Figure \ref{fig:uniaxial_hac_iff_damped_plants_masses_stiff}).
This is due to the interaction between the micro-station (modeled modes at 70Hz, 140Hz and 320Hz) and the nano-hexapod.
This effect will be further explained in Section \ref{sec:uniaxial_support_compliance}.

\begin{figure}[htbp]
\begin{subfigure}{0.37\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_hac_iff_damped_plants_masses_soft.png}
\end{center}
\subcaption{\label{fig:uniaxial_hac_iff_damped_plants_masses_soft}$k_n = 0.01\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_hac_iff_damped_plants_masses_mid.png}
\end{center}
\subcaption{\label{fig:uniaxial_hac_iff_damped_plants_masses_mid}$k_n = 1\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_hac_iff_damped_plants_masses_stiff.png}
\end{center}
\subcaption{\label{fig:uniaxial_hac_iff_damped_plants_masses_stiff}$k_n = 100\,N/\mu m$}
\end{subfigure}
\caption{\label{fig:uniaxial_hac_iff_damped_plants_masses}Obtained damped plant using Integral Force Feedback for three sample masses}
\end{figure}

\section{Position Feedback Controller}
\label{ssec:uniaxial_position_control_design}

The objective is to design high-authority feedback controllers for the three nano-hexapods.
This controller must be robust to the change of sample's mass (from \(1\,\text{kg}\) up to \(50\,\text{kg}\)).

The required feedback bandwidths were estimated in Section \ref{sec:uniaxial_noise_budgeting}:
\begin{itemize}
\item \(f_b \approx 10\,\text{Hz}\) for the soft nano-hexapod (\(k_n = 0.01\,N/\mu m\)).
Near this frequency, the plants (shown in Figure \ref{fig:uniaxial_hac_iff_damped_plants_masses_soft}) are equivalent to a mass line (i.e., slope of \(-40\,dB/\text{dec}\) and a phase of -180 degrees).
The gain of this mass line can vary up to a fact \(\approx 5\) (suspended mass from \(16\,kg\) up to \(65\,kg\)).
This means that the designed controller will need to have \emph{large gain margins} to be robust to the change of sample's mass.
\item \(\approx 50\,\text{Hz}\) for the relatively stiff nano-hexapod (\(k_n = 1\,N/\mu m\)).
Similar to the soft nano-hexapod, the plants near the crossover frequency are equivalent to a mass line (Figure \ref{fig:uniaxial_hac_iff_damped_plants_masses_mid}).
It will probably be easier to have a little bit more bandwidth in this configuration to be further away from the nano-hexapod suspension mode.
\item \(\approx 100\,\text{Hz}\) for the stiff nano-hexapod (\(k_n = 100\,N/\mu m\)).
Contrary to the two first nano-hexapod stiffnesses, here the plants have more complex dynamics near the desired crossover frequency (see Figure \ref{fig:uniaxial_hac_iff_damped_plants_masses_stiff}).
The micro-station is not stiff enough to have a clear stiffness line at this frequency.
Therefore, there is both a change of phase and gain depending on the sample mass.
This makes the robust design of the controller more complicated.
\end{itemize}


Position feedback controllers are designed for each nano-hexapod such that it is stable for all considered sample masses with similar stability margins (see Nyquist plots in Figure \ref{fig:uniaxial_nyquist_hac}).
An arbitrary minimum modulus margin of \(0.25\) was chosen when designing the controllers.
These high authority controllers are generally composed of a lag at low frequency for disturbance rejection, a lead to increase the phase margin near the crossover frequency, and a low pass filter to increase the robustness to high frequency dynamics.
The controllers used for the three nano-hexapod are shown in Equation \eqref{eq:uniaxial_hac_formulas}, and the parameters used are summarized in Table \ref{tab:uniaxial_feedback_controller_parameters}.

\begin{subequations} \label{eq:uniaxial_hac_formulas}
\begin{align}
K_{\text{soft}}(s) &= g \cdot
       \underbrace{\frac{s + \omega_0}{s + \omega_i}}_{\text{lag}} \cdot
       \underbrace{\frac{1 + \frac{s}{\omega_c/\sqrt{a}}}{1 + \frac{s}{\omega_c \sqrt{a}}}}_{\text{lead}} \cdot
       \underbrace{\frac{1}{1 + \frac{s}{\omega_l}}}_{\text{LPF}} \\
K_{\text{mid}}(s) &= g \cdot
       \underbrace{\left(\frac{s + \omega_0}{s + \omega_i}\right)^2}_{\text{2 lags}} \cdot
       \underbrace{\frac{1 + \frac{s}{\omega_c/\sqrt{a}}}{1 + \frac{s}{\omega_c \sqrt{a}}}}_{\text{lead}} \cdot
       \underbrace{\frac{1}{1 + \frac{s}{\omega_l}}}_{\text{LPF}} \\
K_{\text{stiff}}(s) &= g \cdot
       \underbrace{\left(\frac{1}{s + \omega_i}\right)^2}_{\text{2 lags}} \cdot
       \underbrace{\left(\frac{1 + \frac{s}{\omega_c/\sqrt{a}}}{1 + \frac{s}{\omega_c \sqrt{a}}}\right)^2}_{\text{2 leads}} \cdot
       \underbrace{\frac{1}{1 + \frac{s}{\omega_l}}}_{\text{LPF}}
\end{align}
\end{subequations}

\begin{table}[htbp]
\centering
\begin{tabularx}{\linewidth}{lXXX}
\toprule
 & \textbf{Soft} & \textbf{Moderately stiff} & \textbf{Stiff}\\
\midrule
\textbf{Gain} & \(g = 4 \cdot 10^5\) & \(g = 3 \cdot 10^6\) & \(g = 6 \cdot 10^{12}\)\\
\textbf{Lead} & \(a = 5\), \(\omega_c = 20\,Hz\) & \(a = 4\), \(\omega_c = 70\,Hz\) & \(a = 5\), \(\omega_c = 100\,Hz\)\\
\textbf{Lag} & \(\omega_0 = 5\,Hz\), \(\omega_i = 0.01\,Hz\) & \(\omega_0 = 20\,Hz\), \(\omega_i = 0.01\,Hz\) & \(\omega_i = 0.01\,Hz\)\\
\textbf{LPF} & \(\omega_l = 200\,Hz\) & \(\omega_l = 300\,Hz\) & \(\omega_l = 500\,Hz\)\\
\bottomrule
\end{tabularx}
\caption{\label{tab:uniaxial_feedback_controller_parameters}Parameters used for the position feedback controllers}

\end{table}

The loop gains corresponding to the designed high authority controllers for the three nano-hexapod are shown in Figure \ref{fig:uniaxial_loop_gain_hac}.
We can see that for the soft and moderately stiff nano-hexapod (Figures \ref{fig:uniaxial_nyquist_hac_vc} and \ref{fig:uniaxial_nyquist_hac_md}), the crossover frequency varies significantly with the sample mass.
This is because the crossover frequency corresponds to the mass line of the plant (whose gain is inversely proportional to the mass).
For the stiff nano-hexapod (Figure \ref{fig:uniaxial_nyquist_hac_pz}), it was difficult to achieve the desired closed-loop bandwidth of \(\approx 100\,\text{Hz}\).
A crossover frequency of \(\approx 65\,\text{Hz}\) was achieved instead.

Note that these controllers were not designed using any optimization methods.
The goal is to have a first estimation of the attainable performance.

\begin{figure}[htbp]
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_nyquist_hac_vc.png}
\end{center}
\subcaption{\label{fig:uniaxial_nyquist_hac_vc}$k_n = 0.01\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_nyquist_hac_md.png}
\end{center}
\subcaption{\label{fig:uniaxial_nyquist_hac_md}$k_n = 1\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_nyquist_hac_pz.png}
\end{center}
\subcaption{\label{fig:uniaxial_nyquist_hac_pz}$k_n = 100\,N/\mu m$}
\end{subfigure}
\caption{\label{fig:uniaxial_nyquist_hac}Nyquist Plot for the high authority controller. The minimum modulus margin is illustrated by a black circle.}
\end{figure}

\begin{figure}[htbp]
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_loop_gain_hac_vc.png}
\end{center}
\subcaption{\label{fig:uniaxial_loop_gain_hac_vc}$k_n = 0.01\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_loop_gain_hac_md.png}
\end{center}
\subcaption{\label{fig:uniaxial_loop_gain_hac_md}$k_n = 1\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_loop_gain_hac_pz.png}
\end{center}
\subcaption{\label{fig:uniaxial_loop_gain_hac_pz}$k_n = 100\,N/\mu m$}
\end{subfigure}
\caption{\label{fig:uniaxial_loop_gain_hac}Loop gain for the High Authority Controller}
\end{figure}

\section{Closed-Loop Noise Budgeting}
\label{ssec:uniaxial_position_control_cl_noise_budget}

The high authority position feedback controllers are then implemented and the closed-loop sensitivities to disturbances are computed.
These are compared with the open-loop and damped plants cases in Figure \ref{fig:uniaxial_sensitivity_dist_hac_lac} for just one configuration (moderately stiff nano-hexapod with 25kg sample's mass).
As expected, the sensitivity to disturbances decreased in the controller bandwidth and slightly increased outside this bandwidth.

\begin{figure}[htbp]
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_dist_hac_lac_fs.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_dist_hac_lac_fs}Direct forces}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_dist_hac_lac_ft.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_dist_hac_lac_ft}$\mu\text{-station}$ disturbances}
\end{subfigure}
\begin{subfigure}{0.33\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_sensitivity_dist_hac_lac_xf.png}
\end{center}
\subcaption{\label{fig:uniaxial_sensitivity_dist_hac_lac_xf}Floor motion}
\end{subfigure}
\caption{\label{fig:uniaxial_sensitivity_dist_hac_lac}Change of sensitivity to disturbances with LAC and with \acrshort{haclac}. A nano-Hexapod with \(k_n = 1\,N/\mu m\) and a sample mass of \(25\,kg\) is used. \(f_s\) the direct forces applied on the sample (\subref{fig:uniaxial_sensitivity_dist_hac_lac_fs}), \(f_t\) disturbances from the micro-station stages (\subref{fig:uniaxial_sensitivity_dist_hac_lac_ft}) and \(x_f\) the floor motion (\subref{fig:uniaxial_sensitivity_dist_hac_lac_fs})}
\end{figure}

The cumulative amplitude spectrum of the motion \(d\) is computed for all nano-hexapod configurations, all sample masses and in the open-loop (OL), damped (IFF) and position controlled (HAC-IFF) cases.
The results are shown in Figure \ref{fig:uniaxial_cas_hac_lac}.
Obtained root mean square values of the distance \(d\) are better for the soft nano-hexapod (\(\approx 25\,nm\) to \(\approx 35\,nm\) depending on the sample's mass) than for the stiffer nano-hexapod (from \(\approx 30\,nm\) to \(\approx 70\,nm\)).

\begin{figure}[htbp]
\begin{subfigure}{0.37\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_cas_hac_lac_soft.png}
\end{center}
\subcaption{\label{fig:uniaxial_cas_hac_lac_soft}$k_n = 0.01\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_cas_hac_lac_mid.png}
\end{center}
\subcaption{\label{fig:uniaxial_cas_hac_lac_mid}$k_n = 1\,N/\mu m$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_cas_hac_lac_stiff.png}
\end{center}
\subcaption{\label{fig:uniaxial_cas_hac_lac_stiff}$k_n = 100\,N/\mu m$}
\end{subfigure}
\caption{\label{fig:uniaxial_cas_hac_lac}Cumulative Amplitude Spectrum for all three nano-hexapod stiffnesses - Comparison of OL, IFF and \acrshort{haclac} cases}
\end{figure}

\section*{Conclusion}
On the basis of the open-loop noise budgeting made in Section \ref{sec:uniaxial_noise_budgeting}, the closed-loop bandwidth required to obtain a vibration level of \(\approx 20\,nm\,\text{RMS}\) was estimated.
To achieve such bandwidth, the \acrshort{haclac} strategy was followed, which consists of first using an active damping controller (studied in Section \ref{sec:uniaxial_active_damping}) and then adding a high authority position feedback controller.

In this section, feedback controllers were designed in such a way that the required closed-loop bandwidth was reached while being robust to changes in the payload mass.
The attainable vibration control performances were estimated for the three nano-hexapod stiffnesses and were found to be close to the required values.
However, the stiff nano-hexapod (\(k_n = 100\,N/\mu m\)) is requiring the largest feedback bandwidth, which is difficult to achieve while being robust to the change of payload mass.
A slight advantage can be given to the soft nano-hexapod as it requires less feedback bandwidth while providing better stability results.

\chapter{Effect of limited micro-station compliance}
\label{sec:uniaxial_support_compliance}
In this section, the impact of the compliance of the support (i.e., the micro-station) on the dynamics of the plant to control is studied.
This is a critical point because the dynamics of the micro-station is complex, depends on the considered direction (see measurements in Figure \ref{fig:uniaxial_comp_frf_meas_model}) and may vary with position and time.
It would be much better to have a plant dynamics that is not impacted by the micro-station.

Therefore, the objective of this section is to obtain some guidance for the design of a nano-hexapod that will not be impacted by the complex micro-station dynamics.
To study this, two models are used (Figure \ref{fig:uniaxial_support_compliance_models}).
The first one consists of the nano-hexapod directly fixed on top of the granite, thus neglecting any support compliance (Figure \ref{fig:uniaxial_support_compliance_nano_hexapod_only}).
The second one consists of the nano-hexapod fixed on top of the micro-station having some limited compliance (Figure \ref{fig:uniaxial_support_compliance_test_system})

\begin{figure}[htbp]
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_support_compliance_nano_hexapod_only.png}
\end{center}
\subcaption{\label{fig:uniaxial_support_compliance_nano_hexapod_only}Nano-Hexapod fixed directly on the Granite}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_support_compliance_test_system.png}
\end{center}
\subcaption{\label{fig:uniaxial_support_compliance_test_system}Nano-Hexapod fixed on top of the Micro-Station}
\end{subfigure}
\caption{\label{fig:uniaxial_support_compliance_models}Models used to study the effect of limited support compliance}
\end{figure}
\section{Neglected support compliance}

The limited compliance of the micro-station is first neglected and the uniaxial model shown in Figure \ref{fig:uniaxial_support_compliance_nano_hexapod_only} is used.
The nano-hexapod mass (including the payload) is set at \(20\,\text{kg}\) and three hexapod stiffnesses are considered, such that their resonance frequencies are at \(\omega_{n} = 10\,\text{Hz}\), \(\omega_{n} = 70\,\text{Hz}\) and \(\omega_{n} = 400\,\text{Hz}\).
Obtained transfer functions from \(F\) to \(L^\prime\) (shown in Figure \ref{fig:uniaxial_effect_support_compliance_neglected}) are simple second-order low-pass filters.
When neglecting the support compliance, a large feedback bandwidth can be achieved for all three nano-hexapods.

\begin{figure}[htbp]
\begin{subfigure}{0.37\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_effect_support_compliance_neglected_soft.png}
\end{center}
\subcaption{\label{fig:uniaxial_effect_support_compliance_neglected_soft}$\omega_{n} \ll \omega_{\mu}$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_effect_support_compliance_neglected_mid.png}
\end{center}
\subcaption{\label{fig:uniaxial_effect_support_compliance_neglected_mid}$\omega_{n} = \omega_{\mu}$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_effect_support_compliance_neglected_stiff.png}
\end{center}
\subcaption{\label{fig:uniaxial_effect_support_compliance_neglected_stiff}$\omega_{n} \gg \omega_{\mu}$}
\end{subfigure}
\caption{\label{fig:uniaxial_effect_support_compliance_neglected}Obtained transfer functions from \(F\) to \(L^{\prime}\) when neglecting support compliance}
\end{figure}

\section{Effect of support compliance on \(L/F\)}

Some support compliance is now added and the model shown in Figure \ref{fig:uniaxial_support_compliance_test_system} is used.
The parameters of the support (i.e., \(m_{\mu}\), \(c_{\mu}\) and \(k_{\mu}\)) are chosen to match the vertical mode at \(70\,\text{Hz}\) seen on the micro-station (Figure \ref{fig:uniaxial_comp_frf_meas_model}).
The transfer functions from \(F\) to \(L\) (i.e., control of the relative motion of the nano-hexapod) and from \(L\) to \(d\) (i.e., control of the position between the nano-hexapod and the fixed granite) can then be computed.

When the relative displacement of the nano-hexapod \(L\) is controlled (dynamics shown in Figure \ref{fig:uniaxial_effect_support_compliance_dynamics}), having a stiff nano-hexapod (i.e., with a suspension mode at higher frequency than the mode of the support) makes the dynamics less affected by the limited support compliance (Figure \ref{fig:uniaxial_effect_support_compliance_dynamics_stiff}).
This is why it is very common to have stiff piezoelectric stages fixed at the very top of positioning stages.
In such a case, the control of the piezoelectric stage using its integrated metrology (typically capacitive sensors) is quite simple as the plant is not much affected by the dynamics of the support on which it is fixed.
If a soft nano-hexapod is used, the support dynamics appears in the dynamics between \(F\) and \(L\) (see Figure \ref{fig:uniaxial_effect_support_compliance_dynamics_soft}) which will impact the control robustness and performance.

\begin{figure}[htbp]
\begin{subfigure}{0.37\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_effect_support_compliance_dynamics_soft.png}
\end{center}
\subcaption{\label{fig:uniaxial_effect_support_compliance_dynamics_soft}$\omega_{n} \ll \omega_{\mu}$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_effect_support_compliance_dynamics_mid.png}
\end{center}
\subcaption{\label{fig:uniaxial_effect_support_compliance_dynamics_mid}$\omega_{n} = \omega_{\mu}$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_effect_support_compliance_dynamics_stiff.png}
\end{center}
\subcaption{\label{fig:uniaxial_effect_support_compliance_dynamics_stiff}$\omega_{n} \gg \omega_{\mu}$}
\end{subfigure}
\caption{\label{fig:uniaxial_effect_support_compliance_dynamics}Effect of the support compliance on the transfer functions from \(F\) to \(L\)}
\end{figure}

\section{Effect of support compliance on \(d/F\)}

When the motion to be controlled is the relative displacement \(d\) between the granite and the nano-hexapod's top platform (which is the case for the \acrshort{nass}), the effect of the support compliance on the plant dynamics is opposite to that previously observed.
Indeed, using a ``soft'' nano-hexapod (i.e., with a suspension mode at lower frequency than the mode of the support) makes the dynamics less affected by the support dynamics (Figure \ref{fig:uniaxial_effect_support_compliance_dynamics_d_soft}).
Conversely, if a ``stiff'' nano-hexapod is used, the support dynamics appears in the plant dynamics (Figure \ref{fig:uniaxial_effect_support_compliance_dynamics_d_stiff}).

\begin{figure}[htbp]
\begin{subfigure}{0.37\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_effect_support_compliance_dynamics_d_soft.png}
\end{center}
\subcaption{\label{fig:uniaxial_effect_support_compliance_dynamics_d_soft}$\omega_{n} \ll \omega_{\mu}$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_effect_support_compliance_dynamics_d_mid.png}
\end{center}
\subcaption{\label{fig:uniaxial_effect_support_compliance_dynamics_d_mid}$\omega_{n} = \omega_{\mu}$}
\end{subfigure}
\begin{subfigure}{0.31\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_effect_support_compliance_dynamics_d_stiff.png}
\end{center}
\subcaption{\label{fig:uniaxial_effect_support_compliance_dynamics_d_stiff}$\omega_{n} \gg \omega_{\mu}$}
\end{subfigure}
\caption{\label{fig:uniaxial_effect_support_compliance_dynamics_d}Effect of the support compliance on the transfer functions from \(F\) to \(d\)}
\end{figure}

\section*{Conclusion}
To study the impact of support compliance on plant dynamics, simple models shown in Figure \ref{fig:uniaxial_support_compliance_models} were used.
Depending on the quantity to be controlled (\(L\) or \(d\) in Figure \ref{fig:uniaxial_support_compliance_test_system}) and on the relative location of \(\omega_\nu\) (suspension mode of the nano-hexapod) with respect to \(\omega_\mu\) (modes of the support), the interaction between the support and the nano-hexapod dynamics can drastically change (observations made are summarized in Table \ref{tab:uniaxial_effect_compliance}).

For the \acrfull{nass}, having the suspension mode of the nano-hexapod at lower frequencies than the suspension modes of the micro-station would make the plant less dependent on the micro-station dynamics, and therefore easier to control.
Note that the observations made in this section are also affected by the ratio between the support mass \(m_{\mu}\) and the nano-hexapod mass \(m_n\) (the effect is more pronounced when the ratio \(m_n/m_{\mu}\) increases).

\begin{table}[htbp]
\centering
\begin{tabularx}{0.4\linewidth}{Xccc}
\toprule
 & \(\omega_{\nu} \ll \omega_{\mu}\) & \(\omega_{\nu} \approx \omega_{\mu}\) & \(\omega_{\nu} \gg \omega_{\mu}\)\\
\midrule
\(d/F\) & small & large & large\\
\(L/F\) & large & large & small\\
\bottomrule
\end{tabularx}
\caption{\label{tab:uniaxial_effect_compliance}Impact of the support dynamics on the plant dynamics}

\end{table}

\chapter{Effect of Payload Dynamics}
\label{sec:uniaxial_payload_dynamics}
Up to this section, the sample was modeled as a mass rigidly fixed to the nano-hexapod (as shown in Figure \ref{fig:uniaxial_paylaod_dynamics_rigid_schematic}).
However, such a sample may present internal dynamics, and its fixation to the nano-hexapod may have limited stiffness.
To study the effect of the sample dynamics, the models shown in Figure \ref{fig:uniaxial_paylaod_dynamics_schematic} are used.

\begin{figure}[htbp]
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_paylaod_dynamics_rigid_schematic.png}
\end{center}
\subcaption{\label{fig:uniaxial_paylaod_dynamics_rigid_schematic}Rigid payload}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,scale=1]{figs/uniaxial_paylaod_dynamics_schematic.png}
\end{center}
\subcaption{\label{fig:uniaxial_paylaod_dynamics_schematic}Payload with some flexibility}
\end{subfigure}
\caption{\label{fig:uniaxial_payload_dynamics_models}Models used to study the effect of payload dynamics}
\end{figure}

\section{Impact on plant dynamics}
\label{ssec:uniaxial_payload_dynamics_effect_dynamics}

To study the impact of the flexibility between the nano-hexapod and the payload, a first (reference) model with a rigid payload, as shown in Figure \ref{fig:uniaxial_paylaod_dynamics_rigid_schematic} is used.
Then ``flexible'' payload whose model is shown in Figure \ref{fig:uniaxial_paylaod_dynamics_schematic} are considered.
The resonances of the payload are set at \(\omega_s = 20\,\text{Hz}\) and at \(\omega_s = 200\,\text{Hz}\) while its mass is either \(m_s = 1\,\text{kg}\) or \(m_s = 50\,\text{kg}\).

The transfer functions from the nano-hexapod force \(f\) to the motion of the nano-hexapod top platform are computed for all the above configurations and are compared for a soft Nano-Hexapod (\(k_n = 0.01\,N/\mu m\)) in Figure \ref{fig:uniaxial_payload_dynamics_soft_nano_hexapod}.
It can be seen that the mode of the sample adds an anti-resonance followed by a resonance (zero/pole pattern).
The frequency of the anti-resonance corresponds to the ``free'' resonance of the sample \(\omega_s = \sqrt{k_s/m_s}\).
The flexibility of the sample also changes the high frequency gain (the mass line is shifted from \(\frac{1}{(m_n + m_s)s^2}\) to \(\frac{1}{m_ns^2}\)).

\begin{figure}[htbp]
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=\linewidth]{figs/uniaxial_payload_dynamics_soft_nano_hexapod_light.png}
\end{center}
\subcaption{\label{fig:uniaxial_payload_dynamics_soft_nano_hexapod_light}$k_n = 0.01\,N/\mu m$, $m_s = 1\,kg$}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=\linewidth]{figs/uniaxial_payload_dynamics_soft_nano_hexapod_heavy.png}
\end{center}
\subcaption{\label{fig:uniaxial_payload_dynamics_soft_nano_hexapod_heavy}$k_n = 0.01\,N/\mu m$, $m_s = 50\,kg$}
\end{subfigure}
\caption{\label{fig:uniaxial_payload_dynamics_soft_nano_hexapod}Effect of the payload dynamics on the soft Nano-Hexapod. Light sample (\subref{fig:uniaxial_payload_dynamics_soft_nano_hexapod_light}), and heavy sample (\subref{fig:uniaxial_payload_dynamics_soft_nano_hexapod_heavy})}
\end{figure}

The same transfer functions are now compared when using a stiff nano-hexapod (\(k_n = 100\,N/\mu m\)) in Figure \ref{fig:uniaxial_payload_dynamics_stiff_nano_hexapod}.
In this case, the sample's resonance \(\omega_s\) is smaller than the nano-hexapod resonance \(\omega_n\).
This changes the zero/pole pattern to a pole/zero pattern (the frequency of the zero still being equal to \(\omega_s\)).
Even though the added sample's flexibility still shifts the high frequency mass line as for the soft nano-hexapod, the dynamics below the nano-hexapod resonance is much less impacted, even when the sample mass is high and when the sample resonance is at low frequency (see yellow curve in Figure \ref{fig:uniaxial_payload_dynamics_stiff_nano_hexapod_heavy}).

\begin{figure}[htbp]
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=\linewidth]{figs/uniaxial_payload_dynamics_stiff_nano_hexapod_light.png}
\end{center}
\subcaption{\label{fig:uniaxial_payload_dynamics_stiff_nano_hexapod_light}$k_n = 100\,N/\mu m$, $m_s = 1\,kg$}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=\linewidth]{figs/uniaxial_payload_dynamics_stiff_nano_hexapod_heavy.png}
\end{center}
\subcaption{\label{fig:uniaxial_payload_dynamics_stiff_nano_hexapod_heavy}$k_n = 100\,N/\mu m$, $m_s = 50\,kg$}
\end{subfigure}
\caption{\label{fig:uniaxial_payload_dynamics_stiff_nano_hexapod}Effect of the payload dynamics on the stiff Nano-Hexapod. Light sample (\subref{fig:uniaxial_payload_dynamics_stiff_nano_hexapod_light}), and heavy sample (\subref{fig:uniaxial_payload_dynamics_stiff_nano_hexapod_heavy})}
\end{figure}

\section{Impact on close loop performances}
\label{ssec:uniaxial_payload_dynamics_effect_stability}

Having a flexibility between the measured position (i.e., the top platform of the nano-hexapod) and the point-of-interest to be positioned relative to the x-ray may also impact the closed-loop performance (i.e., the remaining sample's vibration).

To estimate whether the sample flexibility is critical for the closed-loop position stability of the sample, the model shown in Figure \ref{fig:uniaxial_sample_flexibility_control} is used.
This is the same model that was used in Section \ref{sec:uniaxial_position_control} but with an added flexibility between the nano-hexapod and the sample (considered sample modes are at \(\omega_s = 20\,\text{Hz}\) and \(\omega_n = 200\,\text{Hz}\)).
In this case, the measured (i.e., controlled) distance \(d\) is no longer equal to the real performance index (the distance \(y\)).

\begin{figure}[htbp]
\centering
\includegraphics[scale=1]{figs/uniaxial_sample_flexibility_control.png}
\caption{\label{fig:uniaxial_sample_flexibility_control}Uniaxial model considering some flexibility between the nano-hexapod top platform and the sample. In this case, the measured and controlled distance \(d\) is different from the distance \(y\) which is the real performance index}
\end{figure}

The system dynamics is computed and IFF is applied using the same gains as those used in Section \ref{sec:uniaxial_active_damping}.
Due to the collocation between the nano-hexapod and the force sensor used for IFF, the damped plants are still stable and similar damping values are obtained than when considering a rigid sample.
The High Authority Controllers used in Section \ref{sec:uniaxial_position_control} are then implemented on the damped plants.
The obtained closed-loop systems are stable, indicating good robustness.

Finally, closed-loop noise budgeting is computed for the obtained closed-loop system, and the cumulative amplitude spectrum of \(d\) and \(y\) are shown in Figure \ref{fig:uniaxial_sample_flexibility_noise_budget_y}.
The cumulative amplitude spectrum of the measured distance \(d\) (Figure \ref{fig:uniaxial_sample_flexibility_noise_budget_d}) shows that the added flexibility at the sample location has very little effect on the control performance.
However, the cumulative amplitude spectrum of the distance \(y\) (Figure \ref{fig:uniaxial_sample_flexibility_noise_budget_y}) shows that the stability of \(y\) is degraded when the sample flexibility is considered and is degraded as \(\omega_s\) is lowered.

What happens is that above \(\omega_s\), even though the motion \(d\) can be controlled perfectly, the sample's mass is ``isolated'' from the motion of the nano-hexapod and the control on \(y\) is not effective.

\begin{figure}[htbp]
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_sample_flexibility_noise_budget_d.png}
\end{center}
\subcaption{\label{fig:uniaxial_sample_flexibility_noise_budget_d}Cumulative Amplitude Spectrum of $d$}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\begin{center}
\includegraphics[scale=1,width=0.95\linewidth]{figs/uniaxial_sample_flexibility_noise_budget_y.png}
\end{center}
\subcaption{\label{fig:uniaxial_sample_flexibility_noise_budget_y}Cumulative Amplitude Spectrum of $y$}
\end{subfigure}
\caption{\label{fig:uniaxial_sample_flexibility_noise_budget}Cumulative Amplitude Spectrum of the distances \(d\) and \(y\). The effect of the sample's flexibility does not affect much \(d\) but is detrimental to the stability of \(y\). A sample mass \(m_s = 1\,\text{kg}\) and a nano-hexapod stiffness of \(100\,N/\mu m\) are used for the simulations.}
\end{figure}

\section*{Conclusion}
Payload dynamics is usually a major concern when designing a positioning system.
In this section, the impact of the sample dynamics on the plant was found to vary with the sample mass and the relative resonance frequency of the sample \(\omega_s\) and of the nano-hexapod \(\omega_n\).
The larger the sample mass, the larger the effect (i.e., change of high frequency gain, appearance of additional resonances and anti-resonances).
A zero/pole pattern is observed if \(\omega_s > \omega_n\) and a pole/zero pattern if \(\omega_s > \omega_n\).
Such additional dynamics can induce stability issues depending on their position relative to the desired feedback bandwidth, as explained in \cite[Section 4.2]{rankers98_machin}.
The general conclusion is that the stiffer the nano-hexapod, the less it is impacted by the payload's dynamics, which would make the feedback controller more robust to a change of payload.
This is why high-bandwidth soft positioning stages are usually restricted to constant and calibrated payloads (CD-player, lithography machines, isolation system for gravitational wave detectors, \ldots{}), whereas stiff positioning systems are usually used when the control must be robust to a change of payload's mass (stiff piezo nano-positioning stages for instance).

Having some flexibility between the measurement point and the point of interest (i.e., the sample point to be position on the x-ray) also degrades the position stability as shown in Section \ref{ssec:uniaxial_payload_dynamics_effect_stability}.
Therefore, it is important to take special care when designing sampling environments, especially if a soft nano-hexapod is used.

\chapter*{Conclusion}
\label{sec:uniaxial_conclusion}

In this study, a uniaxial model of the nano-active-stabilization-system was tuned from both dynamical measurements (Section \ref{sec:micro_station_model}) and from disturbances measurements (Section \ref{sec:uniaxial_disturbances}).

Three active damping techniques can be used to critically damp the nano-hexapod resonances (Section \ref{sec:uniaxial_active_damping}).
However, this model does not allow the determination of which one is most suited to this application (a comparison of the three active damping techniques is done in Table \ref{tab:comp_active_damping}).

Position feedback controllers have been developed for three considered nano-hexapod stiffnesses (Section \ref{sec:uniaxial_position_control}).
These controllers were shown to be robust to the change of sample's masses, and to provide good rejection of disturbances.
Having a soft nano-hexapod makes the plant dynamics easier to control (because its dynamics is decoupled from the micro-station dynamics, see Section \ref{sec:uniaxial_support_compliance}) and requires less position feedback bandwidth to fulfill the requirements.
The moderately stiff nano-hexapod (\(k_n = 1\,N/\mu m\)) is requiring a higher feedback bandwidth, but still gives acceptable results.
However, the stiff nano-hexapod is the most complex to control and gives the worst positioning performance.

\printbibliography[heading=bibintoc,title={Bibliography}]

\printglossaries
\end{document}
