gyz/obsidian_backup
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

vault backup: 2026-01-09 16:17:17

Browse Source
This commit is contained in:
yize 2026-01-09 16:17:19 +08:00
parent 9586af29ed
commit de662f9877
5 changed files with 26 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

35
多体+耦合求解器/多体理论框架及ED求解方法.canvas
View File

@ -41,18 +41,18 @@
{"id":"cccd4482e817fca4","type":"text","text":"$$\n{k'}_{ij}^{B1F} = \\sqrt{FlStTunr^{B1}(i) \\, FlStTunr^{B1}(j)} \\int_{0}^{BldFlexL} EI^{B1F}(r) \\frac{d^2 \\phi_i^{B1F}(r)}{dr^2} \\frac{d^2 \\phi_j^{B1F}(r)}{dr^2} dr \\quad (i, j = 1, 2) \\tag{9}\n$$","x":-180,"y":1855,"width":980,"height":130},
{"id":"722e4b9c310e6d2a","type":"text","text":"$$\nEI^{B1F}(r) = AdjFlSt^{B1} \\cdot FlpStff^{B1}(r) \\tag{10}\n$$","x":-180,"y":2050,"width":980,"height":80},
{"id":"50e241bef2144ea1","type":"text","text":"$$\n\\frac{d^2 \\phi_i^{B1F}(r)}{dr^2} \\frac{d^2 \\phi_j^{B1F}(r)}{dr^2}\n$$","x":200,"y":2190,"width":220,"height":120},
{"id":"f65edb48b74e18d7","type":"text","text":"叶片广义挥舞刚度","x":-540,"y":1904,"width":250,"height":32},
{"id":"f65edb48b74e18d7","type":"text","text":"叶片广义挥舞刚度","x":-540,"y":1904,"width":250,"height":32,"color":"2"},
{"id":"e0d016f4bb106554","type":"text","text":"- 需要叶片挥舞刚度","x":-540,"y":2072,"width":250,"height":36},
{"id":"a9889437875a9263","type":"text","text":"- 模态形状二阶导","x":-540,"y":2232,"width":250,"height":36},
{"id":"da8b0f86919e99d9","type":"text","text":"$$\n\\left[ \\frac{d^2 \\phi_{1}^{B1E}(r)}{dr^2} \\right]^2\n$$","x":10,"y":2710,"width":220,"height":120},
{"id":"6454ea0b5f37aafa","type":"text","text":"- 模态形状二阶导","x":-540,"y":2752,"width":250,"height":36},
{"id":"36c9ab3fd7a0a402","type":"text","text":"- 需要叶片挥舞刚度","x":-540,"y":2592,"width":250,"height":36},
{"id":"ab5c75b7401fb031","type":"text","text":"叶片广义摆振刚度","x":-540,"y":2424,"width":250,"height":32},
{"id":"ab5c75b7401fb031","type":"text","text":"叶片广义摆振刚度","x":-540,"y":2424,"width":250,"height":32,"color":"2"},
{"id":"da00f9548457388f","type":"text","text":"$$\n{k'}_{11}^{B1E} = \\int_{0}^{BldFlexL} EI^{B1E}(r) \\left[ \\frac{d^2 \\phi_{1}^{B1E}(r)}{dr^2} \\right]^2 dr \\tag{11}\n$$","x":-180,"y":2380,"width":600,"height":120},
{"id":"a20bd9b934093730","type":"text","text":"$$\nEI^{B1E}(r) = AdjEdSt^{B1} \\cdot EdgStff^{B1}(r) \\tag{12}\n$$","x":-180,"y":2570,"width":560,"height":80},
{"id":"c477cd2a76cd46fb","x":-540,"y":3101,"width":250,"height":60,"type":"text","text":"- 模态形状一阶导"},
{"id":"2ee598faf0603188","x":-540,"y":2921,"width":250,"height":83,"type":"text","text":"离心刚化效应引起的刚度项"},
{"id":"dcd9236bd7c84bd4","type":"text","text":"叶片质量属性","x":-540,"y":3356,"width":250,"height":32},
{"id":"2ee598faf0603188","x":-540,"y":2921,"width":250,"height":83,"color":"2","type":"text","text":"离心刚化效应引起的刚度项"},
{"id":"dcd9236bd7c84bd4","type":"text","text":"叶片质量属性","x":-540,"y":3356,"width":250,"height":32,"color":"2"},
{"id":"ebc5c6d7451acc3b","type":"text","text":"- 叶片质量","x":-540,"y":3454,"width":250,"height":36},
{"id":"a49465d6a6a3c3c2","type":"text","text":"- 叶片一阶矩","x":-540,"y":3534,"width":250,"height":36},
{"id":"11e198e842aadec4","type":"text","text":"- 叶片质心","x":-540,"y":3853,"width":250,"height":36},
@ -64,22 +64,28 @@
{"id":"253044dc6cebe206","type":"text","text":"```rust\np.bld_cg[k] = p.first_mom[k] / p.bld_mass[k];\n\np.rot_mass += p.bld_mass[k];\n\np.rot_iner += (p.second_mom[k] + p.bld_mass[k] * p.hub_rad * (2.0 * p.bld_cg[k] + p.hub_rad)) * p.cos_pre_c[k].powi(2);\n```\n","x":-156,"y":3853,"width":796,"height":237},
{"id":"9a824939bad661d2","x":-2120,"y":2252,"width":250,"height":36,"type":"text","text":"- 摆振广义质量"},
{"id":"8fa04043b651b57e","x":-2120,"y":2394,"width":250,"height":36,"type":"text","text":"- 模态形状0阶导"},
{"id":"34fed438f96dfa4c","x":-2120,"y":2540,"width":250,"height":60,"type":"text","text":"固有频率"},
{"id":"34fed438f96dfa4c","x":-2120,"y":2540,"width":250,"height":60,"color":"2","type":"text","text":"固有频率"},
{"id":"9005abbfdd8f3246","x":-2120,"y":2660,"width":250,"height":60,"type":"text","text":"挥舞固有频率"},
{"id":"8944328366cdad0c","x":-2120,"y":2780,"width":250,"height":69,"type":"text","text":"考虑叶尖质量、离心刚化"},
{"id":"a96110788e264552","x":-2120,"y":3180,"width":250,"height":60,"type":"text","text":"摆振固有频率"},
{"id":"6f03c475fcf16a01","type":"text","text":"考虑叶尖质量、离心刚化","x":-2120,"y":3320,"width":250,"height":69},
{"id":"f1cd43613a1fdfa5","type":"text","text":"广义质量","x":-2120,"y":1996,"width":250,"height":32},
{"id":"f1cd43613a1fdfa5","type":"text","text":"广义质量","x":-2120,"y":1996,"width":250,"height":32,"color":"2"},
{"id":"5b82d2f5eb26eefc","type":"text","text":"- 挥舞广义质量","x":-2120,"y":2094,"width":250,"height":36},
{"id":"d89b6477389b9a45","x":-2856,"y":2214,"width":476,"height":112,"type":"text","text":"$$\nm_{11}^{\\prime B1E} = \\int_{0}^{BldFlexL} \\mu^{B1}(r) \\left[ \\phi_1^{B1E}(r) \\right]^2 dr \\tag{17}\n$$"},
{"id":"5d4daed7f83a3fa5","x":-2938,"y":2057,"width":640,"height":110,"type":"text","text":"$$\nm_{ij}^{\\prime B1F} = \\int_{0}^{BldFlexL} \\mu^{B1}(r) \\phi_i^{B1F}(r) \\phi_j^{B1F}(r) \\, dr \\quad (i, j = 1, 2) \\tag{16}\n$$"},
{"id":"44300ab8ae0139cc","x":-2618,"y":2630,"width":320,"height":120,"type":"text","text":"$$\nf_i^{\\prime B1F} = \\frac{1}{2\\pi} \\sqrt{\\frac{k_{ii}^{\\prime B1F}}{m_{ii}^{\\prime B1F}}} \\tag{14}\n$$"},
{"id":"0bc6d2d448f012fd","x":-3014,"y":2780,"width":716,"height":315,"type":"text","text":"```rust\n// Natural blade I-flap frequency w/o centrifugal stiffening nor     tip mass effects\n\np.freq_bf[[k, i, 0]] = INV2PI * (p.kbf[[k, i, i]] / (m_bf[[k, i, i]] - p.tip_mass[k])).sqrt();\n\n// Natural blade I-flap frequency w/o centrifugal stiffening, but w/ tip mass effects\n\np.freq_bf[[k, i, 1]] = INV2PI * (p.kbf[[k, i, i]] / m_bf[[k, i, i]]).sqrt();\n\n// Natural blade I-flap frequency w/  centrifugal stiffening and     tip mass effects\n\np.freq_bf[[k, i, 2]] = INV2PI * ((p.kbf[[k, i, i]] + k_bf_cent[[k, i, i]]) / m_bf[[k, i, i]]).sqrt();\n```"},
{"id":"5047212eb42c9e69","type":"text","text":"$$\nf_i^{\\prime B1E} = \\frac{1}{2\\pi} \\sqrt{\\frac{k_{ii}^{\\prime B1E}}{m_{ii}^{\\prime B1E}}} \\tag{15}\n$$","x":-2618,"y":3150,"width":320,"height":120},
{"id":"f5ccbef2e5e250b9","type":"text","text":"```rust\n// Natural blade 1-edge frequency w/o centrifugal stiffening nor      tip mass effects\n\np.freq_be[[k, 0, 0]] = INV2PI * (p.kbe[[k, 0, 0]] / (m_be[[k, 0, 0]] - p.tip_mass[k])).sqrt();\n\n// Natural Blade 1-edge frequency w/o  centrifugal stiffening, but w/ tip mass effects\n\np.freq_be[[k, 0, 1]] = INV2PI * (p.kbe[[k, 0, 0]] / m_be[[k, 0, 0]]).sqrt();\n\n// Natural Blade 1-edge frequency w/  centrifugal stiffening and      tip mass effects\n\np.freq_be[[k, 0, 2]] = INV2PI * ((p.kbe[[k, 0, 0]] + k_be_cent[[k, 0, 0]]) / m_be[[k, 0, 0]]).sqrt();\n```","x":-3014,"y":3320,"width":716,"height":315},
{"id":"b5bee6a3e7c29d31","type":"text","text":"叶片广义参数计算","x":-1440,"y":2916,"width":250,"height":93},
{"id":"204d6c2f3ec03ac0","x":-2120,"y":3740,"width":250,"height":60,"type":"text","text":"广义阻尼"},
{"id":"1684e3f8b6f17829","type":"text","text":"广义阻尼","x":-2120,"y":3860,"width":250,"height":60}
{"id":"b5bee6a3e7c29d31","type":"text","text":"**叶片广义参数计算**","x":-1440,"y":2916,"width":250,"height":93,"color":"1"},
{"id":"d89b6477389b9a45","x":-2774,"y":2196,"width":476,"height":112,"type":"text","text":"$$\nm_{11}^{\\prime B1E} = \\int_{0}^{BldFlexL} \\mu^{B1}(r) \\left[ \\phi_1^{B1E}(r) \\right]^2 dr \\tag{17}\n$$"},
{"id":"d1cb8d99b00c07c8","x":-3022,"y":3699,"width":724,"height":203,"type":"text","text":"$$\n\\left. F_r \\right|_{DampB1} =\n\\begin{cases}\n-\\frac{\\zeta_1^{B1F} {k'}_{11}^{B1F}}{\\pi f^{\\prime B 1 F}_1} \\dot{q}_{B1F1} - \\frac{\\zeta_2^{B1F} {k'}_{12}^{B1F}}{\\pi f^{\\prime B 1 F}_2} \\dot{q}_{B1F2} & \\text{for } r = B1F1 \\\\\n-\\frac{\\zeta_1^{B1E} {k'}_{11}^{B1E}}{\\pi f^{\\prime B 1 E}_1} \\dot{q}_{B1E1} & \\text{for } r = B1E1 \\\\\n-\\frac{\\zeta_1^{B1F} {k'}_{21}^{B1F}}{\\pi f^{\\prime B 1 F}_1} \\dot{q}_{B1F1} - \\frac{\\zeta_2^{B1F} {k'}_{22}^{B1F}}{\\pi f^{\\prime B 1 F}_2} \\dot{q}_{B1F2} & \\text{for } r = B1F2 \\\\\n0 & \\text{otherwise}\n\\end{cases} \\tag{13}\n$$"},
{"id":"204d6c2f3ec03ac0","x":-2120,"y":3771,"width":250,"height":60,"color":"2","type":"text","text":"广义阻尼"},
{"id":"1684e3f8b6f17829","type":"text","text":"挥舞和摆振方向结构阻尼比","x":-2120,"y":3980,"width":250,"height":69},
{"id":"b1dddd3ed142c348","x":-2600,"y":3973,"width":302,"height":84,"type":"text","text":"$\\zeta_{i}^{BIF} = BldFlDmp^{B1}(i)/100$\n$\\zeta_{i}^{BIE} = BldEdDmp^{B1}(i)/100$"},
{"id":"c328cec51da1530b","x":-2963,"y":4080,"width":665,"height":230,"type":"text","text":" ```fortran\n\tp%CBF(K,I,L) = ( 0.01*p%BldFDamp(K,L) )*p%KBF(K,I,L)/( Pi*p%FreqBF(K,L,1) )\n\tp%CBE(K,1,1) = ( 0.01*p%BldEDamp(K,1) )*p%KBE(K,1,1)/( Pi*p%FreqBE(K,1,1) )\n ``` "},
{"id":"169bd99d304f6f85","x":-2120,"y":4165,"width":250,"height":60,"type":"text","text":"广义阻尼系数"},
{"id":"cd865939df274ef7","x":-3291,"y":4420,"width":993,"height":623,"type":"text","text":" ```fortran\n ! Calculate the 2nd derivatives of the twisted shape functions:\n Shape = SHP( p%RNodesNorm(J), p%BldFlexL, p%BldFl1Sh(:,K), 2, ErrStat, ErrMsg )\n p%TwistedSF(K,1,1,J,2) = Shape*p%CThetaS(K,J)\n ...\n ! Integrate to find the 1st derivatives of the twisted shape functions:\n TwstdSF ( I,L, 1) = p%TwistedSF(K,I,L,J,2)*0.5*p%DRNodes(J)\n p%TwistedSF (K,I,L,J,1) = TwstdSF ( I,L, 1)\n ...\n ! Integrate to find the twisted shape functions themselves (i.e., their zeroeth derivative):\n TwstdSF ( I,L, 0) = p%TwistedSF(K,I,L,J,1)*0.5*p%DRNodes(J)\n p%TwistedSF (K,I,L,J,0) = TwstdSF ( I,L, 0)\n ...\n ! Integrate to find the blade axial reduction shape functions:\n AxRdBld ( I,L ) = 0.5*p%DRNodes(J)*( &\n p%TwistedSF(K,1,I,L,1)*p%TwistedSF(K,1,L,J,1) &\n + p%TwistedSF(K,2,I,L,1)*p%TwistedSF(K,2,L,J,1) )\n p%AxRedBld (K,I,L,J) = AxRdBld(I,L)\n ```"},
{"id":"43c52f1d3f4f4337","x":-2120,"y":4580,"width":250,"height":60,"color":"4","type":"text","text":"扭转形状函数和轴向缩短系数"}
],
"edges":[
{"id":"647c1b45edc92b02","fromNode":"9461f7dd96103316","fromSide":"bottom","toNode":"0c8534c8ba68c9a6","toSide":"top"},
@ -142,6 +148,13 @@
{"id":"63ad29914f63862d","fromNode":"34fed438f96dfa4c","fromSide":"right","toNode":"a96110788e264552","toSide":"right"},
{"id":"c58c7b2bd323e581","fromNode":"b5bee6a3e7c29d31","fromSide":"left","toNode":"f1cd43613a1fdfa5","toSide":"right"},
{"id":"b427c81d66172385","fromNode":"b5bee6a3e7c29d31","fromSide":"left","toNode":"34fed438f96dfa4c","toSide":"right"},
{"id":"eaf4a6506bae5f72","fromNode":"b5bee6a3e7c29d31","fromSide":"left","toNode":"204d6c2f3ec03ac0","toSide":"right"}
{"id":"eaf4a6506bae5f72","fromNode":"b5bee6a3e7c29d31","fromSide":"left","toNode":"204d6c2f3ec03ac0","toSide":"right"},
{"id":"ed232ba50bf6c970","fromNode":"1684e3f8b6f17829","fromSide":"left","toNode":"b1dddd3ed142c348","toSide":"right"},
{"id":"a0f82a31cbffe766","fromNode":"204d6c2f3ec03ac0","fromSide":"left","toNode":"d1cb8d99b00c07c8","toSide":"right"},
{"id":"0f934757be6b10ed","fromNode":"204d6c2f3ec03ac0","fromSide":"right","toNode":"1684e3f8b6f17829","toSide":"right"},
{"id":"eb49e3a35fec607d","fromNode":"204d6c2f3ec03ac0","fromSide":"right","toNode":"169bd99d304f6f85","toSide":"right"},
{"id":"4165dc7658a3b437","fromNode":"169bd99d304f6f85","fromSide":"left","toNode":"c328cec51da1530b","toSide":"right"},
{"id":"711a31af774553f9","fromNode":"43c52f1d3f4f4337","fromSide":"left","toNode":"cd865939df274ef7","toSide":"right"},
{"id":"55cbfdf0e50ce7f1","fromNode":"b5bee6a3e7c29d31","fromSide":"left","toNode":"43c52f1d3f4f4337","toSide":"right"}
]
}

