1 : * TRACTUFI PROCEDUR NOBODY 92/01/01 00:00:00 001 2 : 'DEBPROC' TRACTUFI ETAB*'TABLE'; 3 : ************************************************************************ 4 : * * 5 : * PROCEDURE TRANSFORMANT LA COURBE DE TRACTION * 6 : * * 7 : * (SIG,EPS) DU MATERIAU EN COURBE DE TRACTION * 8 : * * 9 : * EXPRIMEE EN TERMES DE VARIABLES GLOBALES * 10 : * * 11 : * MOMENT/ROTATION * 12 : * * 13 : ************************************************************************ 14 : ************************************************************************ 15 : * * 16 : * LA PROCEDURE 'TRACTUFI' UTILISE : * 17 : * * 18 : * SOIT UNE BASE DE DONNEES EXPERIMENTALES * 19 : * - DEFR 20 : * * 21 : * SOIT DES METHODES SIMPLIFIEES * 22 : * - TADA-PARIS * 23 : * - LBB-NRC * 24 : * - LBB1 * 25 : * - LBB2 * 26 : * * 27 : * ELLE FAIT APPEL AUX ARGUMENTS SUIVANTS: * 28 : * * 29 : * ENTREES * 30 : * * 31 : * ETAB.'METHODE' = MOT CLE INDIQUANT LA METHODE EMPLOYEE * 32 : * ETAB.'COUTRA' = COURBE DE TRACTION (SIG,EPS) DU MATERIAU * 33 : * ETAB.'YOUN' = MODULE D'YOUNG * 34 : * ETAB.'SIG1' = CONTRAINTE POUR UNE DEFORMATION DE 0.002 * 35 : * OU ETAB.'SIGF' = CONTRAINTE D'ECOULEMENT * 36 : * ETAB.'REXT' = RAYON EXTERIEUR DU TUYAU FISSURE * 37 : * ETAB.'EPAI' = EPAISSEUR DU TUYAU FISSURE * 38 : * ETAB.'ANGLE' = ANGLE TOTAL DE LA FISSURE (EN DEGRES) * 39 : * ETAB.'COUL' = COULEUR AFFECTEE A LA COURBE * 40 : * ETAB.'ALFA' = VALEURS FACULTATIVES PERMETTANT DE FITER * 41 : * ETAB.'N' = LA COURBE DE TRACTION (METHODE LBBNRC) * 42 : * * 43 : * * 44 : * SORTIES * 45 : * * 46 : * MP = COURBE (M,PHI) * 47 : * CM = COEFFICIENT DE COMPLAISANCE * 48 : * KF = POUR LE CALCUL DU FACTEUR D"INTENSITE DES * 49 : * CONTRAINTES * 50 : * * 51 : ************************************************************************ 52 : *----------------------------------------------------------------------- 53 : * RECEPTION DES ARGUMENTS 54 : 55 : METHODE=ETAB.'METHODE'; 56 : COUTRA=ETAB.'COUTRA';E=ETAB.'YOUN'; 57 : REXT=ETAB.'REXT';EPAI=ETAB.'EPAI';AN=ETAB.'ANGLE'; 58 : 59 : 'SI' ('OU' (AN < 0.) (AN > 360.)); 60 : MESS 'ANGLE DE FISSURE IRREALISTE'; 61 : 'QUITTER' TRACTUFI; 62 : 'FINSI'; 63 : 64 : PI=3.14159D0; 65 : RAYOM=REXT - (EPAI / 2.D0); 66 : THETA=AN / 2.D0; 67 : XMSS0=4.D0 * RAYOM * RAYOM * EPAI * 68 : (('COS' (THETA / 2.D0)) - (0.5D0 * ('SIN' THETA))); 69 : THETAD=THETA; 70 : THETA=THETA * PI / 180.D0; 71 : 72 : LEPS='EXTR' COUTRA ABSC; 73 : LSIG='EXTR' COUTRA ORDO; 74 : NCOURB='DIME' LSIG; 75 : PM='PROG' NCOURB * 0.; 76 : PPHI='PROG' NCOURB * 0.; 77 : SIGM='EXTR' LSIG NCOURB; 78 : 79 : 'SI' ('EXISTE' ETAB 'SIGF'); 80 : SIGF= ETAB . 'SIGF'; 81 : 'SINON'; 82 : SIG1 = ETAB . 'SIG1'; 83 : SIGF= (SIG1 + SIGM) / 2.; 84 : 'FINSI'; 85 : SIGY='EXTR' LSIG 2; 86 : ************************************************************************ 87 : * * 88 : * DETERMINATION DE LA COMPLAISANCE CM DE L'ELEMENT ET DU COEFFICIENT * 89 : * KF * 90 : ************************************************************************ 91 : RST=RAYOM / EPAI; 92 : 'SI' (RST < 10.D0) ; 93 : AXX=(0.125D0 * RST - 0.25D0) ** 0.25D0; 94 : 'SINON'; 95 : AXX=(0.4D0 * RST - 3.D0) ** 0.25D0; 96 : 'FINSI'; 97 : TSP=THETA / PI; 98 : TSP1=TSP ** 1.5D0; 99 : TSP2=TSP ** 4.24D0; 100 : TSP3=TSP ** 3.D0; 101 : TSP4=TSP ** 5.74D0; 102 : TSP5=TSP ** 8.48D0; 103 : FM=AXX * ((5.26D0 * TSP1) + (1.69D0 * TSP2)); 104 : FM=FM + 1.D0; 105 : FM=FM + 106 : (AXX * AXX * ((8.46D0 * TSP3) + (6.28D0 * TSP4) + (1.33D0*TSP5))); 107 : CM=2.D0 * THETA * THETA * FM; 108 : CM=CM / E / PI / RAYOM / RAYOM / EPAI; 109 : 110 : 'SI' (THETAD < 100.); 111 : FF = 1. + (AXX * ((4.6 * (TSP ** 1.5)) + 112 : (2.64 * (TSP ** 4.24)))); 113 : 114 : 'SINON'; 115 : 'MESS' 'ANGLE DE FISSURE TROP GRAND'; 116 : 'FINSI'; 117 : K1 = PI * RAYOM * THETA; 118 : K2 = K1 ** 0.5; 119 : K3 = K2 * FF; 120 : KF = K3 / PI / RAYOM / RAYOM / EPAI; 121 : 122 : 123 : 124 : ************************************************************************ 125 : * * 126 : * DETERMINATION DES COEFFICIENTS DE LA LOI DE RAMBERG - OSGOOD * 127 : * POUR LA METHODE LBB-NRC * 128 : * * 129 : ************************************************************************ 130 : 'SI' (('EXISTE' ETAB 'ALFA') ET ('EXISTE' ETAB 'N')); 131 : ALFA = ETAB.'ALFA'; 132 : N = ETAB.'N'; 133 : 'SINON'; 134 : 'SI' ('EGA' METHODE 'LBBNRC'); 135 : EPSF='IPOL' SIGF LSIG LEPS; 136 : EPSY='EXTR' LEPS 2;SIGY='EXTR' LSIG 2; 137 : ALFA=EPSF / EPSY - 1.D0; 138 : I=1; 139 : 'REPETER' BOUCR (NCOURB - 1); 140 : I=I + 1; 141 : SIG2='EXTR' LSIG I; 142 : 'SI' (SIGF '<EG' SIG2); 143 : SIG1='EXTR' LSIG (I - 1); 144 : EPS2=('EXTR' LEPS I) - (('EXTR' LEPS I) / E); 145 : EPS1=('EXTR' LEPS (I - 1)) - (('EXTR' LEPS (I - 1)) / E); 146 : PENTE=(EPS2 - EPS1) / (SIG2 - SIG1); 147 : 'QUITTER' BOUCR; 148 : 'FINSI'; 149 : 'FIN' BOUCR; 150 : N=(SIGF * PENTE / ALFA / EPSY) + ((SIGF / SIGY - 1.D0) / ALFA); 151 : 'SI' (('VALE' IMPI) 'EGA' 1); 152 : I=1; 153 : 'REPETER' BOUCLE (NCOURB - 1); 154 : I=I + 1; 155 : VSIG='EXTR' LSIG I; 156 : RSIG=VSIG / SIGF; 157 : REPS=RSIG + (ALFA * (RSIG ** N)); 158 : VEPS=EPSY * REPS; 159 : 'REMP' PM I VSIG;'REMP' PPHI I VEPS; 160 : 'FIN' BOUCLE; 161 : COURAM='EVOL' ROUG MANU 'EPS' PPHI 'SIG' PM; 162 : 'DESS' (COUTRA ET COURAM) MARQ CARR ETOI; 163 : 'FINSI'; 164 : 'FINSI'; 165 : 'FINSI'; 166 : ************************************************************************ 167 : * * 168 : * DETERMINATION DU POINT CORRESPONDANT A LA LIMITE ELASTIQUE * 169 : * * 170 : ************************************************************************ 171 : SIGY='EXTR' LSIG 2; 172 : MY=XMSS0 * SIGY; 173 : PHIY=CM * MY; 174 : 'REMP' PM 2 MY;'REMP' PPHI 2 PHIY; 175 : ************************************************************************ 176 : * * 177 : * POUR LA METHODE TADA-PARIS DETERMINATION DE THELIM ET BETA * 178 : * * 179 : ************************************************************************ 180 : 'SI' ('OU' ('OU' ('EGA' METHODE 'TADA') ('EGA' METHODE 'LBBNRC')) 181 : ('EGA' METHODE 'LBB1')); 182 : SP=4. * (('COS' (THETAD / 2.D0)) - (0.5D0 * ('SIN' THETAD))) / PI; 183 : THETA0=THETA; 184 : K=0; 185 : 'REPETER' BOUC1 200; 186 : K=K + 1; 187 : TSP=THETA / PI; 188 : FZM=1.D0 + (AXX * ((4.6D0 * (TSP ** 1.5D0)) 189 : + (2.64D0 * (TSP ** 4.24D0)))); 190 : DFZM=AXX * ((6.9D0 * (TSP ** 0.5D0)) 191 : + (2.64D0 * 4.24D0 * (TSP ** 3.24D0))); 192 : THELIM=THETA0 + ((THETA * FZM) / (FZM + (2.D0 * TSP * DFZM))); 193 : 'SI' (((THELIM - THETA) / THETA) > 1.D-5); 194 : THETA=THELIM; 195 : 'SINON'; 196 : 'QUITTER' BOUC1; 197 : 'FINSI'; 198 : 'FIN' BOUC1; 199 : 'SI' (K 'EGA' 200); 200 : 'MESS' 'PAS DE CONVERGENCE POUR LA DETERMINATION DE THELIM'; 201 : 'QUITTER' TRACTUFI; 202 : 'FINSI'; 203 : THETAY=THELIM - THETA0; 204 : TSP=THELIM / PI; 205 : FZM=1.D0 + (AXX * ((4.6D0 * (TSP ** 1.5D0)) 206 : + (2.64D0 * (TSP ** 4.24D0)))); 207 : GZM=THELIM * FZM * FZM; 208 : DGZM=GZM / THETAY; 209 : BETA=(SP ** 2) * DGZM; 210 : 'FINSI'; 211 : ************************************************************************ 212 : * * 213 : * POUR LA METHODE DEFR DETERMINATION DE LA COURBE PHIPLSEL-EPSPLSEL * 214 : * * 215 : ************************************************************************ 216 : 'SI' ('EGA' METHODE 'DEFR'); 217 : * ---------------------------------------------------------------- 218 : 219 : * COURBES (EPSPL/EL,PHIPL/EL) 220 : 221 : LP15=PROG 0. 0.67 2.67 4.54 5.31 6.06 6.87 8.05 9.14 10.69 11.43 13.37 222 : 15.07 16.93 18.12 20.77 22.14 23.11 25.08 27.78 29.35 30.53 35.14 35.91 223 : 36.51 37.28 39.21 39.57 40.34; 224 : LP15=LP15 * 100.; 225 : 226 : LE15=PROG 0. 0.69 2.02 3.14 3.48 3.79 4.13 4.45 4.75 5.19 5.35 5.69 6.03 227 : 6.30 6.49 6.85 7.03 7.16 7.47 7.75 7.92 8.07 8.39 8.48 8.58 8.76 9.22 228 : 9.43 9.79; 229 : LE15=LE15 * 10.; 230 : 231 : N=DIME LP15; 232 : X1=EXTR LP15 (N-1); 233 : X2=EXTR LP15 N; 234 : Y1=EXTR LE15 (N-1); 235 : Y2=EXTR LE15 N; 236 : P=(Y2 - Y1) / (X2 - X1); 237 : X3=((201. - Y2) / P) + X2; 238 : LP15=LP15 ET (PROG X3); 239 : LE15=LE15 ET (PROG 201.); 240 : 241 : INT15=EVOL MANU 'PHIPLSEL' LP15 'EPSPLSEL' LE15; 242 : * 243 : LP30=PROG 0. 0.27 0.53 0.74 0.96 1.17 1.40 1.73 2.08 2.43 2.78 3.12 3.46 244 : 3.91 4.48 5.07 5.64 6.32 6.99 7.57 8.26 8.82 9.27 9.62 10.78 11.56 245 : 12.12 12.95 13.46 13.80 14.43 14.66 14.89 15.12; 246 : LP30=LP30 * 100.; 247 : 248 : LE30=PROG 0. 1.21 2.24 3.02 3.64 4.25 4.67 5.24 5.64 6.03 6.38 6.67 6.94 249 : 7.36 7.88 8.33 8.73 9.29 9.82 10.08 10.58 11.04 11.45 11.71 12.73 13.34 250 : 13.75 14.61 15.22 15.62 16.59 17.24 17.89 18.54; 251 : LE30=LE30 * 10.; 252 : 253 : N=DIME LP30; 254 : X1=EXTR LP30 (N-1); 255 : X2=EXTR LP30 N; 256 : Y1=EXTR LE30 (N-1); 257 : Y2=EXTR LE30 N; 258 : P=(Y2 - Y1) / (X2 - X1); 259 : X3=((201. - Y2) / P) + X2; 260 : LP30=LP30 ET (PROG X3); 261 : LE30=LE30 ET (PROG 201.); 262 : 263 : INT30=EVOL ROUG MANU LP30 LE30; 264 : * 265 : LP45=PROG 0. 0.15 0.27 0.35 0.43 0.52 0.60 0.69 0.78 0.91 1.04 1.17 1.30 266 : 1.43 1.56 1.69 1.78 2.00 2.13 2.43 2.60 2.77 3.00 3.26 3.43 3.56 3.79 267 : 3.92 4.04 4.14 4.19 4.30 4.36 4.41 4.49 4.54; 268 : LP45=LP45 * 100.; 269 : 270 : LE45=PROG 0. 0.88 1.78 2.41 2.93 3.39 3.84 4.17 4.50 4.90 5.24 5.52 5.83 271 : 6.18 6.43 6.67 6.86 7.26 7.53 7.91 8.20 8.50 8.87 9.19 9.36 9.63 9.88 272 : 9.98 10.03 10.19 10.28 10.47 10.68 10.83 11.20 11.45; 273 : LE45=LE45 * 10.; 274 : 275 : N=DIME LP45; 276 : X1=EXTR LP45 (N-1); 277 : X2=EXTR LP45 N; 278 : Y1=EXTR LE45 (N-1); 279 : Y2=EXTR LE45 N; 280 : P=(Y2 - Y1) / (X2 - X1); 281 : X3=((201. - Y2) / P) + X2; 282 : LP45=LP45 ET (PROG X3); 283 : LE45=LE45 ET (PROG 201.); 284 : 285 : INT45=EVOL ROUG MANU LP45 LE45; 286 : * 287 : LP60=PROG 0. 0.09 0.17 0.56 0.89 1.24 1.60 2.01 2.47 2.86 3.28 3.72 4.14 288 : 4.68 5.07 5.46 5.88 6.34 6.75 7.19 7.63 8.02 8.65 9.04 9.40 10.20 10.79 289 : 11.26 11.96 12.41 12.86 13.26; 290 : LP60=LP60 * 10.; 291 : 292 : LE60=PROG 0. 1.79 2.36 3.55 4.23 4.66 4.90 5.16 5.44 5.65 5.87 6.19 6.33 293 : 6.62 6.84 6.99 7.23 7.54 7.80 8.06 8.26 8.44 8.74 8.94 9.23 10.15 10.79 294 : 11.07 11.55 11.91 12.20 12.49; 295 : LE60=LE60 * 10.; 296 : 297 : N=DIME LP60; 298 : X1=EXTR LP60 (N-1); 299 : X2=EXTR LP60 N; 300 : Y1=EXTR LE60 (N-1); 301 : Y2=EXTR LE60 N; 302 : P=(Y2 - Y1) / (X2 - X1); 303 : X3=((201. - Y2) / P) + X2; 304 : LP60=LP60 ET (PROG X3); 305 : LE60=LE60 ET (PROG 201.); 306 : 307 : INT60=EVOL ROUG MANU LP60 LE60; 308 : * 309 : LP75=PROG 0. 0.2 0.49 0.72 1.05 1.36 1.63 2.03 2.30 2.70 3.09 3.48 3.74 310 : 4.12 4.50 4.87 5.25 5.50 5.96 6.28 6.53 6.79 6.98 7.18; 311 : LP75=LP75 * 10.; 312 : 313 : LE75=PROG 0. 0.70 1.12 1.95 2.53 3.14 3.55 4.06 4.34 4.76 5.16 5.51 5.74 314 : 6.16 6.54 6.99 7.38 7.67 8.14 8.51 8.99 9.53 10.06 11.34; 315 : LE75=LE75 * 10.; 316 : 317 : N=DIME LP75; 318 : X1=EXTR LP75 (N-1); 319 : X2=EXTR LP75 N; 320 : Y1=EXTR LE75 (N-1); 321 : Y2=EXTR LE75 N; 322 : P=(Y2 - Y1) / (X2 - X1); 323 : X3=((201. - Y2) / P) + X2; 324 : LP75=LP75 ET (PROG X3); 325 : LE75=LE75 ET (PROG 201.); 326 : 327 : INT75=EVOL ROUG MANU LP75 LE75; 328 : * 329 : LP90=PROG 0. 0.11 0.18 0.28 0.37 0.54 0.72 0.82 0.94 1.14 1.27 1.47 1.58 330 : 1.68 1.77 1.90 1.99 2.10 2.28 2.37 2.45 2.57 2.79 3.04; 331 : LP90=LP90 * 10.; 332 : 333 : LE90=PROG 0. 1.44 2.29 3.09 3.69 4.35 4.95 5.15 5.38 5.75 5.92 6.23 6.41 334 : 6.65 6.81 7.02 7.17 7.40 7.84 8.06 8.31 8.72 9.44 10.27; 335 : LE90=LE90 * 10.; 336 : 337 : N=DIME LP90; 338 : X1=EXTR LP90 (N-1); 339 : X2=EXTR LP90 N; 340 : Y1=EXTR LE90 (N-1); 341 : Y2=EXTR LE90 N; 342 : P=(Y2 - Y1) / (X2 - X1); 343 : X3=((201. - Y2) / P) + X2; 344 : LP90=LP90 ET (PROG X3); 345 : LE90=LE90 ET (PROG 201.); 346 : 347 : INT90=EVOL ROUG MANU LP90 LE90; 348 : * 349 : LP120=PROG 0. 0.38 0.80 1.34 1.62 2.47 3.65 4.23 4.66 5.86 6.37 7.09 350 : 7.83 8.48 9.11 10.41 11.83 12.32 12.90 13.47 14.70 15.63 16.58; 351 : 352 : LE120=PROG 0. 0.43 1.14 1.82 2.08 3.04 3.96 4.23 4.42 4.86 5.07 5.34 353 : 5.62 5.83 6.04 6.68 7.18 7.44 7.67 8.00 8.60 9.09 10.19; 354 : LE120=LE120 * 10.; 355 : 356 : N=DIME LP120; 357 : X1=EXTR LP120 (N-1); 358 : X2=EXTR LP120 N; 359 : Y1=EXTR LE120 (N-1); 360 : Y2=EXTR LE120 N; 361 : P=(Y2 - Y1) / (X2 - X1); 362 : X3=((201. - Y2) / P) + X2; 363 : LP120=LP120 ET (PROG X3); 364 : LE120=LE120 ET (PROG 201.); 365 : 366 : INT120=EVOL ROUG MANU LP120 LE120; 367 : * 368 : LP150=PROG 0. 0.64 0.85 1.3 1.76 2.16 2.55 2.8 3.13 4.1 4.56 5.33 6.01 369 : 6.43 7.06 7.63 8.08 8.48 8.93 9.78 10.56 10.79 11.24 11.69 12.87 13.75 370 : 14.38; 371 : 372 : LE150=PROG 0. 2.01 3. 3.63 4.23 4.69 5.1 5.38 5.67 6.12 6.34 6.78 7.2 373 : 7.44 7.82 8.17 8.44 8.73 8.98 9.52 10.09 10.32 10.63 10.94 12.13 13.35 374 : 15.20; 375 : LE150=LE150 * 10.; 376 : 377 : N=DIME LP150; 378 : X1=EXTR LP150 (N-1); 379 : X2=EXTR LP150 N; 380 : Y1=EXTR LE150 (N-1); 381 : Y2=EXTR LE150 N; 382 : P=(Y2 - Y1) / (X2 - X1); 383 : X3=((201. - Y2) / P) + X2; 384 : LP150=LP150 ET (PROG X3); 385 : LE150=LE150 ET (PROG 201.); 386 : 387 : INT150=EVOL ROUG MANU LP150 LE150; 388 : * ---------------------------------------------------------------- 389 : * RECHERCHE DE LA COURBE DE TRAVAIL 390 : * ---------------------------------------------------------------- 391 : TAB=TABLE; 392 : J=15.;TAB.J=INT15; 393 : J=30.;TAB.J=INT30; 394 : J=45.;TAB.J=INT45; 395 : J=60.;TAB.J=INT60; 396 : J=75.;TAB.J=INT75; 397 : J=90.;TAB.J=INT90; 398 : J=120.;TAB.J=INT120; 399 : J=150.;TAB.J=INT150; 400 : NIND=8; 401 : PARAM=AN; 402 : COURBE=RESEAU TAB NIND PARAM; 403 : FINSI; 404 : ************************************************************************ 405 : * * 406 : * DETERMINATION DES AUTRES POINTS DE LA COURBE (M,PHI) : * 407 : * UTILISATION DES DIFFERENTES METHODES * 408 : * * 409 : ************************************************************************ 410 : I=2; 411 : 'REPETER' BOUCLE (NCOURB - 2); 412 : I=I + 1; 413 : SIG='EXTR' LSIG I; 414 : SIG0=SIGY + ((SIGF - SIGY) * (SIG - SIGY) / (SIGM - SIGY)); 415 : XM=XMSS0 * SIG0; 416 : * -------------------------------------------------------------------- 417 : * -------------------------------------------------------------------- 418 : 419 : * METHODE DE LA DEFORMATION DE REFERENCE 420 : 421 : 'SI' ('EGA' METHODE 'DEFR'); 422 : 423 : LPAN='EXTR' COURBE ABSC; 424 : LEAN='EXTR' COURBE ORDO; 425 : * ---------------------------------------------------------------- 426 : EPS='EXTR' LEPS I; 427 : PHIELC=CM * XM / 2.D0; 428 : EPSEL=SIG / E; 429 : EPSPLSEL=EPS / EPSEL - 1.; 430 : 'SI' (EPSPLSEL < 0.); 431 : PHI=PHIELC; 432 : 'SINON'; 433 : PHIPLSEL='IPOL' EPSPLSEL LEAN LPAN; 434 : PHI=(PHIPLSEL + 1.) * PHIELC; 435 : 'FINSI'; 436 : PHI=PHI * 2.D0; 437 : * -------------------------------------------------------------------- 438 : 'FINSI'; 439 : * -------------------------------------------------------------------- 440 : * -------------------------------------------------------------------- 441 : * AUTRES METHODES 442 : 443 : EPS='IPOL' SIG0 LSIG LEPS; 444 : 445 : * METHODES TADA - PARIS, LBB-NRC, LBB1 446 : 447 : 'SI' ('OU' ('OU' ('EGA' METHODE 'TADA') ('EGA' METHODE 'LBBNRC')) 448 : ('EGA' METHODE 'LBB1')); 449 : S=SIG0 / SIGF * SP; 450 : THETA=THETA0; 451 : K=0; 452 : 'REPETER' BOUC2 200; 453 : K=K + 1; 454 : TSP=THETA / PI; 455 : FZM=1.D0 + (AXX * ((4.6D0 * (TSP ** 1.5D0)) 456 : + (2.64D0 * (TSP ** 4.24D0)))); 457 : DFZM=AXX * ((6.9D0 * (TSP ** 0.5D0)) 458 : + (2.64D0 * 4.24D0 * (TSP ** 3.24D0))); 459 : HZM=2.D0 * TSP * DFZM; 460 : THETAF=(THETA0 - ((S ** 2) * THETA * FZM * HZM / BETA)) 461 : / (1.D0 - ((S ** 2) * FZM * (FZM + HZM) / BETA)); 462 : 'SI' (ABS(((THETAF-THETA) / THETA)) > 1.D-5); 463 : THETA=THETAF; 464 : 'SINON'; 465 : 'QUITTER' BOUC2; 466 : 'FINSI'; 467 : 'FIN' BOUC2; 468 : 'SI' (K 'EGA' 200); 469 : 'MESS' 'PAS DE CONVERGENCE POUR LA DETERMINATION DE THETAF'; 470 : 'QUITTER' TRACTUFI; 471 : 'FINSI'; 472 : TSP=THETAF / PI; 473 : TSP1=TSP ** 1.5D0; 474 : TSP2=TSP ** 4.24D0; 475 : TSP3=TSP ** 3.D0; 476 : TSP4=TSP ** 5.74D0; 477 : TSP5=TSP ** 8.48D0; 478 : FM=AXX * ((5.26D0 * TSP1) + (1.69D0 * TSP2)); 479 : FM=FM + 1.D0; 480 : FM=FM + 481 : (AXX * AXX * ((8.46D0 * TSP3) + (6.28D0 * TSP4) + (1.33D0 * TSP5))); 482 : CMEFF=2.D0 * THETA * THETA * FM; 483 : CMEFF=CMEFF / E / PI / RAYOM / RAYOM / EPAI; 484 : PHI=CMEFF * XM; 485 : * ----------------------------------------------------------------- 486 : 'SI' ('EGA' METHODE 'LBBNRC'); 487 : COEF=ALFA * (S ** (N - 1.D0)); 488 : PHI=(COEF + 1.D0) * PHI; 489 : 490 : 'FINSI'; 491 : * ----------------------------------------------------------------- 492 : 'SI' ('EGA' METHODE 'LBB1'); 493 : EPSE=SIG0 / E; 494 : EPSP=EPS - EPSE; 495 : COEF=EPSP / EPSE; 496 : PHI=(COEF + 1.D0) * PHI; 497 : 'FINSI'; 498 : * ----------------------------------------------------------------- 499 : 'FINSI'; 500 : * -------------------------------------------------------------------- 501 : * -------------------------------------------------------------------- 502 : 503 : * METHODES LBB2 504 : 505 : 'SI' ('EGA' METHODE 'LBB2'); 506 : PHI=CM * XM; 507 : EPSE=SIG0 / E; 508 : EPSP=EPS - EPSE; 509 : COEF=EPSP / EPSE; 510 : PHI=(COEF + 1.D0) * PHI; 511 : 'FINSI'; 512 : * -------------------------------------------------------------------- 513 : 'REMP' PM I XM; 514 : 'REMP' PPHI I PHI; 515 : 'FIN' BOUCLE; 516 : 'SI' ('EXISTE' ETAB 'COUL');COUL=ETAB.'COUL'; 517 : MP='EVOL' COUL MANU 'ROTATION' PPHI 'MOMENT' PM; 518 : 'SINON'; 519 : MP='EVOL' MANU 'ROTATION' PPHI 'MOMENT' PM; 520 : 'FINSI'; 521 : * 522 : 'FINPROC' MP CM KF;
© Cast3M 2003 - All rights reserved.
Disclaimer