Hello all,
As i know the axial induction factor is determined from the following equation at aerosubs/ SUBROUTINE vindinf:
A(iRadius,iBlade) = A(iRadius,iBlade) + xphi(Rzero,mode) * ( xAlpha(mode) * COS( REAL(MRvector(MODE)) * Windpsi ) + xBeta (mode) * SIN( REAL(MRvector(MODE)) * Windpsi ) )
where windpsi is the azimuth angle.As i understand this azimuth angle is determined from aerodyn using the following relation:
AzimuthAngle = ATAN2( -1.*DOT_PRODUCT( TurbineComponents%Hub%Orientation(3,:), &
TurbineComponents%RotorFurl%Orientation(2, ), &
DOT_PRODUCT( TurbineComponents%Hub%Orientation(3,:), &
TurbineComponents%RotorFurl%Orientation(3, ) ) + pi + (IBlade - 1)*TwoPiNB
while TurbineComponents%Hub%Orientation & TurbineComponents%RotorFurl%Orientation, can be determined from Fast by the following relation:
ADInterfaceComponents%Hub%Orientation(:,1) = (/ e1(1), e2(1), e3(1) /)
ADInterfaceComponents%Hub%Orientation(:,2) = -1.*(/ e1(3), e2(3), e3(3) /)
ADInterfaceComponents%Hub%Orientation(:,3) = (/ e1(2), e2(2), e3(2) /)
ADInterfaceComponents%RotorFurl%Orientation(:,1) = (/ c1(1), -1.c3(1), c2(1) /)
ADInterfaceComponents%RotorFurl%Orientation(:,2) = (/ -1. c1(3), c3(3), -1.*c2(3) /)
ADInterfaceComponents%RotorFurl%Orientation(:,3) = (/ c1(2), -1.*c3(2), c2(2) /)
this point makes me confused. so, is there any schematic that shows the orientation of the vectors e ( ) and c( )?.also , what is QT(DOF-DRTR) & QT(DOF-GEAZ) that are used to determine e2 & e3?.
So finally, i would like to ask how to calculate the azimuth angle that is used in GDW?, and if i assumed a simple rotor with fixed blades and rigid hub, could i calculate the azimuth angle from the simple following equation?
aerodynamic torque-generator torque= j d(rotor speed)/dt
Best regards
Hammam