旋筒风帆的升阻力系数、升阻力力矩系数以及力臂的分析

旋筒帆的升阻力系数、升阻力力矩系数以及力臂系数的数值模拟结果如下：

以上各个物理量的无因次化公式如下：

进行多元非线性回归的表达式如下：

升力系数：

阻力系数：

升力力矩系数：

阻力力矩系数;

基于以上数据对旋筒帆作用于船舶上的横摇力矩以及纵摇力矩进行计算，基于fortran语言编程代码如下：

!-------------------------------------------------------------------
module rotorsail_forces
   
    implicit none
    ! Define parameters
    real(8) :: fn, rlpp
    real(8) :: sail_area, wind_speed, wind_angle,wind_angle0, vel
    real(8) :: cl, cd, Ml, Md, arm_X, arm_Y
    real(8), par