[Debutante] Probleme d'affıchage du resultat d'une fonction de plusieurs variables

**SofC++** · 10/05/2007, 14h45

Bonjour,
pour mon travail je dois realiser un programme qui, a partir de certaines donnees entrees par l'utilisateur, sort le maximum d'une fonction par rapport a trois variables theta1, theta2 et X.
J'ai deja verifie que le main tout seul fonctionne, ainsi que toutes les boucles if et for, elles sont ok. Je pense que le probleme vient des trois dernieres lignes mais je ne sais pas comment faire ! J'espere avoir ete assez claire ...
Merci de votre aide !

Voilou mon petit (!) programme :

Code c :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 
 
/*Function declaration*/
float TMAXF (float,float,float,float,float,float,float,float,float,int,float,float);
 
/*Main function*/
int main()
{
    float qd,Hd,betad,id,fid,cd,rud,gammad,alphad, dholed,Shd,Tmaxd;
    int typed;
 
    /*Initialization*/
   dholed=0;
   Shd=0;
   Tmaxd=0;
 
    /*Input values*/
    printf("Enter q value (kN/m2)\n");
    scanf("%f",&qd);
    printf("Enter H value (m)\n");
    scanf("%f",&Hd);
    printf("Enter beta value (degrees)\n");
    scanf("%f",&betad);
    printf("Enter i value (degrees)\n");
    scanf("%f",&id);
    printf("Enter fi value (degrees)\n");
    scanf("%f",&fid);
    printf("Enter c value (kN/m2)\n");
    scanf("%f",&cd);
    printf("Enter ru value\n");
    scanf("%f",&rud);
    printf("Enter gamma value (kN/m3)\n");
    scanf("%f",&gammad); 
    printf("Enter interface sliding value (Entry 999 for a default value)\n");
    scanf("%f",&alphad);
    printf("Type of reinforcement : 1->Geotextiles/Geogrids, 2->Strip reinforcement, 3->Soil nails\n");
    scanf("%d",&typed);
    if (typed==3)
    {
                printf("Enter the diameter of the grout hole around the nail\n");
                scanf("%f",&dholed);
                printf("Enter the horizontal spacing\n");
                scanf("%f",&Shd);
    }
 
     Tmaxd=TMAXF(qd,Hd,betad,id,fid,cd,rud,gammad,alphad,typed,dholed,Shd);
     printf("Tmax=%f",Tmaxd);
 
    /*Wait before quit*/
    system("pause");
 
}

Code c :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
 
/*TMAXF function*/
 
