vault backup: 2025-09-05 08:17:27
This commit is contained in:
aGYZ
parent
82d9876ffc
commit
a8fd01575b
@ -0,0 +1,248 @@
|
||||
---
|
||||
epoch: 1757030244783
|
||||
modelKey: gemini-2.5-pro|google
|
||||
tags:
|
||||
- copilot-conversation
|
||||
---
|
||||
|
||||
**user**: [[Kallesøe-Equations of motion for a rotor blade]]这篇文章的公式太多了，帮我梳理脉络
|
||||
[Timestamp: 2025/09/05 07:57:24]
|
||||
|
||||
**ai**: 好的，这篇文章 [[Kallesøe-Equations of motion for a rotor blade]] 确实包含了大量复杂的公式。我们来梳理一下它的核心脉络和逻辑结构，让你能更好地理解其内容。
|
||||
|
||||
这篇文章的核心目标是**建立一个能够描述风力机叶片复杂动态行为的数学模型**，这个模型特别考虑了以往模型中被简化的几个重要因素：**重力、变桨动作和风轮转速变化**。
|
||||
|
||||
整个文章的脉络可以分为以下几个部分：
|
||||
|
||||
### 1. 背景与目标 (Introduction)
|
||||
- **问题提出**: 随着风力机尺寸越来越大，叶片自身的运动（挥舞、摆振、扭转）与变桨系统、重力、风轮转速变化之间的耦合作用越来越显著，这会影响机组的疲劳载荷和控制系统设计。
|
||||
- **现有模型不足**: 经典的 Hodges-Dowell 模型虽然很完善，但主要针对直升机旋翼，没有包含重力、变桨和转速变化这些对大型风力机至关重要的因素。
|
||||
- **本文目标**: 扩展 Hodges-Dowell模型，将上述三个因素包含进来，并推导出描述变桨动作和风轮转速的独立方程。最终目标是得到一个既能进行物理解释，又能用于仿真分析的综合模型。
|
||||
|
||||
### 2. 模型的建立与假设 (Model Description)
|
||||
这是所有推导的基础，作者首先定义了研究的物理系统和坐标系。
|
||||
- **物理系统**: 一个旋转的、不可伸长的柔性叶片，具有三个自由度：
|
||||
- **挥舞 (flapwise)**: 垂直于旋转平面的弯曲。
|
||||
- **摆振 (edgewise)**: 在旋转平面内的弯曲。
|
||||
- **扭转 (torsional)**: 绕叶片展向轴线的扭转。
|
||||
- **坐标系**:
|
||||
- **惯性坐标系 (X, Y, Z)**: 固定不动的地面坐标系。
|
||||
- **旋转坐标系 (x̂, ŷ, ẑ)**: 跟随轮毂一起旋转的坐标系。
|
||||
- **叶片坐标系 (x, y, z)**: 在旋转坐标系的基础上，再绕 ẑ 轴旋转一个变桨角 `β` 得到的坐标系，叶片的变形是在这个坐标系下描述的。
|
||||
- **关键假设**: 为了简化问题，模型忽略了塔架运动、偏航、传动链柔性等因素。
|
||||
|
||||
### 3. 核心推导：运动方程的建立 (Derivation of the Equations of Motion)
|
||||
这是文章最核心、公式最密集的部分。作者采用了经典力学中的**扩展哈密顿原理 (extended Hamilton’s principle)** 来推导方程。
|
||||
|
||||
这个方法的思路是：
|
||||
1. **定义系统的能量**:
|
||||
- **动能 (Kinetic Energy)**: 包含叶片平动和转动的能量，以及整个风轮（发电机、齿轮箱等）的转动能量。
|
||||
- **势能 (Potential Energy)**: 包含两部分：
|
||||
- **弹性势能**: 叶片弯曲和扭转时内部产生的应变能。
|
||||
- **重力势能**: 叶片在重力场中的位置决定的能量。
|
||||
- **非保守力功 (Nonconservative Forces)**: 外部作用力（如气动力、发电机力矩、变桨力矩）所做的虚功。
|
||||
|
||||
2. **应用哈密顿原理**: 将上述能量和功的表达式代入作用量积分 `δH`，并令其变分为零。通过复杂的变分运算，就能自然地推导出一系列耦合的**非线性偏微分方程 (PDEs)**。
|
||||
|
||||
3. **方程的分类与解读**: 推导出的方程组可以分为几个部分，分别描述不同的运动：
|
||||
- **叶片弯曲运动方程 (Blade Bending Motion)**: 两个方程，分别描述了 `x` 方向（摆振）和 `y` 方向（挥舞）的运动。作者将方程中的复杂项 (`F_u,1`, `F_u,2` 等) 拆解开，并逐一解释了它们的物理意义，例如：离心力、科里奥利力、重力、变桨动作带来的惯性力等。
|
||||
- **叶片扭转运动方程 (Blade Torsional Motion)**: 一个方程，描述了叶片扭转角 `θ` 的动态。同样，也对各项进行了物理解释。
|
||||
- **变桨动作方程 (Pitch Action)**: 一个独立的方程，描述了整个叶片的变桨角 `β` 的动态。这个方程可以用来计算在特定叶片运动下驱动变桨所需的力矩 `M_pitch`。
|
||||
- **风轮位置方程 (Rotor Position)**: 一个描述风轮方位角 `φ` 的方程，它实际上代表了整个风轮的转速动态。这个方程将单个叶片的运动和受力与整个驱动链的转动惯量联系起来。
|
||||
|
||||
### 4. 模型的简化与应用 (Application Example)
|
||||
直接求解上述复杂的偏微分方程组非常困难。因此，作者采用了**假设模态法 (Assumed Mode Method)** 将其简化。
|
||||
- **核心思想**: 将叶片复杂的连续变形 `u(s, t)` 分解为一个空间形状函数 `u_s(s)` 和一个时间函数 `u_t(t)` 的乘积。空间形状函数（即模态振型）是预先知道的（可以通过有限元等方法计算），而未知量只剩下随时间变化的幅值 `u_t(t)`。
|
||||
- **效果**: 通过这种方法，无限自由度的**偏微分方程 (PDEs)** 被转化为了只有少数几个自由度（本文中为3个）的**常微分方程 (ODEs)**。这个常微分方程组（如文中的公式39）形式更简单，非常适合进行数值仿真和控制分析。
|
||||
|
||||
### 5. 验证与结论 (Discussion & Conclusion)
|
||||
为了证明模型的正确性和有效性，作者进行了两项验证工作：
|
||||
1. **频率对比**: 将本文模型计算出的叶片固有振动频率与成熟的商业软件 HAWCstab 的结果进行对比，两者吻合良好（见 Table I），证明模型准确地捕捉了叶片的基本结构特性。
|
||||
2. **仿真对比**: 模拟了一个快速变桨2度的工况，并将叶尖位移、变桨力矩等结果与高精度的气弹仿真软件 HAWC2 进行对比，结果也表现出很好的一致性（见 Figure 3）。
|
||||
|
||||
**总结一下，这篇文章的脉络就是：**
|
||||
|
||||
**提出问题 (经典模型不足) → 建立物理模型和坐标系 → 运用哈密顿原理推导能量表达式 → 推导出复杂的耦合偏微分运动方程组 → 对各方程的物理意义进行解读 → 采用假设模态法将PDE简化为ODE → 通过与成熟软件对比验证模型的准确性 → 得出结论：模型有效，可用于分析叶片与变桨的耦合作用。**
|
||||
|
||||
希望这个梳理能帮助你更好地理解这篇论文的整体框架和逻辑。
|
||||
[Timestamp: 2025/09/05 07:58:00]
|
||||
|
||||
**user**: 扩展哈密顿原理 (extended Hamilton’s principle)是什么
|
||||
[Timestamp: 2025/09/05 07:58:48]
|
||||
|
||||
**ai**: Of course! Let's break down the **Extended Hamilton's Principle** in a clear and structured way. It's a very powerful concept in advanced mechanics, and understanding it is key to grasping how the equations in that paper were derived.
|
||||
|
||||
The principle is best understood by starting with its simpler, more famous predecessor.
|
||||
|
||||
### 1. The Standard Hamilton's Principle (for Conservative Systems)
|
||||
|
||||
This is a fundamental principle of classical mechanics that offers an alternative, and often more elegant, way to derive the equations of motion than using Newton's laws (F=ma).
|
||||
|
||||
* **The Core Idea**: Nature is "economical." A physical system will always move between two points in time along a path that minimizes (or more accurately, makes *stationary*) a quantity called **Action**.
|
||||
|
||||
* **Key Ingredients**:
|
||||
1. **Kinetic Energy (T)**: The energy of motion.
|
||||
2. **Potential Energy (V)**: The stored energy (e.g., from gravity or elastic deformation).
|
||||
3. **The Lagrangian (L)**: This is the central quantity, defined simply as the difference between kinetic and potential energy.
|
||||
$L = T - V$
|
||||
4. **The Action (S)**: This is the integral of the Lagrangian over a specific time interval, from $t_1$ to $t_2$. It's a single number that represents the "total cost" of a particular path.
|
||||
$S = \int_{t_1}^{t_2} L \, dt = \int_{t_1}^{t_2} (T - V) \, dt$
|
||||
|
||||
* **The Principle Stated**: The actual path taken by the system is the one for which the **variation of the action is zero** ($\delta S = 0$).
|
||||
This means if you imagine slightly altering the path (a "virtual" change), the action value for the true path will be stationary (a minimum, maximum, or saddle point). This mathematical condition, $\delta S = 0$, leads directly to the famous **Euler-Lagrange equations**, which are the system's equations of motion.
|
||||
|
||||
* **Limitation**: The standard principle works perfectly for **conservative systems**, where all forces (like gravity and ideal springs) can be derived from a potential energy function. Energy is conserved in these systems.
|
||||
|
||||
### 2. The "Extended" Hamilton's Principle (for Non-Conservative Systems)
|
||||
|
||||
Now, what happens when you have forces that *don't* conserve energy? These are called **non-conservative forces**.
|
||||
|
||||
* **Examples of Non-Conservative Forces**:
|
||||
* Friction or damping (dissipate energy).
|
||||
* Aerodynamic forces (can add or remove energy).
|
||||
* Applied external forces or torques, like the generator torque ($T_{gen}$) or the pitch actuator moment ($M_{pitch}$) in the paper.
|
||||
|
||||
The standard principle can't handle these because you can't define a simple potential energy `V` for them. This is where the "extension" comes in.
|
||||
|
||||
* **The Extension**: We introduce a new term to account for the work done by these non-conservative forces. This term is called **Virtual Work ($\delta W_{nc}$)**. It represents the work done by all non-conservative forces during an infinitesimal "virtual" displacement.
|
||||
|
||||
* **The Principle Stated**: The principle is modified to say that the variation of the action is no longer zero, but is instead equal to the virtual work done by the non-conservative forces.
|
||||
|
||||
The mathematical formulation is:
|
||||
$$
|
||||
\int_{t_1}^{t_2} (\delta L + \delta W_{nc}) \, dt = 0
|
||||
$$
|
||||
Since $L = T - V$, we can write $\delta L = \delta T - \delta V$. Substituting this in gives the form you see in the paper:
|
||||
$$
|
||||
\int_{t_1}^{t_2} (\delta T - \delta V + \delta W_{nc}) \, dt = 0
|
||||
$$
|
||||
|
||||
### Connecting it to the Paper
|
||||
|
||||
Now, look at the exact expression from the [[Kallesøe-Equations of motion for a rotor blade]] note:
|
||||
$$
|
||||
\delta H \equiv \int_{t_1}^{t_2} (\delta T - \delta V_{ela} - \delta V_{gra} + \delta Q) dt = 0
|
||||
$$
|
||||
This is a direct application of the Extended Hamilton's Principle:
|
||||
* $\delta T$: The variation of the system's total **Kinetic Energy**.
|
||||
* $\delta V_{ela} + \delta V_{gra}$: The variation of the system's total **Potential Energy** (`V`), which has two parts: elastic strain energy and gravitational potential energy.
|
||||
* $\delta Q$: This is exactly the **Virtual Work** term ($\delta W_{nc}$) done by all non-conservative forces (aerodynamics, generator torque, pitch moment).
|
||||
|
||||
### Why is this Principle so Powerful?
|
||||
|
||||
1. **Energy is a Scalar**: It works with energy (a scalar quantity), which is much easier to manage than forces and accelerations (vectors), especially in complex, rotating 3D coordinate systems.
|
||||
2. **System-Level Approach**: It considers the entire system's energy at once. The complex interactions and couplings between different parts (like bending, torsion, and pitching) are automatically handled.
|
||||
3. **Generalizability**: It is extremely effective for deriving the equations of motion for continuous bodies like beams and plates (i.e., a rotor blade), resulting in the correct partial differential equations.
|
||||
4. **Automatic Boundary Conditions**: A key advantage of the variational approach is that the process not only gives you the equations of motion but also naturally derives the system's boundary conditions (e.g., the conditions at the blade tip and root).
|
||||
|
||||
In summary, the **Extended Hamilton's Principle** is a sophisticated and powerful tool that says the motion of any mechanical system, even with non-conservative forces, is governed by a single variational equation based on its kinetic energy, potential energy, and the work done by external forces. By setting up these energy expressions and applying the principle, one can systematically derive the complete and correct equations of motion.
|
||||
[Timestamp: 2025/09/05 07:59:21]
|
||||
|
||||
**user**: [[Kallesøe-Equations of motion for a rotor blade]]方程2-8是怎么得到的
|
||||
[Timestamp: 2025/09/05 08:07:30]
|
||||
|
||||
**ai**: 好的，我们来详细解析一下 [[Kallesøe-Equations of motion for a rotor blade]] 中方程 (2) 到 (8) 的推导过程。
|
||||
|
||||
这些方程都属于**势能 (Potential Energy)**部分的推导，是建立整个系统动力学方程的基础。势能在这里分为两部分：
|
||||
1. **弹性势能 ($V_{ela}$)**：叶片因弯曲和扭转变形而储存的能量。
|
||||
2. **重力势能 ($V_{gra}$)**：叶片在重力场中的位置所具有的能量。
|
||||
|
||||
我们一步步来看：
|
||||
|
||||
---
|
||||
|
||||
### 核心思想：从几何变形到能量
|
||||
|
||||
整个推导的逻辑链是：
|
||||
**几何变形 → 应变 (Strain) → 应力 (Stress) → 弹性势能 (Energy)**
|
||||
|
||||
---
|
||||
|
||||
#### 方程 (2): 格林应变张量 (Green's Strain Tensor)
|
||||
$$
|
||||
2[\mathrm{d}s,\mathrm{d}\eta,\mathrm{d}\xi][\varepsilon_{i j}][\mathrm{d}s,\mathrm{d}\eta,\mathrm{d}\xi]^{\mathrm{T}}=\mathrm{d}\mathbf{r}_{\mathrm{l}}\cdot\mathrm{d}\mathbf{r}_{\mathrm{l}}-\mathrm{d}\mathbf{r}_{\mathrm{0}}\cdot\mathrm{d}\mathbf{r}_{\mathrm{0}}
|
||||
$$
|
||||
* **目的**: 这个方程是**应变**的数学定义。应变是衡量物体变形程度的物理量。
|
||||
* **为什么用这个定义**: 对于风力机叶片这种会发生大位移和转动的结构，简单的“工程应变”定义（长度变化/原始长度）是不够准确的。格林应变张量是一个更精确的定义，它能正确处理大转动下的变形。
|
||||
* **物理意义**:
|
||||
* `d**r**_0`: 想象在**未变形**的叶片上画一个无限小的矢量（一小段线）。
|
||||
* `d**r**_1`: 叶片变形后，原来那段小线段变成了新的矢量 `d**r**_1`。
|
||||
* `d**r**_1 ⋅ d**r**_1 - d**r**_0 ⋅ d**r**_0`: 这个表达式计算的是这段小线段**长度平方的变化量**。如果叶片没有变形，`d**r**_1` 和 `d**r**_0` 只是在空间中转动了一下，长度不变，所以这个差值为零。只有当线段被拉伸、压缩或剪切时，这个差值才不为零。
|
||||
* `ε_ij`: 这就是应变张量，它是一个数学工具，用来系统地描述这个长度平方的变化。
|
||||
|
||||
---
|
||||
|
||||
#### 方程 (3) 和 (4): 变形前后的位置矢量
|
||||
为了使用方程 (2)，我们必须先用数学语言写出 `**r**_0` 和 `**r**_1` 的表达式。
|
||||
|
||||
* **方程 (3): 未变形时的位置矢量 `**r**_0`**
|
||||
$$
|
||||
\mathbf{r}_{0}\!=\![\mathbf{I},\mathbf{J},\mathbf{K}]\mathbf{T}_{\phi}^{\mathrm{T}}\mathbf{T}_{\beta}^{\mathrm{T}}\!\left[[l_{p i},0,w_{0}]^{\mathrm{T}}+(\mathbf{T}_{e}|_{u=\nu=\theta=0})^{\mathrm{T}}[\eta_{0},\xi_{0},0]^{\mathrm{T}}\right]
|
||||
$$
|
||||
这个公式描述了叶片上任意一点在变形前的位置。我们从内向外看：
|
||||
1. `[η_0, ξ_0, 0]^T`: 这是叶片横截面内任意一点相对于**弹性轴**的坐标。
|
||||
2. `(**T**_e|_{u=ν=θ=0})^T`: 将这个截面坐标转换到叶片坐标系 `(x,y,z)` 下。
|
||||
3. `[l_pi, 0, w_0]^T`: 这是弹性轴本身在叶片坐标系 `(x,y,z)` 中的位置。两者相加，就得到了该点在 `(x,y,z)` 坐标系中的完整位置。
|
||||
4. `**T**_φ^T **T**_β^T`: 这一串矩阵是坐标变换，把 `(x,y,z)` 坐标系下的位置转换回**惯性坐标系 `(X,Y,Z)`**。
|
||||
|
||||
* **方程 (4): 变形后的位置矢量 `**r**_1`**
|
||||
$$
|
||||
\mathbf{r}_{\mathrm{1}}\!=\![\mathbf{I},\mathbf{J},\mathbf{K}]\mathbf{T}_{\phi}^{\mathrm{T}}\mathbf{T}_{\beta}^{\mathrm{T}}\big[[l_{p i}+u,\nu,w_{0}+w_{\mathrm{1}}]^{\mathrm{T}}+\mathbf{T}_{\mathrm{c}}^{\mathrm{T}}{[\eta_{\mathrm{1}},\xi_{\mathrm{1}},0]}^{\mathrm{T}}\big]
|
||||
$$
|
||||
这个公式的结构和方程 (3) 完全一样，只是加入了变形量：
|
||||
* 弹性轴的位置从 `(l_pi, 0, w_0)` 变成了 `(l_pi + u, ν, w_0 + w_1)`。
|
||||
* 截面内的点也可能发生微小移动，从 `(η_0, ξ_0)` 变为 `(η_1, ξ_1)`。
|
||||
|
||||
---
|
||||
|
||||
#### 方程 (5): 具体的应变分量
|
||||
$$
|
||||
\begin{array}{l l}{{\displaystyle\varepsilon_{11}\!=\!-u^{\prime\prime}\big(\eta\cos(\overline{{\theta}})-\xi\sin\!\big(\hat{\theta}\big)\big)-\nu^{\prime\prime}\big(\eta\sin\!\big(\hat{\theta}\big)+\xi\cos\!\big(\overline{{\theta}}\big)\big)}}\\ {{\displaystyle\varepsilon_{12}\!=\!-\frac{1}{2}\xi\theta^{\prime},\quad\!\varepsilon_{13}\!=\!\frac{1}{2}\eta\theta^{\prime}}}\end{array}
|
||||
$$
|
||||
* **如何得到**: 这是一个非常繁琐的代数过程，但思路很直接：
|
||||
1. 将方程 (3) 和 (4) 的表达式代入方程 (2)。
|
||||
2. 进行大量的矩阵和矢量运算。
|
||||
3. 应用文中的**量级排序方案 (Order Scheme)**，忽略掉所有高阶小量（例如 `ε^3`, `ε^4` 等）。这是至关重要的一步，它极大地简化了最终的表达式。
|
||||
* **物理意义**:
|
||||
* `ε_11`: **轴向应变**（拉伸或压缩）。这个公式表明，轴向应变主要是由叶片的**弯曲曲率** (`u''` 和 `ν''`) 引起的。这符合梁理论的基本知识：梁弯曲时，外侧被拉伸，内侧被压缩。
|
||||
* `ε_12`, `ε_13`: **剪切应变**。公式表明，剪切应变主要是由叶片的**扭转率** (`θ'`) 引起的。这也符合经典的扭转理论。
|
||||
|
||||
---
|
||||
|
||||
#### 方程 (6): 弹性势能的变分
|
||||
$$
|
||||
\delta V_{e l a}=\int_{r}^{R}\iint_{A}{(\sigma_{s s}\delta\varepsilon_{s s}+\sigma_{s\eta}\delta\varepsilon_{s\eta}+\sigma_{s\xi}\delta\varepsilon_{s\xi})\mathrm{d}\eta\mathrm{d}\xi\mathrm{d}s}
|
||||
$$
|
||||
* **目的**: 将应变和应力联系起来，计算总的弹性势能。
|
||||
* **推导**: 这是弹性力学中计算应变能的标准公式。
|
||||
* `σ * δε`: 应力乘以应变的变分，代表了单位体积内弹性势能的变分。
|
||||
* `∬_A ... dη dξ`: 对整个叶片**横截面**进行积分。
|
||||
* `∫_r^R ... ds`: 沿着叶片**长度**从根部 `r` 到叶尖 `R` 进行积分。
|
||||
* 最终，通过胡克定律（`σ = Eε`）将应力 `σ` 替换为用应变 `ε` 表示的式子，再代入方程 (5) 的结果，就可以得到一个完全由几何变形 (`u`, `ν`, `θ` 及其导数) 表达的弹性势能 `V_ela`。
|
||||
|
||||
---
|
||||
|
||||
#### 方程 (7) 和 (8): 重力势能
|
||||
|
||||
* **方程 (7): 重力势能的定义**
|
||||
$$
|
||||
V_{g r a}=\int_{r}^{R}\mathbf{r}_{c g}^{T}\cdot\mathbf{g}\mathrm{d}s
|
||||
$$
|
||||
* **目的**: 计算整个叶片在重力场中的势能。
|
||||
* **推导**: 这是重力势能的基本公式 `mgh` 的积分形式。
|
||||
* `**g** = [0, 0, -g]^T`: 重力场矢量。
|
||||
* `**r**_cg`: 叶片上某个截面**重心 (center of gravity)** 的位置矢量。
|
||||
* `**r**_cg^T ⋅ **g**`: 计算重力与位置矢量的点积，这与 `-mgh` 成正比。
|
||||
* `∫_r^R ... ds`: 沿着叶片长度积分，把每一段的重力势能加起来。
|
||||
|
||||
* **方程 (8): 重心的位置矢量 `**r**_cg`**
|
||||
$$
|
||||
\mathbf{r}_{c g}\!=\![\mathbf{I},\mathbf{J},\mathbf{K}]\mathbf{T}_{\phi}^{\mathrm{T}}\mathbf{T}_{\beta}^{\mathrm{T}}\big[[l_{p i}+u,\nu,w_{0}+w_{1}]^{\mathrm{T}}+\mathbf{T}_{\mathrm{c}}^{\mathrm{T}}[l_{c g},0,0]^{\mathrm{T}}\big]
|
||||
$$
|
||||
* **目的**: 为了计算方程 (7)，必须先精确地写出重心的位置。
|
||||
* **推导**: 这个公式的结构与方程 (4)（变形后任意一点的位置）几乎完全相同。唯一的区别在于，它描述的不是截面内任意一点 `(η, ξ)`，而是特指**重心 `cg`** 这一点。
|
||||
* 根据文章假设，重心 `cg` 位于弦线上，距离弹性轴 `ea` 为 `l_cg`。因此，它在弦线坐标系下的坐标就是 `[l_cg, 0, 0]^T`。
|
||||
* 所以，方程 (8) 就是将方程 (4) 中的 `[η_1, ξ_1, 0]^T` 替换为 `[l_cg, 0, 0]^T` 得到的。
|
||||
|
||||
总结来说，方程 (2) 到 (8) 是一套完整且逻辑严密的推导链，它从最基本的几何变形定义出发，通过应用连续介质力学和坐标变换，最终得到了计算叶片**弹性势能**和**重力势能**所需的全部数学表达式。这些表达式是后续使用哈密顿原理推导运动方程的必要输入。
|
||||
[Timestamp: 2025/09/05 08:08:16]
|
||||
Loading…
x
Reference in New Issue
Block a user