wecSimInputFile.m

%% Simulation Data
simu = simulationClass();
simu.simMechanicsFile = 'OSWEC_RO.slx';     % Specify Simulink Model File
%simu.mode = 'rapid-accelerator';           % Specify Simulation Mode ('normal','accelerator','rapid-accelerator')
simu.explorer='off';                        % Turn SimMechanics Explorer (on/off)
simu.solver = 'ode4';                       %simu.solver = 'ode4' for fixed step & simu.solver = 'ode45' for variable step 
simu.startTime = 0;                         % Simulation Start Time [s]
simu.rampTime = 50;
simu.endTime = 300;
simu.dt = 0.01;                             %Simulation time-step [s] for a convolution function in the radiation force calculation 
simu.cicEndTime = 30;

%% Wave Information
%% Irregular Waves using PM Spectrum with Convolution Integral Calculation
waves = waveClass('irregular');         % Initialize Wave Class and Specify Type
waves.height = 2.64;                         % Significant Wave Height [m]
waves.period = 9.86;                         % Peak Period [s]
waves.spectrumType = 'PM';              % Specify Wave Spectrum Type
waves.bem.option = 'EqualEnergy';
waves.bem.count = 250;
waves.phaseSeed = 1;

%% Body Data
% Flap
body(1) = bodyClass('../_Common_Input_Files/OSWEC/hydroData/oswec.h5');   % Initialize bodyClass for Flap
body(1).geometryFile = '../_Common_Input_Files/OSWEC/geometry/flap.stl';  % Geometry File
body(1).mass = 127000;                         % User-Defined mass [kg]
body(1).inertia = [1.85e6 1.85e6 1.85e6]; % Moment of Inertia [kg-m^2]
body(1).morisonElement.option = 1;
body(1).morisonElement.cd = ones (5,3);
body(1).morisonElement.ca = zeros(5,3);
body(1).morisonElement.area = zeros(5,3);
body(1).morisonElement.area(:,1) = 18*1.8;
body(1).morisonElement.area(:,3) = 18*1.8;
body(1).morisonElement.VME  = zeros(5,1);
body(1).morisonElement.rgME = [0 0 -3; 0 0 -1.2; 0 0 0.6; 0 0 2.4; 0 0 4.2];

% Base
body(2) = bodyClass('../_Common_Input_Files/OSWEC/hydroData/oswec.h5');   % Initialize bodyClass for Base
body(2).geometryFile = '../_Common_Input_Files/OSWEC/geometry/base.stl';  % Geometry File
body(2).mass = 999;                            % Placeholder mass for fixed body
body(2).inertia = [999 999 999];               % Placeholder inertia for fixed body

%% PTO and Constraint Parameters
constraint(1)= constraintClass('Constraint1'); % Initialize ConstraintClass 
constraint(1).location = [0 0 -10];

constraint(2)= constraintClass('Constraint2'); % Initialize ConstraintClass 
constraint(2).location = [0 0 -8.9];

constraint(3)= constraintClass('Constraint3'); % Initialize ConstraintClass 
constraint(3).location = [4.7021271782+0.9 0 -8.7];

constraint(4)= constraintClass('Constraint4'); % Initialize ConstraintClass 
constraint(4).location = [0+0.9 0 -7];

inertiaBlock = 1e-3; % kg - small mass on follower side of PTO used to resolve degenerate mass errors

pto(1) = ptoClass('PTO1');                     % Initialize ptoClass for PTO1
pto(1).stiffness = 0;                          % PTO Stiffness Coeff [Nm/rad]
pto(1).damping = 0;                            % PTO Damping Coeff [Nsm/rad]
pto(1).location = [2.35106397378+0.9 0 -7.849998936];   % PTO Global Location [m]
pto(1).orientation.z = [-4.7021271782/5 0 1.7/5];       % PTO orientation