Télécharger la plaquette Cast3M


fr en

Télécharger tractufi.procedur

Retour à la liste

Numérotation des lignes : 

  1. * TRACTUFI  PROCEDUR  CB215821  15/12/02    21:15:16     8719           
  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. RAYOM=REXT - (EPAI / 2.D0);                                                     
  65. THETA=AN / 2.D0;                                                                
  66. XMSS0=4.D0 * RAYOM * RAYOM * EPAI *                                             
  67.  (('COS' (THETA / 2.D0)) - (0.5D0 * ('SIN' THETA)));                            
  68. THETAD=THETA;                                                                   
  69. THETA=THETA * PI / 180.D0;                                                      
  70.  
  71. LEPS='EXTR' COUTRA ABSC;                                                        
  72. LSIG='EXTR' COUTRA ORDO;                                                        
  73. NCOURB='DIME' LSIG;                                                             
  74. PM='PROG' NCOURB * 0.;                                                          
  75. PPHI='PROG' NCOURB * 0.;                                                        
  76. SIGM='EXTR' LSIG NCOURB;                                                        
  77.  
  78. 'SI' ('EXISTE' ETAB 'SIGF');                                                    
  79.     SIGF= ETAB . 'SIGF';                                                        
  80.   'SINON';                                                                      
  81.     SIG1 = ETAB . 'SIG1';                                                       
  82.     SIGF= (SIG1 + SIGM) / 2.;                                                   
  83. 'FINSI';                                                                        
  84. SIGY='EXTR' LSIG 2;                                                             
  85. ************************************************************************        
  86. *                                                                      *        
  87. * DETERMINATION DE LA COMPLAISANCE  CM  DE L'ELEMENT ET DU COEFFICIENT *        
  88. * KF                                                                   *        
  89. ************************************************************************        
  90. RST=RAYOM / EPAI;                                                               
  91. 'SI' (RST < 10.D0) ;                                                            
  92.    AXX=(0.125D0 * RST - 0.25D0) ** 0.25D0;                                      
  93.   'SINON';                                                                      
  94.    AXX=(0.4D0 * RST - 3.D0) ** 0.25D0;                                          
  95. 'FINSI';                                                                        
  96. TSP=THETA / PI;                                                                 
  97. TSP1=TSP ** 1.5D0;                                                              
  98. TSP2=TSP ** 4.24D0;                                                             
  99. TSP3=TSP ** 3.D0;                                                               
  100. TSP4=TSP ** 5.74D0;                                                             
  101. TSP5=TSP ** 8.48D0;                                                             
  102. FM=AXX * ((5.26D0 * TSP1) + (1.69D0 * TSP2));                                   
  103. FM=FM + 1.D0;                                                                   
  104. FM=FM +                                                                         
  105.   (AXX * AXX * ((8.46D0 * TSP3) + (6.28D0 * TSP4) + (1.33D0*TSP5)));            
  106. CM=2.D0 * THETA * THETA * FM;                                                   
  107. CM=CM / E / PI / RAYOM / RAYOM / EPAI;                                          
  108.  
  109. 'SI' (THETAD < 100.);                                                           
  110.       FF = 1. + (AXX * ((4.6 * (TSP ** 1.5)) +                                  
  111.                                  (2.64 * (TSP ** 4.24))));                      
  112.  
  113. 'SINON';                                                                        
  114.       'MESS' 'ANGLE DE FISSURE TROP GRAND';                                     
  115. 'FINSI';                                                                        
  116.  K1 = PI * RAYOM * THETA;                                                       
  117.  K2 = K1 ** 0.5;                                                                
  118.  K3 = K2 * FF;                                                                  
  119.  KF = K3 / PI / RAYOM / RAYOM / EPAI;                                           
  120.  
  121.  
  122.  
  123. ************************************************************************        
  124. *                                                                      *        
  125. *  DETERMINATION DES COEFFICIENTS DE LA LOI DE RAMBERG - OSGOOD        *        
  126. *  POUR LA METHODE LBB-NRC                                             *        
  127. *                                                                      *        
  128. ************************************************************************        
  129. 'SI' (('EXISTE' ETAB 'ALFA') ET ('EXISTE' ETAB 'N'));                           
  130.           ALFA = ETAB.'ALFA';                                                   
  131.           N = ETAB.'N';                                                         
  132. 'SINON';                                                                        
  133.    'SI' ('EGA' METHODE 'LBBNRC');                                               
  134.        EPSF='IPOL' SIGF LSIG LEPS;                                              
  135.        EPSY='EXTR' LEPS 2;SIGY='EXTR' LSIG 2;                                   
  136.        ALFA=EPSF / EPSY - 1.D0;                                                 
  137.        I=1;                                                                     
  138.        'REPETER' BOUCR (NCOURB - 1);                                            
  139.          I=I + 1;                                                               
  140.          SIG2='EXTR' LSIG I;                                                    
  141.          'SI' (SIGF '&lt;EG' SIG2);                                                
  142.             SIG1='EXTR' LSIG (I - 1);                                           
  143.             EPS2=('EXTR' LEPS I) - (('EXTR' LEPS I) / E);                       
  144.             EPS1=('EXTR' LEPS (I - 1)) - (('EXTR' LEPS (I - 1)) / E);           
  145.             PENTE=(EPS2 - EPS1) / (SIG2 - SIG1);                                
  146.             'QUITTER' BOUCR;                                                    
  147.          'FINSI';                                                               
  148.        'FIN' BOUCR;                                                             
  149.        N=(SIGF * PENTE / ALFA / EPSY) + ((SIGF / SIGY - 1.D0) / ALFA);          
  150.        'SI' (('VALE' IMPI) 'EGA' 1);                                            
  151.           I=1;                                                                  
  152.           'REPETER' BOUCLE (NCOURB - 1);                                        
  153.               I=I + 1;                                                          
  154.               VSIG='EXTR' LSIG I;                                               
  155.               RSIG=VSIG / SIGF;                                                 
  156.               REPS=RSIG + (ALFA * (RSIG ** N));                                 
  157.               VEPS=EPSY * REPS;                                                 
  158.               'REMP' PM I VSIG;'REMP' PPHI I VEPS;                              
  159.           'FIN' BOUCLE;                                                         
  160.           COURAM='EVOL' ROUG MANU 'EPS' PPHI 'SIG' PM;                          
  161.           'DESS' (COUTRA ET COURAM) MARQ CARR ETOI;                             
  162.        'FINSI';                                                                 
  163.    'FINSI';                                                                     
  164. 'FINSI';                                                                        
  165. ************************************************************************        
  166. *                                                                      *        
  167. *   DETERMINATION DU POINT CORRESPONDANT A LA LIMITE ELASTIQUE         *        
  168. *                                                                      *        
  169. ************************************************************************        
  170. SIGY='EXTR' LSIG 2;                                                             
  171. MY=XMSS0 * SIGY;                                                                
  172. PHIY=CM * MY;                                                                   
  173. 'REMP' PM 2 MY;'REMP' PPHI 2 PHIY;                                              
  174. ************************************************************************        
  175. *                                                                      *        
  176. * POUR LA METHODE TADA-PARIS DETERMINATION DE THELIM ET BETA           *        
  177. *                                                                      *        
  178. ************************************************************************        
  179. 'SI' ('OU' ('OU' ('EGA' METHODE 'TADA') ('EGA' METHODE 'LBBNRC'))               
  180.  ('EGA' METHODE 'LBB1'));                                                       
  181.     SP=4. * (('COS' (THETAD / 2.D0)) - (0.5D0 * ('SIN' THETAD))) / PI;          
  182.     THETA0=THETA;                                                               
  183.     K=0;                                                                        
  184.    'REPETER' BOUC1 200;                                                         
  185.        K=K + 1;                                                                 
  186.        TSP=THETA / PI;                                                          
  187.        FZM=1.D0 + (AXX * ((4.6D0 * (TSP ** 1.5D0))                              
  188.           + (2.64D0 * (TSP ** 4.24D0))));                                       
  189.        DFZM=AXX * ((6.9D0 * (TSP ** 0.5D0))                                     
  190.           + (2.64D0 * 4.24D0 * (TSP ** 3.24D0)));                               
  191.        THELIM=THETA0 + ((THETA * FZM) / (FZM + (2.D0 * TSP * DFZM)));           
  192.       'SI' (((THELIM - THETA) / THETA) > 1.D-5);                                
  193.          THETA=THELIM;                                                          
  194.         'SINON';                                                                
  195.          'QUITTER' BOUC1;                                                       
  196.       'FINSI';                                                                  
  197.    'FIN' BOUC1;                                                                 
  198.    'SI' (K 'EGA' 200);                                                          
  199.       'MESS' 'PAS DE CONVERGENCE POUR LA DETERMINATION DE THELIM';              
  200.       'QUITTER' TRACTUFI;                                                       
  201.    'FINSI';                                                                     
  202.     THETAY=THELIM - THETA0;                                                     
  203.     TSP=THELIM / PI;                                                            
  204.     FZM=1.D0 + (AXX * ((4.6D0 * (TSP ** 1.5D0))                                 
  205.        + (2.64D0 * (TSP ** 4.24D0))));                                          
  206.     GZM=THELIM * FZM * FZM;                                                     
  207.     DGZM=GZM / THETAY;                                                          
  208.     BETA=(SP ** 2) * DGZM;                                                      
  209. 'FINSI';                                                                        
  210. ************************************************************************        
  211. *                                                                      *        
  212. * POUR LA METHODE DEFR    DETERMINATION DE LA COURBE PHIPLSEL-EPSPLSEL *        
  213. *                                                                      *        
  214. ************************************************************************        
  215.    'SI' ('EGA' METHODE 'DEFR');                                                 
  216. *       ----------------------------------------------------------------        
  217.  
  218. *       COURBES (EPSPL/EL,PHIPL/EL)                                             
  219.  
  220. 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          
  221.  15.07 16.93 18.12 20.77 22.14 23.11 25.08 27.78 29.35 30.53 35.14 35.91        
  222.  36.51 37.28 39.21 39.57 40.34;                                                 
  223. LP15=LP15 * 100.;                                                               
  224.  
  225. 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        
  226.  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          
  227.  9.43 9.79;                                                                     
  228. LE15=LE15 * 10.;                                                                
  229.  
  230. N=DIME LP15;                                                                    
  231. X1=EXTR LP15 (N-1);                                                             
  232. X2=EXTR LP15 N;                                                                 
  233. Y1=EXTR LE15 (N-1);                                                             
  234. Y2=EXTR LE15 N;                                                                 
  235. P=(Y2 - Y1) / (X2 - X1);                                                        
  236. X3=((201. - Y2) / P) + X2;                                                      
  237. LP15=LP15 ET (PROG X3);                                                         
  238. LE15=LE15 ET (PROG 201.);                                                       
  239.  
  240. INT15=EVOL MANU 'PHIPLSEL' LP15 'EPSPLSEL' LE15;                                
  241. *                                                                               
  242. 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        
  243.  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             
  244.  12.12 12.95 13.46 13.80 14.43 14.66 14.89 15.12;                               
  245. LP30=LP30 * 100.;                                                               
  246.  
  247. 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        
  248.  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        
  249.  13.75 14.61 15.22 15.62 16.59 17.24 17.89 18.54;                               
  250. LE30=LE30 * 10.;                                                                
  251.  
  252. N=DIME LP30;                                                                    
  253. X1=EXTR LP30 (N-1);                                                             
  254. X2=EXTR LP30 N;                                                                 
  255. Y1=EXTR LE30 (N-1);                                                             
  256. Y2=EXTR LE30 N;                                                                 
  257. P=(Y2 - Y1) / (X2 - X1);                                                        
  258. X3=((201. - Y2) / P) + X2;                                                      
  259. LP30=LP30 ET (PROG X3);                                                         
  260. LE30=LE30 ET (PROG 201.);                                                       
  261.  
  262. INT30=EVOL ROUG MANU LP30 LE30;                                                 
  263. *                                                                               
  264. 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        
  265.  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          
  266.  3.92 4.04 4.14 4.19 4.30 4.36 4.41 4.49 4.54;                                  
  267. LP45=LP45 * 100.;                                                               
  268.  
  269. 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        
  270.  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          
  271.  9.98 10.03 10.19 10.28 10.47 10.68 10.83 11.20 11.45;                          
  272. LE45=LE45 * 10.;                                                                
  273.  
  274. N=DIME LP45;                                                                    
  275. X1=EXTR LP45 (N-1);                                                             
  276. X2=EXTR LP45 N;                                                                 
  277. Y1=EXTR LE45 (N-1);                                                             
  278. Y2=EXTR LE45 N;                                                                 
  279. P=(Y2 - Y1) / (X2 - X1);                                                        
  280. X3=((201. - Y2) / P) + X2;                                                      
  281. LP45=LP45 ET (PROG X3);                                                         
  282. LE45=LE45 ET (PROG 201.);                                                       
  283.  
  284. INT45=EVOL ROUG MANU LP45 LE45;                                                 
  285. *                                                                               
  286. 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        
  287.  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        
  288.  11.26 11.96 12.41 12.86 13.26;                                                 
  289. LP60=LP60 * 10.;                                                                
  290.  
  291. 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        
  292.  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        
  293.  11.07 11.55 11.91 12.20 12.49;                                                 
  294. LE60=LE60 * 10.;                                                                
  295.  
  296. N=DIME LP60;                                                                    
  297. X1=EXTR LP60 (N-1);                                                             
  298. X2=EXTR LP60 N;                                                                 
  299. Y1=EXTR LE60 (N-1);                                                             
  300. Y2=EXTR LE60 N;                                                                 
  301. P=(Y2 - Y1) / (X2 - X1);                                                        
  302. X3=((201. - Y2) / P) + X2;                                                      
  303. LP60=LP60 ET (PROG X3);                                                         
  304. LE60=LE60 ET (PROG 201.);                                                       
  305.  
  306. INT60=EVOL ROUG MANU LP60 LE60;                                                 
  307. *                                                                               
  308. 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         
  309.  4.12 4.50 4.87 5.25 5.50 5.96 6.28 6.53 6.79 6.98 7.18;                        
  310. LP75=LP75 * 10.;                                                                
  311.  
  312. 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        
  313.  6.16 6.54 6.99 7.38 7.67 8.14 8.51 8.99 9.53 10.06 11.34;                      
  314. LE75=LE75 * 10.;                                                                
  315.  
  316. N=DIME LP75;                                                                    
  317. X1=EXTR LP75 (N-1);                                                             
  318. X2=EXTR LP75 N;                                                                 
  319. Y1=EXTR LE75 (N-1);                                                             
  320. Y2=EXTR LE75 N;                                                                 
  321. P=(Y2 - Y1) / (X2 - X1);                                                        
  322. X3=((201. - Y2) / P) + X2;                                                      
  323. LP75=LP75 ET (PROG X3);                                                         
  324. LE75=LE75 ET (PROG 201.);                                                       
  325.  
  326. INT75=EVOL ROUG MANU LP75 LE75;                                                 
  327. *                                                                               
  328. 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        
  329.  1.68 1.77 1.90 1.99 2.10 2.28 2.37 2.45 2.57 2.79 3.04;                        
  330. LP90=LP90 * 10.;                                                                
  331.  
  332. 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        
  333.  6.65 6.81 7.02 7.17 7.40 7.84 8.06 8.31 8.72 9.44 10.27;                       
  334. LE90=LE90 * 10.;                                                                
  335.  
  336. N=DIME LP90;                                                                    
  337. X1=EXTR LP90 (N-1);                                                             
  338. X2=EXTR LP90 N;                                                                 
  339. Y1=EXTR LE90 (N-1);                                                             
  340. Y2=EXTR LE90 N;                                                                 
  341. P=(Y2 - Y1) / (X2 - X1);                                                        
  342. X3=((201. - Y2) / P) + X2;                                                      
  343. LP90=LP90 ET (PROG X3);                                                         
  344. LE90=LE90 ET (PROG 201.);                                                       
  345.  
  346. INT90=EVOL ROUG MANU LP90 LE90;                                                 
  347. *                                                                               
  348. 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            
  349.  7.83 8.48 9.11 10.41 11.83 12.32 12.90 13.47 14.70 15.63 16.58;                
  350.  
  351. 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            
  352.  5.62 5.83 6.04 6.68 7.18 7.44 7.67 8.00 8.60 9.09 10.19;                       
  353. LE120=LE120 * 10.;                                                              
  354.  
  355. N=DIME LP120;                                                                   
  356. X1=EXTR LP120 (N-1);                                                            
  357. X2=EXTR LP120 N;                                                                
  358. Y1=EXTR LE120 (N-1);                                                            
  359. Y2=EXTR LE120 N;                                                                
  360. P=(Y2 - Y1) / (X2 - X1);                                                        
  361. X3=((201. - Y2) / P) + X2;                                                      
  362. LP120=LP120 ET (PROG X3);                                                       
  363. LE120=LE120 ET (PROG 201.);                                                     
  364.  
  365. INT120=EVOL ROUG MANU LP120 LE120;                                              
  366. *                                                                               
  367. 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          
  368.  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         
  369.  14.38;                                                                         
  370.  
  371. 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           
  372.  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         
  373.  15.20;                                                                         
  374.  LE150=LE150 * 10.;                                                             
  375.  
  376. N=DIME LP150;                                                                   
  377. X1=EXTR LP150 (N-1);                                                            
  378. X2=EXTR LP150 N;                                                                
  379. Y1=EXTR LE150 (N-1);                                                            
  380. Y2=EXTR LE150 N;                                                                
  381. P=(Y2 - Y1) / (X2 - X1);                                                        
  382. X3=((201. - Y2) / P) + X2;                                                      
  383. LP150=LP150 ET (PROG X3);                                                       
  384. LE150=LE150 ET (PROG 201.);                                                     
  385.  
  386. INT150=EVOL ROUG MANU LP150 LE150;                                              
  387. *       ----------------------------------------------------------------        
  388. * RECHERCHE DE LA COURBE DE TRAVAIL                                             
  389. *       ----------------------------------------------------------------        
  390. TAB=TABLE;                                                                      
  391. J=15.;TAB.J=INT15;                                                              
  392. J=30.;TAB.J=INT30;                                                              
  393. J=45.;TAB.J=INT45;                                                              
  394. J=60.;TAB.J=INT60;                                                              
  395. J=75.;TAB.J=INT75;                                                              
  396. J=90.;TAB.J=INT90;                                                              
  397. J=120.;TAB.J=INT120;                                                            
  398. J=150.;TAB.J=INT150;                                                            
  399. NIND=8;                                                                         
  400. PARAM=AN;                                                                       
  401. COURBE=RESEAU TAB NIND PARAM;                                                   
  402. FINSI;                                                                          
  403. ************************************************************************        
  404. *                                                                      *        
  405. *   DETERMINATION DES AUTRES POINTS DE LA COURBE (M,PHI) :             *        
  406. *   UTILISATION DES DIFFERENTES METHODES                               *        
  407. *                                                                      *        
  408. ************************************************************************        
  409. I=2;                                                                            
  410. 'REPETER' BOUCLE (NCOURB - 2);                                                  
  411.     I=I + 1;                                                                    
  412.     SIG='EXTR' LSIG I;                                                          
  413.     SIG0=SIGY + ((SIGF - SIGY) * (SIG - SIGY) / (SIGM - SIGY));                 
  414.     XM=XMSS0 * SIG0;                                                            
  415. *   --------------------------------------------------------------------        
  416. *   --------------------------------------------------------------------        
  417.  
  418. *   METHODE DE LA DEFORMATION DE REFERENCE                                      
  419.  
  420.    'SI' ('EGA' METHODE 'DEFR');                                                 
  421.  
  422. LPAN='EXTR' COURBE ABSC;                                                        
  423. LEAN='EXTR' COURBE ORDO;                                                        
  424. *       ----------------------------------------------------------------        
  425.    EPS='EXTR' LEPS I;                                                           
  426.    PHIELC=CM * XM / 2.D0;                                                       
  427.    EPSEL=SIG / E;                                                               
  428.    EPSPLSEL=EPS / EPSEL - 1.;                                                   
  429.    'SI' (EPSPLSEL < 0.);                                                        
  430.         PHI=PHIELC;                                                             
  431.    'SINON';                                                                     
  432.         PHIPLSEL='IPOL' EPSPLSEL LEAN LPAN;                                     
  433.         PHI=(PHIPLSEL + 1.) * PHIELC;                                           
  434.    'FINSI';                                                                     
  435.     PHI=PHI * 2.D0;                                                             
  436. *   --------------------------------------------------------------------        
  437.    'FINSI';                                                                     
  438. *   --------------------------------------------------------------------        
  439. *   --------------------------------------------------------------------        
  440. *   AUTRES METHODES                                                             
  441.  
  442.     EPS='IPOL' SIG0 LSIG LEPS;                                                  
  443.  
  444. *   METHODES TADA - PARIS, LBB-NRC, LBB1                                        
  445.  
  446.    'SI' ('OU' ('OU' ('EGA' METHODE 'TADA') ('EGA' METHODE 'LBBNRC'))            
  447.     ('EGA' METHODE 'LBB1'));                                                    
  448.        S=SIG0 / SIGF * SP;                                                      
  449.        THETA=THETA0;                                                            
  450.        K=0;                                                                     
  451.       'REPETER' BOUC2 200;                                                      
  452.           K=K + 1;                                                              
  453.           TSP=THETA / PI;                                                       
  454.           FZM=1.D0 + (AXX * ((4.6D0 * (TSP ** 1.5D0))                           
  455.              + (2.64D0 * (TSP ** 4.24D0))));                                    
  456.           DFZM=AXX * ((6.9D0 * (TSP ** 0.5D0))                                  
  457.              + (2.64D0 * 4.24D0 * (TSP ** 3.24D0)));                            
  458.           HZM=2.D0 * TSP * DFZM;                                                
  459.           THETAF=(THETA0 - ((S ** 2) * THETA * FZM * HZM / BETA))               
  460.                 / (1.D0 - ((S ** 2) * FZM * (FZM + HZM) / BETA));               
  461.          'SI' (ABS(((THETAF-THETA) / THETA)) > 1.D-5);                          
  462.              THETA=THETAF;                                                      
  463.             'SINON';                                                            
  464.             'QUITTER' BOUC2;                                                    
  465.          'FINSI';                                                               
  466.       'FIN' BOUC2;                                                              
  467.       'SI' (K 'EGA' 200);                                                       
  468.          'MESS' 'PAS DE CONVERGENCE POUR LA DETERMINATION DE THETAF';           
  469.          'QUITTER' TRACTUFI;                                                    
  470.       'FINSI';                                                                  
  471.        TSP=THETAF / PI;                                                         
  472.        TSP1=TSP ** 1.5D0;                                                       
  473.        TSP2=TSP ** 4.24D0;                                                      
  474.        TSP3=TSP ** 3.D0;                                                        
  475.        TSP4=TSP ** 5.74D0;                                                      
  476.        TSP5=TSP ** 8.48D0;                                                      
  477.        FM=AXX * ((5.26D0 * TSP1) + (1.69D0 * TSP2));                            
  478.        FM=FM + 1.D0;                                                            
  479.        FM=FM +                                                                  
  480.    (AXX * AXX * ((8.46D0 * TSP3) + (6.28D0 * TSP4) + (1.33D0 * TSP5)));         
  481.        CMEFF=2.D0 * THETA * THETA * FM;                                         
  482.        CMEFF=CMEFF / E / PI / RAYOM / RAYOM / EPAI;                             
  483.        PHI=CMEFF * XM;                                                          
  484. *      -----------------------------------------------------------------        
  485.       'SI' ('EGA' METHODE 'LBBNRC');                                            
  486.           COEF=ALFA * (S ** (N - 1.D0));                                        
  487.           PHI=(COEF + 1.D0) * PHI;                                              
  488.  
  489.       'FINSI';                                                                  
  490. *      -----------------------------------------------------------------        
  491.       'SI' ('EGA' METHODE 'LBB1');                                              
  492.           EPSE=SIG0 / E;                                                        
  493.           EPSP=EPS - EPSE;                                                      
  494.           COEF=EPSP / EPSE;                                                     
  495.           PHI=(COEF + 1.D0) * PHI;                                              
  496.       'FINSI';                                                                  
  497. *      -----------------------------------------------------------------        
  498.    'FINSI';                                                                     
  499. *   --------------------------------------------------------------------        
  500. *   --------------------------------------------------------------------        
  501.  
  502. *   METHODES LBB2                                                               
  503.  
  504.    'SI' ('EGA' METHODE 'LBB2');                                                 
  505.        PHI=CM * XM;                                                             
  506.        EPSE=SIG0 / E;                                                           
  507.        EPSP=EPS - EPSE;                                                         
  508.        COEF=EPSP / EPSE;                                                        
  509.        PHI=(COEF + 1.D0) * PHI;                                                 
  510.    'FINSI';                                                                     
  511. *   --------------------------------------------------------------------        
  512.    'REMP' PM I XM;                                                              
  513.    'REMP' PPHI I PHI;                                                           
  514. 'FIN' BOUCLE;                                                                   
  515. 'SI' ('EXISTE' ETAB 'COUL');COUL=ETAB.'COUL';                                   
  516.     MP='EVOL' COUL MANU 'ROTATION' PPHI 'MOMENT' PM;                            
  517.    'SINON';                                                                     
  518.     MP='EVOL' MANU 'ROTATION' PPHI 'MOMENT' PM;                                 
  519. 'FINSI';                                                                        
  520. *                                                                               
  521. 'FINPROC' MP CM KF;                                                             
  522.  
  523.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales