diff --git a/wetb/fatigue_tools/tests/test_bearing_damage.py b/wetb/fatigue_tools/tests/test_bearing_damage.py
index 37e33b48af2a1bfc2535ee832822938d5cf98fa9..ae7020864314f3f3fa43a921f915938acf7827c7 100644
--- a/wetb/fatigue_tools/tests/test_bearing_damage.py
+++ b/wetb/fatigue_tools/tests/test_bearing_damage.py
@@ -8,6 +8,7 @@ from __future__ import print_function
from __future__ import division
from __future__ import absolute_import
from future import standard_library
+from wetb import gtsdf
standard_library.install_aliases()
@@ -26,8 +27,8 @@ class TestBearingDamage(unittest.TestCase):
self.tfp = os.path.join(os.path.dirname(__file__), 'test_files/')
def test_bearing_damage_swp(self):
- data = Hawc2io.ReadHawc2(self.tfp + "test_bearing_damage").ReadBinary((np.array([4, 26, 6, 32, 8, 38]) - 1).tolist())
- self.assertAlmostEqual(bearing_damage([(data[:, i], data[:, i + 1]) for i in [0, 2, 4]]), 7.755595081475002e+13)
+ data = Hawc2io.ReadHawc2(self.tfp + "test_bearing_damage").ReadBinary((np.array([1, 4, 2, 5, 3, 6]) ).tolist())
+ self.assertAlmostEqual(bearing_damage([(data[:, i], data[:, i + 1]) for i in [0,2,4]]), 7.755595081475002e+13)
if __name__ == "__main__":
diff --git a/wetb/fatigue_tools/tests/test_files/test_bearing_damage.dat b/wetb/fatigue_tools/tests/test_files/test_bearing_damage.dat
index 7389fa219a0f61778b4108dacf54c78f8d428a30..3162eda66e5f18393355fe716c7de4d015a8c616 100644
Binary files a/wetb/fatigue_tools/tests/test_files/test_bearing_damage.dat and b/wetb/fatigue_tools/tests/test_files/test_bearing_damage.dat differ
diff --git a/wetb/fatigue_tools/tests/test_files/test_bearing_damage.sel b/wetb/fatigue_tools/tests/test_files/test_bearing_damage.sel
index 064eeae15931b1636f04d3677bc8982ad46caa2d..c57f6facaf9851adbd113030342e65f9632b514c 100644
--- a/wetb/fatigue_tools/tests/test_files/test_bearing_damage.sel
+++ b/wetb/fatigue_tools/tests/test_files/test_bearing_damage.sel
@@ -1,328 +1,28 @@
________________________________________________________________________________________________________________________
- Version ID : HAWC2MB 11.8
- Time : 00:05:42
- Date : 11:10.2014
+ Version ID : Pdap 3.12
+ Time : 20:36:25
+ Date : 07:12.2016
________________________________________________________________________________________________________________________
- Result file : ./res/dlc12_iec61400-1ed3/dlc12_wsp12_wdir000_s1005.dat
+ Result file : D:/Dokumenter/Programmering/Python/workspace/WindEnergyToolbox/wetb/fatigue_tools/tests/test_files/test_bearing_damage2.dat
________________________________________________________________________________________________________________________
Scans Channels Time [sec] Format
- 30000 157 600.000 BINARY
+ 30000 7 600.000 BINARY
Channel Variable Description
1 Time s Time
- 2 bea1 angle deg shaft_rot angle
- 3 bea1 angle_speed rpm shaft_rot angle speed
- 4 bea1 angle deg pitch1 angle
- 5 bea1 angle_speed deg/s pitch1 angle speed
- 6 bea1 angle deg pitch2 angle
- 7 bea1 angle_speed deg/s pitch2 angle speed
- 8 bea1 angle deg pitch3 angle
- 9 bea1 angle_speed deg/s pitch3 angle speed
- 10 Omega rad/s Rotor speed
- 11 Ae rot. torque kNm Aero rotor torque
- 12 Ae rot. power kW Aero rotor power
- 13 Ae rot. thrust kN Aero rotor thrust
- 14 WSP gl. coo.,Vx m/s Free wind speed Vx, gl. coo, of gl. pos 0.00, 0.00, -90.00
- 15 WSP gl. coo.,Vy m/s Free wind speed Vy, gl. coo, of gl. pos 0.00, 0.00, -90.00
- 16 WSP gl. coo.,Vz m/s Free wind speed Vz, gl. coo, of gl. pos 0.00, 0.00, -90.00
- 17 Mx coo: tower kNm MomentMx Mbdy:tower nodenr: 1 coo: tower tower base
- 18 My coo: tower kNm MomentMy Mbdy:tower nodenr: 1 coo: tower tower base
- 19 Mz coo: tower kNm MomentMz Mbdy:tower nodenr: 1 coo: tower tower base
- 20 Mx coo: tower kNm MomentMx Mbdy:tower nodenr: 8 coo: tower yaw bearing
- 21 My coo: tower kNm MomentMy Mbdy:tower nodenr: 8 coo: tower yaw bearing
- 22 Mz coo: tower kNm MomentMz Mbdy:tower nodenr: 8 coo: tower yaw bearing
- 23 Mx coo: shaft kNm MomentMx Mbdy:shaft nodenr: 4 coo: shaft main bearing
- 24 My coo: shaft kNm MomentMy Mbdy:shaft nodenr: 4 coo: shaft main bearing
- 25 Mz coo: shaft kNm MomentMz Mbdy:shaft nodenr: 4 coo: shaft main bearing
- 26 Mx coo: blade1 kNm MomentMx Mbdy:blade1 nodenr: 3 coo: blade1 blade 1 root
- 27 My coo: blade1 kNm MomentMy Mbdy:blade1 nodenr: 3 coo: blade1 blade 1 root
- 28 Mz coo: blade1 kNm MomentMz Mbdy:blade1 nodenr: 3 coo: blade1 blade 1 root
- 29 Mx coo: local kNm MomentMx Mbdy:blade1 nodenr: 10 coo: local blade 1 50% local e coo
- 30 My coo: local kNm MomentMy Mbdy:blade1 nodenr: 10 coo: local blade 1 50% local e coo
- 31 Mz coo: local kNm MomentMz Mbdy:blade1 nodenr: 10 coo: local blade 1 50% local e coo
- 32 Mx coo: blade2 kNm MomentMx Mbdy:blade2 nodenr: 3 coo: blade2 blade 2 root
- 33 My coo: blade2 kNm MomentMy Mbdy:blade2 nodenr: 3 coo: blade2 blade 2 root
- 34 Mz coo: blade2 kNm MomentMz Mbdy:blade2 nodenr: 3 coo: blade2 blade 2 root
- 35 Mx coo: local kNm MomentMx Mbdy:blade2 nodenr: 10 coo: local blade 2 50% local e coo
- 36 My coo: local kNm MomentMy Mbdy:blade2 nodenr: 10 coo: local blade 2 50% local e coo
- 37 Mz coo: local kNm MomentMz Mbdy:blade2 nodenr: 10 coo: local blade 2 50% local e coo
- 38 Mx coo: blade3 kNm MomentMx Mbdy:blade3 nodenr: 3 coo: blade3 blade 3 root
- 39 My coo: blade3 kNm MomentMy Mbdy:blade3 nodenr: 3 coo: blade3 blade 3 root
- 40 Mz coo: blade3 kNm MomentMz Mbdy:blade3 nodenr: 3 coo: blade3 blade 3 root
- 41 Mx coo: local kNm MomentMx Mbdy:blade3 nodenr: 10 coo: local blade 3 50% local e coo
- 42 My coo: local kNm MomentMy Mbdy:blade3 nodenr: 10 coo: local blade 3 50% local e coo
- 43 Mz coo: local kNm MomentMz Mbdy:blade3 nodenr: 10 coo: local blade 3 50% local e coo
- 44 State pos x coo: global m State pos x Mbdy:tower E-nr: 7 Z-rel:1.00 coo: global tower top flange position
- 45 State pos y coo: global m State pos y Mbdy:tower E-nr: 7 Z-rel:1.00 coo: global tower top flange position
- 46 State pos z coo: global m State pos z Mbdy:tower E-nr: 7 Z-rel:1.00 coo: global tower top flange position
- 47 State acc x coo: global m/s^2 State acc x Mbdy:tower E-nr: 7 Z-rel:1.00 coo: global tower top flange position
- 48 State acc y coo: global m/s^2 State acc y Mbdy:tower E-nr: 7 Z-rel:1.00 coo: global tower top flange position
- 49 State acc z coo: global m/s^2 State acc z Mbdy:tower E-nr: 7 Z-rel:1.00 coo: global tower top flange position
- 50 State pos x coo: global m State pos x Mbdy:tower E-nr: 3 Z-rel:1.00 coo: global tower approx tip height
- 51 State pos y coo: global m State pos y Mbdy:tower E-nr: 3 Z-rel:1.00 coo: global tower approx tip height
- 52 State pos z coo: global m State pos z Mbdy:tower E-nr: 3 Z-rel:1.00 coo: global tower approx tip height
- 53 State pos x coo: blade1 m State pos x Mbdy:blade1 E-nr: 18 Z-rel:1.00 coo: blade1 blade 1 tip pos
- 54 State pos y coo: blade1 m State pos y Mbdy:blade1 E-nr: 18 Z-rel:1.00 coo: blade1 blade 1 tip pos
- 55 State pos z coo: blade1 m State pos z Mbdy:blade1 E-nr: 18 Z-rel:1.00 coo: blade1 blade 1 tip pos
- 56 State pos x coo: blade2 m State pos x Mbdy:blade2 E-nr: 18 Z-rel:1.00 coo: blade2 blade 2 tip pos
- 57 State pos y coo: blade2 m State pos y Mbdy:blade2 E-nr: 18 Z-rel:1.00 coo: blade2 blade 2 tip pos
- 58 State pos z coo: blade2 m State pos z Mbdy:blade2 E-nr: 18 Z-rel:1.00 coo: blade2 blade 2 tip pos
- 59 State pos x coo: blade3 m State pos x Mbdy:blade3 E-nr: 18 Z-rel:1.00 coo: blade3 blade 3 tip pos
- 60 State pos y coo: blade3 m State pos y Mbdy:blade3 E-nr: 18 Z-rel:1.00 coo: blade3 blade 3 tip pos
- 61 State pos z coo: blade3 m State pos z Mbdy:blade3 E-nr: 18 Z-rel:1.00 coo: blade3 blade 3 tip pos
- 62 State pos x coo: global m State pos x Mbdy:blade1 E-nr: 18 Z-rel:1.00 coo: global blade 1 tip pos
- 63 State pos y coo: global m State pos y Mbdy:blade1 E-nr: 18 Z-rel:1.00 coo: global blade 1 tip pos
- 64 State pos z coo: global m State pos z Mbdy:blade1 E-nr: 18 Z-rel:1.00 coo: global blade 1 tip pos
- 65 State pos x coo: global m State pos x Mbdy:blade2 E-nr: 18 Z-rel:1.00 coo: global blade 2 tip pos
- 66 State pos y coo: global m State pos y Mbdy:blade2 E-nr: 18 Z-rel:1.00 coo: global blade 2 tip pos
- 67 State pos z coo: global m State pos z Mbdy:blade2 E-nr: 18 Z-rel:1.00 coo: global blade 2 tip pos
- 68 State pos x coo: global m State pos x Mbdy:blade3 E-nr: 18 Z-rel:1.00 coo: global blade 3 tip pos
- 69 State pos y coo: global m State pos y Mbdy:blade3 E-nr: 18 Z-rel:1.00 coo: global blade 3 tip pos
- 70 State pos z coo: global m State pos z Mbdy:blade3 E-nr: 18 Z-rel:1.00 coo: global blade 3 tip pos
- 71 WSP Vx, glco, R= 61.5 m/s Wind speed Vx of blade 1 at radius 61.52, global coo.
- 72 WSP Vy, glco, R= 61.5 m/s Wind speed Vy of blade 1 at radius 61.52, global coo.
- 73 WSP Vz, glco, R= 61.5 m/s Wind speed Vz of blade 1 at radius 61.52, global coo.
- 74 Alfa, R= 45.0 deg Angle of attack of blade 1 at radius 45.17
- 75 Alfa, R= 45.0 deg Angle of attack of blade 2 at radius 45.17
- 76 Alfa, R= 45.0 deg Angle of attack of blade 3 at radius 45.17
- 77 Cl, R= 45.0 deg Cl of blade 1 at radius 45.17
- 78 Cl, R= 45.0 deg Cl of blade 2 at radius 45.17
- 79 Cl, R= 45.0 deg Cl of blade 3 at radius 45.17
- 80 DLL out 1: 1 - DLL : 1 outvec : 1 time
- 81 DLL out 1: 2 - DLL : 1 outvec : 2 slow speed shaft rad/s
- 82 DLL out 1: 3 - DLL : 1 outvec : 3 pitch angle 1
- 83 DLL out 1: 4 - DLL : 1 outvec : 4 pitch angle 2
- 84 DLL out 1: 5 - DLL : 1 outvec : 5 pitch angle 3
- 85 DLL out 1: 6 - DLL : 1 outvec : 6 wsp_x_global
- 86 DLL out 1: 7 - DLL : 1 outvec : 7 wsp_y_global
- 87 DLL out 1: 8 - DLL : 1 outvec : 8 wsp_z_global
- 88 DLL out 1: 9 - DLL : 1 outvec : 9 elec. pwr
- 89 DLL out 1: 10 - DLL : 1 outvec : 10 grid flag
- 90 DLL inp 1: 1 - DLL : 1 inpvec : 1 generator torque reference [nm]
- 91 DLL inp 1: 2 - DLL : 1 inpvec : 2 pitch angle reference of blade 1 [rad]
- 92 DLL inp 1: 3 - DLL : 1 inpvec : 3 pitch angle reference of blade 2 [rad]
- 93 DLL inp 1: 4 - DLL : 1 inpvec : 4 pitch angle reference of blade 3 [rad]
- 94 DLL inp 1: 5 - DLL : 1 inpvec : 5 power reference [w]
- 95 DLL inp 1: 6 - DLL : 1 inpvec : 6 filtered wind speed [m/s]
- 96 DLL inp 1: 7 - DLL : 1 inpvec : 7 filtered rotor speed [rad/s]
- 97 DLL inp 1: 8 - DLL : 1 inpvec : 8 filtered rotor speed error for torque [rad/s]
- 98 DLL inp 1: 9 - DLL : 1 inpvec : 9 bandpass filtered rotor speed [rad/s]
- 99 DLL inp 1: 10 - DLL : 1 inpvec : 10 proportional term of torque contr. [nm]
- 100 DLL inp 1: 11 - DLL : 1 inpvec : 11 integral term of torque controller [nm]
- 101 DLL inp 1: 12 - DLL : 1 inpvec : 12 minimum limit of torque [nm]
- 102 DLL inp 1: 13 - DLL : 1 inpvec : 13 maximum limit of torque [nm]
- 103 DLL inp 1: 14 - DLL : 1 inpvec : 14 torque limit switch based on pitch [-]
- 104 DLL inp 1: 15 - DLL : 1 inpvec : 15 filtered rotor speed error for pitch [rad/s]
- 105 DLL inp 1: 16 - DLL : 1 inpvec : 16 power error for pitch [w]
- 106 DLL inp 1: 17 - DLL : 1 inpvec : 17 proportional term of pitch controller [rad]
- 107 DLL inp 1: 18 - DLL : 1 inpvec : 18 integral term of pitch controller [rad]
- 108 DLL inp 1: 19 - DLL : 1 inpvec : 19 minimum limit of pitch [rad]
- 109 DLL inp 1: 20 - DLL : 1 inpvec : 20 maximum limit of pitch [rad]
- 110 DLL inp 1: 21 - DLL : 1 inpvec : 21 torque reference from dt dammper [nm]
- 111 DLL inp 1: 22 - DLL : 1 inpvec : 22 status signal [-]
- 112 DLL inp 1: 23 - DLL : 1 inpvec : 23 total added pitch rate [rad/s]
- 113 DLL inp 1: 25 - DLL : 1 inpvec : 25 flag for mechnical brake [0=off/1=on]
- 114 DLL inp 1: 26 - DLL : 1 inpvec : 26 flag for emergency pitch stop [0=off/1=on]
- 115 DLL out 2: 1 - DLL : 2 outvec : 1 time
- 116 DLL out 2: 2 - DLL : 2 outvec : 2 azimuth angle [rad]
- 117 DLL out 2: 3 - DLL : 2 outvec : 3 azimuth speed [rad/s]
- 118 DLL out 2: 4 - DLL : 2 outvec : 4 flap brm of blade 1 [knm]
- 119 DLL out 2: 5 - DLL : 2 outvec : 5 flap brm of blade 2 [knm]
- 120 DLL out 2: 6 - DLL : 2 outvec : 6 flap brm of blade 3 [knm]
- 121 DLL out 2: 7 - DLL : 2 outvec : 7 pitch reference of blade 1 [rad]
- 122 DLL out 2: 8 - DLL : 2 outvec : 8 pitch reference of blade 2 [rad]
- 123 DLL out 2: 9 - DLL : 2 outvec : 9 pitch reference of blade 3 [rad]
- 124 DLL out 2: 10 - DLL : 2 outvec : 10 status flag from basic ctrl
- 125 DLL inp 2: 1 - DLL : 2 inpvec : 1 pitch reference of blade 1 [rad]
- 126 DLL inp 2: 2 - DLL : 2 inpvec : 2 pitch reference of blade 2 [rad]
- 127 DLL inp 2: 3 - DLL : 2 inpvec : 3 pitch reference of blade 3 [rad]
- 128 DLL inp 2: 4 - DLL : 2 inpvec : 4 cosine moment [knm]
- 129 DLL inp 2: 5 - DLL : 2 inpvec : 5 sine moment [knm]
- 130 DLL inp 2: 6 - DLL : 2 inpvec : 6 low-pass filtered cosine moment [knm]
- 131 DLL inp 2: 7 - DLL : 2 inpvec : 7 low-pass filtered sine moment [knm]
- 132 DLL inp 2: 8 - DLL : 2 inpvec : 8 cosine pitch [rad]
- 133 DLL inp 2: 9 - DLL : 2 inpvec : 9 sine pitch [rad]
- 134 DLL inp 2: 10 - DLL : 2 inpvec : 10 switchfactor [-]
- 135 DLL out 3: 1 - DLL : 3 outvec : 1 time
- 136 DLL out 3: 2 - DLL : 3 outvec : 2 electrical torque reference [nm]
- 137 DLL out 3: 3 - DLL : 3 outvec : 3 omega lss
- 138 DLL inp 3: 1 - DLL : 3 inpvec : 1 mgen lss [nm]
- 139 DLL inp 3: 2 - DLL : 3 inpvec : 2 pelec w
- 140 DLL inp 3: 3 - DLL : 3 inpvec : 3 mframe
- 141 DLL inp 3: 4 - DLL : 3 inpvec : 4 mgen hss
- 142 DLL inp 3: 5 - DLL : 3 inpvec : 5 generator pmech kw
- 143 DLL inp 3: 6 - DLL : 3 inpvec : 6 filtered gen speed
- 144 DLL inp 3: 7 - DLL : 3 inpvec : 7 elec. pwr
- 145 DLL inp 3: 8 - DLL : 3 inpvec : 8 grid flag
- 146 DLL inp 4: 1 - DLL : 4 inpvec : 1 brake torque [nm]
- 147 DLL out 4: 1 - DLL : 4 outvec : 1 time [s]
- 148 DLL out 4: 2 - DLL : 4 outvec : 2 generator lss speed [rad/s]
- 149 DLL out 4: 3 - DLL : 4 outvec : 3 deploy brake
- 150 DLL out 5: 1 - DLL : 5 outvec : 1 time
- 151 DLL out 5: 2 - DLL : 5 outvec : 2 pitchref 1
- 152 DLL out 5: 3 - DLL : 5 outvec : 3 pitchref 2
- 153 DLL out 5: 4 - DLL : 5 outvec : 4 pitchref 3
- 154 DLL out 5: 5 - DLL : 5 outvec : 5 emerg. stop
- 155 DLL inp 5: 1 - DLL : 5 inpvec : 1 pitch 1
- 156 DLL inp 5: 2 - DLL : 5 inpvec : 2 pitch 2
- 157 DLL inp 5: 3 - DLL : 5 inpvec : 3 pitch 3
+ 2 bea1 angle deg pitch1 angle
+ 3 bea1 angle deg pitch2 angle
+ 4 bea1 angle deg pitch3 angle
+ 5 Mx coo: blade1 kNm MomentMx Mbdy:blade1 nodenr: 3 coo: blade1 blade 1 root
+ 6 Mx coo: blade2 kNm MomentMx Mbdy:blade2 nodenr: 3 coo: blade2 blade 2 root
+ 7 Mx coo: blade3 kNm MomentMx Mbdy:blade3 nodenr: 3 coo: blade3 blade 3 root
________________________________________________________________________________________________________________________
Scale factors:
2.18750E-02
- 1.12495E-02
- 4.15512E-04
4.09937E-04
- 1.38743E-04
4.04169E-04
- 1.32760E-04
4.03384E-04
- 1.35139E-04
- 4.34542E-05
- 2.02439E-01
- 2.61031E-01
- 2.65983E-02
- 1.61135E-04
- 6.11529E-04
- 1.42944E-04
- 2.42531E+00
- 3.72165E-01
- 2.84396E-01
- 2.68370E-01
- 1.56493E-01
- 2.83659E-01
- 2.35745E-01
- 2.80691E-01
- 1.58515E-01
3.92345E-01
- 1.61766E-01
- 4.37487E-03
- 1.24812E-01
- 2.80450E-02
- 2.22626E-03
3.96563E-01
- 1.65199E-01
- 3.95572E-03
- 1.23898E-01
- 2.65236E-02
- 2.18839E-03
- 3.76100E-01
- 1.51924E-01
- 4.14227E-03
- 1.14018E-01
- 2.71120E-02
- 2.38153E-03
- 2.83396E-06
- 1.52113E-05
- 2.73736E-03
- 1.03974E-05
- 2.85995E-05
- 1.63177E-06
- 2.95506E-07
- 1.86791E-06
- 9.37454E-04
- 4.76602E-05
- 2.66261E-04
- 1.92210E-03
- 4.95472E-05
- 2.61835E-04
- 1.92207E-03
- 4.64387E-05
- 2.42462E-04
- 1.92216E-03
- 1.97098E-03
- 3.75766E-04
- 4.77808E-03
- 1.96979E-03
- 3.65905E-04
- 4.77677E-03
- 1.97069E-03
- 4.01181E-04
- 4.77924E-03
- 2.09823E-04
- 7.06817E-04
- 1.59703E-04
- 3.40315E-04
- 3.51060E-04
- 2.96693E-04
- 4.55834E-05
- 4.69429E-05
- 4.27909E-05
- 2.18750E-02
- 4.35124E-05
- 7.15475E-06
- 7.05409E-06
- 7.04039E-06
- 1.61135E-04
- 6.11529E-04
- 1.42944E-04
- 1.56661E+02
- 1.00000E+00
- 1.33565E+02
- 6.23526E-06
- 6.23526E-06
- 6.23526E-06
- 1.56268E+02
- 4.92568E-04
- 4.33530E-05
- 5.71198E-06
- 1.00000E+00
- 1.10332E+02
- 1.89505E+02
- 1.33565E+02
- 1.33565E+02
- 3.12500E-05
- 5.71198E-06
- 6.98173E+01
- 4.51520E-06
- 5.88678E-06
- 1.00000E+00
- 4.58149E-05
- 1.00000E+00
- 1.00000E+00
- 1.00000E+00
- 1.00000E+00
- 1.00000E+00
- 2.18750E-02
- 9.81722E-05
- 4.35124E-05
- 3.92970E-01
- 3.97220E-01
- 3.76856E-01
- 6.23526E-06
- 6.23526E-06
- 6.23526E-06
- 1.00000E+00
- 7.15371E-06
- 7.05734E-06
- 7.03710E-06
- 1.58750E-01
- 1.70334E-01
- 9.60481E-02
- 9.77392E-02
- 2.18166E-06
- 2.18166E-06
- 3.12344E-05
- 2.18750E-02
- 1.33565E+02
- 4.35124E-05
- 1.41487E+02
- 1.56661E+02
- 1.40028E+02
- 1.45863E+00
- 1.65955E-01
- 1.58515E+02
- 2.95000E-05
- 1.00000E+00
- 1.00000E+00
- 2.18750E-02
- 4.35124E-05
- 1.00000E+00
- 2.18750E-02
- 7.15371E-06
- 7.05734E-06
- 7.03710E-06
- 1.00000E+00
- 7.15475E-06
- 7.05409E-06
- 7.04039E-06
+ 3.76100E-01
\ No newline at end of file