obsidian_backup/书籍/力学书籍/OpenFast/FASTKinetics/auto/Drivertrain.md

The inertia of the drivetrain brings about generalized inertia forces and the load in the generator, high-speed shaft brake, gearbox (friction forces resulting from nonzero GBoxEff ) and the flexibility of the low speed shaft bring about generalized active forces. Note that all of these equations assume that the rotor is spinning about the positive c_{I} axis (they assume that the rotor can’t be forced to rotate in the opposite direction). This model works for any gearbox arrangement (including no gearbox, single stage, or multi-stage) as long as the generator rotates about the shaft axis (it may not be skewed relative to the shaft, even though it may rotate in the opposite direction of the low-speed shaft due to the gearbox stages). If there is no gearbox, simply set G B R a t i o\,=\,G B o x E f f=G e n D i r=1 (GB R e\nu e r s e={\mathrm{False}} ).
传动系惯性带来广义惯性力和发电机负载、高速轴制动器、齿轮箱（由于GBoxEff非零导致的摩擦力）以及低速轴挠曲性，都会产生广义主动力。需要注意的是，这些方程都假设转子绕正$c_{I}$轴旋转（即假设转子不能被强制旋转到相反方向）。该模型适用于任何齿轮箱布置（包括无齿轮箱、单级或多级），只要发电机绕轴旋转即可（即使它可能由于齿轮箱级数而以低速轴相反的方向旋转，但不能相对于轴倾斜）。如果不存在齿轮箱，只需将G B R a t i o\,=\,G B o x E f f=G e n D i r=1 (GB R e\nu e r s e={\mathrm{False}} )。

The mechanical torque within the generator is applied to the high speed shaft and equally and oppositely to the structure that furls with the rotor as follows:
发电机内部的机械扭矩施加到高速轴上，并以相等且相反的方式作用于与转子一起摆动的结构，具体如下：

\left.F_{r}\right|_{G e n}=\left(^{E}\omega_{r}^{G}-{^{E}\!\omega_{r}^{R}}\right)\cdot M_{G e n}^{G}\quad\left(r=l,2,...,22\right)

Thus,

