Skip to content
Snippets Groups Projects
Commit 60947a7d authored by Jenni Rinker's avatar Jenni Rinker
Browse files

adding version 1.0 files

parent b2162fe2
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
nrel-5mw-rwt/data/NREL_5MW_ae.txt 0 → 100644
+ 21
− 0
View file @ 60947a7d
1 Chord[m] T/C[%] Set no.
1 18 Blade
0.000 3.542 100 1
1.367 3.542 100 1
4.100 3.854 100 1
6.833 4.167 90 1
10.250 4.557 40 1
14.350 4.652 35 1
18.450 4.458 35 1
22.550 4.249 30 1
26.650 4.007 25 1
30.750 3.748 25 1
34.850 3.502 21 1
38.950 3.256 21 1
43.050 3.01 17 1
47.150 2.764 17 1
51.250 2.518 17 1
54.667 2.313 17 1
57.400 2.086 17 1
60.133 1.419 17 1
61.500 0.961 17 1
nrel-5mw-rwt/data/NREL_5MW_pc.txt 0 → 100644
+ 840
− 0
View file @ 60947a7d
This diff is collapsed.
nrel-5mw-rwt/data/NREL_5MW_st.txt 0 → 100644
+ 272
− 0
View file @ 60947a7d
This diff is collapsed.
nrel-5mw-rwt/htc/NREL_5MW_reference_wind_turbine.htc 0 → 100644
+ 468
− 0
View file @ 60947a7d
begin Simulation;
time_stop 100;
solvertype 1 ; (newmark)
on_no_convergence continue ;
logfile ./log/NREL_5MW_reference_wind_turbine.log ;
animation ./animation/NREL_5MW_reference_wind_turbine.dat;
;
begin newmark;
deltat 0.02;
end newmark;
end simulation;
;
begin new_htc_structure;
beam_output_file_name ./log/NREL_5MW_reference_wind_turbine_beam.dat; Optional - Calculated beam properties of the bodies are written to file
body_output_file_name ./log/NREL_5MW_reference_wind_turbine_body.dat; Optional - Body initial position and orientation are written to file
; body_eigenanalysis_file_name ./eigenfrq/NREL_5MW_reference_wind_turbine_body_eigen.dat;
; structure_eigenanalysis_file_name ./eigenfrq/NREL_5MW_reference_wind_turbine_strc_eigen.dat ;
;-------------------------------------------------------------------------------------------------------------------------------
begin main_body; tower 90m
name tower ;
type timoschenko ;
nbodies 1 ;
node_distribution c2_def ;
damping_posdef 8.142E-04 8.14E-04 4.0E-03 4.3E-04 4.3E-04 4.3E-03 ; Mx My Mz Kx Ky Kz , Ms raises overall level, Ks raises high freguency level
begin timoschenko_input;
filename ./data/NREL_5MW_st.txt ;
set 1 1 ;
end timoschenko_input;
begin c2_def; Definition of centerline (main_body coordinates)
nsec 8;
sec 1 0.0 0.0 0.0 0.0 ; x,y,z,twist
sec 2 0.0 0.0 -10.0 0.0 ;
sec 3 0.0 0.0 -20.0 0.0 ;
sec 4 0.0 0.0 -30.0 0.0 ;
sec 5 0.0 0.0 -45.0 0.0 ;
sec 6 0.0 0.0 -60.0 0.0 ;
sec 7 0.0 0.0 -70.0 0.0 ;
sec 8 0.0 0.0 -87.6 0.0 ;
end c2_def ;
end main_body;
;
begin main_body;
name towertop ;
type timoschenko ;
nbodies 1 ;
node_distribution c2_def ;
; damping_posdef 9.025E-06 9.025E-06 8.0E-05 8.3E-06 8.3E-06 8.5E-05 ;
damping 2.50E-04 1.40E-04 2.00E-03 3.00E-05 3.00E-05 2.00E-04 ;
concentrated_mass 2 0.0 1.9 0.21256 2.4E5 1741490.0 1.7E5 1741490.0 ; Nacelle mass and inertia
begin timoschenko_input;
filename ./data/NREL_5MW_st.txt ;
set 2 1 ;
end timoschenko_input;
begin c2_def; Definition of centerline (main_body coordinates)
nsec 2;
sec 1 0.0 0.0 0.0 0.0 ; x,y,z,twist
sec 2 0.0 0.0 -1.96256 0.0 ;
end c2_def ;
end main_body;
;
begin main_body;
name shaft ;
type timoschenko ;
nbodies 1 ;
node_distribution c2_def ;
; damping_posdef 7.00E-3 7.00E-03 7.00E-02 3.48E-04 3.48E-04 1.156E-03 ;
damping_posdef 7.00E-3 7.00E-03 7.00E-02 6.5E-04 6.5E-04 1.84E-02 ;
concentrated_mass 1 0.0 0.0 0.0 0.0 0.0 0.0 5025497.444 ;generator equivalent slow shaft
concentrated_mass 5 0.0 0.0 0.0 56780 0.0 0.0 115926 ; hub mass and inertia;
begin timoschenko_input;
filename ./data/NREL_5MW_st.txt ;
set 3 1 ;
end timoschenko_input;
begin c2_def; Definition of centerline (main_body coordinates)
nsec 5;
sec 1 0.0 0.0 0.0 0.0 ; Tower top x,y,z,twist
sec 2 0.0 0.0 1.0 0.0 ;
sec 3 0.0 0.0 2.0 0.0 ;
sec 4 0.0 0.0 3.1071 0.0 ; Main bearing
sec 5 0.0 0.0 5.0191 0.0 ; Rotor centre
end c2_def ;
end main_body;
;
begin main_body;
name hub1 ;
type timoschenko ;
nbodies 1 ;
node_distribution c2_def ;
damping_posdef 2.00E-05 2.00E-05 2.00E-04 3.00E-06 3.00E-06 2.00E-05;
begin timoschenko_input;
filename ./data/NREL_5MW_st.txt ;
set 4 1 ;
end timoschenko_input;
begin c2_def; Definition of centerline (main_body coordinates)
nsec 2;
sec 1 0.0 0.0 0.0 0.0 ; x,y,z,twist
sec 2 0.0 0.0 1.5 0.0 ;
end c2_def ;
end main_body;
;
begin main_body;
name hub2 ;
copy_main_body hub1;
end main_body;
;
begin main_body;
name hub3 ;
copy_main_body hub1 ;
end main_body;
;
begin main_body;
name blade1 ;
type timoschenko ;
nbodies 1 ;
node_distribution c2_def;
; damping 3.5e-2 5.5e-4 5.0e-4 3.0e-4 0.5e-3 5.5e-3 ;
damping_posdef 1.16e-4 5.75e-5 6.1e-6 6.5e-4 5.1e-4 6.4e-4 ;
begin timoschenko_input ;
filename ./data/NREL_5MW_st.txt ;
set 5 1 ; set subset
end timoschenko_input;
begin c2_def; Definition of centerline (main_body coordinates)
nsec 19 ;
sec 1 0.0000 0.0000 0.000 0.000 ; x.y.z. twist
sec 2 -0.0041 0.0010 1.367 -13.308 ;
sec 3 -0.1058 0.0250 4.100 -13.308 ;
sec 4 -0.2502 0.0592 6.833 -13.308 ;
sec 5 -0.4594 0.1087 10.250 -13.308 ;
sec 6 -0.5699 0.1157 14.350 -11.480 ;
sec 7 -0.5485 0.0983 18.450 -10.162 ;
sec 8 -0.5246 0.0832 22.550 -9.011 ;
sec 9 -0.4962 0.0679 26.650 -7.795 ;
sec 10 -0.4654 0.0534 30.750 -6.544 ; 50% blade radius
sec 11 -0.4358 0.0409 34.850 -5.361 ;
sec 12 -0.4059 0.0297 38.950 -4.188 ;
sec 13 -0.3757 0.0205 43.050 -3.125 ;
sec 14 -0.3452 0.0140 47.150 -2.319 ;
sec 15 -0.3146 0.0084 51.250 -1.526 ;
sec 16 -0.2891 0.0044 54.667 -0.863 ;
sec 17 -0.2607 0.0017 57.400 -0.370 ;
sec 18 -0.1774 0.0003 60.133 -0.106 ;
sec 19 -0.1201 0.0000 61.500 -0.000 ;
end c2_def ;
end main_body;
;
begin main_body;
name blade2 ;
copy_main_body blade1;
end main_body;
;
begin main_body;
name blade3 ;
copy_main_body blade1 ;
end main_body;
;-------------------------------------------------------------------------------------------------------------------------------
;
begin orientation;
begin base;
body tower;
inipos 0.0 0.0 0.0 ; initial position of node 1
body_eulerang 0.0 0.0 0.0;
end base;
;
begin relative;
body1 tower last;
body2 towertop 1;
body2_eulerang 0.0 0.0 0.0;
end relative;
;
begin relative;
body1 towertop last;
body2 shaft 1;
body2_eulerang 90.0 0.0 0.0;
body2_eulerang 5.0 0.0 0.0; 5 deg tilt angle
body2_ini_rotvec_d1 0.0 0.0 -1.0 0.5 ; body initial rotation velocity x.y.z.angle velocity[rad/s] (body 2 coordinates)
; body2_ini_rotvec_d1 0.0 0.0 -1.0 0.9424 ; body initial rotation velocity x.y.z.angle velocity[rad/s] (body 2 coordinates)
end relative;
;
begin relative;
body1 shaft last;
body2 hub1 1;
body2_eulerang -90.0 0.0 0.0;
body2_eulerang 0.0 180.0 0.0;
body2_eulerang 2.5 0.0 0.0; 2.5deg cone angle
end relative;
;
begin relative;
body1 shaft last;
body2 hub2 1;
body2_eulerang -90.0 0.0 0.0;
body2_eulerang 0.0 60.0 0.0;
body2_eulerang 2.5 0.0 0.0; 2.5deg cone angle
end relative;
;
begin relative;
body1 shaft last;
body2 hub3 1;
body2_eulerang -90.0 0.0 0.0;
body2_eulerang 0.0 -60.0 0.0;
body2_eulerang 2.5 0.0 0.0; 2.5deg cone angle
end relative;
;
begin relative;
body1 hub1 last;
body2 blade1 1;
body2_eulerang 0.0 0.0 0;
end relative;
;
begin relative;
body1 hub2 last;
body2 blade2 1;
body2_eulerang 0.0 0.0 0.0;
end relative;
;
begin relative;
body1 hub3 last;
body2 blade3 1;
body2_eulerang 0.0 0.0 0.0;
end relative;
;
end orientation;
;-------------------------------------------------------------------------------------------------------------------------------
begin constraint;
;
begin fix0; fixed to ground in translation and rotation of node 1
body tower;
end fix0;
;
begin fix1;
body1 tower last ;
body2 towertop 1;
end fix1;
;
begin bearing1; free bearing
name shaft_rot;
body1 towertop last;
body2 shaft 1;
bearing_vector 2 0.0 0.0 -1.0; x=coo (0=global.1=body1.2=body2) vector in body2 coordinates where the free rotation is present
end bearing1;
;
begin fix1;
body1 shaft last ;
body2 hub1 1;
end fix1;
;
begin fix1;
body1 shaft last ;
body2 hub2 1;
end fix1;
;
begin fix1;
body1 shaft last ;
body2 hub3 1;
end fix1;
;
begin bearing2;
name pitch1;
body1 hub1 last;
body2 blade1 1;
bearing_vector 2 0.0 0.0 -1.0;
end bearing2;
;
begin bearing2;
name pitch2;
body1 hub2 last;
body2 blade2 1;
bearing_vector 2 0.0 0.0 -1.0;
end bearing2;
;
begin bearing2;
name pitch3;
body1 hub3 last;
body2 blade3 1;
bearing_vector 2 0.0 0.0 -1.0;
end bearing2;
end constraint;
;
end new_htc_structure;
;----------------------------------------------------------------------------------------------------------------------------------------------------------------
begin wind ;
density 1.225 ; to be checked
wsp 8.0 ;
tint 0.0 ;
horizontal_input 1 ; 0=false, 1=true
windfield_rotations 0.0 0.0 0.0 ; yaw, tilt, rotation
center_pos0 0.0 0.0 -90.0 ;
shear_format 1 0.14 ;0=none,1=constant,2=log,3=power,4=linear
turb_format 0 ; 0=none, 1=mann,2=flex
tower_shadow_method 3 ; 0=none, 1=potential flow, 2=jet
; scale_time_start 0 ;
; wind_ramp_factor 0.0 [t0] [wsp factor] 1.0 ;
; [gust] iec_gust [gust_type] [gust_A] [gust_phi0] [gust_t0] [gust_T] ;
;
begin mann;
filename_u ./turb/dummyu.bin ;
filename_v ./turb/dummyv.bin ;
filename_w ./turb/dummyw.bin ;
box_dim_u 8192 0.0 ;
box_dim_v 32 4.6875;
box_dim_w 32 4.6875;
std_scaling 1.0 0.7 0.5 ;
end mann;
;
begin tower_shadow_potential_2;
tower_mbdy_link tower;
nsec 2;
radius 0.0 3.0 ;
radius 90.0 1.935 ;
end tower_shadow_potential_2;
end wind;
;
begin aero ;
nblades 3;
hub_vec shaft -3 ; rotor rotation vector (normally shaft composant directed from pressure to sustion side)
link 1 mbdy_c2_def blade1;
link 2 mbdy_c2_def blade2;
link 3 mbdy_c2_def blade3;
ae_filename ./data/NREL_5MW_ae.txt;
pc_filename ./data/NREL_5MW_pc.txt;
induction_method 1 ; 0=none, 1=normal
aerocalc_method 1 ; 0=ingen aerodynamic, 1=med aerodynamic
aerosections 30 ;
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 dll;
begin hawc_dll;
filename ./control/bladed2hawc.dll ;
dll_subroutine regulation ;
arraysizes 15 15 ;
deltat 0.02;
begin output;
general time ; 1
constraint bearing2 pitch1 1; angle and angle velocity written to dll 2,3
constraint bearing2 pitch2 1; angle and angle velocity written to dll 4,5
constraint bearing2 pitch3 1; angle and angle velocity written to dll 6,7
constraint bearing2 shaft_rot 1; angle and angle velocity written to dll (slow speed shaft) 8,9
wind free_wind 1 0.0 0.0 -90.0; local wind at fixed position: coo (1=global,2=non-rotation rotor coo.), pos x, pos y, pos z 10,11,12
general constant 97.0 ; generator exchange ratio 13
end output;
;
begin actions;
body moment_int shaft 1 3 towertop 2 ;
end actions;
end hawc_dll;
;
begin hawc_dll;
filename ./control/pitchservo_pos.dll ;
dll_subroutine servo ;
arraysizes 15 15 ;
deltat 0.02 ;
begin output;
general time ; 1
dll inpvec 1 2; 2
dll inpvec 1 3; 3
dll inpvec 1 4; 4
constraint bearing2 pitch1 1; angle and angle velocity written to dll 5,6
constraint bearing2 pitch2 1; angle and angle velocity written to dll 7,8
constraint bearing2 pitch3 1; angle and angle velocity written to dll 9,10
end output;
;
begin actions;
body bearing_angle pitch1;
body bearing_angle pitch2;
body bearing_angle pitch3;
end actions;
end hawc_dll;
;
begin hawc_dll;
filename ./control/damper.dll ;
dll_subroutine damp ;
arraysizes 15 15 ;
begin output;
general time ; 1
general constant 5.0;
general constant 10.0;
general constant -1.0E1 ;
mbdy state vel towertop 1 1.0 tower;
end output;
;
begin actions;
mbdy force_ext towertop 2 1 towertop;
mbdy force_ext towertop 2 2 towertop;
end actions;
end hawc_dll;
end dll;
;
;----------------------------------------------------------------------------------------------------------------------------------------------------------------
;
begin output;
filename ./res/NREL_5MW_reference_wind_turbine ;
; time 390.0 450.0 ;
buffer 1 ;
general time;
data_format hawc_binary;
;
constraint bearing1 shaft_rot 2; angle and angle velocity
constraint bearing2 pitch1 5; angle and angle velocity
constraint bearing2 pitch2 5; angle and angle velocity
constraint bearing2 pitch3 5; angle and angle velocity
aero omega ;
aero torque;
aero power;
aero thrust;
wind free_wind 1 0.0 0.0 -90.0; local wind at fixed position: coo (1=global,2=non-rotation rotor coo.), pos x, pos y, pos z
mbdy momentvec towertop 1 2 towertop # yaw bearing ;
mbdy forcevec towertop 1 2 towertop # yaw bering ;
mbdy momentvec shaft 4 1 shaft # main bearing ;
mbdy momentvec blade1 3 1 blade1 # blade 1 root ;
mbdy momentvec blade1 10 1 local # blade 1 50% local e coo ;
mbdy momentvec hub1 1 2 hub1 # blade 1 root ;
mbdy momentvec hub2 1 2 hub2 # blade 2 root ;
mbdy momentvec hub3 1 2 hub3 # blade 3 root ;
mbdy state pos towertop 1 1.0 global # tower top flange position ;
mbdy state pos blade1 18 1.0 blade1 # blade 1 tip pos ;
mbdy state pos blade2 18 1.0 blade2 # blade 2 tip pos ;
mbdy state pos blade3 18 1.0 blade3 # blade 3 tip pos ;
mbdy state pos blade1 18 1.0 global # blade 1 tip pos ;
aero windspeed 3 1 1 63.0; wind seen from the blade: coo(1=local ae,2=blade,3=global,4=rotor polar), blade number, component xyz, blade radius
aero windspeed 3 1 2 63.0;
aero windspeed 3 1 3 63.0;
aero alfa 1 45.0;
aero alfa 2 45.0;
aero alfa 3 45.0;
mbdy momentvec towertop 1 1 tower # tower top -1: below top mass ;
mbdy forcevec towertop 1 1 tower # tower top -1: below top mass ;
mbdy momentvec tower 1 1 tower # tower base ;
mbdy forcevec tower 1 1 tower # tower base ;
;
dll outvec 1 1 # time;
dll outvec 1 2 # pitch angle 1;
dll outvec 1 3 # pitch vel 1;
dll outvec 1 4 # pitch angle 2;
dll outvec 1 5 # pitch vel 2;
dll outvec 1 6 # pitch angle 3;
dll outvec 1 7 # pitch vel 3;
dll outvec 1 8 # gen. azi slow;
dll outvec 1 9 # gen. speed slow;
dll outvec 1 10 # free wind x;
dll outvec 1 11 # free wind y;
dll outvec 1 12 # free wind z;
dll outvec 1 13 # gear ratio;
dll inpvec 1 1 # Mgen slow;
dll inpvec 1 2 # pitchref 1;
dll inpvec 1 3 # pitchref 2;
dll inpvec 1 4 # pitchref 3;
dll inpvec 1 7 # F;
dll inpvec 1 8 # Mechanical power generator [kW];
dll inpvec 1 10 # Pitch rate [rad/s];
dll inpvec 2 1 # pitch 1;
dll inpvec 2 2 # pitch 2;
dll inpvec 2 3 # pitch 3;
dll outvec 2 1 # time;
dll outvec 2 2 # pitchref 1;
dll outvec 2 3 # pitchref 2;
dll outvec 2 4 # pitchref 3;
dll outvec 2 5 # pitch angle 1;
dll outvec 2 6 # pitch speed 1;
dll outvec 2 7 # pitch angle 2;
dll outvec 2 8 # pitch speed 2;
dll outvec 2 9 # pitch angle 3;
dll outvec 2 10 # pitch speed 3;
end output;
;
exit;
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