Commit aa7bf93f authored by Enzo Vergnaud's avatar Enzo Vergnaud
Browse files

Update u4-fonctions.cpp

parent 695312ab
......@@ -87,23 +87,6 @@ void InitialiserDonnees()
gDonnees.corde.t_boum = 0;
//gDonnees.corde.Oke[0] = 0;
gInterface.SliderK1->value(gDonnees.corde.k1);
gInterface.SliderK2->value(gDonnees.corde.k2);
gInterface.SliderKc->value(gDonnees.corde.kc);
gInterface.SliderM->value(gDonnees.corde.m*10e3);
gInterface.SliderR1->value(gDonnees.corde.r1*10e2);
gInterface.SliderR2->value(gDonnees.corde.r2*10e2);
gInterface.SliderX3->value(gDonnees.corde.x[3]*10e4);
gInterface.SliderX2->value(gDonnees.corde.x[2]*10e4);
gInterface.SliderX1->value(gDonnees.corde.x[1]*10e4);
gInterface.SliderH0->value(gDonnees.corde.h[0]*10e2);
gInterface.SliderH1->value(gDonnees.corde.h[1]*10e3);
gInterface.SliderH2->value(gDonnees.corde.h[2]*10e3);
// Exemple son
//JouerSon("media/r2d2.mp3");
......@@ -361,7 +344,7 @@ double calculUg_Matlab()
if ((gDonnees.corde.h[1]<=0) || (gDonnees.corde.h[2]<=0)) //cas où on a fermeture de la glottique, donc abscence d'écoulement
{
Ug = 0; //on fixe alors le débit à zéro
dUg = (Ug-gDonnees.corde.Ugm1)*gDonnees.fe; //calcul de la dérivée du débit
//dUg = (Ug-gDonnees.corde.Ugm1)*gDonnees.fe; //calcul de la dérivée du débit
Delta=0; //par défaut on met Delta à Zéro, par NR
GA = 0;
}
......@@ -370,7 +353,7 @@ double calculUg_Matlab()
GA = 0.5*Rho/(gDonnees.corde.Lg*gDonnees.corde.Lg)*(1/(gDonnees.corde.hs*gDonnees.corde.hs)-1/(gDonnees.corde.h[0]*gDonnees.corde.h[0])); //terme de bernoulli
// if (h1~=h2) // cas ou la plaque 2 est non horizontale
printf("xs : %.10lf et hs : %.10lf et GA : %lf\n", gDonnees.corde.xs, gDonnees.corde.hs, GA);
//calcul du discriminant de l'équation du second degré écrit en débit
Delta = 4*GA*(DeltaPtotal);
if (Delta>=0) //>= au lieu de >
......@@ -383,8 +366,9 @@ double calculUg_Matlab()
Ug = 0 ;
}
//if (Ugi <0), Ugi = 0;end
dUg = (Ug-gDonnees.corde.Ugm1)*gDonnees.fe;
// dUg = (Ug-gDonnees.corde.Ugm1)*gDonnees.fe;
}
printf("xs : %.10lf et hs : %.10lf et GA : %lf\n", gDonnees.corde.xs, gDonnees.corde.hs, GA);
gDonnees.corde.Ugm1=Ug;
return Ug;
......@@ -440,6 +424,7 @@ double X_Matlab(int i) // terminé
else
{
X = 1/(pow(gDonnees.geometrie.A[i],2))*log(gDonnees.corde.h[i+1]/gDonnees.corde.h[i])+1/gDonnees.geometrie.A[i]*(gDonnees.corde.x[i]-gDonnees.corde.h[i]/gDonnees.geometrie.A[i])*(1/gDonnees.corde.h[i]-1/gDonnees.corde.h[i+1]);
}
return X;
......@@ -515,6 +500,7 @@ void calculForceMatlab()
// int i2=2;
//terme de bernoulli
bern = tmp1*(1/(pow(gDonnees.corde.h[0],2))+2*(X_Matlab(1)-gDonnees.corde.x[2]*W_Matlab(1))/(pow(gDonnees.corde.x[2]-gDonnees.corde.x[1],2)));
}
else //glotte formant un canal convergent ou divergent avec h2<hs
{
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment