Télécharger pasapas.procedur

Retour à la liste

Numérotation des lignes :

  1. * PASAPAS PROCEDUR PASCAL 22/04/27 21:15:09 11356
  2.  
  3. 'DEBP' PASAPAS PRECED*'TABLE';
  4. 'SAUT' 1 'LIGNE';
  5. '------------------ DEBUT DE LA PROCEDURE PASAPAS --------------------';
  6. ************************************************************************
  7. * *
  8. * *
  9. *** PRECED objet de type TABLE *
  10. * C'est une table donnant , les arguments d'entrees, les *
  11. * options, puis les resultats des pas precedents *
  12. * *
  13. * Etat initial au debut du pas de temps decrit par le MCHELM WTAB.ETAT1*
  14. * champs de materiaux : *
  15. * WTAB.'MAT_REF' materiau etat a la temperature de REFERENCE *
  16. * WTAB.'MAT0' materiau etat debut de calcul *
  17. * WTAB.'MAT1' materiau fin du pas de temps *
  18. * *
  19. ************************************************************************
  20.  
  21. **************** I PARTIE ******************************************
  22.  
  23. *---------- attribution des valeurs par defaut et verifs elementaires
  24. PAS_DEFA PRECED ;
  25. *-------- initialisation du pas precedent si besoin et creation des tables de sorties.
  26. PAS_INIT PRECED ;
  27.  
  28. WTAB = PRECED.'WTABLE' ;
  29. conti = PRECED.'CONTINUATION';
  30. estim = PRECED.'ESTIMATION' ;
  31. *
  32. WTAB.'GE0_DEB'=WTAB.'FOR0';
  33. 'SI'( WTAB.'GRANDS_DEPLACEMENTS' 'OU' WTAB.'FEFP_FORMULATION');
  34. 'SI' ( 'EXIS' conti 'DEPLACEMENTS');
  35. WTAB.'GE0_DEB'= 'FORME' conti . 'DEPLACEMENTS';
  36. 'FINSI';
  37. 'FINSI';
  38. 'FORM' WTAB.'FOR0';
  39.  
  40. *--------- quel est le delta-temps du dernier pas calcule ? ---------
  41. 'SI' ('EGA' WTAB.'DT' ('MOT' 'INCONNU'));
  42. DELTAN = 0.D0 ;
  43. WTAB.'DT' = DELTAN ;
  44. 'SINON';
  45. DELTAN= WTAB.'DT';
  46. 'FINSI';
  47.  
  48. *------------Appel a une procedure utilisateur PERSO1
  49. 'SI' WTAB.'PROCEDURE_PERSO1';
  50. PRECED.'PERSO1_APPEL' = 1 ;
  51. PERSO1 PRECED ;
  52. 'SI' ('EXIS' PRECED 'ARRET') ;
  53. WTAB.'ARRET' = PRECED.'ARRET' ;
  54. 'FINSI' ;
  55. 'FINSI';
  56.  
  57. * ------------Appel a une procedure utilisateur REEV_MEC
  58. 'SI' WTAB.'PROCEDURE_REEV_MEC' ;
  59. REEV_MEC PRECED 0 ;
  60. 'FINSI' ;
  61.  
  62. * ------------Appel a une procedure utilisateur REEV_THE
  63. 'SI' WTAB.'PROCEDURE_REEV_THE' ;
  64. REEV_THE PRECED 0 ;
  65. 'FINSI' ;
  66.  
  67. *--------- initialisation ETAT ET MATERIAUX
  68. WTAB.'ETAT0' = PAS_ETAT PRECED WTAB.'TEMPS0' ;
  69. WTAB.'ETAT1' = WTAB.'ETAT0' ;
  70. WTAB.'MAT0' = PAS_MATE PRECED WTAB.'ETAT0' ;
  71. WTAB.'MAT1' = WTAB.'MAT0' ;
  72.  
  73. *--------- initialisation deformation thermique de reference
  74. 'SI' WTAB.'FOR_MECA' ;
  75. 'SI' (WTAB.'CHAR_THE' 'OU' WTAB.'FOR_THER' ) ;
  76. 'SI' ('EGA' WTAB.'ETREF' 'INCONNU');
  77. WTAB.'ETAT_REF'= PAS_ETAT PRECED WTAB.'TEMPS0' 'REFERENCE';
  78. WTAB.'MAT_REF' = PAS_MATE PRECED WTAB.'ETAT_REF' ;
  79. WTAB.'ETREF' ='EPTH' WTAB.'MO_TOT' WTAB.'MAT_REF' WTAB.'TEMPERATURE_REFERENCE' ;
  80. 'FINS' ;
  81.  
  82. 'SI' WTAB.'POR1';
  83. 'SI' ('EGA' ('TYPE' WTAB.'TET1') 'CHPOINT');
  84. TET1 ='CHAN' 'CHAM' WTAB.'TET1' WTAB.'MO_TOT' 'STRESSES' 'TEMPERATURES' ;
  85. 'FINS';
  86. WTAB.'MSRTHER1'= PAS_EPTH PRECED WTAB.'MO_TOT' WTAB.'MAT0' TET1 ;
  87. 'FINS' ;
  88. 'FINS' ;
  89. 'FINS' ;
  90.  
  91. ************************** II PARTIE *************************
  92. ******************************************************************
  93. ******************* boucle sur les pas de temps *****************
  94. ******************************************************************
  95.  
  96. * ------------quelques initialisations de l'algorithme
  97. PASFINAL = 0 ;
  98. IPAS = 0 ;
  99. ICHG = 0 ;
  100. IAUTO = 1 ;
  101. IREREDU = VRAI;
  102. TEMP0 = WTAB.'TEMPS0';
  103. ISAUV = VRAI;
  104. DTV =('EXTR' WTAB.'TEMPS_CALCULES' 1) - TEMP0 ;
  105. 'SI' ( DTV '<' 1.e-15 );
  106. 'SI' WTAB.'FREQUENTIEL' ;
  107. DTV = 1.e-14 ;
  108. 'SINON' ;
  109. DTV =('EXTR' WTAB.'TEMPS_CALCULES' 2) - TEMP0 ;
  110. 'FINSI';
  111. 'FINSI';
  112. WTAB.'ARRET' = FAUX;
  113. WTAB.'AUTODEUX'= FAUX;
  114. WTAB.'CONV' = VRAI;
  115. DT_AVANT = DTV ;
  116. WTAB.'NOMBRE_ITERATIONS'=12;
  117.  
  118. * Attention Bizarre : RFROT est initialise ici et utilise dans UNPAS ????
  119. RFROT=0;
  120.  
  121. *---------- debut boucle sur pas de temps
  122. 'REPETER' BEXTERN ;
  123.  
  124. *---recherche du temps a calculer TI
  125. NPAS='DIME' WTAB.'TEMPS_CALCULES';
  126. ICHG=ICHG + 1; WTAB.'ICHG'= ICHG;
  127. 'SI' ( ICHG '>' NPAS) ; 'QUITTER' BEXTERN; 'FINSI';
  128. TI = 'EXTRAIRE' WTAB.'TEMPS_CALCULES' ICHG ;
  129. TIV=TI;
  130. 'SI' (WTAB.'FREQUENTIEL') ;
  131. 'SI' (TI '<' TEMP0 ) ; 'ITERER' BEXTERN ; 'FINSI' ;
  132. 'SINON' ;
  133. 'SI' ((TI '<' TEMP0 ) 'OU' (TI 'EGA' PRECED.'TEMPS'. 0));
  134. 'ITERER' BEXTERN ;
  135. 'FINSI'; 'FINSI';
  136.  
  137. *---ajustement du pas en fonction du nombre d iteration du pas precedent
  138. PASFINAL=PASFINAL+1;
  139. 'SI' WTAB.'PAS_AJUSTE';
  140. DTV=DT_AVANT;
  141. 'SI' WTAB.'CONV';
  142. 'SI' (WTAB.'NOMBRE_ITERATIONS' > 20);
  143. DTV=DTV'/' (2. ** 0.5);
  144. 'FINSI';
  145. 'SI' (WTAB.'NOMBRE_ITERATIONS' < 7);
  146. 'SI' IREREDU;
  147. DTV=DTV '*' (2. ** 0.5);
  148. 'FINSI';
  149. 'FINSI';
  150. DT_AVANT= DTV;
  151. IREREDU=VRAI;
  152. 'SINON';
  153. DTV=0.0000000001D0* DT_AVANT;
  154. DT_AVANT= DT_AVANT /2.;
  155. IREREDU=FAUX;
  156. 'FINSI';
  157. TTI = DTV* 1.0000000001 + TEMP0;
  158. 'SI' ( TTI '&lt;EG' TIV );
  159. TI=TTI;
  160. PASFINAL=0;
  161. 'FINSI';
  162. 'FINSI';
  163. estim.'TEMPS'=TI;
  164. 'SI' (PASFINAL '>EG' 1.5) ; DT_AVANT= TI-TEMP0; 'FINSI';
  165. WTAB.'T_FINAL'= TI;
  166. WTAB.'DT_INIT' = TI-TEMP0;
  167. 'SI'((WTAB.'CONSOLIDATION' 'OU' WTAB.'DYNAMIQUE')
  168. 'ET' (WTAB.'DT_INIT' 'EGA' 0.));
  169. 'ITERER' BEXTERN;
  170. 'FINSI';
  171. WTAB.'DT'= WTAB.'DT_INIT';
  172.  
  173. *---compteur IPAS des pas de temps calcules
  174. IPAS = IPAS + 1;
  175. WTAB.'PAS'=IPAS;
  176.  
  177. *---Mise a jour du modele si chargement de type MODE :
  178. 'SI' (WTAB.'CHAR_MODE') ;
  179. MOD1P = WTAB.'MODELE' ;
  180. MODEVOL1 = 'TIRE' PRECED.'MODELE' 'MODE' TI ;
  181. WTAB.'MODELE' = MODEVOL1 ;
  182.  
  183. 'SI' (WTAB.'CHAR_MATE') ;
  184. MATEVOL1 = 'TIRE' PRECED.'CARACTERISTIQUES' 'MATE' TI ;
  185. WTAB.'CARACTERISTIQUES' = MATEVOL1 ;
  186. 'FINSI' ;
  187.  
  188. PAS_MODL PRECED ;
  189.  
  190. * MaJ table rayonnement :
  191. 'SI' WTAB.'RAYO' ;
  192. WTAB.'RAYONNEMENT' = 'PRRA' WTAB.'MOD_RAY' ;
  193. 'FINSI' ;
  194.  
  195. * On sauvegarde les modeles successifs pour analyse resultats :
  196. ESTIM.'MODELES' = WTAB.'MODELE' ;
  197.  
  198. * Initialisation TINI si thermique :
  199. 'SI' (WTAB.'THE1') ;
  200. MODX = WTAB.'MOD_COND' ;
  201. MODXP = MOD1P 'EXTR' 'MATE' 'CONDUCTION' ;
  202. GEO1 = MODX 'EXTR' 'MAIL' ;
  203. GEO1P = MODXP 'EXTR' 'MAIL' ;
  204. DGEO1 = GEO1 'DIFF' GEO1P ;
  205.  
  206. * Adapation du modele : transport de la solution
  207. 'SI' ((WTAB.'ADAP_MODE') 'ET' ('NON' (TI '<' (WTAB.'TPS1_ADAP')))) ;
  208. mess '***** Transport solution adaptation' ;
  209. IADAP1 = (WTAB.'INDI_ADAP') + 1 ;
  210. 'SI' (IADAP1 '&lt;EG' ('DIME' (WTAB.'LTPS_ADAP'))) ;
  211. WTAB.'TPS1_ADAP' = WTAB.'LTPS_ADAP' 'EXTR' IADAP1 ;
  212. WTAB.'INDI_ADAP' = IADAP1 ;
  213. 'SINON' ;
  214. WTAB.'TPS1_ADAP' = 'MAXI' WTAB.'TEMPS_CALCULES' + 1. ;
  215. 'FINSI' ;
  216. *erre 21 ;
  217. * Projection du champ de temperature sur la partie adaptee :
  218. DGEO1P = DGEO1 'INTE' GEO1P ;
  219. DGEO2P = (DGEO1 'INTE' GEO1) 'INCLU' DGEO1P 'VOLU' 'BARY' ;
  220. CHT1P = (WTAB.'TET1') 'EXCO' 'T' 'T' ;
  221. CHAMT1P = 'CHAN' 'CHAM' DGEO1P CHT1P ;
  222. CHT2P = 'PROI' CHAMT1P DGEO2P ;
  223. CHT1P = ('REDU' CHT1P (GEO1P 'DIFF' DGEO1P)) 'ET' CHT2P ;
  224. WTAB.'TET1' = CHT1P ;
  225.  
  226. * Calcul de la capacite courante sur la partie adaptee :
  227. DMODX = 'REDU' MODX DGEO2P ;
  228. DMATX = 'REDU' (WTAB.'CARACTERISTIQUES') DMODX ;
  229. DMATX = 'VARI' 'NUAG' DMODX DMATX CHT2P ;
  230. CAPAX = 'CAPA' DMODX DMATX ;
  231. CAPA1P = ('REDU' (WTAB.'CAPA_COURANTE') (GEO1P 'DIFF' DGEO1P)) 'ET' CAPAX ;
  232. WTAB.'CAPA_COURANTE' = CAPA1P ;
  233.  
  234. DGEO1 = (DGEO1 'INTE' GEO1) 'DIFF' DGEO2P ;
  235. 'FINSI' ;
  236.  
  237. * Calcul temperature dans partie ajoutee du modele :
  238. 'SI' ('NON' (('VIDE' DGEO1) 'OU' ('VIDE' (DGEO1 'INTE' GEO1)))) ;
  239. * Quantite de chaleur sur geometrie GEO1P :
  240. CHT1P = (WTAB.'TET1') 'EXCO' 'T' 'T' ;
  241. CHQ1P = (WTAB.'CAPA_COURANTE') * CHT1P ;
  242. * Quantite de chaleur dans ajout :
  243. DMODX = 'REDU' MODX DGEO1 ;
  244. DMATX = 'REDU' (WTAB.'CARACTERISTIQUES') DMODX ;
  245. CHTINIX = DMATX 'EXCO' 'TINI' 'T' ;
  246. CHTINIX = 'CHAN' 'CHPO' CHTINIX DMODX 'MOYE' ;
  247. DMATX = 'VARI' 'NUAG' DMODX DMATX CHTINIX ;
  248. CAPAX = 'CAPA' DMODX DMATX ;
  249. CHQX = CAPAX * CHTINIX ;
  250. * Calcul T sur GEO1 :
  251. CAPA1 = (WTAB.'CAPA_COURANTE') et CAPAX ;
  252. CHTX1 = 'MANU' 'CHPO' GEO1 1 'T' 1. 'NATURE' 'DIFFUS' ;
  253. CHQX1 = CAPA1 * CHTX1;
  254. CHT1 = (CHQ1P + CHQX) * ('INVE' CHQX1) ('MOTS' 'Q') ('MOTS' 'Q') ('MOTS' 'T') ;
  255. WTAB.'TET1' = CHT1 ;
  256. *List ((resu (CHQ1P + CHQX)) maxi abs) ;
  257. *List ((resu (CAPA1*CHT1)) maxi abs) ;
  258. *trac CHT1 GEO1 ;
  259. 'FINSI' ;
  260.  
  261. 'FINSI' ;
  262.  
  263. * Initialisation PROPORTIONS_PHASES si metallurgie :
  264. 'SI' WTAB.'FOR_METALLU' ;
  265. MODX = (WTAB.'MODELE') 'EXTR' 'FORM' 'METALLURGIE' ;
  266. MODXP = MOD1P 'EXTR' 'FORM' 'METALLURGIE' ;
  267. GEO1 = MODX 'EXTR' 'MAIL' ;
  268. GEO1P = MODXP 'EXTR' 'MAIL' ;
  269. DGEO1 = GEO1 'DIFF' GEO1P ;
  270. 'SI' ('NON' (('VIDE' DGEO1) 'OU' ('VIDE' (DGEO1 'INTE' GEO1)))) ;
  271. CHMET0 = ESTIM.'PROPORTIONS_PHASES' ;
  272. DMODX = 'REDU' MODX DGEO1 ;
  273. CHMET1 = 'ZERO' DMODX 'VARINTER' ;
  274. LISCO = 'EXTR' DMODX 'VARI' ;
  275. CMP1 = 'EXCO' CHMET1 ('EXTR' LISCO 1) ('EXTR' LISCO 1) ;
  276. CHMET1 = CHMET1 + (CMP1 + 1.D0);
  277. CHMET1 = CHMET0 'ET' CHMET1 ;
  278. ESTIM.'PROPORTIONS_PHASES' = CHMET1 'REDU' MODX ;
  279. *trac CHMET1 MODX ;
  280. 'FINSI' ;
  281. 'FINS' ;
  282.  
  283. 'FINSI' ;
  284.  
  285. *--- Mise a jour blocages passes dans chargement
  286. * Cas BLOCAGES_MECANIQUES (BLOM) :
  287. 'SI' (WTAB.'CHAR_BLOM') ;
  288. BLOM1 = 'TIRE' PRECED.'BLOCAGES_MECANIQUES' 'BLOM' TI ;
  289. WTAB.'BLOCAGES_MECANIQUES' = BLOM1 ;
  290.  
  291. * MAJ Operateur Resolution :
  292. 'SI' ('EXIS' WTAB 'RRRR') ;
  293. RRRR = WTAB.'RRRR' 'EXTR' 'RIGI' 'NOMU' ;
  294. WTAB.'RRRR' = RRRR 'ET' BLOM1 ;
  295. 'FINSI' ;
  296. 'FINSI' ;
  297.  
  298. * Cas BLOCAGES_THERMIQUES (BLOT) :
  299. 'SI' (WTAB.'CHAR_BLOT') ;
  300. BLOT1 = 'TIRE' PRECED.'BLOCAGES_THERMIQUES' 'BLOT' TI ;
  301. WTAB.'BLOCAGES_THERMIQUES' = BLOT1 ;
  302. 'FINSI' ;
  303.  
  304. * Cas BLOCAGES_DIFFUSIONS (BLOD) :
  305. 'SI' (WTAB.'CHAR_BLOD') ;
  306. BLOD1 = 'TIRE' PRECED.'BLOCAGES_DIFFUSIONS' 'BLOD' TI ;
  307. WTAB.'BLOCAGES_DIFFUSIONS' = BLOD1 ;
  308. 'FINSI' ;
  309.  
  310. *---on se place dans la configuration deformee
  311. 'SI'( WTAB.'GRANDS_DEPLACEMENTS' 'OU' WTAB.'FEFP_FORMULATION');
  312. 'SI' (WTAB.'RECALCUL');
  313. 'SI' ('EGA' ( 'TYPE' WTAB.'GE0_DEB') 'CONFIGUR');
  314. 'FORM' WTAB.'GE0_DEB';
  315. 'FINSI';
  316. 'FINSI';
  317. 'FINSI';
  318.  
  319. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  320. *++++++++++++++++++ BOUCLE interne de bouclage mec_ther ++++++++++++++
  321. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  322. I_BOTH = 0 ;
  323. 'REPETER' BO_BOTH WTAB.'BOU_METH' ;
  324. I_BOTH = I_BOTH + 1 ;
  325. *
  326. * re-initialisation de WTAB.'CHARGEMENT' en tete de boucle mec-ther
  327. *
  328. 'SI' WTAB.'CONVERGENCE_MEC_THE' ;
  329. WTAB.'CHARGEMENT' = WTAB . 'CHARGEMENT_I' ;
  330. 'FINSI';
  331.  
  332. *######CALCUL D'UN PAS SI IL Y A DE LA THERMIQUE ####################
  333. 'SI' (WTAB.'FOR_THER' 'OU' WTAB.'FOR_DIFF') ;
  334. 'SI' WTAB.'FOR_MECA' ;
  335. *---------- en cas de convergence mecanique thermique
  336. 'SI' (I_BOTH > 1) ;
  337. *---------- reprise de la configuration meca de la derniere estimation
  338. * du pas present si coherence mecanique/thermique
  339. 'SI' ('EXIS' estim 'FOR' );
  340. WTAB .'RECARI' = VRAI ;
  341. WTAB .'RECADET' = VRAI ;
  342. WTAB.'REA_GEOM' = VRAI ;
  343. 'FORM' estim.'FOR' ;
  344. 'FINSI';
  345. 'FINSI' ;
  346. 'FINSI' ;
  347.  
  348. 'SI' (TI 'NEG' TEMP0);
  349. 'SI' ('NON' ('EXISTE' PRECED 'PAS_MAX'));
  350. WTAB.'TEM_CALC' = 'PROG' TI;
  351. 'SINON';
  352. 'SI' ((TEMP0 + PRECED.'PAS_MAX') '>EG' TI);
  353. WTAB.'TEM_CALC' = 'PROG' TI;
  354. 'SINON';
  355. WTAB.'TEM_CALC' = 'PROG' TEMP0
  356. 'PAS' PRECED.'PAS_MAX' TI;
  357. WTAB.'TEM_CALC' = 'ENLE' WTAB.'TEM_CALC' 1;
  358. 'FINSI';
  359. 'FINSI';
  360. *----- Appel a la procedure thermique (lineaire ou nonlineaire) ----
  361. 'SI' ('NON' WTAB.'COUPLAGE'); 'MENAGE'; 'FINSI';
  362. 'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'NONLINEAIRE');
  363. CHTER = TRANSNON PRECED;
  364. 'FINSI';
  365. 'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'LINEAIRE');
  366. CHTER = TRANSLIN PRECED;
  367. 'FINSI';
  368. 'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'DUPONT');
  369. WTAB.'NOUV_PAS' = ('EXTR' WTAB.'TEM_CALC' 1) - TEMP0;
  370. CHTER = DUPONT2 PRECED;
  371. 'FINSI';
  372.  
  373. * Verification d'une erreur en sortie
  374. 'SI' CHTER.'ERREUR';
  375. 'ERREUR' ' ';
  376. 'FINSI';
  377.  
  378. * Modele Metallurgique (T.L. en juin 2018)
  379. 'SI' WTAB.'FOR_METALLU' ;
  380. estim.'PROPORTIONS_PHASES' = CHTER.'PROPORTIONS_PHASES' ;
  381. 'FINSI';
  382.  
  383. 'SI' WTAB.'PHASE';
  384. WTAB.'PROPPHAS' = CHTER.'PROPPHAS' ;
  385. estim.'PROPORTIONS_PHASE'= CHTER.'PROPPHAS' ;
  386. 'FINSI' ;
  387.  
  388. 'SI' WTAB.'FOR_THER' ;
  389. WTAB.'TET2'=CHTER.'INITIAL(2)' ;
  390. WTAB.'CLIT'=CHTER.'CLIT' ;
  391. 'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'NONLINEAIRE') ;
  392. 'SI' ('NEG' WTAB.'RIBLO_T' 'INCONNU');
  393. WTAB.'RIBLO_T' = CHTER.'RIBLO_T' ;
  394. WTAB.'LISEA_T' = CHTER.'LISEA_T' ;
  395. 'FINSI' ;
  396. 'FINSI' ;
  397. 'SI' WTAB.'RAYO' ;
  398. WTAB.'RAYONNEMENT' = CHTER.'RAYONNEMENT' ;
  399. 'FINSI' ;
  400. 'FINSI' ;
  401. 'SI' WTAB.'FOR_DIFF' ;
  402. WTAB.'CO2' =CHTER.'CO2';
  403. WTAB.'CLID'=CHTER.'CLID' ;
  404. 'FINSI' ;
  405. 'FINS';
  406.  
  407. 'SI' WTAB.'FOR_THER' ;
  408. estim.'TEMPERATURES' = WTAB.'TET2' ;
  409. estim.'REACTIONS_THERMIQUES'='REAC' WTAB.'CLIT' WTAB.'TET2';
  410. 'FINSI';
  411.  
  412. 'SI' WTAB.'FOR_DIFF' ;
  413. estim.'CONCENTRATIONS' = WTAB.'CO2' ;
  414. estim.'REACTIONS_DIFFUSIONS' ='REAC' WTAB.'CLID' WTAB.'CO2';
  415. 'FINSI' ;
  416. *** Reevaluation des donnees du calcul dependant de la temperature
  417. 'SI' WTAB.'PROCEDURE_REEV_THE' ;
  418. REEV_THE PRECED 1 ;
  419. 'FINSI' ;
  420.  
  421. *-----convergence mecanique thermique : reprise configuration ------
  422. *---------- du pas precedent pour refaire la mecanique -------------
  423. 'SI'( WTAB.'GRANDS_DEPLACEMENTS' 'OU' WTAB.'FEFP_FORMULATION');
  424. 'FORM' WTAB.'GE0_DEB';
  425. 'FINSI' ;
  426. 'FINSI' ;
  427. *################## fin du calcul d'un pas thermique #################
  428.  
  429.  
  430. *#######################################################################
  431. *#################### CALCUL D'UN PAS SI MECANIQUE #####################
  432. *#######################################################################
  433. 'SI' WTAB.'FOR_MECA' ;
  434. 'SAUT' 1 'LIGN';
  435. 'MESSAGE' ' Numero du pas : ' IPAS ' Indice d evolution : ' ICHG ' -> temps :' TI;
  436.  
  437. *--------------------- calcul de la temperature ---------------------
  438. * CB215821 : On met la thermique dans le Chargement pour 'TIRER' le chargement
  439. * au bon instant dans UNPAS pour la MECANIQUE
  440. * Il faudra faire de meme avec la MECANIQUE si la THERMIQUE en depend
  441. * 'SI' WTAB.'FOR_THER';
  442. * THTA1='TABLE';
  443. * THTA2='TABLE';
  444. * THTA1.0=TEMP0;
  445. * THTA1.1=TI;
  446. * THTA2.0= conti.'TEMPERATURES' 'ENLE' 'LX' ;
  447. * THTA2.1= estim.'TEMPERATURES' 'ENLE' 'LX' ;
  448. * CHTH='CHAR' 'T' THTA1 THTA2;
  449. * CHASANST= WTAB.'CHARGEMENT' 'ENLE' 'T';
  450. * WTAB.'CHARGEMENT'= CHASANST 'ET' CHTH;
  451. * 'FINSI';
  452. 'SI' WTAB.'CHAR_THE';
  453. WTAB.'TET2' ='TIRE' WTAB.'CHARGEMENT' 'T ' TI;
  454. 'FINSI';
  455.  
  456. *------- Caracteristiques initiales en cas de grands deplacements -----
  457. 'SI' WTAB.'ITCAR';
  458. 'SI' ('EGA' WTAB.'CARA' 'INCONNU');
  459. WTAB.'CARA' = WTAB.'MAT0';
  460. 'FINSI';
  461. 'FINSI';
  462.  
  463. ************************************************************************
  464. **************** APPEL A UNPAS ********************************
  465. ************************************************************************
  466. 'MENA' ;
  467. TT = UNPAS PRECED;
  468.  
  469. *---------- Mise a jour du MODELE du PAS precedent ---------------
  470. WTAB.'MO_TOT_PREC' = TT.'MO_TOT_PREC' ;
  471.  
  472. *---------- retour a la configuration initiale -------------------
  473. 'SI' (WTAB.'GRANDS_DEPLACEMENTS' 'OU' WTAB.'FEFP_FORMULATION');
  474. 'FORM' WTAB.'FOR0';
  475. 'FINSI';
  476.  
  477. *------- test si on a converge
  478. WTAB.'CONV' = TT.'CONV';
  479. 'SI' ('NON' TT.'CONV' 'ET' ('NON' WTAB.'CONVERGENCE_FORCEE'));
  480. * en dynamique et en poreux on garde le PRECED.FREA1 du dernier pas
  481. * converge (dans le cas ou on voudrait continuer le calcul)
  482. * 'SI' WTAB.'DYNAMIQUE' ; WTAB.'FREA1' = F1 ; 'FINSI';
  483. * 'SI' WTAB.'CONSOLIDATION' ;WTAB.'FREA1' = FF1 ; 'FINSI';
  484. 'ERREUR' 997 ;
  485. 'FINSI';
  486.  
  487. * Remplissage de la table ESTIMATION : champs en fin de pas
  488. DEP0= conti.'DEPLACEMENTS';
  489. * DEP = ENLEVER DEP0 'LX';
  490. estim.'DEPLACEMENTS'= TT.'ZU1' ;
  491. * 'DETRUIRE' DEP;
  492. estim.'DEFORMATIONS'= TT.'DEFF';
  493. estim.'CONTRAINTES' = TT.'SIGF' ;
  494. estim.'DFGRAD' = TT.'DFGRAD';
  495. 'SI' WTAB.'IPLAVI' ;
  496. estim.'VARIABLES_INTERNES' = TT.'VARF';
  497. estim.'DEFORMATIONS_INELASTIQUES' = TT.'DFPF';
  498. 'FINSI' ;
  499. 'SI' WTAB.'GRANDS_DEPLACEMENTS' ;
  500. estim.'FOR' = TT.'FOR' ;
  501. 'FINSI' ;
  502. 'SI' WTAB . 'DYNAMIQUE' ;
  503. estim.'VITESSES' = TT.'VITI' ;
  504. estim.'ACCELERATIONS' = TT.'ACCEI' ;
  505. 'FINSI';
  506. * 'SI' WTAB .'ITCAR';
  507. * estim.'CARA_ACTU' = TT.'CARA' ;
  508. * 'FINSI';
  509. *** Reevaluation des donnees du calcul dependant de la mecanique
  510. 'SI' WTAB.'PROCEDURE_REEV_MEC' ;
  511. REEV_MEC PRECED 1;
  512. 'FINSI' ;
  513. 'FINSI';
  514. *#######################################################################
  515. *################### FIN DU CALCUL D'UN PAS MECANIQUE ##################
  516. *#######################################################################
  517.  
  518.  
  519.  
  520. ***
  521. *** Test de la convergence mecanique thermique
  522. ***
  523. 'SI' WTAB.'CONVERGENCE_MEC_THE' ;
  524. * je cherche la convergence thermo-mecanique : je dois au moins
  525. * calculer deux fois mon probleme
  526. 'SI' (I_BOTH > 1) ;
  527. * Test sur toutes les composantes (notamment, si DIFFUSION)
  528. CHT1 = WTAB . 'THER_COURANT' ;
  529. DUABS2 = 'ABS' (WTAB.'THER_COURANTP' - CHT1) ;
  530. CRICOH1 = 0. ;
  531. NBINCO1 = 'DIME' WTAB . 'PRIM_TOT' ;
  532. 'REPE' BINCO1 NBINCO1 ;
  533. LINCOI1 = 'MOTS' ('EXTR' WTAB.'PRIM_TOT' &BINCO1) ;
  534. MDUABSI1 = 'MAXI' DUABS2 'AVEC' LINCOI1 ;
  535. NORMI1 = 'MAXI' 'ABS' CHT1 'AVEC' LINCOI1 ;
  536. 'SI' (NORMI1 '&lt;EG' (1.E-14 * MDUABSI1)) ;
  537. NORMI1 = 1. ;
  538. 'FINS' ;
  539. CRICOHI1 = MDUABSI1 / NORMI1 ;
  540. 'SI' (CRICOHI1 '>' CRICOH1) ;
  541. CRICOH1 = CRICOHI1 ;
  542. 'FINS' ;
  543. 'FIN' BINCO1 ;
  544.  
  545. WTAB.'CONV_THMEC' = CRICOH1 '<' (WTAB . 'CRITERE_COHERENCE') ;
  546. * Fin de boucle Si ne converge Pas
  547. 'SI'('NON' WTAB.'CONV_THMEC') ;
  548. 'SI' (I_BOTH >EG WTAB.'BOU_METH');
  549. 'MESS' 'PAS DE COHERENCE TROUVEE MECANIQUE THERMIQUE' ;
  550. 'ERREUR' 999;
  551. 'FINSI' ;
  552. 'FINSI' ;
  553. 'SINON' ;
  554. * je suis sur la premiere boucle de convergence thermomecanique
  555. WTAB.'CONV_THMEC' = FAUX ;
  556. 'FINSI' ;
  557. * THER_COURANTPrecedent dans WTAB pour test CONVERGENCE_MEC_THE
  558. WTAB . 'THER_COURANTP' = WTAB . 'THER_COURANT' ;
  559. 'SINON' ;
  560. * ici je ne cherche pas la convergence thermomecanique
  561. WTAB.'CONV_THMEC' = VRAI ;
  562. 'FINSI' ;
  563.  
  564. 'SI' (WTAB.'CONV_THMEC') ;
  565. 'QUITTER' BO_BOTH ;
  566. 'FINSI' ;
  567.  
  568. 'SI' ('NEG' ('TYPE' WTAB.'GE0_DEB') 'MOT') ;
  569. FORME WTAB.'GE0_DEB';
  570. 'FINSI' ;
  571.  
  572. 'FIN' BO_BOTH ;
  573.  
  574. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  575. *++++++++++++++++ fin de boucle interne de convergence mec_ther ++++++++
  576. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  577.  
  578. *----- Preparation du pas suivant : Thermique et Diffusion
  579. 'SI' WTAB.'FOR_THER';
  580. WTAB.'TET2' = estim.'TEMPERATURES' ;
  581. WTAB.'TET1'=WTAB.'TET2';
  582. 'SI' ('EGA' WTAB.'PROCEDURE_THERMIQUE' 'DUPONT');
  583. WTAB.'TET0' = WTAB.'TET1' ;
  584. WTAB.'VIEUXPAS' = WTAB.'NOUV_PAS' ;
  585. WTAB.'RELAXATION_DUPONT' = CHTER.'RELAXATION_DUPONT';
  586. 'SI' ('EXIS' CHTER 'INITIAL(1)') ;
  587. WTAB.'TET1' = CHTER.'INITIAL(1)';
  588. 'FINSI' ;
  589. 'FINSI';
  590. 'FINSI' ;
  591.  
  592. 'SI' WTAB.'FOR_DIFF';
  593. WTAB.'CO1' = WTAB.'CO2' ;
  594. WTAB.'CO2' = estim.'CONCENTRATIONS' ;
  595. 'FINSI' ;
  596.  
  597. *----- Preparation du pas suivant : Mecanique
  598. 'SI' WTAB.'FOR_MECA' ;
  599. 'SI' ('EXISTE' TT 'SECOND_MEMBRE');
  600. WTAB . 'SECOND_MEMBRE' = TT . 'SECOND_MEMBRE' ;
  601. 'FINSI' ;
  602. 'SI' ('EXISTE' TT 'SUCCES');
  603. WTAB . 'SUCCES' = TT . 'SUCCES' ;
  604. 'FINSI' ;
  605. 'SI' ('EXISTE' TT 'AUTOCOEF');
  606. WTAB . 'AUTOCOEF' = TT . 'AUTOCOEF' ;
  607. 'FINSI' ;
  608. 'SI' ('EXISTE' TT 'AUTOREDU');
  609. WTAB . 'AUTOREDU' = TT . 'AUTOREDU' ;
  610. 'FINSI' ;
  611. 'SI' ('EXISTE' TT 'AUTORED1');
  612. WTAB . 'AUTORED1' = TT . 'AUTORED1' ;
  613. 'FINSI' ;
  614. 'SI' ('EXISTE' TT 'RIBLO_M');
  615. WTAB . 'RIBLO_M' = TT . 'RIBLO_M' ;
  616. WTAB . 'LISEA_M' = TT . 'LISEA_M' ;
  617. 'FINSI' ;
  618. 'SI' ('EXISTE' TT 'INCREMENT');
  619. WTAB . 'INCREMENT' = TT . 'INCREMENT' ;
  620. 'FINSI' ;
  621. 'SI' ('EXISTE' TT 'INITEMPS');
  622. WTAB . 'INITEMPS' = TT . 'INITEMPS' ;
  623. 'FINSI' ;
  624. 'SI' ('EXISTE' TT 'DT');
  625. WTAB . 'DT' = TT . 'DT' ;
  626. 'FINSI' ;
  627. 'SI' ('EXISTE' TT 'NBPLAS');
  628. WTAB . 'NBPLAS' = TT . 'NBPLAS' ;
  629. 'FINSI' ;
  630. 'SI' ('EXISTE' TT 'NOMBRE_ITERATIONS');
  631. WTAB . 'NOMBRE_ITERATIONS' = TT . 'NOMBRE_ITERATIONS' ;
  632. 'FINSI' ;
  633. 'SI' ('EXISTE' TT 'FFROT');
  634. WTAB . 'FFROT' = TT . 'FFROT' ;
  635. 'FINSI' ;
  636. 'SI' ('EXISTE' TT 'RESIDU');
  637. WTAB . 'RESIDU' = TT . 'RESIDU' ;
  638. 'FINSI' ;
  639. 'SI' ('EXISTE' TT 'FNONL');
  640. WTAB . 'FNONL' = TT . 'FNONL' ;
  641. 'FINSI' ;
  642. 'SI' ('EXISTE' TT 'DTPREC');
  643. WTAB . 'DTPREC' = TT . 'DTPREC' ;
  644. 'FINSI' ;
  645. * SP : Modification de WTAB.'TET2' suite a appel a UNPAS
  646. * CB : SI (IPILOT) ==> TET2 est TIRE entre TEMPS0 et TI
  647. 'SI' ('EXISTE' TT 'TET2');
  648. WTAB . 'TET2' = TT . 'TET2' ;
  649. 'FINSI' ;
  650. 'SI' ('EXISTE' TT 'XDENO');
  651. WTAB . 'XDENO' = TT . 'XDENO' ;
  652. 'FINSI' ;
  653. 'SI' ('EXISTE' TT 'XDENOM');
  654. WTAB . 'XDENOM' = TT . 'XDENOM' ;
  655. 'FINSI' ;
  656. 'SI' ('EXISTE' TT 'LASTKTAN');
  657. WTAB . 'LASTKTAN' = TT . 'LASTKTAN' ;
  658. 'FINSI' ;
  659. 'SI' ('EXISTE' TT 'DFGRAD');
  660. WTAB . 'DFGRAD'= TT .'DFGRAD';
  661. 'FINSI' ;
  662. 'SI' ('EXIS' TT 'ETHER2');
  663. WTAB.'ETHER1'= TT.'ETHER2';
  664. WTAB.'ETHER2'= TT.'ETHER2';
  665. 'FINSI';
  666. 'SI' ('EXIS' TT 'MSRTHER2');
  667. WTAB.'MSRTHER1'= TT.'MSRTHER2';
  668. WTAB.'MSRTHER2'= TT.'MSRTHER2';
  669. 'FINSI';
  670. 'SI' ('EXIS' TT 'DEFOR2');
  671. WTAB.'DEFOR1'= TT.'DEFOR2';
  672. WTAB.'DEFOR2'= TT.'DEFOR2';
  673. 'FINSI';
  674. 'SI'('EXIS' TT 'ETAT2');
  675. WTAB.'ETAT1'= TT.'ETAT2';
  676. 'FINSI';
  677.  
  678. * fin des reactualisations.
  679. *--------------------------- deplacement ------------------------------
  680. TT.'DEPTOT' = estim.'DEPLACEMENTS' ;
  681. *-------------------- forces de reaction ------------------------------
  682. TT.'REAC'= 'REAC' TT.'DEPTOT' WTAB.'CLIM';
  683. estim.'REACTIONS' = TT.'REAC' ;
  684. *------------------ deformation d'origine thermique ------------------
  685. 'SI' WTAB.'CHAR_THE';
  686. WTAB.'ETHER2' = TT.'ETHER2' ;
  687. 'SI' WTAB.'POR1' ;
  688. WTAB.'MSRTHER2' = TT.'MSRTHER2';
  689. 'FINSI';
  690. 'FINSI';
  691. *------- en dynamique : vitesse, forces de reactions, operateur -------
  692.  
  693. 'SI' WTAB . 'DYNAMIQUE' ;
  694. 'SI' WTAB . 'LIAISON_PERSISTANTE' ;
  695. * corection des vitesses pour avoir des vitesses relatives nulles
  696. * aux points qui sont en contact
  697. TT.'VITI' = TT.'VITI' + TT.'VAUX' ;
  698. 'SINON' ;
  699. * forces d'acceleration et d'amortissement a la fin du pas *
  700. WTAB.'FREA1' = TT.'FOPL' - TT.'REAC';
  701. * en cas d'appuis unilateraux avec choc elastique essai de corriger
  702. * les vitesses fournies par le schema
  703.  
  704. TT.'VITI'= VITEUNIL TT.'ZRAIDV' WTAB.'MASSE' TT.'VITI'
  705. tt.'DEPT' DEP0 WTAB.'DT' WTAB.'SECOND_MEMBRE' WTAB;
  706. 'SI' ( 'EXIS' WTAB 'RATE_VITEUNIL');
  707. 'ERREUR' 'Probleme dans viteunil' ; 'FINSI';
  708. 'FINSI' ;
  709. 'SI' ( 'NEG' WTAB.'REAPREC' 'INCONNU');
  710. 'SI' ('EGA' ('TYPE' (EXTR TT.'REAC' 'MAIL')) 'MAILLAGE');
  711. reacdif='REDU' WTAB.'REAPREC' (EXTR TT.'REAC' 'MAIL');
  712. TT.'REAC'= TT.'REAC' - REACDIF;
  713. 'FINSI' ;
  714. WTAB.'REAPREC'=TT.'REAC';
  715. 'SINON';
  716. WTAB . 'REAPREC'= TT.'REAC';
  717. 'FINSI';
  718. estim.'VITESSES' = TT.'VITI' ;
  719. estim.'REACTIONS' = TT.'REAC' ;
  720. 'FINSI';
  721.  
  722. *---------- en consolidation : operateur ------------------------------
  723. 'SI' WTAB . 'CONSOLIDATION' ;
  724. WTAB.'OPERATEUR' = ROP;
  725. 'FINSI';
  726.  
  727. *-------------------------- configuration -----------------------------
  728.  
  729. 'SI' WTAB.'GRANDS_DEPLACEMENTS' ;
  730. WTAB.'RECALCUL'= VRAI ;
  731. WTAB.'GE0_DEB' = TT.'FOR' ;
  732. 'SI' WTAB.'ITCAR' ;
  733. WTAB.'CARA' = TT.'CARA' ;
  734. 'FINSI' ;
  735. 'FINSI';
  736. 'FORM' WTAB.'FOR0';
  737.  
  738. ***********************************************************************
  739. 'SI' ('ET' (WTAB.'FEFP_FORMULATION')
  740. (WTAB.'UPDATE_LAGRANGIAN')) ;
  741. * tc 'SI' ( 'NEG' WTAB.'FOR' 'INCONNU'); DETR GEOM1; 'FINSI';
  742. WTAB.'GE0_DEB' = TT.'FOR' ; WTAB.'RECALCUL'=VRAI;
  743. 'FINSI';
  744. ************************************************************************
  745.  
  746. *--------------------- destruction de l operateur ----------------------
  747.  
  748. 'SI' (WTAB . 'RECADET') ;
  749. * 'DETR' RH 'ELEM';
  750. 'OUBLIER' ROP;
  751. 'FINSI';
  752.  
  753. *----------------------------- temps initial --------------------------
  754. DT2 =( TI - TEMP0)/10000.;
  755. 'SI' (WTAB . 'AUTOMATIQUE' 'OU' ('NON' WTAB.CONV));
  756. TEMP0 = ( TI - TEMP0) * WTAB.'AUTOCOEF' + TEMP0;
  757. 'MESS' ' Temps effectivement calcule et non converge : ' TEMP0;
  758. WTAB.'TEMPS0' = TEMP0;
  759. 'SINON';
  760. TEMP0=TI;
  761. WTAB.'TEMPS0' = TEMP0;
  762. 'FINSI';
  763.  
  764. *----------------- Appel a une procedure personnelle -------------------
  765.  
  766. 'SI' (WTAB.'PROCEDURE_PERSO1');
  767. PRECED.'PERSO1_APPEL' = 2 ;
  768. PERSO1 PRECED;
  769. 'SI' ('EXIS' PRECED 'ARRET') ;
  770. WTAB . 'ARRET' = PRECED . 'ARRET' ;
  771. 'FINSI' ;
  772. TEMP0 = WTAB.'TEMPS0' ;
  773. 'FINSI';
  774. *----------------- On redefinit WTAB.'CHARGEMENT_I' ------------------
  775. 'SI' WTAB.'CONVERGENCE_MEC_THE' ;
  776. WTAB.'CHARGEMENT_I' = WTAB.'CHARGEMENT' ;
  777. 'FINSI';
  778. 'FINSI';
  779.  
  780. *++++++++++++++++++ fin preparation pas suivant mecanique +++++
  781. *################## Sauvegarde ###############################
  782. 'SI' ('NON' WTAB.'FOR_MECA') ;
  783. DT2 =( TI - TEMP0)/10000.;
  784. TEMP0=TI;
  785. WTAB.'TEMPS0' = TEMP0;
  786. TT = 'TABLE' ;
  787. 'FINSI';
  788. * pour ne pas sauvegarder la matrice de raideur si on doit la recalculer
  789. 'SI' ('EXIS' WTAB 'RRRR') ;
  790. 'SI' WTAB.'RECARI' ;
  791. WTAB ='ENLE' WTAB 'RRRR' ;
  792. PRECED.'WTABLE'= WTAB ;
  793. 'SI' ('EXIS' estim 'FOR') ;
  794. estim ='ENLE' estim 'FOR' ;
  795. PRECED.'ESTIMATION'=estim ;
  796. 'FINSI';
  797. 'FINSI';
  798. 'FINSI';
  799.  
  800. *-- remplissage de la table de continuation et des tables de resultats
  801. *-- s'il faut garder les resultats, sauvegarde incrementale et
  802. *-- fantomisation si economique
  803. ISAUV = PAS_RESU TEMP0 iauto PRECED ;
  804.  
  805. *******************************************
  806. ***** cas du pilotage automatique *******
  807. *******************************************
  808.  
  809. 'SI' WTAB . 'AUTOMATIQUE'; IAUTO=IAUTO+1;
  810. 'SI' ((WTAB . 'AUTOCOEF' 'EGA' 0.9999999D0) 'OU'
  811. (TEMP0 '>EG' TI)) ;
  812. TEMP0=TI;
  813. 'MESS' ' Objectif atteint ';
  814. 'MESS' ' ';
  815. WTAB.'AUTODEUX' = FAUX ;
  816. WTAB.AUTORED1 = 0 ;
  817. WTAB.AUTOREDU = 1. ;
  818. IAUTO=1;
  819. 'SINON';
  820. 'MESS' ' Objectif non atteint: On recommence';
  821. 'MESS' ' ';
  822. 'SI' ( IAUTO >EG WTAB.'AUTOPAS');
  823. 'ERREUR' 1000;
  824. * 'MESS' 'Pilotage non converge apres ' IAUTO ' pas.';
  825. 'FINSI';
  826. WTAB.'AUTODEUX'=VRAI;ICHG=ICHG - 1;
  827. 'FINSI';
  828. 'FINSI';
  829. *
  830. 'SI'(WTAB.'CHAR_THE' 'OU' WTAB.'FOR_THER');
  831. WTAB.'TET1'=WTAB.'TET2';
  832. 'SI' WTAB.'FOR_MECA';
  833. WTAB.'ETHER1' =WTAB.'ETHER2';
  834. 'FINSI';
  835. 'SI' WTAB.'POR1';
  836. WTAB.'MSRTHER1'=WTAB.'MSRTHER2';
  837. 'FINSI';
  838. 'FINSI' ;
  839. 'SI' WTAB.'CHAR_DEFI' ; WTAB.'DEFOR1'=WTAB.'DEFOR2';'FINSI';
  840. 'SI' WTAB.'FOR_DIFF' ; WTAB.'CO1' =WTAB.'CO2' ;'FINSI' ;
  841. ************************************************************************
  842. ****** cas de la convergence forcee
  843. ************************************************************************
  844. 'SI' WTAB.'PAS_AJUSTE';
  845. 'SI' ((TI < TIV) 'OU' ('NON' WTAB.'CONV')); ICHG=ICHG-1;'FINSI';
  846. 'FINSI';
  847. 'SI' WTAB.'ARRET';
  848. 'QUITTER' BEXTERN;
  849. 'FINSI';
  850. ************************************************************************
  851. *************** FIN DE LA BOUCLE SUR LES PAS DE TEMPS ******************
  852. ************************************************************************
  853. 'FIN' BEXTERN;
  854. ************************************************************************
  855.  
  856. 'SI' ('NEG' WTAB.'AMORTISSEMENT' 'INCONNU');
  857. WTAB.'AMORTISSEMENT' = TOTO ;
  858. 'FINSI';
  859.  
  860. 'SI' ('NEG' ('TYPE' WTAB.'GE0_DEB') 'MOT') ;
  861. 'OUBL' WTAB 'GE0_DEB' ;
  862. 'FINSI' ;
  863.  
  864. * on tue WTAB.'MAT1' s'il existe car il peut etre tres gourmand en place
  865. WTAB.'MAT1' = 'MOT' 'INCONNU';
  866.  
  867. 'MESS' ' ';
  868. 'MESS'
  869. '------------------- FIN DE LA PROCEDURE PASAPAS ---------------------';
  870. 'SAUT' 1 'LIGNE';
  871. *
  872. ************************************************************************
  873. *
  874. 'SI' ('EXIS' WTAB 'RRRR') ;
  875. * pour forcer le calcul de la rigidite si on rappele pasapas
  876. WTAB1= 'ENLEVER' WTAB 'RRRR';
  877. PRECED .'WTABLE'= WTAB1;
  878. 'FINSI';
  879.  
  880. * Typage de la table en sortie de PASAPAS ;
  881. PRECED.'SOUSTYPE' = 'MOT' 'PASAPAS' ;
  882.  
  883. 'FINPROC' PRECED ;
  884.  
  885.  
  886.  

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