Télécharger pas_init.procedur

Retour à la liste

Numérotation des lignes :

  1. * PAS_INIT PROCEDUR MB234859 24/06/10 21:15:02 11939
  2. 'DEBPROC' PAS_INIT PRECED*'TABLE' ;
  3. ************************************************************************
  4. * Procedure qui renseigne les champs dans la table CONTINUATION.
  5. * et qui creee les tables ou seront sauves les resultats.
  6. *
  7. * Trois cas sont a distinguer :
  8. * - en cas d'un premier appel a PASAPAS, CONTINUATION n'existe pas
  9. * et elle est creee et completee de facon a contenir les champs
  10. * a l'instant TEMPS0 = 0.
  11. * - en cas de reprise d'un calcul, la procedure PAS_DEFA s'est chargee
  12. * de verifier que cet instant est bien present dans la liste des temps
  13. * calcules. Les champs resultats pour les instants posterieurs a cet
  14. * instant de reprise ont ete supprimes des tables de sauvegardes.
  15. * TEMPS0 = PRECED.'TEMPS'. NSOR et il faut renseigner CONTINUATION de
  16. * facon a ce qu'elle contienne les champs a cet instant.
  17. * - en cas de poursuite d'un calcul, on repart du dernier pas calcule.
  18. * La table CONTINUATION est redefinie dans le cas ou le modele
  19. * aurait change (cas de la fabrication additive ou des elements ont pu
  20. * etre ajoutes entre deux appels de PASAPAS).
  21. *
  22. * Entree : Table PRECED transmises a PASAPAS
  23. * Sortie : Aucune mais la table de PASAPAS contient desormais
  24. * - des tables pour stocker les resultats de calculs
  25. * - la table CONTINUATION qui contient les champs a TEMPS0
  26. * - la table ESTIMATION qui contiendra les champs a TI
  27. ************************************************************************
  28. WTAB = PRECED.'WTABLE';
  29. TABSAUV = WTAB.'TABSAUV';
  30. inouv = 'NON' ('EXIS' PRECED 'CONTINUATION') ;
  31. 'SI' inouv ;
  32. PRECED.'CONTINUATION' = 'TABLE';
  33. 'FINSI';
  34. conti = PRECED.'CONTINUATION';
  35. ************************************************************************
  36. * Initialisations de l'instant de calcul
  37. ************************************************************************
  38. 'SI' ('EXIS' PRECED 'TEMPS') ;
  39. NSOR = ('DIME' (PRECED.'TEMPS')) - 1;
  40. 'SI' inouv ;
  41. TEMP0 = PRECED.'TEMPS'. NSOR;
  42. 'SINON';
  43. TEMP0 = conti.'TEMPS';
  44. 'FINSI';
  45. 'SINON';
  46. NSOR = 0;
  47. 'SI' WTAB.'FREQUENTIEL' ;
  48. TEMP0 = WTAB.'TEMPS0' ;
  49. 'SINON' ;
  50. TEMP0 = 0.D0 ;
  51. 'FINSI' ;
  52. PRECED.'TEMPS' = 'TABLE';
  53. PRECED.'TEMPS'. 0 = TEMP0 ;
  54. 'FINSI';
  55. WTAB.'TEMPS0' = TEMP0;
  56. WTAB.'NSOR' = NSOR;
  57. conti.'TEMPS' = WTAB.'TEMPS0';
  58. ************************************************************************
  59. * Initialisation indices 'MODELES' de la table de sortie
  60. ************************************************************************
  61. 'SI' WTAB.'CHAR_MODE' ;
  62.  
  63. * Indices MODELES de la table de sortie :
  64. 'SI' ('NON' ('EXIS' PRECED 'MODELES')) ;
  65. PRECED.'MODELES' = 'TABLE' ;
  66. PRECED.'MODELES' . 0 = WTAB.'MODELE' ;
  67. * PRECED.'MAILLAGES' = 'TABLE' ;
  68. * PRECED.'MAILLAGES' . 0 = 'EXTR' (PRECED.'MODELES' . 0) 'MAIL' ;
  69. 'FINSI' ;
  70. conti.'MODELES' = PRECED.'MODELES'. NSOR ;
  71. * conti.'MAILLAGES' = PRECED.'MAILLAGES' . NSOR ;
  72.  
  73. 'FINSI' ;
  74. ************************************************************************
  75. * Initialisations pour la THERMIQUE
  76. ************************************************************************
  77. 'SI' WTAB.'FOR_THER';
  78. *
  79. * - Initialisation de la temperature
  80. 'SI' ('EXISTE' PRECED 'TEMPERATURES');
  81. TH_COUR = PRECED.'TEMPERATURES'. NSOR;
  82. 'SI' (NSOR '>EG' 2) ;
  83. WTAB.'TET0' = PRECED.'TEMPERATURES'. (NSOR - 2);
  84. 'FINSI';
  85. 'SINON';
  86. 'SI' ('EXIS' WTAB.'MAT_T' 'TINI') ;
  87. MOTSTEMP = 'EXTR' WTAB.'MOD_T' 'DEPL' ;
  88. DDDD = VIDE 'MCHAML' ;
  89. 'REPE' BMO ('DIME' MOTSTEMP) ;
  90. MOTEMPI = 'EXTR' MOTSTEMP &BMO ;
  91. DDDD = DDDD ET ('EXCO' WTAB.'MAT_T' 'TINI' MOTEMPI) ;
  92. 'FIN' BMO ;
  93. 'SINO' ;
  94. DDDD = 'ZERO' WTAB.'MOD_T' 'TEMPERAT' ;
  95. 'FINS' ;
  96. TH_COUR = 'CHANGER' 'CHPO' WTAB.'MOD_T' DDDD;
  97. PRECED.'TEMPERATURES' = 'TABLE';
  98. PRECED.'TEMPERATURES'. 0 = TH_COUR;
  99. 'FINSI';
  100. CONTI.'TEMPERATURES' = 'CHAN' 'ATTRIBUT' TH_COUR 'NATURE' 'DIFFUS';
  101. *
  102. * - Initialisation des reactions
  103. 'SI' ('NON' ('EXIS' PRECED 'REACTIONS_THERMIQUES')) ;
  104. 'SI' WTAB.'CHAR_BLOT' ;
  105. BLOTHER = 'TIRE' WTAB.'BLOCAGES_THERMIQUES' 'BLOT' TEMP0 ;
  106. 'SINON' ;
  107. BLOTHER = WTAB.'BLOCAGES_THERMIQUES' ;
  108. 'FINSI' ;
  109. QTH_COUR = 'REAC' BLOTHER CONTI.'TEMPERATURES' ;
  110. PRECED.'REACTIONS_THERMIQUES' = 'TABLE' ;
  111. PRECED.'REACTIONS_THERMIQUES' . 0 = QTH_COUR ;
  112. 'FINS' ;
  113. conti.'REACTIONS_THERMIQUES' = PRECED.'REACTIONS_THERMIQUES'. NSOR ;
  114. *
  115. * Indice 'TEMPERATURES_APPORT' de la table de sortie :
  116. 'SI' WTAB.'CHAR_MODE' ;
  117. 'SI' ('NON' ('EXIS' PRECED 'TEMPERATURES_APPORT')) ;
  118. PRECED.'TEMPERATURES_APPORT' = TABLE ;
  119. NBT1 = 'DIME' PRECED.'TEMPERATURES' ;
  120. 'REPE' BTAPP1 NBT1 ;
  121. IBT1 = &BTAPP1 - 1 ;
  122. PRECED.'TEMPERATURES_APPORT' . IBT1 = PRECED.'TEMPERATURES' . IBT1 ;
  123. 'FIN' BTAPP1 ;
  124. 'FINSI' ;
  125. conti.'TEMPERATURES_APPORT' = PRECED.'TEMPERATURES_APPORT'. NSOR ;
  126. 'FINSI' ;
  127. *
  128. 'FINS';
  129. ************************************************************************
  130. * Initialisations pour la MECANIQUE
  131. ************************************************************************
  132. 'SI' WTAB.'FOR_MECA' ;
  133. MO_TOT = WTAB.'MO_TOT';
  134. *
  135. * - Initialisation du deplacement
  136. 'SI' ('NON' ('EXIS' PRECED 'DEPLACEMENTS'));
  137. DDDD = 'ZERO' MO_TOT 'DEPLACEM';
  138. DEP_COUR = 'CHANGER' 'CHPO' MO_TOT DDDD;
  139. 'DETRUIRE' DDDD;
  140. PRECED.'DEPLACEMENTS' = 'TABLE';
  141. PRECED.'DEPLACEMENTS'. 0 = DEP_COUR;
  142. 'FINSI';
  143. conti.'DEPLACEMENTS' = PRECED.'DEPLACEMENTS'. NSOR;
  144. *
  145. * - Initialisation des contraintes
  146. 'SI' ('NON' ('EXIS' PRECED 'CONTRAINTES'));
  147. SIG_COUR = 'ZERO' MO_TOT 'CONTRAIN';
  148. PRECED.'CONTRAINTES' = 'TABLE';
  149. PRECED.'CONTRAINTES'. 0 = SIG_COUR;
  150. 'FINSI';
  151. conti.'CONTRAINTES' = PRECED.'CONTRAINTES'. NSOR;
  152. *
  153. * - Initialisation des vitesses
  154. 'SI' (WTAB.'DYNAMIQUE' 'OU' ('EXISTE' WTAB 'MOD_LIA'));
  155. 'SI' ('NON' ('EXISTE' PRECED 'VITESSES')) ;
  156. DDDD = 'ZERO' MO_TOT 'DEPLACEMENTS';
  157. VIT_COUR = 'CHANGER' 'CHPO' MO_TOT DDDD;
  158. 'DETRUIRE' DDDD;
  159. PRECED.'VITESSES' = 'TABLE';
  160. PRECED.'VITESSES'. 0 = VIT_COUR;
  161. 'FINSI';
  162. conti.'VITESSES' = PRECED.'VITESSES'. NSOR;
  163. 'FINSI';
  164. *
  165. * - Initialisation des accelerations
  166. 'SI' WTAB.'DYNAMIQUE';
  167. 'SI' ('NON' ('EXISTE' PRECED 'ACCELERATIONS')) ;
  168. DDDD = 'ZERO' MO_TOT 'DEPLACEMENTS';
  169. ACC_COUR = 'CHANGER' 'CHPO' MO_TOT DDDD;
  170. 'DETRUIRE' DDDD;
  171. PRECED.'ACCELERATIONS' = 'TABLE';
  172. PRECED.'ACCELERATIONS'. 0 = ACC_COUR;
  173. 'FINSI';
  174. conti.'ACCELERATIONS' = PRECED.'ACCELERATIONS'. NSOR;
  175. 'FINSI';
  176. *
  177. * - Initialisation des reactions
  178. 'SI' ('NON' ('EXIS' PRECED 'REACTIONS'));
  179. 'SI' WTAB.'CHAR_BLOM' ;
  180. BLOMECA = 'TIRE' WTAB.'BLOCAGES_MECANIQUES' 'BLOM' TEMP0 ;
  181. 'SINON' ;
  182. BLOMECA = WTAB.'BLOCAGES_MECANIQUES' ;
  183. 'FINSI' ;
  184. REA_COUR = 'REAC' BLOMECA conti.'DEPLACEMENTS';
  185. PRECED.'REACTIONS' = 'TABLE' ;
  186. PRECED.'REACTIONS' . 0 = REA_COUR;
  187. 'FINSI';
  188. conti.'REACTIONS' = PRECED.'REACTIONS'. NSOR;
  189. *
  190. 'SI' WTAB.'IPLAVI';
  191. *
  192. * - Initialisation des variables internes
  193. 'SI' ('NON' ('EXIS' PRECED 'VARIABLES_INTERNES')) ;
  194. VAR_COUR = 'ZERO' MO_TOT 'VARINTER';
  195. PRECED.'VARIABLES_INTERNES' = 'TABLE';
  196. PRECED.'VARIABLES_INTERNES'. 0 = VAR_COUR;
  197. 'FINSI';
  198. conti.'VARIABLES_INTERNES' = PRECED.'VARIABLES_INTERNES'. NSOR ;
  199. *
  200. * - Initialisation des deformations inelastiques
  201. 'SI' ('EXISTE' PRECED 'DEFORMATIONS_INELASTIQUES');
  202. DEI_COUR = PRECED.'DEFORMATIONS_INELASTIQUES' . NSOR;
  203. 'SINON';
  204. 'SI' inouv ;
  205. DEI_COUR = 'ZERO' MO_TOT 'DEFINELA';
  206. 'SI' ('EXIS' TABSAUV 'DEFIN');
  207. 'SI' TABSAUV.'DEFIN';
  208. PRECED.'DEFORMATIONS_INELASTIQUES' = 'TABLE';
  209. PRECED.'DEFORMATIONS_INELASTIQUES'. 0 = DEI_COUR;
  210. 'FINSI';
  211. 'FINSI';
  212. 'SINON' ;
  213. DEI_COUR = conti.'DEFORMATIONS_INELASTIQUES' ;
  214. 'FINSI' ;
  215. 'FINSI';
  216. conti.'DEFORMATIONS_INELASTIQUES' = DEI_COUR;
  217. *
  218. 'FINSI';
  219. *
  220. * - Initialisation des deformations totales
  221. 'SI' ('EXISTE' PRECED 'DEFORMATIONS');
  222. DEF_COUR = PRECED.'DEFORMATIONS'. NSOR;
  223. 'SINON';
  224. 'SI' inouv ;
  225. DEF_COUR = 'ZERO' MO_TOT 'DEFORMATIONS';
  226. 'SI' ('EXIS' TABSAUV 'DEFTO');
  227. 'SI' TABSAUV.'DEFTO';
  228. PRECED.'DEFORMATIONS' = 'TABLE';
  229. PRECED.'DEFORMATIONS' . 0 = DEF_COUR;
  230. 'FINSI';
  231. 'FINSI';
  232. 'SINON' ;
  233. DEF_COUR = conti.'DEFORMATIONS' ;
  234. 'FINSI' ;
  235. 'FINSI';
  236. conti.'DEFORMATIONS' = DEF_COUR;
  237. *
  238. * - Initialisation des deformations apport
  239. 'SI' ('EXISTE' PRECED 'DEFORMATIONS_APPORT');
  240. DEF_APPO = PRECED.'DEFORMATIONS_APPORT'. NSOR ;
  241. 'SINON';
  242. 'SI' inouv;
  243. DEF_APPO = 'ZERO' MO_TOT 'DEFORMATIONS';
  244. 'SI' ('EXIS' TABSAUV 'DEFAP');
  245. 'SI' TABSAUV.'DEFAP';
  246. PRECED.'DEFORMATIONS_APPORT' = 'TABLE';
  247. PRECED.'DEFORMATIONS_APPORT' . 0 = DEF_APPO;
  248. 'FINSI';
  249. 'FINSI' ;
  250. 'SINON' ;
  251. DEF_APPO = conti.'DEFORMATIONS_APPORT' ;
  252. 'FINSI';
  253. 'FINSI';
  254. conti.'DEFORMATIONS_APPORT' = DEF_APPO;
  255.  
  256. * - Initialisation pour un tracking
  257. 'SI' ('EXIS' WTAB.'MODELE' 'MATE' 'EFEM');
  258. 'SI' ('EXIS' PRECED 'TRACKING');
  259. 'SI' PRECED.'TRACKING' ;
  260. 'SI' ('NON' ('EXIS' PRECED 'PROOTS'));
  261. PRECED.'PROOTS' = 'TABLE';
  262. PRECED.'PROOTS'. 0 = 0.;
  263. 'FINSI';
  264. * 'SI' inouv; conti.'PROOTS' = PRECED.'PROOTS'.NSOR;'FINS';
  265. 'SI' ('NON' ('EXIS' PRECED 'EROOTS'));
  266. PRECED.'EROOTS' = 'TABLE';
  267. PRECED.'EROOTS'. 0 = 0.;
  268. 'FINSI';
  269. * 'SI' inouv; conti.'PROOTS' = PRECED.'PROOTS'.NSOR;'FINS';
  270. 'SI' ('NON' ('EXIS' PRECED 'NROOTS'));
  271. PRECED.'NROOTS' = 0;
  272. 'FINSI';
  273. * 'SI' inouv; conti.'EROOTS' = PRECED.'EROOTS'.NSOR;'FINS';
  274. 'SI' ('NON' ('EXIS' PRECED 'TABISO'));
  275. PRECED.'TABISO' = 'TABLE';
  276. PRECED.'TABISO'. 0 = 0.;
  277. 'FINSI';
  278. * 'SI' inouv; conti.'TABISO' = PRECED.'TABISO'.NSOR;'FINS';
  279. 'SI' ('NON' ('EXIS' PRECED 'RESTH'));
  280. PRECED.'RESTH' = 'TABLE';
  281. PRECED.'RESTH'. 0 = 0.;
  282. 'FINSI';
  283. * 'SI' inouv; conti.'RESTH' = PRECED.'RESTH'.NSOR;'FINS';
  284. 'SI' ('NON' ('EXIS' PRECED 'ISOTEST'));
  285. PRECED.'ISOTEST' = 'TABLE';
  286. PRECED.'ISOTEST'. 0 = 0.;
  287. 'FINSI';
  288. * 'SI' inouv; conti.'ISOTEST' = PRECED.'ISOTEST'.NSOR;'FINS';
  289. 'SI' ('NON' ('EXIS' PRECED 'ISOTOT'));
  290. PRECED.'ISOTOT' = 'TABLE';
  291. PRECED.'ISOTOT'. 0 = 0.;
  292. 'FINSI';
  293. * 'SI' inouv; conti.'ISOTOT' = PRECED.'ISOTOT'.NSOR;'FINS';
  294. 'SI' ('NON' ('EXIS' PRECED 'MESH_ISO'));
  295. PRECED.'MESH_ISO' = 'TABLE';
  296. PRECED.'MESH_ISO'. 0 = 0.;
  297. 'FINSI';
  298. * 'SI' inouv; conti.'MESH_ISO' = PRECED.'MESH_ISO'.NSOR;'FINS';
  299. 'SI' ('NON' ('EXIS' PRECED 'MESH'));
  300. PRECED.'MESH' = 'TABLE';
  301. PRECED.'MESH'. 0 = 0.;
  302. 'FINSI';
  303. * 'SI' inouv; conti.'MESH' = PRECED.'MESH'.NSOR;'FINS';
  304. 'FINSI';
  305. 'FINSI';
  306. 'FINSI';
  307. *
  308. * - Initialisation pour le non-local HELM
  309. 'SI' WTAB.'NLOC' ;
  310. 'SI' ('EGA' WTAB.'NON_LOCAL' 'HELM');
  311. WTAB.'HELMHOLTZ' = 'TABLE' ;
  312. MOLVIA = 'EXTR' WTAB.'NLOC_MODL' 'NLOC' ;
  313. NHELM = 'DIME' MOLVIA ;
  314. WTAB.'HELMHOLTZ' . 'N_VARI_NL'= NHELM ;
  315. GEOHELM = 'EXTR' WTAB.'NLOC_MODL' 'MAIL' ;
  316. MAHELM = 'REDU' WTAB.'CARACTERISTIQUES' WTAB.'NLOC_MODL' ;
  317. * liste des parametres variables
  318. DDVA = 'EXTR' MAHELM 'COVA' ;
  319. WTAB.'HELMHOLTZ' . 'NLOC_COVA'= DDVA ;
  320. 'REPE' BH NHELM ;
  321. WTAB.'HELMHOLTZ'. &BH = 'TABLE' ;
  322. LEMO = 'EXTR' MOLVIA &BH ;
  323. chaih1=chai 'pas_init: Helmholtz pour :' LEMO;
  324. * mess chaih1;
  325. WTAB.'HELMHOLTZ'. &BH. 'NOM' = LEMO ;
  326. 'SI' ('EXIS' DDVA LEMO) ;
  327. WTAB.'HELMHOLTZ'. &BH. 'CTE' = FAUX ;
  328. 'SINON' ;
  329. WTAB.'HELMHOLTZ'. &BH. 'CTE' = VRAI ;
  330. 'FINSI' ;
  331. LEMOD = 'MODE' GEOHELM 'THERMIQUE' 'ANISOTROPE' ;
  332. WTAB.'HELMHOLTZ'. &BH . 'H_MODELE' = LEMOD ;
  333. * non linearite par defaut des helmholtz (sellier )
  334. WTAB.'HELMHOLTZ'. &BH . 'LINEAIRE' = FAUX ;
  335. 'FIN' BH ;
  336. 'FINSI';
  337. 'FINSI';
  338. *
  339. * En cas de chargement de deformations initialisation
  340. 'SI' WTAB.'CHAR_DEFI';
  341. WTAB.'DEFOR1' = 'TIRE' WTAB.'CHARGEMENT' 'DEFI' TEMP0 ;
  342. 'FINSI';
  343. *
  344. * En presence de THERMIQUE (CHARGEMENT ou FORMULATION) initialisation de TEMPERATURE_REFERENCE
  345. 'SI' WTAB.'ITHER' ;
  346. *
  347. * En cas de chargement thermique initialisation
  348. 'SI' WTAB.'CHAR_THE';
  349. WTAB.'TET1' = 'TIRE' WTAB.'CHARGEMENT' 'T' TEMP0;
  350. 'SINON';
  351. WTAB.'TET1' = CONTI.'TEMPERATURES' ;
  352. 'FINSI';
  353. *
  354. * On prend la temperature de reference dans le MATERIAU
  355. 'SI' ('EXIS' WTAB.'MA_TOT' 'TREF');
  356. TREF = 'ZERO' WTAB.'MO_TOT' 'TEMPERAT' ;
  357. TSZ = 'EXTR' WTAB.'MO_TOT' 'ZONE' ;
  358. ISZ = 1 ;
  359. 'REPE' BSZ (('DIME' TSZ) / 2) ;
  360. MOSZI = TSZ . ISZ ;
  361. MOTSTEMP = 'EXTR' MOSZI 'TEMP' ;
  362. 'REPE' BMO ('DIME' MOTSTEMP) ;
  363. MOTEMPI = 'EXTR' MOTSTEMP &BMO ;
  364. TREFI = 'CHAN' 'STRESSES' MOSZI ('EXCO' ('REDU' WTAB.'MA_TOT' MOSZI) 'TREF' MOTEMPI) ;
  365. TREF = TREF + TREFI ;
  366. 'FIN' BMO ;
  367. ISZ = ISZ + 2 ;
  368. 'FIN' BSZ ;
  369. WTAB.'TEMPERATURE_REFERENCE' = 'CHAN' 'TYPE' TREF 'TEMPERATURES' ;
  370. 'SINO';
  371. ERREUR 77 'AVEC' ('MOT' 'TREF') ('MOT' 'de MATERIAU');
  372. 'FINS';
  373. *
  374. *** Desormais initialise dans PASAPAS, au meme endroit que ETREF
  375. *** 'SI' ('EGA' WTAB.'ETHER1' 'INCONNU');
  376. *** On initialise les deformations thermiques a ZERO
  377. *** WTAB.'ETHER1' = 'ZERO' MO_TOT 'DEFORMATIONS' ;
  378. *** CB215821 : Pour une reprise propre, il faudrait mettre DEFORMATION_THERMIQUE dans CONTI
  379. *** 'FINS' ;
  380. 'FINS';
  381. 'FINS';
  382. ************************************************************************
  383. * Initialisation des differentes configurations (reference, debut pas)
  384. ************************************************************************
  385. WTAB.'FOR0' = 'FORM' ;
  386. WTAB.'GE0_DEB' = WTAB.'FOR0';
  387. 'SI' WTAB.'FOR_MECA' ;
  388. 'SI'(WTAB.'GRANDS_DEPLACEMENTS' 'OU' WTAB.'FEFP_FORMULATION');
  389. WTAB.'GE0_DEB' = 'FORM' conti.'DEPLACEMENTS';
  390. 'FINS';
  391. 'FORM' WTAB.'FOR0' ;
  392. 'FINS';
  393. ************************************************************************
  394. * Initialisations pour la MECANIQUE DES FLUIDES
  395. ************************************************************************
  396. 'SI' WTAB.'NVSTNL' ;
  397. 'SI' ('NON' ('EXISTE' PRECED 'VITESSES_FLUIDE')) ;
  398. DDDD = 'ZERO' WTAB.'MOD_NSL' 'DEPLACEMENTS' ;
  399. VFL_COUR = 'CHANGER' 'CHPO' WTAB.'MOD_NSL' DDDD ;
  400. 'DETRUIRE' DDDD ;
  401. PRECED.'VITESSES_FLUIDE' = 'TABLE' ;
  402. PRECED.'VITESSES_FLUIDE'. 0 = VFL_COUR ;
  403. 'FINSI' ;
  404. conti.'VITESSES_FLUIDE' = PRECED.'VITESSES_FLUIDE'. NSOR ;
  405. conti.'VITESSES_FLUIDE_0' = PRECED.'VITESSES_FLUIDE'. NSOR;
  406. 'FINSI' ;
  407. ************************************************************************
  408. * Initialisations pour la METALLURGIE (T.L. en juin 2018)
  409. ************************************************************************
  410. 'SI' WTAB.'FOR_METALLU' ;
  411. *
  412. * - Initialisation des proportions de phase
  413. 'SI' ('NON' ('EXIS' PRECED 'PROPORTIONS_PHASES')) ;
  414. * Valeurs par defaut : ZERO pour toutes les phases sauf la premiere : 1.D0
  415. CHMET0 ='ZERO' WTAB.'MOD_MET' 'VARINTER';
  416. LISCO ='EXTR' WTAB.'MOD_MET' 'VARI';
  417. CMP1 ='EXCO' CHMET0 ('EXTR' LISCO 1) ('EXTR' LISCO 1);
  418. CHMET0 = CHMET0 + (CMP1 + 1.D0);
  419. PRECED.'PROPORTIONS_PHASES' = 'TABLE' ;
  420. PRECED.'PROPORTIONS_PHASES'. 0 = CHMET0;
  421. 'FINS' ;
  422. conti.'PROPORTIONS_PHASES' = PRECED.'PROPORTIONS_PHASES'. NSOR ;
  423. CHTPOI0 = ('ZERO' WTAB.'MOD_MET' 'TEMPERATURES') 'NOMC' 'TPOI' ;
  424. CHTPOI0 = 'CHAN' 'CHPO' WTAB.'MOD_MET' CHTPOI0 ;
  425. WTAB . 'T_POINT0' = CHTPOI0 ;
  426. * 'SI' ('EXIS' PRECED 'T_POINT') ;
  427. * TAB_TPOI = 'INDE' (PRECED.'T_POINT') ;
  428. * CHTPOI0 = PRECED.'T_POINT'.NSOR ;
  429. * 'SINO' ;
  430. * CHTPOI0 = ('ZERO' WTAB.'MOD_MET' 'TEMPERATURES') 'NOMC' 'TPOI' ;
  431. * CHTPOI0 = 'CHAN' 'CHPO' WTAB.'MOD_MET' CHTPOI0 ;
  432. * PRECED.'T_POINT' = 'TABLE' ;
  433. * PRECED.'T_POINT' . 0 = CHTPOI0 ;
  434. * 'FINS' ;
  435. * conti . 'T_POINT' = CHTPOI0 ;
  436. 'FINSI' ;
  437. ************************************************************************
  438. * Initialisations pour la DIFFUSION
  439. ************************************************************************
  440. 'SI' WTAB.'FOR_DIFF' ;
  441. *
  442. * - Initialisation des concentrations
  443. 'SI' ('NON' ('EXIS' PRECED 'CONCENTRATIONS')) ;
  444. CHCO0 = 'ZERO' WTAB. 'MOD_DIF' 'DEPLACEMENTS' ;
  445. CHCO0 = 'CHAN' 'CHPO' WTAB. 'MOD_DIF' CHCO0 ;
  446. CHCO0 = 'CHAN' 'ATTRIBUT' CHCO0 'NATURE' 'DIFFUS';
  447. PRECED.'CONCENTRATIONS' = 'TABLE' ;
  448. PRECED.'CONCENTRATIONS'. 0 = CHCO0 ;
  449. 'FINS' ;
  450. CONTI.'CONCENTRATIONS' = PRECED.'CONCENTRATIONS'. NSOR ;
  451. *
  452. * - Initialisation des reactions
  453. 'SI' ('NON' ('EXIS' PRECED 'REACTIONS_DIFFUSIONS')) ;
  454. 'SI' WTAB.'CHAR_BLOD' ;
  455. BLODIFF = 'TIRE' WTAB.'BLOCAGES_DIFFUSIONS' 'BLOD' TEMP0 ;
  456. 'SINON' ;
  457. BLODIFF = WTAB.'BLOCAGES_DIFFUSIONS' ;
  458. 'FINSI' ;
  459. QDI_COUR = 'REAC' BLODIFF CONTI.'CONCENTRATIONS' ;
  460. PRECED.'REACTIONS_DIFFUSIONS' = 'TABLE' ;
  461. PRECED.'REACTIONS_DIFFUSIONS' . 0 = QDI_COUR ;
  462. 'FINS' ;
  463. CONTI.'REACTIONS_DIFFUSIONS' = PRECED.'REACTIONS_DIFFUSIONS'. NSOR ;
  464. 'FINS' ;
  465. ************************************************************************
  466. * Initialisations pour le CHANGEMENT DE PHASE (CB@!%*@! 2019)
  467. *-----------Doit etre place en dernier car peut necessiter les autres
  468. *-----------champs
  469. ************************************************************************
  470. 'SI' WTAB.'PHASE';
  471. *
  472. MO_COVA ='EXTR' WTAB.'MAT_PHA' 'COVA';
  473. 'SI' (('DIME' MO_COVA) > 0);
  474. CHAI1 ='CHAI' ' Le materiau du modele CHANGEMENT_PHASE ne peut pas varier';
  475. 'ERRE' CHAI1 ;
  476. 'FINS';
  477.  
  478. 'SI' ( 'EXISTE' PRECED 'PROPORTIONS_PHASE') ;
  479. CHPHA0 ='REDU' WTAB.'MOD_PHA' PRECED.'PROPORTIONS_PHASE'.NSOR ;
  480. 'SINON';
  481. * On met un etat qui a du sens physique pour la proportion de phases
  482. TABPHA ='EXTR' WTAB.'MOD_PHA' 'ZONE';
  483. CHPHA0 ='VIDE' 'MCHAML';
  484. 'REPE' SURi (('DIME' TABPHA) / 2);
  485. * Boucle sur les modeles CHANGEMENT_PHASE
  486. MODEi = TABPHA.((&SURi-1)*2 + 1) ;
  487. MATEi =('REDU' MODEi WTAB.'MAT_PHA' ) 'EXCO' 'PRIM';
  488. LMOPRIM='EXTR' MODEi 'DEPL';
  489. MOPRIM ='EXTR' LMOPRIM 1;
  490. 'SI' ('EXIS' WTAB.'CHARGEMENT' MOPRIM);
  491. * On regarde dans le CHARGEMENT si la PRIMALE existe
  492. CHPRIM ='TIRE' WTAB.'CHARGEMENT' MOPRIM TEMP0;
  493. 'SI' ('EXIS' ('EXTR' CHPRIM 'COMP') MOPRIM);
  494. CHPRIM='CHAN' 'CHAM' CHPRIM MODEi 'NOEUD';
  495. CHPRIM= CHPRIM 'NOMC' 'PRIM';
  496. CHPHAi=('MASQ' (MATEi-CHPRIM) 'EGINFE' 0.D0) 'NOMC' 'PPHA';
  497. CHPHA0= CHPHAi 'ET' CHPHA0;
  498. 'ITER' SURi ;
  499. 'FINS';
  500. 'FINS';
  501.  
  502. * On regarde dans les indices connus de la TABLE PRECED
  503. 'SI' ('EXIS' PRECED 'TEMPERATURES');
  504. CHPRIM = TH_COUR;
  505. 'SI' ('EXIS' ('EXTR' CHPRIM 'COMP') MOPRIM);
  506. CHPRIM='CHAN' 'CHAM' CHPRIM MODEi 'NOEUD';
  507. CHPRIM= CHPRIM 'NOMC' 'PRIM';
  508. CHPHAi=('MASQ' (MATEi-CHPRIM) 'EGINFE' 0.D0) 'NOMC' 'PPHA';
  509. CHPHA0= CHPHAi 'ET' CHPHA0;
  510. 'ITER' SURi ;
  511. 'FINS';
  512. 'FINS';
  513.  
  514. 'SI' ('EXIS' PRECED 'CONCENTRATIONS');
  515. CHPRIM = CHCO0;
  516. 'SI' ('EXIS' ('EXTR' CHPRIM 'COMP') MOPRIM);
  517. CHPRIM='CHAN' 'CHAM' CHPRIM MODEi 'NOEUD';
  518. CHPRIM= CHPRIM 'NOMC' 'PRIM';
  519. CHPHAi=('MASQ' (MATEi-CHPRIM) 'EGINFE' 0.D0) 'NOMC' 'PPHA';
  520. CHPHA0= CHPHAi 'ET' CHPHA0;
  521. 'ITER' SURi ;
  522. 'FINS';
  523. 'FINS';
  524.  
  525. * Si on arrive ici c'est qu'on n'a pas trouve le champ primal : on met zero
  526. CHPHAi='ZERO' MODEi 'VARINTER';
  527. CHPHA0= CHPHAi 'ET' CHPHA0;
  528. 'FIN' SURi;
  529.  
  530. PRECED.'PROPORTIONS_PHASE' = 'TABLE' ;
  531. PRECED.'PROPORTIONS_PHASE'. 0 = CHPHA0;
  532. 'FINSI';
  533. conti.'PROPORTIONS_PHASE' = CHPHA0;
  534. WTAB.'BLO_PHASE'='PMAT' WTAB.'MOD_PHA' ;
  535. 'FINSI';
  536. *
  537. 'SI' WTAB.'SOLUBILITE';
  538. WTAB.'BLO_SOL' = 'PMAT' WTAB.'MOD_SOL' ;
  539. 'FINSI';
  540. *
  541. PRECED.'ESTIMATION' = 'COPIER' conti;
  542. 'FINP';
  543.  
  544.  
  545.  
  546.  

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