Active Damping applied on the Simscape Model
Table of Contents
First, in section 1, we will looked at the undamped system.
Then, we will compare three active damping techniques:
- In section 2: the integral force feedback is used
- In section 3: the direct velocity feedback is used
- In section 4: inertial control is used
For each of the active damping technique, we will:
- Look at the damped plant
- Simulate tomography experiments
- Compare the sensitivity from disturbances
The disturbances are:
- Ground motion
- Motion errors of all the stages
1 Undamped System
All the files (data and Matlab scripts) are accessible here.
We first look at the undamped system. The performance of this undamped system will be compared with the damped system using various techniques.
1.1 Identification of the dynamics for Active Damping
1.1.1 Initialize the Simulation
We initialize all the stages with the default parameters.
initializeGround(); initializeGranite(); initializeTy(); initializeRy(); initializeRz(); initializeMicroHexapod(); initializeAxisc(); initializeMirror();
The nano-hexapod is a piezoelectric hexapod and the sample has a mass of 50kg.
initializeNanoHexapod('actuator', 'piezo'); initializeSample('mass', 50);
We set the references to zero.
initializeReferences('Rz_type', 'rotating', 'Ry_period', 1);
And all the controllers are set to 0.
K = tf(zeros(6)); save('./mat/controllers.mat', 'K', '-append'); K_ine = tf(zeros(6)); save('./mat/controllers.mat', 'K_ine', '-append'); K_iff = tf(zeros(6)); save('./mat/controllers.mat', 'K_iff', '-append'); K_dvf = tf(zeros(6)); save('./mat/controllers.mat', 'K_dvf', '-append');
1.1.2 Identification
First, we identify the dynamics of the system using the linearize
function.
%% Options for Linearized options = linearizeOptions; options.SampleTime = 0; %% Name of the Simulink File mdl = 'sim_nass_active_damping'; %% Input/Output definition clear io; io_i = 1; io(io_i) = linio([mdl, '/Fnl'], 1, 'openinput'); io_i = io_i + 1; io(io_i) = linio([mdl, '/Micro-Station'], 3, 'openoutput', [], 'Dnlm'); io_i = io_i + 1; io(io_i) = linio([mdl, '/Micro-Station'], 3, 'openoutput', [], 'Fnlm'); io_i = io_i + 1; io(io_i) = linio([mdl, '/Micro-Station'], 3, 'openoutput', [], 'Vlm'); io_i = io_i + 1; %% Run the linearization G = linearize(mdl, io, 1, options); G.InputName = {'Fnl1', 'Fnl2', 'Fnl3', 'Fnl4', 'Fnl5', 'Fnl6'}; G.OutputName = {'Dnlm1', 'Dnlm2', 'Dnlm3', 'Dnlm4', 'Dnlm5', 'Dnlm6', ... 'Fnlm1', 'Fnlm2', 'Fnlm3', 'Fnlm4', 'Fnlm5', 'Fnlm6', ... 'Vnlm1', 'Vnlm2', 'Vnlm3', 'Vnlm4', 'Vnlm5', 'Vnlm6'};
We then create transfer functions corresponding to the active damping plants.
G_iff = minreal(G({'Fnlm1', 'Fnlm2', 'Fnlm3', 'Fnlm4', 'Fnlm5', 'Fnlm6'}, {'Fnl1', 'Fnl2', 'Fnl3', 'Fnl4', 'Fnl5', 'Fnl6'})); G_dvf = minreal(G({'Dnlm1', 'Dnlm2', 'Dnlm3', 'Dnlm4', 'Dnlm5', 'Dnlm6'}, {'Fnl1', 'Fnl2', 'Fnl3', 'Fnl4', 'Fnl5', 'Fnl6'})); G_ine = minreal(G({'Vnlm1', 'Vnlm2', 'Vnlm3', 'Vnlm4', 'Vnlm5', 'Vnlm6'}, {'Fnl1', 'Fnl2', 'Fnl3', 'Fnl4', 'Fnl5', 'Fnl6'}));
And we save them for further analysis.
save('./active_damping/mat/undamped_plants.mat', 'G_iff', 'G_dvf', 'G_ine');
1.2 Tomography Experiment
1.2.1 Simulation
We initialize elements for the tomography experiment.
prepareTomographyExperiment();
We change the simulation stop time.
load('mat/conf_simscape.mat'); set_param(conf_simscape, 'StopTime', '3');
And we simulate the system.
sim('sim_nass_active_damping');
Finally, we save the simulation results for further analysis
save('./active_damping/mat/tomo_exp.mat', 'En', 'Eg', '-append');
1.2.2 Results
2 Integral Force Feedback
All the files (data and Matlab scripts) are accessible here.
Integral Force Feedback is applied on the simscape model.
2.1 Control Design
2.1.1 Plant
Let’s load the previously indentified undamped plant:
load('./active_damping/mat/undamped_plants.mat', 'G_iff');
Let’s look at the transfer function from actuator forces in the nano-hexapod to the force sensor in the nano-hexapod legs for all 6 pairs of actuator/sensor (figure 6).
2.1.2 Control Design
2.1.3 Diagonal Controller
We create the diagonal controller and we add a minus sign as we have a positive feedback architecture.
K_iff = -K_iff*eye(6);
We save the controller for further analysis.
save('./active_damping/mat/K_iff.mat', 'K_iff');
2.1.4 IFF with High Pass Filter
w_hpf = 2*pi*10; % Cut-off frequency for the high pass filter [rad/s] K_iff = 2*pi*200/s * (s/w_hpf)/(s/w_hpf + 1);
The corresponding loop gains are shown in figure 8.
We create the diagonal controller and we add a minus sign as we have a positive feedback architecture.
K_iff = -K_iff*eye(6);
We save the controller for further analysis.
save('./active_damping/mat/K_iff_hpf.mat', 'K_iff');
2.2 Tomography Experiment
2.2.1 Simulation with IFF Controller
We initialize elements for the tomography experiment.
prepareTomographyExperiment();
We set the IFF controller.
load('./active_damping/mat/K_iff.mat', 'K_iff'); save('./mat/controllers.mat', 'K_iff', '-append');
We change the simulation stop time.
load('mat/conf_simscape.mat'); set_param(conf_simscape, 'StopTime', '3');
And we simulate the system.
sim('sim_nass_active_damping');
Finally, we save the simulation results for further analysis
En_iff = En; Eg_iff = Eg; save('./active_damping/mat/tomo_exp.mat', 'En_iff', 'Eg_iff', '-append');
2.2.2 Simulation with IFF Controller with added High Pass Filter
We initialize elements for the tomography experiment.
prepareTomographyExperiment();
We set the IFF controller with the High Pass Filter.
load('./active_damping/mat/K_iff_hpf.mat', 'K_iff'); save('./mat/controllers.mat', 'K_iff', '-append');
We change the simulation stop time.
load('mat/conf_simscape.mat'); set_param(conf_simscape, 'StopTime', '3');
And we simulate the system.
sim('sim_nass_active_damping');
Finally, we save the simulation results for further analysis
En_iff_hpf = En; Eg_iff_hpf = Eg; save('./active_damping/mat/tomo_exp.mat', 'En_iff_hpf', 'Eg_iff_hpf', '-append');
2.2.3 Compare with Undamped system
We load the results of tomography experiments.
load('./active_damping/mat/tomo_exp.mat', 'En', 'En_iff', 'En_iff_hpf'); t = linspace(0, 3, length(En(:,1)));
2.3 Conclusion
Integral Force Feedback:
- Robust (guaranteed stability)
- Acceptable Damping
- Increase the sensitivity to disturbances at low frequencies
3 Direct Velocity Feedback
All the files (data and Matlab scripts) are accessible here.
In the Direct Velocity Feedback (DVF), a derivative feedback is applied between the measured actuator displacement to the actuator force input. The actuator displacement can be measured with a capacitive sensor for instance.
3.1 Control Design
3.1.1 Plant
Let’s load the undamped plant:
load('./active_damping/mat/undamped_plants.mat', 'G_dvf');
Let’s look at the transfer function from actuator forces in the nano-hexapod to the measured displacement of the actuator for all 6 pairs of actuator/sensor (figure 12).
3.1.2 Control Design
3.1.3 Diagonal Controller
We create the diagonal controller and we add a minus sign as we have a positive feedback architecture.
K_dvf = -K_dvf*eye(6);
We save the controller for further analysis.
save('./active_damping/mat/K_dvf.mat', 'K_dvf');
3.2 Tomography Experiment
3.2.1 Initialize the Simulation
We initialize elements for the tomography experiment.
prepareTomographyExperiment();
We set the DVF controller.
load('./active_damping/mat/K_dvf.mat', 'K_dvf'); save('./mat/controllers.mat', 'K_dvf', '-append');
3.2.2 Simulation
We change the simulation stop time.
load('mat/conf_simscape.mat'); set_param(conf_simscape, 'StopTime', '3');
And we simulate the system.
sim('sim_nass_active_damping');
Finally, we save the simulation results for further analysis
En_dvf = En; Eg_dvf = Eg; save('./active_damping/mat/tomo_exp.mat', 'En_dvf', 'Eg_dvf', '-append');
3.2.3 Compare with Undamped system
We load the results of tomography experiments.
load('./active_damping/mat/tomo_exp.mat', 'En', 'En_dvf'); t = linspace(0, 3, length(En(:,1)));
3.3 Conclusion
Direct Velocity Feedback:
4 Inertial Control
All the files (data and Matlab scripts) are accessible here.
In Inertial Control, a feedback is applied between the measured absolute motion (velocity or acceleration) of the platform to the actuator force input.
4.1 Control Design
4.1.1 Plant
Let’s load the undamped plant:
load('./active_damping/mat/undamped_plants.mat', 'G_ine');
Let’s look at the transfer function from actuator forces in the nano-hexapod to the measured velocity of the nano-hexapod platform in the direction of the corresponding actuator for all 6 pairs of actuator/sensor (figure 17).
4.1.2 Control Design
4.1.3 Diagonal Controller
We create the diagonal controller and we add a minus sign as we have a positive feedback architecture.
K_ine = -K_ine*eye(6);
We save the controller for further analysis.
save('./active_damping/mat/K_ine.mat', 'K_ine');
4.2 Tomography Experiment
4.2.1 Initialize the Simulation
We initialize elements for the tomography experiment.
prepareTomographyExperiment();
We set the Inertial controller.
load('./active_damping/mat/K_ine.mat', 'K_ine'); save('./mat/controllers.mat', 'K_ine', '-append');
4.2.2 Simulation
We change the simulation stop time.
load('mat/conf_simscape.mat'); set_param(conf_simscape, 'StopTime', '3');
And we simulate the system.
sim('sim_nass_active_damping');
Finally, we save the simulation results for further analysis
En_ine = En; Eg_ine = Eg; save('./active_damping/mat/tomo_exp.mat', 'En_ine', 'Eg_ine', '-append');
4.2.3 Compare with Undamped system
We load the results of tomography experiments.
load('./active_damping/mat/tomo_exp.mat', 'En', 'En_ine'); t = linspace(0, 3, length(En_ine(:,1)));
4.3 Conclusion
Inertial Control:
5 Variability of the system dynamics for Active Damping
All the files (data and Matlab scripts) are accessible here.
5.1 Variation of the Sample Mass
For all the identifications, the disturbances are disabled and no controller are used.
We identify the dynamics for the following sample mass.
masses = [1, 10, 50]; % [kg]
5.2 Variation of the Spindle Angle
We identify the dynamics for the following Spindle angles.
Rz_amplitudes = [0, pi/4, pi/2, pi]; % [rad]
5.3 Variation of the Spindle Rotation Speed
We identify the dynamics for the following Spindle rotation periods.
Rz_periods = [60, 10, 1]; % [s]
5.4 Variation of the Tilt Angle
We identify the dynamics for the following Tilt stage angles.
Ry_amplitudes = [0, 3*pi/180]; % [rad]
5.5 Conclusion
6 Comparison
6.1 Load the plants
load('./active_damping/mat/plants.mat', 'G', 'G_iff', 'G_ine', 'G_dvf');
6.2 Sensitivity to Disturbance
Figure 38: Compliance in the Z direction: Sensitivity of direct forces applied on the sample in the Z direction on the Z motion error (png, pdf)
Figure 39: Sensitivity to forces applied in the Z direction by the Spindle on the Z motion error (png, pdf)
6.3 Damped Plant
6.4 Tomography Experiment
6.4.1 Load the Simulation Data
load('./active_damping/mat/tomo_exp.mat', 'En', 'En_iff_hpf', 'En_dvf', 'En_ine'); En_iff = En_iff_hpf; t = linspace(0, 3, length(En(:,1)));
6.4.2 Frequency Domain Analysis
Window used for pwelch
function.
n_av = 8; han_win = hanning(ceil(length(En(:, 1))/n_av));
Figure 45: PSD of the translation errors in the X direction for applied Active Damping techniques (png, pdf)
Figure 46: PSD of the rotation errors in the X direction for applied Active Damping techniques (png, pdf)
7 Useful Functions
7.1 prepareTomographyExperiment
This Matlab function is accessible here.
Function Description
function [] = prepareTomographyExperiment(args)
Optional Parameters
arguments args.nass_actuator char {mustBeMember(args.nass_actuator,{'piezo', 'lorentz'})} = 'piezo' args.sample_mass (1,1) double {mustBeNumeric, mustBePositive} = 50 args.Ry_period (1,1) double {mustBeNumeric, mustBePositive} = 1 end
Initialize the Simulation
We initialize all the stages with the default parameters.
initializeGround(); initializeGranite(); initializeTy(); initializeRy(); initializeRz(); initializeMicroHexapod(); initializeAxisc(); initializeMirror();
The nano-hexapod is a piezoelectric hexapod and the sample has a mass of 50kg.
initializeNanoHexapod('actuator', args.nass_actuator); initializeSample('mass', args.sample_mass);
We set the references to zero.
initializeReferences('Rz_type', 'rotating', 'Rz_period', args.Ry_period);
And all the controllers are set to 0.
K = tf(zeros(6)); save('./mat/controllers.mat', 'K', '-append'); K_ine = tf(zeros(6)); save('./mat/controllers.mat', 'K_ine', '-append'); K_iff = tf(zeros(6)); save('./mat/controllers.mat', 'K_iff', '-append'); K_dvf = tf(zeros(6)); save('./mat/controllers.mat', 'K_dvf', '-append');