# This file defines the const varibles

const ρ = 1.225 #空气密度 量纲kg*m^-3
const ν = 1.46e-4 #运动黏性系数 无量纲
const g = 9.8 #重力加速度  量纲kg*m*s^-2
const v_sound = 343.2 #音速 量纲m/s
const NR = 2  #桨叶片数
const Nbe = 10 #旋翼分段数量
const R = 1.05  #旋翼半径 量纲m
const ecut = 0.1 #桨叶根切比例 无量纲
const m_ = 0.35/R #桨叶质量密度 量纲kg/m
const Ω = 100  #旋翼转速 量纲rad/s
const αs = -2*π/180.0  #旋翼轴倾角  量纲rad
const Kβ = 0.0  #桨叶根部挥舞弹簧刚度 量纲？？？
const Iβ = m_/3*(1-ecut^3)*R^3  #桨叶挥舞惯量 量纲kg*m^2
# const R0 = 0.382  #lenght of root cut 量纲m
# const R1 = 0.136  #length of tip loss 量纲m
const vair = 50 # 来流速度 量纲 m/s
const v_air = [vair*cos(αs),0.0,vair*sin(αs)] #forward wind speed  量纲m/s
const μ_air = v_air[1]/(Ω*R)  #来流在桨盘edgewise方向分量  无量纲
const λ_air = v_air[3]/(Ω*R)  #来流在桨盘轴向分量  无量纲
const A = π*R^2 #参考面积 量纲为m^2
const fnonc = ρ*A*Ω*R^2*Ω #力的无量纲化参数 量纲kg*m*s^-2
const mnonc = ρ*A*Ω^2*R^3 #力矩的无量纲化参数 量纲kg*m^2/s^2
const T = 250 #飞行器重量 (量纲为kg*m*s^-2)
const CT = T/fnonc  #无量纲重量系数
const dpsideg = 10  # 方位角不仅长度（量纲为deg）
const dψ = dpsideg*π/180 #方位角步进步长 (量纲为rad)
const npsi = 360/dpsideg # 周向分割步数
const βang0 = 0.0 # 挥舞角初值
const dβ0 = 0.0 # 挥舞角速度初值
const θ7 = 2.0/180*π # 旋翼总距
const θtw = -5./180*π # 桨叶扭转角
const thelat = -0/180*π # 横向周期变距
const thelon = 0/180*π # 纵向周期变距

ch = Array{Float64}(NR,Nbe) #桨叶弦长
rb = Array{Float64}(NR,Nbe)  #桨叶叶素点分段径向长度值
dr = Array{Float64}(NR,Nbe)  #桨叶叶素分段长度值

for k in 1:NR
  for i in 1:Nbe
    ch[k,i] = 0.35
    dr[k,i] = R/Nbe
    rb[k,i] = R/Nbe*(i-1)+dr[k,i]/2
  end
end
