Commit d11d6820 authored by Jenni Rinker's avatar Jenni Rinker
Browse files

notrack step/hs2, auto from ci

parent 341fd73e
; DTU 10MW RWT. https://gitlab.windenergy.dtu.dk/hawc-reference-models/dtu-10-mw
;
begin new_htc_structure; material properties, orientations and constraints of bodies
; beam_output_file_name ./log/DTU_10MW_RWT_beam.dat; Optional - beam properties of bodies to file
; body_output_file_name ./log/DTU_10MW_RWT_body.dat; Optional - initial position/orientation of bodies to file
; body_eigenanalysis_file_name ./eig/DTU_10MW_RWT_body_eigen.dat; Optional - eigenvalues of bodies to file
; structure_eigenanalysis_file_name ./eig/DTU_10MW_RWT_strc_eigen.dat ; Optional - eigenvalues of structure to file
;-------------------------------------------------------------------------------------------------------------------------------
begin main_body; tower
name tower;
type timoschenko;
nbodies 1;
node_distribution c2_def;
damping_posdef 0 0 0 0.00412 0.00412 0.00045; tuned by Larh
begin timoschenko_input;
filename ./data/DTU_10MW_RWT_Tower_st.dat;
set 1 1; 1 1=flexibe, 1 2=stiff
end timoschenko_input;
begin c2_def; definition of centerline in main_body coordinates
nsec 11;
sec 1 0 0 0 0; x, y, z, twist
sec 2 0 0 -11.5 0;
sec 3 0 0 -23 0;
sec 4 0 0 -34.5 0;
sec 5 0 0 -46 0;
sec 6 0 0 -57.5 0;
sec 7 0 0 -69 0;
sec 8 0 0 -80.5 0;
sec 9 0 0 -92 0;
sec 10 0 0 -103.5 0;
sec 11 0 0 -115.63 0;
end c2_def;
end main_body;
;
begin main_body; towertop
name towertop;
type timoschenko;
nbodies 1;
node_distribution c2_def;
damping_posdef 0 0 0 0.007 0.007 0.007; tuned by Larh
concentrated_mass 2 0 2.687 0.30061 446040 4106000 410600 4106000; nacelle mass and inertia corrected by Anyd 25/4/13
begin timoschenko_input;
filename ./data/DTU_10MW_RWT_Towertop_st.dat;
set 1 2; 1 1=flexibe, 1 2=stiff
end timoschenko_input;
begin c2_def; Definition of centerline (main_body coordinates)
nsec 2;
sec 1 0 0 0 0; x,y,z,twist
sec 2 0 0 -2.75 0;
end c2_def;
end main_body;
;
begin main_body; shaft
name shaft;
type timoschenko;
nbodies 1;
node_distribution c2_def;
damping_posdef 0 0 0 0.000465 0.000465 0.003983; tuned by Anyd 23/5/13 to 31.45 log decr. damping for free free with stiff rotor and tower
concentrated_mass 1 0 0 0 0 0 0 3751000; generator equivalent slow shaft re-tuned by Anyd 20/2/13
concentrated_mass 5 0 0 0 105520 0 0 325700; hub mass and inertia re-tuned by Anyd 20/2/13
begin timoschenko_input;
filename ./data/DTU_10MW_RWT_Shaft_st.dat;
set 1 1; 1 1=flexibe, 1 2=stiff
end timoschenko_input;
begin c2_def; definition of centerline in main_body coordinates
nsec 5;
sec 1 0 0 0 0; x, y, z, twist
sec 2 0 0 1.5 0;
sec 3 0 0 3 0;
sec 4 0 0 4.4 0; main bearing
sec 5 0 0 7.1 0; rotor centre
end c2_def;
end main_body;
;
begin main_body; hub1
name hub1;
type timoschenko;
nbodies 1;
node_distribution c2_def;
damping_posdef 0 0 0 3e-06 3e-06 2e-05; changed by Larh
begin timoschenko_input;
filename ./data/DTU_10MW_RWT_Hub_st.dat;
set 1 2; 1 1=flexibe, 1 2=stiff
end timoschenko_input;
begin c2_def; definition of centerline in main_body coordinates
nsec 2;
sec 1 0 0 0 0; x, y, z, twist
sec 2 0 0 2.8 0;
end c2_def;
end main_body;
;
begin main_body; hub2
name hub2;
copy_main_body hub1;
end main_body;
;
begin main_body; hub3
name hub3;
copy_main_body hub1;
end main_body;
;
begin main_body; blade1
name blade1;
type timoschenko;
nbodies 10;
node_distribution c2_def;
damping_posdef 0 0 0 0.00153 0.00255 0.00033; 3% damping tuned by tkim 23/03/13, unable to fit 3rd and higher mode
begin timoschenko_input;
filename ./data/DTU_10MW_RWT_Blade_st.dat;
set 1 1; 1 1=flexibe, 1 2=stiff
end timoschenko_input;
begin c2_def; definition of centerline in main_body coordinates
nsec 27;
sec 1 0 7.006e-05 4.44089e-16 -14.5;
sec 2 -2.06477e-05 -0.0122119 3 -14.5;
sec 3 -0.0072881 -0.0249251 6 -14.4851;
sec 4 -0.0189235 -0.0273351 7.00004 -14.461;
sec 5 -0.0541282 -0.0282163 8.70051 -14.3388;
sec 6 -0.126633 -0.021321 10.402 -14.0201;
sec 7 -0.225666 -0.0128378 12.2046 -13.3904;
sec 8 -0.288563 -0.00770659 13.2065 -12.9371;
sec 9 -0.399194 -0.00488317 15.01 -11.9445;
sec 10 -0.576634 -0.0180296 18.2151 -9.98243;
sec 11 -0.707136 -0.0501772 21.4178 -8.45147;
sec 12 -0.791081 -0.0941228 24.6189 -7.46417;
sec 13 -0.837195 -0.14888 27.8193 -6.72916;
sec 14 -0.853948 -0.214514 31.0194 -6.08842;
sec 15 -0.849367 -0.290618 34.2197 -5.49322;
sec 16 -0.79392 -0.462574 40.2204 -4.39222;
sec 17 -0.716284 -0.688437 46.6217 -3.09315;
sec 18 -0.634358 -0.960017 53.0232 -1.75629;
sec 19 -0.553179 -1.28424 59.4245 -0.50065;
sec 20 -0.475422 -1.66402 65.8255 0.601964;
sec 21 -0.40318 -2.10743 72.2261 1.5556;
sec 22 -0.330085 -2.6563 79.0266 2.51935;
sec 23 -0.31014 -2.78882 80.5267 2.7295;
sec 24 -0.286719 -2.92517 82.0271 2.93201;
sec 25 -0.255823 -3.06577 83.5274 3.11874;
sec 26 -0.207891 -3.20952 85.0277 3.28847;
sec 27 -0.089894 -3.33685 86.3655 3.42796;
end c2_def;
end main_body;
;
begin main_body; blade2
name blade2;
copy_main_body blade1;
end main_body;
;
begin main_body; blade3
name blade3;
copy_main_body blade1;
end main_body;
;
;-------------------------------------------------------------------------------------------------------------------------------
begin orientation; orientations of body coordinate systems
begin base; tower to ground
body tower;
inipos 0 0 0; initial position of node 1
body_eulerang 0 0 0; same as global; zT down and yT downwind.
end base;
;
begin relative; towertop to tower
body1 tower last;
body2 towertop 1;
body2_eulerang 0 0 0; same as tower; zTT down and yTT downwind.
end relative;
;
begin relative; shaft to towertop
body1 towertop last;
body2 shaft 1;
body2_eulerang 90 0 0;
body2_eulerang 5 0 0; 5 deg tilt angle; zS directed upwind along shaft, xS parallel to ground
mbdy2_ini_rotvec_d1 0 0 -1 0.2; initial shaft rotation
end relative;
;
begin relative; hub1 to shaft
body1 shaft last;
body2 hub1 1;
body2_eulerang -90 0 0;
body2_eulerang 0 180 0;
body2_eulerang 2.5 0 0; 2.5deg cone angle; yH downwind along shaft, zH outward towards blade tip
end relative;
;
begin relative; hub2 to shaft
body1 shaft last;
body2 hub2 1;
body2_eulerang -90 0 0;
body2_eulerang 0 60 0;
body2_eulerang 2.5 0 0; 2.5deg cone angle; yH downwind along shaft, zH outward towards blade tip
end relative;
;
begin relative; hub3 to shaft
body1 shaft last;
body2 hub3 1;
body2_eulerang -90 0 0;
body2_eulerang 0 -60 0;
body2_eulerang 2.5 0 0; 2.5deg cone angle; yH downwind along shaft, zH outward towards blade tip
end relative;
;
begin relative; blade1 to hub1
body1 hub1 last;
body2 blade1 1;
body2_eulerang 0 0 0; same as hub1; yB downwind along shaft, zB outward towards blade tip
end relative;
;
begin relative; blade2 to hub2
body1 hub2 last;
body2 blade2 1;
body2_eulerang 0 0 0; same as hub2; yB downwind along shaft, zB outward towards blade tip
end relative;
;
begin relative; blade3 to hub3
body1 hub3 last;
body2 blade3 1;
body2_eulerang 0 0 0; same as hub3; yB downwind along shaft, zB outward towards blade tip
end relative;
;
end orientation;
;
;-------------------------------------------------------------------------------------------------------------------------------
begin constraint; constraints of each body to another
;
begin fix0; tower fixed to ground in translation and rotation of node 1
body tower;
end fix0;
;
begin fix1; bottom of towertop rigidly fixed to top of tower
body1 tower last;
body2 towertop 1;
end fix1;
;
begin bearing1; shaft_rot: shaft rotates about its -z axis with respect to towertop
name shaft_rot;
body1 towertop last;
body2 shaft 1;
bearing_vector 2 0 0 -1;
end bearing1;
;
begin fix1; bottom of hub1 rigidly fixed to shaft
body1 shaft last;
body2 hub1 1;
end fix1;
;
begin fix1; bottom of hub2 rigidly fixed to shaft
body1 shaft last;
body2 hub2 1;
end fix1;
;
begin fix1; bottom of hub3 rigidly fixed to shaft
body1 shaft last;
body2 hub3 1;
end fix1;
;
begin bearing2; pitch1: blade1 rotates about its -z axis with respect to hub
name pitch1;
body1 hub1 last;
body2 blade1 1;
bearing_vector 2 0 0 -1;
end bearing2;
;
begin bearing2; pitch2: blade2 rotates about its -z axis with respect to hub
name pitch2;
body1 hub2 last;
body2 blade2 1;
bearing_vector 2 0 0 -1;
end bearing2;
;
begin bearing2; pitch3: blade3 rotates about its -z axis with respect to hub
name pitch3;
body1 hub3 last;
body2 blade3 1;
bearing_vector 2 0 0 -1;
end bearing2;
end constraint;
;
end new_htc_structure;
;
;----------------------------------------------------------------------------------------------------------------------------------------------------------------
begin wind; wind input to structure
density 1.225;
wsp 24;
tint 0.157; iec: i_ref*(0.75*U+5.6)/U
horizontal_input 1;
windfield_rotations 0 0 0; yaw, tilt, rotation
center_pos0 0 0 -119; hub height
shear_format 3 0.2; 0=no mean wind, 1=constant, 2=log, 3=power, 4=linear
turb_format 0;
tower_shadow_method 0;
wind_ramp_factor 0 40 0.5 1; ramp wind from 50% to 100% starting at 0s and ending at 40s
;
end wind;
;
begin aero; aerodynamic forces
nblades 3;
hub_vec shaft -3; vector of rotor rotation direction
link 1 mbdy_c2_def blade1;
link 2 mbdy_c2_def blade2;
link 3 mbdy_c2_def blade3;
ae_filename ./data/DTU_10MW_RWT_ae.dat;
pc_filename ./data/DTU_10MW_RWT_pc.dat;
induction_method 1; 0=none, 1=normal
aerocalc_method 1; 0=no aero calculations, 1=with aero calculations
aerosections 50;
ae_sets 1 1 1;
tiploss_method 1; 0=none, 1=prandtl
dynstall_method 2; 0=none, 1=stig øye method, 2=mhh method
;
end aero;
;-------------------------------------------------------------------------------------------------
begin hawcstab2;
; define structure;
begin ground_fixed_substructure;
main_body tower;
main_body towertop;
end ground_fixed_substructure;
begin rotating_axissym_substructure;
main_body shaft;
end rotating_axissym_substructure;
begin rotating_threebladed_substructure;
main_body hub1;
main_body blade1;
second_order_actuator pitch1 100 0.7;
end rotating_threebladed_substructure;
operational_data_filename ./data/operation.dat;
; inputs for finding optimal operational data;
begin operational_data;
windspeed 4 25 22; cut-in [m/s], cut-out [m/s], points [-]
genspeed 299.85 479.56567; [rpm]
gearratio 50; [-]
minpitch 0; [deg]
opt_lambda 7.5; [-]
maxpow 10638.297872340427; [kW]
prvs_turbine 1; [-]
include_torsiondeform 1; [-]
end operational_data;
; define controller couplings;
begin controller;
begin input;
constraint bearing1 shaft_rot;
constraint bearing2 pitch1 collective;
end input;
begin output;
constraint bearing1 shaft_rot 1 only 2;
constraint bearing2 pitch1 1 only 1 collective;
end output;
end controller;
; inputs for controller tuning;
begin controller_tuning;
partial_load 0.05 0.7; fn [hz], zeta [-]
full_load 0.06 0.7; fn [hz], zeta [-]
gain_scheduling 2; 1 linear, 2 quadratic
constant_power 1; 0 constant torque, 1 constant power at full load
end controller_tuning;
; HAWC2S commands (uncomment as needed);
;compute_optimal_pitch_angle use_operational_data;
compute_steady_states bladedeform tipcorrect induction nogradients; compute steady states using hawcstab2 (need for other commands)
save_power; save steady-state values to input_htc_file.pwr
compute_structural_modal_analysis nobladeonly 12; struct modal analysis saved to <htc>_struc.cmb
compute_stability_analysis windturbine 12; aeroelastic modal analysis saved to <htc>.cmb
;compute_controller_input; tune controller parameters
end hawcstab2;
exit;
\ No newline at end of file
This diff is collapsed.
Supports Markdown
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