load params from file

turbie/MATLAB_model/main.m

 clc; clear; close all
 
+%%
+
 % Figure settings    
 set(0,'DefaultAxesFontSize',14)
 set(0,'DefaultTextFontSize',14)
@@ -10,25 +12,26 @@ global t f M C K rho Dr CT V
 
 %% Turbine parameters
 
-mb = 41e3;          % Blade mass [kg]
-mn = 446e3;         % Nacelle mass [kg]
-mh = 105e3;         % Hub mass [kg]
-mt = 628e3;         % Tower mass [kg]
+parameters = readvars('../turbie_parameters.txt', 'CommentStyle','%',  'NumHeaderLines',1, 'ExpectedNumVariables',14);
+mb = parameters(1);         % Blade mass [kg]
+mn = parameters(2);         % Nacelle mass [kg]
+mh = parameters(3);         % Hub mass [kg]
+mt = parameters(4);         % Tower mass [kg]
 
-c1 = 4.208e3;       % Blade flap structural damping [N/(m/s)]
-c2 = 1.273e4;       % Tower fore-aft structural damping [N/(m/s)]
+c1 = parameters(5);       % Blade flap structural damping [N/(m/s)]
+c2 = parameters(6);       % Tower fore-aft structural damping [N/(m/s)]
 
-k1 = 1.711e6;       % Blade flap stiffness [N/m]
-k2 = 3.278e6;       % Tower fore-aft stiffness [N/m]
+k1 = parameters(7);       % Blade flap stiffness [N/m]
+k2 = parameters(8);       % Tower fore-aft stiffness [N/m]
 
-fb = 0.63;          % First flap blade modal frequency [Hz]
-ft = 0.25;          % First fore-aft tower modal frequency [Hz]
+fb = parameters(9);          % First flap blade modal frequency [Hz]
+ft = parameters(10);         % First fore-aft tower modal frequency [Hz]
 
-drb = 0.031/2/pi;   % First flap blade modal damping ratio [-]
-drt = 0.019/2/pi;   % First fore-aft tower modal damping ratio [-]
+drb = parameters(11);   % First flap blade modal damping ratio [-]
+drt = parameters(12);   % First fore-aft tower modal damping ratio [-]
 
-Dr = 178;           % Rotor diameter (m)
-rho = 1.22;         % Air density (kg/m3)
+Dr = parameters(13);          % Rotor diameter (m)
+rho = parameters(14);         % Air density (kg/m3)
 
 %% System matrices
 
@@ -79,15 +82,11 @@ table(fntar,fnu,fnd,fnerr,drtar,dr,drerr)
 
 %% Thrust coefficient
 
-cd ..\
-
 % Load table with CT as function of mean wind speed
-temp = importdata('CT.txt');
+temp = importdata('../CT.txt');
 Vtab = temp.data(:,1);
 CTtab = temp.data(:,2);
 
-cd MATLAB_model
-
 figure()
 plot(Vtab,CTtab,'ko'), hold on, grid minor, axis tight
 xlabel(temp.colheaders{1})

turbie/turbie_parameters.txt

-%% Turbine parameters
-
-mb = 41e3;          % Blade mass [kg]
-mn = 446e3;         % Nacelle mass [kg]
-mh = 105e3;         % Hub mass [kg]
-mt = 628e3;         % Tower mass [kg]
-
-c1 = 4.208e3;       % Blade flap structural damping [N/(m/s)]
-c2 = 1.273e4;       % Tower fore-aft structural damping [N/(m/s)]
-
-k1 = 1.711e6;       % Blade flap stiffness [N/m]
-k2 = 3.278e6;       % Tower fore-aft stiffness [N/m]
-
-fb = 0.63;          % First flap blade modal frequency [Hz]
-ft = 0.25;          % First fore-aft tower modal frequency [Hz]
-
-drb = 0.031/2/pi;   % First flap blade modal damping ratio [-]
-drt = 0.019/2/pi;   % First fore-aft tower modal damping ratio [-]
-
-Dr = 178;           % Rotor diameter (m)
-rho = 1.22;         % Air density (kg/m3)
\ No newline at end of file
+% Turbie parameters
+41e3     % mb [kg]
+446e3    % mn [kg]
+105e3    % mh [kg]
+628e3    % mt [kg]
+4.208e3  % c1 [N/(m/s)]
+1.273e4  % c2 [N/(m/s)]
+1.711e6  % k1 [N/(m/s)]
+3.278e6  % k2 [N/(m/s)]
+0.63     % fb [Hz]
+0.25     % ft [Hz]
+0.004933803235848756  % drb [-]
+0.0030239439187460114  % drt [-]
+178   % Dr [m]
+1.22  % rho [kg/m3]