Télécharger la plaquette Cast3M


fr en

Télécharger pas_defa.procedur

Retour à la liste

Numérotation des lignes : 

  1. * PAS_DEFA  PROCEDUR  MB234859  24/11/04    21:15:04     12064          
  2. DEBPROC PAS_DEFA PRECED*'TABLE';
  3. ************************************************************************
  4. * Procedure qui se charge de :
  5. * - verifier si d'anciens mots-cles ne sont pas utilises
  6. * - s'assurer de la presence de certaines donnees obligatoires
  7. * - verifier la coherence des donnees.
  8. * - initialiser les valeurs par defaut des options de calcul et les 
  9. *   surcharger par les valeurs fournies
  10. * - afficher un recapitulatif du calcul
  11. *
  12. * En cas de reprise, la configuration peut etre erronnee 
  13. * (reprise/poursuite d'un calcul suite a une erreur)
  14. ************************************************************************
  15. 'SI' (('EXISTE' PRECED 'GRANDES_ROTATIONS'   ) 'OU'
  16.       ('EXISTE' PRECED 'GRANDES_DEFORMATIONS'));
  17.   'MESS' ' ATTENTION :' ;
  18.   'MESS' ' Les donnees concernant les GRANDES ROTATIONS et GRANDES DEFORMATIONS ont ete modifiees :';
  19.   'MESS' ' Utilisez l''indice HYPOTHESE_DEFORMATIONS';
  20.   'ERREUR' 21 ;
  21. 'FINSI';
  22. *
  23. 'SI' ('EXISTE' PRECED 'CONTACT') ;
  24.   'MESS' ' ATTENTION :' ;
  25.   'MESS' ' Les donnees concernant le CONTACT dans PASAPAS ont ete modifiees :' ;
  26.   'MESS' ' - L''indice CONTACT de la table de PASAPAS n''existe plus.';
  27.   'MESS' ' - Il faut definir un MODELE de CONTACT.' ;
  28.   'ERREUR' 21 ;
  29. 'FINSI';
  30. *
  31. 'SI' (('EXISTE' PRECED 'PROCEDURE_REAC_MECA') 'OU'
  32.       ('EXISTE' PRECED 'PROCEDURE_CHAR_MECA') 'OU'
  33.       ('EXISTE' PRECED 'PROCEDURE_REAC_THER') 'OU'
  34.       ('EXISTE' PRECED 'PROCEDURE_CHAR_THER') 'OU'
  35.       ('EXISTE' PRECED 'PROCEDURE_PARA_THER')) ;
  36.   'MESS' ' ATTENTION :' ;
  37.   'MESS' ' Les donnees concernant les procedures utilisateur ont ete modifiees :' ;
  38.   'MESS' ' - L''indice PROCEDURE_REAC_MECA est remplace par l''indice PROCEDURE_REEV_MEC' ;
  39.   'MESS' ' - L''indice PROCEDURE_CHAR_MECA est remplace par l''indice PROCEDURE_CHARMECA' ;
  40.   'MESS' ' - L''indice PROCEDURE_REAC_THER est remplace par l''indice PROCEDURE_REEV_THE' ;
  41.   'MESS' ' - L''indice PROCEDURE_CHAR_THER est remplace par l''indice PROCEDURE_CHARTHER' ;
  42.   'MESS' ' - L''indice PROCEDURE_PARA_THER est remplace par l''indice PROCEDURE_PARATHER' ;
  43.   'ERREUR' 21 ;
  44. 'FINSI';
  45.  
  46. * Modifs CB215821 Juillet 2020
  47. 'SI' ('EXISTE' PRECED 'TEMPERATURE_REFERENCE');
  48.   'MESS' ' ATTENTION :' ;
  49.   'MESS' ' L''indice ''TEMPERATURE_REFERENCE'' de PASAPAS a ete supprime';
  50.   'MESS' ' Il faut desormais fournir la composante ''TREF'' dans le materiau';
  51.   'ERREUR' 21 ;
  52. 'FINSI';
  53. 'SI' ('EXISTE' PRECED 'TALPHA_REFERENCE');
  54.   'MESS' ' ATTENTION :' ;
  55.   'MESS' ' L''indice ''TALPHA_REFERENCE'' de PASAPAS a ete supprime';
  56.   'MESS' ' Il faut desormais fournir la composante ''TALP'' dans le materiau';
  57.   'ERREUR' 21 ;
  58. 'FINSI';
  59. ************************************************************************
  60. * Verifier la presence de certains indices de PRECED, obligatoires
  61. ************************************************************************
  62. 'SI' ('NON' ('EXISTE' PRECED 'MODELE')) ;
  63.   'ERREUR' 'Il manque la donnee du modele.';
  64. 'SINON' ;
  65.   MO    = PRECED.'MODELE';
  66.   TYPMO = 'TYPE' MO ;
  67.   'SI' (('NEG' TYPMO 'MMODEL  ') 'ET' ('NEG' TYPMO 'CHARGEME')) ;
  68.     'ERREUR' 'Le type du MODELE est incorrect';
  69.   'FINSI';
  70. 'FINSI' ;
  71. *
  72. 'SI' ('NON' ('EXISTE' PRECED 'CARACTERISTIQUES')) ;
  73.   'ERREUR' 'Il manque la donnee du materiau.';
  74. 'SINON' ;
  75.   MA    = PRECED.'CARACTERISTIQUES'  ;
  76.   TYPMA = 'TYPE' MA ;
  77.   'SI' (('NEG' TYPMA 'MCHAML  ') 'ET' ('NEG' TYPMA 'CHARGEME')) ;
  78.     'ERREUR' 'Le type du champ de materiau est incorrect';
  79.   'FINSI';
  80. 'FINSI' ;
  81. *
  82. 'SI' ('EXISTE' PRECED 'TEMPS_CALCULES');
  83.   TYP1 = 'TYPE' PRECED.'TEMPS_CALCULES';
  84.   'SI' ('NEG' TYP1 'LISTREEL');
  85.     'ERREUR' 'Le type de TEMPS_CALCULES est incorrect';
  86.   'FINSI';
  87. 'FINSI' ;
  88. *
  89. 'SI' ('EXISTE' PRECED 'CHARGEMENT');
  90.   TYP1 = 'TYPE' PRECED.'CHARGEMENT' ;
  91.   'SI' ('NEG' TYP1 'CHARGEME');
  92.     'ERREUR' 'Le type du CHARGEMENT est incorrect';
  93.   'FINSI';
  94. 'FINSI';
  95. *
  96. * Typage de la table en sortie de PASAPAS
  97. PRECED.'SOUSTYPE' = 'MOT' 'PASAPAS' ;
  98. ************************************************************************
  99. * Initialiser les indices de WTABLE a leur valeur par defaut
  100. ************************************************************************
  101. TWORK = 'TABLE' ;
  102. * adap_mode : adaptation du modele au cours du calcul
  103. TWORK.'ADAP_MODE'              = FAUX;
  104. * Adherence presence de la composante ADHE (modele contact ou cont/fro)
  105. TWORK.'ADHERENCE'              = FAUX ;
  106. * Advection presence d'un modele thermique advection ou diffusion advection
  107. TWORK.'ADVECTION'              = FAUX;
  108. * matrice d'amortissement en dynamique
  109. TWORK.'AMORTISSEMENT'          = 'MOT' 'INCONNU';
  110. * arret: indicateur pour quitter la boucle temporelle de PASAPAS
  111. TWORK.'ARRET'                  = FAUX;
  112. * autocrit: valeur du critere a utiliser dans autopilo(option automatique)
  113. TWORK.'AUTOCRIT'               = 0.001;
  114. * autocoef: coefficient de pilotage
  115. TWORK.'AUTOCOEF'               = 'MOT' 'INCONNU';
  116. * autodeux=vrai si pilotage et pas premiere tentative : pasapas
  117. TWORK.'AUTODEUX'               = FAUX;
  118. * automatique =vrai si demande de pilotage automatique
  119. TWORK.'AUTOMATIQUE'            = FAUX;
  120. * autopas : nb maxi de tentative en pilotage automatique
  121. TWORK.'AUTOPAS'                = 250;
  122. * autoreca : indique quand recalculer le coefficient de pilotage
  123. TWORK.'AUTORECA'               = 1;
  124. * beton_ht ( ne marche pas!)
  125. TWORK.'BETON_HT'               = FAUX;
  126. * blocages_diffusions; matrices
  127. TWORK.'BLOCAGES_DIFFUSIONS'    = 'VIDE' 'RIGIDITE'/'RIGIDITE';
  128. * blocages_mecaniques; matrices
  129. TWORK.'BLOCAGES_MECANIQUES'    = 'VIDE' 'RIGIDITE'/'RIGIDITE';
  130. * blocages_thermiques; matrices
  131. TWORK.'BLOCAGES_THERMIQUES'    = 'VIDE' 'RIGIDITE'/'RIGIDITE';
  132. * bou_meth : valeur par defaut du nombre de bouclage mecanique-thermique
  133. TWORK.'BOU_METH'               = 1;
  134. * cafrotte : veut dire il existe modele contact-frottement
  135. TWORK.'CAFROTTE'               = FAUX ;
  136. * capaconst logique il existe une capacite constante
  137. TWORK.'CAPACONST'              = FAUX;
  138. * caracteristiques :
  139. TWORK.'CARACTERISTIQUES'       = 'MOT' 'INCONNU';
  140. * ceramique presence de modele ceramique
  141. TWORK.'CERAMIQUE'              = FAUX;
  142. * indice celsius
  143. TWORK.'CELSIUS'                = FAUX;
  144. * chargement idem celui de preced ou creation d'un chargement farfelu
  145. TWORK.'CHARGEMENT'             = 'VIDE' 'CHARGEMENT' ;
  146. * chargement_i idem chargement  DFIS et FI si certain materiau pas_defa
  147. TWORK.'CHARGEMENT_I'           = 'MOT' 'INCONNU';
  148. * char_BLOD presence d'un chargement de BLOcages Diffusion
  149. TWORK.'CHAR_BLOD'              = FAUX;
  150. * char_blom presence d'un chargement de BLOcages Mecaniques
  151. TWORK.'CHAR_BLOM'              = FAUX;
  152. * char_blot presence d'un chargement de BLOcages Thermiques
  153. TWORK.'CHAR_BLOT'              = FAUX;
  154. * char_defi presence de chargement DEFI deformation initiale
  155. TWORK.'CHAR_DEFI'              = FAUX;
  156. * char_mate presence d'un chargement 'MATE'
  157. TWORK.'CHAR_MATE'              = FAUX;
  158. * char_mode : le modele est de type chargement
  159. TWORK.'CHAR_MODE'              = FAUX;
  160. * char_pres presence de chargement 'P   ' de pression suiveuse
  161. TWORK.'CHAR_PRES'              = FAUX;
  162. * char_the presence de chargement thermique pour la mecanique
  163. TWORK.'CHAR_THE'               = FAUX;
  164. * char_pilo presence de chargement de pilotage indirect
  165. TWORK.'CHAR_PILO'              = FAUX;
  166. TWORK.'DEPLACEMENTS_PILOTES'   = 'VIDE' 'CHPOINT' / 'DISCRET' ;
  167. TWORK.'FORCES_PILOTEES'        = 'VIDE' 'CHPOINT' / 'DISCRET' ;
  168. TWORK.'ZONE_DE_PILOTAGE'       = FAUX ;
  169. TWORK.'PARAMETRE_DE_PILOTAGE'  = 'MOT' 'INCONNU' ;
  170. * clim matrice des blocagesmecaniques y compris contact et frottement pasapas
  171. TWORK.'CLIM'                   = 'MOT' 'INCONNU';
  172. * condconst logique il existe une conductivite constante
  173. TWORK.'CONDCONST'              = FAUX;
  174. * conduction presence d'un modele de conduction
  175. TWORK.'CONDUCTION'             = FAUX;
  176. * consolidation presence de consolidation
  177. TWORK.'CONSOLIDATION'          = FAUX;
  178. * contact presence de contact
  179. TWORK.'CONTACT'                = FAUX ;
  180. * conv logique precisant si unpas a bien converge
  181. TWORK.'CONV'                   = VRAI;
  182. * convection presence d'un modele de convection
  183. TWORK.'CONVECTION'             = FAUX;
  184. * convergence_forcee logique demandant la convergence forcee
  185. TWORK.'CONVERGENCE_FORCEE'     = VRAI;
  186. * demande de convergence mecanique-thermique en cas de bouclage
  187. TWORK.'CONVERGENCE_MEC_THE'    = FAUX;
  188. * logique pour imposer au residu de decroitre de facon monotone
  189. TWORK.'CONVERGENCE_MONOTONE'   = FAUX;
  190. * conv_thmec on a convergee le bouclage mecanique-thermique  pasapas
  191. TWORK.'CONV_THMEC'             = VRAI;
  192. * couplage le calcul contient des modeles thermiques et des mecaniques
  193. TWORK.'COUPLAGE'               = FAUX;
  194. * cte_stefan_boltzmann valeur de la constante
  195. TWORK.'CTE_STEFAN_BOLTZMANN'   = 5.673E-8 ;
  196. * deja_calc_charther indique le dernier temps pour lequel on a calcule le
  197. * chargement thermique dans charther : charther
  198. TWORK.'DEJA_CALC_CHARTHER'     = 0.d0;
  199. * deltaiter sert au test de non convergence : unpas
  200. TWORK.'DELTAITER'              = 10;
  201. * dt: valeur du delta de temps entre le dernier pas et le pas en cours:pasapas
  202. TWORK.'DT'                     = 'MOT' 'INCONNU';
  203. * dtprec : dt du pas precedent: pasapas
  204. TWORK.'DTPREC'                 = 0 ;
  205. * dynamique calcul en dynamique
  206. TWORK.'DYNAMIQUE'              = FAUX;
  207. * endommagement il esiste un tel modele de comportement
  208. TWORK.'ENDOMMAGEMENT'          = FAUX;
  209. * Enregistrement de la date du depart de PASAPAS
  210. TWORK.'EPOQUE'                 = DATE 'EPOCH';
  211. * etat0 etat au temps 0 : pasapas pas_etat
  212. TWORK.'ETAT0'                  = 'MOT' 'INCONNU';
  213. * etat1 etat au temps courant : pasapas pas_etat
  214. TWORK.'ETAT1'                  = 'MOT' 'INCONNU';
  215. * etat_ref etat de reference : pasapas pas_etat
  216. TWORK.'ETAT_REF'               = 'MOT' 'INCONNU';
  217. * etref deformation d'origine thermique de reference : pas_epth
  218. TWORK.'ETREF'                  = 'MOT' 'INCONNU';
  219. * ether1 stocke des deformations thermiques au temps precedent : pasapas
  220. TWORK.'ETHER1'                 = 'MOT' 'INCONNU';
  221. * fnonl permet forces de reequilibrage du pas precedent : unpas
  222. TWORK.'FNONL'                  = 'MOT' 'INCONNU';
  223. * fopl : forces d'acceleration et d'amortissement 
  224. TWORK.'FOPL'                   = 'MOT' 'INCONNU';
  225. * configuration finale  : pasapas unpas
  226. TWORK.'FOR'                    = 'MOT' 'INCONNU';
  227. * configuration en debut de pas : pasapas
  228. TWORK.'FOR0'                   = 'MOT' 'INCONNU';
  229. * FOR_DIFF le calcul inclus des modeles DIFFUSION
  230. TWORK.'FOR_DIFF'               = FAUX;
  231. * FOR_METALLU le calcul inclus des modeles METALLURGIE
  232. TWORK.'FOR_METALLU'            = FAUX;
  233. * utilisation de la frormulation fefp
  234. TWORK.'FEFP_FORMULATION'       = FAUX;
  235. * frea1 partie constante du second membre  dynamique,consolidation: pasapas
  236. TWORK.'FREA1'                  = 'MOT' 'INCONNU';
  237. * frequentiel calcul en frequentiel
  238. TWORK.'FREQUENTIEL'            = 'EGA' ('VALE' 'MODE') 'FREQ' ;
  239. * ftol desequilibre tolere
  240. TWORK.'FTOL'                   = 'MOT' 'INCONNU';
  241. * Configuration courante
  242. TWORK.'GE0_DEB'                = 'MOT' 'INCONNU';
  243. * grands_deplacement
  244. TWORK.'GRANDS_DEPLACEMENTS'    = FAUX;
  245. * hypothese de calcul pour les deformations
  246. TWORK.'HYPOTHESE_DEFORMATIONS' = 'MOT' 'LINEAIRE';
  247. * indice d'evolution de la boucle temporelle de pasapas
  248. TWORK.'ICHG'                   = 0;
  249. * increment de charge du pas precedent : unpas
  250. TWORK.'INCREMENT'              ='MOT' 'INCONNU';
  251. * logique pour autoriser l'Initialisation a partir du pas precedent : unpas
  252. TWORK.'INITIALISATION'         = VRAI;
  253. * iplavi il existe un modele devant passer par COMP
  254. TWORK.'IPLAVI'                 = FAUX;
  255. * ipresu : listreel des temps a mettre dans les tables resultats
  256. TWORK.'IPRESU'                 = 'MOT' 'INCONNU';
  257. * iptime  liste des temps a calculer
  258. TWORK.'IPTIME'                 = 'MOT' 'INCONNU';
  259. * isauvega liste des temps sauvegardes (appel SAUV dans pasapas)
  260. TWORK.'ISAUVEGA'               = 'MOT' 'INCONNU';
  261. * logique la sauvegarde dans pasapas est activee
  262. TWORK.'ISAUVEGL'               = FAUX;
  263. * nombre de sous-pas realise dans unpas
  264. TWORK.'ISOUSPAS'               = 0 ;
  265. * itcar  des elements ont des caracteristiques geometriques
  266. TWORK.'ITCAR'                  = FAUX;
  267. * ither : pour un calcul meca presence de thermique (via mode ou char)
  268. TWORK.'ITHER'                  = FAUX;
  269. * k_sigma il faut utiliser l'operateur KSIGMA
  270. TWORK.'K_SIGMA'                = FAUX;
  271. * k_tangent il faut utiliser l'operateur KTANG
  272. TWORK.'K_TANGENT'              = FAUX;
  273. * k_tangent_iter0 matrice tangente a utiliser pour iteration 0 : unpas
  274. TWORK.'K_TANGENT_ITER0'        = 'MOT' 'INCONNU';
  275. * k_tangent_pert operateur KTAN matrice tangente par perturbation : unpas
  276. TWORK.'K_TANGENT_PERT'         = FAUX;
  277. * k_tangent_syme matrice tangente symetrisee : unpas
  278. TWORK.'K_TANGENT_SYME'         = FAUX;
  279. * k_tang_pert_c1 matrice tangente par perturbation coefficient C1 : unpas
  280. TWORK.'K_TANG_PERT_C1'         = 1.D-3;
  281. * k_tang_pert_c2 matrice tangente par perturbation coefficient C2 : unpas
  282. TWORK.'K_TANG_PERT_C2'         = 'MOT' 'INCONNU';
  283. * lastktan matrice donnant la derniere matrice tangente : pasapas unpas
  284. TWORK.'LASTKTAN'               = 'MOT' 'INCONNU';
  285. * langrangien schema numerique lagrangien reactualise par defaut
  286. TWORK.'LAG_TOT'                = 3 ;
  287. * ldeva liste des variables dont dependent les materiaux
  288. TWORK.'LDEVA'                  = 'MOT' 'INCONNU';
  289. * leconomi l'option economique est activee : pas_sauv
  290. TWORK.'LECONOMI'               = FAUX;
  291. * liaison_persistante si il existe des liaisons persistantes
  292. TWORK.'LIAISON_PERSISTANTE'    = FAUX;
  293. * linesearch en grand deplacement
  294. TWORK.'LINESEARCH'             = FAUX;
  295. * lisea_m  pour les contacts unilateraux en mecaniques :pasapas unpas
  296. TWORK.'LISEA_M'                = 'MOT' 'INCONNU';
  297. * lisea_t pour les changements de phase en thermique : pasapas transnon
  298. TWORK.'LISEA_T'                = 'MOT' 'INCONNU';
  299. * lpar0 liste des parametres externes a passer a COMP
  300. TWORK.'LPAR0'                  = 'MOT' 'INCONNU';
  301. * man methode asymptotique numerique
  302. TWORK.'MAN'                   =  FAUX;
  303. * massconst logique il existe des masses constantes
  304. TWORK.'MASSCONST'              = FAUX;
  305. * mat0 materiau au temps 0 : pasapas pas_epth
  306. TWORK.'MAT0'                   = 'MOT' 'INCONNU';
  307. * mat1  materiaux pour le pas en cours : pasapas unpas
  308. TWORK.'MAT1'                   = 'MOT' 'INCONNU';
  309. * matvar le materiau est variable  :  pas_verm
  310. TWORK.'MATVAR'                 = FAUX;
  311. * mat_con materiau des modeles de convection
  312. TWORK.'MAT_CON'                = 'MOT' 'INCONNU';
  313. * mat_dif materiau des modeles de diffusion
  314. TWORK.'MAT_DIF'                = 'MOT' 'INCONNU';
  315. * mat_mec materiau des modeles de mecanique
  316. TWORK.'MAT_MEC'                = 'MOT' 'INCONNU';
  317. * mat_pha materiau des modeles de changement de phase 'PARFAIT'
  318. TWORK.'MAT_PHA'                = 'MOT' 'INCONNU';
  319. * mat_por materiau des modeles poreux 
  320. TWORK.'MAT_POR'                = 'MOT' 'INCONNU';
  321. * mat_ray materiau des modeles de rayonnement 
  322. TWORK.'MAT_RAY'                = 'MOT' 'INCONNU';
  323. * mat1  materiaux de reference
  324. TWORK.'MAT_REF'                = 'MOT' 'INCONNU';
  325. * mat_sol materiau des modeles de changement de phase 'SOLUBILITE' 
  326. TWORK.'MAT_SOL'                = 'MOT' 'INCONNU';
  327. * mat_T materiau thermique + thermohydrique 
  328. TWORK.'MAT_T'                  = 'MOT' 'INCONNU';
  329. * mat_the materiau des modeles thermiques 
  330. TWORK.'MAT_THE'                = 'MOT' 'INCONNU';
  331. * mat_thm materiau thermohydrique 
  332. TWORK.'MAT_THM'                = 'MOT' 'INCONNU';
  333. * materiau des modeles convection + thermiques 
  334. TWORK.'MAT_TOT'                = 'MOT' 'INCONNU';
  335. * maxdefor valeur limitant la deformation en cas de non convergence
  336. * ou de RENORMALISATION souhaitee de ZDEP1 
  337. TWORK.'MAXDEFOR'               = 1.e-3;
  338. * maxiteration nombre maximum d'iteration pour un pas de calcul 
  339. TWORK.'MAXITERATION'           = 49;
  340. * maxisouspas nombre maximum de sous pas d'integration du fluage
  341. TWORK.'MAXISOUSPAS'            = 2000;
  342. * maxsousiterations nombre maximum de sous iterations pour un calcul non-local
  343. TWORK.'MAXSOUSITERATION'       = 500;  
  344. * maxsouspas nombre maximum de sous pas en cas de non convergence 
  345. TWORK.'MAXSOUSPAS'             = 100;
  346. * ma_tot materiau des modeles mecaniques + poreux +liaison 
  347. TWORK.'MA_TOT'                 = 'MOT' 'INCONNU';
  348. * mec1 il existe au moins un modele mecanique  
  349. TWORK.'MEC1'                   = FAUX;
  350. * affichage messages options calcul entete procedure resolution THERMIQUE
  351. TWORK.'MESSTHER'               = VRAI ;
  352. * modal ou statique 
  353. TWORK.'MODAL'                  = FAUX;
  354. * mod_adv modele de thermique, advection
  355. TWORK.'MOD_ADV'                = 'MOT' 'INCONNU';
  356. * mod_con modele de thermique convection  
  357. TWORK.'MOD_CON'                = 'MOT' 'INCONNU';
  358. * mod_dif modele de diffusion  
  359. TWORK.'MOD_DIF'                = 'MOT' 'INCONNU';
  360. * mod_mec modele mecanique
  361. TWORK.'MOD_MEC'                = 'MOT' 'INCONNU';
  362. * mod_pha modele de changement de phase 'PARFAIT' 
  363. TWORK.'MOD_PHA'                = 'MOT' 'INCONNU';
  364. * mod_por modele poreux
  365. TWORK.'MOD_POR'                = 'MOT' 'INCONNU';
  366. * mod_ray modele rayonnement
  367. TWORK.'MOD_RAY'                = 'MOT' 'INCONNU';
  368. * mod_sol modele de changement de phase 'SOLUBILITE' 
  369. TWORK.'MOD_SOL'                = 'MOT' 'INCONNU';
  370. * mod_soq : modele source de chaleur
  371. TWORK.'MOD_SOQ'                = 'MOT' 'INCONNU';
  372. * mod_T modele thermique + thermohydrique (contient inconnue T)
  373. TWORK.'MOD_T'                  = 'MOT' 'INCONNU';
  374. * mod_the modele thermique 
  375. TWORK.'MOD_THE'                = 'MOT' 'INCONNU';
  376. * mod_thm modele thermohydrique 
  377. TWORK.'MOD_THM'                = 'MOT' 'INCONNU';
  378. * mod_tot modele thermique 
  379. TWORK.'MOD_TOT'                = 'MOT' 'INCONNU';
  380. * mo_tot modele mecanique + modele poreux 
  381. TWORK.'MO_TOT'                 = 'MOT' 'INCONNU';
  382. * mo_tot_prec modele mecanique + modele poreux au pas precedent
  383. TWORK.'MO_TOT_PREC'            = 'MOT' 'INCONNU';
  384. * MODELE ce qui est dans preced.modele  
  385. TWORK.'MODELE'                 = 'MOT' 'INCONNU' ;
  386. * mova nom de la variables internes dont le max est imprime 
  387. TWORK.'MOVA'                   = 'MOT' 'EPSE';
  388. * msrther1 deformations thermiques pour les poreux
  389. TWORK.'MSRTHER1'               = 'MOT' 'INCONNU';
  390. * msrther2 deformations thermiques pour les poreux
  391. TWORK.'MSRTHER2'               = 'MOT' 'INCONNU';
  392. * mtol valeur autorisee du desequilibre (en terme de moment)
  393. TWORK.'MTOL'                   = 'MOT' 'INCONNU';
  394. * naures valeur de preced .autoresu   
  395. TWORK.'NAURES'                 = 100000;
  396. * premier pas de sauvegarde encore present (option economique de pasapas)
  397. TWORK.'NFSOR'                  = 0;
  398. * nombre maxi d'iteration internes si k_tangent 
  399. TWORK.'NITERINTER_MAX'         = 15;
  400. * niter_ktangent recalcule de la matrice tangente toutes les .. iterations
  401. TWORK.'NITER_KTANGENT'         = 1;
  402. * nloc s'il existe un modele de comportement non_local (preced.non_local)
  403. TWORK.'NLOC'                   = FAUX ;
  404. * nmaxsubsteps nb maximun de substepping si fefp  
  405. TWORK.'NMAXSUBSTEPS'           = 1;
  406. * nombre_iterations nb d'iterations du pas precedent : pasapas unpas
  407. TWORK.'NOMBRE_ITERATIONS'      = 12 ;
  408. * non_lineaire il existe un materiau utilisateur non-lineaire  
  409. TWORK.'NON_LINEAIRE'           = FAUX;
  410. * non_local decrit le modele de comportement non_local (preced.non_local)
  411. TWORK.'NON_LOCAL'              = 'MOT' 'INCONNU' ;
  412. * delta temps du pas en cours pour dupont2 : pasapas
  413. TWORK.'NOUV_PAS'               = 0.D0;
  414. *  nsor numero a ecrire dans les tables de resultats : pas_ init
  415. TWORK.'NSOR'                   = 0;
  416. * formulation NAVIER_STOKES materiau NLIN
  417. TWORK.'NVSTNL'                 = FAUX ; 
  418. * operateur matrice contenant l'operateur d'iterations : pasapas
  419. TWORK.'OPERATEUR'              = 'MOT' 'INCONNU';
  420. * Ordre est l'ordre de developpement de la MAN
  421. TWORK.'ORDRE'                  = 12;
  422. * para_exte il existe des parametres externes  : pas_verm
  423. TWORK.'PARA_EXTE'              = FAUX;
  424. * pas numero du pas de temps utile a transnon : pasapas
  425. TWORK.'PAS'                    = 0;
  426. * pas_ajuste si ajustement automatique de longueur de pas 
  427. TWORK.'PAS_AJUSTE'             = FAUX;
  428. * pas_redu  en automatiqu il a fallu reduire le pas  : pasapas unpas
  429. TWORK.'PAS_REDU'               = FAUX;
  430. * phase  il existe du changement de phase 'PARFAIT' 
  431. TWORK.'PHASE'                  = FAUX;
  432. * plastique il exite un comportement plastique 
  433. TWORK.'PLASTIQUE'              = FAUX;
  434. * por1 il existe un model poreux  
  435. TWORK.'POR1'                   = FAUX;
  436. * precflex critere de desequilibre sur moment autorise 
  437. TWORK.'PRECFLEX'               = 1.E-4;
  438. * precisinter precision pour les iterations locales (fefp)
  439. TWORK.'PRECISINTER'            = 1.D-8;
  440. * precision critere de convergence Thermique & Mecanique 
  441. TWORK.'PRECISION'              = 1.E-4 ;
  442. * precsousiteration precision pour les sous iterations en non local
  443. TWORK.'PRECSOUSITERATION'      = 1.0E-6;
  444. * predicteur
  445. TWORK.'PREDICTEUR'             = 'MOT' 'INCONNU';
  446. * procedure initialisation pour g_theta
  447. TWORK.'PROCEDURE'              = 'MOT' 'PASAPAS';
  448. * procedure_charmeca il faut appele la procedure charmeca  
  449. TWORK.'PROCEDURE_CHARMECA'     = FAUX;
  450. * procedure_charther il faut appele la procedure charther  
  451. TWORK.'PROCEDURE_CHARTHER'     = FAUX;
  452. * procedure_parather il faut appeler la procedure parather
  453. TWORK.'PROCEDURE_PARATHER'     = FAUX ;
  454. * procedure_perso1 il faut appeler la procedur perso1
  455. TWORK.'PROCEDURE_PERSO1'       = FAUX;
  456. * procedure_perso2 il faut appeler la procedur perso2
  457. TWORK.'PROCEDURE_PERSO2'       = FAUX;
  458. * procedure_reev_mec  il faut appeler la procedur reev_mec 
  459. TWORK.'PROCEDURE_REEV_MEC'     = FAUX ;
  460. * procedure_reev_the  il faut appeler la procedur reev_the 
  461. TWORK.'PROCEDURE_REEV_THE'     = FAUX ;
  462. * procedur_thermique non de la procedure thermique a appeler 
  463. TWORK.'PROCEDURE_THERMIQUE'    = 'MOT' 'NONLINEAIRE';
  464. * processeurs pour imposer l'option de parallelisme
  465. TWORK.'PROCESSEURS'            = 'MOT' 'INCONNU';
  466. * projection le champ de temperature clcule doit etere projete
  467. TWORK.'PROJECTION'             = FAUX;
  468. * raidconst il existe des raideurs constantes 
  469. TWORK.'RAIDCONST'              = FAUX;
  470. * raidaugm diffaugm theraugm il existe des raideurs augmentees 
  471. TWORK.'RAIDAUGM'               = FAUX;
  472. TWORK.'DIFFAUGM'               = FAUX;
  473. TWORK.'THERAUGM'               = FAUX;
  474. * autaug  augmentation automatique 
  475. TWORK.'AUTOAUGM'               = FAUX;
  476. * rayo il existe du rayonnement en thermique 
  477. TWORK.'RAYO'                   = FAUX;
  478. * rayonnement on y met des pre-calculs pour le rayonnement  : transnon charther
  479. TWORK.'RAYONNEMENT'            = 'MOT' 'INCONNU';
  480. * Reactualisation de la matrice en grands deplacements pour 0,1 de defor
  481. TWORK.'REAC_GRANDS'            = 0.1 ;
  482. * logique disant pour charther si geometrie a change : pasapas
  483. TWORK.'REA_GEOM'               = FAUX ;
  484. * reaprec utile pour les contacts unilateraux en dynamiques
  485. TWORK.'REAPREC'                = 'MOT' 'INCONNU';
  486. * recaop  = recalcul operateur
  487. TWORK.'RECAOP'                 = VRAI;
  488. * recari : recalcul de la rigidite
  489. TWORK.'RECARI'                 = FAUX;
  490. * reequilibrage : recalcul du dernier pas de temps lors d'une reprise
  491. TWORK.'REEQUILIBRAGE'          = FAUX ;
  492. * relaxation_dupont2 coeff de relaxation pour dupont2  
  493. TWORK.'RELAXATION_DUPONT'      = 0.25;
  494. * parametre de relaxation du pas de temps en non convergence
  495. TWORK.'RELAXATION_NONCONV'     = 1.;
  496. * relaxation_thet   coeff de relaxation pour transnon 
  497. TWORK.'RELAXATION_THETA'       = 1. ;
  498. * renormalisation de l'itere pour le limiter 
  499. TWORK.'RENORMALISATION'        = FAUX;
  500. * residu chpoint donnant le dernier residu d'equilibre : pasapas unpas
  501. TWORK.'RESIDU'                 = 'MOT' 'INCONNU';
  502. * riblo_m  matrice de blocages des appuis unilateraux : unpas
  503. TWORK.'RIBLO_M'                = 'MOT' 'INCONNU';
  504. * riblo_t matrice de blocages des changement de phase : transnon
  505. TWORK.'RIBLO_T'                = 'MOT' 'INCONNU';
  506. * valeurr des rigidites constantes
  507. TWORK.'RIGIDITE_CONSTANTE'     = 'MOT' 'INCONNU';
  508. * solubilite  il existe du changement de phase 'SOLUBILITE' 
  509. TWORK.'SOLUBILITE'             = FAUX;
  510. * sous_increment nombre de sous increment
  511. TWORK.'SOUS_INCREMENT'         = 1 ;
  512. TWORK.'RECALCUL_SOUSINC'       = VRAI ;
  513. * SOURCE_Q presence d'un modele thermique source de chaleur
  514. TWORK.'SOURCE_Q'               = FAUX;
  515. * test stabilite de la solution
  516. TWORK.'STABILITE'              = VRAI;
  517. * substepping  pour integrer en petites deformation 
  518. TWORK.'SUBSTEPPING'            = FAUX;
  519. * tabsauv table indiquant les options de sauv des resultats 
  520. TWORK.'TABSAUV'                = 'TABL' ;
  521. TWORK.'TABSAUV'.'DEFIN'        = FAUX;
  522. TWORK.'TABSAUV'.'DEFTO'        = FAUX;
  523. TWORK.'TABSAUV'.'DEFAP'        = FAUX;
  524. * temps_calcules liste des tempsd a calculer  
  525. TWORK.'TEMPS_CALCULES'         = 'MOT' 'INCONNU';
  526. * tem_calc liste des temps a calculer pour les proc thermique : pasapas
  527. TWORK.'TEM_CALC'               = 'MOT' 'INCONNU';
  528. * teta facteur pour la theta methode en consolidation  
  529. TWORK.'TETA'                   = 1.0D0;
  530. * tet0 champ de temperature 2  pas avant   : pasapas dupont2
  531. TWORK.'TET0'                   = 0;
  532. * tet1 champ de temperature 1  pas avant   : pasapas
  533. TWORK.'TET1'                   = 0;
  534. * tet2 champ de temperature du pas courant : pasapas
  535. TWORK.'TET2'                   = 0;
  536. * the1 le calcul inclus des modeles thermique 
  537. TWORK.'THE1'                   = FAUX;
  538. * ther_courant champ de temperature propose pour charther : dupont2 transnon
  539. TWORK.'THER_COURANT'           = 'MOT' 'INCONNU';
  540. * thm1 le calcul inclus des modeles thermohydriques 
  541. TWORK.'THM1'                   = FAUX;
  542. * tref temperature de reference 0 ou 273 : transnon dupont2
  543. TWORK.'TREF'                   = 0.D0;
  544. * ttol desequilibre absolu dans TRANSNON
  545. TWORK.'TTOL'                   = 'MOT' 'INCONNU';
  546. * update_lagrangian  on veut faire du lagrangien reactualise fefp
  547. TWORK.'UPDATE_LAGRANGIAN'      = FAUX;
  548. * presence d'un modele dit utilisateur
  549. TWORK.'UTILISATEUR'            = FAUX ;
  550. * vieuxpas dt precedent  pour dupont2
  551. TWORK.'VIEUXPAS'               = 0.D0;
  552. * viscodommage il existe un modele de comportement viscodommage 
  553. TWORK.'VISCODOMMAGE'           = FAUX;
  554. * viscoplastique il existe un modele de comportement viscoplastique  
  555. TWORK.'VISCOPLASTIQUE'         = FAUX;
  556. * visco_externe il existe un  modele de comportementvisco_externe  
  557. TWORK.'VISCO_EXTERNE'          = FAUX;
  558. * xdeno valeur pour calculer le critere de convergence force : unpas
  559. TWORK.'XDENO'                  = 0.d0;
  560. * xdenom valeur pour calculer le critere de convergence moment : unpas
  561. TWORK.'XDENOM'                 = 0.d0;
  562. * zraidv : rigidite K+4M/H2 (VITEUNIL)
  563. TWORK.'ZRAIDV'                 = 'MOT' 'INCONNU' ;
  564. ************************************************************************
  565. * Reprise/poursuite d'un calcul
  566. ************************************************************************
  567. * - Reprise d'un calcul : supprimer les infos des instants t > treprise
  568. 'SI' ('EXIS' PRECED 'REPRISE');
  569.   'SI' (PRECED.'REPRISE' > 0.);
  570.      PAS_REPR PRECED;
  571.   'FINSI';
  572. 'FINSI';
  573. *
  574. * - Poursuite d'un calcul : repartir du dernier instant calcule
  575. 'SI' ('EXIS' PRECED 'CONTINUATION');
  576. *
  577. * Si le dernier pas calcule n'a pas ete sauvegarde, l'ajouter aux tables 
  578. * de resultats en recuperant le contenu de CONTINUATION
  579.   NN = ('DIME' (PRECED.'TEMPS')) '-' 1 ;
  580.   'SI' ('NEG' PRECED.'CONTINUATION'.'TEMPS' PRECED.'TEMPS'. NN) ;
  581.     PRECED.'ESTIMATION' = 'TABLE' ;
  582.     ISAUV = PAS_RESU 0 PRECED ;
  583.   'FINSI' ;
  584. *
  585. * En cas de poursuite de calcul, indices eventuellement utiles a :
  586. * UNPAS, TRANSNON, TRANSLIN, DUPONT2 et autres procedures.
  587. * (a ameliorer...)
  588.   'SI' ('EXIS' PRECED 'WTABLE') ;
  589. *
  590. *   Pour initialisation a partir du pas precedent
  591.     'SI' ('EXIS' PRECED.'WTABLE' 'ZDEP1') ;
  592.       TWORK.'ZDEP1'     = PRECED.'WTABLE'.'ZDEP1' ;
  593.       TWORK.'FNONL'     = PRECED.'WTABLE'.'FNONL' ;
  594.       TWORK.'INCREMENT' = PRECED.'WTABLE'.'INCREMENT' ;
  595.       TWORK.'RESIDU'    = PRECED.'WTABLE'.'RESIDU' ;
  596.       TWORK.'ISOUSPAS'  = PRECED.'WTABLE'.'ISOUSPAS' ;
  597.     'FINSI' ;
  598. *
  599. *   Pour le pilotage automatique
  600.     TWORK.'AUTOCOEF' = PRECED.'WTABLE'.'AUTOCOEF';
  601.     TWORK.'XDENO'    = PRECED.'WTABLE'.'XDENO' ;
  602.     TWORK.'XDENOM'   = PRECED.'WTABLE'.'XDENOM' ;
  603. *
  604. *   En dynamique
  605.     TWORK.'FREA1'     = PRECED.'WTABLE'.'FREA1' ;
  606.     TWORK.'OPERATEUR' = PRECED.'WTABLE'.'OPERATEUR' ;
  607.     TWORK.'DTPREC'    = PRECED.'WTABLE'.'DTPREC' ;
  608. *
  609. *   Rigidite tangente
  610.     TWORK.'LASTKTAN' = PRECED.'WTABLE'.'LASTKTAN' ;
  611. *
  612. *   Pour PAS_AJUSTE
  613.     TWORK.'NOMBRE_ITERATIONS' = PRECED.'WTABLE'.'NOMBRE_ITERATIONS' ;
  614.  
  615. *   THERMIQUE ajout de matiere
  616.     'SI' ('EXIS' PRECED.'WTABLE' 'CAPA_COURANTE') ;
  617.       TWORK.'CAPA_COURANTE' = PRECED.'WTABLE'.'CAPA_COURANTE' ;
  618.     'FINSI' ;
  619.   'FINSI' ;
  620. *
  621. 'FINSI';
  622. *
  623. 'SI' (('EXIS' PRECED 'REPRISE') 'OU' ('EXIS' PRECED 'CONTINUATION')) ;
  624.   TWORK.'NSOR' = ('DIME' (PRECED.'TEMPS')) '-' 1 ;
  625.   'SI' ('EXISTE' PRECED 'REEQUILIBRAGE') ;
  626.     TWORK.'REEQUILIBRAGE' = PRECED.'REEQUILIBRAGE' ;
  627.   'FINSI' ;
  628. 'FINSI' ;
  629. TWORK.'RECAOP' = 'EGA' TWORK.'OPERATEUR' 'INCONNU' ;
  630. ************************************************************************
  631. * Stocker TWORK a l'indice WTABLE de PRECED
  632. ************************************************************************
  633. PRECED.'WTABLE' = TWORK ;
  634. ************************************************************************
  635. * Renseigner WTABLE a partir des informations de PRECED
  636. ************************************************************************
  637. 'SI' ('EGA' TYPMO 'CHARGEME') ;
  638.   TWORK.'CHAR_MODE' = VRAI ;
  639.   TWORK.'MODELE'    = 'TIRE' MO 'MODE' 0.D0 ;
  640. 'SINON' ;
  641.   TWORK.'MODELE' = MO ;
  642. 'FINSI' ;
  643. *-----------------------------------------------------------------------
  644. 'SI' ('EGA' TYPMA 'CHARGEME') ;
  645.   TWORK.'CHAR_MATE' = VRAI ;
  646.   TWORK.'CARACTERISTIQUES' = 'TIRE' MA 'MATE' 0.D0 ;
  647. 'SINON' ;
  648.   TWORK.'CARACTERISTIQUES' = MA ;
  649. 'FINSI';
  650. *-----------------------------------------------------------------------
  651. 'SI' ('EXIS' PRECED 'CHARGEMENT') ;
  652.   TWORK.'CHARGEMENT' = PRECED.'CHARGEMENT';
  653. 'FINSI';
  654. *-----------------------------------------------------------------------
  655. 'SI' ('EXISTE' PRECED 'PAS_AJUSTE');
  656.   TWORK.'PAS_AJUSTE' = PRECED.'PAS_AJUSTE';
  657.   'SI' TWORK.'PAS_AJUSTE' ;
  658.     TWORK.'CONVERGENCE_FORCEE' = FAUX ;
  659.   'FINSI';
  660. 'FINSI';
  661. 'SI' ('EXISTE' PRECED 'AUTOMATIQUE');
  662.   TWORK.'AUTOMATIQUE' = PRECED.'AUTOMATIQUE';
  663. * par defaut pas de controle de stabilite en pilotage
  664.   TWORK.'STABILITE'   = FAUX;
  665. 'FINSI' ;
  666. 'SI' TWORK.'AUTOMATIQUE'; 
  667.   TWORK.'PAS_AJUSTE' = FAUX; 
  668. 'FINSI';
  669. *
  670. 'SI' (('EXIS' PRECED 'TEMPS_CALCULES') 'ET' ('NON' TWORK.'PAS_AJUSTE'));
  671.   TWORK.'TEMPS_CALCULES'= PRECED.'TEMPS_CALCULES';
  672. 'SINON' ;
  673.   'SI' TWORK.'FREQUENTIEL' ;
  674.     l_zf = 'PROG' ; l_ii = 'PROG' 1. 'PAS' 0.3 4. 'PAS' 1. 10. ;
  675.     'SI' ('EXISTE' TWORK.'CARACTERISTIQUES' 'FREQ') ;
  676.       ccfr1 = 'EXCO' TWORK.'CARACTERISTIQUES' 'FREQ' 'NOID' ;
  677.       'SI' ('EGA' ('TYPE' ccfr1) 'MCHAML  ') ;
  678.         ncfr1 = 'EXTR' ccfr1 'NBZO' ;
  679.         'REPE' vfr1 ncfr1 ;
  680. ** hypothese 1 modele = 1 mode !!
  681.           f1 = 'EXTR' ccfr1 'FREQ' &vfr1 1 1 ;
  682.           log_f1 = ('LOG' f1) / ('LOG' 10.) ;
  683.           e10_f1 = 'ENTI' log_f1 ;
  684.           'SI' (log_f1 < 0.) ; e10_f1 = e10_f1 - 1 ; 'FINSI' ;
  685.           'SI'((log_f1 - e10_f1) 'EGA' 0.) ; e10_f1 = e10_f1 - 1 ; 'FINSI' ;
  686.           p10_f1 = 10. ** e10_f1 ;
  687.           l_zf1 = 'PROG' (f1 - (p10_f1/10.)) (f1 - (p10_f1/100.))
  688.                          (f1 - (p10_f1/1000.)) (f1 + (p10_f1/1000.)) ;
  689.           l_i10 = l_ii * p10_f1 ;
  690.           'SI' (('DIME' l_zf) 'EGA' 0) ;
  691.             e10_min = e10_f1 ; e10_max = e10_f1 + 1 ;
  692.             l_zf = 'ORDO' (l_i10 'ET' l_zf1) 'UNIQ' ;
  693.           'SINON' ;
  694.             l_zf0 = l_zf ;
  695.             'SI' ((e10_f1 '>EG' e10_min) 'ET' (e10_f1 < e10_max)) ;
  696.               l_zf = 'ORDO' (l_zf 'ET' l_zf1) 'UNIQ' ;
  697.             'SINON' ;
  698.               'SI' (e10_f1 < e10_min) ;
  699.                 e10_u = e10_f1 ;
  700.                 'REPE'  b_lz ;
  701.                   l_zf0 = l_i10 'ET' l_zf0 ; e10_u = e10_u + 1 ;
  702.                   'SI' (e10_u 'EGA' e10_min) ; 'QUIT' b_lz ; 'FINSI'  ;
  703.                   l_i10 = l_i10 * 10 ;
  704.                 'FIN'  b_lz ;
  705.                 e10_min = e10_f1 ;
  706.                 l_zf = 'ORDO' (l_zf0 'ET' l_zf1) 'UNIQ' ;
  707.               'FINSI'  ;
  708.               'SI' (e10_f1 '>EG' e10_max) ;
  709.                 e10_u = e10_f1 ;
  710.                 'REPE'  b_lz ;
  711.                   l_zf0 = l_i10 'ET' l_zf0 ; e10_u = e10_u - 1 ;
  712.                   'SI' (e10_u < e10_max) ; 'QUIT' b_lz ; 'FINSI'  ;
  713.                   l_i10 = l_i10 / 10 ;
  714.                 'FIN'  b_lz ;
  715.                 e10_max = e10_f1 + 1;
  716.                 l_zf = 'ORDO' (l_zf0 'ET' l_zf1) 'UNIQ' ;
  717.               'FINSI' ;
  718.             'FINSI' ;
  719.           'FINSI' ;
  720.         'FIN'  vfr1 ;
  721.       'SINON' ;
  722.         l_zf = l_ii ;
  723.       'FINSI' ;
  724.       TWORK.'TEMPS_CALCULES' = l_zf ;
  725.     'SINON' ;
  726.       TWORK.'TEMPS_CALCULES' = l_ii ;
  727.     'FINSI' ;
  728.   'SINON';
  729.     'SI' ('NON' ('EXISTE' PRECED 'TEMPS_SAUVES'));
  730.       'ERREUR' 'Il manque la donnee des temps a sauver ou a calculer';
  731.     'FINSI';
  732.     TWORK.'PAS_AJUSTE'     = VRAI;
  733.     TWORK.'TEMPS_CALCULES' = PRECED.'TEMPS_SAUVES';
  734.     TWORK.'TEMPS_SAUVES'   = PRECED.'TEMPS_SAUVES';
  735.   'FINSI';
  736. 'FINSI';
  737. TWORK.'IPTIME' = TWORK.'TEMPS_CALCULES';
  738. TWORK.'IPRESU' = TWORK.'TEMPS_CALCULES';
  739. 'SI' ('EXISTE' PRECED 'TEMPS_SAUVES' );
  740.   TWORK.'IPRESU' = PRECED.'TEMPS_SAUVES';
  741. 'FINSI';
  742. 'SI' ('EXISTE' PRECED 'TEMPS_SAUVEGARDES');
  743.    TWORK.'ISAUVEGL' = VRAI; 
  744.    TWORK.'ISAUVEGA' = PRECED.'TEMPS_SAUVEGARDES';
  745. 'FINSI';
  746. 'SI' ('EXISTE' PRECED 'ECONOMIQUE');
  747.     TWORK.'LECONOMI' = PRECED.'ECONOMIQUE';
  748. 'FINSI';
  749. 'SI' ('EXIS' PRECED 'MES_SAUVEGARDES');
  750.   TABSAUV = PRECED .'MES_SAUVEGARDES';
  751.   TWORK.'TABSAUV' = TABSAUV ;
  752.   'SI' ('NON' ( 'EXIS' TABSAUV 'DEFIN'));
  753.     TWORK.'TABSAUV'.'DEFIN' = FAUX;
  754.   'FINSI';
  755.   'SI' ('NON' ( 'EXIS' TABSAUV 'DEFTO'));
  756.     TWORK.'TABSAUV'.'DEFTO' = FAUX;
  757.   'FINSI';
  758.   'SI' ('NON' ( 'EXIS' TABSAUV 'DEFAP'));
  759.     TWORK.'TABSAUV'.'DEFAP' = FAUX;
  760.   'FINSI';
  761. 'FINSI';
  762. LOGERR1 = 'DANS' 'QUEL' TWORK.'IPRESU' TWORK.'IPTIME';
  763. LOGERR2 = VRAI;
  764. 'SI' TWORK.'ISAUVEGL';
  765.   LOGERR2= 'DANS' 'QUEL' TWORK.'ISAUVEGA' TWORK.'IPTIME';
  766. 'FINSI';
  767. 'SI' ('NON' (LOGERR1 'ET' LOGERR2));
  768.   'ERRE' 1013;
  769. 'FINSI';
  770. *
  771. *----- Indicateur ADAP_MODE d'adaptation du modele :
  772. 'SI' ('EXIS' PRECED 'TEMPS_ADAPTATION_MODELE') ;
  773.   TWORK.'ADAP_MODE' = VRAI ;
  774.   TWORK.'LTPS_ADAP' = PRECED.'TEMPS_ADAPTATION_MODELE' ;
  775. * TPS1_ADAP : prochain temps d'adaptation
  776.   TWORK.'TPS1_ADAP' = TWORK.'LTPS_ADAP' 'EXTR' 1 ;
  777.   TWORK.'INDI_ADAP' = 1 ;
  778. 'SINON' ;
  779.   TWORK.'TPS1_ADAP' = ('MAXI' TWORK.'TEMPS_CALCULES') + 1. ;
  780. 'FINSI' ;
  781. *-----------------------------------------------------------------------
  782. 'SI' ('EXISTE' PRECED 'CONSOLIDATION')  ;
  783.   TWORK.'CONSOLIDATION' = PRECED.'CONSOLIDATION'  ;
  784. 'FINSI';
  785. *
  786. 'SI' ('EXISTE' PRECED 'CONVERGENCE_MONOTONE');
  787.   TWORK.'CONVERGENCE_MONOTONE' = PRECED.'CONVERGENCE_MONOTONE';
  788. 'FINSI';
  789. *
  790. 'SI' ('EXISTE' PRECED 'DELTAITER');
  791.   TWORK.'DELTAITER' = PRECED.'DELTAITER';
  792. 'FINSI' ;
  793. *
  794. 'SI' ('EXISTE' PRECED 'DYNAMIQUE') ;
  795.   TWORK.'DYNAMIQUE' = PRECED.'DYNAMIQUE';
  796.   'SI' ('EXISTE' PRECED 'AMORTISSEMENT' );
  797.     TWORK.'AMORTISSEMENT' = 'EXTR' PRECED.'AMORTISSEMENT' 'RIGI' 'NOMU';
  798.   'FINSI';
  799. 'FINSI';
  800. *
  801. 'SI' ('EXISTE' PRECED 'FTOL' );
  802.   TWORK.'FTOL' = PRECED.'FTOL';
  803. 'FINSI';
  804. *
  805. 'SI' ('EXISTE' PRECED 'INITIALISATION');
  806.   TWORK.'INITIALISATION' = PRECED.'INITIALISATION';
  807. 'FINSI';
  808. *
  809. 'SI' ('EXISTE' PRECED 'LIAISON_PERSISTANTE' ) ;
  810.   TWORK.'LIAISON_PERSISTANTE' = PRECED.'LIAISON_PERSISTANTE';
  811. 'FINSI';
  812. *
  813. 'SI' ('EXISTE' PRECED 'MAXSOUSITERATION') ;
  814.   TWORK.'MAXSOUSITERATION' = PRECED.'MAXSOUSITERATION';
  815. 'FINSI'; 
  816. *
  817. 'SI' ('EXISTE' PRECED 'MAXSOUSPAS');
  818.   TWORK.'MAXSOUSPAS' = PRECED.'MAXSOUSPAS' ;
  819. 'FINSI';
  820. *
  821. 'SI' ('EXISTE' PRECED 'LINESEARCH');
  822.   'SI' ('EGA' PRECED.'LINESEARCH' VRAI);
  823.     TWORK.'LINESEARCH'          = VRAI;
  824.     TWORK.'GRANDS_DEPLACEMENTS' = VRAI;
  825.   'FINSI';
  826. 'FINSI';
  827. *
  828. 'SI' ('EXISTE' PRECED 'MTOL' );
  829.   TWORK.'MTOL' = PRECED.'MTOL';
  830. 'FINSI';
  831. *
  832. 'SI' ('EXISTE' PRECED 'PRECISION') ;
  833.   TWORK.'PRECISION' = PRECED.'PRECISION' ;
  834. 'FINS' ;
  835. *
  836. 'SI' ('EXISTE' PRECED 'PRECSOUSITERATION') ;
  837.   TWORK.'PRECSOUSITERATION' = PRECED.'PRECSOUSITERATION';
  838. 'FINSI'; 
  839. *
  840. 'SI' ('EXISTE' PRECED 'PREDICTEUR');
  841.   TWORK.'PREDICTEUR'          = 'MOT' PRECED.'PREDICTEUR';
  842.   TWORK.'GRANDS_DEPLACEMENTS' = VRAI;
  843. 'FINSI';
  844. *
  845. 'SI' ('EXISTE' PRECED 'PROCEDURE_CHARMECA');
  846.   TWORK.'PROCEDURE_CHARMECA' = PRECED.'PROCEDURE_CHARMECA';
  847. 'FINSI';
  848. *
  849. 'SI' ('EXISTE' PRECED 'PROCEDURE_CHARTHER');
  850.   TWORK.'PROCEDURE_CHARTHER' = PRECED.'PROCEDURE_CHARTHER';
  851. 'FINSI';
  852. *
  853. 'SI' ('EXISTE' PRECED 'PROCEDURE_PARATHER') ;
  854.   TWORK.'PROCEDURE_PARATHER' = PRECED.'PROCEDURE_PARATHER' ;
  855. 'FINSI';
  856. *
  857. 'SI' ('EXISTE' PRECED 'PROCEDURE_PERSO1');
  858.   TWORK.'PROCEDURE_PERSO1' = PRECED.'PROCEDURE_PERSO1';
  859. 'FINSI';
  860. *
  861. 'SI' ('EXISTE' PRECED 'PROCEDURE_PERSO2');
  862.   TWORK.'PROCEDURE_PERSO2' = PRECED.'PROCEDURE_PERSO2';
  863. 'FINSI';
  864. *
  865. 'SI' ('EXISTE' PRECED 'PROCEDURE_REEV_MEC');
  866.   TWORK.'PROCEDURE_REEV_MEC' = PRECED.'PROCEDURE_REEV_MEC' ;
  867. 'FINSI';
  868. *
  869. 'SI' ('EXISTE' PRECED 'PROCEDURE_REEV_THE');
  870.   TWORK.'PROCEDURE_REEV_THE' = PRECED.'PROCEDURE_REEV_THE' ;
  871. 'FINSI';
  872. *
  873. 'SI' ('EXISTE' PRECED 'PROCEDURE_THERMIQUE');
  874.   TWORK.'PROCEDURE_THERMIQUE' = PRECED.'PROCEDURE_THERMIQUE';
  875. 'FINSI';
  876. *
  877. 'SI' ('EXIS' PRECED 'REAC_GRANDS');
  878.   'SI' ('EGA' ('TYPE' PRECED . 'REAC_GRANDS') 'FLOTTANT') ;
  879.     TWORK.'REAC_GRANDS' = PRECED.'REAC_GRANDS';
  880.   'FINSI';
  881. 'FINSI';
  882. *
  883. 'SI' ('EXISTE' PRECED 'RENORMALISATION');
  884.   TWORK.'RENORMALISATION' = PRECED.'RENORMALISATION';
  885. 'FINSI';
  886. *
  887. 'SI' ('EXIS' PRECED 'SOUS-INCREMENT' );
  888.   TWORK.'SOUS-INCREMENT'=PRECED.'SOUS-INCREMENT';
  889. 'FINSI';
  890. *
  891. 'SI' ('EXISTE' PRECED 'STABILITE');
  892.   'SI' ('EGA' PRECED.'STABILITE' FAUX);
  893.     TWORK.'STABILITE'= FAUX;
  894.   'SINON';
  895.     TWORK.'STABILITE'= VRAI;
  896.   'FINSI';
  897. 'FINSI';
  898. *
  899. 'SI' ('EXISTE' PRECED 'CONVERGENCE_FORCEE');
  900.   TWORK.'CONVERGENCE_FORCEE' = PRECED.'CONVERGENCE_FORCEE';
  901. 'SINON';
  902.   'SI' ('NON' TWORK.'AUTOMATIQUE');
  903.     TWORK.'CONVERGENCE_FORCEE' = VRAI;
  904.   'FINSI';
  905. 'FINSI';
  906. *
  907. 'SI' ('EXISTE' PRECED 'CONVERGENCE_MEC_THE');
  908.   TWORK.'CONVERGENCE_MEC_THE' = PRECED.'CONVERGENCE_MEC_THE';
  909. 'FINSI';
  910. 'SI' TWORK.'CONVERGENCE_MEC_THE';
  911.   TWORK.'AUTOMATIQUE'        = FAUX;
  912.   TWORK.'CONVERGENCE_FORCEE' = FAUX;
  913.   'SI' TWORK.'PAS_AJUSTE' ;
  914.     'ERREUR' 'Options CONVERGENCE_MEC_THE et PAS_AJUSTE incompatibles';
  915.   'FINSI';
  916.   TWORK.'CRITERE_COHERENCE' = TWORK.'PRECISION' ;
  917.   'SI' ('EXISTE' PRECED 'CRITERE_COHERENCE') ;
  918.     TWORK.'CRITERE_COHERENCE' = PRECED.'CRITERE_COHERENCE' ;
  919.   'FINSI' ;
  920.   TWORK.'BOU_METH' = 10;
  921.   'SI' ('EXISTE' PRECED 'NB_BOTH');
  922.     TWORK.'BOU_METH'=PRECED.'NB_BOTH';
  923.   'FINSI';
  924. 'FINSI';
  925. *
  926. 'SI' ('EXISTE' PRECED 'MAXITERATION') ;
  927.   TWORK.'MAXITERATION' = PRECED.'MAXITERATION';
  928. 'SINON';
  929.   'SI' ((TWORK.'AUTOMATIQUE') 'OU' (TWORK.'PAS_AJUSTE')) ;
  930.     TWORK.'MAXITERATION' = 250;
  931.   'FINSI';
  932. 'FINSI';
  933. *
  934. 'SI' (TWORK.'CONVERGENCE_FORCEE' 'OU' TWORK.'RENORMALISATION');
  935.   TWORK.'MAXITERATION' = 250;
  936.   'SI' ('EXISTE' PRECED 'MAXDEFOR');
  937.     TWORK.'MAXDEFOR' = PRECED.'MAXDEFOR' ;
  938.   'FINSI';
  939.   'SI' TWORK.'CONVERGENCE_FORCEE';
  940.     'SI' ('EXISTE' PRECED 'RELAXATION_NONCONV');
  941.       TWORK.'RELAXATION_NONCONV' = PRECED.'RELAXATION_NONCONV' ;
  942.     'FINSI';
  943.   'FINSI';
  944. 'FINSI';
  945. *-----------------------------------------------------------------------
  946. * Note : Il faudrait separer la thermique de la mecanique pour la //.
  947. NBPART = 'VALEUR' 'ASSI' ;
  948. TWORK.'NBPART' = NBPART ;
  949. MPARA = 'MOT' TWORK.'PROCESSEURS' ;
  950. 'SI' ('EXISTE' PRECED 'PROCESSEURS');
  951.   MPARA = 'MOT' PRECED.'PROCESSEURS';
  952. 'FINSI';
  953. 'SI' ('&lt;EG' NBPART 1);
  954.   MPARA = 'MOT' 'MONO_PROCESSEUR' ;
  955. 'SINON' ;
  956.   'SI' ('EGA' MPARA 'INCONNU') ;
  957. *   anciennes valeurs: 2000 200
  958.     NBELMIN = 200 ; NBELMM = 400;
  959.     toutmo = 'EXTR' TWORK.'MODELE' 'ZONE' ;
  960.     IMOD = ('DIME' toutmo) / 2 ;
  961.     MINELE = 1000000 ; MAXELE = 0 ;
  962.     'REPE' BOUCMO IMOD ;
  963.       BBB = 'NBEL' TOUTMO . (&BOUCMO * 2);
  964.       'SI' (MINELE > BBB); MINELE = BBB; 'FINSI';
  965.       'SI' (MAXELE < BBB); MAXELE = BBB; 'FINSI';
  966.     'FIN' BOUCMO ;
  967.     NBTMP = (MINELE / NBELMM) ;
  968. **  mess 'nbtmp' ' ' nbtmp;
  969.     SI (NBTMP < NBPART); NBPART = NBTMP; 'FINSI';
  970.     TWORK.'NBPART' = NBPART ;
  971.     'SI' ('>' NBPART 1);
  972.       MINELE=MINELE/NBPART ; MAXELE=MAXELE/NBPART;
  973. **    mess 'minele maxele' ' ' minele ' ' maxele;
  974.       'SI' (MINELE > NBELMIN) ;
  975.         MPARA = 'MOT' 'AUTOMATIQUE' ;
  976.       'SINON';
  977.         'SI'( (MINELE > NBELMM) 'OU' (MAXELE > NBELMM)) ;
  978.           MPARA = 'MOT' 'COMPORTEMENT' ;
  979.         'FINSI';
  980.       'FINSI';
  981.     'SINON';
  982.       MPARA = 'MOT' 'MONO_PROCESSEUR' ;
  983.     'FINSI';
  984.   'FINSI';
  985. 'FINSI';
  986. LMPARA = 'MOTS' 'MONO_PROCESSEUR' 'AUTOMATIQUE' 'COMPORTEMENT' ;
  987. 'SI' ('EXIS' LMPARA MPARA) ;
  988.   TWORK.'PROCESSEURS' = 'MOT' MPARA ;
  989. 'SINON' ;
  990.   'MESS' ' *** Cas // non prevu' ; 'ERREUR' 5 ;
  991. 'FINSI' ;
  992. ************************************************************************
  993. * Decortiquer ce que contient le MODELE dans PAS_MODL
  994. ************************************************************************
  995. * Oublier l'indice MODEL_COURANT pour re-initialiser la table
  996. si (exis TWORK 'MODELE_COURANT') ;
  997.   oubl TWORK 'MODELE_COURANT' ;
  998. fins ;
  999. *
  1000. * Oublier l'indice 'MO_TOT_PREC' pour le reinitialiser dans UNPAS
  1001. si (exis TWORK 'MO_TOT_PREC') ;
  1002.   oubl TWORK 'MO_TOT_PREC' ;
  1003. fins ;
  1004.  
  1005. PAS_MODL PRECED ;
  1006. *
  1007. 'SI' ('NON' ('OU' TWORK.'FOR_MECA' TWORK.'FOR_THER' TWORK.'FOR_DIFF'
  1008.                   TWORK.'NVSTNL')) ;
  1009.   'ERREUR' 'Le modele ne correspond a aucune formulation valable';
  1010. 'FINSI';
  1011. *
  1012. 'SI' TWORK.'ADAP_MODE' ;
  1013.   'SI' ('OU' TWORK.'FOR_METALLU' TWORK.'PHASE' TWORK.'SOLUBILITE'
  1014.              TWORK.'FOR_DIFF' TWORK.'THM1' TWORK.'FOR_MECA') ;
  1015.     'MESS' ' > Adaptation du modele pas disponible pour cette formulation.' ;
  1016.     'ERRE' 251 ;
  1017.   'FINSI' ;
  1018. 'FINSI' ;
  1019. *
  1020. 'SI' (TWORK .'NVSTNL' 'ET' TWORK .'FOR_MECA') ;
  1021.   TWORK .'DYNAMIQUE' = VRAI ;
  1022. 'FINSI' ;
  1023. *
  1024. * Noms des inconnues passees a TRANSNON :
  1025. 'SI' ('NEG' TWORK.'MOD_TOT' 'INCONNU') ;
  1026.   TWORK.'PRIM_TOT' ='EXTR' TWORK.'MOD_TOT' 'DEPL' ;
  1027.   TWORK.'DUAL_TOT' ='EXTR' TWORK.'MOD_TOT' 'FORC' ;
  1028. 'FINS' ;
  1029. *
  1030. * GRANDS_DEPLACEMENTS si CONTACT :
  1031. 'SI' TWORK.'CONTACT' ;
  1032.   TWORK.'GRANDS_DEPLACEMENTS' = VRAI;
  1033. 'FINSI';
  1034. *
  1035. *---------- Calcul purement thermique / diffusion ? ou couplage? -------
  1036. 'SI' (TWORK.'FOR_THER' 'OU' TWORK.'FOR_DIFF');
  1037.   'SI' TWORK.'FOR_MECA' ;
  1038.     TWORK.'COUPLAGE' = VRAI;
  1039.     'SI' ('EXIS' PRECED 'PROJECTION');
  1040.       TWORK.'PROJECTION' = PRECED.'PROJECTION';
  1041.     'FINSI';
  1042.   'FINSI';
  1043. 'FINSI';
  1044. *
  1045. 'SI' TWORK.'POR1';
  1046.   TWORK.'MAI_POR' = 'EXTR' TWORK.'MOD_POR' 'MAIL';
  1047.   TWORK.'MOT_POR' = 'MOTS' 'FP' 'FPQ' 'FTP' ;
  1048. *
  1049. * poreux+consolidation -> il faut donner la matrice de permeabilite
  1050. *                         et ne pas etre en dynamique
  1051.   'SI' TWORK.'CONSOLIDATION' ;
  1052.     'SI' TWORK.'DYNAMIQUE' ;
  1053.       'ERREUR' ' On ne peut pas etre en dynamique.' ;
  1054.     'FINSI';
  1055.     'SI' ('NON' ('EXISTE' PRECED 'PERMEABILITE'));
  1056.       'ERREUR' 'Il faut donner la matrice de permeabilite';
  1057.     'SINON';
  1058.       TWORK.'PERMEABILITE'=PRECED.'PERMEABILITE';
  1059.     'FINSI';
  1060.     'SI' ('EXISTE' PRECED 'TETA' ) ;
  1061.       TWORK.'TETA' = PRECED.'TETA';
  1062.     'FINSI';
  1063.   'FINSI';
  1064. 'FINSI';
  1065. *
  1066. * Verif. CHARGEMENT si formulation CHARGEMENT :
  1067. 'SI' ('EXIS' TWORK 'MOD_CHA') ;
  1068.   'SI' ('NON' ('EXIS' (TWORK.'CHARGEMENT') 'PRES')) ;
  1069.     'ERREUR' 'Il manque la donnee d un chargement de pression';
  1070.     'ERREUR' '-> voir la notice de l operateur CHAR';
  1071.   'FINSI' ;
  1072. 'FINSI' ;
  1073.  
  1074. TWORK.'CHAR_THE'  = 'EXISTE' TWORK.'CHARGEMENT' 'T' ;
  1075. TWORK.'CHAR_DEFI' = 'EXISTE' TWORK.'CHARGEMENT' 'DEFI' ;
  1076. TWORK.'CHAR_MATE' = TWORK.'CHAR_MATE' 'OU' ('EXIS' TWORK.'CHARGEMENT' 'MATE') ;
  1077. TWORK.'CHAR_PRES' = 'EXISTE' TWORK 'MOD_PRE' ;
  1078. TWORK.'ITHER'     = 'OU' TWORK.'FOR_THER' TWORK.'CHAR_THE' ;
  1079.  
  1080. *------------ existe-t-il un chargement en pilotage indirect ?  --------
  1081. 'SI' ('EXISTE' PRECED 'PILOTAGE_INDIRECT') ;
  1082.    TWORK.'CHAR_PILO' = VRAI;
  1083.  
  1084. *--------- et des deplacements pilotes
  1085.   'SI' ('EXISTE' PRECED 'DEPLACEMENTS_PILOTES' );
  1086.      TWORK.'DEPLACEMENTS_PILOTES' = PRECED.'DEPLACEMENTS_PILOTES' ;
  1087.   'FINSI';
  1088.  
  1089. *--------- et des forces pilotees
  1090.   'SI' ('EXISTE' PRECED 'FORCES_PILOTEES' );
  1091.      TWORK.'FORCES_PILOTEES' = PRECED.'FORCES_PILOTEES';
  1092.   'FINSI';
  1093.  
  1094. *--------- et une region de recherche
  1095.   'SI' ('EXISTE' PRECED 'ZONE_DE_PILOTAGE' );
  1096.      TWORK.'ZONE_DE_PILOTAGE' = PRECED.'ZONE_DE_PILOTAGE';
  1097.   'FINSI';
  1098.  
  1099. *--------- et le parametre de pilotage
  1100.   'SI' ('EXISTE' PRECED 'PARAMETRE_DE_PILOTAGE' );
  1101.      TWORK.'PARAMETRE_DE_PILOTAGE' = PRECED.'PARAMETRE_DE_PILOTAGE';
  1102.   'FINSI';
  1103. * 
  1104. *----------  Parametre de pilotage indirect et table de sortie
  1105.   'SI' ('NON' ('EXIS' PRECED 'COEFFICIENT_DE_PILOTAGE'));
  1106.     PRECED.'COEFFICIENT_DE_PILOTAGE' = PROG 0. ;
  1107.   'FINSI';
  1108. 'FINSI';
  1109. *
  1110. 'SI' ('EXISTE' PRECED 'CAPACITE_CONSTANTE');
  1111.   TWORK.'CAPACONST' = VRAI;
  1112.   TWORK.'CAPACITE_CONSTANTE' = PRECED.'CAPACITE_CONSTANTE';
  1113. 'FINSI';
  1114. *
  1115. 'SI' ('EXISTE' PRECED 'CONDUCTIVITE_CONSTANTE');
  1116.   TWORK.'CONDCONST' = VRAI;
  1117.   TWORK.'CONDUCTIVITE_CONSTANTE' = PRECED.'CONDUCTIVITE_CONSTANTE';
  1118. 'FINSI';
  1119. ************************************************************************
  1120. * Verifications/initialisations propres a la THERMIQUE
  1121. ************************************************************************
  1122. 'SI' TWORK.'FOR_THER' ;
  1123. *
  1124.   'SI' ('OU' TWORK.'PHASE' TWORK.'SOLUBILITE' TWORK.'FOR_DIFF' TWORK.'THM1');
  1125.     TWORK.'PROCEDURE_THERMIQUE' = 'MOT' 'NONLINEAIRE';
  1126.   'FINSI';
  1127. *
  1128.   MOD_T  = TWORK.'MOD_T' ;
  1129.   PRIM_T ='EXTR' MOD_T 'DEPL' ;
  1130.  
  1131.  'SI' ('EXIS' PRECED 'TTOL');
  1132.     TTOL=PRECED.'TTOL' ;
  1133.    'SI' ('EGA' ('TYPE' TTOL) ('MOT' 'FLOTTANT'));
  1134.       TTOL ='MANU' 'CHPO' ('EXTR' MOD_T 'MAIL') PRIM_T ('PROG' ('DIME' PRIM_T) '*' TTOL) 'NATURE' 'DIFFUS';
  1135.    'SINO';
  1136.      'ERRE' 'PAS_DEFA : l''indice ''TTOL'' ne peut etre qu''un ''FLOTTANT''' ;
  1137.    'FINS';
  1138.    TWORK.'TTOL' = TTOL ;
  1139.  'FINS';
  1140.   TWORK.'PRIM_T' =  PRIM_T              'ET' 'LX'  ;
  1141.   TWORK.'DUAL_T' =('EXTR' MOD_T 'FORC') 'ET' 'FLX' ;
  1142. *
  1143.   'SI'('EXIS' PRECED 'CELSIUS');
  1144.     TWORK.'CELSIUS' = PRECED.'CELSIUS';
  1145.     TWORK.'TREF'    = 273.15 ;
  1146.   'FINSI';
  1147. *
  1148.   'SI' ('EXIS' PRECED 'BLOCAGES_THERMIQUES');
  1149.     BLOTHE1 = PRECED.'BLOCAGES_THERMIQUES' ;
  1150.     TWORK.'BLOCAGES_THERMIQUES' = BLOTHE1 ;
  1151.     TWORK.'CHAR_BLOT' = 'EGA' ('TYPE' BLOTHE1) 'CHARGEME' ;
  1152.   'FINSI';
  1153. *
  1154.   'SI'  TWORK.'RAYO' ;
  1155.     TWORK.'RAYONNEMENT' = 'PRRA' TWORK.'MOD_RAY';
  1156.     'SI' ('NON' ('EXISTE' (TWORK.'CARACTERISTIQUES') 'EMIS'));
  1157.        'ERREUR' ' Il manque la donnee des facteurs d emissivite';
  1158.     'FINSI';
  1159.     'SI'('EXISTE' PRECED 'CTE_STEFAN_BOLTZMANN') ;
  1160.        TWORK.'CTE_STEFAN_BOLTZMANN'=PRECED.'CTE_STEFAN_BOLTZMANN';
  1161.       'SAUT' 1 'LIGN' ;
  1162.       'MESS' 'ATTENTION ! Constante de Stefan-Boltzmann en unites SI';
  1163.       'MESS' '            soit ' (TWORK . 'CTE_STEFAN_BOLTZMANN') ;
  1164.     'FINSI' ;
  1165.   'FINSI';
  1166. *
  1167.   'SI' ('EXIS' PRECED 'RELAXATION_THETA');
  1168.     TWORK.'RELAXATION_THETA'=PRECED.'RELAXATION_THETA' ;
  1169.   'FINSI';
  1170. *
  1171.   'SI' ('EXIS' PRECED 'RELAXATION_DUPONT');
  1172.     TWORK.'RELAXATION_DUPONT'=  PRECED . 'RELAXATION_DUPONT' ;
  1173.   'FINSI';
  1174. *
  1175.   'SI' ('EXISTE' PRECED 'THERMIQUE_AUGMENTEE');
  1176.     TWORK.'THERAUGM' = VRAI;
  1177.     TWORK.'THERMIQUE_AUGMENTEE' = PRECED.'THERMIQUE_AUGMENTEE';
  1178.     RR= PRECED.'THERMIQUE_AUGMENTEE';
  1179.     MNOMU = ('EXTR' RR 'MAILLAGE' 'NOMU') 'NBNO';
  1180.     MAVMU = ('EXTR' RR 'MAILLAGE'       ) 'NBNO';
  1181.     'SI' (MNOMU 'NEG' MAVMU) ;
  1182.       'ERREUR' 'Il ne faut pas de multiplicateur dans THERMIQUE_AUGMENTEE';
  1183.     'FINSI';
  1184.   'FINSI';
  1185. 'FINSI';
  1186. ************************************************************************
  1187. * Verifications/initialisations propres a la DIFFUSION
  1188. ************************************************************************
  1189. 'SI' TWORK.'FOR_DIFF' ;
  1190. *
  1191.   MOD_DIF  = TWORK.'MOD_DIF' ;
  1192.   PRIM_DIF ='EXTR' MOD_DIF 'DEPL' ;
  1193.   TWORK.'PRIM_DIF' = PRIM_DIF ;
  1194.  
  1195.   'SI' ('EXIS' PRECED 'CTOL');
  1196.     CTOL=PRECED.'CTOL' ;
  1197.     'SI' ('EGA' ('TYPE' CTOL) ('MOT' 'FLOTTANT'));
  1198.       CTOL ='MANU' 'CHPO' ('EXTR' MOD_DIF 'MAIL') PRIM_DIF ('PROG' ('DIME' PRIM_DIF) '*' CTOL) 'NATURE' 'DIFFUS';
  1199.     'SINO';
  1200.       'SI' ('EGA' ('TYPE' CTOL) ('MOT' 'CHPOINT'));
  1201.         'SI' ('NON' ('EXIS' ('EXTR' CTOL 'COMP') PRIM_DIF 'ET'));
  1202.          'ERRE' 'CTOL ne contient pas toutes les composantes de la DIFFUSION';
  1203.         'FINS';
  1204.       'SINO';
  1205.        'ERRE' 'PAS_DEFA : l''indice ''CTOL'' ne peut etre qu''un CHPOINT ou un ''FLOTTANT''';
  1206.       'FINS';
  1207.     'FINS';
  1208.  
  1209.     'SI' ('EGA' ('TYPE' TWORK.'TTOL') ('MOT' 'CHPOINT'));
  1210.       TWORK.'TTOL' = TWORK.'TTOL' 'ET' CTOL ;
  1211.     'SINO';
  1212.       TWORK.'TTOL' = CTOL ;
  1213.     'FINS';
  1214.   'FINS';
  1215.   TWORK.'PRIM_DIF' =  PRIM_DIF              'ET' 'LX'  ;
  1216.   TWORK.'DUAL_DIF' =('EXTR' MOD_DIF 'FORC') 'ET' 'FLX' ;
  1217. *
  1218.   'SI' ('EXIS' PRECED 'BLOCAGES_DIFFUSIONS') ;
  1219.     BLODIF1 = PRECED.'BLOCAGES_DIFFUSIONS' ;
  1220.     TWORK.'BLOCAGES_DIFFUSIONS' = BLODIF1 ;
  1221.     TWORK.'CHAR_BLOD' = 'EGA' ('TYPE' BLODIF1) 'CHARGEME' ;
  1222.   'FINSI' ;
  1223. *
  1224.   'SI' ('EXIS' PRECED 'RELAXATION_THETA' );
  1225.     TWORK.'RELAXATION_THETA'=PRECED.'RELAXATION_THETA' ;
  1226.   'FINSI';
  1227. *
  1228.   'SI' ('EXIS'  PRECED  'RELAXATION_DUPONT');
  1229.     TWORK.'RELAXATION_DUPONT'=  PRECED . 'RELAXATION_DUPONT' ;
  1230.   'FINSI';
  1231. *
  1232.   'SI' ('EXISTE' PRECED 'DIFFUSION_AUGMENTEE');
  1233.     TWORK.'DIFFAUGM' = VRAI;
  1234.     TWORK.'DIFFUSION_AUGMENTEE' = PRECED.'DIFFUSION_AUGMENTEE';
  1235.     RR= PRECED.'DIFFUSION_AUGMENTEE';
  1236.     MNOMU = ('EXTR' RR 'MAILLAGE' 'NOMU') 'NBNO';
  1237.     MAVMU = ('EXTR' RR 'MAILLAGE'       ) 'NBNO';
  1238.     'SI' (MNOMU 'NEG' MAVMU) ;
  1239.       'ERREUR' 'Il ne faut pas de multiplicateur dans DIFFUSION_AUGMENTEE';
  1240.     'FINSI';
  1241.   'FINSI';
  1242. 'FINSI';
  1243. ************************************************************************
  1244. * Verifications/initialisations propres a la MECANIQUE
  1245. ************************************************************************
  1246. 'SI' ('EXIS' PRECED 'BLOCAGES_MECANIQUES');
  1247.   BLOMEC1 = PRECED.'BLOCAGES_MECANIQUES' ;
  1248.   TWORK.'BLOCAGES_MECANIQUES' = BLOMEC1 ;
  1249.   TWORK.'CHAR_BLOM' = 'EGA' ('TYPE' BLOMEC1) 'CHARGEME' ;
  1250. 'FINSI';
  1251. *
  1252. 'SI' TWORK.'FOR_MECA' ;
  1253. *
  1254.    'SI' ('EXISTE' PRECED 'PRECISION') ;
  1255.      TWORK.'PRECFLEX' = PRECED.'PRECISION' ;
  1256.    'FINSI';
  1257.    'SI' ('EXISTE' PRECED 'PRECFLEX') ;
  1258.      TWORK.'PRECFLEX' = PRECED.'PRECFLEX';
  1259.    'FINSI';
  1260. *
  1261.    'SI' ('EXISTE' PRECED 'K_TANGENT');
  1262.      TWORK.'K_TANGENT' = PRECED.'K_TANGENT';
  1263.    'FINSI';
  1264.    'SI' ('EXISTE' PRECED 'K_TANGENT_PERT') ;
  1265.      TWORK.'K_TANGENT_PERT' = PRECED.'K_TANGENT_PERT' ;
  1266.      'SI' TWORK.'K_TANGENT_PERT' ;
  1267.        TWORK.'K_TANGENT' = VRAI ;
  1268.      'FINSI' ;
  1269.    'FINSI' ;
  1270.    'SI' ('EXISTE' PRECED 'K_TANG_PERT_C1') ;
  1271.      TWORK.'K_TANG_PERT_C1' = PRECED.'K_TANG_PERT_C1' ;
  1272.    'FINSI' ;
  1273.    'SI' ('EXISTE' PRECED 'K_TANG_PERT_C2') ;
  1274.      TWORK.'K_TANG_PERT_C2' = PRECED.'K_TANG_PERT_C2' ;
  1275.    'FINSI' ;
  1276. *  Par defaut, C2 = C1/100 (C1 donne ou par defaut 1.D-3)
  1277.    'SI' ('EGA' TWORK.'K_TANG_PERT_C2' 'INCONNU') ;
  1278.      TWORK.'K_TANG_PERT_C2' = TWORK.'K_TANG_PERT_C1' / 1.D+02 ;
  1279.    'FINSI' ;
  1280. * Matrice tangente a utiliser pour iteration 0 dans UNPAS
  1281. * Par defaut, matrice tangente precedemment calculee si existe
  1282.    'SI' ('EXISTE' PRECED 'K_TANGENT_ITER0') ;
  1283.      TWORK.'K_TANGENT_ITER0' = PRECED.'K_TANGENT_ITER0' ;
  1284.    'FINSI' ;
  1285. * Utilisation matrice tangente symetrisee
  1286.    'SI' ('EXISTE' PRECED 'K_TANGENT_SYME') ;
  1287.      TWORK.'K_TANGENT_SYME' = PRECED.'K_TANGENT_SYME' ;
  1288.    'FINSI' ;
  1289.    'SI' (  'EXISTE' PRECED 'NITER_KTANGENT');
  1290.      TWORK.'NITER_KTANGENT' =PRECED.'NITER_KTANGENT' ;
  1291.    'FINSI';
  1292. *----------MAN? ordre de la MAN?
  1293.    'SI' ('EXISTE' PRECED 'MAN' );
  1294.      TWORK.'MAN'= PRECED.'MAN';
  1295.    'FINSI';
  1296.    'SI' ( 'EXISTE' PRECED 'ORDRE');
  1297.      TWORK.'ORDRE'=PRECED.'ORDRE';
  1298.    'FINSI';
  1299. *-------------y-a-t-il du substepping?-----------
  1300.    'SI' ('EXISTE' PRECED 'SUBSTEPPING');
  1301.      TWORK.'SUBSTEPPING'=PRECED.'SUBSTEPPING' ;
  1302.    'FINSI';
  1303.    'SI' ('EXISTE' PRECED 'PRECISINTER');
  1304.      TWORK.'PRECISINTER' = PRECED.'PRECISINTER';
  1305.    'FINSI';
  1306.    'SI' ('EXISTE' PRECED 'NITERINTER_MAX');
  1307.      TWORK.'NITERINTER_MAX' =PRECED.'NITERINTER_MAX';
  1308.    'FINSI';
  1309.    'SI' ('EXISTE' PRECED  'NMAXSUBSTEPS');
  1310.      TWORK.'NMAXSUBSTEPS'=PRECED.'NMAXSUBSTEPS';
  1311.    'FINSI';
  1312. *------------- quelle est l'hypothese de calcul des deformations ?
  1313.    'SI' ('EXISTE' PRECED 'HYPOTHESE_DEFORMATIONS' );
  1314.      HYPDEF = MOT PRECED.'HYPOTHESE_DEFORMATIONS';
  1315.      'SI' (('EGA' HYPDEF ('MOT' 'LINEAIRE'))  'OU' ('EGA' HYPDEF ('MOT' 'QUADRATIQUE')) 'OU'
  1316.            ('EGA' HYPDEF ('MOT' 'TRUESDELL')) 'OU' ('EGA' HYPDEF ('MOT' 'JAUMANN'))     'OU'
  1317.            ('EGA' HYPDEF ('MOT' 'UTILISATEUR')));
  1318.        TWORK.'HYPOTHESE_DEFORMATIONS'='MOT' HYPDEF;
  1319.      'SINON';
  1320.        'MESS' 'Les valeurs permises pour HYPOTHESE_DEFORMATIONS sont :';
  1321.        'MESS' 'LINEAIRE, QUADRATIQUE, TRUESDELL, JAUMANN, UTILISATEUR';
  1322.        'ERRE' 308;
  1323.      'FINSI';
  1324.    'FINSI';
  1325. *------------- y-a-t-il du grand deplacement ?
  1326.    'SI' ('EXISTE' PRECED 'GRANDS_DEPLACEMENTS' );
  1327.      'SI' (TWORK.'CONTACT' 'ET' ('NON' PRECED.'GRANDS_DEPLACEMENTS'));
  1328.        'MESS' 'Un calcul avec un modele de CONTACT doit necessairement'
  1329.        'MESS' 'etre realise en GRANDS_DEPLACEMENTS !';
  1330.        'MESS' 'Arret de pasapas';
  1331.        'ERRE' 21;
  1332.      'FINSI';
  1333.      TWORK.'GRANDS_DEPLACEMENTS' = PRECED.'GRANDS_DEPLACEMENTS';
  1334.    'FINSI';
  1335. *------------- quelques reglages forces en grands deplacements
  1336.    'SI' TWORK.'GRANDS_DEPLACEMENTS' ;
  1337.      TWORK.'RECARI'   =  VRAI ;
  1338. *
  1339. *  PV Quel interet d'avoir les deformations?
  1340. **   par defaut, on sauvegarde des def. totales en grands deplacements
  1341. **   'SI' ('NON' ('EXIS' PRECED 'MES_SAUVEGARDES'));
  1342. **     TWORK.'TABSAUV'.'DEFTO'= VRAI;
  1343. **   'FINS';
  1344. *    par defaut, on calcule avec une deformation quadratique en grands deplacements
  1345. *    mais on laisse faire l'utilisateur s'il souhaite autrement
  1346.      'SI' ('NON' ('EXIS' PRECED 'HYPOTHESE_DEFORMATIONS'));
  1347.        TWORK.'HYPOTHESE_DEFORMATIONS'='MOT' 'QUADRATIQUE';
  1348.      'FINSI';
  1349. *    par defaut, on active K_SIGMA en grands deplacements si pas de controle de stabilite
  1350.      'SI' ('NON' TWORK.'STABILITE'); 
  1351.        TWORK.'K_SIGMA'=VRAI;
  1352.      'FINSI';
  1353. *    par defaut, on active K_SIGMA en grands deplacements si deformations quadratiques    
  1354.      'SI' ('EGA' TWORK.'HYPOTHESE_DEFORMATIONS' 'QUADRATIQUE');
  1355.        TWORK.'K_SIGMA'=VRAI;
  1356.      'FINSI';
  1357.  
  1358. *------Y-a-t-il des elements dont il faudra reactualiser les caracteristiques?
  1359.      MOCA = 'MOTS' 'VECT' 'VX  ' 'VY  ' 'VZ  ' 'VXF ' 'VYF ' 'VZF '
  1360.                    'V1X ' 'V1Y ' 'V1Z ' 'V2X ' 'V2Y ' 'V2Z ';
  1361.      'REPETER' BOU_CAR ('DIME' MOCA);
  1362.        ITE2 = &BOU_CAR;
  1363.        MOT1 = 'EXTR' MOCA ITE2;
  1364.        'SI' ('EXIS' TWORK.'MA_TOT' MOT1);
  1365.          TWORK.'ITCAR' = VRAI;
  1366.          'QUITTER' BOU_CAR;
  1367.        'FINSI';
  1368.      'FIN' BOU_CAR;
  1369.    'FINSI';
  1370. * ---------- y-a-t-il une demande sur K_SIGMA ?
  1371.    'SI' ('EXISTE' PRECED 'K_SIGMA');
  1372.      TWORK.'K_SIGMA' = PRECED.'K_SIGMA';
  1373.    'FINSI';
  1374. * ---------- y-a-t-il une demande sur LAGRANGIEN ?
  1375.    'SI' ('EXISTE' PRECED 'LAGRANGIEN');
  1376.      'SI' ('EGA' preced.'LAGRANGIEN' 'TOTAL');
  1377.        TWORK.'LAG_TOT' = 1;
  1378.      'SINON';
  1379.        'SI' ('EGA' preced.'LAGRANGIEN' 'FIN_PAS');
  1380.          TWORK.'LAG_TOT' = 2;
  1381.        'SINON';
  1382.          'SI' (('EGA' preced.'LAGRANGIEN' 'REACTUALISE') 'OU'
  1383.                ('EGA' preced.'LAGRANGIEN' 'MI_PAS'));
  1384.            TWORK.'LAG_TOT' = 3;
  1385.          'SINON';
  1386.            'ERRE' 'Seules les valeurs TOTAL, REACTUALISE (a MI_PAS) et FIN_PAS sont permises';
  1387.          'FINSI';
  1388.        'FINSI';
  1389.      'FINSI';
  1390.    'FINSI';
  1391. *
  1392. *  MOdele TOTal si FOR_MECA :
  1393.    MO_TOT = TWORK.'MO_TOT' ;
  1394. *--------- initialisation des comportements d'apres le modele
  1395.    'SI' ('EXISTE' MO_TOT 'MATE' 'PLASTIQUE');
  1396.      TWORK.'PLASTIQUE'      = VRAI;
  1397.      TWORK.'IPLAVI'         = VRAI;
  1398.    'FINSI';
  1399.    'SI' ('EXISTE' MO_TOT 'MATE' 'PLASTIQUE_ENDOM');
  1400.      TWORK.'PLASTIQUE'      = VRAI;
  1401.      TWORK.'IPLAVI'         = VRAI;
  1402.      TWORK.'ENDOMMAGEMENT'  = VRAI;
  1403.    'FINSI';
  1404.    'SI' (('EXISTE' MO_TOT 'MATE' 'VISCOPLASTIQUE') 'OU'
  1405.          ('EXISTE' MO_TOT 'MATE' 'FLUAGE'        ));
  1406.      TWORK.'VISCOPLASTIQUE' = VRAI;
  1407.      TWORK.'IPLAVI'         = VRAI;
  1408.    'FINSI';
  1409.    'SI' ('EXISTE' MO_TOT 'MATE' 'VISCODOMMAGE');
  1410.      TWORK.'VISCODOMMAGE'   = VRAI;
  1411.      TWORK.'IPLAVI'         = VRAI;
  1412.    'FINSI';
  1413.    'SI' (('EXISTE' MO_TOT 'MATE' 'ENDOMMAGEMENT') 'OU'
  1414.         ('EXISTE' MO_TOT 'MATE' 'ENDOMMAGEABLE'));
  1415.      TWORK.'ENDOMMAGEMENT'  = VRAI;
  1416.      TWORK.'IPLAVI'         = VRAI;
  1417.    'FINSI';
  1418.    'SI' ('EXISTE' MO_TOT 'MATE' 'NON_LINEAIRE') ;
  1419.      TWORK.'NON_LINEAIRE'   = VRAI;
  1420.      TWORK.'IPLAVI'         = VRAI;
  1421.      TWORK.'UTILISATEUR'    = 'EXISTE' MO_TOT 'MATE' 'UTILISATEUR' ;
  1422.    'FINSI';
  1423.    'SI' ('EXISTE' MO_TOT 'MATE' 'VISCO_EXTERNE') ;
  1424.      TWORK.'VISCO_EXTERNE'  = VRAI;
  1425.      TWORK.'IPLAVI'         = VRAI;
  1426.    'FINSI';
  1427.    'SI' ('EXISTE' MO_TOT 'MATE' 'CERAMIQUE');
  1428.      TWORK.'CERAMIQUE'      = VRAI;
  1429.      TWORK.'IPLAVI'         = VRAI;
  1430.    'FINSI';
  1431.    'SI' (('EXISTE' MO_TOT 'MATE' 'MODAL') 'OU'
  1432.              ('EXISTE' MO_TOT 'MATE' 'STATIQUE' ));
  1433.      TWORK.'MODAL'          = VRAI;
  1434.      TWORK.'IPLAVI'         = VRAI;
  1435.    'FINSI';
  1436.    'SI' ('EXISTE' MO_TOT 'MATE' 'IMPEDANCE');
  1437.      TWORK.'IMPEDANCE'      = VRAI;
  1438.      TWORK.'IPLAVI'         = VRAI;
  1439.    'FINSI';
  1440. *
  1441.    'SI' ('OU' TWORK.'CERAMIQUE' TWORK.'ENDOMMAGEMENT' TWORK.'VISCODOMMAGE');
  1442.      TWORK.'RECARI'  = VRAI;
  1443.      TWORK.'LAG_TOT' = 0 ;
  1444.    'FINS';
  1445. *
  1446. *- existe-t-il un chargement de flux neutronique pour le modele mistral
  1447. *- existe-t-il un chargement de densite de fissions pour gatt_monerie & uo2
  1448. *----  si oui on sort si non on le cree a zero
  1449.    'SI' ('EXIS'  MO_TOT  'MATE' 'MISTRAL');
  1450.      'SI' ('NON' ('EXISTE' (TWORK.'CHARGEMENT') 'FI  ')) ;
  1451.         FIZER = 'MANU' 'CHPO' MAI_TOT 1 'FI  ' 0. 'NATURE' 'DIFFUS' ;
  1452.         TWORK.'CHARGEMENT' = TWORK.'CHARGEMENT' 'ET' ('CHAR' 'FI' FIZER
  1453.       ('EVOL' 'MANU'  ('PROG' 0. 10.e10 ) ('PROG' 1. 1.)));
  1454.      'FINSI';
  1455.    'FINSI';
  1456.    'SI' (('EXIS'  MO_TOT  'MATE' 'UO2') 'OU'
  1457.          ('EXIS'  MO_TOT  'MATE' 'GATT_MONERIE') ) ;
  1458.      'SI' ('NON' ('EXISTE' (TWORK.'CHARGEMENT') 'DFIS')) ;
  1459.        DFISZER = 'MANU' 'CHPO'  MAI_TOT 1 'DFIS' 0. 'NATURE' 'DIFFUS' ;
  1460.        TWORK.'CHARGEMENT' = TWORK.'CHARGEMENT' 'ET' ('CHAR' 'DFIS' DFISZER
  1461.          ('EVOL' 'MANU'   ('PROG' 0. 10.e10 ) ('PROG' 1. 1.))) ;
  1462.      'FINSI';
  1463.    'FINSI';
  1464. *
  1465.    MAI_TOT             = 'EXTR' MO_TOT 'MAIL'       ;
  1466.    TWORK.'LDEVA'       = 'EXTR' TWORK.'MA_TOT' 'DEVA';
  1467.    TWORK.'LPAR0'       = 'EXTR' MO_TOT 'PARA'       ;
  1468.    'SI' ('NEG' 0 ('DIME'  TWORK.'LPAR0'))     ;
  1469.      TWORK.'PARA_EXTE'=VRAI;
  1470.    'FINSI';
  1471. *--------------y-a-t-il une formulation feFP
  1472.    'SI' ('EXISTE' PRECED 'FEFP_FORMULATION');
  1473.      TWORK.'FEFP_FORMULATION' = PRECED.'FEFP_FORMULATION';
  1474.    'FINSI';
  1475.    'SI' TWORK.'FEFP_FORMULATION';
  1476.      TWORK.'GRANDS_DEPLACEMENTS'  = FAUX ;
  1477.      TWORK.'K_SIGMA'              = VRAI ;
  1478.      TWORK.'K_TANGENT'            = VRAI ;
  1479.      TWORK.'SUBSTEPPING'          = FAUX ;
  1480.      TWORK.'RECARI'               = VRAI ;
  1481.    'FINSI';
  1482.    'SI' ( 'EXISTE' PRECED 'UPDATE_LAGRANGIAN'  );
  1483.      TWORK.'UPDATE_LAGRANGIAN' = PRECED.'UPDATE_LAGRANGIAN' ;
  1484.    'FINSI';
  1485.    'SI' ('NON' TWORK.'GRANDS_DEPLACEMENTS');
  1486.      TWORK.'LAG_TOT' = 0 ;
  1487.    'FINS';
  1488. *---------- en cas de pilotage automatique ?
  1489.    'SI' TWORK.'AUTOMATIQUE' ;
  1490.      'SI' ('EXISTE' PRECED 'AUTOCRIT');
  1491.        TWORK.'AUTOCRIT'= PRECED.'AUTOCRIT';
  1492.      'FINSI';
  1493.      'SI' ('EXISTE' PRECED 'AUTOPAS');
  1494.        TWORK.'AUTOPAS'= PRECED.'AUTOPAS';
  1495.      'FINSI';
  1496.      'SI' ('EXISTE' PRECED 'AUTORECA');
  1497.        TWORK.'AUTORECA'= 'MAXI' ('ENTI' PRECED.'AUTORECA') 1 ;
  1498.      'FINSI';
  1499.      'SI' ('EXISTE' PRECED 'AUTORESU');
  1500.        TWORK.'NAURES' = PRECED.'AUTORESU';
  1501.      'FINSI';
  1502.      TWORK.'DELTAITER' = TWORK.'MAXITERATION';
  1503.    'FINSI';
  1504. *------ y-a-t-il un nom de variable interne specifique ?
  1505.    'SI' ('EXISTE' PRECED 'MOVA');
  1506.      TWORK.'MOVA' =PRECED.'MOVA';
  1507.    'SINON';
  1508.      'SI'(TWORK.'VISCO_EXTERNE' 'OU' TWORK.'UTILISATEUR');
  1509.        TWORK.MOVA = 'MOT' 'RIEN';
  1510.      'FINSI';
  1511.      'SI' ('NON' ('EXISTE' ('EXTR' TWORK.'MODELE' 'VARI') 'EPSE'));
  1512.        TWORK.'MOVA' = 'MOT' 'RIEN';
  1513.      'FINSI'  ;
  1514.    'FINSI';
  1515. *---------- y-a-t-il des elements de rigidite constante ? --------
  1516.    'SI' ( 'EXISTE' PRECED 'RIGIDITE_CONSTANTE');
  1517.      TWORK.'RAIDCONST' =VRAI;
  1518.      TWORK.'RIGIDITE_CONSTANTE'=PRECED.'RIGIDITE_CONSTANTE';
  1519.      RR= PRECED.'RIGIDITE_CONSTANTE';
  1520.      MNOMU = ('EXTR' RR 'MAILLAGE' 'NOMU') 'NBNO';
  1521.      MAVMU = ('EXTR' RR 'MAILLAGE' ) 'NBNO';
  1522.      'SI' ( MNOMU 'NEG' MAVMU) ;
  1523.        'ERREUR' 'Il ne faut pas de multiplicateur dans RIGIDITE_CONSTANTE';
  1524.      'FINSI';
  1525.    'FINSI';
  1526. *---------- y-a-t-il des elements de rigidite augmentee ? --------
  1527.    'SI' ('EXISTE' PRECED 'RIGIDITE_AUGMENTEE');
  1528.      TWORK.'RAIDAUGM' =VRAI;
  1529.      TWORK.'RIGIDITE_AUGMENTEE'=PRECED.'RIGIDITE_AUGMENTEE';
  1530.      RR= PRECED.'RIGIDITE_AUGMENTEE';
  1531.      MNOMU = ('EXTR' RR 'MAILLAGE' 'NOMU') 'NBNO';
  1532.      MAVMU = ('EXTR' RR 'MAILLAGE' ) 'NBNO';
  1533.      'SI' ( MNOMU 'NEG' MAVMU) ;
  1534.        'ERREUR' 'Il ne faut pas de multiplicateur dans RIGIDITE_AUGMENTEE';
  1535.      'FINSI';
  1536.    'SINON';
  1537.      'SI' ('EXISTE' PRECED 'AUGMENTATION_AUTOMATIQUE');
  1538.         TWORK.'RIGIDITE_AUGMENTEE'='MASSE' TWORK.'MODELE' TWORK.'CARACTERISTIQUES';
  1539.      'FINSI';
  1540.    'FINSI';
  1541.  
  1542.    'SI' ('EXISTE' PRECED 'AUGMENTATION_AUTOMATIQUE');
  1543.      'SI' PRECED.'AUGMENTATION_AUTOMATIQUE';
  1544.        TWORK.'AUTOAUGM'  = VRAI;
  1545.        TWORK.'STABILITE' = VRAI ;
  1546. ***  'SI' ('NON' TWORK.'RAIDAUGM');
  1547. ***    'ERREUR' 'il faut donner la matrice d augmentation dans RIGIDITE_AUGMENTEE';
  1548. ***  'FINSI';
  1549.      'FINSI';
  1550.    'FINSI';
  1551. *---------- y-a-t-il des elements de masse constante ? ---------
  1552.    'SI' ('EXISTE' PRECED 'MASSE_CONSTANTE');
  1553.      TWORK.'MASSE_CONSTANTE'=PRECED.'MASSE_CONSTANTE';
  1554.      TWORK.'MASSCONST'=VRAI;
  1555.    'FINSI';
  1556. *---------- y-a-t-il un sous_increment ? ---------
  1557.    'SI' ('EXISTE' PRECED 'SOUS_INCREMENT');
  1558.      TWORK.'SOUS_INCREMENT'   = PRECED.'SOUS_INCREMENT';
  1559.      TWORK.'RECALCUL_SOUSINC' = FAUX ;
  1560.      'SI' ('&lt;EG' TWORK.'SOUS_INCREMENT' 0) ;
  1561.        'MESS' 'Indice SOUS_INCREMENT doit etre > 0' ;
  1562.        'ERREUR' 21 ;
  1563.      'FINSI';
  1564.    'FINSI';
  1565. *------------------- Cas du non local :
  1566.    'SI' TWORK.'NLOC' ;
  1567.      'SI' ( 'EXISTE'  PRECED 'CONN' );
  1568.        'MESS' 'On garde les connectivites dans TAB CONN';
  1569.        TWORK.'CONN'=PRECED . 'CONN' ;
  1570.      'SINON';
  1571.        'SI' ('EGA'  (TWORK.'NON_LOCAL')  'HELM' ); 
  1572.          'MESS' 'Pas besoin de la connectivite pour Helmholtz';
  1573.        'SINON';  
  1574.          'MESS' 'Calcul des connectivites en mode normal.';
  1575.          TWORK.'CONN'='CONN' TWORK.'NLOC_MODL'
  1576.             ('REDU' TWORK.'MA_TOT' TWORK.'NLOC_MODL') 'NORMAL' ;
  1577.        'FINSI' ;
  1578.      'FINSI' ;
  1579.    'FINSI' ;
  1580. 'FINSI';
  1581. TWORK.'CHARGEMENT_I' = TWORK.'CHARGEMENT';
  1582. ************************************************************************
  1583. * Recapitulatif des options choisies
  1584. ************************************************************************
  1585. 'SAUTER' 1 'LIGNE';
  1586. 'SI' TWORK.'FOR_MECA' ;
  1587.   'SI' TWORK.'FOR_THER';
  1588.     'MESS' 'Couplage THERMIQUE-MECANIQUE';
  1589.   'SINON';
  1590.     'MESS' 'Calcul MECANIQUE';
  1591.   'FINSI';
  1592.   'MESS' ' *** Hypothese calcul des deformations :' ' ' TWORK.'HYPOTHESE_DEFORMATIONS' ' ***';
  1593.   'SI' ('NEG' TWORK.'PREDICTEUR' 'INCONNU');
  1594.     'MESS'   ' *** Predicteur :' ' ' TWORK.'PREDICTEUR' ' ***' ;
  1595.   'FINSI';
  1596.   'SI' TWORK.'K_SIGMA';
  1597.     'MESS' ' *** K_SIGMA *** ';
  1598.   'FINSI';
  1599.   'SI' TWORK.'GRANDS_DEPLACEMENTS';
  1600.     'SI' (TWORK.'LAG_TOT' 'EGA' 1) ;
  1601.       'MESS' ' *** Calcul en lagrangien total *** ';
  1602.     'FINSI';
  1603.     'SI' (TWORK.'LAG_TOT' 'EGA' 2) ;
  1604.       'MESS' ' *** Calcul en lagrangien fin pas *** ';
  1605.     'FINSI';
  1606.     'SI' (TWORK.'LAG_TOT' 'EGA' 3) ;
  1607.       'MESS' ' *** Calcul en configuration mi-pas *** ';
  1608.     'FINSI';
  1609.   'FINSI';
  1610.   'SI' TWORK.'STABILITE';
  1611.     'MESS' ' *** Controle de la stabilite *** ';
  1612.   'FINSI';
  1613.   'SI' TWORK.'LINESEARCH';
  1614.     'MESS' ' *** Linesearch  *** ';
  1615.   'FINSI';
  1616.   'SI' (non TWORK.'INITIALISATION');
  1617.     'MESS' ' *** Initialisation a partir du pas precedent interdite *** ';
  1618.   'FINSI';
  1619.   'SI' TWORK.'CONVERGENCE_MONOTONE' ;
  1620.     'MESS' ' *** Convergence monotone imposee *** ';
  1621.   'FINSI';
  1622.   'SI' TWORK.'RENORMALISATION' ;
  1623.     'MESS' ' *** Renormalisation des iteres *** ';
  1624.   'FINSI';
  1625.   'SI' TWORK.'POR1';
  1626.     'MESS' ' *** POREUX ***';
  1627.   'FINSI';
  1628.   'SI' TWORK.'DYNAMIQUE' ;
  1629.     'MESS' ' *** DYNAMIQUE PAS A PAS ***';
  1630.   'FINSI' ;
  1631.   'SI' TWORK.'FREQUENTIEL'  ;
  1632.     'MESS' ' *** FREQUENTIEL PAS A PAS ***' ;
  1633.   'FINSI' ;
  1634.   'SI' TWORK.'LIAISON_PERSISTANTE' ;
  1635.     'MESS' ' *** LIAISON PERSISTANTE ***';
  1636.   'FINSI' ;
  1637.   'SI' TWORK.'CONSOLIDATION' ;
  1638.     'MESS' ' *** CONSOLIDATION ***';
  1639.   'FINSI' ;
  1640.   'SI' TWORK.'PLASTIQUE' ;
  1641.     'MESS' ' *** PLASTICITE ***';
  1642.   'FINSI';
  1643.   'SI' TWORK.'VISCOPLASTIQUE' ;
  1644.     'MESS' ' *** VISCO - PLASTICITE ***';
  1645.   'FINSI' ;
  1646.   'SI' TWORK.'NON_LINEAIRE';
  1647.     'SI' TWORK.'UTILISATEUR' ;
  1648.       'MESS' ' *** MODELE DE COMPORTEMENT NON LINEAIRE EXTERNE ***';
  1649.     'SINON';
  1650.       'MESS' ' *** ELASTICITE NON_LINEAIRE ***';
  1651.     'FINSI';
  1652.   'FINSI' ;
  1653.   'SI' TWORK.'VISCO_EXTERNE';
  1654.     'MESS' ' *** VISCO - PLASTICITE, MODELE EXTERNE ***';
  1655.   'FINSI' ;
  1656.   'SI' TWORK.'CHAR_THE';
  1657.     'MESS' ' *** CHARGEMENT THERMIQUE ***';
  1658.   'FINSI';
  1659.   'SI' TWORK.'CHAR_DEFI';
  1660.     'MESS' ' *** CHARGEMENT PAR DES DEFORMATIONS ***';
  1661.   'FINSI';
  1662.   'SI' TWORK.'CHAR_PRES';
  1663.     'MESS' ' *** CHARGEMENT PAR UNE PRESSION SUIVEUSE ***';
  1664.   'FINSI';
  1665.   'SI' TWORK.'CHAR_PILO';
  1666.     'MESS' ' *** CHARGEMENT PAR PILOTAGE INDIRECT ***';
  1667.   'FINSI';
  1668.   'SI' TWORK.'PARA_EXTE';
  1669.     'MESS' ' *** LE MODELE DEPEND DE PARAMETRES EXTERNES ***';
  1670.   'FINSI';
  1671.   'SI' TWORK.'GRANDS_DEPLACEMENTS' ;
  1672.     'MESS' ' *** GRANDS DEPLACEMENTS ***';
  1673.     'SI' ('NEG' TWORK.'HYPOTHESE_DEFORMATIONS' 'QUADRATIQUE');
  1674.       'MESS' '     ATTENTION : l''hypothese de deformations choisie est' ' ' PRECED.'HYPOTHESE_DEFORMATIONS';
  1675.       'MESS' '                 on recommande plutot l''hypothese QUADRATIQUE';
  1676.     'FINSI';
  1677.   'FINSI';
  1678.   'SI' TWORK.'CONVERGENCE_FORCEE' ;
  1679.     'MESS' ' *** CONVERGENCE FORCEE   ***';
  1680.   'FINSI';
  1681.   'SI' TWORK.'ENDOMMAGEMENT' ;
  1682.     'MESS' ' *** PLASTICITE + ENDOMMAGEMENT ***';
  1683.   'FINSI';
  1684.   'SI' TWORK.'VISCODOMMAGE' ;
  1685.     'MESS' ' *** VISCO-PLASTICITE + ENDOMMAGEMENT ***';
  1686.   'FINSI';
  1687.   'SI' TWORK.'CERAMIQUE' ;
  1688.     'MESS' ' ***CERAMIQUE FLUAGE  + ENDOMMAGEMENT ***';
  1689.     'MESS' ' ***             +OTTOSEN             ***';
  1690.   'FINSI';
  1691. *   'SI' IUO2;
  1692. *      'MESS' ' *** VISCO-PLASTICITE + OTTOSEN ***';
  1693. *   'FINSI';
  1694.   'SI' TWORK.'NLOC' ;
  1695.     'MESS' ' *** EN NON-LOCAL *** ' ;
  1696.   'FINSI';
  1697.   'SI' TWORK.'CONTACT';
  1698.     'MESS' ' *** CONTACT ***';
  1699.   'FINSI';
  1700.   'SI' (TWORK.'SOUS_INCREMENT' > 1);
  1701.     'MESS' ' *** SOUS_INCREMENT :'TWORK.'SOUS_INCREMENT';
  1702.   'FINSI';
  1703.   'SI' TWORK.'CAFROTTE';
  1704.     'MESS' ' *** FROTTEMENT ***';
  1705.   'FINSI';
  1706.   'SI' TWORK.'K_TANGENT' ;
  1707.     'SI' TWORK.'K_TANGENT_PERT' ;
  1708.       'MESS' ' *** MATRICE TANGENTE PAR PERTURBATION ***' ;
  1709.     'SINON' ;
  1710.       'MESS' ' *** MATRICE TANGENTE (COHERENTE) ***' ;
  1711.     'FINSI' ;
  1712.     'MESS' '     MATRICE EVALUEE TOUTES LES' TWORK.'NITER_KTANGENT'
  1713.             'ITERATIONS' ;
  1714.     'SI' TWORK.'K_TANGENT_SYME' ;
  1715.       'MESS' '     PARTIE SYMETRIQUE DE LA MATRICE TANGENTE UTILISEE';
  1716.     'FINSI' ;
  1717.   'FINSI';
  1718.   'SI' TWORK.'SUBSTEPPING';
  1719.     'MESS' ' *** SUBSTEPPING *** NUM_MAXI_SUBSTEPS ='TWORK.'NMAXSUBSTEPS' ;
  1720.     'MESS' ' *** SUBSTEPPING *** NUM_ITER_INTERNAS ='TWORK.'NITERINTER_MAX' ;
  1721.   'FINSI';
  1722.    'SI' TWORK.'UPDATE_LAGRANGIAN';
  1723.      'MESS' ' *** UPDATE LAGRANGIAN ***' ;
  1724.    'FINSI';
  1725.    'SI' TWORK.'FEFP_FORMULATION';
  1726.      'MESS' ' *** FORMULATION FeFp ***' ;
  1727.    'FINSI';
  1728. 'FINSI';
  1729. *
  1730. 'SI' TWORK.'FOR_THER';
  1731.   'MESS' 'Calcul THERMIQUE ';
  1732. 'FINSI';
  1733. *
  1734. 'SI' TWORK.'NVSTNL' ;
  1735.   'MESS' 'Calcul NAVIER-STOKES NLIN' ;
  1736. 'FINSI' ;
  1737. 'SAUTER' 1 'LIGNE';
  1738. 'SI' ('EGA' TWORK.'PROCESSEURS' 'MONO_PROCESSEUR') ;
  1739.   'MESS' ' *** Pas de parallelisation ***' ;
  1740. 'SINON' ;
  1741.   'SI' ('EGA' TWORK.'PROCESSEURS' 'COMPORTEMENT') ;
  1742.     'MESS' ' *** Parallelisation du comportement sur ' nbpart ' assistants ***';
  1743.   'SINON' ;
  1744.     'SI' ('EGA' TWORK.'PROCESSEURS' 'AUTOMATIQUE') ;
  1745.       'MESS' ' *** Parallelisation automatique sur ' nbpart ' assistants ***' ;
  1746.     'FINSI' ;
  1747.   'FINSI' ;
  1748. 'FINSI' ;
  1749. 'SAUTER' 1 'LIGNE';
  1750. 'FINPROC' ;
  1751.  
  1752.  

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