\left.F_{r}\right|_{G e n}=\left\{\!\!\!\begin{array}{l l}{{^{E}\pmb{\omega}_{G e d z}^{G}\cdot M_{G e n}^{G}}}&{{f o r~~r=G e A z}}\\ {{0}}&{{o t h e r w i s e}}\end{array}\ \right.

\boldsymbol{M_{G e n}^{G}}=-G e n D i\boldsymbol{r}\cdot\boldsymbol{T^{G e n}}\left(G B R a t i o\cdot\dot{q}_{G e d z},t\right)\boldsymbol{c}_{I}

Note that a positive T^{G e n} represents a load (positive power extracted) whereas a negative T^{G e n} represents a motoring-up situation (negative power extracted, or power input). Thus,

F_{r}\Big|_{G e n}=\left\{\begin{array}{l l}{\left(G e n D i r\cdot G B R a t i o\cdot c_{I}\right)\cdot\left[-G e n D i r\cdot T^{G e n}\left(G B R a t i o\cdot\dot{q}_{G e d z},t\right)c_{I}\right]}&{f o r\;\;\;r=G e n D i r}\\ {0}&{o t h e r w i s e}\end{array}\right.

Or since G e n D i r^{2}=1 ,

F_r \Big|_{Gen} = \begin{cases}

-GRRatio \cdot T^{Gen} (GRRatio \cdot \dot{q}_{GeAz}, t) & \text{for } r = GeAz \\

0 & \text{otherwise}

\end{cases}

Thus,

[C(q,t)]|_{\text{Gen}} = 0

\{-f(\dot{q},q,t)\}|_{\text{Gen}} = \left\{
\begin{array}{l}
\ldots \\
\ldots \\
\ldots \\
\ldots \\
\ldots \\
-GBRatio \cdot T^{\text{Gen}} (GBRatio \cdot \dot{q}_{\text{GeAz}}, t) \\
\ldots \\
\ldots \\
\ldots \\
\ldots \\
\ldots
\end{array}
\right\}

Similarly, the mechanical torque applied to the high-speed shaft from the high-speed shaft brake is applied equally and oppositely to the structure that furls with the rotor. Thus,
同样地，从高速轴制动器施加到高速轴上的机械扭矩，也会以相等且相反的方式作用于与转子一同回转的结构上。因此，

F_{r}\big|_{B r a k e}=\left\{\!\!\!\begin{array}{l l}{{^{E}\pmb{\omega}_{G e A z}^{G}\cdot M_{B r a k e}^{G}}}&{{f o r\;\;\;r=G e A z}}\\ {{0}}&{{o t h e r w i s e}}\end{array}\!\!\right. where M_{B r a k e}^{G}=-G e n D i r\cdot T^{B r a k e}\left(t\right)c_{I}\qquad\mathrm{~and~where}\qquad T^{B r a k e}\left(t\right)=H S S B r k T\left(t\right)

which is assumed to be positive in value always. Thus,

F_{r}\bigr|_{B r a k e}=\left\{\begin{array}{l l}{-G B R a t i o\cdot T^{B r a k e}\left(t\right)}&{f o r\;\;\;r=G e A z}\\ {O}&{o t h e r w i s e}\end{array}\right.

Thus,

[C(q,t)]|_{\text{Brake}} = 0

\{-f(\dot{q},q,t)\}|_{\text{Brake}} = \left\{
\begin{array}{l}
\ldots \\
\ldots \\
\ldots \\
\ldots \\
\ldots \\
\ldots \\
-GBRatio \cdot T^{\text{Brake}}(t) \\
\ldots \\
\ldots \\
\ldots \\
\ldots \\
\ldots \\
\ldots
\end{array}
\right\}

If the translational inertia of the drivetrain is assumed to be incorporated into that of the structure that furls with the rotor, then the high-speed shaft generator inertia generalized force is as follows: 如果假设动能惯性由驱动系传递至随转子摆动的结构部分，则高速轴发电机惯性概括力如下：

F_r^*|_{\text{G}} = {}^E \omega_r^{\text{G}} \cdot (-\overline{\overline{{I}}}^{\text{G}} \cdot {}^E {a}^{\text{G}} - {}^E \boldsymbol{\omega}^{\text{G}} \times \overline{\overline{{I}}}^{\text{G}} \cdot {}^E \boldsymbol{\omega}^{\text{G}}) \quad \text{where} \quad \overline{\overline{{I}}}^{\text{G}} = \text{GenInerc}_{1c1} \tag{r=1,2,...,22}

or,

F_r^*|_{\text{G}} = {}^E \boldsymbol{\omega}_r^{\text{G}} \cdot \left\{ -\overline{\overline{{I}}}^{\text{G}} \cdot \left\{ \left( \sum_{i=4}^{13} {}^E \omega_i^{\text{G}} \ddot{q}_i \right) + \left[ \sum_{i=7}^{13} \frac{d}{dt} \left( {}^E \boldsymbol{\omega}_i^{\text{G}} \right) \dot{q}_i \right] \right\} - {}^E \boldsymbol{\omega}^{\text{G}} \times \overline{\overline{{I}}}^{\text{G}} \cdot {}^E \boldsymbol{\omega}^{\text{G}} \right\} \tag{r=1,2,...,22}

However, since c_{1}\cdot\frac{d}{d t}\Big(\{^{E}\omega_{G e A z}^{G}\Big)\propto c_{1}\cdot(^{E}\omega^{R}\times c_{1}\Big)={^{E}\omega^{R}}\cdot\Big(c_{1}\times c_{1}\Big)=0 (the first c_{1} coming from \overline{\overline{{I}}}^{\text{G}} ), this simplifies as follows:

F_r^*|_{\text{G}} = {}^E \boldsymbol{\omega}_r^{\text{G}} \cdot \left\{ -\overline{\overline{{I}}}^{\text{G}} \cdot \left\{ \left( \sum_{i=4}^{13} {}^E \boldsymbol{\omega}_i^{\text{G}} \ddot{q}_i \right) + \left[ \sum_{i=7}^{12} \frac{d}{dt} \left( {}^E \boldsymbol{\omega}_i^{\text{G}} \right) \dot{q}_i \right] \right\} - {}^E \boldsymbol{\omega}^{\text{G}} \times \overline{\overline{{I}}}^{\text{G}} \cdot {}^E \boldsymbol{\omega}^{\text{G}} \right\} \tag{r=1,2,...,22}
$$Or, 

F_r^*|_{G} = \begin{cases}

• {}^E \boldsymbol{\omega}r^{\text{R}} \cdot \overline{\overline{{I}}}^{\text{G}} \cdot \left{ \left( \sum{i=4}^{13} {}^E \boldsymbol{\omega}_i^{\text{G}} \ddot{q}i \right) + \left[ \sum{i=7}^{12} \frac{d}{dt} \left( {}^E \boldsymbol{\omega}_i^{\text{R}} \right) \dot{q}_i \right] \right} - ({}^E \boldsymbol{\omega}_r^{\text{R}} \cdot ({}^E \boldsymbol{\omega}^{\text{G}} \times \overline{\overline{{I}}}^{\text{G}} \cdot {}^E \boldsymbol{\omega}^{\text{G}})) & \text{for } r = 4,5,...,12 \
• \text{GenDir} \cdot \text{GenIner} \cdot \text{GBRatio} \left{ \left( \left( \sum_{i=4}^{12} {}^E \boldsymbol{\omega}_i^{\text{R}} \ddot{q}i \right) + \left[ \sum{i=7}^{12} \frac{d}{dt} \left( {}^E \boldsymbol{\omega}_i^{\text{R}} \right) \dot{q}i \right] \right) \cdot c_1 - \text{GenIner} \cdot \text{GBRatio}^2 \cdot \ddot{q}{\text{GeAz}} \right} - \text{GenDir} \cdot \text{GBRatioc}_1 \cdot ({}^E \boldsymbol{\omega}^{\text{G}} \times \overline{\overline{{I}}}^{\text{G}} \cdot {}^E \boldsymbol{\omega}^{\text{G}}) & \text{for } r = \text{GeAz} \ 0 & \text{otherwise} \end{cases}

However since $c_{1}\cdot\left({\,^{E}}{\omega^{G}}\times c_{1}\,\right)={^{E}}{\omega^{G}}\cdot\left(c_{1}\times c_{1}\right)=0$ (the first $c_{I}$ coming from $\overline{{\overline{{I}}}}^{G}$ ), this simplifies again as follows:  

F_r^*|_{G} = \begin{cases}

• {}^E \boldsymbol{\omega}r^{\text{R}} \cdot \overline{\overline{{I}}}^{\text{G}} \cdot \left{ \left( \sum{i=4}^{13} {}^E \boldsymbol{\omega}_i^{\text{G}} \ddot{q}i \right) + \left[ \sum{i=7}^{12} \frac{d}{dt} \left( {}^E \boldsymbol{\omega}_i^{\text{R}} \right) \dot{q}_i \right] \right} - ({}^E \boldsymbol{\omega}_r^{\text{R}} \cdot ({}^E \boldsymbol{\omega}^{\text{G}} \times \overline{\overline{{I}}}^{\text{G}} \cdot {}^E \boldsymbol{\omega}^{\text{G}})) & \text{for } r = 4,5,...,12 \
• \text{GenDir} \cdot \text{GenIner} \cdot \text{GBRatio} \left{ \left( \left( \sum_{i=4}^{12} {}^E \boldsymbol{\omega}_i^{\text{R}} \ddot{q}i \right) + \left[ \sum{i=7}^{12} \frac{d}{dt} \left( {}^E \boldsymbol{\omega}_i^{\text{R}} \right) \dot{q}i \right] \right) \cdot c_1 - \text{GenIner} \cdot \text{GBRatio}^2 \cdot \ddot{q}{\text{GeAz}} \right} - \text{GenDir} \cdot \text{GBRatioc}_1 \cdot ({}^E \boldsymbol{\omega}^{\text{G}} \times \overline{\overline{{I}}}^{\text{G}} \cdot {}^E \boldsymbol{\omega}^{\text{G}}) & \text{for } r = \text{GeAz} \ 0 & \text{otherwise} \end{cases}

The terms associated with DOFs 4,5,…,12 represent the fact that the rate of change of angular momentum of the generator can be considered as an additional torque on the structure that furls with the rotor (i.e., in addition to the torques on the structure transmitted directly from the low-speed shaft).  
与DOF 4、5、…、12相关的项，表示发电机角动量变化率可以被视为一个附加于随转子展开的结构上的扭矩（即，除了直接从低速轴传递到结构上的扭矩之外）。

[C(q,t)]G = \begin{cases} {}^E \boldsymbol{\omega}{\text{Row}}^{\text{G}} \cdot \overline{\overline{{I}}}^{\text{G}} \cdot {}^E \boldsymbol{\omega}_{\text{Col}}^{\text{G}} & \text{for (Row, Col = 4,5,...,13)} \ 0 & \text{otherwise} \end{cases}

{-f(\dot{q},q,t)}_G = \begin{cases}

• {}^E \boldsymbol{\omega}{\text{Row}}^{\text{G}} \cdot \left{ \overline{\overline{{I}}}^{\text{G}} \cdot \left[ \sum{i=7}^{12} \frac{d}{dt} \left( {}^E \boldsymbol{\omega}_i^{\text{R}} \right) \dot{q}_i \right] + {}^E \boldsymbol{\omega}^{\text{G}} \times \overline{\overline{{I}}}^{\text{G}} \cdot {}^E \boldsymbol{\omega}^{\text{G}} \right} & \text{for (Row = 4,5,...,12)} \
• {}^E \boldsymbol{\omega}{\text{Row}}^{\text{G}} \cdot \overline{\overline{{I}}}^{\text{G}} \cdot \left[ \sum{i=7}^{12} \frac{d}{dt} \left( {}^E \boldsymbol{\omega}_i^{\text{R}} \right) \dot{q}_i \right] & \text{for (Row = 13)} \ 0 & \text{otherwise} \end{cases}

F_r |{\text{ElasticDrive}} = -\frac{\partial V{\text{Drive}}}{\partial q_r} \quad (r = 1,2,...,22) \quad \text{where} \quad V_{\text{Drive}} = \frac{1}{2} \text{DTTorSpr} \cdot q_{\text{DrTr}}^2

So, 

F_r |{\text{ElasticDrive}} = \begin{cases} -DTTorSpr \cdot q{DrTr} & \text{for } r = DrTr \ 0 & \text{otherwise} \end{cases}

and likewise

F_r |{\text{DampDrive}} = \begin{cases} -DTTorDmp \cdot \dot{q}{DrTr} & \text{for } r = DrTr \ 0 & \text{otherwise} \end{cases}

Thus, 

[C(q,t)]_{\text{ElasticDrive}} = 0

{-f(\dot{q},q,t)}{\text{ElasticDrive}} = \begin{Bmatrix} \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ -DTTorSpr \cdot q{DrTr} \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \end{Bmatrix}

[C(q,t)]_{\text{DampDrive}} = 0

{-f(\dot{q},q,t)}{\text{DampDrive}} = \begin{Bmatrix} \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ -DTTorDmp \cdot \dot{q}{DrTr} \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \end{Bmatrix}

Similar to the generator and HSS brake, the mechanical friction torque applied to the high speed shaft is applied equally and oppositely to the structure that furls with the rotor.  Thus,  
类似于发电机和HSS制动器，施加到高速轴上的机械摩擦扭矩，以相等且相反的方式作用于与转子一同旋转的结构上。因此，

F_r |{\text{GBFric}} = \begin{cases} {}^E\omega{\text{GeAz}}^G \cdot M_{\text{GBFric}}^G & \text{for } r = \text{GeAz} \ 0 & \text{otherwise} \end{cases} \quad \text{where} \quad M_{\text{GBFric}}^G = \frac{T_{\text{GBFric}}^G(\ddot{q}, \dot{q}, q, t) c_1}{\text{GBRatio} \cdot \text{GenDir}}

where, from a free-body diagram of the high and low-speed shafts, it is easily seen that the friction torque applied on the LSS upon the gearbox entrance, $T^{G B F r i c}\left(\ddot{q},\dot{q},q,t\right)$ , is always positive in value and equal to: 
从高速轴和低速轴的受力图上可以清楚地看到，作用于低速轴进入变速箱处的摩擦转矩，$T^{G B F r i c}\left(\ddot{q},\dot{q},q,t\right)$，其值始终为正，且等于：

T^{\text{GBFric}}(\ddot{q}, \dot{q}, q, t) = \left[ \text{GenIner} \cdot \text{GBRatio}^2 \cdot \ddot{q}{\text{GeAz}} + \text{GenDir} \cdot \text{GenIner} \cdot \text{GBRatio} \cdot a^E \cdot c_1^R + \text{GBRatio} \cdot T^{\text{Gen}}(\text{GBRatio} \cdot \dot{q}{\text{GeAz}}, t) + \text{GBRatio} \cdot T^{\text{Brake}}(t) \right] \cdot \left[ \frac{1}{\text{GBBoxEff}^{\text{SIGN(LSShftTq)}}} - 1 \right]

Thus, 

F_r |_{\text{GBFric}} = \begin{cases} -T^{\text{GBFric}}(\ddot{q}, \dot{q}, q, t) & \text{for } r = \text{GeAz} \ 0 & \text{otherwise} \end{cases}

F_r |_{\text{GBFric}} = \begin{cases}

• \left[ \text{GenIner} \cdot \text{GBRatio}^2 \cdot \ddot{q}{\text{GeAz}} + \text{GenDir} \cdot \text{GenIner} \cdot \text{GBRatio}^E a^R \cdot c_1 + \text{GBRatio} \cdot T^{\text{Gen}}(\text{GBRatio} \cdot \dot{q}{\text{GeAz}}, t) + \text{GBRatio} \cdot T^{\text{Brake}}(t) \right] \cdot \left[ \frac{1}{\text{GBoxEff}^{\text{SIGN(LSShftTq)}}} - 1 \right] & \text{for } r = \text{GeAz} \ 0 & \text{otherwise} \end{cases}

Or, 

F_r |_{\text{GBFric}} = \begin{cases}

• \left[ \text{GenIner} \cdot \text{GBRatio}^2 \cdot \ddot{q}{\text{GeAz}} \right. \ \quad \left. + \text{GenDir} \cdot \text{GenIner} \cdot \text{GBRatio} \left{ \sum{i=4}^{12} {^E \omega_i^R \ddot{q}i} + \left[ \sum{i=7}^{12} \frac{d}{dt} ({^E \omega_i^R}) \dot{q}i \right] \right} \cdot c_1 \right. \ \quad \left. + \text{GBRatio} \cdot T^{\text{Gen}}(\text{GBRatio} \cdot \dot{q}{\text{GeAz}}, t) + \text{GBRatio} \cdot T^{\text{Brake}}(t) \right] \cdot \left[ \frac{1}{\text{GBoxEff}^{\text{SIGN(LSShftTq)}}} - 1 \right] & \text{for } r = \text{GeAz} \ 0 & \text{otherwise} \end{cases}

Thus, 

[C(q,t)]{\text{GBFric}} = \begin{cases} {^E\omega{\text{Row}}^G \cdot \overline{\overline{{I}}}^{\text{G}} \cdot {^E\omega_{\text{Col}}^G}} \left[ \frac{1}{\text{GBoxEff}^{\text{SIGN(LSShftTq)}}} - 1 \right] & \text{for } (\text{Row}=13, \text{Col}=4,5,...,13) \ 0 & \text{otherwise} \end{cases}

{-f(\dot{q},q,t)}_{\text{GBFric}} = \begin{cases} \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \ \ldots \

• \left{ {^E\omega_{\text{GeAz}}^G \cdot \overline{\overline{{I}}}^{\text{G}} \cdot \left[ \sum_{i=7}^{12} \frac{d}{dt} ({^E\omega_i^R}) \dot{q}i \right]} + \text{GBRatio} \cdot T^{\text{Gen}}(\text{GBRatio} \cdot \dot{q}{\text{GeAz}}, t) + \text{GBRatio} \cdot T^{\text{Brake}}(t) \right} \cdot \left[ \frac{1}{\text{GBoxEff}^{\text{SIGN(LSShftTq)}}} - 1 \right] \ \ldots \ \ldots \ \ldots \ \ldots \ \end{cases}

It is noted that the gearbox friction generalized active force effects both the mass matrix and the forcing function.  Its effect on the mass matrix can be thought of as an apparent mass in the system (i.e., an active friction force which is a function of the generator acceleration).  The gearbox friction causes the mass matrix to become unsymmetric.  Note that the equation for DOF GeAz is an implicit equation (since the gearbox friction is a function of DOF GeAz), which has no closed-form solution.  To avoid the complications resulting from this implicit characterization, the value of the LSShftTq  from the previous time step is used in the SIGN() function in place of the value of the LSShftTq  in the current time step.  This may be done since it may be assumed that LSShftTq  will not change much between time steps if the time step is considered small enough.  Note that gearbox friction is the only term in the equations of motion that cause the mass matrix to be unsymmetrical.  The mass matrix will only be unsymmetric if GBoxEff  is not $100\%$ —if GBoxEff  is $100\%$ , then the mass matrix is completely symmetric.  
需要注意的是，齿箱摩擦产生的广义主动力效应既影响质量矩阵，也影响强制函数。它对质量矩阵的影响可以被认为是系统中的视质量（即，一个随发电机加速度变化的活动摩擦力）。齿箱摩擦导致质量矩阵变得非对称。需要注意的是，DOF GeAz 的方程是一个隐式方程（因为齿箱摩擦是 DOF GeAz 的函数），因此没有闭合解。为了避免由此隐式特性带来的复杂性，在 SIGN() 函数中，使用上一个时间步长的 LSShftTq 值代替当前时间步长的 LSShftTq 值。这是可以做到的，因为如果时间步长足够小，可以假设 LSShftTq 在时间步长之间不会发生太大变化。需要注意的是，齿箱摩擦是导致质量矩阵非对称的唯一项。只有当 GBoxEff 不是 100% 时，质量矩阵才会是非对称的——如果 GBoxEff 是 100%，那么质量矩阵将完全对称。