Commit 1a73b48a authored by Alan Wai Hou Lio's avatar Alan Wai Hou Lio

tidy up the code

parent 07e82d2d
;DTU_10MW_RWT, version 9, 25-09-2017, mhha
;
begin simulation;
time_stop 150;
time_stop 600;
solvertype 1 ; (newmark)
on_no_convergence continue ;
convergence_limits 1E3 1.0 1E-7 ;
......@@ -289,7 +289,7 @@ end new_htc_structure;
;----------------------------------------------------------------------------------------------------------------------------------------------------------------
begin wind ;
density 1.225 ;
wsp 14.0 ;
wsp 4.0 ;
tint 0.0 ;
horizontal_input 1 ;
windfield_rotations 0.0 0.0 0.0 ; yaw, tilt, rotation
......@@ -517,6 +517,7 @@ begin type2_dll;
mbdy momentvec blade1 2 2 blade1 only 1; # blade 1 root [kNm] ;
mbdy momentvec blade2 2 2 blade2 only 1; # blade 2 root [kNm];
mbdy momentvec blade3 2 2 blade3 only 1; # blade 3 root [kNm];
dll inpvec 1 45;
end output;
;
begin actions;
......@@ -637,7 +638,7 @@ end dll;
;
begin output;
filename ./res/DTU_10MW_RWT ;
data_format hawc_ascii;
data_format hawc_binary;
buffer 1 ;
;
general time;
......
......@@ -19,6 +19,7 @@ while time <= Tend
Mb1 = ArrayFromHawc(5);% 8: blade1 root 1 Mx [kNm]
Mb2 = ArrayFromHawc(6);% 11: blade1 root 2 Mx [kNm]
Mb3 = ArrayFromHawc(7);% 14: blade1 root 3 Mx [kNm]
dr = ArrayFromHawc(8);
% check simulation status
% clc;
% fprintf('----------------------------------'); fprintf('\n');
......@@ -31,7 +32,7 @@ while time <= Tend
theta1 = 0; % additional pithc angle
theta2 = 0;
theta3 = 0;
dr_update = 0.6; % de-rated at 0.6*Prated.
dr_update = dr; % Derate percentage (e.g. 0.7 means 70% of nominal power)
%%%%%%%%%%%%%%%%%%% Signals to HAWC2 output %%%%%%%%%%%%%%%%%
ArrayToHawc = zeros(40,1);
......
......@@ -406,7 +406,8 @@ subroutine init_regulation_advanced(array1, array2) bind(c,name='init_regulation
! Initialization
TimerExcl = -0.02_mk
! Derating parameters
Deratevar%strat = array1(79)
Deratevar%strat = array1(79)
write(6,*) "Warning: This controller only works if the power set-point is provided as a controller input ! "
if (Deratevar%strat > 0) then
Deratevar%dr = array1(80)/100.0_mk
endif
......@@ -683,7 +684,6 @@ subroutine update_regulation(array1, array2) bind(c,name='update_regulation')
! Real-time derating
if (time .gt. 10.0_mk) then
Deratevar%dr = array1(13)
print*, Deratevar%dr,CtrlStatus
endif
!***********************************************************************************************
! Safety system
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment