1

力学书籍/OpenFast/modeling of the uae wind turbine for refinement of fast_ad/auto/modeling of the uae wind turbine for refinement of fast_ad.md

@@ -1091,63 +1091,42 @@ u(z,t)\!=\sum_{b=p}^{N+p-I}\!\varphi_{b}(z)c_{b}(t)
 $$  
 
 where $c_{b}(t)$ is the generalized coordinate associated with the function $\varphi_{b}(z)$ .  The $\varphi_{b}(z)$ ’s are known as shape functions and the parameter $p$ is chosen for convenience.  
+其中，$c_{b}(t)$ 是与函数 $\varphi_{b}(z)$ 相关的广义坐标。$\varphi_{b}(z)$ 被称为形函数，参数 $p$ 为了方便起见而选择。
 
 Since the shape functions are not unique to each normal mode, meaning that each normal mode is related to all shape functions, there is a relationship such that the normal mode shapes form a linear combination of the shape functions:  
 
+
+由于形函数并非每个正模特有的，意味着每个正模都与所有形函数相关，因此存在一种关系，使得正模形函数构成形函数的线性组合：
 $$
 \phi_{a}(z)\!=\!\sum_{b=p}^{N+p-l}\!C_{a,b}\varphi_{b}(z)\ \ \left(a={1,2,...,N}\right)
 $$  
 
 where $C_{a,b}$ is the constant proportionality coefficient associated with the $b^{\mathrm{th}}$ shape function and the $a^{\mathrm{th}}$ normal mode.  This is known as the Rayleigh-Ritz method.  
+其中，$C_{a,b}$ 是与第 $b$ 个形函数和第 $a$ 个正模相关的比例常数。 这被称为瑞利-里兹法。
 
 In FAST_AD, each normal mode shape is assumed to be expressible as a polynomial; thus, the $b^{\mathrm{th}}$ shape function is defined as:  
+在 FAST_AD 中，假设每个固有振型都可以用多项式表达；因此，第 $b^{\mathrm{th}}$ 形变函数定义如下：
 
 $$
 \varphi_{b}(z)\!=\!\left(\frac{z}{Z}\right)^{b}
 $$  
 
 Since the slope of a cantilevered beam must be zero at the fixed end, $p$ must be no smaller than two if the shape functions are to accurately represent the normal mode shapes.  Thus, the FAST_AD design code requires that $p$ equal two.  FAST_AD allows $N$ to be as high as five.  In FAST_AD, the constant proportionality coefficients associated with each shape function and normal mode are parameters requested in the input file.  A preprocessor code entitled Modes enables users of FAST_AD to obtain these parameters.  The theory employed by Modes is also developed here.  
-FAST_AD 的结构模型将叶片和塔架视为具有连续质量和刚度分布的灵活悬臂梁。理论上，这些体系拥有无穷多个自由度（DOFs），因为需要无数坐标来指定体系中每一点的位置。实际上，这些体系被建模为主要正常振动模式已知形状的线性和。此技术称为正常模态求和法，将自由度数从无穷减少到 $N_{\ast}$，即考虑为主导的正常模态数量。使用该方法，灵活梁在任何时间 $t$ 和位置 $z$ 的侧向挠度（垂直于未变形梁） $u(z,t)$ 可以表示为各个正常模态形状 $\phi_{a}(z)$ 与其相关的广义坐标 $q_{a}(t)$ 乘积之和：
+由于悬臂梁的slope必须在其固定端为零，如果希望形状函数能够准确地表示模态形状，则 $p$ 必须不小于二。因此，FAST_AD 设计代码要求 $p$ 等于二。FAST_AD 允许 $N$ 高达五。在 FAST_AD 中，与每个形状函数和模态相关的比例常数是输入文件中要求的参数。名为 Modes 的预处理器代码允许 FAST_AD 的用户获得这些参数。Modes 所采用的理论也在此进行阐述。
 
 $$
 u(z,t)\!=\!\sum_{a=l}^{N}\phi_{a}\!\left(z\right)\!q_{a}\!\left(t\right)
 $$  
 
-固有模态形状 $\phi_{a}(z)$ 对于模式 $a$，仅是沿梁的距离 $z$ 的函数（其中 $z=0$ 在固定端，$z=Z$ 在自由端），而与固有模态 $a$ 相关的广义坐标 $q_{a}(t)$ 仅是时间 $t$ 的函数。每个正常模态都有一个相关的固有频率 $\omega_{a}$ 和相位 $\psi_{a}$。通常，**与正常模态相关联的广义坐标被允许为悬臂梁自由端的挠度**；**因此，每个正常模态形状是无量纲的，并且在自由端归一化**。
-
-当各个固有模态形状已知时，需要 $N$ 个参数来指定任何时间点上灵活体系的挠度。因此，替代地，灵活体系的侧向挠度可以用 $N$ 个其他函数 $\varphi_{b}(z)$ 表示，这些函数不是每个正常模态独特的：
-
-$$
-u(z,t)\!=\sum_{b=p}^{N+p-I}\!\varphi_{b}(z)c_{b}(t)
-$$  
-
-其中 $c_{b}(t)$ 是与函数 $\varphi_{b}(z)$ 相关联的广义坐标。$\varphi_{b}(z)$ 称为形状函数，参数 $p$ 由便利性决定。
-
-由于形状函数不是每个正常模态独特的，即每个正常模态与所有形状函数相关联，因此存在关系使得正常模态形状可以表示为形状函数的线性组合：
-
-$$
-\phi_{a}(z)\!=\!\sum_{b=p}^{N+p-l}\!C_{a,b}\varphi_{b}(z)\ \ \left(a={1,2,...,N}\right)
-$$  
-
-其中 $C_{a,b}$ 是与第 $b$ 形状函数和第 $a$ 固有模态相关联的比例系数。这称为雷利-里兹方法。
-
-在 FAST_AD 中，假设每个正常模态形状可以表示为多项式；因此，第 $b$ 形状函数定义为：
-
-$$
-\varphi_{b}(z)\!=\!\left(\frac{z}{Z}\right)^{b}
-$$  
-
-由于悬臂梁在固定端的斜率必须为零，如果形状函数要准确表示正常模态形状，则 $p$ 不得小于二。因此，FAST_AD 设计代码要求 $p$ 等于二。FAST_AD 允许 $N$ 为高达五。在 FAST_AD 中，与每个形状函数和正常模态相关联的比例系数是输入文件中请求的参数。一个名为 Modes 的预处理代码使 FAST_AD 用户能够获得这些参数。Modes 所采用的理论也在此开发。
-
-
 Using Lagrange’s equations for a conservative, scleronomic system13, the equations of motion for an $N\!\cdot$ -DOF system are equivalent to:  
+使用拉格朗日方程，对于一个保守、刚体系统，N-DOF系统运动方程等同于：
 
 $$
-\sum_{j=p}^{N+p-I}m_{i j}\ddot{c}_{j}(t)+\sum_{j=p}^{N+p-I}k_{i j}c_{j}\big(t\big)=O\quad\big(i=p,p+I,...,N+p-I\big)
+\sum_{j=p}^{N+p-1}m_{i j}\ddot{c}_{j}(t)+\sum_{j=p}^{N+p-1}k_{i j}c_{j}\big(t\big)=O\quad\big(i=p,p+1,...,N+p-1\big)
 $$  
 
 where the generalized mass and stiffness, $m_{i j}$ and $k_{i j}$ respectively, are defined in terms of the kinetic energy, $T,$ and potential energy, $V$ :  
-
+其中，广义质量和刚度，分别表示为 $m_{i j}$ 和 $k_{i j}$，由动能 $T$ 和势能 $V$ 定义：
 $$
 T=\frac{I}{2}\sum_{i=p}^{N+p-l}\sum_{j=p}^{N+p-l}m_{i j}\dot{c}_{i}(t)\dot{c}_{j}(t)
 $$  
@@ -1157,21 +1136,24 @@ V=\frac{I}{2}\sum_{i=p}^{N+p-I}\sum_{j=p}^{N+p-I}k_{i j}c_{i}(t)c_{j}(t)
 $$  
 
 Now, when the flexible beam is vibrating at a specific natural mode, say $a=m$ , the following conditions result:  
-
+现在，当柔性梁以特定的固有振动模式振动，例如 $a=m$ 时，会产生以下条件：
 $$
