[实验流体力学][Matlab] 船模阻力实验

例1

% C_total_model = R_total_model/(1/2*rho_model*S_model*V_model^2)
% 式中：
% C_total_model 船模总阻力系数
% R_total_model 船模阻力 N
% rho_model 实验水池中的水的密度 kg*m^(-3)
% S_model 船模湿面积 m^2
% V_model 船模速度 m*s^(-1)

% C_friction_model = 0.075/(lg(Re_model)-2)^2
% 式中：
% C_friction_model 船模摩擦阻力系数
% Re_model = V_model*L_model/mu_model 船模雷诺数
% L_model 船模水线长 m
% mu_model 实验水池中的水的运动粘性系数 可从水温表查得

% C_rest_model = C_total_model-C_friction_model
% 式中：
% C_rm 船模的剩余阻力系数
% C_total_model 船模总阻力系数
% C_friction_model 船模摩擦阻力系数

% Fr = V^2/(g*L) 佛鲁德数
% 根据相似定律 Fr_ruler = Fr_model/Fr_simulation
% => V_simulation = V_model*lambda^(1/2)
% 式中：
% lambda = L_simulation/L_model
% 根据相似定律 C_rest_simulation = C_rest_model

% C_friction_simulation = 0.075/(lg(Re_simulation-2)^2)
% 式中：
% C_friction_simulation 实船摩擦阻力系数
% Re_simulation = V_simulation*L_simulation/mu_simulation 实船雷诺数
% L_simulation 实船水线长 m
% mu_simulation 实船的水的运动粘性系数 一般取 15℃ 时的值

% 通常取实船粗糙度补贴 Delta_C_friction = 0.4*10^(-3)

% C_total_simulation = C_rest_model + C_friction_simulation + Delta_C_friction
% 式中：
% C_total_simulation 实船总阻力系数

% R_total_simulation =
% 1/2*rho_simulation*S_simulation*V_simulation^2*C_total_simulation N
% 式中：
% R_total_simulation 实船总阻力
% rho_simulation 实船的水的密度 kg*m^(-3) 一般取 15℃ 时的值
% S_simulation 实船湿面积 m^2

% P_effective_simulation = P_total_simulation * V_simulation/1000
% 式中：
% P_effective_simulation 有效功率 kW
% P_total_simulation 总功率 W
% V_simulation 实船速度 m*s(-1)

clear;

% 设计数据

% 尺度比
lambda = 50;

% 实船水线长 m
L_simulation = 152;

% 实船湿面积 m^2
S_simulation = 4075;

% 实验数据

% 船模速度 m*s^(-1)
V_model = [
    0.7
    0.8
    0.9
    1
    1.1
    1.2
    1.2
    1.3
    1.4
    1.5
    1.6];

% 船模总阻力 N
R_total_model = [
    0.18042
    0.22504
    0.27963
    0.35303
    0.43209
    0.50693
    0.49687
    0.6238
    0.78914
    1.15075
    1.59732]/