vault backup: 2025-02-10 08:53:02

.obsidian/core-plugins.json

@@ -26,5 +26,6 @@
   "workspaces": false,
   "file-recovery": true,
   "publish": false,
-  "sync": false
+  "sync": false,
+  "webviewer": false
 }

多体求解器debug/多体+气动 yaw debug.md

@@ -5,6 +5,9 @@
 
 ## yaw求解过程
 
+## 思路
+yaw_mom 到 augmat - 再到 qdt dof_yaw
+哪里发生了比较大的变化导致 发散  外推？drtr？
 
 1 输入yaw_mom from 控制
 
@@ -27,10 +30,76 @@ Partial moment at the base plate (body B) / yaw bearing (point O) due the nacell
 pmom_ngn_rt 就有问题
 Partial moment at the nacelle (body N) / selected point on rotor-furl axis (point V) due the structure that furls with the rotor, generator, and rotor [-]
 
-pmom_lprot 有问题 （2，10）
+pmom_lprot 有问题 （2，10）PMomLPRot-（3，11）
+pmom_lprot - Partial moment at the teeter pin (point P) on the low-speed shaft (body L) due to the rotor [-]
+pmom_h0b - Partial moment at the hub (body H) / blade root (point S(0)) due to the blade [-]
+pmom_lprot[.., dof] += rt_hs.pmom_h0b[.., K, dof] + tmp_vec
+
+n = 0
+rust
+
+| pmom_lprot[.., 10]  | [10065.169025312112, -880.58815982021201, 0]                     |
+| ------------------- | ---------------------------------------------------------------- |
+| tmp_vec             | [-106008.52186436673, 4628.43188287714, -1.0682859215676195e-10] |
+| pmom_h0b[.., 0, 10] | [-1273126.2724884835, 55585.89184040811, -2.254358687222148e-9]  |
+| pmom_lprot[.., 10]  | [-1369069.625327538, 59333.735563465045, -2.36118727937891e-9]   |
+|                     |                                                                  |
+| tmp_vec             | [89979.12573566216, -452011.2188342586, -248942.84544012396]     |
+| pmom_h0b[.., 1, 10] | [1422330.1493642807, -9454670.27374079, -5253335.786885787]      |
+| pmom_lprot[.., 10]  | [143239.64977240493, -9847347.757011585, -5502278.632325914]     |
+|                     |                                                                  |
+| tmp_vec             | [89979.12573566235, -452011.21883425827, 248942.8454401244]      |
+| pmom_h0b[.., 2, 10] | [1422330.149364284, -9454670.27374078, 5253335.786885795]        |
+| pmom_lprot[.., 10]  | [1655548.9248723513, -19754029.249586623, 5.820766091346741e-9]  |
+
+fortran 
+
+| pmom_lprot[.., 10]  | [10065.17, -880.5881, 0.0]              |
+| ------------------- | --------------------------------------- |
+| tmp_vec             | [-106008.5, 4628.433, -1.0682862E-10]   |
+| pmom_h0b[.., 0, 10] | [-1273126., 55585.91, -2.2543594E-09]   |
+| pmom_lprot[.., 10]  | [-1369070., 59333.75, -2.3611879E-09]   |
+|                     |                                         |
+| tmp_vec             | [89979.14, -452011.3, -248942.9]        |
+| pmom_h0b[.., 1, 10] | [1422330., -9454672., -5253337.]        |
+| pmom_lprot[.., 10]  | [143239.8, -9847349., -5502280.]        |
+|                     |                                         |
+| tmp_vec             | [89979.14, -452011.3, 248942.9]         |
+| pmom_h0b[.., 1, 10] | [1422330., -9454672., 5253337.]         |
+| pmom_lprot[.., 10]  | [1655549., -1.9754032E+07., -0.1093750] |
 
 
 
 pfrc_prot 有问题
-plin_vel_ec 10, 0, 0有问题
+pfrc_prot - Partial force at the teeter pin (point P) due to the rotor [-]
-Partial linear velocity (and its 1st time derivative) of the hub center of mass (point C) in the inertia frame (body E for earth) [-]
+pfrc_s0b - Partial force at the blade root (point S(0)) due to the blade [-]
+pfrc_prot[.., dof] += rt_hs.pfrc_s0b[.., K, dof]
+
+
+## 切换双精度后：依然出现
+
+plin_vel_ec 10, 0, 0有问题 
+Partial linear velocity (and its 1st time derivative) of the hub center of mass (point C) in the inertia frame (body E for earth) [-] 
+
+rt_hs.pmom_lprot 2, 10 差两倍
+
+```fortran
+RtHSdat%PMomLPRot(:,p%DOFs%PCE(I)) = TmpVec2 - p%Hubg1Iner*CoordSys%g1*DOT_PRODUCT( CoordSys%g1, RtHSdat%PAngVelEH(p%DOFs%PCE(I),0,:) ) & - p%Hubg2Iner*CoordSys%g2*DOT_PRODUCT( CoordSys%g2, RtHSdat%PAngVelEH(p%DOFs%PCE(I),0,:) )
+```
+
+g1 [2]差2倍
+
+MomBNcRtt有差距，导致右端项算的不对 
+
+frc_vgn_rtt
+mom_ngn_rtt
+
+frc_prot_t - Portion of the force at the teeter pin (point P) due to the rotor associated with everything but the QD2T()'s [-]
+
+frc_s0bt - Force at the blade root (point S(0)) due to the blade (excluding QD2T components) [-]
+
+
+## 问题
+
+- Fortran 与rust加和 趋于0的结果 精度不同 Fortran e-5 rust e-11
+- Fortran写到txt文本的数值，读进rust有少许差异

杂项/中国移动会员.md

@@ -0,0 +1,8 @@
+plus会员 ：
+兑换月卡，价值较高  优酷爱奇艺腾讯+咪咕会员
+电费立减 百度网盘月卡 加油券 腾讯体育vip月卡  饿了吗季卡
+
+黄金会员
+爱奇艺日卡 咪咕月卡 优酷周卡 美团饿了么5元红包 音乐7天
+
+定向流量