空化器阻力系数

  • 函数：DragCone(beta)
  • 计算公式：

  1. 将角度 beta 转换为弧度：

    μ=β180\mu = \frac{\beta}{180}

    其中，μ\mu 是 锥形空化器的锥角beta 转换为弧度的值。

  2. 根据 beta 的不同范围计算 Cx0：

    • 如果 0≤β<30∘0 \leq \beta < 30^\circ：

      Cx0=μ2(0.915+9.52μ)C_{x0} = \frac{\mu}{2} \left( 0.915 + \frac{9.5}{2} \mu \right)

    • 如果 30∘≤β≤180∘30^\circ \leq \beta \leq 180^\circ：

      Cx0=0.5+1.81(μ2−0.25)−2(μ2−0.25)2C_{x0} = 0.5 + 1.81 \left( \frac{\mu}{2} - 0.25 \right) - 2 \left( \frac{\mu}{2} - 0.25 \right)^2

    • 如果 β>180∘\beta > 180^\circ：

      Cx0=0.8275+0.00225⋅(μ⋅180−180)−0.000007176⋅(μ⋅180−180)2C_{x0} = 0.8275 + 0.00225 \cdot (\mu \cdot 180 - 180) - 0.000007176 \cdot (\mu \cdot 180 - 180)^2

  • 代码
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
def DragCone(beta):
    # 计算空化器的阻力系数
    # 输入为
    # Beta: 锥形空化器的锥角，单位为deg
    # 输出为Cx0
    # Cx0: 空化器的阻力系数（无量纲）
    # beta角度转弧度
    miu = beta / 180
    # 根据beta的范围计算阻力系数
    if beta < 30:
        # 阻力系数
        cx0 = miu / 2 * (0.915 + 9.5 / 2 * miu)
    elif 30 <= beta <= 180:
        cx0 = 0.5 + 1.81 * (miu / 2 - 0.25) - 2 * (miu / 2 - 0.25)**2
    else:
        cx0 = 0.8275 + 0.00225 * (miu * 180 - 180) - 0.000007176 * (miu * 180 - 180)**2
    return cx0