Ajout de la calculette "passe à enrochement simple"

Cette calculette correspond au module de précalcul d'une passe à enrochement développée par Ludovic Cassan.

Il s'agit d'une passe à enrochement sans déversoir latéral (la totalité du débit transite de l'amont vers l'aval de la passe sans dérivation).

Les paramètres de l'équation sont les suivants :

• La côte amont (m) (n'est utilisée que pour calculer la côte du bas de la rampe à partir de la longueur et de la pente);
• La longueur de la rampe (m);
• La rugosité de fond (m);
• Trois paramètres de blocs :
  • Diamètre (m);
  • Hauteur (m);
  • Forme Cd (1 pour la forme ronde, 2 pour la forme carrée).
• Pente (m/m);
• Largeur (m);
• Profondeur (m);
• Concentration (-);
• Débit (m3/s);

La version présentée par L. Cassan permet de calculer certains paramètres : La pente, la largeur, la profondeur, la concentration et le débit.

Elle permet de faire varier ces mêmes paramètres et dans ce cas, permet de calculer la profondeur ou alors le débit si c'est la profondeur qui varie.

Des résultats complémentaires sont calculés :

• La cote de bas de rampe (m);
• La vitesse débitante (m/s);
• La vitesse max (m/s);
• La vitesse max guide technique (m/s);
• P/V (W/m3);
• Le Froude (-);
• Le type d'écoulement (émergent, quasi-émergent, immergé);
• Le débit du guide technique (m3/s).

Les résultats graphiques proposés quand un paramètres varie sont les vitesses en fonction du paramètre qui varie et la variable calculée en fonction du paramètre qui varie.

Le code source de l'équation au format Matlab : find_Q_nat.m

TAF

Programmer la calculette de la passe à enrochement :

• Définir les paramètres de l'équation;
• Programmer l'équation à partir du code source matlab;
• Calculer les résultats complémentaires.

Le code source Matlab pour le calcul des résultats complémentaires est le suivant :

cote_bas=z_amont-S*long;
set(handles.edit11,'string',num2str(cote_bas))

Vdeb=Q/L/pf;
set(handles.edit12,'string',num2str(Vdeb))

Vg=Q/L/pf/(1-C^0.5);
Fr=Vg./(g.*pf).^0.5;
if Cd==2
    Vmax=Vg.*(min(1.6./(1-(Fr.^2)/4),Fr.^(-2/3)));%
else
    Vmax=Vg.*(min(1.2./(1-(Fr.^2)/4),Fr.^(-2/3)));
end

set(handles.edit16,'string',num2str(Vmax))

P=1000*g*Q/L*S;
set(handles.edit13,'string',num2str(P))


set(handles.edit14,'string',num2str(Fr))

if pf/h<1
    set(handles.edit15,'string','emergent')
elseif pf/h<1.1 & pf/h>=1
    set(handles.edit15,'string','quasi emergent')
elseif pf/h>1.2
    set(handles.edit15,'string','immerge')
end



if pf/h>1.1

    q_technique=0.955*(pf/h)^2.282*S^0.466*C^(-0.23)*(9.81*D)^0.5.*h*L;
else
    if Cd==2
        
        q_technique= 0.648*(pf/D)^1.084*S^0.56*C^(-0.456)*(9.81*D)^0.5.*D;
        V_technique=3.35*(pf/D)^0.27*S^0.53*(9.81*D)^0.5;
    else
        q_technique=0.815*(pf/D)^1.45*S^0.557*C^(-0.456)*(9.81*D)^0.5.*D*L;
        V_technique=4.54*(pf/D)^0.32*S^0.56*(9.81*D)^0.5;
        
    end
    
   set(handles.edit19,'string',num2str(V_technique)) ; 
end
set(handles.edit18,'string',num2str(q_technique)) ;

added specification and removed ~62 labels

changed the description

mentioned in issue nghyd#44 (closed)

changed title from Ajout de la calculette "passe à enrochement" to Ajout de la calculette "passe à enrochement simple"

changed the description

Pour le coefficient Cd proposer une liste déroulante pour la forme des plots (rond, carré) avec la possibilit de modifier le coefficient Cd sous la même forme qua Lechapt et Calmon

Etendre "Varier" à tous les paramètres. Voir pour les paramètres à calculer.

Voir pour fusionner avec la calculette #35 (closed)

removed assignee

removed ~63 label

changed milestone to %"Passes à macro-rugosités"

added To Do label

assigned to @david.dorchies

added Doing and removed To Do labels

created branch 22-ajout-de-la-calculette-passe-a-enrochement-simple

mentioned in commit d35f198d

mentioned in commit d418d0fb

mentioned in commit de7a95a7

mentioned in commit dbae67af

mentioned in commit 50c28cee

mentioned in commit 77512fbb

mentioned in commit 5403bfe0

mentioned in commit e6cabb2a

mentioned in commit 4e6145c3

mentioned in commit b90e0bfe

mentioned in commit d6b482ad

mentioned in commit 6f3cd0ab

Il y a une erreur dans le calcul du débit du guide technique dans le cas émergent.

r.extraResults.Q_GuideTech = cQ[0] * Math.pow(this.prms.Y.v / hdk, cQ[1]) *
            Math.pow(this.prms.If.v, cQ[2]) * Math.pow(this.prms.C.v, cQ[3]) *
            Math.sqrt(MacroRugo.g * this.prms.PBD.v) * this.prms.PBD.v * this.prms.B.v;
        if (this.prms.Y.v / this.prms.PBH.v <= MacroRugo.limitSubmerg) {
            r.extraResults.V_GuideTech = cV[0] * Math.pow(this.prms.Y.v / this.prms.PBD.v, cV[1]) *
                Math.pow(this.prms.If.v, cQ[2]) * Math.sqrt(MacroRugo.g * this.prms.PBD.v);
        }

Il faut remplacer tous les this.prms.PBD.v par hdk.

mentioned in commit c6f3597f

mentioned in commit 8eb14f35

mentioned in commit 89f6e2fe

mentioned in commit 9d3abff3

mentioned in commit 76b5631a