Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found

Target

Select target project
  • toolbox/WindEnergyToolbox
  • tlbl/WindEnergyToolbox
  • cpav/WindEnergyToolbox
  • frza/WindEnergyToolbox
  • borg/WindEnergyToolbox
  • mmpe/WindEnergyToolbox
  • ozgo/WindEnergyToolbox
  • dave/WindEnergyToolbox
  • mmir/WindEnergyToolbox
  • wluo/WindEnergyToolbox
  • welad/WindEnergyToolbox
  • chpav/WindEnergyToolbox
  • rink/WindEnergyToolbox
  • shfe/WindEnergyToolbox
  • shfe1/WindEnergyToolbox
  • acdi/WindEnergyToolbox
  • angl/WindEnergyToolbox
  • wliang/WindEnergyToolbox
  • mimc/WindEnergyToolbox
  • wtlib/WindEnergyToolbox
  • cmos/WindEnergyToolbox
  • fabpi/WindEnergyToolbox
22 results
Show changes
Hide whitespace changes
Inline Side-by-side
Showing
with 2677 additions and 378 deletions
wetb/hawc2/tests/test_files/simulation_setup/DTU10MWRef6.0_IOS/input/data/DTU_10MW_RWT_Shaft_st.dat 0 → 100644
View file @ 62f0f2e4
#1 Shaft
r m x_cg y_cg ri_x ri_y x_sh y_sh E G I_x I_y I_p k_x k_y A pitch x_e y_e
$1 2 shaft with hub mass on last element. flexible
0.0 1.00E+00 0 0 0.2 0.2 0 0 2.10E+11 8.08E+10 1.00E+00 1.00E+00 0.2036 0.5 0.5 10 0 0 0
7.1 1.00E+00 0 0 0.2 0.2 0 0 2.10E+11 8.08E+10 1.00E+00 1.00E+00 0.2036 0.5 0.5 10 0 0 0
r m x_cg y_cg ri_x ri_y x_sh y_sh E G I_x I_y I_p k_x k_y A pitch x_e y_e
$2 2 shaft with hub mass on last element. stiff
0.0 1.00E+00 0 0 0.2 0.2 0 0 2.10E+16 8.08E+17 1.00E+00 1.00E+00 0.2036 0.5 0.5 10 0 0 0
7.1 1.00E+00 0 0 0.2 0.2 0 0 2.10E+16 8.08E+17 1.00E+00 1.00E+00 0.2036 0.5 0.5 10 0 0 0
wetb/hawc2/tests/test_files/simulation_setup/DTU10MWRef6.0_IOS/input/data/DTU_10MW_RWT_Tower_st.dat 0 → 100644
View file @ 62f0f2e4
#1 Tower made by anyd 25.02.2013
r m x_cg y_cg ri_x ri_y x_sh y_sh E G I_x I_y I_p k_x k_y A pitch x_e y_e
$1 20
0.000 8.3837E+03 0.0000E+00 0.0000E+00 2.9211E+00 2.9211E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 8.4160E+00 8.4160E+00 1.6832E+01 5.0000E-01 5.0000E-01 9.8632E-01 0.0000E+00 0.0000E+00 0.0000E+00
11.500 8.1012E+03 0.0000E+00 0.0000E+00 2.8226E+00 2.8226E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 7.5934E+00 7.5934E+00 1.5187E+01 5.0000E-01 5.0000E-01 9.5308E-01 0.0000E+00 0.0000E+00 0.0000E+00
11.501 7.6767E+03 0.0000E+00 0.0000E+00 2.8233E+00 2.8233E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 7.1991E+00 7.1991E+00 1.4398E+01 5.0000E-01 5.0000E-01 9.0314E-01 0.0000E+00 0.0000E+00 0.0000E+00
23.000 7.4090E+03 0.0000E+00 0.0000E+00 2.7249E+00 2.7249E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 6.4720E+00 6.4720E+00 1.2944E+01 5.0000E-01 5.0000E-01 8.7165E-01 0.0000E+00 0.0000E+00 0.0000E+00
23.001 6.9992E+03 0.0000E+00 0.0000E+00 2.7256E+00 2.7256E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 6.1171E+00 6.1171E+00 1.2234E+01 5.0000E-01 5.0000E-01 8.2343E-01 0.0000E+00 0.0000E+00 0.0000E+00
34.500 6.7464E+03 0.0000E+00 0.0000E+00 2.6271E+00 2.6271E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 5.4779E+00 5.4779E+00 1.0956E+01 5.0000E-01 5.0000E-01 7.9369E-01 0.0000E+00 0.0000E+00 0.0000E+00
34.501 6.3512E+03 0.0000E+00 0.0000E+00 2.6278E+00 2.6278E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 5.1598E+00 5.1598E+00 1.0320E+01 5.0000E-01 5.0000E-01 7.4720E-01 0.0000E+00 0.0000E+00 0.0000E+00
46.000 6.1133E+03 0.0000E+00 0.0000E+00 2.5294E+00 2.5294E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 4.6013E+00 4.6013E+00 9.2027E+00 5.0000E-01 5.0000E-01 7.1921E-01 0.0000E+00 0.0000E+00 0.0000E+00
46.001 5.7328E+03 0.0000E+00 0.0000E+00 2.5301E+00 2.5301E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 4.3173E+00 4.3173E+00 8.6346E+00 5.0000E-01 5.0000E-01 6.7444E-01 0.0000E+00 0.0000E+00 0.0000E+00
57.500 5.5097E+03 0.0000E+00 0.0000E+00 2.4316E+00 2.4316E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 3.8327E+00 3.8327E+00 7.6654E+00 5.0000E-01 5.0000E-01 6.4820E-01 0.0000E+00 0.0000E+00 0.0000E+00
57.501 5.1439E+03 0.0000E+00 0.0000E+00 2.4323E+00 2.4323E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 3.5803E+00 3.5803E+00 7.1605E+00 5.0000E-01 5.0000E-01 6.0516E-01 0.0000E+00 0.0000E+00 0.0000E+00
69.000 4.9357E+03 0.0000E+00 0.0000E+00 2.3339E+00 2.3339E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 3.1629E+00 3.1629E+00 6.3258E+00 5.0000E-01 5.0000E-01 5.8067E-01 0.0000E+00 0.0000E+00 0.0000E+00
69.001 4.5845E+03 0.0000E+00 0.0000E+00 2.3346E+00 2.3346E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 2.9396E+00 2.9396E+00 5.8792E+00 5.0000E-01 5.0000E-01 5.3935E-01 0.0000E+00 0.0000E+00 0.0000E+00
80.500 4.3912E+03 0.0000E+00 0.0000E+00 2.2361E+00 2.2361E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 2.5832E+00 2.5832E+00 5.1664E+00 5.0000E-01 5.0000E-01 5.1661E-01 0.0000E+00 0.0000E+00 0.0000E+00
80.501 4.0547E+03 0.0000E+00 0.0000E+00 2.2368E+00 2.2368E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 2.3867E+00 2.3867E+00 4.7734E+00 5.0000E-01 5.0000E-01 4.7702E-01 0.0000E+00 0.0000E+00 0.0000E+00
92.000 3.8762E+03 0.0000E+00 0.0000E+00 2.1384E+00 2.1384E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 2.0852E+00 2.0852E+00 4.1705E+00 5.0000E-01 5.0000E-01 4.5602E-01 0.0000E+00 0.0000E+00 0.0000E+00
92.001 3.5543E+03 0.0000E+00 0.0000E+00 2.1391E+00 2.1391E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 1.9133E+00 1.9133E+00 3.8267E+00 5.0000E-01 5.0000E-01 4.1816E-01 0.0000E+00 0.0000E+00 0.0000E+00
103.500 3.3908E+03 0.0000E+00 0.0000E+00 2.0406E+00 2.0406E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 1.6611E+00 1.6611E+00 3.3223E+00 5.0000E-01 5.0000E-01 3.9891E-01 0.0000E+00 0.0000E+00 0.0000E+00
103.501 3.0836E+03 0.0000E+00 0.0000E+00 2.0413E+00 2.0413E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 1.5117E+00 1.5117E+00 3.0234E+00 5.0000E-01 5.0000E-01 3.6277E-01 0.0000E+00 0.0000E+00 0.0000E+00
115.630 2.9267E+03 0.0000E+00 0.0000E+00 1.9375E+00 1.9375E+00 0.0000E+00 0.0000E+00 2.1000E+11 8.0769E+10 1.2925E+00 1.2925E+00 2.5850E+00 5.0000E-01 5.0000E-01 3.4432E-01 0.0000E+00 0.0000E+00 0.0000E+00
$2 20
0.000 8.3837E+03 0.0000E+00 0.0000E+00 2.9211E+00 2.9211E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 8.4160E+00 8.4160E+00 1.6832E+01 5.0000E-01 5.0000E-01 9.8632E-01 0.0000E+00 0.0000E+00 0.0000E+00
11.500 8.1012E+03 0.0000E+00 0.0000E+00 2.8226E+00 2.8226E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 7.5934E+00 7.5934E+00 1.5187E+01 5.0000E-01 5.0000E-01 9.5308E-01 0.0000E+00 0.0000E+00 0.0000E+00
11.501 7.6767E+03 0.0000E+00 0.0000E+00 2.8233E+00 2.8233E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 7.1991E+00 7.1991E+00 1.4398E+01 5.0000E-01 5.0000E-01 9.0314E-01 0.0000E+00 0.0000E+00 0.0000E+00
23.000 7.4090E+03 0.0000E+00 0.0000E+00 2.7249E+00 2.7249E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 6.4720E+00 6.4720E+00 1.2944E+01 5.0000E-01 5.0000E-01 8.7165E-01 0.0000E+00 0.0000E+00 0.0000E+00
23.001 6.9992E+03 0.0000E+00 0.0000E+00 2.7256E+00 2.7256E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 6.1171E+00 6.1171E+00 1.2234E+01 5.0000E-01 5.0000E-01 8.2343E-01 0.0000E+00 0.0000E+00 0.0000E+00
34.500 6.7464E+03 0.0000E+00 0.0000E+00 2.6271E+00 2.6271E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 5.4779E+00 5.4779E+00 1.0956E+01 5.0000E-01 5.0000E-01 7.9369E-01 0.0000E+00 0.0000E+00 0.0000E+00
34.501 6.3512E+03 0.0000E+00 0.0000E+00 2.6278E+00 2.6278E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 5.1598E+00 5.1598E+00 1.0320E+01 5.0000E-01 5.0000E-01 7.4720E-01 0.0000E+00 0.0000E+00 0.0000E+00
46.000 6.1133E+03 0.0000E+00 0.0000E+00 2.5294E+00 2.5294E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 4.6013E+00 4.6013E+00 9.2027E+00 5.0000E-01 5.0000E-01 7.1921E-01 0.0000E+00 0.0000E+00 0.0000E+00
46.001 5.7328E+03 0.0000E+00 0.0000E+00 2.5301E+00 2.5301E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 4.3173E+00 4.3173E+00 8.6346E+00 5.0000E-01 5.0000E-01 6.7444E-01 0.0000E+00 0.0000E+00 0.0000E+00
57.500 5.5097E+03 0.0000E+00 0.0000E+00 2.4316E+00 2.4316E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 3.8327E+00 3.8327E+00 7.6654E+00 5.0000E-01 5.0000E-01 6.4820E-01 0.0000E+00 0.0000E+00 0.0000E+00
57.501 5.1439E+03 0.0000E+00 0.0000E+00 2.4323E+00 2.4323E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 3.5803E+00 3.5803E+00 7.1605E+00 5.0000E-01 5.0000E-01 6.0516E-01 0.0000E+00 0.0000E+00 0.0000E+00
69.000 4.9357E+03 0.0000E+00 0.0000E+00 2.3339E+00 2.3339E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 3.1629E+00 3.1629E+00 6.3258E+00 5.0000E-01 5.0000E-01 5.8067E-01 0.0000E+00 0.0000E+00 0.0000E+00
69.001 4.5845E+03 0.0000E+00 0.0000E+00 2.3346E+00 2.3346E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 2.9396E+00 2.9396E+00 5.8792E+00 5.0000E-01 5.0000E-01 5.3935E-01 0.0000E+00 0.0000E+00 0.0000E+00
80.500 4.3912E+03 0.0000E+00 0.0000E+00 2.2361E+00 2.2361E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 2.5832E+00 2.5832E+00 5.1664E+00 5.0000E-01 5.0000E-01 5.1661E-01 0.0000E+00 0.0000E+00 0.0000E+00
80.501 4.0547E+03 0.0000E+00 0.0000E+00 2.2368E+00 2.2368E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 2.3867E+00 2.3867E+00 4.7734E+00 5.0000E-01 5.0000E-01 4.7702E-01 0.0000E+00 0.0000E+00 0.0000E+00
92.000 3.8762E+03 0.0000E+00 0.0000E+00 2.1384E+00 2.1384E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 2.0852E+00 2.0852E+00 4.1705E+00 5.0000E-01 5.0000E-01 4.5602E-01 0.0000E+00 0.0000E+00 0.0000E+00
92.001 3.5543E+03 0.0000E+00 0.0000E+00 2.1391E+00 2.1391E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 1.9133E+00 1.9133E+00 3.8267E+00 5.0000E-01 5.0000E-01 4.1816E-01 0.0000E+00 0.0000E+00 0.0000E+00
103.500 3.3908E+03 0.0000E+00 0.0000E+00 2.0406E+00 2.0406E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 1.6611E+00 1.6611E+00 3.3223E+00 5.0000E-01 5.0000E-01 3.9891E-01 0.0000E+00 0.0000E+00 0.0000E+00
103.501 3.0836E+03 0.0000E+00 0.0000E+00 2.0413E+00 2.0413E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 1.5117E+00 1.5117E+00 3.0234E+00 5.0000E-01 5.0000E-01 3.6277E-01 0.0000E+00 0.0000E+00 0.0000E+00
115.630 2.9267E+03 0.0000E+00 0.0000E+00 1.9375E+00 1.9375E+00 0.0000E+00 0.0000E+00 2.1000E+17 8.0769E+17 1.2925E+00 1.2925E+00 2.5850E+00 5.0000E-01 5.0000E-01 3.4432E-01 0.0000E+00 0.0000E+00 0.0000E+00
\ No newline at end of file
wetb/hawc2/tests/test_files/simulation_setup/DTU10MWRef6.0_IOS/input/data/DTU_10MW_RWT_Towertop_st.dat 0 → 100644
View file @ 62f0f2e4
#1 Tower top
r m x_cg y_cg ri_x ri_y x_sh y_sh E G I_x I_y I_p k_x k_y A pitch x_e y_e
$1 2 towertop with nacelle mass on bottom element. flexible
0.0 1.00E-04 0 0 1.36E+00 1.36E+00 0 0 2.10E+11 8.08E+10 5.52E-01 5.52E-01 1.10E+00 0.5 0.5 2.98E-01 0 0 0
2.75 1.00E-04 0 0 1.36E+00 1.36E+00 0 0 2.10E+11 8.08E+10 5.52E-01 5.52E-01 1.10E+00 0.5 0.5 2.98E-01 0 0 0
r m x_cg y_cg ri_x ri_y x_sh y_sh E G I_x I_y I_p k_x k_y A pitch x_e y_e
$2 2 towertop with nacelle mass on bottom element. Stiff
0.0 1.00E-04 0 0 1.36E+00 1.36E+00 0 0 2.10E+17 8.08E+17 5.52E-01 5.52E-01 1.10E+00 0.5 0.5 2.98E-01 0 0 0
2.75 1.00E-04 0 0 1.36E+00 1.36E+00 0 0 2.10E+17 8.08E+17 5.52E-01 5.52E-01 1.10E+00 0.5 0.5 2.98E-01 0 0 0
\ No newline at end of file
wetb/hawc2/tests/test_files/simulation_setup/DTU10MWRef6.0_IOS/input/data/DTU_10MW_RWT_ae.dat 0 → 100644
View file @ 62f0f2e4
1
1 40
0.0000E+00 5.3800E+00 1.0000E+02 1 ;
2.0000E+00 5.3800E+00 1.0000E+02 1 ;
4.7120E+00 5.3800E+00 9.6915E+01 1 ;
5.3960E+00 5.3800E+00 9.5291E+01 1 ;
6.3190E+00 5.3886E+00 9.2664E+01 1 ;
7.4750E+00 5.4212E+00 8.8775E+01 1 ;
8.8580E+00 5.4865E+00 8.3448E+01 1 ;
1.0458E+01 5.5887E+00 7.6689E+01 1 ;
1.2268E+01 5.7247E+00 6.8754E+01 1 ;
1.4275E+01 5.8817E+00 6.0273E+01 1 ;
1.6467E+01 6.0346E+00 5.2291E+01 1 ;
1.8833E+01 6.1478E+00 4.5826E+01 1 ;
2.1356E+01 6.2020E+00 4.0950E+01 1 ;
2.4023E+01 6.1950E+00 3.7343E+01 1 ;
2.6817E+01 6.1292E+00 3.4518E+01 1 ;
2.9721E+01 6.0096E+00 3.2270E+01 1 ;
3.2719E+01 5.8432E+00 3.0488E+01 1 ;
3.5791E+01 5.6400E+00 2.9020E+01 1 ;
3.8920E+01 5.4107E+00 2.7756E+01 1 ;
4.2086E+01 5.1613E+00 2.6693E+01 1 ;
4.5272E+01 4.8974E+00 2.5829E+01 1 ;
4.8457E+01 4.6255E+00 2.5157E+01 1 ;
5.1623E+01 4.3519E+00 2.4665E+01 1 ;
5.4750E+01 4.0827E+00 2.4338E+01 1 ;
5.7820E+01 3.8220E+00 2.4156E+01 1 ;
6.0815E+01 3.5724E+00 2.4100E+01 1 ;
6.3716E+01 3.3364E+00 2.4100E+01 1 ;
6.6506E+01 3.1161E+00 2.4100E+01 1 ;
6.9168E+01 2.9130E+00 2.4100E+01 1 ;
7.1687E+01 2.7275E+00 2.4100E+01 1 ;
7.4047E+01 2.5595E+00 2.4100E+01 1 ;
7.6234E+01 2.4087E+00 2.4100E+01 1 ;
7.8234E+01 2.2660E+00 2.4100E+01 1 ;
8.0037E+01 2.1175E+00 2.4100E+01 1 ;
8.1631E+01 1.9588E+00 2.4100E+01 1 ;
8.3006E+01 1.7913E+00 2.4100E+01 1 ;
8.4155E+01 1.6013E+00 2.4100E+01 1 ;
8.5070E+01 1.3858E+00 2.4100E+01 1 ;
8.5746E+01 1.1384E+00 2.4100E+01 1 ;
8.6366E+01 8.3354E-01 2.4100E+01 1 ;
wetb/hawc2/tests/test_files/simulation_setup/DTU10MWRef6.0_IOS/input/data/DTU_10MW_RWT_pc.dat 0 → 100644
View file @ 62f0f2e4
1 DTU 10MW RWT, DTU Wind Energy. Data generated using EllipSys2D with free transition TI=0.1 and 3D corrected using the model by Bak et al.
6
1 105 24.10000 FFA-W3-241 (Re=12x10^6)
-180.00000 0.00000 0.00000 0.00000
-175.00000 0.17360 0.01140 0.02180
-170.00000 0.34200 0.04520 0.04340
-165.00000 0.50000 0.10050 0.06470
-160.00000 0.64280 0.17550 0.08550
-155.00000 0.76600 0.26790 0.10570
-150.00000 0.86600 0.37500 0.12500
-145.00000 0.93970 0.49350 0.14340
-140.00000 0.98480 0.61970 0.16070
-135.00000 1.00000 0.75000 0.17680
-130.00000 0.98480 0.88030 0.19150
-125.00000 0.93970 1.00650 0.20480
-120.00000 0.86600 1.12500 0.21650
-115.00000 0.76600 1.23210 0.22660
-110.00000 0.64280 1.32450 0.23490
-105.00000 0.50000 1.39950 0.24150
-100.00000 0.34200 1.45480 0.24620
-95.00000 0.17360 1.48860 0.24900
-90.00000 0.00000 1.50000 0.25000
-85.00000 -0.17360 1.48860 0.24900
-80.00000 -0.34200 1.45480 0.24620
-75.00000 -0.50000 1.39950 0.24150
-70.00000 -0.64280 1.32450 0.23490
-65.00000 -0.76600 1.23210 0.22660
-60.00000 -0.86600 1.12500 0.21650
-55.00000 -0.93970 1.00650 0.20480
-50.00000 -0.98480 0.86030 0.19150
-45.00000 -1.01200 0.71200 0.17080
-40.00000 -1.03760 0.54750 0.14160
-39.00000 -1.04190 0.51650 0.13460
-38.00000 -1.04620 0.48160 0.12760
-37.00000 -1.05250 0.44870 0.11850
-36.00000 -1.05680 0.41770 0.10950
-35.00000 -1.06110 0.38480 0.10050
-34.00000 -1.06540 0.35580 0.08940
-33.00000 -1.07170 0.32890 0.08240
-32.00000 -1.07650 0.30310 0.06780
-30.00000 -1.08890 0.25600 0.05080
-28.00000 -1.09930 0.20900 0.03370
-26.00000 -1.11680 0.17560 0.02060
-24.00000 -1.12820 0.14230 0.00750
-22.00000 -1.12150 0.11830 -0.00080
-20.00000 -1.11480 0.09430 -0.00910
-18.00000 -1.09190 0.07650 -0.01230
-16.00000 -1.06910 0.05870 -0.01560
-14.00000 -1.03790 0.04540 -0.01560
-12.00000 -1.00670 0.03210 -0.01550
-10.00000 -0.84790 0.02300 -0.03180
-8.00000 -0.68920 0.01380 -0.04800
-6.00000 -0.42780 0.01180 -0.06110
-4.00000 -0.16650 0.00980 -0.07420
-2.00000 0.08630 0.00950 -0.08110
0.00000 0.33910 0.00920 -0.08800
2.00000 0.58670 0.00940 -0.09330
4.00000 0.83010 0.00990 -0.09770
6.00000 1.06560 0.01090 -0.10080
8.00000 1.29140 0.01240 -0.10260
10.00000 1.50120 0.01440 -0.10240
12.00000 1.68860 0.01730 -0.09980
14.00000 1.81030 0.02260 -0.09410
16.00000 1.81390 0.03540 -0.08740
18.00000 1.75450 0.06470 -0.08500
20.00000 1.60710 0.10350 -0.09130
22.00000 1.52570 0.14370 -0.10260
24.00000 1.44280 0.18410 -0.11400
26.00000 1.38260 0.22900 -0.12820
28.00000 1.32180 0.27380 -0.14230
30.00000 1.25830 0.32780 -0.15820
32.00000 1.19440 0.38160 -0.17390
33.00000 1.17340 0.42090 -0.17900
34.00000 1.15680 0.45040 -0.18310
35.00000 1.13790 0.48150 -0.18730
36.00000 1.12480 0.51050 -0.19140
37.00000 1.11770 0.54500 -0.19530
38.00000 1.10240 0.57370 -0.19930
39.00000 1.08910 0.59990 -0.20320
40.00000 1.07550 0.62800 -0.20710
45.00000 1.01750 0.75780 -0.21610
50.00000 0.97160 0.88200 -0.21840
55.00000 0.92680 1.01040 -0.22140
60.00000 0.86600 1.12500 -0.22550
65.00000 0.76600 1.23210 -0.23060
70.00000 0.64280 1.32450 -0.23490
75.00000 0.50000 1.39950 -0.24150
80.00000 0.34200 1.45480 -0.24620
85.00000 0.17360 1.48860 -0.24900
90.00000 0.00000 1.50000 -0.25000
95.00000 -0.17360 1.48860 -0.24900
100.00000 -0.34200 1.45480 -0.24620
105.00000 -0.50000 1.39950 -0.24150
110.00000 -0.64280 1.32450 -0.23490
115.00000 -0.76600 1.23210 -0.22660
120.00000 -0.86600 1.12500 -0.21650
125.00000 -0.93970 1.00650 -0.20480
130.00000 -0.98480 0.88030 -0.19150
135.00000 -1.00000 0.75000 -0.17680
140.00000 -0.98480 0.61970 -0.16070
145.00000 -0.93970 0.49350 -0.14340
150.00000 -0.86600 0.37500 -0.12500
155.00000 -0.76600 0.26790 -0.10570
160.00000 -0.64280 0.17550 -0.08550
165.00000 -0.50000 0.10050 -0.06470
170.00000 -0.34200 0.04520 -0.04340
175.00000 -0.17360 0.01140 -0.02180
180.00000 0.00000 0.00000 0.00000
2 105 30.10000 FFA-W3-301 (Re=10x10^6)
-180.00000 0.00000 0.00000 0.00000
-175.00000 0.17360 0.00990 0.02180
-170.00000 0.34200 0.03920 0.04340
-165.00000 0.50000 0.08710 0.06470
-160.00000 0.64280 0.15210 0.08550
-155.00000 0.76600 0.23220 0.10570
-150.00000 0.86600 0.32500 0.12500
-145.00000 0.93970 0.42770 0.14340
-140.00000 0.98480 0.53710 0.16070
-135.00000 1.00000 0.65000 0.17680
-130.00000 0.98480 0.76290 0.19150
-125.00000 0.93970 0.87230 0.20480
-120.00000 0.86600 0.97500 0.21650
-115.00000 0.76600 1.06780 0.22660
-110.00000 0.64280 1.14790 0.23490
-105.00000 0.50000 1.21290 0.24150
-100.00000 0.34200 1.26080 0.24620
-95.00000 0.17360 1.29010 0.24900
-90.00000 0.00000 1.30000 0.25000
-85.00000 -0.17360 1.29010 0.24900
-80.00000 -0.34200 1.26080 0.24620
-75.00000 -0.50000 1.21290 0.24150
-70.00000 -0.64280 1.14790 0.23490
-65.00000 -0.76600 1.06780 0.22660
-60.00000 -0.86600 0.97500 0.21650
-55.00000 -0.93970 0.87230 0.20480
-50.00000 -0.98480 0.76290 0.19150
-45.00000 -1.00000 0.65000 0.17680
-40.00000 -1.02160 0.53520 0.14160
-39.00000 -1.02390 0.50660 0.13460
-38.00000 -1.01620 0.48510 0.12760
-37.00000 -1.01850 0.45650 0.12050
-36.00000 -1.01080 0.42790 0.11350
-35.00000 -1.00310 0.39940 0.09650
-34.00000 -0.99540 0.37080 0.07940
-33.00000 -0.98770 0.33530 0.06240
-32.00000 -0.98350 0.30970 0.05150
-30.00000 -0.97030 0.26630 0.03890
-28.00000 -0.96720 0.22290 0.02630
-26.00000 -0.94410 0.19410 0.01540
-24.00000 -0.93100 0.16540 0.00450
-22.00000 -0.90470 0.14170 -0.00390
-20.00000 -0.87840 0.11810 -0.01230
-18.00000 -0.84590 0.09860 -0.01750
-16.00000 -0.81340 0.07920 -0.02270
-14.00000 -0.77280 0.06430 -0.02350
-12.00000 -0.73220 0.04950 -0.02440
-10.00000 -0.69350 0.03810 -0.02270
-8.00000 -0.65470 0.02670 -0.02100
-6.00000 -0.45070 0.02040 -0.03890
-4.00000 -0.24670 0.01400 -0.05690
-2.00000 0.02950 0.01290 -0.07170
0.00000 0.30560 0.01180 -0.08650
2.00000 0.56700 0.01190 -0.09540
4.00000 0.81990 0.01250 -0.10240
6.00000 1.06140 0.01360 -0.10710
8.00000 1.28740 0.01520 -0.10940
10.00000 1.48400 0.01800 -0.10796
12.00000 1.63880 0.02240 -0.10256
14.00000 1.73270 0.03030 -0.09524
16.00000 1.71420 0.05390 -0.09104
18.00000 1.68280 0.09540 -0.09794
20.00000 1.65670 0.14350 -0.11030
22.00000 1.64440 0.22800 -0.12996
24.00000 1.63290 0.31480 -0.15332
26.00000 1.63330 0.39260 -0.17600
28.00000 1.61750 0.46230 -0.19658
30.00000 1.59750 0.52660 -0.21490
32.00000 1.57080 0.58550 -0.23130
33.00000 1.57580 0.62150 -0.23660
34.00000 1.55710 0.65620 -0.24060
35.00000 1.54750 0.68970 -0.24470
36.00000 1.52680 0.71480 -0.24870
37.00000 1.52590 0.74640 -0.24960
38.00000 1.51390 0.77620 -0.24550
39.00000 1.49180 0.79820 -0.24440
40.00000 1.45960 0.81970 -0.24230
45.00000 1.29510 0.91220 -0.23730
50.00000 1.12560 0.96650 -0.23537
55.00000 0.98080 0.99570 -0.23243
60.00000 0.86600 1.03500 -0.23350
65.00000 0.76600 1.08280 -0.23360
70.00000 0.64280 1.14790 -0.23490
75.00000 0.50000 1.21290 -0.24150
80.00000 0.34200 1.26080 -0.24620
85.00000 0.17360 1.29010 -0.24900
90.00000 0.00000 1.30000 -0.25000
95.00000 -0.17360 1.29010 -0.24900
100.00000 -0.34200 1.26080 -0.24620
105.00000 -0.50000 1.21290 -0.24150
110.00000 -0.64280 1.14790 -0.23490
115.00000 -0.76600 1.06780 -0.22660
120.00000 -0.86600 0.97500 -0.21650
125.00000 -0.93970 0.87230 -0.20480
130.00000 -0.98480 0.76290 -0.19150
135.00000 -1.00000 0.65000 -0.17680
140.00000 -0.98480 0.53710 -0.16070
145.00000 -0.93970 0.42770 -0.14340
150.00000 -0.86600 0.32500 -0.12500
155.00000 -0.76600 0.23220 -0.10570
160.00000 -0.64280 0.15210 -0.08550
165.00000 -0.50000 0.08710 -0.06470
170.00000 -0.34200 0.03920 -0.04340
175.00000 -0.17360 0.00990 -0.02180
180.00000 0.00000 0.00000 0.00000
3 105 36.00000 FFA-W3-360 (Re=10x10^6)
-180.00000 0.00000 0.00000 0.00000
-175.00000 0.17360 0.00990 0.02180
-170.00000 0.34200 0.03920 0.04340
-165.00000 0.50000 0.08710 0.06470
-160.00000 0.64280 0.15210 0.08550
-155.00000 0.76600 0.23220 0.10570
-150.00000 0.86600 0.32500 0.12500
-145.00000 0.93970 0.42770 0.14340
-140.00000 0.98480 0.53710 0.16070
-135.00000 1.00000 0.65000 0.17680
-130.00000 0.98480 0.76290 0.19150
-125.00000 0.93970 0.87230 0.20480
-120.00000 0.86600 0.97500 0.21650
-115.00000 0.76600 1.06780 0.22660
-110.00000 0.64280 1.14790 0.23490
-105.00000 0.50000 1.21290 0.24150
-100.00000 0.34200 1.26080 0.24620
-95.00000 0.17360 1.29010 0.24900
-90.00000 0.00000 1.30000 0.25000
-85.00000 -0.17360 1.29010 0.24900
-80.00000 -0.34200 1.26080 0.24620
-75.00000 -0.50000 1.21290 0.24150
-70.00000 -0.64280 1.14790 0.23490
-65.00000 -0.76600 1.06780 0.22660
-60.00000 -0.86600 0.97500 0.21650
-55.00000 -0.93970 0.87230 0.20480
-50.00000 -0.98480 0.76290 0.19150
-45.00000 -1.00000 0.65000 0.17680
-40.00000 -0.97160 0.53520 0.14160
-39.00000 -0.96390 0.51360 0.13460
-38.00000 -0.94620 0.48510 0.12760
-37.00000 -0.92850 0.46350 0.12050
-36.00000 -0.91080 0.43490 0.11350
-35.00000 -0.89310 0.40640 0.09650
-34.00000 -0.88540 0.37780 0.07940
-33.00000 -0.86770 0.35230 0.06240
-32.00000 -0.85640 0.33150 0.04470
-30.00000 -0.83780 0.27770 0.03270
-28.00000 -0.81910 0.23380 0.02080
-26.00000 -0.79240 0.20450 0.01020
-24.00000 -0.77560 0.17510 -0.00040
-22.00000 -0.74410 0.15130 -0.00890
-20.00000 -0.71260 0.12740 -0.01740
-18.00000 -0.66780 0.10850 -0.02290
-16.00000 -0.62310 0.08960 -0.02850
-14.00000 -0.57420 0.07480 -0.03000
-12.00000 -0.52520 0.06000 -0.03140
-10.00000 -0.48270 0.04850 -0.02790
-8.00000 -0.44020 0.03700 -0.02440
-6.00000 -0.29830 0.02940 -0.03590
-4.00000 -0.15640 0.02190 -0.04740
-2.00000 0.17440 0.02030 -0.07820
0.00000 0.50530 0.01870 -0.10900
2.00000 0.82410 0.01880 -0.13290
4.00000 1.12090 0.01960 -0.15100
6.00000 1.38970 0.02130 -0.16298
8.00000 1.62540 0.02400 -0.16908
10.00000 1.81090 0.02790 -0.16832
12.00000 1.85890 0.03650 -0.15848
14.00000 1.81590 0.07600 -0.15576
16.00000 1.77860 0.11650 -0.16696
18.00000 1.75600 0.15710 -0.18094
20.00000 1.76300 0.20630 -0.19598
22.00000 1.80020 0.30570 -0.21428
24.00000 1.84950 0.41530 -0.23648
26.00000 1.87750 0.51630 -0.25882
28.00000 1.88280 0.60690 -0.27966
30.00000 1.86890 0.68920 -0.29830
32.00000 1.84390 0.76250 -0.31466
33.00000 1.83490 0.79540 -0.31660
34.00000 1.82970 0.82320 -0.31060
35.00000 1.81360 0.84460 -0.30470
36.00000 1.78540 0.86840 -0.29870
37.00000 1.76690 0.89120 -0.29260
38.00000 1.72690 0.90660 -0.28650
39.00000 1.68640 0.92550 -0.28440
40.00000 1.64580 0.93850 -0.28030
45.00000 1.44730 1.02900 -0.26330
50.00000 1.13560 1.07150 -0.25437
55.00000 0.98080 1.08570 -0.24543
60.00000 0.86600 1.10000 -0.23650
65.00000 0.76600 1.12780 -0.23660
70.00000 0.64280 1.16290 -0.23890
75.00000 0.50000 1.21290 -0.24150
80.00000 0.34200 1.26080 -0.24620
85.00000 0.17360 1.29010 -0.24900
90.00000 0.00000 1.30000 -0.25000
95.00000 -0.17360 1.29010 -0.24900
100.00000 -0.34200 1.26080 -0.24620
105.00000 -0.50000 1.21290 -0.24150
110.00000 -0.64280 1.14790 -0.23490
115.00000 -0.76600 1.06780 -0.22660
120.00000 -0.86600 0.97500 -0.21650
125.00000 -0.93970 0.87230 -0.20480
130.00000 -0.98480 0.76290 -0.19150
135.00000 -1.00000 0.65000 -0.17680
140.00000 -0.98480 0.53710 -0.16070
145.00000 -0.93970 0.42770 -0.14340
150.00000 -0.86600 0.32500 -0.12500
155.00000 -0.76600 0.23220 -0.10570
160.00000 -0.64280 0.15210 -0.08550
165.00000 -0.50000 0.08710 -0.06470
170.00000 -0.34200 0.03920 -0.04340
175.00000 -0.17360 0.00990 -0.02180
180.00000 0.00000 0.00000 0.00000
4 105 48.00000 FFA-W3-"480" (Re=10x10^6)
-180.00000 0.00000 0.00000 0.00000
-175.00000 0.17360 0.00990 0.02180
-170.00000 0.34200 0.03920 0.04340
-165.00000 0.50000 0.08710 0.06470
-160.00000 0.64280 0.15210 0.08550
-155.00000 0.76600 0.23220 0.10570
-150.00000 0.86600 0.32500 0.12500
-145.00000 0.93970 0.42770 0.14340
-140.00000 0.98480 0.53710 0.16070
-135.00000 1.00000 0.65000 0.17680
-130.00000 0.98480 0.76290 0.19150
-125.00000 0.93970 0.87230 0.20480
-120.00000 0.86600 0.97500 0.21650
-115.00000 0.76600 1.06780 0.22660
-110.00000 0.64280 1.14790 0.23490
-105.00000 0.50000 1.21290 0.24150
-100.00000 0.34200 1.26080 0.24620
-95.00000 0.17360 1.29010 0.24900
-90.00000 0.00000 1.30000 0.25000
-85.00000 -0.17360 1.29010 0.24900
-80.00000 -0.34200 1.26080 0.24620
-75.00000 -0.50000 1.21290 0.24150
-70.00000 -0.64280 1.14790 0.23490
-65.00000 -0.76600 1.06780 0.22660
-60.00000 -0.86600 0.97500 0.21650
-55.00000 -0.93970 0.87230 0.19780
-50.00000 -0.98480 0.76290 0.17750
-45.00000 -1.00000 0.65000 0.15580
-40.00000 -0.98160 0.53520 0.12460
-39.00000 -0.95390 0.51360 0.11060
-38.00000 -0.92620 0.48510 0.09660
-37.00000 -0.88850 0.45650 0.08250
-36.00000 -0.85080 0.42790 0.06550
-35.00000 -0.82310 0.39240 0.05850
-34.00000 -0.78540 0.36380 0.04840
-33.00000 -0.74770 0.33830 0.03840
-32.00000 -0.70110 0.31230 0.03050
-30.00000 -0.62080 0.26630 0.02000
-28.00000 -0.54060 0.24020 0.00940
-26.00000 -0.46940 0.21690 -0.00060
-24.00000 -0.38810 0.19360 -0.01060
-22.00000 -0.31610 0.17320 -0.01970
-20.00000 -0.24420 0.15290 -0.02890
-18.00000 -0.16410 0.13550 -0.03680
-16.00000 -0.08410 0.11800 -0.04480
-14.00000 0.00210 0.10350 -0.05100
-12.00000 0.08830 0.08900 -0.05720
-10.00000 0.17220 0.07730 -0.06040
-8.00000 0.25610 0.06560 -0.06360
-6.00000 0.31790 0.05630 -0.06050
-4.00000 0.37980 0.04700 -0.05740
-2.00000 0.35730 0.04050 -0.03510
0.00000 0.33480 0.03410 -0.01280
2.00000 0.56520 0.03160 -0.04940
4.00000 0.87690 0.03430 -0.08940
6.00000 1.04250 0.04510 -0.11168
8.00000 0.94870 0.07000 -0.12084
10.00000 0.90880 0.08860 -0.13760
12.00000 0.97610 0.09930 -0.15942
14.00000 1.11300 0.10700 -0.18234
16.00000 1.30650 0.11630 -0.20664
18.00000 1.54140 0.13170 -0.23150
20.00000 1.80490 0.15700 -0.25662
22.00000 2.00200 0.27570 -0.28058
24.00000 2.12160 0.42240 -0.30574
26.00000 2.19160 0.56100 -0.32940
28.00000 2.22910 0.68610 -0.35198
30.00000 2.23220 0.79580 -0.37100
32.00000 2.21880 0.89150 -0.38784
33.00000 2.20920 0.92290 -0.38450
34.00000 2.19270 0.94720 -0.38060
35.00000 2.16250 0.96640 -0.37470
36.00000 2.13130 0.98890 -0.36870
37.00000 2.08780 1.00730 -0.36260
38.00000 2.04350 1.02840 -0.35650
39.00000 1.99840 1.05210 -0.35040
40.00000 1.94180 1.07690 -0.34430
45.00000 1.65350 1.16970 -0.32330
50.00000 1.30160 1.19650 -0.30437
55.00000 1.03080 1.19570 -0.29543
60.00000 0.86600 1.19000 -0.28650
65.00000 0.76600 1.19280 -0.27660
70.00000 0.64280 1.20290 -0.26490
75.00000 0.50000 1.21290 -0.25150
80.00000 0.34200 1.26080 -0.24620
85.00000 0.17360 1.29010 -0.24900
90.00000 0.00000 1.30000 -0.25000
95.00000 -0.17360 1.29010 -0.24900
100.00000 -0.34200 1.26080 -0.24620
105.00000 -0.50000 1.21290 -0.24150
110.00000 -0.64280 1.14790 -0.23490
115.00000 -0.76600 1.06780 -0.22660
120.00000 -0.86600 0.97500 -0.21650
125.00000 -0.93970 0.87230 -0.20480
130.00000 -0.98480 0.76290 -0.19150
135.00000 -1.00000 0.65000 -0.17680
140.00000 -0.98480 0.53710 -0.16070
145.00000 -0.93970 0.42770 -0.14340
150.00000 -0.86600 0.32500 -0.12500
155.00000 -0.76600 0.23220 -0.10570
160.00000 -0.64280 0.15210 -0.08550
165.00000 -0.50000 0.08710 -0.06470
170.00000 -0.34200 0.03920 -0.04340
175.00000 -0.17360 0.00990 -0.02180
180.00000 0.00000 0.00000 0.00000
5 105 60.00000 FFA-W3-"600" (Re=6x10^6)
-180.00000 0.00000 0.00000 0.00000
-175.00000 0.17360 0.00990 0.02180
-170.00000 0.34200 0.03920 0.04340
-165.00000 0.50000 0.08710 0.06470
-160.00000 0.64280 0.15210 0.08550
-155.00000 0.76600 0.23220 0.10570
-150.00000 0.86600 0.32500 0.12500
-145.00000 0.93970 0.42770 0.14340
-140.00000 0.98480 0.53710 0.16070
-135.00000 1.00000 0.65000 0.17680
-130.00000 0.98480 0.76290 0.19150
-125.00000 0.93970 0.87230 0.20480
-120.00000 0.86600 0.97500 0.21650
-115.00000 0.76600 1.06780 0.22660
-110.00000 0.64280 1.14790 0.23490
-105.00000 0.50000 1.21290 0.24150
-100.00000 0.34200 1.26080 0.24620
-95.00000 0.17360 1.29010 0.24900
-90.00000 0.00000 1.30000 0.25000
-85.00000 -0.17360 1.29010 0.24900
-80.00000 -0.34200 1.26080 0.24620
-75.00000 -0.50000 1.21290 0.24150
-70.00000 -0.64280 1.14790 0.23490
-65.00000 -0.76600 1.06780 0.22660
-60.00000 -0.86600 0.97500 0.21650
-55.00000 -0.93970 0.87230 0.19780
-50.00000 -0.98480 0.76290 0.17750
-45.00000 -1.00000 0.65000 0.14880
-40.00000 -0.95160 0.52820 0.11760
-39.00000 -0.90390 0.50260 0.10760
-38.00000 -0.82620 0.46710 0.09060
-37.00000 -0.71850 0.43450 0.08050
-36.00000 -0.62080 0.39890 0.06650
-35.00000 -0.52310 0.36340 0.05950
-34.00000 -0.44540 0.32780 0.05240
-33.00000 -0.34770 0.30130 0.04540
-32.00000 -0.25710 0.27290 0.04699
-30.00000 -0.15740 0.23320 0.06060
-28.00000 -0.07760 0.20340 0.07420
-26.00000 0.01150 0.18580 0.08749
-24.00000 0.10070 0.16810 0.10079
-22.00000 0.19620 0.15260 0.11329
-20.00000 0.29180 0.13700 0.12579
-18.00000 0.39320 0.12400 0.13700
-16.00000 0.49470 0.11090 0.14822
-14.00000 0.60040 0.10090 0.15711
-12.00000 0.70610 0.09080 0.16600
-10.00000 0.78800 0.08360 0.16771
-8.00000 0.87000 0.07650 0.16942
-6.00000 0.83580 0.07480 -0.14700
-4.00000 0.81620 0.07311 -0.12588
-2.00000 0.66030 0.07553 -0.08764
0.00000 0.51990 0.07795 -0.05062
2.00000 0.26360 0.07795 -0.00078
4.00000 -0.09060 0.07082 0.06962
6.00000 0.16560 0.07483 0.01974
8.00000 0.44750 0.08587 -0.02718
10.00000 0.70710 0.09990 -0.06712
12.00000 0.95400 0.11148 -0.10334
14.00000 1.18910 0.12072 -0.13694
16.00000 1.41830 0.13678 -0.16866
18.00000 1.63920 0.14998 -0.19886
20.00000 1.82160 0.16894 -0.22524
22.00000 1.98370 0.20657 -0.25094
24.00000 2.09850 0.30919 -0.27474
26.00000 2.19040 0.43440 -0.29658
28.00000 2.25410 0.55120 -0.31712
30.00000 2.27310 0.67320 -0.33352
32.00000 2.27400 0.77560 -0.34854
33.00000 2.27800 0.82090 -0.34660
34.00000 2.26580 0.85510 -0.34760
35.00000 2.24080 0.88620 -0.34870
36.00000 2.23310 0.91410 -0.34270
37.00000 2.21340 0.93900 -0.33660
38.00000 2.18200 0.97090 -0.33050
39.00000 2.13910 0.99960 -0.32440
40.00000 2.09510 1.02540 -0.31830
45.00000 1.83700 1.14190 -0.30030
50.00000 1.51560 1.16150 -0.29137
55.00000 1.20080 1.16570 -0.28243
60.00000 0.96600 1.17500 -0.27350
65.00000 0.76600 1.17780 -0.26360
70.00000 0.64280 1.19790 -0.25490
75.00000 0.50000 1.21290 -0.25150
80.00000 0.34200 1.26080 -0.24620
85.00000 0.17360 1.29010 -0.24900
90.00000 0.00000 1.30000 -0.25000
95.00000 -0.17360 1.29010 -0.24900
100.00000 -0.34200 1.26080 -0.24620
105.00000 -0.50000 1.21290 -0.24150
110.00000 -0.64280 1.14790 -0.23490
115.00000 -0.76600 1.06780 -0.22660
120.00000 -0.86600 0.97500 -0.21650
125.00000 -0.93970 0.87230 -0.20480
130.00000 -0.98480 0.76290 -0.19150
135.00000 -1.00000 0.65000 -0.17680
140.00000 -0.98480 0.53710 -0.16070
145.00000 -0.93970 0.42770 -0.14340
150.00000 -0.86600 0.32500 -0.12500
155.00000 -0.76600 0.23220 -0.10570
160.00000 -0.64280 0.15210 -0.08550
165.00000 -0.50000 0.08710 -0.06470
170.00000 -0.34200 0.03920 -0.04340
175.00000 -0.17360 0.00990 -0.02180
180.00000 0.00000 0.00000 0.00000
6 105 100.00000 Cylinder
-180.00000 0.00000 0.60000 0.00000
-175.00000 0.00000 0.60000 0.00000
-170.00000 0.00000 0.60000 0.00000
-165.00000 0.00000 0.60000 0.00000
-160.00000 0.00000 0.60000 0.00000
-155.00000 0.00000 0.60000 0.00000
-150.00000 0.00000 0.60000 0.00000
-145.00000 0.00000 0.60000 0.00000
-140.00000 0.00000 0.60000 0.00000
-135.00000 0.00000 0.60000 0.00000
-130.00000 0.00000 0.60000 0.00000
-125.00000 0.00000 0.60000 0.00000
-120.00000 0.00000 0.60000 0.00000
-115.00000 0.00000 0.60000 0.00000
-110.00000 0.00000 0.60000 0.00000
-105.00000 0.00000 0.60000 0.00000
-100.00000 0.00000 0.60000 0.00000
-95.00000 0.00000 0.60000 0.00000
-90.00000 0.00000 0.60000 0.00000
-85.00000 0.00000 0.60000 0.00000
-80.00000 0.00000 0.60000 0.00000
-75.00000 0.00000 0.60000 0.00000
-70.00000 0.00000 0.60000 0.00000
-65.00000 0.00000 0.60000 0.00000
-60.00000 0.00000 0.60000 0.00000
-55.00000 0.00000 0.60000 0.00000
-50.00000 0.00000 0.60000 0.00000
-45.00000 0.00000 0.60000 0.00000
-40.00000 0.00000 0.60000 0.00000
-39.00000 0.00000 0.60000 0.00000
-38.00000 0.00000 0.60000 0.00000
-37.00000 0.00000 0.60000 0.00000
-36.00000 0.00000 0.60000 0.00000
-35.00000 0.00000 0.60000 0.00000
-34.00000 0.00000 0.60000 0.00000
-33.00000 0.00000 0.60000 0.00000
-32.00000 0.00000 0.60000 0.00000
-30.00000 0.00000 0.60000 0.00000
-28.00000 0.00000 0.60000 0.00000
-26.00000 0.00000 0.60000 0.00000
-24.00000 0.00000 0.60000 0.00000
-22.00000 0.00000 0.60000 0.00000
-20.00000 0.00000 0.60000 0.00000
-18.00000 0.00000 0.60000 0.00000
-16.00000 0.00000 0.60000 0.00000
-14.00000 0.00000 0.60000 0.00000
-12.00000 0.00000 0.60000 0.00000
-10.00000 0.00000 0.60000 0.00000
-8.00000 0.00000 0.60000 0.00000
-6.00000 0.00000 0.60000 0.00000
-4.00000 0.00000 0.60000 0.00000
-2.00000 0.00000 0.60000 0.00000
0.00000 0.00000 0.60000 0.00000
2.00000 0.00000 0.60000 0.00000
4.00000 0.00000 0.60000 0.00000
6.00000 0.00000 0.60000 0.00000
8.00000 0.00000 0.60000 0.00000
10.00000 0.00000 0.60000 0.00000
12.00000 0.00000 0.60000 0.00000
14.00000 0.00000 0.60000 0.00000
16.00000 0.00000 0.60000 0.00000
18.00000 0.00000 0.60000 0.00000
20.00000 0.00000 0.60000 0.00000
22.00000 0.00000 0.60000 0.00000
24.00000 0.00000 0.60000 0.00000
26.00000 0.00000 0.60000 0.00000
28.00000 0.00000 0.60000 0.00000
30.00000 0.00000 0.60000 0.00000
32.00000 0.00000 0.60000 0.00000
33.00000 0.00000 0.60000 0.00000
34.00000 0.00000 0.60000 0.00000
35.00000 0.00000 0.60000 0.00000
36.00000 0.00000 0.60000 0.00000
37.00000 0.00000 0.60000 0.00000
38.00000 0.00000 0.60000 0.00000
39.00000 0.00000 0.60000 0.00000
40.00000 0.00000 0.60000 0.00000
45.00000 0.00000 0.60000 0.00000
50.00000 0.00000 0.60000 0.00000
55.00000 0.00000 0.60000 0.00000
60.00000 0.00000 0.60000 0.00000
65.00000 0.00000 0.60000 0.00000
70.00000 0.00000 0.60000 0.00000
75.00000 0.00000 0.60000 0.00000
80.00000 0.00000 0.60000 0.00000
85.00000 0.00000 0.60000 0.00000
90.00000 0.00000 0.60000 0.00000
95.00000 0.00000 0.60000 0.00000
100.00000 0.00000 0.60000 0.00000
105.00000 0.00000 0.60000 0.00000
110.00000 0.00000 0.60000 0.00000
115.00000 0.00000 0.60000 0.00000
120.00000 0.00000 0.60000 0.00000
125.00000 0.00000 0.60000 0.00000
130.00000 0.00000 0.60000 0.00000
135.00000 0.00000 0.60000 0.00000
140.00000 0.00000 0.60000 0.00000
145.00000 0.00000 0.60000 0.00000
150.00000 0.00000 0.60000 0.00000
155.00000 0.00000 0.60000 0.00000
160.00000 0.00000 0.60000 0.00000
165.00000 0.00000 0.60000 0.00000
170.00000 0.00000 0.60000 0.00000
175.00000 0.00000 0.60000 0.00000
180.00000 0.00000 0.60000 0.00000
wetb/hawc2/tests/test_files/simulation_setup/DTU10MWRef6.0_IOS/input/htc/DTU_10MW_RWT.htc 0 → 100644
View file @ 62f0f2e4
;DTU_10MW_RWT, version 5, 04-21-15, anyd
;
begin simulation;
time_stop 1;
solvertype 1; (newmark)
on_no_convergence continue;
convergence_limits 1000 1 1e-07;
logfile ./log/dtu_10mw_rwt_ver4.log;
begin newmark;
deltat 0.02;
end newmark;
end simulation; ;
;----------------------------------------------------------------------------------------------------------------------------------------------------------------
begin new_htc_structure; ; beam_output_file_name ./log/DTU_10MW_RWT_beam.dat; Optional - Calculated beam properties of the bodies are written to file
; body_output_file_name ./log/DTU_10MW_RWT_body.dat; Optional - Body initial position and orientation are written to file
; body_eigenanalysis_file_name ./eig/DTU_10MW_RWT_body_eigen.dat;
; structure_eigenanalysis_file_name ./eig/DTU_10MW_RWT_strc_eigen.dat ;
;-------------------------------------------------------------------------------------------------------------------------------
;-------------------------------------------------------------------------------------------------------------------------------
begin main_body; tower 115m
name tower;
type timoschenko;
nbodies 1;
node_distribution c2_def;
damping_posdef 0 0 0 0.00412 0.00412 0.00045; Mx My Mz Kx Ky Kz , M´s raises overall level, K´s raises high freguency level "tuned by Larh"
begin timoschenko_input;
filename ./data/dtu_10mw_rwt_tower_st.dat;
set 1 1;
end timoschenko_input;
begin c2_def; Definition of centerline (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;
name towertop;
type timoschenko;
nbodies 1;
node_distribution c2_def;
damping_posdef 0 0 0 0.007 0.007 0.007; "changed 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;
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;
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;
end timoschenko_input;
begin c2_def; Definition of centerline (main_body coordinates)
nsec 5;
sec 1 0 0 0 0; Tower top 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;
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;
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.8 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 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; set subset
end timoschenko_input;
begin c2_def; Definition of centerline (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;
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; initial position of node 1
body_eulerang 0 0 0;
end base; ;
begin relative;
body1 tower last;
body2 towertop 1;
body2_eulerang 0 0 0;
end relative; ;
begin relative;
body1 towertop last;
body2 shaft 1;
body2_eulerang 90 0 0;
body2_eulerang 5 0 0; 5 deg tilt angle
body2_eulerang 0 0 0;
mbdy2_ini_rotvec_d1 0 0 -1 0.7; mbdy2_ini_rotvec_d1 0.0 0.0 -1.0 [init_wr];
end relative; ;
begin relative;
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
end relative; ;
begin relative;
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
end relative; ;
begin relative;
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
end relative; ;
begin relative;
body1 hub1 last;
body2 blade1 1;
body2_eulerang 0 0 0;
end relative; ;
begin relative;
body1 hub2 last;
body2 blade2 1;
body2_eulerang 0 0 0;
end relative; ;
begin relative;
body1 hub3 last;
body2 blade3 1;
body2_eulerang 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 -1; 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 -1;
end bearing2; ;
begin bearing2;
name pitch2;
body1 hub2 last;
body2 blade2 1;
bearing_vector 2 0 0 -1;
end bearing2; ;
begin bearing2;
name pitch3;
body1 hub3 last;
body2 blade3 1;
bearing_vector 2 0 0 -1;
end bearing2;
end constraint; ;
end new_htc_structure; ;----------------------------------------------------------------------------------------------------------------------------------------------------------------
begin wind;
density 1.225;
wsp 4;
tint 0;
horizontal_input 1;
windfield_rotations 0 0 0; yaw, tilt, rotation
center_pos0 0 0 -119; hub heigth
shear_format 1 0.2;
turb_format 1; 0=none, 1=mann,2=flex
tower_shadow_method 0; 0=none, 1=potential flow, 2=jet
scale_time_start 0;
wind_ramp_factor 0 40 0.6 1; ; Steps ;
wind_ramp_abs 140 141 0 1; wsp. after the step: 5.0
wind_ramp_abs 181 182 0 1; wsp. after the step: 6.0
wind_ramp_abs 222 223 0 1; wsp. after the step: 7.0
wind_ramp_abs 263 264 0 1; wsp. after the step: 8.0
wind_ramp_abs 304 305 0 1; wsp. after the step: 9.0
wind_ramp_abs 345 346 0 1; wsp. after the step: 10.0
wind_ramp_abs 386 387 0 1; wsp. after the step: 11.0
wind_ramp_abs 427 428 0 1; wsp. after the step: 12.0
wind_ramp_abs 468 469 0 1; wsp. after the step: 13.0
wind_ramp_abs 509 510 0 1; wsp. after the step: 14.0
wind_ramp_abs 550 551 0 1; wsp. after the step: 15.0
wind_ramp_abs 591 592 0 1; wsp. after the step: 16.0
wind_ramp_abs 632 633 0 1; wsp. after the step: 17.0
wind_ramp_abs 673 674 0 1; wsp. after the step: 18.0
wind_ramp_abs 714 715 0 1; wsp. after the step: 19.0
wind_ramp_abs 755 756 0 1; wsp. after the step: 20.0
wind_ramp_abs 796 797 0 1; wsp. after the step: 21.0
wind_ramp_abs 837 838 0 1; wsp. after the step: 22.0
wind_ramp_abs 878 879 0 1; wsp. after the step: 23.0
wind_ramp_abs 919 920 0 1; wsp. after the step: 24.0
wind_ramp_abs 960 961 0 1; wsp. after the step: 25.0
;
begin tower_shadow_potential_2;
tower_mbdy_link tower;
nsec 2;
radius 0 4.15;
radius 115.63 2.75;
end tower_shadow_potential_2;
end wind; ;
begin aerodrag;
begin aerodrag_element;
mbdy_name tower;
aerodrag_sections uniform 10;
nsec 2;
sec 0 0.6 8.3; tower bottom
sec 115.63 0.6 5.5; tower top
end aerodrag_element; ;
begin aerodrag_element; Nacelle drag side
mbdy_name shaft;
aerodrag_sections uniform 2;
nsec 2;
sec 0 0.8 10;
sec 7.01 0.8 10;
end aerodrag_element;
end aerodrag; ;
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/dtu_10mw_rwt_ae.dat;
pc_filename ./data/dtu_10mw_rwt_pc.dat;
induction_method 1; 0=none, 1=normal
aerocalc_method 1; 0=ingen aerodynamic, 1=med aerodynamic
aerosections 50; def. 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 dll; ;
begin type2_dll;
name risoe_controller;
filename ./control/dtu_we_controller.dll;
dll_subroutine_init init_regulation;
dll_subroutine_update update_regulation;
arraysizes_init 52 1;
arraysizes_update 12 100;
begin init; ; Overall parameters
constant 1 10000; Rated power [kW]
constant 2 0.5236; Minimum rotor speed [rad/s]
constant 3 1.005; Rated rotor speed [rad/s]
constant 4 15600000; Maximum allowable generator torque [Nm]
constant 5 100; Minimum pitch angle, theta_min [deg],
; if |theta_min|>90, then a table of <wsp,theta_min> is read ;
; from a file named 'wptable.n', where n=int(theta_min)
constant 6 82; Maximum pitch angle [deg]
constant 7 10; Maximum pitch velocity operation [deg/s]
constant 8 0.4; Frequency of generator speed filter [Hz]
constant 9 0.7; Damping ratio of speed filter [-]
constant 10 1.92; Frequency of free-free DT torsion mode [Hz], if zero no notch filter used
; Partial load control parameters
constant 11 11750000; Optimal Cp tracking K factor [Nm/(rad/s)^2], ;
; Qg=K*Omega^2, K=eta*0.5*rho*A*Cp_opt*R^3/lambda_opt^3
constant 12 70840000; Proportional gain of torque controller [Nm/(rad/s)]
constant 13 15900000; Integral gain of torque controller [Nm/rad]
constant 14 0; Differential gain of torque controller [Nm/(rad/s^2)]
; Full load control parameters
constant 15 2; Generator control switch [1=constant power, 2=constant torque]
constant 16 1.304; Proportional gain of pitch controller [rad/(rad/s)]
constant 17 0.3511; Integral gain of pitch controller [rad/rad]
constant 18 0; Differential gain of pitch controller [rad/(rad/s^2)]
constant 19 4e-09; Proportional power error gain [rad/W]
constant 20 4e-09; Integral power error gain [rad/(Ws)]
constant 21 11.35; Coefficient of linear term in aerodynamic gain scheduling, KK1 [deg]
constant 22 400.7; Coefficient of quadratic term in aerodynamic gain scheduling, KK2 [deg^2] &
; (if zero, KK1 = pitch angle at double gain)
constant 23 1.3; Relative speed for double nonlinear gain [-]
; Cut-in simulation parameters
constant 24 -1; Cut-in time [s]
constant 25 1; Time delay for soft start of torque [1/1P]
; Cut-out simulation parameters
constant 26 1500; Cut-out time [s]
constant 27 5; Time constant for linear torque cut-out [s]
constant 28 1; Stop type [1=normal, 2=emergency]
constant 29 1; Time delay for pitch stop after shut-down signal [s]
constant 30 3; Maximum pitch velocity during initial period of stop [deg/s]
constant 31 3; Time period of initial pitch stop phase [s] (maintains pitch speed specified in constant 30)
constant 32 4; Maximum pitch velocity during final phase of stop [deg/s]
; Expert parameters (keep default values unless otherwise given)
constant 33 2; Lower angle above lowest minimum pitch angle for switch [deg]
constant 34 2; Upper angle above lowest minimum pitch angle for switch [deg], if equal then hard switch
constant 35 95; Ratio between filtered speed and reference speed for fully open torque limits [%]
constant 36 2; Time constant of 1st order filter on wind speed used for minimum pitch [1/1P]
constant 37 1; Time constant of 1st order filter on pitch angle used for gain scheduling [1/1P]
; Drivetrain damper
constant 38 0; Proportional gain of active DT damper [Nm/(rad/s)], requires frequency in input 10
; Over speed
constant 39 25; Overspeed percentage before initiating turbine controller alarm (shut-down) [%]
; Additional non-linear pitch control term (not used when all zero)
constant 40 0; Err0 [rad/s]
constant 41 0; ErrDot0 [rad/s^2]
constant 42 0; PitNonLin1 [rad/s]
; Storm control command
constant 43 28; Wind speed 'Vstorm' above which derating of rotor speed is used [m/s]
constant 44 28; Cut-out wind speed (only used for derating of rotor speed in storm) [m/s]
; Safety system parameters
constant 45 300; Overspeed percentage before initiating safety system alarm (shut-down) [%]
constant 46 1.5; Max low-pass filtered tower top acceleration level [m/s^2]
; Turbine parameter
constant 47 178; Nominal rotor diameter [m]
; Parameters for rotor inertia reduction in variable speed region
constant 48 0; Proportional gain on rotor acceleration in variable speed region [Nm/(rad/s^2)] (not used when zero)
; Parameters for alternative partial load controller with PI regulated TSR tracking
constant 49 0; Optimal tip speed ratio [-] (only used when K=constant 11 = 0 otherwise Qg=K*Omega^2 is used)
; Parameters for adding aerodynamic drivetrain damping on gain scheduling
constant 50 0; Proportional gain of aerodynamic DT damping [Nm/(rad/s)]
constant 51 0; Coefficient of linear term in aerodynamic DT damping scheduling, KK1 [deg]
constant 52 0; Coefficient of quadratic term in aerodynamic DT damping scheduling, KK2 [deg^2]
end init; ;
begin output;
general time; [s]
constraint bearing1 shaft_rot 1 only 2; Drivetrain speed [rad/s]
constraint bearing2 pitch1 1 only 1; [rad]
constraint bearing2 pitch2 1 only 1; [rad]
constraint bearing2 pitch3 1 only 1; [rad]
wind free_wind 1 0 0 -119; Global coordinates at hub height
dll inpvec 2 2; Elec. power from generator servo .dll
dll inpvec 2 8; Grid state flag from generator servo .dll
mbdy state acc tower 10 1 global only 1; Tower top x-acceleration [m/s^2]
mbdy state acc tower 10 1 global only 2; Tower top y-acceleration [m/s^2]
end output;
end type2_dll; ;
begin type2_dll;
name generator_servo;
filename ./control/generator_servo.dll;
dll_subroutine_init init_generator_servo;
dll_subroutine_update update_generator_servo;
arraysizes_init 7 1;
arraysizes_update 4 8;
begin init;
constant 1 20; Frequency of 2nd order servo model of generator-converter system [Hz]
constant 2 0.9; Damping ratio 2nd order servo model of generator-converter system [-]
constant 3 15600000; Maximum allowable LSS torque (pull-out torque) [Nm]
constant 4 0.94; Generator efficiency [-]
constant 5 1; Gearratio [-]
constant 6 0; Time for half value in softstart of torque [s]
constant 7 1500; Time for grid loss [s]
end init; ;
begin output;
general time; Time [s]
dll inpvec 1 1; Electrical torque reference [Nm]
constraint bearing1 shaft_rot 1 only 2; Generator LSS speed [rad/s]
mbdy momentvec shaft 1 1 shaft only 3; Shaft moment [kNm] (Qshaft)
end output;
;
begin actions;
mbdy moment_int shaft 1 -3 shaft towertop 2; Generator LSS torque [Nm]
end actions;
end type2_dll; ;
begin type2_dll;
name mech_brake;
filename ./control/mech_brake.dll;
dll_subroutine_init init_mech_brake;
dll_subroutine_update update_mech_brake;
arraysizes_init 7 1;
arraysizes_update 4 6;
begin init;
constant 1 2727252; Fully deployed maximum brake torque [Nm]
constant 2 100; Parameter alpha used in Q = tanh(omega*alpha), typically 1e2/Omega_nom
constant 3 0.625; Delay time for before brake starts to deploy [s] - from 5MW*1P_5/1P_10
constant 4 0.75; Time for brake to become fully deployed [s]
end init; ;
begin output;
general time; Time [s]
constraint bearing1 shaft_rot 1 only 2; Generator LSS speed [rad/s]
dll inpvec 1 25; Command to deploy mechanical disc brake [0,1]
end output;
;
begin actions;
mbdy moment_int shaft 1 3 shaft towertop 2; Brake LSS torque [Nm]
end actions;
end type2_dll; ;
begin type2_dll;
name servo_with_limits;
filename ./control/servo_with_limits.dll;
dll_subroutine_init init_servo_with_limits;
dll_subroutine_update update_servo_with_limits;
arraysizes_init 10 1;
arraysizes_update 5 9;
begin init;
constant 1 3; Number of blades [-]
constant 2 1; Frequency of 2nd order servo model of pitch system [Hz]
constant 3 0.7; Damping ratio 2nd order servo model of pitch system [-]
constant 4 10; Max. pitch speed [deg/s]
constant 5 15; Max. pitch acceleration [deg/s^2]
constant 6 -5; Min. pitch angle [deg]
constant 7 90; Max. pitch angle [deg]
constant 8 1500; Time for pitch runaway [s]
constant 9 1500; Time for stuck blade 1 [s]
constant 10 0; Angle of stuck blade 1 [deg]
end init;
begin output;
general time; Time [s]
dll inpvec 1 2; Pitch1 demand angle [rad]
dll inpvec 1 3; Pitch2 demand angle [rad]
dll inpvec 1 4; Pitch3 demand angle [rad]
dll inpvec 1 26; Flag for emergency pitch stop [0=off/1=on]
end output;
;
begin actions;
constraint bearing2 angle pitch1; Angle pitch1 bearing [rad]
constraint bearing2 angle pitch2; Angle pitch2 bearing [rad]
constraint bearing2 angle pitch3; Angle pitch3 bearing [rad]
end actions;
end type2_dll; ;
; --- DLL for tower-blade tip distance -- ;
begin type2_dll;
name disttowtip;
filename ./control/towclearsens.dll;
dll_subroutine_init initialize;
dll_subroutine_update update;
arraysizes_init 1 1;
arraysizes_update 12 4;
begin init;
constant 1 2.66; Tower radius close to downward blade tip [m]
end init;
begin output;
mbdy state pos tower 3 0.62 global; [1,2,3]. Tower position: 30.18 m
mbdy state pos blade1 26 1 global; [4,5,6]
mbdy state pos blade2 26 1 global; [7,8,9]
mbdy state pos blade3 26 1 global; [10,11,12]
end output;
end type2_dll;
end dll; ;----------------------------------------------------------------------------------------------------------------------------------------------------------------
;
begin output;
filename ./res/dtu_10mw_rwt; ; time 99.0 1000.0 ;
; data_format hawc_ascii;
data_format hawc_binary;
buffer 1; ;
general time;
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 -119; local wind at fixed position: coo (1=global,2=non-rotation rotor coo.), pos x, pos y, pos z
; Moments:
mbdy momentvec tower 1 1 tower # tower base;
mbdy momentvec tower 10 2 tower # tower yaw bearing;
mbdy momentvec shaft 4 1 shaft # main bearing;
mbdy momentvec blade1 2 2 blade1 # blade 1 root;
mbdy momentvec blade2 2 2 blade2 # blade 2 root;
mbdy momentvec blade3 2 2 blade3 # blade 3 root;
mbdy momentvec blade1 13 1 local # blade 1 50% local e coo;
mbdy momentvec blade2 13 1 local # blade 2 50% local e coo;
mbdy momentvec blade3 13 1 local # blade 3 50% local e coo;
; Displacements and accellerations
mbdy state pos tower 10 1 global only 1 # tower top fa displ;
mbdy state pos tower 10 1 global only 2 # tower top ss displ;
mbdy state acc tower 10 1 global only 1 # tower top fa acc;
mbdy state acc tower 10 1 global only 2 # tower top ss acc;
;
mbdy state pos blade1 26 1 blade1 # blade 1 tip pos;
mbdy state pos blade2 26 1 blade2 # blade 2 tip pos;
mbdy state pos blade3 26 1 blade3 # blade 3 tip pos;
mbdy state pos blade1 26 1 global # gl blade 1 tip pos;
; - Monitor Aerodynamics - ;
aero windspeed 3 1 1 72.5;
aero alfa 1 72.5;
aero alfa 2 72.5;
aero alfa 3 72.5;
aero cl 1 72.5;
aero cl 2 72.5;
aero cl 3 72.5;
aero cd 1 72.5;
aero cd 2 72.5;
aero cd 3 72.5;
; - Main Controller -
; Output to controller
; dll outvec 1 1 # time;
; dll outvec 1 2 # slow speed shaft rad/s;
; dll outvec 1 3 # pitch angle 1;
; dll outvec 1 4 # pitch angle 2;
; dll outvec 1 5 # pitch angle 3;
; dll outvec 1 6 # WSP_x_global;
; dll outvec 1 7 # WSP_y_global;
; dll outvec 1 8 # WSP_z_global;
; dll outvec 1 9 # Elec. pwr ;
; dll outvec 1 10 # Grid flag ;
; Input from controller
dll inpvec 1 1 # generator torque reference [nm];
dll inpvec 1 2 # pitch angle reference of blade 1 [rad];
dll inpvec 1 3 # pitch angle reference of blade 2 [rad];
dll inpvec 1 4 # pitch angle reference of blade 3 [rad];
dll inpvec 1 5 # power reference [w];
dll inpvec 1 6 # filtered wind speed [m/s];
dll inpvec 1 7 # filtered rotor speed [rad/s];
dll inpvec 1 8 # filtered rotor speed error for torque [rad/s];
dll inpvec 1 9 # bandpass filtered rotor speed [rad/s];
dll inpvec 1 10 # proportional term of torque contr. [nm];
dll inpvec 1 11 # integral term of torque controller [nm];
dll inpvec 1 12 # minimum limit of torque [nm];
dll inpvec 1 13 # maximum limit of torque [nm];
dll inpvec 1 14 # torque limit switch based on pitch [-];
dll inpvec 1 15 # filtered rotor speed error for pitch [rad/s];
dll inpvec 1 16 # power error for pitch [w];
dll inpvec 1 17 # proportional term of pitch controller [rad];
dll inpvec 1 18 # integral term of pitch controller [rad];
dll inpvec 1 19 # minimum limit of pitch [rad];
dll inpvec 1 20 # maximum limit of pitch [rad];
dll inpvec 1 21 # torque reference from dt dammper [nm];
dll inpvec 1 22 # status signal [-];
dll inpvec 1 23 # total added pitch rate [rad/s];
dll inpvec 1 24 # filtered mean pitch for gain sch [rad];
dll inpvec 1 25 # flag for mechnical brake [0=off/1=on];
dll inpvec 1 26 # flag for emergency pitch stop [0=off/1=on];
dll inpvec 1 27 # lp filtered acceleration level [m/s^2];
; ; Output to generator model
; dll outvec 2 1 # time ;
; dll outvec 2 2 # Electrical torque reference [Nm] ;
; dll outvec 2 3 # omega LSS ;
; Input from generator model
dll inpvec 2 1 # mgen lss [nm];
dll inpvec 2 2 # pelec w;
dll inpvec 2 3 # mframe;
dll inpvec 2 4 # mgen hss;
dll inpvec 2 5 # generator pmech kw;
dll inpvec 2 6 # filtered gen speed;
dll inpvec 2 7 # elec. pwr;
dll inpvec 2 8 # grid flag;
; Output to mechanical brake
dll inpvec 3 1 # brake torque [nm];
; ; Input from mechanical brake
; dll outvec 3 1 # Time [s] ;
; dll outvec 3 2 # Generator LSS speed [rad/s] ;
; dll outvec 3 3 # Deploy brake ;
; ; Output to pitch servo
; dll outvec 4 1 # time;
; dll outvec 4 2 # pitchref 1;
; dll outvec 4 3 # pitchref 2;
; dll outvec 4 4 # pitchref 3;
; dll outvec 4 5 # Emerg. stop;
; Input from pitch servo
dll inpvec 4 1 # pitch 1;
dll inpvec 4 2 # pitch 2;
dll inpvec 4 3 # pitch 3;
; Check tower clearence
dll inpvec 5 1 # bltip tow min d [m];
end output;
;
begin output_at_time aero 0.1;
filename res/at;
twist 1;
chord 1;
end output_at_time;
exit;
wetb/hawc2/tests/test_files/simulation_setup/DTU10MWRef6.0_IOS/input/htc/test_no_turb.htc 0 → 100644
View file @ 62f0f2e4
;DTU_10MW_RWT, version 5, 04-21-15, anyd
;
begin simulation;
time_stop 1;
solvertype 1; (newmark)
on_no_convergence continue;
convergence_limits 1000 1 1e-07;
logfile ../output/log/test_no_turb.log;
begin newmark;
deltat 0.02;
end newmark;
end simulation; ;
;----------------------------------------------------------------------------------------------------------------------------------------------------------------
begin new_htc_structure; ; beam_output_file_name ./log/DTU_10MW_RWT_beam.dat; Optional - Calculated beam properties of the bodies are written to file
; body_output_file_name ./log/DTU_10MW_RWT_body.dat; Optional - Body initial position and orientation are written to file
; body_eigenanalysis_file_name ./eig/DTU_10MW_RWT_body_eigen.dat;
; structure_eigenanalysis_file_name ./eig/DTU_10MW_RWT_strc_eigen.dat ;
;-------------------------------------------------------------------------------------------------------------------------------
;-------------------------------------------------------------------------------------------------------------------------------
begin main_body; tower 115m
name tower;
type timoschenko;
nbodies 1;
node_distribution c2_def;
damping_posdef 0 0 0 0.00412 0.00412 0.00045; Mx My Mz Kx Ky Kz , M´s raises overall level, K´s raises high freguency level "tuned by Larh"
begin timoschenko_input;
filename ./data/dtu_10mw_rwt_tower_st.dat;
set 1 1;
end timoschenko_input;
begin c2_def; Definition of centerline (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;
name towertop;
type timoschenko;
nbodies 1;
node_distribution c2_def;
damping_posdef 0 0 0 0.007 0.007 0.007; "changed 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;
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;
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;
end timoschenko_input;
begin c2_def; Definition of centerline (main_body coordinates)
nsec 5;
sec 1 0 0 0 0; Tower top 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;
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;
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.8 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 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; set subset
end timoschenko_input;
begin c2_def; Definition of centerline (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;
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; initial position of node 1
body_eulerang 0 0 0;
end base; ;
begin relative;
body1 tower last;
body2 towertop 1;
body2_eulerang 0 0 0;
end relative; ;
begin relative;
body1 towertop last;
body2 shaft 1;
body2_eulerang 90 0 0;
body2_eulerang 5 0 0; 5 deg tilt angle
body2_eulerang 0 0 0;
mbdy2_ini_rotvec_d1 0 0 -1 0.7; mbdy2_ini_rotvec_d1 0.0 0.0 -1.0 [init_wr];
end relative; ;
begin relative;
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
end relative; ;
begin relative;
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
end relative; ;
begin relative;
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
end relative; ;
begin relative;
body1 hub1 last;
body2 blade1 1;
body2_eulerang 0 0 0;
end relative; ;
begin relative;
body1 hub2 last;
body2 blade2 1;
body2_eulerang 0 0 0;
end relative; ;
begin relative;
body1 hub3 last;
body2 blade3 1;
body2_eulerang 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 -1; 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 -1;
end bearing2; ;
begin bearing2;
name pitch2;
body1 hub2 last;
body2 blade2 1;
bearing_vector 2 0 0 -1;
end bearing2; ;
begin bearing2;
name pitch3;
body1 hub3 last;
body2 blade3 1;
bearing_vector 2 0 0 -1;
end bearing2;
end constraint; ;
end new_htc_structure; ;----------------------------------------------------------------------------------------------------------------------------------------------------------------
begin wind;
density 1.225;
wsp 4;
tint 0;
horizontal_input 1;
windfield_rotations 0 0 0; yaw, tilt, rotation
center_pos0 0 0 -119; hub heigth
shear_format 1 0.2;
turb_format 0; 0=none, 1=mann,2=flex
tower_shadow_method 0; 0=none, 1=potential flow, 2=jet
scale_time_start 0;
wind_ramp_factor 0 40 0.6 1; ; Steps ;
wind_ramp_abs 140 141 0 1; wsp. after the step: 5.0
wind_ramp_abs 181 182 0 1; wsp. after the step: 6.0
wind_ramp_abs 222 223 0 1; wsp. after the step: 7.0
wind_ramp_abs 263 264 0 1; wsp. after the step: 8.0
wind_ramp_abs 304 305 0 1; wsp. after the step: 9.0
wind_ramp_abs 345 346 0 1; wsp. after the step: 10.0
wind_ramp_abs 386 387 0 1; wsp. after the step: 11.0
wind_ramp_abs 427 428 0 1; wsp. after the step: 12.0
wind_ramp_abs 468 469 0 1; wsp. after the step: 13.0
wind_ramp_abs 509 510 0 1; wsp. after the step: 14.0
wind_ramp_abs 550 551 0 1; wsp. after the step: 15.0
wind_ramp_abs 591 592 0 1; wsp. after the step: 16.0
wind_ramp_abs 632 633 0 1; wsp. after the step: 17.0
wind_ramp_abs 673 674 0 1; wsp. after the step: 18.0
wind_ramp_abs 714 715 0 1; wsp. after the step: 19.0
wind_ramp_abs 755 756 0 1; wsp. after the step: 20.0
wind_ramp_abs 796 797 0 1; wsp. after the step: 21.0
wind_ramp_abs 837 838 0 1; wsp. after the step: 22.0
wind_ramp_abs 878 879 0 1; wsp. after the step: 23.0
wind_ramp_abs 919 920 0 1; wsp. after the step: 24.0
wind_ramp_abs 960 961 0 1; wsp. after the step: 25.0
;
begin tower_shadow_potential_2;
tower_mbdy_link tower;
nsec 2;
radius 0 4.15;
radius 115.63 2.75;
end tower_shadow_potential_2;
end wind; ;
begin aerodrag;
begin aerodrag_element;
mbdy_name tower;
aerodrag_sections uniform 10;
nsec 2;
sec 0 0.6 8.3; tower bottom
sec 115.63 0.6 5.5; tower top
end aerodrag_element; ;
begin aerodrag_element; Nacelle drag side
mbdy_name shaft;
aerodrag_sections uniform 2;
nsec 2;
sec 0 0.8 10;
sec 7.01 0.8 10;
end aerodrag_element;
end aerodrag; ;
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/dtu_10mw_rwt_ae.dat;
pc_filename ./data/dtu_10mw_rwt_pc.dat;
induction_method 1; 0=none, 1=normal
aerocalc_method 1; 0=ingen aerodynamic, 1=med aerodynamic
aerosections 50; def. 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 dll; ;
; begin type2_dll;
; name risoe_controller;
; filename ./contro_controller.dll;
; dll_subroutine_init init_regulation;
; dll_subroutine_update update_regulation;
; arraysizes_init 52 1;
; arraysizes_update 12 100;
; begin init; ; Overall parameters
; constant 1 10000; Rated power [kW]
; constant 2 0.5236; Minimum rotor speed [rad/s]
; constant 3 1.005; Rated rotor speed [rad/s]
; constant 4 15600000; Maximum allowable generator torque [Nm]
; constant 5 100; Minimum pitch angle, theta_min [deg],
; ; if |theta_min|>90, then a table of <wsp,theta_min> is read ;
; ; from a file named 'wptable.n', where n=int(theta_min)
; constant 6 82; Maximum pitch angle [deg]
; constant 7 10; Maximum pitch velocity operation [deg/s]
; constant 8 0.4; Frequency of generator speed filter [Hz]
; constant 9 0.7; Damping ratio of speed filter [-]
; constant 10 1.92; Frequency of free-free DT torsion mode [Hz], if zero no notch filter used
; ; Partial load control parameters
; constant 11 11750000; Optimal Cp tracking K factor [Nm/(rad/s)^2], ;
; ; Qg=K*Omega^2, K=eta*0.5*rho*A*Cp_opt*R^3/lambda_opt^3
; constant 12 70840000; Proportional gain of torque controller [Nm/(rad/s)]
; constant 13 15900000; Integral gain of torque controller [Nm/rad]
; constant 14 0; Differential gain of torque controller [Nm/(rad/s^2)]
; ; Full load control parameters
; constant 15 2; Generator control switch [1=constant power, 2=constant torque]
; constant 16 1.304; Proportional gain of pitch controller [rad/(rad/s)]
; constant 17 0.3511; Integral gain of pitch controller [rad/rad]
; constant 18 0; Differential gain of pitch controller [rad/(rad/s^2)]
; constant 19 4e-09; Proportional power error gain [rad/W]
; constant 20 4e-09; Integral power error gain [rad/(Ws)]
; constant 21 11.35; Coefficient of linear term in aerodynamic gain scheduling, KK1 [deg]
; constant 22 400.7; Coefficient of quadratic term in aerodynamic gain scheduling, KK2 [deg^2] &
; ; (if zero, KK1 = pitch angle at double gain)
; constant 23 1.3; Relative speed for double nonlinear gain [-]
; ; Cut-in simulation parameters
; constant 24 -1; Cut-in time [s]
; constant 25 1; Time delay for soft start of torque [1/1P]
; ; Cut-out simulation parameters
; constant 26 1500; Cut-out time [s]
; constant 27 5; Time constant for linear torque cut-out [s]
; constant 28 1; Stop type [1=normal, 2=emergency]
; constant 29 1; Time delay for pitch stop after shut-down signal [s]
; constant 30 3; Maximum pitch velocity during initial period of stop [deg/s]
; constant 31 3; Time period of initial pitch stop phase [s] (maintains pitch speed specified in constant 30)
; constant 32 4; Maximum pitch velocity during final phase of stop [deg/s]
; ; Expert parameters (keep default values unless otherwise given)
; constant 33 2; Lower angle above lowest minimum pitch angle for switch [deg]
; constant 34 2; Upper angle above lowest minimum pitch angle for switch [deg], if equal then hard switch
; constant 35 95; Ratio between filtered speed and reference speed for fully open torque limits [%]
; constant 36 2; Time constant of 1st order filter on wind speed used for minimum pitch [1/1P]
; constant 37 1; Time constant of 1st order filter on pitch angle used for gain scheduling [1/1P]
; ; Drivetrain damper
; constant 38 0; Proportional gain of active DT damper [Nm/(rad/s)], requires frequency in input 10
; ; Over speed
; constant 39 25; Overspeed percentage before initiating turbine controller alarm (shut-down) [%]
; ; Additional non-linear pitch control term (not used when all zero)
; constant 40 0; Err0 [rad/s]
; constant 41 0; ErrDot0 [rad/s^2]
; constant 42 0; PitNonLin1 [rad/s]
; ; Storm control command
; constant 43 28; Wind speed 'Vstorm' above which derating of rotor speed is used [m/s]
; constant 44 28; Cut-out wind speed (only used for derating of rotor speed in storm) [m/s]
; ; Safety system parameters
; constant 45 300; Overspeed percentage before initiating safety system alarm (shut-down) [%]
; constant 46 1.5; Max low-pass filtered tower top acceleration level [m/s^2]
; ; Turbine parameter
; constant 47 178; Nominal rotor diameter [m]
; ; Parameters for rotor inertia reduction in variable speed region
; constant 48 0; Proportional gain on rotor acceleration in variable speed region [Nm/(rad/s^2)] (not used when zero)
; ; Parameters for alternative partial load controller with PI regulated TSR tracking
; constant 49 0; Optimal tip speed ratio [-] (only used when K=constant 11 = 0 otherwise Qg=K*Omega^2 is used)
; ; Parameters for adding aerodynamic drivetrain damping on gain scheduling
; constant 50 0; Proportional gain of aerodynamic DT damping [Nm/(rad/s)]
; constant 51 0; Coefficient of linear term in aerodynamic DT damping scheduling, KK1 [deg]
; constant 52 0; Coefficient of quadratic term in aerodynamic DT damping scheduling, KK2 [deg^2]
; end init; ;
; begin output;
; general time; [s]
; constraint bearing1 shaft_rot 1 only 2; Drivetrain speed [rad/s]
; constraint bearing2 pitch1 1 only 1; [rad]
; constraint bearing2 pitch2 1 only 1; [rad]
; constraint bearing2 pitch3 1 only 1; [rad]
; wind free_wind 1 0 0 -119; Global coordinates at hub height
; dll inpvec 2 2; Elec. power from generator servo .dll
; dll inpvec 2 8; Grid state flag from generator servo .dll
; mbdy state acc tower 10 1 global only 1; Tower top x-acceleration [m/s^2]
; mbdy state acc tower 10 1 global only 2; Tower top y-acceleration [m/s^2]
; end output;
; end type2_dll; ;
; begin type2_dll;
; name generator_servo;
; filename ./control/generator_servo.dll;
; dll_subroutine_init init_generator_servo;
; dll_subroutine_update update_generator_servo;
; arraysizes_init 7 1;
; arraysizes_update 4 8;
; begin init;
; constant 1 20; Frequency of 2nd order servo model of generator-converter system [Hz]
; constant 2 0.9; Damping ratio 2nd order servo model of generator-converter system [-]
; constant 3 15600000; Maximum allowable LSS torque (pull-out torque) [Nm]
; constant 4 0.94; Generator efficiency [-]
; constant 5 1; Gearratio [-]
; constant 6 0; Time for half value in softstart of torque [s]
; constant 7 1500; Time for grid loss [s]
; end init; ;
; begin output;
; general time; Time [s]
; dll inpvec 1 1; Electrical torque reference [Nm]
; constraint bearing1 shaft_rot 1 only 2; Generator LSS speed [rad/s]
; mbdy momentvec shaft 1 1 shaft only 3; Shaft moment [kNm] (Qshaft)
; end output;
; ;
; begin actions;
; mbdy moment_int shaft 1 -3 shaft towertop 2; Generator LSS torque [Nm]
; end actions;
; end type2_dll; ;
; begin type2_dll;
; name mech_brake;
; filename ./control/mech_brake.dll;
; dll_subroutine_init init_mech_brake;
; dll_subroutine_update update_mech_brake;
; arraysizes_init 7 1;
; arraysizes_update 4 6;
; begin init;
; constant 1 2727252; Fully deployed maximum brake torque [Nm]
; constant 2 100; Parameter alpha used in Q = tanh(omega*alpha), typically 1e2/Omega_nom
; constant 3 0.625; Delay time for before brake starts to deploy [s] - from 5MW*1P_5/1P_10
; constant 4 0.75; Time for brake to become fully deployed [s]
; end init; ;
; begin output;
; general time; Time [s]
; constraint bearing1 shaft_rot 1 only 2; Generator LSS speed [rad/s]
; dll inpvec 1 25; Command to deploy mechanical disc brake [0,1]
; end output;
; ;
; begin actions;
; mbdy moment_int shaft 1 3 shaft towertop 2; Brake LSS torque [Nm]
; end actions;
; end type2_dll; ;
; begin type2_dll;
; name servo_with_limits;
; filename ./control/servo_with_limits.dll;
; dll_subroutine_init init_servo_with_limits;
; dll_subroutine_update update_servo_with_limits;
; arraysizes_init 10 1;
; arraysizes_update 5 9;
; begin init;
; constant 1 3; Number of blades [-]
; constant 2 1; Frequency of 2nd order servo model of pitch system [Hz]
; constant 3 0.7; Damping ratio 2nd order servo model of pitch system [-]
; constant 4 10; Max. pitch speed [deg/s]
; constant 5 15; Max. pitch acceleration [deg/s^2]
; constant 6 -5; Min. pitch angle [deg]
; constant 7 90; Max. pitch angle [deg]
; constant 8 1500; Time for pitch runaway [s]
; constant 9 1500; Time for stuck blade 1 [s]
; constant 10 0; Angle of stuck blade 1 [deg]
; end init;
; begin output;
; general time; Time [s]
; dll inpvec 1 2; Pitch1 demand angle [rad]
; dll inpvec 1 3; Pitch2 demand angle [rad]
; dll inpvec 1 4; Pitch3 demand angle [rad]
; dll inpvec 1 26; Flag for emergency pitch stop [0=off/1=on]
; end output;
; ;
; begin actions;
; constraint bearing2 angle pitch1; Angle pitch1 bearing [rad]
; constraint bearing2 angle pitch2; Angle pitch2 bearing [rad]
; constraint bearing2 angle pitch3; Angle pitch3 bearing [rad]
; end actions;
; end type2_dll; ;
; ; --- DLL for tower-blade tip distance -- ;
; begin type2_dll;
; name disttowtip;
; filename ./control/towclearsens.dll;
; dll_subroutine_init initialize;
; dll_subroutine_update update;
; arraysizes_init 1 1;
; arraysizes_update 12 4;
; begin init;
; constant 1 2.66; Tower radius close to downward blade tip [m]
; end init;
; begin output;
; mbdy state pos tower 3 0.62 global; [1,2,3]. Tower position: 30.18 m
; mbdy state pos blade1 26 1 global; [4,5,6]
; mbdy state pos blade2 26 1 global; [7,8,9]
; mbdy state pos blade3 26 1 global; [10,11,12]
; end output;
; end type2_dll;
; end dll; ;----------------------------------------------------------------------------------------------------------------------------------------------------------------
;
begin output;
filename ../output/res/test_no_turb; ; time 99.0 1000.0 ;
; data_format hawc_ascii;
data_format hawc_binary;
buffer 1; ;
general time;
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 -119; local wind at fixed position: coo (1=global,2=non-rotation rotor coo.), pos x, pos y, pos z
; Moments:
mbdy momentvec tower 1 1 tower # tower base;
mbdy momentvec tower 10 2 tower # tower yaw bearing;
mbdy momentvec shaft 4 1 shaft # main bearing;
mbdy momentvec blade1 2 2 blade1 # blade 1 root;
mbdy momentvec blade2 2 2 blade2 # blade 2 root;
mbdy momentvec blade3 2 2 blade3 # blade 3 root;
mbdy momentvec blade1 13 1 local # blade 1 50% local e coo;
mbdy momentvec blade2 13 1 local # blade 2 50% local e coo;
mbdy momentvec blade3 13 1 local # blade 3 50% local e coo;
; Displacements and accellerations
mbdy state pos tower 10 1 global only 1 # tower top fa displ;
mbdy state pos tower 10 1 global only 2 # tower top ss displ;
mbdy state acc tower 10 1 global only 1 # tower top fa acc;
mbdy state acc tower 10 1 global only 2 # tower top ss acc;
;
mbdy state pos blade1 26 1 blade1 # blade 1 tip pos;
mbdy state pos blade2 26 1 blade2 # blade 2 tip pos;
mbdy state pos blade3 26 1 blade3 # blade 3 tip pos;
mbdy state pos blade1 26 1 global # gl blade 1 tip pos;
; - Monitor Aerodynamics - ;
aero windspeed 3 1 1 72.5;
aero alfa 1 72.5;
aero alfa 2 72.5;
aero alfa 3 72.5;
aero cl 1 72.5;
aero cl 2 72.5;
aero cl 3 72.5;
aero cd 1 72.5;
aero cd 2 72.5;
aero cd 3 72.5;
; - Main Controller -
; Output to controller
; dll outvec 1 1 # time;
; dll outvec 1 2 # slow speed shaft rad/s;
; dll outvec 1 3 # pitch angle 1;
; dll outvec 1 4 # pitch angle 2;
; dll outvec 1 5 # pitch angle 3;
; dll outvec 1 6 # WSP_x_global;
; dll outvec 1 7 # WSP_y_global;
; dll outvec 1 8 # WSP_z_global;
; dll outvec 1 9 # Elec. pwr ;
; dll outvec 1 10 # Grid flag ;
; Input from controller
dll inpvec 1 1 # generator torque reference [nm];
dll inpvec 1 2 # pitch angle reference of blade 1 [rad];
dll inpvec 1 3 # pitch angle reference of blade 2 [rad];
dll inpvec 1 4 # pitch angle reference of blade 3 [rad];
dll inpvec 1 5 # power reference [w];
dll inpvec 1 6 # filtered wind speed [m/s];
dll inpvec 1 7 # filtered rotor speed [rad/s];
dll inpvec 1 8 # filtered rotor speed error for torque [rad/s];
dll inpvec 1 9 # bandpass filtered rotor speed [rad/s];
dll inpvec 1 10 # proportional term of torque contr. [nm];
dll inpvec 1 11 # integral term of torque controller [nm];
dll inpvec 1 12 # minimum limit of torque [nm];
dll inpvec 1 13 # maximum limit of torque [nm];
dll inpvec 1 14 # torque limit switch based on pitch [-];
dll inpvec 1 15 # filtered rotor speed error for pitch [rad/s];
dll inpvec 1 16 # power error for pitch [w];
dll inpvec 1 17 # proportional term of pitch controller [rad];
dll inpvec 1 18 # integral term of pitch controller [rad];
dll inpvec 1 19 # minimum limit of pitch [rad];
dll inpvec 1 20 # maximum limit of pitch [rad];
dll inpvec 1 21 # torque reference from dt dammper [nm];
dll inpvec 1 22 # status signal [-];
dll inpvec 1 23 # total added pitch rate [rad/s];
dll inpvec 1 24 # filtered mean pitch for gain sch [rad];
dll inpvec 1 25 # flag for mechnical brake [0=off/1=on];
dll inpvec 1 26 # flag for emergency pitch stop [0=off/1=on];
dll inpvec 1 27 # lp filtered acceleration level [m/s^2];
; ; Output to generator model
; dll outvec 2 1 # time ;
; dll outvec 2 2 # Electrical torque reference [Nm] ;
; dll outvec 2 3 # omega LSS ;
; Input from generator model
dll inpvec 2 1 # mgen lss [nm];
dll inpvec 2 2 # pelec w;
dll inpvec 2 3 # mframe;
dll inpvec 2 4 # mgen hss;
dll inpvec 2 5 # generator pmech kw;
dll inpvec 2 6 # filtered gen speed;
dll inpvec 2 7 # elec. pwr;
dll inpvec 2 8 # grid flag;
; Output to mechanical brake
dll inpvec 3 1 # brake torque [nm];
; ; Input from mechanical brake
; dll outvec 3 1 # Time [s] ;
; dll outvec 3 2 # Generator LSS speed [rad/s] ;
; dll outvec 3 3 # Deploy brake ;
; ; Output to pitch servo
; dll outvec 4 1 # time;
; dll outvec 4 2 # pitchref 1;
; dll outvec 4 3 # pitchref 2;
; dll outvec 4 4 # pitchref 3;
; dll outvec 4 5 # Emerg. stop;
; Input from pitch servo
dll inpvec 4 1 # pitch 1;
dll inpvec 4 2 # pitch 2;
dll inpvec 4 3 # pitch 3;
; Check tower clearence
dll inpvec 5 1 # bltip tow min d [m];
end output;
;
exit;
wetb/hawc2/tests/test_hawc2_input_writer.py 0 → 100644
View file @ 62f0f2e4
import os
import pandas as pd
from pathlib import Path
import pytest
from wetb.hawc2.tests import test_files
from wetb.hawc2.hawc2_input_writer import HAWC2InputWriter
from wetb.hawc2.htc_file import HTCFile
from wetb.hawc2 import hawc2_input_writer
from tests.run_main import run_module_main
htc_base_dir = os.path.dirname(test_files.__file__) + '/simulation_setup/DTU10MWRef6.0/htc/'
@pytest.fixture
def h2writer():
"""HAWC2InputWriter with DTU 10 MW as base file, used in several tests"""
htc_base_file = htc_base_dir + 'DTU_10MW_RWT.htc'
return HAWC2InputWriter(htc_base_file)
def test_hawc2_writer_frompath(tmp_path):
"""Load from a path base, and write to path and string"""
path = Path(htc_base_dir + 'DTU_10MW_RWT.htc')
h2writer = HAWC2InputWriter(path)
h2writer.write(tmp_path / 'test.htc') # path
h2writer.write((tmp_path / 'test.htc').as_posix()) # string
def test_hawc2_writer_noname_error(h2writer, tmp_path):
"""Should throw error if set_name is called without name"""
wsp_lst = [4, 6, 8]
df = pd.DataFrame({'wind.wsp': wsp_lst, 'Folder': ['tmp'] * 3})
h2writer.from_pandas(df)
with pytest.raises(KeyError): # keyerror in write_all
h2writer.write_all(tmp_path)
def test_pandas2htc_dottag(h2writer, tmp_path):
"""Load dlc from dict + pandas dataframe (wind speed tag with dot, name), write htc"""
wsp_lst = [4, 6, 8]
df = pd.DataFrame({'wind.wsp': wsp_lst, 'Name': ['c1', 'c2', 'c3'], 'Folder': ['tmp'] * 3})
h2writer.from_pandas(df)
h2writer.write_all(tmp_path)
for i, wsp in enumerate(wsp_lst, 1):
htc_path = tmp_path / ('tmp/c%d.htc' % i)
htc = HTCFile(htc_path, modelpath=tmp_path.as_posix())
assert htc.wind.wsp[0] == wsp
assert htc.simulation.logfile[0] == ('./log/tmp/c%d.log' % i)
assert htc.output.filename[0] == ('./res/tmp/c%d' % i)
def test_excel2htc(h2writer, tmp_path):
"""Load dlc from Excel (wind speed tag with dot, name), write to htc.
Both string and pathlib.Path"""
excel_path = os.path.dirname(test_files.__file__) + '/htc_input_table.xlsx'
for i in range(2):
h2writer.from_excel(excel_path if i else Path(excel_path))
# print(h2writer.contents)
h2writer.write_all(tmp_path)
for i, wsp in enumerate([4, 6], 1):
htc = HTCFile(tmp_path / ('tmp/a%d.htc' % i), modelpath=tmp_path.as_posix())
assert htc.wind.wsp[0] == wsp
assert htc.simulation.logfile[0] == ('./log/tmp/a%d.log' % i)
assert htc.output.filename[0] == ('./res/tmp/a%d' % i)
def test_hawc2_writer_main():
"""Run the example in the hawc2_input_writer module"""
run_module_main(hawc2_input_writer)
def test_hawc2_writer_custom(tmp_path):
"""Input writer with a custom 'set_'-style attribute"""
htc_base_file = htc_base_dir + 'DTU_10MW_RWT.htc'
class MyWriter(HAWC2InputWriter):
def set_time(self, htc, **kwargs):
htc.set_time(self.time_start, self.time_start + kwargs['time'])
myWriter = MyWriter(htc_base_file, time_start=100)
myWriter.write(tmp_path / 'w1.htc', **{'Name': 'w1', 'wind.wsp': 10, 'time': 600})
htc = HTCFile(tmp_path / 'w1.htc', modelpath=tmp_path.as_posix())
assert htc.simulation.time_stop[0] == 700
assert htc.wind.wsp[0] == 10
assert htc.output.filename[0] == './res/w1'
def test_CVF2pandas(h2writer):
"""Load dlc with constants, variables, functions format"""
constants = {'simulation.time_stop': 100}
variables = {'wind.wsp': [4, 6, 8],
'wind.tint': [0.1, 0.15, 0.2]}
functions = {'Name': lambda x: 'sim_wsp' + str(x['wind.wsp']) + '_ti' + str(x['wind.tint'])}
h2writer.from_CVF(constants, variables, functions)
assert len(h2writer.contents) == 9
assert set(list(h2writer.contents)) == set(['simulation.time_stop', 'wind.wsp', 'wind.tint', 'Name'])
wetb/hawc2/tests/test_hawc2io.py
View file @ 62f0f2e4
......@@ -3,12 +3,6 @@ Created on 05/11/2015
@author: MMPE
'''
from __future__ import unicode_literals
from __future__ import print_function
from __future__ import division
from __future__ import absolute_import
from future import standard_library
standard_library.install_aliases()
import unittest
import numpy as np
from wetb.hawc2.Hawc2io import ReadHawc2
......
wetb/hawc2/tests/test_htc_extensions.py
View file @ 62f0f2e4
......@@ -3,15 +3,6 @@ Created on 17/07/2014
@author: MMPE
'''
from __future__ import print_function
from __future__ import unicode_literals
from __future__ import division
from __future__ import absolute_import
from io import open
from builtins import str
from builtins import zip
from future import standard_library
standard_library.install_aliases()
import os
import unittest
......
wetb/hawc2/tests/test_htc_file.py
View file @ 62f0f2e4
'''
Created on 17/07/2014
@author: MMPE
'''
from __future__ import print_function
from __future__ import unicode_literals
from __future__ import division
from __future__ import absolute_import
from io import open
from builtins import str
from builtins import zip
from future import standard_library
standard_library.install_aliases()
import os
import unittest
from datetime import datetime
from wetb.hawc2.htc_file import HTCFile, HTCLine
import numpy as np
class TestHtcFile(unittest.TestCase):
def setUp(self):
unittest.TestCase.setUp(self)
self.testfilepath = os.path.join(os.path.dirname(__file__), 'test_files/htcfiles/') # test file path
def check_htc_file(self, f):
with open(f) as fid:
orglines = fid.read().strip().split("\n")
htcfile = HTCFile(f,"../")
newlines = str(htcfile).split("\n")
htcfile.save(self.testfilepath + 'tmp.htc')
#with open(self.testfilepath + 'tmp.htc') as fid:
# newlines = fid.readlines()
for i, (org, new) in enumerate(zip(orglines, newlines), 1):
fmt = lambda x : x.strip().replace("\t", " ").replace(" ", " ").replace(" ", " ").replace(" ", " ").replace(" ", " ")
if fmt(org) != fmt(new):
print ("----------%d-------------" % i)
print (fmt(org))
print (fmt(new))
self.assertEqual(fmt(org), fmt(new))
break
print ()
assert len(orglines) == len(newlines)
def test_htc_files(self):
for f in ['test3.htc']:
self.check_htc_file(self.testfilepath + f)
def test_htc_file_get(self):
htcfile = HTCFile(self.testfilepath + "test3.htc",'../')
self.assertEqual(htcfile['simulation']['time_stop'][0], 200)
self.assertEqual(htcfile['simulation/time_stop'][0], 200)
self.assertEqual(htcfile['simulation.time_stop'][0], 200)
self.assertEqual(htcfile.simulation.time_stop[0], 200)
self.assertEqual(htcfile.dll.type2_dll.name[0], "risoe_controller")
self.assertEqual(htcfile.dll.type2_dll__2.name[0], "risoe_controller2")
s = """begin simulation;\n time_stop\t200;"""
self.assertEqual(str(htcfile.simulation)[:len(s)], s)
def test_htc_file_get2(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
self.assertEqual(htcfile['simulation']['logfile'][0], './logfiles/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.log')
self.assertEqual(htcfile['simulation/logfile'][0], './logfiles/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.log')
self.assertEqual(htcfile['simulation.logfile'][0], './logfiles/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.log')
self.assertEqual(htcfile.simulation.logfile[0], './logfiles/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.log')
self.assertEqual(htcfile.simulation.newmark.deltat[0], 0.02)
def test_htc_file_set(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
time_stop = htcfile.simulation.time_stop[0]
htcfile.simulation.time_stop = time_stop * 2
self.assertEqual(htcfile.simulation.time_stop[0], 2 * time_stop)
self.assertEqual(htcfile.simulation.time_stop.__class__, HTCLine)
htcfile.output.time = 10, 20
self.assertEqual(htcfile.output.time[:2], [10, 20])
self.assertEqual(str(htcfile.output.time), "time\t10 20;\n")
htcfile.output.time = [11, 21]
self.assertEqual(htcfile.output.time[:2], [11, 21])
htcfile.output.time = "12 22"
self.assertEqual(htcfile.output.time[:2], [12, 22])
def test_htc_file_set_key(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
htcfile.simulation.name = "value"
self.assertEqual(htcfile.simulation.name[0], "value")
def test_htc_file_del_key(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
del htcfile.simulation.logfile
self.assertTrue("logfile" not in str(htcfile.simulation))
try:
del htcfile.hydro.water_properties.water_kinematics_dll
except KeyError:
pass
def test_htcfile_setname(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
htcfile.set_name("mytest")
self.assertEqual(os.path.relpath(htcfile.filename, self.testfilepath).replace("\\","/"), r'../htc/mytest.htc')
self.assertEqual(htcfile.simulation.logfile[0], './log/mytest.log')
self.assertEqual(htcfile.output.filename[0], './res/mytest')
htcfile.set_name("mytest", 'subfolder')
self.assertEqual(os.path.relpath(htcfile.filename, self.testfilepath).replace("\\","/"), r'../htc/subfolder/mytest.htc')
self.assertEqual(htcfile.simulation.logfile[0], './log/subfolder/mytest.log')
self.assertEqual(htcfile.output.filename[0], './res/subfolder/mytest')
def test_set_time(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
htcfile.set_time(10, 20, 0.2)
self.assertEqual(htcfile.simulation.time_stop[0], 20)
self.assertEqual(htcfile.simulation.newmark.deltat[0], 0.2)
self.assertEqual(htcfile.wind.scale_time_start[0], 10)
self.assertEqual(htcfile.output.time[:2], [10, 20])
def test_add_section(self):
htcfile = HTCFile()
htcfile.wind.add_section('mann')
htcfile.wind.mann.add_line("create_turb_parameters", [29.4, 1.0, 3.9, 1004, 1.0], "L, alfaeps, gamma, seed, highfrq compensation")
self.assertEqual(htcfile.wind.mann.create_turb_parameters[0], 29.4)
self.assertEqual(htcfile.wind.mann.create_turb_parameters[3], 1004)
self.assertEqual(htcfile.wind.mann.create_turb_parameters.comments, "L, alfaeps, gamma, seed, highfrq compensation")
def test_add_section2(self):
htcfile = HTCFile()
htcfile.add_section('hydro')
#self.assertEqual(str(htcfile).strip()[-5:], "exit;")
htcfile = HTCFile(self.testfilepath + "test.htc")
htcfile.add_section('hydro')
self.assertEqual(str(htcfile).strip()[-5:], "exit;")
def test_add_mann(self):
htcfile = HTCFile()
htcfile.add_mann_turbulence(30.1, 1.1, 3.3, 102, False)
s = """begin mann;
create_turb_parameters\t30.1 1.1 3.3 102 0;\tL, alfaeps, gamma, seed, highfrq compensation
filename_u\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.12x3.12_s0102u.turb;
filename_v\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.12x3.12_s0102v.turb;
filename_w\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.12x3.12_s0102w.turb;
box_dim_u\t16384 0.3662;
box_dim_v\t32 3.2258;
box_dim_w\t32 3.2258;"""
for a, b in zip(s.split("\n"), str(htcfile.wind.mann).split("\n")):
self.assertEqual(a.strip(), b.strip())
self.assertEqual(htcfile.wind.turb_format[0], 1)
self.assertEqual(htcfile.wind.turb_format.comments, "")
def test_add_turb_export(self):
htc = HTCFile()
htc.add_mann_turbulence(30.1, 1.1, 3.3, 102, False)
htc.set_time(100,700,0.01)
htc.add_turb_export()
s = """begin turb_export;
filename_u\texport_u.turb;
filename_v\texport_v.turb;
filename_w\texport_w.turb;
samplefrq\t3;
time_start\t100;
nsteps\t60000.0;
box_dim_v\t32 3.2258;
box_dim_w\t32 3.2258;
end turb_export;"""
for a, b in zip(s.split("\n"), str(htc.wind.turb_export).split("\n")):
self.assertEqual(a.strip(), b.strip())
def test_sensors(self):
htcfile = HTCFile()
htcfile.set_name("test")
htcfile.output.add_sensor('wind', 'free_wind', [1, 0, 0, -30])
s = """begin output;
filename\t./res/test;
general time;
wind free_wind\t1 0 0 -30;"""
for a, b in zip(s.split("\n"), str(htcfile.output).split("\n")):
self.assertEqual(a.strip(), b.strip())
#print (htcfile)
def test_output_at_time(self):
htcfile = HTCFile(self.testfilepath + "test2.htc",'../')
self.assertTrue('begin output_at_time aero 15.0;' in str(htcfile))
def test_output_files(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
output_files = htcfile.output_files()
#print (htcfile.output)
for f in ['./logfiles/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.log',
'./visualization/dlc12_wsp10_wdir000_s1004.hdf5',
'./animation/structure_aero_control_turb.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_wsp10_wdir000_s1004_beam.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_wsp10_wdir000_s1004_body.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_wsp10_wdir000_s1004_struct.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_wsp10_wdir000_s1004_body_eigen.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_wsp10_wdir000_s1004_strc_eigen.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/mode*.dat',
'./launcher_test/ssystem_eigenanalysis.dat', './launcher_test/mode*.dat',
'./res/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.sel',
'./res/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.dat',
'./res/rotor_check_inipos.dat',
'./res/rotor_check_inipos2.dat']:
try:
output_files.remove(f)
except ValueError:
raise ValueError(f + " is not in list")
self.assertFalse(output_files)
def test_turbulence_files(self):
htcfile = HTCFile(self.testfilepath + "dlc14_wsp10_wdir000_s0000.htc",'../')
self.assertEqual(htcfile.turbulence_files(), ['./turb/turb_wsp10_s0000u.bin', './turb/turb_wsp10_s0000v.bin', './turb/turb_wsp10_s0000w.bin'])
def test_input_files(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
input_files = htcfile.input_files()
#print (htcfile.output)
for f in ['./data/DTU_10MW_RWT_Tower_st.dat',
'./data/DTU_10MW_RWT_Towertop_st.dat',
'./data/DTU_10MW_RWT_Shaft_st.dat',
'./data/DTU_10MW_RWT_Hub_st.dat',
'./data/DTU_10MW_RWT_Blade_st.dat',
'./data/DTU_10MW_RWT_ae.dat',
'./data/DTU_10MW_RWT_pc.dat',
'./control/risoe_controller.dll',
'./control/generator_servo.dll',
'./control/mech_brake.dll',
'./control/servo_with_limits.dll',
'./control/towclearsens.dll',
'./data/user_shear.dat',
self.testfilepath.replace("\\","/") + 'test.htc'
]:
try:
input_files.remove(f)
except ValueError:
raise ValueError(f + " is not in list")
self.assertFalse(input_files)
htcfile = HTCFile(self.testfilepath + "DTU_10MW_RWT.htc")
self.assertTrue('./control/wpdata.100' in htcfile.input_files())
def test_input_files2(self):
htcfile = HTCFile(self.testfilepath + "ansi.htc",'../')
input_files = htcfile.input_files()
self.assertTrue('./htc_hydro/ireg_airy_h6_t10.inp' in input_files)
#
def test_continue_in_files(self):
htcfile = HTCFile(self.testfilepath + "continue_in_file.htc", ".")
self.assertIn('main_body__31', htcfile.new_htc_structure.keys())
self.assertIn(os.path.abspath(self.testfilepath + 'orientation.dat'), [os.path.abspath(f) for f in htcfile.input_files()])
self.assertIn('./data/NREL_5MW_st1.txt', htcfile.input_files())
self.assertEqual(str(htcfile).count("exit"), 1)
self.assertIn('filename\t./res/oc4_p2_load_case_eq;', str(htcfile))
def test_continue_in_files_autodetect_path(self):
htcfile = HTCFile(self.testfilepath + "sub/continue_in_file.htc")
self.assertIn('main_body__31', htcfile.new_htc_structure.keys())
self.assertIn(os.path.abspath(self.testfilepath + 'orientation.dat'), [os.path.abspath(f) for f in htcfile.input_files()])
self.assertIn('./data/NREL_5MW_st1.txt', htcfile.input_files())
self.assertEqual(str(htcfile).count("exit"), 1)
self.assertIn('filename\t./res/oc4_p2_load_case_eq;', str(htcfile))
def test_tjul_example(self):
htcfile = HTCFile(self.testfilepath + "./tjul.htc", ".")
htcfile.save("./temp.htc")
def test_ansi(self):
htcfile = HTCFile(self.testfilepath + "./ansi.htc",'../')
def test_file_with_BOM(self):
htcfile = HTCFile(self.testfilepath + 'DLC15_wsp11_wdir000_s0000_phi000_Free_v2_visual.htc')
self.assertEqual(str(htcfile)[0], ";")
def test_htc_reset(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
self.assertEqual(htcfile.wind.wsp[0], 10)
def test_htc_model_autodetect(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
self.assertEqual(os.path.relpath(htcfile.modelpath,os.path.dirname(htcfile.filename)), "..")
htcfile = HTCFile(self.testfilepath + "sub/test.htc")
self.assertEqual(os.path.relpath(htcfile.modelpath,os.path.dirname(htcfile.filename)).replace("\\","/"), "../..")
self.assertRaisesRegex(ValueError, "Modelpath cannot be autodetected", HTCFile, self.testfilepath + "test2.htc")
def test_open_eq_save(self):
HTCFile(self.testfilepath + "test3.htc","../").save(self.testfilepath + "tmp.htc")
htcfile = HTCFile(self.testfilepath + "tmp.htc","../")
htcfile.save(self.testfilepath + "tmp.htc")
self.assertEqual(str(htcfile).count("\t"), str(HTCFile(self.testfilepath + "tmp.htc", "../")).count("\t"))
self.assertEqual(str(htcfile), str(HTCFile(self.testfilepath + "tmp.htc", "../")))
def test_2xoutput(self):
htc = HTCFile(self.testfilepath + "test_2xoutput.htc","../")
self.assertEqual(len(htc.res_file_lst()), 4)
if __name__ == "__main__":
#import sys;sys.argv = ['', 'Test.testName']
unittest.main()
'''
Created on 17/07/2014
@author: MMPE
'''
import os
import unittest
from wetb.hawc2.htc_file import HTCFile, HTCLine
class TestHtcFile(unittest.TestCase):
def setUp(self):
unittest.TestCase.setUp(self)
self.testfilepath = os.path.join(os.path.dirname(__file__), 'test_files/htcfiles/') # test file path
def check_htc_file(self, f):
with open(f) as fid:
orglines = fid.read().strip().split("\n")
htcfile = HTCFile(f, "../")
newlines = str(htcfile).split("\n")
htcfile.save(self.testfilepath + 'tmp.htc')
# with open(self.testfilepath + 'tmp.htc') as fid:
# newlines = fid.readlines()
for i, (org, new) in enumerate(zip(orglines, newlines), 1):
def fmt(x): return x.strip().replace("\t", " ").replace(
" ", " ").replace(" ", " ").replace(" ", " ").replace(" ", " ")
if fmt(org) != fmt(new):
print("----------%d-------------" % i)
print(fmt(org))
print(fmt(new))
self.assertEqual(fmt(org), fmt(new))
break
print()
assert len(orglines) == len(newlines)
def test_htc_files(self):
for f in ['test3.htc']:
self.check_htc_file(self.testfilepath + f)
def test_htc_file_get(self):
htcfile = HTCFile(self.testfilepath + "test3.htc", '../')
self.assertEqual(htcfile['simulation']['time_stop'][0], 200)
self.assertEqual(htcfile['simulation/time_stop'][0], 200)
self.assertEqual(htcfile['simulation.time_stop'][0], 200)
self.assertEqual(htcfile.simulation.time_stop[0], 200)
self.assertEqual(htcfile.dll.type2_dll.name[0], "risoe_controller")
self.assertEqual(htcfile.dll.type2_dll__2.name[0], "risoe_controller2")
s = """begin simulation;\n time_stop\t200;"""
self.assertEqual(str(htcfile.simulation)[:len(s)], s)
def test_htc_file_get2(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
self.assertEqual(htcfile['simulation']['logfile'][0],
'./logfiles/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.log')
self.assertEqual(htcfile['simulation/logfile'][0],
'./logfiles/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.log')
self.assertEqual(htcfile['simulation.logfile'][0],
'./logfiles/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.log')
self.assertEqual(htcfile.simulation.logfile[0], './logfiles/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.log')
self.assertEqual(htcfile.simulation.newmark.deltat[0], 0.02)
def test_htc_file_set(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
time_stop = htcfile.simulation.time_stop[0]
htcfile.simulation.time_stop = time_stop * 2
self.assertEqual(htcfile.simulation.time_stop[0], 2 * time_stop)
self.assertEqual(htcfile.simulation.time_stop.__class__, HTCLine)
htcfile.output.time = 10, 20
self.assertEqual(htcfile.output.time[:2], [10, 20])
self.assertEqual(str(htcfile.output.time), "time\t10 20;\n")
htcfile.output.time = [11, 21]
self.assertEqual(htcfile.output.time[:2], [11, 21])
htcfile.output.time = "12 22"
self.assertEqual(htcfile.output.time[:2], [12, 22])
def test_htc_file_set_key(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
htcfile.simulation.name = "value"
self.assertEqual(htcfile.simulation.name[0], "value")
htcfile.simulation.name2 = ("value", 1)
self.assertEqual(htcfile.simulation.name2[0], "value")
self.assertEqual(htcfile.simulation.name2[1], 1)
def test_htc_file_set_key2(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
htcfile.simulation['name'] = "value"
self.assertEqual(htcfile.simulation.name[0], "value")
htcfile.simulation['name2'] = ("value", 1)
self.assertEqual(htcfile.simulation.name2[0], "value")
self.assertEqual(htcfile.simulation.name2[1], 1)
def test_htc_file_del_key(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
del htcfile.simulation.logfile
self.assertTrue("logfile" not in str(htcfile.simulation))
try:
del htcfile.hydro.water_properties.water_kinematics_dll
except KeyError:
pass
def test_htc_file_delete(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
self.assertTrue("logfile" in str(htcfile.simulation))
htcfile.simulation.logfile.delete()
self.assertTrue("logfile" not in str(htcfile.simulation))
self.assertTrue('newmark' in str(htcfile.simulation))
htcfile.simulation.newmark.delete()
with self.assertRaises(KeyError):
htcfile.simulation.newmark
def test_htcfile_setname(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
htcfile.set_name("mytest")
self.assertEqual(os.path.relpath(htcfile.filename, self.testfilepath).replace("\\", "/"), r'../htc/mytest.htc')
self.assertEqual(htcfile.simulation.logfile[0], './log/mytest.log')
self.assertEqual(htcfile.output.filename[0], './res/mytest')
htcfile.set_name("mytest", 'subfolder')
self.assertEqual(os.path.relpath(htcfile.filename, self.testfilepath).replace(
"\\", "/"), r'../htc/subfolder/mytest.htc')
self.assertEqual(htcfile.simulation.logfile[0], './log/subfolder/mytest.log')
self.assertEqual(htcfile.output.filename[0], './res/subfolder/mytest')
def test_set_time(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
htcfile.set_time(10, 20, 0.2)
self.assertEqual(htcfile.simulation.time_stop[0], 20)
self.assertEqual(htcfile.simulation.newmark.deltat[0], 0.2)
self.assertEqual(htcfile.wind.scale_time_start[0], 10)
self.assertEqual(htcfile.output.time[:2], [10, 20])
def test_add_section(self):
htcfile = HTCFile()
htcfile.wind.add_section('mann')
htcfile.wind.mann.add_line("create_turb_parameters", [
29.4, 1.0, 3.9, 1004, 1.0], "L, alfaeps, gamma, seed, highfrq compensation")
self.assertEqual(htcfile.wind.mann.create_turb_parameters[0], 29.4)
self.assertEqual(htcfile.wind.mann.create_turb_parameters[3], 1004)
self.assertEqual(htcfile.wind.mann.create_turb_parameters.comments,
"L, alfaeps, gamma, seed, highfrq compensation")
def test_add_section2(self):
htcfile = HTCFile()
htcfile.add_section('hydro')
#self.assertEqual(str(htcfile).strip()[-5:], "exit;")
htcfile = HTCFile(self.testfilepath + "test.htc")
htcfile.add_section('hydro')
self.assertEqual(str(htcfile).strip()[-5:], "exit;")
def test_add_mann(self):
htcfile = HTCFile()
htcfile.add_mann_turbulence(30.1, 1.1, 3.3, 102, False)
s = """begin mann;
create_turb_parameters\t30.1 1.1 3.3 102 0;\tL, alfaeps, gamma, seed, highfrq compensation
filename_u\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.12x3.12_s0102u.turb;
filename_v\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.12x3.12_s0102v.turb;
filename_w\t./turb/mann_l30.1_ae1.1000_g3.3_h0_16384x32x32_0.366x3.12x3.12_s0102w.turb;
box_dim_u\t16384 0.3662;
box_dim_v\t32 3.125;
box_dim_w\t32 3.125;"""
for a, b in zip(s.split("\n"), str(htcfile.wind.mann).split("\n")):
self.assertEqual(a.strip(), b.strip())
self.assertEqual(htcfile.wind.turb_format[0], 1)
self.assertEqual(htcfile.wind.turb_format.comments, "0=none, 1=mann,2=flex")
def test_add_turb_export(self):
htc = HTCFile()
htc.add_mann_turbulence(30.1, 1.1, 3.3, 102, False)
htc.set_time(100, 700, 0.01)
htc.add_turb_export()
s = """begin turb_export;
filename_u\texport_u.turb;
filename_v\texport_v.turb;
filename_w\texport_w.turb;
samplefrq\t3;
time_start\t100;
nsteps\t60000.0;
box_dim_v\t32 3.125;
box_dim_w\t32 3.125;
end turb_export;"""
for a, b in zip(s.split("\n"), str(htc.wind.turb_export).split("\n")):
self.assertEqual(a.strip(), b.strip())
def test_sensors(self):
htcfile = HTCFile()
htcfile.set_name("test")
htcfile.output.add_sensor('wind', 'free_wind', [1, 0, 0, -30])
s = """begin output;
filename\t./res/test;
data_format\tgtsdf;
general time;
wind free_wind\t1 0 0 -30;"""
for a, b in zip(s.split("\n"), str(htcfile.output).split("\n")):
self.assertEqual(a.strip(), b.strip())
#print (htcfile)
def test_output_at_time(self):
htcfile = HTCFile(self.testfilepath + "test2.htc", '../')
self.assertTrue('begin output_at_time aero 15.0;' in str(htcfile))
def test_output_files(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
output_files = htcfile.output_files()
#print (htcfile.output)
for f in ['./logfiles/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.log',
'./visualization/dlc12_wsp10_wdir000_s1004.hdf5',
'./animation/structure_aero_control_turb.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_wsp10_wdir000_s1004_beam.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_wsp10_wdir000_s1004_body.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_wsp10_wdir000_s1004_struct.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_wsp10_wdir000_s1004_body_eigen.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_wsp10_wdir000_s1004_strc_eigen.dat',
'./res_eigen/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004/mode*.dat',
'./launcher_test/ssystem_eigenanalysis.dat', './launcher_test/mode*.dat',
'./res/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.sel',
'./res/dlc12_iec61400-1ed3/dlc12_wsp10_wdir000_s1004.dat',
'./res/rotor_check_inipos.dat',
'./res/rotor_check_inipos2.dat']:
try:
output_files.remove(f)
except ValueError:
raise ValueError(f + " is not in list")
self.assertFalse(output_files)
def test_turbulence_files(self):
htcfile = HTCFile(self.testfilepath + "dlc14_wsp10_wdir000_s0000.htc", '../')
self.assertEqual(htcfile.turbulence_files(), [
'./turb/turb_wsp10_s0000u.bin', './turb/turb_wsp10_s0000v.bin', './turb/turb_wsp10_s0000w.bin'])
def test_input_files(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
input_files = htcfile.input_files()
#print (htcfile.output)
for f in ['./data/DTU_10MW_RWT_Tower_st.dat',
'./data/DTU_10MW_RWT_Towertop_st.dat',
'./data/DTU_10MW_RWT_Shaft_st.dat',
'./data/DTU_10MW_RWT_Hub_st.dat',
'./data/DTU_10MW_RWT_Blade_st.dat',
'./data/DTU_10MW_RWT_ae.dat',
'./data/DTU_10MW_RWT_pc.dat',
'./control/risoe_controller.dll',
'./control/risoe_controller_64.dll',
'./control/risoe_controller.so',
'./control/generator_servo.dll',
'./control/generator_servo_64.dll',
'./control/generator_servo.so',
'./control/mech_brake.dll',
'./control/mech_brake_64.dll',
'./control/mech_brake.so',
'./control/servo_with_limits.dll',
'./control/servo_with_limits_64.dll',
'./control/servo_with_limits.so',
'./control/towclearsens.dll',
'./control/towclearsens_64.dll',
'./control/towclearsens.so',
'./data/user_shear.dat',
self.testfilepath.replace("\\", "/") + 'test.htc'
]:
try:
input_files.remove(f)
except ValueError:
raise ValueError(f + " is not in list")
self.assertFalse(input_files)
htcfile = HTCFile(self.testfilepath + "DTU_10MW_RWT.htc")
self.assertTrue('./control/wpdata.100' in htcfile.input_files())
def test_input_files2(self):
htcfile = HTCFile(self.testfilepath + "ansi.htc", '../')
input_files = htcfile.input_files()
self.assertTrue('./htc_hydro/ireg_airy_h6_t10.inp' in input_files)
#
def test_continue_in_files(self):
htcfile = HTCFile(self.testfilepath + "continue_in_file.htc", ".")
self.assertIn('main_body__31', htcfile.new_htc_structure.keys())
self.assertIn(os.path.abspath(self.testfilepath + 'orientation.dat'),
[os.path.abspath(f) for f in htcfile.input_files()])
self.assertIn('./data/NREL_5MW_st1.txt', htcfile.input_files())
self.assertEqual(str(htcfile).count("exit"), 1)
self.assertIn('filename\t./res/oc4_p2_load_case_eq;', str(htcfile).lower())
def test_continue_in_files_autodetect_path(self):
htcfile = HTCFile(self.testfilepath + "sub/continue_in_file.htc")
self.assertIn('main_body__31', htcfile.new_htc_structure.keys())
self.assertIn(os.path.abspath(self.testfilepath + 'orientation.dat'),
[os.path.abspath(f) for f in htcfile.input_files()])
self.assertIn('./data/NREL_5MW_st1.txt', htcfile.input_files())
self.assertEqual(str(htcfile).count("exit"), 1)
self.assertIn('filename\t./res/oc4_p2_load_case_eq;', str(htcfile).lower())
def test_tjul_example(self):
htcfile = HTCFile(self.testfilepath + "./tjul.htc", ".")
htcfile.save("./temp.htc")
def test_ansi(self):
htcfile = HTCFile(self.testfilepath + "./ansi.htc", '../')
def test_file_with_BOM(self):
htcfile = HTCFile(self.testfilepath + 'DLC15_wsp11_wdir000_s0000_phi000_Free_v2_visual.htc')
self.assertEqual(str(htcfile)[0], ";")
def test_htc_reset(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
self.assertEqual(htcfile.wind.wsp[0], 10)
def test_htc_model_autodetect(self):
htcfile = HTCFile(self.testfilepath + "test.htc")
self.assertEqual(os.path.relpath(htcfile.modelpath, os.path.dirname(htcfile.filename)), "..")
htcfile = HTCFile(self.testfilepath + "sub/test.htc")
self.assertEqual(os.path.relpath(htcfile.modelpath, os.path.dirname(
htcfile.filename)).replace("\\", "/"), "../..")
self.assertRaisesRegex(ValueError, "Modelpath cannot be autodetected",
HTCFile, self.testfilepath + "missing_input_files.htc")
def test_htc_model_autodetect_upper_case_files(self):
htcfile = HTCFile(self.testfilepath + "../simulation_setup/DTU10MWRef6.0/htc/DTU_10MW_RWT.htc")
self.assertEqual(os.path.relpath(htcfile.modelpath, os.path.dirname(htcfile.filename)), "..")
def test_open_eq_save(self):
HTCFile(self.testfilepath + "test3.htc", "../").save(self.testfilepath + "tmp.htc")
htcfile = HTCFile(self.testfilepath + "tmp.htc", "../")
htcfile.save(self.testfilepath + "tmp.htc")
self.assertEqual(str(htcfile).count("\t"), str(HTCFile(self.testfilepath + "tmp.htc", "../")).count("\t"))
self.assertEqual(str(htcfile), str(HTCFile(self.testfilepath + "tmp.htc", "../")))
def test_2xoutput(self):
htc = HTCFile(self.testfilepath + "test_2xoutput.htc", "../")
self.assertEqual(len(htc.res_file_lst()), 4)
def test_access_section__1(self):
htc = HTCFile(self.testfilepath + "test_2xoutput.htc", "../")
assert htc.output__1.name_ == "output"
def test_compare(self):
htc = HTCFile(self.testfilepath + "test_cmp1.htc", "../")
s = htc.compare(self.testfilepath + 'test_cmp2.htc')
ref = """- begin subsection1;
- end subsection1;
+ begin subsection2;
+ end subsection2;
- begin section1;
- end section1;
- ;
+ alfa 2;
- alfa 1;
+ sensor1 1;
- sensor2 2;
+ begin section2;
+ end section2;
+ ;"""
assert s.strip() == ref
def test_pbs_file(self):
htc = HTCFile(self.testfilepath + "../simulation_setup/DTU10MWRef6.0/htc/DTU_10MW_RWT.htc")
assert os.path.relpath(htc.modelpath, self.testfilepath) == os.path.relpath(
"../simulation_setup/DTU10MWRef6.0/")
from wetb.hawc2.hawc2_pbs_file import JESS_WINE32_HAWC2MB
htc.pbs_file(r"R:\HAWC2_tests\v12.6_mmpe3\hawc2\win32", JESS_WINE32_HAWC2MB)
def test_pbs_file_inout(self):
htc = HTCFile(self.testfilepath + "../simulation_setup/DTU10MWRef6.0_IOS/input/htc/DTU_10MW_RWT.htc")
assert os.path.relpath(htc.modelpath, self.testfilepath) == os.path.relpath(
"../simulation_setup/DTU10MWRef6.0_IOS/input")
from wetb.hawc2.hawc2_pbs_file import JESS_WINE32_HAWC2MB
# print(htc.pbs_file(r"R:\HAWC2_tests\v12.6_mmpe3\hawc2\win32",
# JESS_WINE32_HAWC2MB, input_files=["./input/*"], output_files=['./output/*']))
def test_htc_file_Path_object(self):
from pathlib import Path
htcfile = HTCFile(Path(self.testfilepath) / "test.htc")
def test_htc_copy(self):
htc = HTCFile(self.testfilepath + "test.htc")
tower2 = htc.new_htc_structure.main_body.copy()
tower2.name = "tower2"
htc.new_htc_structure.add_section(tower2, allow_duplicate=True)
assert htc.new_htc_structure.main_body__8 is tower2
assert htc.new_htc_structure.main_body.name[0] == 'tower'
assert htc.new_htc_structure.main_body__8.name[0] == 'tower2'
ti2 = tower2.add_section(section_name='timoschenko_input',
section=tower2.timoschenko_input.copy(), allow_duplicate=True)
ti2.set = 3, 3
assert tower2.timoschenko_input.set.values == [1, 2]
assert tower2.timoschenko_input__2.set.values == [3, 3]
def test_htc_copy2(self):
htc1 = HTCFile(self.testfilepath + "test.htc")
htc2 = htc1.copy()
test1 = len(htc2.compare(htc1)) == 0
test2 = id(htc1) != id(htc2)
assert test1 and test2
assert htc1.modelpath == htc2.modelpath
assert htc2.filename == self.testfilepath.replace("\\", "/") + "test_copy.htc"
assert str(htc1.output) == str(htc2.output)
assert str(htc1.output_at_time) == str(htc2.output_at_time)
def test_location(self):
htc = HTCFile(self.testfilepath + "test.htc")
assert htc.new_htc_structure.main_body__3.location() == 'test.htc/new_htc_structure/main_body__3'
def test__call__(self):
htc = HTCFile(self.testfilepath + "test.htc")
assert htc.new_htc_structure.main_body(name='shaft').name.values[0] == 'shaft'
assert htc.new_htc_structure.main_body.c2_def.sec(v0=3).values[0] == 3
def test_jinja_tags(self):
htc = HTCFile(self.testfilepath + "jinja.htc",
jinja_tags={'wsp': 12, 'logfilename': None, 'begin_step': 100})
assert htc.wind.wsp[0] == 12
assert 'logfile' not in htc.simulation
assert htc.wind.wind_ramp_abs.values == [100, 101, 0, 1]
assert htc.wind.wind_ramp_abs__2.values == [150, 151, 0, 1]
htc = HTCFile(self.testfilepath + "jinja.htc",
jinja_tags={'wsp': 12, 'logfilename': 'test.log', 'begin_step': 100})
assert htc.simulation.logfile[0] == 'test.log'
if __name__ == "__main__":
#import sys;sys.argv = ['', 'Test.testName']
unittest.main()
wetb/hawc2/tests/test_logfile.py
View file @ 62f0f2e4
......@@ -3,13 +3,6 @@ Created on 18/11/2015
@author: MMPE
'''
from __future__ import unicode_literals
from __future__ import print_function
from __future__ import division
from __future__ import absolute_import
from io import open
from future import standard_library
standard_library.install_aliases()
import unittest
from wetb.hawc2.log_file import LogFile, \
INITIALIZATION, SIMULATING, DONE, PENDING, UNKNOWN
......
wetb/hawc2/tests/test_pc_file.py
View file @ 62f0f2e4
'''
Created on 17/07/2014
@author: MMPE
'''
from __future__ import unicode_literals
from __future__ import print_function
from __future__ import division
from __future__ import absolute_import
from future import standard_library
"""Test wetb.hawc2.pc_file
"""
from wetb.hawc2.ae_file import AEFile
standard_library.install_aliases()
import os
import unittest
from wetb.hawc2.pc_file import PCFile
import tempfile
import numpy as np
from wetb.hawc2.pc_file import PCFile
class TestPCFile(unittest.TestCase):
......@@ -26,6 +18,7 @@ class TestPCFile(unittest.TestCase):
self.testfilepath = os.path.join(os.path.dirname(__file__), 'test_files/') # test file path
def test_PCFile_ae(self):
"""Verify correct values in values of object loaded from pc file"""
pc = PCFile(self.testfilepath + "NREL_5MW_pc.txt")
ae = AEFile(self.testfilepath + "NREL_5MW_ae.txt")
thickness = ae.thickness(36)
......@@ -33,6 +26,22 @@ class TestPCFile(unittest.TestCase):
self.assertEqual(pc.CD(thickness, 10), 0.0255)
self.assertEqual(pc.CM(thickness, 10), -0.1103)
def test_write_PCFile(self):
"""Round-trip loading and saving a pc file
"""
pc1 = PCFile(self.testfilepath + 'NREL_5MW_pc.txt')
with tempfile.TemporaryDirectory() as tdir:
pc1.save(tdir + '/test_pc.txt')
pc2 = PCFile(pc1.filename)
self.assertEqual(str(pc1), str(pc2))
tc1, pcs1 = pc1.pc_sets[1]
tc2, pcs2 = pc2.pc_sets[1]
np.testing.assert_array_almost_equal(tc1, tc2)
for pc1, pc2 in zip(pcs1, pcs2):
np.testing.assert_array_almost_equal(pc1, pc2)
if __name__ == "__main__":
#import sys;sys.argv = ['', 'Test.testName']
unittest.main()
wetb/hawc2/tests/test_sel_file.py
View file @ 62f0f2e4
......@@ -3,12 +3,6 @@ Created on 17/07/2014
@author: MMPE
'''
from __future__ import unicode_literals
from __future__ import print_function
from __future__ import division
from __future__ import absolute_import
from future import standard_library
standard_library.install_aliases()
import unittest
from wetb.hawc2.sel_file import SelFile, BINARY, ASCII
from datetime import datetime
......
wetb/hawc2/tests/test_sensor_names.py 0 → 100644
View file @ 62f0f2e4
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Mon Mar 10 11:11:48 2025
@author: dave
"""
import unittest
import numpy as np
import pandas as pd
from wetb.hawc2.sensor_names import unified_channel_names
from wetb.hawc2.Hawc2io import ReadHawc2
import os
testfilepath = os.path.join(os.path.dirname(__file__), 'test_files/hawc2io/') # test file path
class TestHAWC2SensorNames(unittest.TestCase):
def test_unified_channel_names(self):
# only read the sel file (not the actual data)
fname = os.path.join(testfilepath, 'hawc2bin_chantest_3.sel')
res = ReadHawc2(fname)
ch_dict, ch_df = unified_channel_names(res.ChInfo)
# when loading the reference result from csv everything will be a string
# so this makes the comparison easier. What is missing in unified_channel_names
# is also casting each column to the correct data types.
ch_df = ch_df.astype(str)
ch_df_ref = pd.read_csv(fname.replace('.sel', '.ch_df.csv'), sep=',',
dtype=str, keep_default_na=False, index_col=0)
np.testing.assert_equal(ch_df.columns.values, ch_df_ref.columns.values)
np.testing.assert_equal(ch_df.values, ch_df_ref.values)
if __name__ == "__main__":
unittest.main()
wetb/hawc2/tests/test_shear_file.py
View file @ 62f0f2e4
......@@ -3,14 +3,7 @@ Created on 05/11/2015
@author: MMPE
'''
from __future__ import unicode_literals
from __future__ import print_function
from __future__ import division
from __future__ import absolute_import
from io import open
from future import standard_library
from wetb.hawc2.shear_file import ShearFile
standard_library.install_aliases()
import unittest
from wetb.hawc2 import shear_file
import numpy as np
......
wetb/hawc2/tests/test_st_file.py
View file @ 62f0f2e4
......@@ -3,34 +3,65 @@ Created on 05/11/2015
@author: MMPE
'''
from __future__ import unicode_literals
from __future__ import print_function
from __future__ import division
from __future__ import absolute_import
from future import standard_library
standard_library.install_aliases()
import unittest
from wetb.hawc2.st_file import StFile
import os
from numpy import testing
from wetb.hawc2.st_file import StFile
testfilepath = os.path.join(os.path.dirname(__file__), 'test_files/') # test file path
class TestStFile(unittest.TestCase):
def test_stfile(self):
st = StFile(testfilepath + "DTU_10MW_RWT_Blade_st.dat")
st = StFile(testfilepath + 'DTU_10MW_RWT_Blade_st.dat')
self.assertEqual(st.radius_st()[2], 3.74238)
self.assertEqual(st.radius_st(3), 3.74238)
self.assertEqual(st.x_e(67.7351), 4.4320990737400E-01)
self.assertEqual(st.E(3.74238, 1, 1), 1.2511695058500E+10)
self.assertEqual(st.E(3.74238, 1, 2), 1.2511695058500E+27)
def test_fpmfile(self):
st = StFile(os.path.join(testfilepath, 'IEA_15MW_RWT_Blade_fpm_st.dat'))
self.assertEqual(st.radius_st()[2], 5.857445433495505)
self.assertEqual(st.radius_st(10), 11.71489728616879)
self.assertEqual(st.x_e(67.7351), 0.7150636117520498)
self.assertEqual(st.K_11(3.5, 1, 1), 5572685723.993099)
self.assertEqual(st.K_11(3.5, 1, 2), 5572685723.993099)
def test_stfile_interpolate(self):
st = StFile(testfilepath + "DTU_10MW_RWT_Blade_st.dat")
st = StFile(testfilepath + 'DTU_10MW_RWT_Blade_st.dat')
self.assertAlmostEqual(st.x_e(72.2261), 0.381148048)
self.assertAlmostEqual(st.y_e(72.2261), 0.016692967)
def test_save(self):
fname = os.path.join(testfilepath, 'DTU_10MW_RWT_Blade_st.dat')
fname2 = os.path.join(testfilepath, 'DTU_10MW_RWT_Blade_st2.dat')
st = StFile(fname)
st.save(fname2, encoding='utf-8', precision='%20.12e')
st2 = StFile(fname2)
self.assertEqual(len(st.main_data_sets), len(st2.main_data_sets))
self.assertEqual(len(st.main_data_sets[1]), len(st2.main_data_sets[1]))
for k in st.main_data_sets[1]:
testing.assert_almost_equal(st.main_data_sets[1][k],
st2.main_data_sets[1][k], decimal=12)
os.remove(fname2)
def test_save_fpm(self):
fname = os.path.join(testfilepath, 'IEA_15MW_RWT_Blade_fpm_st.dat')
fname2 = os.path.join(testfilepath, 'IEA_15MW_RWT_Blade_fpm_st2.dat')
st = StFile(fname)
st.save(fname2, encoding='utf-8', precision='% 24.15e')
st2 = StFile(fname2)
self.assertEqual(len(st.main_data_sets), len(st2.main_data_sets))
self.assertEqual(len(st.main_data_sets[1]), len(st2.main_data_sets[1]))
for k in st.main_data_sets[1]:
testing.assert_almost_equal(st.main_data_sets[1][k],
st2.main_data_sets[1][k], decimal=12)
os.remove(fname2)
if __name__ == "__main__":
#import sys;sys.argv = ['', 'Test.testName']
......