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