float TMAXF(float q,float H,float beta,float i,float fi,float c,float ru, float gamma, float alpha,int type,float dhole, float Sh)
{
 
    float Tmax, H1,L,Xlim,theta1,theta2,X,Tmax1,Y,a,b,k,d,e,f,g,s,t,j,m,u,v,w,z,lambdas,W1,W2,U1,U2,K1,K2;
    X=0;
    Tmax=0;
    Tmax1=0;
 
    /*alpha default values*/
    if (alpha==999)
    {
       switch (type)
       {
       case 1:
       alpha=0.85;
       break;
       case 2:
       alpha=0.9;
       break;
       case 3:
       alpha=0.9;
       break;
       }
    }
 
    /*H1=Height with surcharge*/
    H1=H+(q/gamma);
    /*L=length from the toe to the crest*/
    L=H1/tan(beta*M_PI/180);
   /*Imposed limit value for X*/  
    Xlim=2.5*H1;   
 
    for (theta1=0;theta1<=90;theta1=theta1+0.1)
    {
        for (theta2=0;theta2<=90;theta2=theta2+0.1)
        {
            for (X=0;X<Xlim;X=X+0.1)
            {
                Y=X*tan(theta1*M_PI/180);
                a=(H1-X*tan(beta*M_PI/180));
                b=(H1-Y-a);
                k=((a+b)/tan(theta1*M_PI/180))-a/tan(beta*M_PI/180);
                d=k*tan(theta1*M_PI/180);
                e=k/cos(theta1*M_PI/180);
                f=((a+b)/sin(theta1*M_PI/180))-e;
                g=(a+b)/sin(theta1*M_PI/180);
                s=(a+b)/tan(theta1*M_PI/180);
                t=k-s;
                j=t*tan(theta2*M_PI/180);
                m=sqrt(X+Y);
                u=b/(tan(theta1*M_PI/180)-tan(beta*M_PI/180));
                   if (i==0)
                   v=0;                
                   else 
                   v=k/((1/tan(i*M_PI/180))-(1/tan(theta1*M_PI/180)));
                w=v/sin(theta1*M_PI/180);
                z=t*tan(i*M_PI/180);
 
                /*Different expressions of lambdas depending on the 
                reinforcement material*/
                switch (type)
                {
                case 1:
                lambdas=alpha;
                break;
                case 2:
                lambdas=alpha*b+(1-b);
                break;
                case 3:
                lambdas=(alpha*dhole/Sh)+(1-dhole/Sh);
                break;
                }
 
                if (X<=(H1/tan(beta*M_PI/180)))
                {
                   if (X+(H1-L*tan(theta1*M_PI/180))/(tan(theta2*M_PI/180)-tan(theta1*M_PI/180))<L)
                   /*CASE 3*/
                   {
                   W1=0.5*gamma*b*u;
                   W2=0.5*gamma*b*X;
                   U1=0.5*ru*gamma*b*u/cos(theta1*M_PI/180);
                   U2=0.5*ru*gamma*b*m;
                   K1=c*u/cos(theta1*M_PI/180);
                   K2=c*m;
                   }
 
 
                   else 
                   /*CASE 1*/
                   {
                   W1=0.5*gamma*((a+b)*(a+b)/tan(theta1*M_PI/180)-a*a/tan(beta*M_PI/180)+v*k);
                   W2=0.5*gamma*b*X;
                   U1=0.5*ru*gamma*(d*(e+w)+f*(d+b));
                   U2=0.5*ru*gamma*b*m;
                   K1=c*(g+w);
                   K2=c*m;
                   }
 
                }
 
 
                else
                /*CASE 2*/
                {
                W1=0.5*gamma*(H1-Y)*(H1-Y)/tan(theta1*M_PI/180)+0.5*gamma*(v*k-t*z);
                W2=0.5*gamma*(2*X*H1-X*tan(theta2*M_PI/180)-H1*H1/tan(beta*M_PI/180)+t*t*z);
                U1=0.5*ru*gamma*(g+w)*(H1-Y+z);
                U2=0.5*ru*((H1-Y)*(X+t)/cos(theta2*M_PI/180)+j*m+z*t/cos(theta2*M_PI/180));
                K1=c*(g+w);
                K2=c*m;
                }   
 
                  Tmax1=((W1*(tan(theta1*M_PI/180)-tan(fi*M_PI/180))+(U1*tan(fi*M_PI/180)-K1)/cos(theta1*M_PI/180))/(1+tan(theta1*M_PI/180)*tan(fi*M_PI/180)))+((W2*(tan(theta2*M_PI/180)-lambdas*tan(fi*M_PI/180))+lambdas*(U2*tan(fi*M_PI/180)-K2)/cos(theta2*M_PI/180))/(1+lambdas*tan(theta2*M_PI/180)*tan(fi*M_PI/180)));
 
 
                   /*Take the maximum value of Tmax*/
                   if (Tmax1>Tmax) Tmax=Tmax1;
 
                     }
 
             }
        }
        return (Tmax);
    }

**ZeRiL** · 10/05/2007, 14h55

Euh c'est quoi le probleme exactement ? j'avoue ne pas avoir tout compris la

la valeur retournée n'est pas bonne ?

PS : bouuuuh printf !!! bouuuuuh scanf !!!

**juls64** · 10/05/2007, 14h57

Au fait pour que le code soit plus lisible utilise le # et tu mets ton code au lieu il y aura la coloration syntaxique sinon ca donne pas envie

**Médinoc** · 10/05/2007, 15h45

C'est plus du C que du C++, ce code...

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

float Tmax, H1,L,Xlim,theta1,theta2,X,Tmax1,Y,a,b,k,d,e,f,g,s,t,j,m,u,v,w,z,lambdas,W1,W2,U1,U2,K1,K2;

Tu es sûr que tu ne peux pas faire plus court, comme noms de variables ?

Ne compte pas sur nous pour t'aider si tu produis un code incompréhensible (et pour ainsi dire non commenté)...

**HanLee** · 10/05/2007, 21h01

Ben disons que c'est plutôt le schéma qui nous manque... Les noms de variables désignent clairement des trucs de physique...

**SofC++** · 11/05/2007, 09h47

Désolée pour les couleurs je suis nouvelle sur le site je n'avais pas vu cette possibilité.