-q_{a}(t)\!=\!\left\{\!\!\begin{array}{l l}{\!Q_{a}\,s i n(\omega_{a}t\!+\!\psi_{a})}&{f\!o r\,a=m}\\ {\!0}&{\!o t h e r w i s e}\end{array}\!\!\right.
+q_{a}(t)\!=\!\left\{\begin{array}{l l}{Q_{a}\,s i n(\omega_{a}t\!+\!\psi_{a})}&{f\!o r\,a=m}\\ {0}&{\!o t h e r w i s e}\end{array}\!\!\right.
 $$  
 
 $$
-\begin{array}{r l}{c_{b}(t)\!=\!C_{m,b}q_{m}(t)}&{{}\big(b=p,p+l,...,N+p-l\big)}\end{array}
+\begin{array}{r l}{c_{b}(t)\!=\!C_{m,b}q_{m}(t)}&{{}\big(b=p,p+1,...,N+p-1\big)}\end{array}
 $$  
 
 where $Q_{a}$ is the amplitude of the deflection of the tip of the flexible beam associated with natural mode $a$ .  
 
 Substituting Eq. (3.45) into Eq. (3.46), then substituting the resulting equation into Eq. (3.42), results in (the subscript has been dropped from the specific natural mode):  
 
+其中，$Q_{a}$ 为与自然模式 $a$ 相关的柔性梁尖端挠曲振幅。
+
+将公式 (3.45) 代入公式 (3.46)，然后将所得方程代入公式 (3.42)，可得（此处已省略特定自然模式的下标）：
 $$
-\sum_{j=p}^{N+p-l}(-\omega^{2}m_{i j}+k_{i j})C_{j}=O\quad\displaystyle\big(i=p,p+l,...,N+p-l\big)
+\sum_{j=p}^{N+p-1}(-\omega^{2}m_{i j}+k_{i j})C_{j}=O\quad\displaystyle\big(i=p,p+l,...,N+p-l\big)
 $$  
 
 Equation (3.47) can be written in matrix form:  
@@ -1182,7 +1164,10 @@ $$
 
 where the generalized mass matrix, $[M]$ , and generalized stiffness matrix, $[K]$ , are both $N\times N$ matrices and the coefficient vector, $\{C\}$ , is an $N\times1$ vector.  The determinant of the matrix premultiplying the coefficient vector results in an $N^{\mathrm{th}}$ -degree algebraic equation in $\omega^{2}$ , which is called the characteristic equation.  The $N$ roots, ${\omega_{a}}^{2}$ , are the eigenvalues, each being the square of the natural frequency associated with normal mode $a$ .  The eigenvector associated with each eigenvalue, $\{C\}_{a}$ , defines the constant proportionality coefficients associated with normal mode $a$ (the $C_{a,b}$ ’s).  
 
+其中，广义质量矩阵 $[M]$ 和广义刚度矩阵 $[K]$ 均为 $N\times N$ 阶矩阵，系数向量 $\{C\}$ 是一个 $N\times1$ 阶向量。系数向量前乘的矩阵的行列式会得到一个关于 $\omega^{2}$ 的 $N$ 次代数方程，称为特征方程。该方程的 $N$ 个根，${\omega_{a}}^{2}$ ，是特征值，每个特征值都是与模态 $a$ 相关的固有频率的平方。与每个特征值相关的特征向量，$\{C\}_{a}$ ，定义了与模态 $a$ 相关的常数比例系数（即 $C_{a,b}$）。
+
 In FAST_AD and Modes, the tower is modeled as an inverted cantilever beam with a point mass affixed to its free end.  The point mass, $M_{T o p}.$ , represents the combined mass of the base plate, nacelle, hub, and blades.  The tower is assumed to deflect in the longitudinal and lateral directions independently.  The stiffness distributions in each direction are assumed to be identical; consequently, the associated natural mode shapes and frequencies are assumed to be identical in each direction.  
+在 FAST_AD 和 Modes 中，塔架被建模为一个倒置悬臂梁，其自由端附有一个质点。该质点，$M_{T o p}.$，代表了底板、舱壳、轮毂和叶片的联合质量。塔架被假定在纵向和横向独立挠曲。每个方向的刚度分布被假定相同；因此，相关的固有振型和频率也被假定在每个方向上相同。
 
 The kinetic energy of the tower has a component associated with the distributed mass of the beam and a component associated with the point mass:  
 

多体+耦合求解器/V0.5版本结果对比/对比步骤.md

@@ -15,3 +15,23 @@ q 叶片、塔架的振型如何看？
 2、V0.5软件与fast 15mw对比
 
 3、气动改用AeroDyn ED对比
+
+
+V4.0 vs V3.5 新增
+.fst
+- seast
+ElastoDyn
+- 增加ptfmxyIner
+- 增加YAW-FRICTION
+
+AD 15
+-  删去 AFAeroMod
+-  删去 FrozenWake
+-  增加NacelleDrag
+- 增加BEM_Mod
+
+潜在问题：
+1 刚度、模态曲线是否一致
+   bladed 模态振型导出
+   bladed 模态与fast一致后，计算结果如何
+2 气动结果差异

多体+耦合求解器/rust 经验.md

@@ -19,3 +19,6 @@ aug_mat矩阵大小600
     ];
 ```
 ![[Pasted image 20250123103848.png]]
+
+
+# 创建lib

工作OKRs/25.2-5 OKR.canvas

@@ -1,11 +1,11 @@
 {
 	"nodes":[
 		{"id":"b698e88ca5fb9c51","type":"text","text":"状态指标：\n推进OKR的时候也要关注这些事情，它们是完成OKR的保障。\n\n\n效率状态  green","x":-96,"y":80,"width":456,"height":347},
-		{"id":"2b068bfe5df15a72","type":"text","text":"# 目标：多体动力学模块完善\n### 每周盘点一下它们\n\n\n关键结果：建模原理、建模方法掌握    （8/10）\n\n关键结果：对标Bladed模块完成  （8.5/10）\n\n关键结果：风机多体动力学文献调研情况完成 （5.5/10）","x":-96,"y":-307,"width":456,"height":347},
+		{"id":"2b068bfe5df15a72","type":"text","text":"# 目标：多体动力学模块完善\n### 每周盘点一下它们\n\n\n关键结果：建模原理、建模方法掌握    （8.5/10）\n\n关键结果：对标Bladed模块完成  （9/10）\n\n关键结果：风机多体动力学文献调研情况完成 （5.5/10）","x":-96,"y":-307,"width":456,"height":347},
 		{"id":"01ee5c157d0deeae","type":"text","text":"# 推进计划\n未来四周计划推进的重要事情\n\n文献调研启动\n\n建模重新推导\n\n\n","x":-620,"y":80,"width":456,"height":347},
-		{"id":"cb9319d24c3e70e3","type":"text","text":"# 5月已完成\n\nP1 Steady Operational Loads求解器测试 \n- 变桨算法测试完成\n- 转速算法基本完成\n- 两个结合点测试 完成\n\nP1 simpack多体对CAE的需求梳理 分成塔架、叶片、传动链模型 完成","x":-240,"y":520,"width":440,"height":560},
+		{"id":"cb9319d24c3e70e3","type":"text","text":"# 5月已完成\n\nP1 Steady Operational Loads求解器编写测试 \n- 变桨算法测试完成\n- 转速算法基本完成\n- 两个结合点测试 完成\n\nP1 Steady Parked Loads求解器编写及测试\n\nP1 simpack多体对CAE的需求梳理 分成塔架、叶片、传动链模型 完成\n\nP1 建立IEA 15yaml文件 完成\nP1 结果对比\n-  完成 bladed、fast模型建立，工况设置，对比","x":-240,"y":520,"width":440,"height":560},
 		{"id":"1ebeabaf5c73ddbb","type":"text","text":"# 4月已完成\n\n多体原理学习 YouTube课程 018\n\n气动模块联合调试，跑通\n\n使用python搭建风电机组多体模型  刚性部件主动力、惯性力计算 \n\n编写Steady Operational Loads求解器\n- 稳态工况多体动力学求解方法 --龙格库塔+ed_caloutput 初步方案完成\n- 遍历风速框架 完成\n- 不同风速下转速、变桨角度算法 完成\n- 多体设置参数 完成\n- 每个风速直接是否需要重新初始化 需要 完成\n\n\n\ngenerator torque计算 简单了解，确定方案","x":-720,"y":520,"width":440,"height":560},
-		{"id":"58be7961ae7275a7","type":"text","text":"# 计划\n这周要做的3~5件重要的事情，这些事情能有效推进实现OKR。\n\nP1 必须做。P2 应该做\n\n\nP1 柔性部件 叶片、塔架主动力惯性力算法 主线\n- 变形体动力学 简略看看ing\n- 浮动坐标系方法 如何用于梁模型 \n\t- Q 问孟航  不用浮动坐标系\n- 柔性梁弯曲变形振动学习，主线 \n\t- 广义质量 刚度矩阵及含义\n- 如何静力学求解\n\t- 基于本构方程 读孟的论文\n\t\nP1 建立IEA 15yaml文件\nP1 结合yaml解析代码，联合气动更新对yaml文件的支持\nP1 结果对比\n- 优先瞬态，**如何设置计算？**\n\nP1 如何优雅的存储、输出结果。\nP1 国产化适配交给甲子营，对接\nP1 推进气动、控制、多体、水动\nP2 yaw 自由度再bug确认 已知原理了\n\n","x":-614,"y":-307,"width":450,"height":347}
+		{"id":"58be7961ae7275a7","type":"text","text":"# 计划\n这周要做的3~5件重要的事情，这些事情能有效推进实现OKR。\n\nP1 必须做。P2 应该做\n\n\nP1 柔性部件 叶片、塔架主动力惯性力算法 主线\n- 变形体动力学 简略看看ing\n- 浮动坐标系方法 如何用于梁模型 \n\t- Q 问孟航  不用浮动坐标系\n- 柔性梁弯曲变形振动学习，主线 \n\t- 广义质量 刚度矩阵及含义\n- 如何静力学求解\n\t- 基于本构方程 读孟的论文\n\t\n\nP1 结合yaml解析代码，联合气动更新对yaml文件的支持 \nP1 结果对比\n- Herowind 不带气动与fast3.5对比\n- Herowind 不带气动与fast4.0对比\n- Herowind 带气动与fast对比\n\nP1 如何优雅的存储、输出结果。\nP1 国产化适配交给甲子营，对接\nP1 推进气动、控制、多体、水动 耦合计算\nP2 yaw 自由度再bug确认 已知原理了\n\n","x":-614,"y":-307,"width":450,"height":347}
 	],
 	"edges":[]
 }