瞬态：
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