Effectivement mes variables ont des noms longs mais ce programme n'est que le début d'une longue suite, et il servira a dimensionner des murs de soutènement. C'est donc un problème physique, et je garde des noms proches des paramètres physiques pour m'y retrouver...
Quant a mon problème, je m'explique :
déjà je pense que le pas que j'utilisais pour mes variables theta1, theta2 et X était trop petit, si bien que les calculs duraient des heures et je n'attendais jamais assez longtemps... J'ai diminué le pas.
Ensuite la j'ai introduit une valeur p pour tester le programme pour une valeur particulière je veux qu'il me sorte les valeurs de W1, W2, U1, U2, K1 et K2. Les valeurs de W1 et W2 sont surprenantes ....En plus, au lieu de me sortir p=1,2 ou 3 il me sort une valeur aberrante !
Enfin, et c'est la la problème principal, la valeur de Tmax en sortie est complètement aberrante ...
J'ai essaye de rajouter des commentaires, j'espère que j'ai été suffisamment claire cette fois-ci merci de votre aide !

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 
 
/*Function declaration*/
float TMAXF (float,float,float,float,float,float,float,float,float,int,float,float);
 
/*Main function*/
int main()
{
    float qd,Hd,betad,id,fid,cd,rud,gammad,alphad, dholed,Shd,Tmaxd;
    int typed;
 
    /*Initialization*/
   dholed=0;
   Shd=0;
   Tmaxd=0;
 
    /*Input values*/
    printf("Enter q value (kN/m2)\n");
    scanf("%f",&qd);
    printf("Enter H value (m)\n");
    scanf("%f",&Hd);
    printf("Enter beta value (degrees)\n");
    scanf("%f",&betad);
    printf("Enter i value (degrees)\n");
    scanf("%f",&id);
    printf("Enter fi value (degrees)\n");
    scanf("%f",&fid);
    printf("Enter c value (kN/m2)\n");
    scanf("%f",&cd);
    printf("Enter ru value\n");
    scanf("%f",&rud);
    printf("Enter gamma value (kN/m3)\n");
    scanf("%f",&gammad); 
    printf("Enter interface sliding value (Entry 000 for a default value)\n");
    scanf("%f",&alphad);
    printf("Type of reinforcement : 1->Geotextiles/Geogrids, 2->Strip reinforcement, 3->Soil nails\n");
    scanf("%d",&typed);
    if (typed==3)
    {
                printf("Enter the diameter of the grout hole around the nail\n");
                scanf("%f",&dholed);
                printf("Enter the horizontal spacing\n");
                scanf("%f",&Shd);
    }
 
     /*Return Tmax value*/
     Tmaxd=TMAXF(qd,Hd,betad,id,fid,cd,rud,gammad,alphad,typed,dholed,Shd);
     printf("Tmax=%f",Tmaxd);
 
    /*Wait before quit*/
    system("pause");
 
} 
 
 
 
 
 
 
 
 
 
 
 
/*TMAXF function*/
 
