clear all;
close all;
clc;
%-----------------------------------------------------------------------------
-
% Disturbance variables
%-----------------------------------------------------------------------------
-
% Compressed supply pressure
p_sup = 9.5e5; % [Pa]
% Ambient pressure
p_amb = 1e5; % [Pa]
% Ambient temperature
T_amb = 293; % [K]
%-----------------------------------------------------------------------------
-
% Parameter setup
%-----------------------------------------------------------------------------
-
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Physical constants
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specific gas constant of air
R_air = 287.14; % [J/kgK]
% Adiabatic exponent of air
kappa_air = 1.4; % [-]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Magnet valve parameters - 1.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Intake cross section diameter
d_in_1 = 0.0028; % [m]
% Intake cross section
A_in_1 = (d_in_1^2*pi)/4; % [m^2]
% Intake contraction coefficient
alpha_in_1 = 0.55; % [-]
% Exhaust cross section diameter
d_exh_1 = 0.0030; % [m]
% Exhaust cross section
A_exh_1 = (d_exh_1^2*pi)/4; % [m^2]
% Exhaust contraction coefficient
alpha_exh_1 = 0.55; % [-]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Magnet valve parameters - 2.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Intake cross section diameter
d_in_2 = 0.0028; % [m]
% Intake cross section
A_in_2 = (d_in_2^2*pi)/4; % [m^2]
% Intake contraction coefficient
alpha_in_2 = 0.55; % [-]
% Exhaust cross section diameter
d_exh_2 = 0.0030; % [m]
% Exhaust cross section
A_exh_2 = (d_exh_2^2*pi)/4; % [m^2]
% Exhaust contraction coefficient
alpha_exh_2 = 0.55; % [-]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Chamber parameters
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Piston rod diameter
d_rod = 0.016; % [m]
% Piston diameter
d_pst = 0.063; % [m]
% Cross section area of piston
A_rod = (d_rod^2*pi)/4; % [m^2]
A_pst = (d_pst^2*pi)/4; % [m^2]
% Dead volumes of chambers
V_ch_d_1 = 53e-6; % [m^3]
V_ch_d_2 = 49e-6
; % [m^3]
% Stroke limiting of piston
x_pst_max = 0.0307; % [m]
x_pst_min = 0.0025; % [m]
% Mass of piston
m_pst = 0.65; % [kg]
% Damping coefficient
k_pst = 500; % [Nm/s]
%-----------------------------------------------------------------------------
-
% Initial conditions
%-----------------------------------------------------------------------------
-
% Piston velocity
v_pst_0 = 0; % [m/s]
% Piston displacement
x_pst_0 = x_pst_min; % [m]
% Chambers pressures
p_ch1_0 = p_amb; % [Pa]
p_ch2_0 = p_amb; % [Pa]
load('KOGJA512_SzimuHF_v31');
y_ch1 = max(meas_res.signals(3).values)-min(meas_res.signals(3).values);
y_ch2 = max(meas_res.signals(4).values)-min(meas_res.signals(4).values);
y_pst = max(meas_res.signals(5).values)-min(meas_res.signals(5).values);