4.2.7.3. ElastoDyn Theory
Note this document is work in progress and is greatly incomplete. This documentation was started to document some code changes to the the tail furl and rotor furl part of ElastoDyn. Please refer to the different ressources provided in Section 4.2.7 for additional documents.
4.2.7.3.1. Notations
Points
The following (partial) list of points are defined by ElastoDyn:
Z
: the platform reference pointO
: the tower-top/base plate pointW
: the specified point on the tail-furl axisI
: the tail boom center of massJ
: the tail fin center of mass
Bodies
The following (partial) list of bodies are defined by ElastoDyn:
E
: the earth/inertial frameX
: the platform bodyN
: the nacelle bodyA
: the tail-furl body
4.2.7.3.2. Kinematics
ElastoDyn computes the position, velocity and accelerations of key points of the structure, starting from the platform reference point Z and going up in the structure.
The different position vectors are available in the data stucture RtHSdat
.
For instance, the global position of point J is given by:
The translational displacement vector (how much a point has moved compared to its reference position) is calculated as follows: \(\boldsymbol{r}_J-\boldsymbol{r}_{J,\text{ref}}\).
The coordinate systems of ElastoDyn are stored in the variable CoordSys
.
The orientation matrix of a given coordinate system can be formed using the unit vectors (assumed to be column vectors) of a given coordinate system expressed in the inertial frame.
For instance for the tailfin coordinate system:
Angular velocities are stored in variables RtHSdat%AngVelE*
with respect to the initial frame (“Earth”, E).
For instance, the angular velocity of the tail-furl body (body A) is:
where \(\boldsymbol{\omega}_{N/X}=\boldsymbol{\omega}_{B/X}+\boldsymbol{\omega}_{N/B}\)
Linear (translational) velocities of the different points are found in the variables RtHSdat%LinVelE*
, and are computed based on Kane’s partial velocities (which are Jacobians of the velocity with respect to the time derivatives of the degrees of freedom).
For instance, the linear velocity of point J is computed as:
where the Jacobians \(\frac{\partial v_J}{\partial \dot{q}_j}\) are stored in RtHSdat%PLinVelEJ(:,0)
Translational accelerations are computed as the sum of contribution from the first and second time derivatives of the degrees of freedom. For instance, the acceleration of point J is computed as:
where \(\frac{\partial a_J}{\partial \dot{q}_j}\) are stored in RtHSdat%PLinVelEJ(:,1)
Angular accelerations requires similar computations currently not documented.
4.2.7.3.3. Rotor and tail furl
The user can select linear spring and damper models, together with up- and down-stop springs, and up- and down-stop dampers.
The torque applied from the linear spring and damper is:
where \(\theta\) is the degree of freedom (rotor or tail furl),
\(k\) is the linear spring constant (RFrlSpr
or TFrlSpr
)
\(d\) is the linear damping constant (RFrlDmp
or TFrlDmp
).
The up-/down- stop spring torque is defined as:
where
\(k_{US}\) is the up-stop spring constant (RFrlUSSpr
or TFrlUSSpr
),
\(\theta_{k_{US}}\) is the up-stop spring angle (RFrlUSSP
or TFrlUSSP
),
and similar notations are used for the down-stop spring.
The up-/down- stop damping torque is defined as:
where similar nnotations are used. The total moment on the given degree of freedom is:
4.2.7.4. Developer notes
Internal coordinate systems
The different coordinate systems of ElastoDyn are stored in the variable CoordSys
.
The coordinate systems used internally by ElastoDyn are using a different convention than the OpenFAST input/output coordinate system.
For instance, for the coordinate system of the nacelle, with unit axes noted \(x_n,y_n,z_n\) in OpenFAST, and \(d_1,d_2,d_3\) in ElastoDyn, the following conversions apply: \(d_1 = x_n\), \(d_2 =z_n\) and \(d_3 =-y_n\).
The following (partial) list of coordinate systems are defined internally by ElastoDyn:
z : inertial coordinate system
a : tower base coordinate system
t : tower-node coordinate system (one per node)
d : nacelle coordinate system
c : shaft-tilted coordinate system
rf : rotor furl coordinate system
tf : tail furl coordinate system
g : hub coordinate system