圆形断面临界水深

% -- 计算  抛物线 形渠道的 临界水深计算
clc; clear all
beta=0.001:0.01:10;  % 无量纲参数的取值范围  β=αQ/b^2.5;
g=9.81; % 重力加速度
Mk=8*(1/g*beta.^2).^(1/3);
 
 
hk025=0.207*(Mk).^0.7704; % 单位 m；
plot(beta,hk025,'b-','linewidth',2); hold on
  
  axis([0.001,3,0.01,1])
 

xlabel('无量纲参数  Q/d^2^.^5 ');
 ylabel('无量纲临界水深 h_k /d');
  
% legend('地基加固后2#桩检测值','本文预估计算值（1#桩）', '原状土1#桩检测值','本文计算值（1#桩）') 
 
title('圆形渠道 临界水深查算图');