diff --git a/wetb/fatigue_tools/bearing_damage.py b/wetb/fatigue_tools/bearing_damage.py index 145def9a58727e151e552cbf1b953eca663ea6b8..62b99d4906d1cf2a12afad9a8801b90d932a61de 100644 --- a/wetb/fatigue_tools/bearing_damage.py +++ b/wetb/fatigue_tools/bearing_damage.py @@ -4,35 +4,9 @@ Created on 13/10/2014 @author: MMPE ''' -def bearing_damage_SWP(angle_moment_lst, m=3, thresshold=0.1): +def bearing_damage(angle_moment_lst, m=3, thresshold=0.1): """Function ported from Matlab. - - Matlab documentation: - "Compute bearing damage according to SWP procedure - The basic equation is a damage equivalent loading calculated as - - Damage = sum(dPitch*Mr^3.0) - - Where - dPitch = pitch change in a given timestep [deg] - Mr = load component in a given sector [kNm] - - We evaluate Mr (blade root moment resultant) in 10deg sectors and take - the largest damage. For a delta loading as you are investigating it is - accurate enough just to take the flap loading and ignore the changes - coming from edge load variation. - - The delta pitch is a bit tricky. You need to filter you signal so minor - quick variation in pitch does not give an actual delta pitch. The - variation in pitch reference is filtered. So effectively a variation - below 0.1deg is ignored and the pitch is fixed with no variation as - long as the change is below this value. - - So basically for each time step calulate dPitch*Mr^3 and sum it up with - Mr=blade flap moment. I’m sure your model will not have the details for - capacity comparison, but for changes from one configuration to another - should be ok. We do not use the sum(dPitch) for design evaluation directly." - + Parameters ---------- angle_moment_lst : ((angle_A_vector, moment_A_vector),(angle_B_vector, moment_B_vector),...)