1.2 KiB
瞬态： 0时刻 多体 变形 -> 气动 气动 载荷 -> 多体
t=1 多体 龙格库塔 -> 变形 -> 气动
初步构想：
沿用 多体 init -> 气动 -> 多体 龙格库塔 -> 变形
或者 气动刚体计算 载荷 -> 多体 多体 龙格库塔 -> 变形 -> 气动 气动 载荷 -> 多体 多体 变形 -> 气动 ...至收敛，作为稳态解
Bladed
结果方面：
叶片：叶片1载荷 Mx My Mxy Mz Fx Fy Fxy Fz Hub: 载荷 rotating GL coordinate Hub: 载荷 fixed GL coordinate yaw bearing : 载荷 GL coordinate Tower member : 载荷 local coordinate Nacelle Motion: fore-after\side-side disp / vel / acc Supporting structure : deflection / vel / acc
叶片没有motion，只有塔顶motion
blade tip to tower centre
问题
如何解决 转速设定、变桨角度设定
转速
multibody_input_file.rot_speed = 1.267109; pub fn ed_set_parameters p.rot_speed = multibody_input_file.rot_speed;
变桨角度
init _u
初始化阶段 pitch_com等于0状态下，设置坐标系
set_position后 给u.bl_pitch_com赋值
u.bl_pitch_com = multibody_input_file.bl_pitch.clone();
t_0时刻ed_calc_const_state_deriv函数中 重新计算set_coord_sy