Commit c11cf7c1 authored by phelma's avatar phelma
Browse files

voilaa

parent b0aed1d5
......@@ -625,7 +625,7 @@
1330559674 /users/phelma/phelma2020/vergnaue/Documents/kitc-phelma/fltklinux64/FL/Fl_Tabs.H
"Fl_Group.H"
1615237772 source:/home/phelma/Documents/kitc-phelma/4projetcordesvocales/main.cpp
1615675130 source:/home/phelma/Documents/kitc-phelma/4projetcordesvocales/main.cpp
<iostream>
<FL/Fl.H>
<FL/fl_draw.H>
......@@ -645,13 +645,13 @@
1330559674 /home/phelma/Documents/kitc-phelma/fltklinux64/FL/fl_draw.H
"Enumerations.H"
1615237772 /home/phelma/Documents/kitc-phelma/4projetcordesvocales/main.h
1615675130 /home/phelma/Documents/kitc-phelma/4projetcordesvocales/main.h
<FL/Fl_Widget.H>
1330559674 /home/phelma/Documents/kitc-phelma/fltklinux64/FL/Fl_Widget.H
"Enumerations.H"
1615237772 /home/phelma/Documents/kitc-phelma/4projetcordesvocales/u1-interface.h
1615675130 /home/phelma/Documents/kitc-phelma/4projetcordesvocales/u1-interface.h
"main.h"
<FL/Fl_Double_Window.H>
<FL/Fl_Button.H>
......@@ -694,13 +694,13 @@
1330559674 /home/phelma/Documents/kitc-phelma/fltklinux64/FL/fl_ask.H
"Enumerations.H"
1615237772 /home/phelma/Documents/kitc-phelma/4projetcordesvocales/u3-callbacks.h
1615675130 /home/phelma/Documents/kitc-phelma/4projetcordesvocales/u3-callbacks.h
<FL/Fl_Widget.H>
1615492742 /home/phelma/Documents/kitc-phelma/4projetcordesvocales/u4-fonctions.h
1615684080 /home/phelma/Documents/kitc-phelma/4projetcordesvocales/u4-fonctions.h
<math.h>
1615237772 source:/home/phelma/Documents/kitc-phelma/4projetcordesvocales/u1-interface.cpp
1615675651 source:/home/phelma/Documents/kitc-phelma/4projetcordesvocales/u1-interface.cpp
<iostream>
<stdio.h>
"u1-interface.h"
......@@ -708,9 +708,9 @@
"u3-callbacks.h"
"u4-fonctions.h"
1615237772 /home/phelma/Documents/kitc-phelma/4projetcordesvocales/u2-dessin.h
1615675130 /home/phelma/Documents/kitc-phelma/4projetcordesvocales/u2-dessin.h
1615237772 source:/home/phelma/Documents/kitc-phelma/4projetcordesvocales/u2-dessin.cpp
1615675130 source:/home/phelma/Documents/kitc-phelma/4projetcordesvocales/u2-dessin.cpp
<iostream>
<FL/Fl.H>
<FL/fl_draw.H>
......@@ -718,7 +718,7 @@
"u2-dessin.h"
"u4-fonctions.h"
1615492046 source:/home/phelma/Documents/kitc-phelma/4projetcordesvocales/u3-callbacks.cpp
1615764528 source:/home/phelma/Documents/kitc-phelma/4projetcordesvocales/u3-callbacks.cpp
<iostream>
<FL/Fl.H>
<FL/fl_ask.H>
......@@ -813,7 +813,7 @@
1330559674 /home/phelma/Documents/kitc-phelma/fltklinux64/FL/Fl_Return_Button.H
"Fl_Button.H"
1615492176 source:/home/phelma/Documents/kitc-phelma/4projetcordesvocales/u4-fonctions.cpp
1615862892 source:/home/phelma/Documents/kitc-phelma/4projetcordesvocales/u4-fonctions.cpp
<iostream>
<stdlib.h>
<stdio.h>
......
<?xml version="1.0" encoding="UTF-8" standalone="yes" ?>
<CodeBlocks_layout_file>
<FileVersion major="1" minor="0" />
<ActiveTarget name="Debug" />
<File name="u3-callbacks.h" open="1" top="0" tabpos="4" split="0" active="1" splitpos="0" zoom_1="0" zoom_2="0">
<Cursor>
<Cursor1 position="672" topLine="0" />
</Cursor>
</File>
<File name="main.cpp" open="1" top="0" tabpos="9" split="0" active="1" splitpos="0" zoom_1="0" zoom_2="0">
<Cursor>
<Cursor1 position="1276" topLine="0" />
</Cursor>
</File>
<File name="u4-fonctions.cpp" open="1" top="1" tabpos="8" split="0" active="1" splitpos="0" zoom_1="0" zoom_2="0">
<Cursor>
<Cursor1 position="12661" topLine="293" />
</Cursor>
</File>
<File name="u1-interface.cpp" open="1" top="0" tabpos="1" split="0" active="1" splitpos="0" zoom_1="0" zoom_2="0">
<Cursor>
<Cursor1 position="3085" topLine="11" />
</Cursor>
</File>
<File name="u4-fonctions.h" open="1" top="0" tabpos="7" split="0" active="1" splitpos="0" zoom_1="0" zoom_2="0">
<Cursor>
<Cursor1 position="3233" topLine="70" />
</Cursor>
</File>
<File name="u2-dessin.cpp" open="1" top="0" tabpos="6" split="0" active="1" splitpos="0" zoom_1="0" zoom_2="0">
<Cursor>
<Cursor1 position="4023" topLine="112" />
</Cursor>
</File>
<File name="u2-dessin.h" open="1" top="0" tabpos="5" split="0" active="1" splitpos="0" zoom_1="0" zoom_2="0">
<Cursor>
<Cursor1 position="223" topLine="0" />
</Cursor>
</File>
<File name="u3-callbacks.cpp" open="1" top="0" tabpos="3" split="0" active="1" splitpos="0" zoom_1="0" zoom_2="0">
<Cursor>
<Cursor1 position="697" topLine="6" />
</Cursor>
<Folding>
<Collapse line="107" />
<Collapse line="113" />
<Collapse line="125" />
</Folding>
</File>
<File name="u1-interface.h" open="1" top="0" tabpos="2" split="0" active="1" splitpos="0" zoom_1="0" zoom_2="0">
<Cursor>
<Cursor1 position="616" topLine="0" />
</Cursor>
</File>
<File name="main.h" open="0" top="0" tabpos="6" split="0" active="1" splitpos="0" zoom_1="0" zoom_2="0">
<Cursor>
<Cursor1 position="0" topLine="0" />
</Cursor>
</File>
</CodeBlocks_layout_file>
......@@ -34,15 +34,16 @@ void TraiterCycleCB()
printf("Tour: %lf\n",10000*gDonnees.t + 1);
printf("Rho= %lf\n", Rho);
collision();
calculGeometrie();
printf("A:");
afficherTableau(gDonnees.geometrie.A,3);
printf("B:");
afficherTableau(gDonnees.geometrie.B,3);
printf("h:");
afficherTableau(gDonnees.corde.h,4);
// printf("h:");
// afficherTableau(gDonnees.corde.h,4);
Pseparation_Matlab();
collision();
// collision();
printf("k1 : %.10lf k2 : %.10lf\n", gDonnees.corde.k1, gDonnees.corde.k2);
printf("r1 : %.10lf r2 : %.10lf\n", gDonnees.corde.r1, gDonnees.corde.r2);
printf("Ug: %.12lf\n", calculUg_Matlab());
......@@ -56,6 +57,8 @@ void TraiterCycleCB()
printf("Forces:");
afficherTableau(gDonnees.corde.f, 4);
calculdeplacement();
printf("h:");
afficherTableau(gDonnees.corde.h,4);
//printf("Déplacements: y1 : %.18lf et y2 : %.18lf\n",gDonnees.corde.y1, gDonnees.corde.y2);
//printf("Déplacements: y1 t-1 : %.18lf et y2 t-1 : %.18lf\n",gDonnees.corde.ytm1[1], gDonnees.corde.ytm1[2]);
//printf("Déplacements: y1 t-2 : %.18lf et y2 t-2 : %.18lf\n",gDonnees.corde.ytm2[1], gDonnees.corde.ytm2[2]);
......@@ -184,9 +187,9 @@ void SliderCallback2(Fl_Widget* w, void* data)
k = gInterface.Slider2->value() ;
}
void SliderKCallback( Fl_Widget* w, void* data )
void SliderK1Callback( Fl_Widget* w, void* data )
{
gDonnees.corde.k1 = gInterface.SliderK->value();
gDonnees.corde.k1 = gInterface.SliderK1->value();
}
void SliderKcCallback( Fl_Widget* w, void* data )
......@@ -218,8 +221,24 @@ void SliderMCallback( Fl_Widget* w, void* data )
{
}
//j'ai modifier ça
void SliderR1Callback( Fl_Widget* w, void* data )
{
void SliderRCallback( Fl_Widget* w, void* data )
}
// et ça
void SliderR2Callback( Fl_Widget* w, void* data )
{
}
// j'ai ajouter ça
void SliderH0Callback( Fl_Widget* w, void* data )
{
}
// ça aussi
void SliderH1Callback( Fl_Widget* w, void* data )
{
}
......@@ -368,7 +368,7 @@ double calculUg_Matlab()
//if (Ugi <0), Ugi = 0;end
// dUg = (Ug-gDonnees.corde.Ugm1)*gDonnees.fe;
}
printf("xs : %.10lf et hs : %.10lf et GA : %lf\n", gDonnees.corde.xs, gDonnees.corde.hs, GA);
printf("xs : %.18lf et hs : %.18lf et GA : %lf\n", gDonnees.corde.xs, gDonnees.corde.hs, GA);
gDonnees.corde.Ugm1=Ug;
return Ug;
......@@ -414,7 +414,6 @@ double X_Matlab(int i) // terminé
//{
// X = (gDonnees.corde.x[i+1]-gDonnees.corde.x[i])/(pow(gDonnees.corde.h[i],2)); //cas où la plaque est horizontale
//}
if (fabs(gDonnees.geometrie.A[i]*gDonnees.corde.x[i+1]/gDonnees.geometrie.B[i])<pres) //cas de presque horizontalité, avec dvpt limité
{
X2 = 1/(pow(gDonnees.geometrie.B[i],2))*(pow(gDonnees.corde.x[i+1],2)/2-2*gDonnees.geometrie.A[i]*pow(gDonnees.corde.x[i+1],3)/3/gDonnees.geometrie.B[i]+3*(pow(gDonnees.geometrie.A[i],2))*pow(gDonnees.corde.x[i+1],4)/4/(pow(gDonnees.geometrie.B[i],2)));
......@@ -428,7 +427,6 @@ double X_Matlab(int i) // terminé
}
return X;
}
void X_W_s_Matlab(double* X,double* W) // entre x1 et xS, terminé
......@@ -443,8 +441,8 @@ void X_W_s_Matlab(double* X,double* W) // entre x1 et xS, terminé
// Xs = (gDonnees.corde.xs-gDonnees.corde.x[1])/(pow(gDonnees.corde.h[1],2)); //cas où la plaque est horizontale
//Ws = (gDonnees.corde.xs-gDonnees.corde.x[1])/(pow(gDonnees.corde.h[1],2)); //cas où la plaque est horizontale
//}
if (fabs(gDonnees.geometrie.A[1]*gDonnees.corde.xs/gDonnees.geometrie.B[1])<pres) //cas de presque horizontalité, avec dvpt limité
double Qs=fabs(gDonnees.geometrie.A[1]*gDonnees.corde.xs/gDonnees.geometrie.B[1]);
if ( Qs <pres) //cas de presque horizontalité, avec dvpt limité
{
X2 = 1/(pow(gDonnees.geometrie.B[1],2))*(pow(gDonnees.corde.xs,2)/2-2*gDonnees.geometrie.A[1]*pow(gDonnees.corde.xs,3)/3/gDonnees.geometrie.B[1]+3*(pow(gDonnees.geometrie.A[1],2))*pow(gDonnees.corde.xs,4)/4/(pow(gDonnees.geometrie.B[1],2)));
X1 = 1/(pow(gDonnees.geometrie.B[1],2))*(pow(gDonnees.corde.x[1],2)/2-2*gDonnees.geometrie.A[1]*pow(gDonnees.corde.x[1],3)/3/gDonnees.geometrie.B[1]+3*(pow(gDonnees.geometrie.A[1],2))*pow(gDonnees.corde.x[1],4)/4/(pow(gDonnees.geometrie.B[1],2)));
......@@ -462,7 +460,7 @@ void X_W_s_Matlab(double* X,double* W) // entre x1 et xS, terminé
//Ws = 1/gDonnees.geometrie.A[1]*(1/gDonnees.corde.h[1]-1/gDonnees.corde.hs);
*W = 1/gDonnees.geometrie.A[1]*(1/gDonnees.corde.h[1]-1/gDonnees.corde.hs);
}
printf("Qs=%.10lf\n", Qs);
}
void calculForceMatlab()
......@@ -474,6 +472,9 @@ void calculForceMatlab()
double Ug, P, bern, tmp1, tmpG, tmpH;
Ug = calculUg_Matlab();
P = calculPsub();
double Q1=fabs(gDonnees.geometrie.A[0]*gDonnees.corde.x[1]/gDonnees.geometrie.B[0]);
double Q2=fabs(gDonnees.geometrie.A[1]*gDonnees.corde.x[2]/gDonnees.geometrie.B[1]);
double Qs=fabs(gDonnees.geometrie.A[1]*gDonnees.corde.xs/gDonnees.geometrie.B[1]);
// quelques variables temporaires
tmp1 = 0.5*Rho*(Ug*Ug)/(gDonnees.corde.Lg*(gDonnees.corde.Lg)); //variable temporaire associée au terme de Bernoulli
......@@ -488,6 +489,7 @@ void calculForceMatlab()
bern = tmp1*(1/(pow(gDonnees.corde.h[0],2))-2*(X_Matlab(0)-gDonnees.corde.x[0]*W_Matlab(0))/pow((gDonnees.corde.x[1]-gDonnees.corde.x[0]),2));
printf("bern: %.10lf \n", bern);
printf("Q1= %.10lf\n", Q1);
gDonnees.corde.f[0] = 0.5*(gDonnees.corde.x[1]-gDonnees.corde.x[0])*gDonnees.corde.Lg*(P+bern);
// *******************************************************
......@@ -510,7 +512,6 @@ void calculForceMatlab()
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)));
}
gDonnees.corde.f[1] = 0.5*(gDonnees.corde.x[2]-gDonnees.corde.x[1])*gDonnees.corde.Lg*(P+bern);
}
else //cas où la glotte forme un canal divergent avec h2>hs
......@@ -535,6 +536,7 @@ void calculForceMatlab()
gDonnees.corde.f[1] = gDonnees.corde.Lg*tmpG*(P+bern)+gDonnees.corde.Lg*P_supra*(gDonnees.corde.x[2]-gDonnees.corde.xs)/(gDonnees.corde.x[2]-gDonnees.corde.x[1])*(gDonnees.corde.x[2]-(gDonnees.corde.x[2]+gDonnees.corde.xs)/2);
}
printf("Q2= %.10lf\n",Q2);
}
// *******************************************************
// ********************* Calcul de Fl_h2 *****************
......@@ -574,8 +576,8 @@ void calculForceMatlab()
bern = tmp1*(1/(pow(gDonnees.corde.h[0],2))-(Xs-gDonnees.corde.x[1]*Ws)/(tmpH*(gDonnees.corde.x[2]-gDonnees.corde.x[1])));
}
gDonnees.corde.f[2] = gDonnees.corde.Lg*tmpH*(P+bern)+ gDonnees.corde.Lg*P_supra*(gDonnees.corde.x[2]-gDonnees.corde.xs)/(gDonnees.corde.x[2]-gDonnees.corde.x[1])*((gDonnees.corde.x[2]+gDonnees.corde.xs)/2-gDonnees.corde.x[1]);
// printf("Qs=%.10lf\n", Qs );
gDonnees.corde.f[2] = gDonnees.corde.Lg*tmpH*(P+bern)+ gDonnees.corde.Lg*P_supra*(gDonnees.corde.x[2]-gDonnees.corde.xs)/(gDonnees.corde.x[2]-gDonnees.corde.x[1])*((gDonnees.corde.x[2]+gDonnees.corde.xs)/2-gDonnees.corde.x[1]);
}
//AJOUT DU PROF
if ((gDonnees.corde.h[1]>0)&&(gDonnees.corde.h[2]<=0)) //ajout pas NR (04/04/2007), pour calculer les forces de pression
......@@ -620,13 +622,16 @@ void calculdeplacement(){
printf("Déplacements: y1 t-1 : %.18lf et y2 t-1 : %.18lf\n",gDonnees.corde.ytm1[1], gDonnees.corde.ytm1[2]);
printf("Déplacements: y1 t-2 : %.18lf et y2 t-2 : %.18lf\n",gDonnees.corde.ytm2[1], gDonnees.corde.ytm2[2]);
gDonnees.corde.h[1]= 2*(gDonnees.corde.y1 + gDonnees.corde.y1_u_0); //2*gDonnees.corde.y1;
gDonnees.corde.h[2]= 2*(gDonnees.corde.y2 + gDonnees.corde.y2_u_0);
gDonnees.corde.ytm2[1]=gDonnees.corde.ytm1[1];
gDonnees.corde.ytm2[2]=gDonnees.corde.ytm1[2];
gDonnees.corde.ytm1[1]=gDonnees.corde.y1;
gDonnees.corde.ytm1[2]=gDonnees.corde.y2;
gDonnees.corde.h[1]= 2*(gDonnees.corde.y1 + gDonnees.corde.y1_u_0); //2*gDonnees.corde.y1;
gDonnees.corde.h[2]= 2*(gDonnees.corde.y2 + gDonnees.corde.y2_u_0); //2*gDonnees.corde.y2;
// gDonnees.corde.h[1]= 2*(gDonnees.corde.y1 + gDonnees.corde.y1_u_0); //2*gDonnees.corde.y1;
// gDonnees.corde.h[2]= 2*(gDonnees.corde.y2 + gDonnees.corde.y2_u_0); //2*gDonnees.corde.y2;
//gDonnees.corde.h[1]=gDonnees.geometrie.A[0]*gDonnees.corde.x[1] + gDonnees.geometrie.B[0];
......
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