Resolve normal depth from Manning's equation

2020 8 1
1 个回答
回答已采纳
更新时间：2020 8 1
2 次查看（30 天）

0 个投票

Hello,
I aim to obtain the normal depth of a channel using Mannig's equation. Somehow I don't manage to resolve its value. Here it's the pieco of code that I'm using:
riverSlope=0.0114; % [m/m] - inletSlope, outletSlope or riverSlope
bottom_width=33.5937; % [m] - inlet or outlet bottom width
slope_Rbank=1.1336; % [m/m] - slope_Rbank_in or slope_Rbank_out
slope_Lbank=0.3334; % [m/m] - slope_Lbank_in or slope_Lbank_out
q=10; % [m3/s] - Flow discharge
n=0.04; % [-] - Manning's roughness coefficient
syms y
area=(bottom_width+(y/(2*slope_Rbank))+(y/(2*slope_Lbank)))*y;
wetted_perimeter=bottom_width+y*(sqrt(1+(1/slope_Rbank)^2)+sqrt(1+(1/slope_Lbank)^2));
manning_eqn=@(y)(1/n)*((area/wetted_perimeter)^(2/3))*(riverSlope^(1/2))*area==q;
soly=solve(manning_eqn,y)
I would really appreciate if someone can help to fix it in order to obtain the desired values and avoid the coding of an iteration loop for the manual calculation. Thanks in advance!!
Álvaro

请先登录，再回答此问题。

 采纳的回答

Alan Stevens
Alan Stevens 2020-8-1
编辑：Alan Stevens 2020-8-1

0 个投票

This shoud do it:
depth0 = 1; % Initial guess
depth = fzero(@manningfn, depth0);
function manning = manningfn(y)
riverSlope=0.0114; % [m/m] - inletSlope, outletSlope or riverSlope
bottom_width=33.5937; % [m] - inlet or outlet bottom width
slope_Rbank=1.1336; % [m/m] - slope_Rbank_in or slope_Rbank_out
slope_Lbank=0.3334; % [m/m] - slope_Lbank_in or slope_Lbank_out
q=10; % [m3/s] - Flow discharge
n=0.04; % [-] - Manning's roughness coefficient
area=(bottom_width+(y/(2*slope_Rbank))+(y/(2*slope_Lbank)))*y;
wetted_perimeter=bottom_width+y*(sqrt(1+(1/slope_Rbank)^2)+sqrt(1+(1/slope_Lbank)^2));
manning = (1/n)*((area/wetted_perimeter)^(2/3))*(riverSlope^(1/2))*area-q;
end

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Álvaro Pardo
Álvaro Pardo 2020-8-1
Is there any way I can define the parameters out of the function? That would be nice because I aim to use them for other applications. Many thanks!!
Alan Stevens
Alan Stevens 2020-8-1
Yes, you could do this:
riverSlope=0.0114; % [m/m] - inletSlope, outletSlope or riverSlope
bottom_width=33.5937; % [m] - inlet or outlet bottom width
slope_Rbank=1.1336; % [m/m] - slope_Rbank_in or slope_Rbank_out
slope_Lbank=0.3334; % [m/m] - slope_Lbank_in or slope_Lbank_out
q=10; % [m3/s] - Flow discharge
n=0.04; % [-] - Manning's roughness coefficient
data =[riverSlope; bottom_width; slope_Rbank; slope_Lbank; q; n];
depth0 = 1; % Initial guess
depth = fzero(@manningfn, depth0,[],data);
function manning = manningfn(y, data)
riverSlope=data(1);
bottom_width=data(2);
slope_Rbank=data(3);
slope_Lbank=data(4);
q=data(5);
n=data(6);
area=(bottom_width+(y/(2*slope_Rbank))+(y/(2*slope_Lbank)))*y;
wetted_perimeter=bottom_width+y*(sqrt(1+(1/slope_Rbank)^2)+sqrt(1+(1/slope_Lbank)^2));
manning = (1/n)*((area/wetted_perimeter)^(2/3))*(riverSlope^(1/2))*area-q;
end

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Mathematics and Optimization 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by