float TMAXF(float q,float H,float beta,float i,float fi,float c,float ru, float gamma, float alpha,int type,float dhole, float Sh)
{
 
    float Tmax, H1,L,Xlim,theta1,theta2,X,Tmax1,Y,a,b,k,d,e,f,g,s,t,j,m,u,v,w,z,lambdas,W1,W2,U1,U2,K1,K2;
    int p;
    X=0;
    Tmax=0;
    Tmax1=0;
    p=0;
 
    /*alpha default values*/
    if (alpha==000)
    {
       switch (type)
       {
       case 1:
       alpha=0.85;
       break;
       case 2:
       alpha=0.9;
       break;
       case 3:
       alpha=0.9;
       break;
       }
    }
 
    /*H1=Height with surcharge*/
    H1=H+(q/gamma);
    /*L=length from the toe to the crest*/
    L=H1/tan(beta*M_PI/180);
   /*Imposed limit value for X*/  
    Xlim=2.5*H1;   
 
    /*Scanning of all possible combinations for theta1, theta2 and X*/
    for (theta1=1;theta1<=90;theta1=theta1+1)
    {
        for (theta2=1;theta2<=90;theta2=theta2+1)
        {
            for (X=0;X<Xlim;X=X+0.2)
            {
                /*Definition of parameters for all cases*/
                Y=X*tan(theta1*M_PI/180);
                a=(H1-X*tan(beta*M_PI/180));
                b=(H1-Y-a);
                k=((a+b)/tan(theta1*M_PI/180))-a/tan(beta*M_PI/180);
                d=k*tan(theta1*M_PI/180);
                e=k/cos(theta1*M_PI/180);
                f=((a+b)/sin(theta1*M_PI/180))-e;
                g=(a+b)/sin(theta1*M_PI/180);
                s=(a+b)/tan(theta1*M_PI/180);
                t=k-s;
                j=t*tan(theta2*M_PI/180);
                m=sqrt(X+Y);
                u=b/(tan(theta1*M_PI/180)-tan(beta*M_PI/180));
                   if (i==0)
                   v=0;                
                   else 
                   v=k/((1/tan(i*M_PI/180))-(1/tan(theta1*M_PI/180)));
                w=v/sin(theta1*M_PI/180);
                z=t*tan(i*M_PI/180);
 
                /*Different expressions of lambdas depending on the 
                reinforcement material*/
                switch (type)
                {
                case 1:
                lambdas=alpha;
                break;
                case 2:
                lambdas=alpha*b+(1-b);
                break;
                case 3:
                lambdas=(alpha*dhole/Sh)+(1-dhole/Sh);
                break;
                }
 
 
                /*Different calculations of forces depending on the case*/
                if (X<=(H1/tan(beta*M_PI/180)))
                {
                   if ((X<((H1-L*tan(theta1*M_PI/180))/(tan(theta2*M_PI/180)-tan(theta1*M_PI/180))))&&(X+(H-X*tan(theta2*M_PI/180)/tan(theta1*M_PI/180))))
                   /*CASE 3*/
                   {
                   /*Weight of wedge 1*/
                   W1=0.5*gamma*b*u;
                   /*Weight of wedge 2*/
                   W2=0.5*gamma*b*X;
                   /*Pore water force acting on base of wedge 1*/
                   U1=0.5*ru*gamma*b*u/cos(theta1*M_PI/180);
                   /*Pore water force acting on base of wedge 2*/
                   U2=0.5*ru*gamma*b*m;
                   /*Cohesion force acting on base of wedge 1*/
                   K1=c*u/cos(theta1*M_PI/180);
                   /*Cohesion force acting on base of wedge 2*/
                   K2=c*m;
                   p=3;
                   }
 
 
                   else 
                   /*CASE 1*/
                   {
                   W1=0.5*gamma*((a+b)*(a+b)/tan(theta1*M_PI/180)-a*a/tan(beta*M_PI/180)+v*k);
                   W2=0.5*gamma*b*X;
                   U1=0.5*ru*gamma*(d*(e+w)+f*(d+b));
                   U2=0.5*ru*gamma*b*m;
                   K1=c*(g+w);
                   K2=c*m;
                   p=1;
                   }
 
                }
 
 
                else
                /*CASE 2*/
                {
                W1=0.5*gamma*(H1-Y)*(H1-Y)/tan(theta1*M_PI/180)+0.5*gamma*(v*k-t*z);
                W2=0.5*gamma*(2*X*H1-X*tan(theta2*M_PI/180)-H1*H1/tan(beta*M_PI/180)+t*t*z);
                U1=0.5*ru*gamma*(g+w)*(H1-Y+z);
                U2=0.5*ru*((H1-Y)*(X+t)/cos(theta2*M_PI/180)+j*m+z*t/cos(theta2*M_PI/180));
                K1=c*(g+w);
                K2=c*m;
                p=2;
                }   
 
                  if (theta1==50&&theta2==10&&X==1) 
                  {
                  printf("p=%f\n",p);
                  printf("W1=%f\n",W1);
                  printf("W2=%f\n",W2);
                  printf("U1=%f\n",U1);
                  printf("U2=%f\n",U2);
                  printf("K1=%f\n",K1);
                  printf("K2=%f\n",K2);
                  }
 
                  /*Total quantity of horizontal reinforcement force required*/
                  Tmax1=((W1*(tan(theta1*M_PI/180)-tan(fi*M_PI/180))+(U1*tan(fi*M_PI/180)-K1)/cos(theta1*M_PI/180))/(1+tan(theta1*M_PI/180)*tan(fi*M_PI/180)))+((W2*(tan(theta2*M_PI/180)-lambdas*tan(fi*M_PI/180))+lambdas*(U2*tan(fi*M_PI/180)-K2)/cos(theta2*M_PI/180))/(1+lambdas*tan(theta2*M_PI/180)*tan(fi*M_PI/180)));
 
                   /*Take the maximum value of Tmax*/
                   if (Tmax1>Tmax) Tmax=Tmax1;
 
                     }
 
             }
        }
        return (Tmax);
    }
 
////////////////////////////////////////////////////////////////////////////////

**SofC++** · 11/05/2007, 10h51

Alors j'ai change mes bornes pour la variation de mes angles parce que la variations de 0 a 90 pour les deux n'était pas physique.

J'ai donc des valeurs de Tmax beaucoup moins aberrantes .... mais toujours fausses !

J'ai bien vérifié mes expressions des différents termes, elles sont juste, j'ai même ajoute des parenthèses pour être sur des priorités sur les opérations.
Je demande toujours p et il me sort toujours une valeur bizarre ...

Merci de votre aide !

[Debutante] Probleme d'affıchage du resultat d'une fonction de plusieurs variables

C++

Vue hybride

Discussions similaires

Partager

Partager