2
工作OKRs/25.12-26.2 OKR.canvas
View File

@ -6,7 +6,7 @@
{"id":"82708a439812fdc7","type":"text","text":"# 1月已完成\n\n\n\n","x":-220,"y":134,"width":440,"height":560},
{"id":"505acb3e6b119076","type":"text","text":"# 12月已完成\n\nP1 明阳机型验证\n\n- 商业机型建模 done\n- 正常发电工况对故障工况支持 故障建模 done\n\t- 超速n4 普通超速 多体设fault结构体 \n\t- 卡桨、 控制设fault结构体\n\n- 故障工况检查 done\n- 批量计算更新配置文件，风文件，程序版本，再计算 done\n\nP1 稳态工况init_with_yaml检查 done\n\nP1 前端\n- 所有simulation功能测试及对接 done\n\nP1 演示ppt补充内容 再补充\n- 补充steady operational loads / steady parked loads 缺结果 done\n- 6个算例的跑通 找一个与bladed对比 缺结果 done \n- 内部控制器 pass\n- batch 并行计算 单个工况是否快 多工况并行 暂时不做\n\nP1 前端\n- steady输出的名字改掉 done\n- 批量计算 改成并行计算 功能界面\n- 瞬态计算更新控制器 测试 done\n- 简单内控\n- 比较Bladed与正常发电工况速度，总时间短一点 multicase done\n- 是否需要增加相对路径问题 done\n\nP1 dlc 71 done\n\nP1 输出量更新到Bladed相同命名 done\n\nP1 专利 done\n- 做出solidworks模型，写专利\n\n","x":-700,"y":134,"width":440,"height":560},
{"id":"30cb7486dc4e224c","type":"text","text":"# 2月已完成\n\n\n","x":260,"y":134,"width":440,"height":560},
{"id":"c18d25521d773705","type":"text","text":"# 计划\n这周要做的3~5件重要的事情，这些事情能有效推进实现OKR。\n\nP1 必须做。P2 应该做\n\n\nP2 柔性部件 叶片、塔架变形算法 主线\n- 变形体动力学 简略看看ing\n- 柔性梁弯曲变形振动学习，主线 \n\t- 广义质量 刚度矩阵及含义\n\t\n- 梳理bladed动力学框架\n\t- 子结构文献阅读\n\t- 叶片模型建模 done\n- 共旋方法学习\n- DTU 变形量计算方法学习\n\n\nP1 线性化方法编写 ing\n\n- 开始编写扰动代码\n- 形成系统矩阵-输出矩阵\n\n\nP1 明阳机型验证 \n\n- 商业机型模型验证 气动建模有问题\n- 动态失速计算慢 done\n- 叶片发散 叶片往前位移\n\n\nP1 服务器网络申请 ing\n- 等待自托管机架信息\n\n\nP1 报告更新到0.6版本 ing\n\nP2 bladed对比--稳态运行载荷，产出报告\n- 气动参与模块对比\n- 模态对比 两种描述方法不同，bladed方向更多，x y z deflection, x y z rotation，不好对比\n- 气动对比 aerodynamic info 轴向切向诱导因子，根部，尖部差距较大\n\nP2 yaw 自由度再bug确认 已知原理了\n","x":-597,"y":-803,"width":453,"height":457},
{"id":"c18d25521d773705","type":"text","text":"# 计划\n这周要做的3~5件重要的事情，这些事情能有效推进实现OKR。\n\nP1 必须做。P2 应该做\n\n\nP2 柔性部件 叶片、塔架变形算法 主线\n- 变形体动力学 简略看看ing\n- 柔性梁弯曲变形振动学习，主线 \n\t- 广义质量 刚度矩阵及含义\n\t\n- 梳理bladed动力学框架\n\t- 子结构文献阅读\n\t- 叶片模型建模 done\n- 共旋方法学习\n- DTU 变形量计算方法学习\n\n\nP1 线性化方法编写 ing\n\n- 开始编写扰动代码\n- 形成系统矩阵-输出矩阵\n\n\nP1 明阳机型验证 \n\n- 商业机型模型验证 气动建模有问题\n- 叶片发散 叶片往前位移 原因查找 done -- 刘\n\nP1 检查foundation fz\n\nP1 检查有没有塔筒上所有部件合并计算质量、质心、惯性参数\n\n\nP1 服务器网络申请 ing\n- 等待自托管机架信息\n- 机架有门\n- 电源 - 高\n- 散热\n\nP2 初始化方法，找到稳态结果开始\n\nP1 叶素、塔架station输出结果\n\n\nP1 报告更新到0.6版本 ing\n\nP2 bladed对比--稳态运行载荷，产出报告\n- 气动参与模块对比\n- 模态对比 两种描述方法不同，bladed方向更多，x y z deflection, x y z rotation，不好对比\n- 气动对比 aerodynamic info 轴向切向诱导因子，根部，尖部差距较大\n\nP2 yaw 自由度再bug确认 已知原理了\n","x":-597,"y":-803,"width":453,"height":457},
{"id":"86ab96a25a3bf82e","type":"text","text":" 湍流风+ 控制的联调，bladed也算一个算例\n- 加水动的联调\n- 8月份底完成这两个\n- 9月份完成停机等工况测试\n- 10月份明阳实际机型测试","x":580,"y":-803,"width":480,"height":220},
{"id":"e355f33c92cf18ea","type":"text","text":"9月份定常计算对接前端\n非定常测试完也对接前端","x":580,"y":-500,"width":480,"height":100},
{"id":"859e6853b7f1b92b","type":"text","text":"年底考核：\n专利\n线性化模块","x":1200,"y":-803,"width":320,"height":110},

BIN
工作总结/周报/周报109-郭翼泽.docx Normal file
View File

Binary file not shown.

3
未命名 1.base
View File

@ -1,3 +0,0 @@
views:
- type: table
name: 表格