Télécharger unpas.procedur

Retour à la liste

Numérotation des lignes :

  1. * UNPAS PROCEDUR PV090527 24/06/10 21:15:02 11936
  2. 'DEBPROC' UNPAS PRECED*'TABLE';
  3. *----------------------------------------------------------------------*
  4. * PROCEDURE UNPAS *
  5. * *
  6. * Calcul d'un increment de solution en grand deplacement plastique par *
  7. * la methode des residus. *
  8. * *
  9. * Les differentes configurations qui interviennent sont : *
  10. * WTAB.'FOR0' : configuration de debut de calcul *
  11. * WTAB.'GE0_DEB' : *
  12. * -> en GRANDS_DEPLACEMENTS ou FEFP_FORMULATION (1) : *
  13. * configuration de debut de pas *
  14. * -> dans les autres cas (2) : *
  15. * configuration de debut de calcul (= WTAB.'FOR0) *
  16. * Dans le cas (1), la configuration WTAB.'GE0_DEB' est actualisee *
  17. * dans PASAPAS avec GEOM2. *
  18. * Dans le cas (2), il n'y a pas d'actualisation de la configuration *
  19. * au cours du calcul. *
  20. * *
  21. * Pour les calculs DYNAMIQUE, un schema de Newmark (implicite) *
  22. * avec gamma = 1/2 et beta = 1/4 est utilise. *
  23. * *
  24. *----------------------------------------------------------------------*
  25. *CB215821 : Recuperation de XPETIT (07/12/2016)
  26. XPETIT = 'VALE' 'PETI' ;
  27. XGRAND = 'VALE' 'GRAND' ;
  28. *
  29. WTAB = PRECED.'WTABLE' ;
  30. conti = PRECED.'CONTINUATION';
  31. estim = PRECED.'ESTIMATION' ;
  32. *
  33. * Liste de composantes utiles
  34. MXMYMZ = 'MOTS' 'MX' 'MY' 'MZ' 'MT' 'FP' 'FPQ' 'FTP'
  35. 'IMX' 'IMY' 'IMZ' 'IMT' 'IFP' 'IFPQ' 'IFTP' ;
  36. MXMFLX = 'MOTS' 'MX' 'MY' 'MZ' 'MT' 'FLX' 'FP' 'FPQ' 'FTP'
  37. 'IMX' 'IMY' 'IMZ' 'IMT' 'FLX' 'IFP' 'IFPQ' 'IFTP' ;
  38. MNPRIM = 'MOTS' 'UX' 'UY' 'UZ' 'UR' 'UT' 'RX' 'RY' 'RZ' 'RT'
  39. 'P' 'PQ' 'TP' 'ALFA' 'BETA'
  40. 'IUX' 'IUY' 'IUZ' 'IUR' 'IUT' 'IRX' 'IRY' 'IRZ' 'IRT'
  41. 'IP' 'IPQ' 'ITP' 'IALF' 'IBET' ;
  42. MLPRIM = MNPRIM 'ET' ('MOTS' 'LX');
  43. MNDUAL = 'MOTS' 'FX' 'FY' 'FZ' 'FR' 'FT' 'MX' 'MY' 'MZ' 'MT'
  44. 'FP' 'FPQ' 'FTP' 'FALF' 'FBET'
  45. 'IFX' 'IFY' 'IFZ' 'IFR' 'IFT' 'IMX' 'IMY' 'IMZ' 'IMT'
  46. 'IFP' 'IFPQ' 'IFTP' 'IFAL' 'IFBE' ;
  47. MLDUAL = MNDUAL 'ET' ('MOTS' 'FLX');
  48. MLDEPL = 'MOTS' 'UX' 'UY' 'UZ' 'UR' 'UT' 'IUX' 'IUY' 'IUZ' 'IUR' 'IUT'
  49. 'ALFA' 'BETA' 'IALF' 'IBET' ;
  50. MLROTA = 'MOTS' 'RX' 'RY' 'RZ' 'RT' 'P' 'PQ' 'TP'
  51. 'IRX' 'IRY' 'IRZ' 'IRT' 'IP' 'IPQ' 'ITP';
  52. MLDEFOR = 'MOTS' 'EPXX' 'EPYY' 'EPZZ' 'EPSS' 'EPTT' 'EPRR'
  53. 'GAXY' 'GAXZ' 'GAYZ' 'GAST' 'GASN' 'GATN'
  54. 'GARZ' 'GART' 'GAZT' 'GXY '
  55. * 'CX ' 'CY ' 'CZ '
  56. 'EPSE' 'EPS' ;
  57. * Pour MOD_LIAISON :
  58. MLPRIM_LIA = MLPRIM ;
  59. MVPRIM_LIA = 'MOTS' 'VTX' 'VTY' 'VTZ' 'VTR' 'VTT' 'VWX' 'VWY' 'VWZ' 'VWT'
  60. 'VVP' 'VVPQ' 'VVTP' 'VALF' 'VBET'
  61. 'IVTX' 'IVTY' 'IVTZ' 'IVTR' 'IVTT' 'IVWX' 'IVWY' 'IVWZ' 'IVWT'
  62. 'IVVP' 'IVPQ' 'IVTP' 'IVAL' 'IVBE' 'VLX';
  63. MLDUAL_LIA = MLDUAL ;
  64.  
  65. *=DEB===== Formulation HHO = Definition des DDLs ===============================
  66. * Definition des DDLs primaux : inconnues de deplacements des cellules et des Faces
  67. MNPRIM_HHO = 'MOTS'
  68. * 'UXC0' 'UXC1' 'UXC2' 'UXC3' 'UXC4' 'UXC5' 'UXC6' 'UXC7' 'UXC8' 'UXC9'
  69. * 'UYC0' 'UYC1' 'UYC2' 'UYC3' 'UYC4' 'UYC5' 'UYC6' 'UYC7' 'UYC8' 'UYC9'
  70. * 'UZC0' 'UZC1' 'UZC2' 'UZC3' 'UZC4' 'UZC5' 'UZC6' 'UZC7' 'UZC8' 'UZC9'
  71. 'UXF0' 'UXF1' 'UXF2' 'UXF3' 'UXF4' 'UXF5' 'UXF6' 'UXF7' 'UXF8' 'UXF9'
  72. 'UYF0' 'UYF1' 'UYF2' 'UYF3' 'UYF4' 'UYF5' 'UYF6' 'UYF7' 'UYF8' 'UYF9'
  73. 'UZF0' 'UZF1' 'UZF2' 'UZF3' 'UZF4' 'UZF5' 'UZF6' 'UZF7' 'UZF8' 'UZF9' ;
  74. MNDUAL_HHO = 'MOTS'
  75. * 'FXC0' 'FXC1' 'FXC2' 'FXC3' 'FXC4' 'FXC5' 'FXC6' 'FXC7' 'FXC8' 'FXC9'
  76. * 'FYC0' 'FYC1' 'FYC2' 'FYC3' 'FYC4' 'FYC5' 'FYC6' 'FYC7' 'FYC8' 'FYC9'
  77. * 'FZC0' 'FZC1' 'FZC2' 'FZC3' 'FZC4' 'FZC5' 'FZC6' 'FZC7' 'FZC8' 'FZC9'
  78. 'FXF0' 'FXF1' 'FXF2' 'FXF3' 'FXF4' 'FXF5' 'FXF6' 'FXF7' 'FXF8' 'FXF9'
  79. 'FYF0' 'FYF1' 'FYF2' 'FYF3' 'FYF4' 'FYF5' 'FYF6' 'FYF7' 'FYF8' 'FYF9'
  80. 'FZF0' 'FZF1' 'FZF2' 'FZF3' 'FZF4' 'FZF5' 'FZF6' 'FZF7' 'FZF8' 'FZF9' ;
  81. * Pour l'istant on ne considere que les inconnues (DEPL et FORC) de face pour
  82. * le calcul de l'equilibre et l'acceleration de convergence
  83. * D'ou les lignes commentees ci-dessus dans la definition de MNPRIM_HHO et MNDUAL_HHO
  84. MNPRIM = MNPRIM 'ET' MNPRIM_HHO ;
  85. MLPRIM = MLPRIM 'ET' MNPRIM_HHO ;
  86. MNDUAL = MNDUAL 'ET' MNDUAL_HHO ;
  87. MLDUAL = MLDUAL 'ET' MNDUAL_HHO ;
  88. MLDEPL = MLDEPL 'ET' MNPRIM_HHO ;
  89. *A Ameliorer : ne faire que si presence HHO et n'ajouter que les DDLs utilises par HHO !
  90. *=FIN===== Formulation HHO =====================================================
  91. *
  92. * Pour stocker des informations necessaires aux calculs d'usure
  93. WTAB.'POST_COFR' = 'TABL' ;
  94. *----------------------------------------------------------------------*
  95. * Options de pasapas - parametres du calcul *
  96. *----------------------------------------------------------------------*
  97. * 1- rigidite
  98. AUTAUG = WTAB.'AUTOAUGM';
  99. IKSIA = WTAB.'K_SIGMA';
  100. IKTAN = WTAB.'K_TANGENT' ;
  101. ZKTASYM = 'TEXTE' ' ' ;
  102. 'SI' WTAB.'K_TANGENT_SYME' ;
  103. ZKTASYM = 'MOT' 'SYME' ;
  104. 'FINS';
  105. IMPLP = WTAB.'LIAISON_PERSISTANTE';
  106. IRAUGLU = WTAB.'RAIDAUGM';
  107. IRAUG = IRAUGLU ;
  108. 'SI' IRAUG;
  109. RIG_AUG = WTAB.'RIGIDITE_AUGMENTEE';
  110. 'FINS';
  111. IRCON = WTAB.'RAIDCONST';
  112. 'SI' IRCON;
  113. RIG_CONS = WTAB.'RIGIDITE_CONSTANTE';
  114. MAI_CONS ='EXTR' RIG_CONS 'MAIL' ;
  115. 'FINS';
  116. * declenchement du recalcul de la matrice
  117. DLTAIT = WTAB.'DELTAITER' ;
  118. ITRCLC = 1 '*' DLTAIT ;
  119. 'SI' (ITRCLC < 20) ; ITRCLC = 20; 'FINSI';
  120. *
  121. * 2- Type de formulation
  122. IFEFP = WTAB.'FEFP_FORMULATION' ;
  123. LAG_TOT= WTAB.'LAG_TOT';
  124. ISSTE = WTAB.'SUBSTEPPING';
  125. IFEFPUL= WTAB.'UPDATE_LAGRANGIAN';
  126. LNLOC = WTAB.'NLOC';
  127. 'SI' (LNLOC 'ET' ('EGA' WTAB.'NON_LOCAL' 'HELM'));
  128. TAHELM = WTAB.'HELMHOLTZ' ;
  129. NHELM = TAHELM . 'N_VARI_NL' ;
  130. 'FINSI' ;
  131. *
  132. * 3- Type de materiau
  133. ICERAM = WTAB.'CERAMIQUE' ;
  134. IENDOM = WTAB.'ENDOMMAGEMENT';
  135. IPLAVI = WTAB.'IPLAVI';
  136. POR1 = WTAB.'POR1' ;
  137. IVIEXT = WTAB.'VISCO_EXTERNE';
  138. IVIDOM = WTAB.'VISCODOMMAGE';
  139. IVISCO = WTAB.'VISCOPLASTIQUE';
  140. *
  141. * 4- Critere de convergence/non-convergence
  142. ZMAXIT = WTAB.'MAXITERATION' ;
  143. NITMA = WTAB.'NITERINTER_MAX';
  144. ZPRECM = WTAB.'PRECFLEX' ;
  145. ZPREK = WTAB.'PRECISINTER' ;
  146. ZPREC = WTAB.'PRECISION' ;
  147. ZPRECHPP = ZPREC ;
  148. IFTOL = 'NEG' WTAB.'FTOL' 'INCONNU' ;
  149. 'SI' IFTOL ;
  150. ZFTOL = 'ABS' WTAB.'FTOL' ;
  151. 'FINS';
  152. IMTOL = 'NEG' WTAB.'MTOL' 'INCONNU';
  153. 'SI' IMTOL ;
  154. ZMTOL = 'ABS' WTAB.'MTOL' ;
  155. 'FINS';
  156. 'SI' ('OU' IVISCO IVIDOM IVIEXT); ZPREK = 5.E-7 ; 'FINS';
  157. 'SI' IENDOM; ZPREK = ZPREC ; 'FINS';
  158. NSOINCR = WTAB.'SOUS_INCREMENT' ;
  159. *
  160. * 5- Type de calcul
  161. IPILOT = WTAB.'AUTOMATIQUE';
  162. ISOL = WTAB.'CONSOLIDATION';
  163. IDYN = WTAB.'DYNAMIQUE';
  164. IGRD = WTAB.'GRANDS_DEPLACEMENTS';
  165. HYPDEF = 'MOT' WTAB.'HYPOTHESE_DEFORMATIONS';
  166. HPP_EPS = 'EGA' WTAB.'PREDICTEUR' 'HPP';
  167. *
  168. * 6- Chargement particulier
  169. LOGDEF = WTAB.'CHAR_DEFI';
  170. LOGPRE = WTAB.'CHAR_PRES' ;
  171. LOGPIL = WTAB.'CHAR_PILO' ;
  172. ITHER = WTAB.'ITHER' ;
  173. *
  174. * 7- Acceleration de convergence
  175. ZNACCE = VRAI;
  176. * Option a n'utiliser que par les utilisateurs avertis
  177. 'SI' ('EXIS' PRECED 'ACCELERATION') ;
  178. III = PRECED.'ACCELERATION' ;
  179. 'SI' ('EGA' ('TYPE' III) 'ENTIER') ;
  180. ZNACCE = III ;
  181. 'FINS' ;
  182. 'FINS' ;
  183. *
  184. * 8- Instant de debut et de fin de pas
  185. TEMPS0 = WTAB.'TEMPS0';
  186. TI = WTAB.'T_FINAL';
  187. *
  188. * 9- Blocages mecaniques (BLOM)
  189. 'SI' WTAB.'CHAR_BLOM' ;
  190. BLOM1 = 'TIRE' PRECED.'BLOCAGES_MECANIQUES' 'BLOM' TI ;
  191. WTAB.'BLOCAGES_MECANIQUES' = BLOM1 ;
  192. 'FINSI' ;
  193. ZCLIM0 = WTAB.'BLOCAGES_MECANIQUES';
  194. ZCLIM = ZCLIM0 ;
  195. *----------------------------------------------------------------------*
  196. * Verification de la compatibilite entre options *
  197. *----------------------------------------------------------------------*
  198. 'SI' (IKTAN 'ET' ('NON' IPLAVI)) ;
  199. 'MESS' 'IPLAVI faux : pas de matrice tangente ->'
  200. ' on utilise la rigidite elastique' ;
  201. IKTAN = FAUX ;
  202. 'FINS' ;
  203. IPERT = WTAB.'K_TANGENT_PERT' 'ET' ('NON' LNLOC) 'ET' IPLAVI ;
  204. * Matrice tangente : pas d'acceleration en cas de modele FEFP ou SSTE
  205. 'SI' IKTAN ;
  206. 'SI' (IFEFP 'OU' ISSTE) ; ZNACCE = FAUX; 'FINS';
  207. 'FINS' ;
  208. *
  209. * Configurations de reference et de debut de pas
  210. GEOREF0 = WTAB.'FOR0' ;
  211. GEOM1 = WTAB.'GE0_DEB' ;
  212. *----------------------------------------------------------------------*
  213. * Modele complet *
  214. * *
  215. * Notations utilisees : *
  216. * ZMODLI : modeles mecanique + poreux (non //) *
  217. * -> aucune parallelisation alors ZMODL = ZMODLI *
  218. * -> parallelisation comportement alors ZMODL = ZMODLI *
  219. * MODRELOC est // *
  220. * -> parallelisation automatique alors ZMODL est // *
  221. *----------------------------------------------------------------------*
  222. * PAS_MODL : mise a jour des indices de PRECED.WTABLE relatifs aux
  223. * modeles si PRECED.WTABLE.MODELE a ete modifie.
  224. PAS_MODL PRECED ;
  225. 'SI' WTAB.'FOR_MECA' ;
  226. PARALLEL = FAUX ; PARTLOCA = FAUX ;
  227. ZMODLI = WTAB.'MO_TOT' ;
  228. ZMODL = ZMODLI ;
  229. 'SI' ('EGA' WTAB.'PROCESSEURS' 'COMPORTEMENT') ;
  230. PARALLEL = VRAI ; PARTLOCA = VRAI ;
  231. MODRELOC = 'PART' 'ARLE' ZMODLI WTAB.'NBPART' ;
  232. 'FINS' ;
  233. 'SI' ('EGA' WTAB.'PROCESSEURS' 'AUTOMATIQUE') ;
  234. PARALLEL = VRAI ; PARTLOCA = FAUX ;
  235. ZMODL = 'PART' 'ARLE' ZMODLI WTAB.'NBPART' ;
  236. 'OPTI' 'PARA' VRAI ;
  237. 'FINS' ;
  238. *
  239. * Si le modele a change, adapter egalement les MCHAMLS
  240. 'SI' ('NEG' ZMODLI WTAB.'MO_TOT_PREC') ;
  241. *
  242. CONTI.'CONTRAINTES' =('REDU' CONTI. 'CONTRAINTES' ZMODLI) '+' ('ZERO' ZMODLI 'CONTRAINTES') ;
  243. CONTI.'DEFORMATIONS'=('REDU' CONTI.'DEFORMATIONS' ZMODLI) '+' ('ZERO' ZMODLI 'DEFORMAT' ) ;
  244. 'SI' IPLAVI ;
  245. CONTI.'DEFORMATIONS_INELASTIQUES'=('REDU' CONTI.'DEFORMATIONS_INELASTIQUES' ZMODLI) '+' ('ZERO' ZMODLI 'DEFINELA') ;
  246. CONTI.'VARIABLES_INTERNES' =('REDU' CONTI.'VARIABLES_INTERNES' ZMODLI) '+' ('ZERO' ZMODLI 'VARINTER') ;
  247. 'FINS';
  248. *
  249. WTAB.'ETAT1' = PAS_ETAT PRECED TEMPS0 ;
  250. WTAB.'MAT1' = PAS_MATE PRECED WTAB.'ETAT1';
  251. *
  252. 'SI' ITHER ;
  253. WTAB.'ETHER1' =('REDU' WTAB.'ETHER1' ZMODLI) '+' ('ZERO' ZMODLI 'DEFORMAT') ;
  254. 'SI' POR1;
  255. WTAB.'MSRTHER1'=('REDU' WTAB.'MSRTHER1' ZMODLI) '+' ('ZERO' ZMODLI 'DEFORMAT') ;
  256. 'FINS';
  257. 'FINS';
  258. *
  259. 'SI' LOGDEF ;
  260. WTAB.'DEFOR1' = ('REDU' WTAB.'DEFOR1' ZMODLI) '+' ('ZERO' ZMODLI 'DEFORMAT') ;
  261. WTAB.'DEFOR2' = ('REDU' WTAB.'DEFOR2' ZMODLI) '+' ('ZERO' ZMODLI 'DEFORMAT') ;
  262. 'FINS';
  263. *
  264. 'FINS';
  265. 'FINSI';
  266. *----------------------------------------------------------------------*
  267. * Champ materiau et caracteristiques a l'instant TEMPS0 *
  268. * *
  269. * Notations utilisees : *
  270. * ZMAT11 : materiau a l'instant TEMPS0 (non //) *
  271. * -> aucune parallelisation alors ZMAT1 = ZMAT11 *
  272. * -> parallelisation comportement alors ZMAT1 = ZMAT11 *
  273. * -> parallelisation automatique alors ZMAT1 est // *
  274. *----------------------------------------------------------------------*
  275. ZMAT11 = WTAB.'MAT1' ;
  276. ZMAT1 = ZMAT11;
  277. 'SI' WTAB.'FOR_MECA' ;
  278. 'SI' ('EGA' WTAB.'PROCESSEURS' 'AUTOMATIQUE') ;
  279. ZMAT1 = 'REDU' ZMAT11 ZMODL ;
  280. 'FINS';
  281. *
  282. * Distinguer les modeles poreux des modeles mecaniques
  283. 'SI' POR1;
  284. MO_PORI = 'EXTR' ZMODLI 'FORM' 'POREUX';
  285. MO_POR = 'EXTR' ZMODL 'FORM' 'POREUX';
  286. *
  287. ** kich ma_por0 intervient si ISOL
  288. ** initialement MA_POR0= 'REDU' MO_POR MAT1 ;
  289. MA_POR0 = 'REDU' ZMAT11 MO_POR ;
  290. 'FINS';
  291. 'FINSI';
  292. *----------------------------------------------------------------------*
  293. * Champ materiau et caracteristiques a l'instant TI *
  294. * *
  295. * Notations utilisees : *
  296. * ZMAT22 : materiau a l'instant TI (non //) *
  297. * -> aucune parallelisation alors ZMAT2 = ZMAT22 *
  298. * -> parallelisation comportement alors ZMAT2 = ZMAT22 *
  299. * -> parallelisation automatique alors ZMAT2 est // *
  300. *----------------------------------------------------------------------*
  301. ZETAT2 = PAS_ETAT PRECED TI ;
  302. ZMAT2I = PAS_MATE PRECED ZETAT2;
  303. *
  304. 'SI' WTAB.'FOR_MECA' ;
  305. ZMAT22 = 'REDU' ZMAT2I ZMODLI ;
  306. 'SI' WTAB.'ITCAR';
  307. 'FORM' GEOREF0 ;
  308. GEOBID ZMAT22 = 'FORM' CONTI.'DEPLACEMENTS' ZMODLI ZMAT22;
  309. 'FORM' GEOM1 ;
  310. 'FINS';
  311. ZMAT2 = 'REDU' ZMAT22 ZMODL ;
  312. 'SI' POR1;
  313. MA_POR = 'REDU' ZMAT2I MO_POR ;
  314. 'FINS';
  315. 'FINSI' ;
  316. *----------------------------------------------------------------------*
  317. * Chargement a imposer a l'instant TI (F^ext_n+1) *
  318. *----------------------------------------------------------------------*
  319. ZFEXT2= 'VIDE' 'CHPOINT' / 'DISCRET';
  320. 'SI' ('EXIS' WTAB.'CHARGEMENT' 'MECA');
  321. ZFEXT2 = 'TIRE' WTAB.'CHARGEMENT' 'MECA' TI;
  322. TYP_2 = 'TYPE' ZFEXT2;
  323. 'SI' ('NEG' TYP_2 'CHPOINT ');
  324. 'MESS' '*** Erreur dans la definition du chargement ( MECA ) ***';
  325. 'ERRE' '*** Le type du champ n est pas CHPOINT ***';
  326. 'FINS';
  327. 'FINS';
  328. *
  329. 'SI' ('EXIS' WTAB.'CHARGEMENT' 'FORC');
  330. F2_FOR = 'TIRE' WTAB.'CHARGEMENT' 'FORC' TI;
  331. TYP_2 = 'TYPE' F2_FOR;
  332. 'SI' ('NEG' TYP_2 'CHPOINT ');
  333. 'MESS' '*** Erreur dans la definition du chargement ( FORC ) ***';
  334. 'ERRE' '*** Le type du champ n est pas CHPOINT ***';
  335. 'FINS';
  336. ZFEXT2 = ZFEXT2 '+' F2_FOR ;
  337. 'FINS';
  338. *
  339. 'SI' (LOGPRE 'ET' ('NON' IGRD)) ;
  340. MOP = WTAB.'MOD_PRE' ;
  341. ZPEXT ='TIRE' WTAB.'CHARGEMENT' 'PRES' TI ;
  342. 'SI' ('EXIS' WTAB 'MAT_PRE') ;
  343. ZFPEXT = 'BSIG' MOP ZPEXT WTAB.'MAT_PRE' ;
  344. 'SINON' ;
  345. ZFPEXT = 'BSIG' MOP ZPEXT ;
  346. 'FINS' ;
  347. ZFEXT2 = ZFEXT2 '+' ZFPEXT ;
  348. 'FINS' ;
  349. *
  350. 'SI' ('EXIS' WTAB.'CHARGEMENT' 'DIMP');
  351. F2_DEP = 'TIRE' WTAB.'CHARGEMENT' 'DIMP' TI;
  352. ZFEXT2 = ZFEXT2 + F2_DEP;
  353. 'FINS';
  354. *
  355. 'SI' ('EXIS' WTAB.'CHARGEMENT' 'DINC');
  356. F2_INC = 'TIRE' WTAB.'CHARGEMENT' 'DINC' TI;
  357. F2_mail = extr f2_inc maillage;
  358. zclim_mail = zclim extrai 'MAILLAGE';
  359. F2_Mcmp = zclim_mail elem appu larg f2_mail;
  360. zclim_inc = zclim redu f2_mcmp;
  361. f2_base = (zclim_inc * conti.'DEPLACEMENTS') exco 'FLX' 'NOID' 'FLX';
  362. ZFEXT2 = ZFEXT2 + F2_INC + F2_base; ;
  363. 'FINS';
  364. *
  365. 'SI' LOGDEF ;
  366. WTAB.'DEFOR2' = 'TIRE' WTAB.'CHARGEMENT' 'DEFI' TI;
  367. 'FINS';
  368. *
  369. 'SI' WTAB.'FOR_MECA' ;
  370. 'SI' ITHER ;
  371. 'SI' WTAB.'CHAR_THE';
  372. WTAB.'TET2' = 'TIRE' WTAB.'CHARGEMENT' 'T ' TI;
  373. 'FINSI';
  374. 'SI' WTAB.'FOR_THER';
  375. WTAB.'TET2' = estim.'TEMPERATURES' ;
  376. 'FINSI';
  377.  
  378. * Champs de temperature en debut (ZTEMP1) et fin de pas (ZTEMP2)
  379. * Si les maillages mecanique/thermique sont differents,
  380. * on projette le champ de temperature sur le modele mecanique
  381. 'SI' WTAB.'PROJECTION';
  382. 'FORM' GEOREF0 ;
  383. 'SI' ('EGA' ('TYPE' WTAB.'TET1') 'CHPOINT') ;
  384. ZTEMP1 = 'PROI' ('EXTR' ZMODL 'MAILLAGE') ('CHAN' 'CHAM' WTAB.'TET1' WTAB.'MOD_THE' 'NOEUD') ;
  385. ZTEMP2 = 'PROI' ('EXTR' ZMODL 'MAILLAGE') ('CHAN' 'CHAM' WTAB.'TET2' WTAB.'MOD_THE' 'NOEUD') ;
  386. 'SINON' ;
  387. ZTEMP1 = 'PROI' ('EXTR' ZMODL 'MAILLAGE') WTAB.'TET1' ;
  388. ZTEMP2 = 'PROI' ('EXTR' ZMODL 'MAILLAGE') WTAB.'TET2' ;
  389. 'FINSI' ;
  390. 'FORM' GEOM1 ;
  391. * Sinon, on vient seulement recuperer les champs de temperature
  392. 'SINON' ;
  393. ZTEMP1 = WTAB.'TET1';
  394. ZTEMP2 = WTAB.'TET2';
  395. 'FINSI' ;
  396. ETT0 ='REDU' WTAB.'ETHER1' ZMODL; 'COMM' 'REDU au cas ou le modele soit parallele';
  397. 'SI' POR1;
  398. MSRTT0='REDU' WTAB.'MSRTHER1' ZMODL; 'COMM' 'REDU au cas ou le modele soit parallele';
  399. 'FINS';
  400. 'SINON' ;
  401. che3 = 'MANU' 'CHML' ZMODL 'T' 20. 'STRESSES' ;
  402. che4 = 'MANU' 'CHML' ZMODL 'T' 20. 'STRESSES' ;
  403. 'FINSI';
  404. 'FINSI';
  405. *----------------------------------------------------------------------*
  406. * Dynamique : ajout au second membre de *
  407. * F0 + 4/DT*M*V0 - B0*SIG0 *
  408. *----------------------------------------------------------------------*
  409. * Matrice de masse
  410. 'SI' (IDYN 'OU' WTAB.'FREQUENTIEL') ;
  411. ICALMAS = FAUX ;
  412. 'SI' ('NON' ('EXIS' WTAB 'MASSE'));
  413. ICALMAS = VRAI ;
  414. 'FINSI' ;
  415. 'SI' (IGRD 'ET' IDYN) ;
  416. ICALMAS = VRAI ;
  417. 'FINSI' ;
  418. 'SI' ICALMAS ;
  419. WTAB.'MASSE' = 'MASS' WTAB.'MO_TOT' ('CHAN' 'MASSE'
  420. ZMAT22 WTAB.'MO_TOT' ) ;
  421. 'SI' (IDYN 'ET' WTAB.'MASSCONST') ;
  422. WTAB.'MASSE' = WTAB.'MASSE' 'ET' WTAB.'MASSE_CONSTANTE';
  423. 'FINS';
  424. 'FINS';
  425. 'FINS';
  426. *
  427. 'SI' IDYN ;
  428. UNSURH = 1.D0 '/' WTAB.'DT' ;
  429. 'SI' ('EGA' WTAB.'FREA1' 'INCONNU');
  430. F1 = 'VIDE' 'CHPOINT' / 'DISCRET';
  431. *
  432. * Forces exterieures a TEMPS0
  433. 'SI' ('EXIS' WTAB.'CHARGEMENT' 'MECA');
  434. F1M = 'TIRE' WTAB.'CHARGEMENT' 'MECA' WTAB.'TEMPS0' ;
  435. F1 = F1 '+' F1M ;
  436. 'FINS';
  437. *
  438. 'SI' ('EXIS' WTAB.'CHARGEMENT' 'FORC');
  439. F1F = 'TIRE' WTAB.'CHARGEMENT' 'FORC' WTAB.'TEMPS0' ;
  440. F1 = F1 '+' F1F ;
  441. 'FINS';
  442. *
  443. 'SI' WTAB.'PROCEDURE_CHARMECA';
  444. TFF1 = CHARMECA PRECED TEMPS0 ;
  445. 'SI' ('EXIS' TFF1 'ADDI_SECOND');
  446. FF1=TFF1 .'ADDI_SECOND';
  447. 'FINS';
  448. F1 = F1 '+' FF1;
  449. 'FINS';
  450. *
  451. 'SI' LOGPRE ;
  452. MOP = WTAB.'MOD_PRE' ;
  453. ZPEXT0 = 'TIRE' WTAB.'CHARGEMENT' 'PRES' TEMPS0 ;
  454. 'SI' ('EXIS' WTAB 'MAT_PRE') ;
  455. FF1 = 'BSIG' MOP ZPEXT0 ('REDU' WTAB.'MAT1' MOP) ;
  456. 'SINON' ;
  457. FF1 = 'BSIG' MOP ZPEXT0 ;
  458. 'FINS' ;
  459. F1 = F1 '+' FF1;
  460. 'FINS';
  461. *
  462. * Forces interieures a TEMPS0
  463. *HHO : Modifications appel a BSIGMA (ajout du champ de deplacements necessaire en HHO)
  464. LAF0 = 'BSIG' WTAB.'MO_TOT' conti.'CONTRAINTES' ZMAT22 conti.'DEPLACEMENTS' ;
  465. 'SI' IRCON ;
  466. LAF0 = LAF0 'ET'
  467. ('REDU' MAI_CONS (RIG_CONS '*' conti.'DEPLACEMENTS'));
  468. 'FINS';
  469. *
  470. * FREA1 : (masse*acceleration initiale)+(amortissement*vitesse initiale)
  471. WTAB.'FREA1' = F1 '-' LAF0 ;
  472. 'FINS';
  473. *
  474. 'SI' IMPLP ;
  475. * forces d'acceleration au debut du pas
  476. 'SI' ('NEG' WTAB.'AMORTISSEMENT' 'INCONNU');
  477. FF4 = WTAB.'AMORTISSEMENT' '*' conti.'VITESSES';
  478. WTAB.'FMAN' = WTAB.'FREA1' - FF4 ;
  479. 'SINON' ;
  480. WTAB.'FMAN' = WTAB.'FREA1' ;
  481. 'FINS';
  482. 'FINS';
  483. *
  484. FF = WTAB.'MASSE' '*' conti.'VITESSES';
  485. FF4 = 4. '/' WTAB.'DT' '*' FF; 'DETR' FF;
  486. *
  487. * partie du second membre qui ne depend que des informations du pas prec
  488. WTAB.'FREA1' = FF4 '+' WTAB.'FREA1';
  489. 'FINS';
  490. *
  491. * Advection
  492. 'SI' WTAB.'NVSTNL' ;
  493. 'SI' ('EXISTE' WTAB.'CHARGEMENT' 'ADVE') ;
  494. UADV = 'TIRE' WTAB.'CHARGEMENT' 'ADVE' TI ;
  495. LOGADV = VRAI ;
  496. 'SINON' ;
  497. LOGADV = FAUX ;
  498. 'FINSI' ;
  499. 'FINSI' ;
  500. *----------------------------------------------------------------------*
  501. * Consolidation : ajout au second membre de *
  502. * -B0*SIG0 + DT*(1-TETA)*FI0 + DT*H*P *
  503. *----------------------------------------------------------------------*
  504. 'SI' ISOL ;
  505. FF = 'BSIG' WTAB.'MOD_POR' conti.'CONTRAINTES' ;
  506. FF4 = 'EXCO' WTAB.'MOT_POR' FF
  507. WTAB.'MOT_POR' 'NOID' 'NATURE' 'DISCRET' ;
  508. ZFEXT2 = ZFEXT2 '+' FF4 ; 'DETR' FF ;
  509. *
  510. * ---- traitement des flux si besoin ----
  511. 'SI' ( 'EXIS' WTAB.'CHARGEMENT' 'FLUX' ) ;
  512. FLUXT0= 'TIRE' WTAB.'CHARGEMENT' 'FLUX' TEMPS0;
  513. FLUXTI= 'TIRE' WTAB.'CHARGEMENT' 'FLUX' TI ;
  514. FACFLU = -1. * WTAB.'DT';
  515. FLUXT = ( FACFLU * (1 - WTAB.'TETA') * FLUXT0 )
  516. + ( FACFLU * WTAB.'TETA' *FLUXTI ) ;
  517. ZFEXT2 = ZFEXT2 '+' FLUXT ;
  518. 'DETR' FLUXT ; 'DETR' FLUXT0; 'DETR' FLUXTI;
  519. 'FINS' ;
  520. 'FINS';
  521. *----------------------------------------------------------------------*
  522. * Pilotage indirect : preparation *
  523. *----------------------------------------------------------------------*
  524. 'SI' LOGPIL ;
  525. ZFPILIN = WTAB.'FORCES_PILOTEES' 'ET' WTAB.'DEPLACEMENTS_PILOTES' ;
  526. DTAU = ('IPOL' WTAB.'PARAMETRE_DE_PILOTAGE' TI) -
  527. ('IPOL' WTAB.'PARAMETRE_DE_PILOTAGE' TEMPS0) ;
  528. ETA0 = 'EXTR' PRECED.'COEFFICIENT_DE_PILOTAGE'
  529. ('DIME' PRECED.'COEFFICIENT_DE_PILOTAGE') ;
  530. D_ETA=0.;
  531. 'FINS' ;
  532. *----------------------------------------------------------------------*
  533. * Partie constante du second membre sur le pas de temps *
  534. *----------------------------------------------------------------------*
  535. ZFCONST1 = ZFEXT2 ;
  536. 'SI' IDYN ;
  537. ZDYFEXT = 'ENLE' ZFEXT2 'FLX';
  538. ZFEXT2 = ZFEXT2 '+' WTAB.'FREA1';
  539. 'FINSI' ;
  540. ZFLXB = 'EXCO' ZFEXT2 'FLX' 'NOID' 'FLX' ;
  541. *----------------------------------------------------------------------*
  542. * Recuperation de valeurs de WTAB pour initialisations eventuelles *
  543. *----------------------------------------------------------------------*
  544. ZRIBLO_M = 'MOT' 'INCONNU' ;
  545. 'SI' WTAB.'FOR_MECA' ;
  546. *
  547. ZXDENO = WTAB.'XDENO' ;
  548. ZXDENOM = WTAB.'XDENOM' ;
  549. * XDENOo = WTAB.'XDENO';
  550. * XDENOMo = WTAB.'XDENOM';
  551. ZETAT1 = WTAB.'ETAT1';
  552. 'SI' LOGDEF ;
  553. ZDEFOR1 = WTAB.'DEFOR1' ;
  554. ZDEFOR2 = WTAB.'DEFOR2' ;
  555. 'FINS';
  556. 'SI' ITHER ;
  557. ZTET1 = WTAB.'TET1' ;
  558. ZTET2 = WTAB.'TET2' ;
  559. 'FINS';
  560. 'SI' IPILOT ;
  561. ZAUTOCOEF = WTAB.'AUTOCOEF' ;
  562. ZAUTOREDU = WTAB.'AUTOREDU' ;
  563. ZAUTORED1 = WTAB.'AUTORED1' ;
  564. 'FINS';
  565. ZLASTKTAN = WTAB.'LASTKTAN' ;
  566. *
  567. 'SI' ('NEG' WTAB.'RIBLO_M' 'INCONNU') ;
  568. ZRIBLO_M = WTAB.'RIBLO_M' ;
  569. ZLISEA_M = WTAB.'LISEA_M' ;
  570. 'FINS' ;
  571. *
  572. ZFNONL = WTAB.'FNONL' ;
  573. 'SI' ('NEG' WTAB.'INCREMENT' 'INCONNU');
  574. INCRPREC = WTAB.'INCREMENT';
  575. 'FINS';
  576. *
  577. * teste t'on les moments ?
  578. TSTMOM = 'EGA' ('VALE' 'MODE') 'PLANGENE' ;
  579. TSTMOM = TSTMOM 'OU' ('EXIS' ZMAT22 'EPAI') 'OU' ('EXIS' ZMAT22 'INRY')
  580. 'OU' ('EXIS' ZMAT22 'MODS') 'OU'
  581. (('EXIS' ZMODLI 'ELEM' 'JOI1') 'ET' ('EXIS' ZMAT22 'V1X ')
  582. 'ET' ('EXIS' ZMAT22 'V1Y '));
  583. * teste t'on les POREUX ?
  584. 'SI' POR1 ; TSTMOM=VRAI ; 'FINS';
  585. *
  586. IKLFFF=VRAI;
  587. 'SI'TSTMOM;
  588. 'SI' IFTOL; 'SI' IMTOL;
  589. IKLFFF=FAUX;
  590. 'FINS'; 'FINS';
  591. 'SINON' ;
  592. 'SI' IFTOL; IKLFFF=FAUX; 'FINS';
  593. 'FINS';
  594. *
  595. * Jacobien du modele pour ponderer les champs de act3
  596. CH1 = 'MANU' 'CHML' zmodli 'SCAL' 1 'GRAVITE';
  597. zjac = 'INTG' 'ELEM' CH1 zmodli;
  598. zjac = zjac zmodli chan 'CHPO' 'SOMM' + xpetit;
  599. zjac = zjac **-0.5;
  600. *
  601. * Initialisation CHArgement SANS T :
  602. CHASANST= WTAB.'CHARGEMENT' 'ENLE' 'T' ;
  603. *
  604. * isoucomp : Indicateur SOUci dans COMPortement, FAUX par defaut
  605. isoucomp = faux ;
  606. PASUNIL = FAUX;
  607. residc = 'VIDE' 'CHPOINT' / 'DISCRET';
  608. RFNS = 'VIDE' 'RIGIDITE' / 'RIGIDITE' ;
  609. RTRSF = 'VIDE' 'RIGIDITE' / 'RIGIDITE' ;
  610. ZNITE = 0 ;
  611. ISOUSPPP = WTAB.'ISOUSPAS' ;
  612. WTAB.'ISOUSPAS' = 0;
  613. 'FINSI' ;
  614. *
  615. RED_URG = 0;
  616. 'SI' autaug; iraug = faux; 'FINSI';
  617. resmul = 1;
  618. augmult = 0.30000000;
  619. XCONV = 0.;
  620. DEPSTDM = 0.;
  621. WTAB.'CONV' = VRAI;
  622. KNOCONV = 0 ;
  623. *----------------------------------------------------------------------*
  624. * Valeurs des champs en debut de pas TEMPS0 *
  625. *----------------------------------------------------------------------*
  626. 'SI' WTAB.'FOR_MECA' ;
  627. ZDEP0 = CONTI.'DEPLACEMENTS' ;
  628. 'SI' (LAG_TOT 'EGA' 1); ZDEP0N = ZDEP0 * -1; 'FINSI';
  629. ZDEF0 = 'REDU' CONTI.'DEFORMATIONS' ZMODL ;
  630. ZSIG0 = 'REDU' CONTI.'CONTRAINTES' ZMODL ;
  631. com_sig = 'EXTR' ZMODLI 'CONT';
  632. 'SI' IPLAVI ;
  633. ZDEI0 = 'REDU' CONTI.'DEFORMATIONS_INELASTIQUES' ZMODL ;
  634. ZVAR0 = 'REDU' CONTI.'VARIABLES_INTERNES' ZMODL ;
  635. 'SI' ('NEG' WTAB.'MOVA' 'RIEN') ;
  636. ACC0 = 'EXCO' ZVAR0 WTAB.'MOVA' ;
  637. 'FINS' ;
  638. com_var = 'EXTR' ZMODLI 'VARI' ;
  639. com_dei = 'EXTR' ZMODLI 'DEIN' ;
  640. 'SI' ISOL ;
  641. com_maa = 'EXTR' ZMODLI 'MATE' ;
  642. 'FINS' ;
  643. 'FINS' ;
  644. ZGRDU0 = 'MOT' 'INCONNU' ;
  645. 'SI' ('NEG' WTAB.'MOD_MEC' 'INCONNU') ;
  646. 'SI' (IGRD 'ET' WTAB.'UTILISATEUR');
  647. 'FORM' GEOREF0 ;
  648. ZGRDU0 ='GRAD' ZMODL ZMAT2 ZDEP0 ;
  649. 'FORM' GEOM1 ;
  650. 'FINS' ;
  651. 'FINS' ;
  652. 'FINSI' ;
  653. dpsmax = 0;
  654. *######################################################################*
  655. *----------------------------------------------------------------------*
  656. * Boucle de non convergence BONOCONV (Configuration GEOM1) *
  657. *----------------------------------------------------------------------*
  658. 'REPETER' BONOCONV WTAB.'MAXSOUSPAS' ;
  659. **assi 2 trac cach v1 nclk;
  660. *
  661. * Initialisation objectif de fin de pas. Forces et deplacements
  662. ZFCONSTA = ZFEXT2 ;
  663. 'SI' IGRD ;
  664. 'FORM' GEOM1;
  665. 'FINSI' ;
  666. *
  667. excfconv = 0;
  668. 'SI' ('EGA' RED_URG 0);
  669. 'SI' ((XCONV < ZPREC) 'ET' (DEPSTDM < ZPREC) 'ET' (resmul > 0.99) 'ET' (AUGMULT < 100.));
  670. augmult = 0.30000000;
  671. *pv resmul = 1;
  672. augauto = augmult;
  673. augk = augmult;
  674.  
  675. 'SI' autaug ;
  676. RE_CNTRL = 'VIDE' 'RIGIDITE' / 'RIGIDITE';
  677. DE_CNTRL = FAUX;
  678. 'SI' (dpsmax < zprec);
  679. iraug = faux;
  680. 'FINSI';
  681. 'FINSI';
  682. znacce=VRAI;
  683. 'SI' (IFEFP 'OU' ISSTE) ; ZNACCE = FAUX; 'FINS';
  684. 'SINON';
  685. augmult = augmult / 1.8;
  686. si (augmult < 1d-3); augmult = 1d-3; finsi;
  687. 'FINSI';
  688. 'SINON';
  689. augmult = augmult * 1.5;
  690. 'SI' autaug; IRAUG = VRAI; 'FINSI';
  691. 'FINSI';
  692. *
  693. 'SI' WTAB.'FOR_MECA' ;
  694. KNOCONV = KNOCONV+1 ;
  695. *
  696. ZD0SLX = 'ENLE' ZDEP0 'LX';
  697. ZDEPT = ZDEP0 '-' ZD0SLX ;
  698. DEPST = 'ZERO' ZMODL 'DEFORMATIONS' ;
  699. ZSIGF = ZSIG0 ;
  700. *
  701. 'SI' ITHER ;
  702. DTETD = ZTEMP2 '-' ZTEMP1 ;
  703. 'FINS' ;
  704. *
  705. * --------------------------------------------------------------------
  706. * Reevaluer les materiaux en convergence forcee (TEMPS0 =TI)
  707. 'SI' (knoconv '>' 1) ;
  708. ZMAT11 = PAS_MATE PRECED ZETAT1 ;
  709. 'SI' WTAB.'ITCAR';
  710. 'FORM' GEOREF0 ;
  711. GEOBID ZMAT11 = 'FORM' ZDEP0 ZMODLI ZMAT11;
  712. 'FORM' GEOM1 ;
  713. 'FINS';
  714. ZMAT1 = 'REDU' ZMAT11 ZMODL ;
  715. *
  716. MCHC = PAS_ETAT PRECED TI ;
  717. ZMAT22 = PAS_MATE PRECED MCHC;
  718. 'SI' (IGRD 'ET' ('EGA' LAG_TOT 1)) ;
  719. ZMAT2I = ZMAT22 ;
  720. 'FINSI' ;
  721. 'SI' WTAB.'ITCAR';
  722. 'FORM' GEOREF0 ;
  723. GEOBID ZMAT22 = 'FORM' ZDEP0 ZMODLI ZMAT22;
  724. 'FORM' GEOM1 ;
  725. 'FINS';
  726. ZMAT2 = 'REDU' ZMAT22 ZMODL ;
  727. 'FINS' ;
  728. *
  729. * --------------------------------------------------------------------
  730. * Rigidite a la fin du pas
  731. 'SI' ('OU' WTAB.'RECARI' IRAUG ('NON' ('EXIS' WTAB 'RRRR')));
  732. *
  733. 'SI' ('OU' IENDOM IVIDOM ICERAM);
  734. HOOKENDO = 'HOOK' ZMODL ZMAT2 ZVAR0 ;
  735. RH = 'RIGI' ZMODL HOOKENDO ZMAT2 ;
  736. 'DETR' HOOKENDO;
  737. 'SINON';
  738. 'SI' (IGRD 'ET' ('EGA' LAG_TOT 1)) ;
  739. 'FORM' GEOREF0;
  740. HOOKRH = 'HOOK' ZMODL ('REDU' ZMAT2I ZMODL) ;
  741. HOOKRH2 = 'PICA' HYPDEF HOOKRH ZMODL ZDEP0;
  742. 'FORM' GEOM1;
  743. RH = 'RIGI' HOOKRH2 ZMODL ZMAT2;
  744. 'DETR' HOOKRH; 'DETR' HOOKRH2;
  745. 'SINON';
  746. RH = 'RIGI' ZMODL ZMAT2 ;
  747. 'FINSI';
  748. 'FINS';
  749. *
  750. * RH peut contenir des CL
  751. ZCL = 'EXTR' RH 'RIGI' 'MULT' ;
  752. 'SI' ('NEG' ('DIME' ZCL) 0) ;
  753. ZCLIM = ZCLIM0 'ET' ZCL ;
  754. 'FINSI' ;
  755. RRRR = RH 'ET' ZCLIM0 ;
  756. *
  757. * Option 'RIGIDITE_CONSTANTE'
  758. 'SI' IRCON;
  759. RRRR = RRRR 'ET' RIG_CONS;
  760. 'FINS';
  761. ZRAIDNA = RRRR;
  762. *
  763. * Option 'RIGIDITE_AUGMENTEE'
  764. 'SI' (IRAUG 'ET' ('NON' AUTAUG));
  765. RRRR = RRRR 'ET' RIG_AUG ;
  766. 'FINSI';
  767. *
  768. * Option 'AUGMENTATION_AUTOMATIQUE'
  769. 'SI' (AUTAUG 'ET' ('NON' IRAUGLU));
  770. RIG_AUG = 'MASSE' ZMODL ZMAT2 ;
  771. *** mess 'actualisation rig_aug';
  772. 'FINSI';
  773. *
  774. 'SI' (AUTAUG 'ET' IRAUG);
  775. * recalcul masse pour l'augmentation
  776. RRRR = RRRR 'ET' (RIG_AUG * augauto) 'ET' (RH * augk) ;
  777. 'FINS';
  778. *
  779. 'SI' (('NON' IRAUG) 'ET' IRAUGLU) ;
  780. *** 'MESS' ' augmentation residuelle ';
  781. RRRR = RRRR 'ET' RIG_AUG ;
  782. 'FINSI';
  783. *
  784. * Stockage de la rigidite pour eviter de la recalculer
  785. 'SI' ('EGA' ('DIME' ZCL) 0) ;
  786. WTAB.'RRRR'=RRRR;
  787. 'FINSI';
  788. *
  789. 'FINS';
  790. ZRAID = WTAB.'RRRR' ;
  791.  
  792. 'FINSI' ;
  793. *FOR_MECA
  794. *
  795. 'SI' WTAB.'NVSTNL' ;
  796. 'SI' ('NON' ('EXISTE' WTAB 'RNSL')) ;
  797. 'SI' (WTAB.'DT' '&lt;EG' 0.) ; 'MESS' 'PAS DE TEMPS !! ' WTAB.'DT' ;
  798. 'ERRE' 21 ; 'FINSI' ;
  799. WTAB.'MAT_NSL' = WTAB.'MAT_NSL' 'ET'
  800. ('MATE' WTAB.'MOD_NSL' 'FREQ' (WTAB.'DT'**(-1))) ;
  801. RINSL = 'RIGI' WTAB.'MOD_NSL' WTAB.'MAT_NSL' ;
  802. ndiv = 'RIGI' 'DIVE' WTAB.'MOD_NSL' WTAB.'MAT_NSL' ;
  803. WTAB.'NSMA' = 'MASSE' WTAB.'MOD_NSL' WTAB.'MAT_NSL' ;
  804. WTAB.'RNSL' = RINSL 'ET' ZCLIM0 'ET' ndiv ;
  805. 'FINSI' ;
  806. 'FINSI' ;
  807. *
  808. *------------ consolidation ou dynamique faut-il recalculer l'operateur?
  809. 'SI' (ISOL 'OU' IDYN);
  810. 'SI' ( '>' (WTAB.'DTPREC' '*' 0.9999) WTAB.'DT') ;
  811. WTAB.'RECAOP' = VRAI;
  812. 'FINS';
  813. 'SI' ( '<' (WTAB.'DTPREC' '*' 1.0001) WTAB.'DT') ;
  814. WTAB.'RECAOP' = VRAI;
  815. 'FINS';
  816. 'SI' WTAB.'MATVAR';
  817. WTAB.'RECAOP' = VRAI;
  818. 'FINS';
  819. 'SI' ('NON' WTAB.'RECAOP') ;
  820. 'SI' ('NEG' WTAB.'OPERATEUR' 'INCONNU');
  821. ZRAID = WTAB.'OPERATEUR';
  822. 'FINS';
  823. 'FINS';
  824. 'FINS';
  825. *
  826. *------------ operateur amortissement en frequentiel
  827. 'SI' WTAB.'FREQUENTIEL' ;
  828. RRR2 = 'AMOR' ZMODL ZMAT2 ;
  829. RR2 = 'CHAN' 'INCO' RRR2
  830. ('MOTS' 'ALFA' 'BETA') ('MOTS' 'IALF' 'IBET')
  831. ('MOTS' 'FALF' 'FBET') ('MOTS' 'FALF' 'FBET') 'QUEL' ;
  832. RR3 = 'CHAN' 'INCO' RRR2
  833. ('MOTS' 'ALFA' 'BETA') ('MOTS' 'ALFA' 'BETA')
  834. ('MOTS' 'FALF' 'FBET') ('MOTS' 'IFAL' 'IFBE') 'QUEL' ;
  835. RRR2 = RR2 'ET' RR3 ;
  836.  
  837. RR1 = ZRAID ;
  838. OMEGI= 2.* PI * TI ;
  839. RRR1 = OMEGI * OMEGI * (-1.) * WTAB.'MASSE' ;
  840. RR1 = ZRAID 'ET' RRR1 ;
  841. RR4 = 'CHAN' 'INCO' (RR1 '*' (-1.D0))
  842. ('MOTS' 'ALFA' 'BETA') ('MOTS' 'IALF' 'IBET')
  843. ('MOTS' 'FALF' 'FBET') ('MOTS' 'IFAL' 'IFBE') 'QUEL';
  844. ZRAID= RR1 'ET' RR4 ;
  845. RR5 = OMEGI '*' RRR2 ;
  846. ZRAID= ZRAID 'ET' RR5 ;
  847. 'FINS' ;
  848.  
  849. *--------------- et la perméabilité ----------------------------------
  850. 'SI' ISOL ;
  851. 'SI' (IGRD 'OU' WTAB.'MATVAR') ;
  852. WTAB.'PERMEABILITE' = 'PERM' WTAB.'MOD_POR' ZMAT2;
  853. WTAB.'RECAOP' = VRAI ;
  854. 'FINS';
  855. 'FINS';
  856.  
  857. *------------- Cas de la consolidation ou de la dynamique -------------
  858. *------------- il faut recalculer l'operateur d'iteration -------------
  859. 'SI' WTAB.'RECAOP' ;
  860. *
  861. 'SI' IDYN;
  862. ZRMAS = 4.D0 '/' (WTAB.'DT' '**' 2) '*' WTAB.'MASSE' ;
  863. ZRAID = ZRMAS 'ET' ZRAID;
  864. 'SI' ('NEG' WTAB.'AMORTISSEMENT' 'INCONNU');
  865. ZRAID = WTAB.'AMORTISSEMENT' '*' (2.D0 '/' WTAB.'DT') 'ET' ZRAID;
  866. 'FINS';
  867. 'FINS' ;
  868. *
  869. 'SI' ISOL ;
  870. ZRAID =-1.* WTAB.'DT'* WTAB.'TETA'* WTAB.'PERMEABILITE' 'ET' ZRAID ;
  871. 'FINS' ;
  872. *
  873. WTAB.'OPERATEUR'= ZRAID ;
  874. 'FINS';
  875. ZRAIDNA = ZRAID;
  876. *
  877. 'SI' WTAB.'FOR_MECA' ;
  878. *
  879. * --------------------------------------------------------------------
  880. * Conditions de contact-frottement
  881. 'SI' WTAB.'CONTACT';
  882. ITACC = 4;
  883. MODCON = WTAB.'MODCONTA';
  884. 'SI' ('EXIS' WTAB 'MATCONTA') ;
  885. WTAB.'MATCONTA' = 'REDU' MODCON WTAB.'MATCONTA';
  886. CJEU CRR = 'RFCO' MODCON WTAB.'CONV' WTAB.'MATCONTA' ;
  887. 'SINON' ;
  888. CJEU CRR = 'RFCO' MODCON WTAB.'CONV' ;
  889. 'FINS';
  890. *
  891. * Separer ce qui concerne l'adherence (FADH) et les jeux (FLX)
  892. 'SI' ('NEG' CJEU 0) ;
  893. FADHE = 'EXCO' CJEU 'FADH' 'NOID' 'FADH' ;
  894. CDAP = 'EXCO' CJEU 'FLX' 'NOID' 'FLX' ;
  895. 'FINSI' ;
  896. *
  897. 'SI' WTAB.'MODAL' ;
  898. CRR = 'PJBA' CRR ZMODL ZMAT2 ;
  899. MCRR = 'EXTR' CRR 'MAIL' 'MULT' ;
  900. MCDAP= 'EXTR' CDAP 'MAIL' ;
  901. 'REPETER' BCDA ('NBNO' MCDAP) ;
  902. PBCDA ='POIN' MCDAP &BCDA ;
  903. PCRR ='POIN' MCRR &BCDA ;
  904. CHCR ='MANU' 'CHPO' PCRR 1 'FLX' ('EXTR' CDAP 'FLX' PBCDA)
  905. 'NATURE' 'DISCRETE' ;
  906. 'SI' ('EGA' 1 &BCDA) ;
  907. CCDA = CHCR ;
  908. 'SINON' ;
  909. CCDA = CHCR 'ET' CCDA ;
  910. 'FINS' ;
  911. 'FIN' BCDA ;
  912. CDAP = CCDA ;
  913. 'FINS' ;
  914. *
  915. ZFCONT = 'VIDE' 'CHPOINT' / 'DISCRET' ;
  916. 'SI' ('NEG' CRR 0) ;
  917. ZFCONT = CRR '*' ZD0SLX 'EXCO' 'FLX' 'NOID' 'FLX';
  918. ZCLIM = CRR 'ET' ZCLIM ;
  919. ZRAID = CRR 'ET' ZRAID ;
  920. 'FINS';
  921. *
  922. 'SI' ('NEG' CJEU 0) ;
  923. ZFCONT = ZFCONT '+' CDAP ;
  924. 'FINS';
  925. *
  926. 'SI' (('NEG' CRR 0) 'OU' ('NEG' CJEU 0)) ;
  927. ZFCONSTA = ZFCONSTA '+' ZFCONT ;
  928. 'FINS';
  929. *
  930. 'FINS';
  931. *
  932. * --------------------------------------------------------------------
  933. * Matrice tangente et fefp
  934. 'SI' IKTAN ;
  935. 'SI' IFEFP ;
  936. IKT_SAUV = VRAI ;
  937. 'SI' ('NEG' WTAB.'LASTKTAN' 'INCONNU') ;
  938. 'MESS' 'FEFP: Start with LASTKTAN' ;
  939. ZRIKTA = ZLASTKTAN ;
  940. ZRAID = ZCLIM 'ET' ZRIKTA ;
  941. 'SINON' ;
  942. 'MESS' 'FEFP: Previous KTAN not available' ;
  943. ZRAID = ZRAID 'ET' ('KSIGMA' ZMODL ZSIG0 ZMAT2) ;
  944. 'FINS' ;
  945. 'SINON' ;
  946. 'SI' ('NEG' WTAB.'LASTKTAN' 'INCONNU') ;
  947. IKT_SAUV = VRAI ;
  948. 'SI' IPERT ;
  949. 'MESS' 'Matrice tangente par perturbation - '
  950. 'Demarrage avec KTAN = LASTKTAN' ;
  951. 'SINON' ;
  952. 'MESS' 'Matrice tangente "coherente" - '
  953. 'Demarrage avec KTAN = LASTKTAN' ;
  954. 'FINS' ;
  955. ZRIKTA = ZLASTKTAN ;
  956. ZRAID = ZCLIM 'ET' ZRIKTA ;
  957. 'SINON' ;
  958. 'SI' IPERT ;
  959. IKT_SAUV='NEG' WTAB.'K_TANGENT_ITER0'
  960. 'MAT_ELASTIQUE';
  961. 'MESS' 'Matrice tangente par perturbation - '
  962. 'Demarrage avec KTAN = rigidite elastique' ;
  963. 'SINON' ;
  964. IKT_SAUV = ('NEG' WTAB.'K_TANGENT_ITER0' 'MAT_ELASTIQUE')
  965. 'ET' ('NEG' WTAB.'K_TANGENT_ITER0' 'MAT_TANGENTE') ;
  966. 'SI' ('EGA' WTAB.'K_TANGENT_ITER0' 'MAT_ELASTIQUE') ;
  967. 'MESS' 'Matrice tangente "coherente" - '
  968. 'Demarrage avec KTAN = rigidite elastique' ;
  969. 'SINON' ;
  970. 'MESS' 'Matrice tangente "coherente" - '
  971. 'Demarrage avec KTAN (DTTAN = 0.)' ;
  972. DTTAN = 0. ;
  973. ZRIKTA = 'KTAN' ZMODL ZSIG0 ZVAR0 ZMAT2
  974. 'PREC' ZPREK 'DT ' DTTAN ZKTASYM ;
  975. ZRAID = ZCLIM 'ET' ZRIKTA ;
  976. 'FINS' ;
  977. 'FINS' ;
  978. 'FINS' ;
  979. 'FINS' ;
  980. 'FINS' ;
  981. *
  982. * --------------------------------------------------------------------
  983. * Raideur geometrique (option K_SIGMA)
  984. 'SI' IGRD ;
  985. 'SI' (IKSIA 'ET' ('NON' IFEFP) 'ET' ('NON' HPP_EPS)) ;
  986. KSIGTC = 'KSIGMA' ZMODL ZSIG0 ZMAT2;
  987. ZRAIDINI= ZRAID ;
  988. ZRAID = ZRAID 'ET' KSIGTC;
  989. 'FINS' ;
  990. 'FINS' ;
  991. *----------------------------------------------------------------------*
  992. * Pilotage automatique
  993. *----------------------------------------------------------------------*
  994. COEPI = 1.D0 ; CORPREC = 1. ;
  995. *
  996. 'SI' IPILOT ;
  997. AL1 = 1. ; COEPI0 = 1.d0;
  998. 'SI' WTAB.'AUTODEUX' ;
  999. COEPI = 'ABS' ZAUTOCOEF ;
  1000. COEPI = COEPI / (1.-COEPI);
  1001. 'SI' (COEPI > 1.D0) ; COEPI=1.D0;'FINS';
  1002. COEPI0=COEPI;ZAUTOCOEF =COEPI;
  1003. 'SINON';
  1004. ZAUTOCOEF = 1.D0;
  1005. 'FINS' ;
  1006. RED1 = 1. ;
  1007. RED2 = 0 ;
  1008. 'SI' ('NEG' ZAUTORED1 'INCONNU') ;
  1009. 'SI' (ZAUTORED1 > 0);
  1010. ZAUTORED1 = ZAUTORED1 - 1;
  1011. 'SI' (ZAUTORED1 'EGA' 0) ;
  1012. COEPI = 3 * COEPI ;
  1013. ZAUTOREDU = ZAUTOREDU / 3.;
  1014. 'SI' (ZAUTOREDU > 1.1 );
  1015. * on travaille encore avec un critere reduit
  1016. ZAUTORED1 = 4 ;
  1017. RED1 = 3. ;
  1018. 'FINS';
  1019. 'MESS' 'On multiplie le critere de pilotage par 3';
  1020. 'FINS';
  1021. 'FINS';
  1022. 'FINS';
  1023. 'SI' (COEPI > 1d0);
  1024. RED1 = RED1 / COEPI ;COEPI =1d0;
  1025. 'FINS';
  1026. 'SI' ('NEG' WTAB.'NBPLAS' 'INCONNU') ;
  1027. 'SI' (WTAB.'NBPLAS' 'EGA' 0) ;
  1028. 'SI' (COEPI < 0.) ;COEPI = COEPI * -2.;'FINS';
  1029. 'FINS';
  1030. 'FINS';
  1031. COEPI = 'ABS' COEPI ;
  1032. COEPI0 = COEPI;
  1033. 'SINON';
  1034. ZAUTOCOEF = 1.D0 ;
  1035. 'FINS';
  1036. *
  1037. 'FINSI' ;
  1038. *FOR_MECA
  1039. *----------------------------------------------------------------------*
  1040. * Partie variable du chargement a imposer a TI *
  1041. * (geometrie, materiaux, chargement en deformations) *
  1042. *----------------------------------------------------------------------*
  1043. XDENDEF = 0. ; FDENDEF = 0. ;
  1044. DEPST0 = 0. ; DMSRT0 = 0. ;
  1045. DFTHE = 0. ; DFDEF = 0. ;
  1046. ZFSUIV = 0. ;
  1047. *
  1048. * - Corrections en cas de materiaux variables
  1049. 'SI' WTAB.'MATVAR' ;
  1050. 'SI' ('OU' IENDOM IVIDOM ICERAM);
  1051. XXX3 ='ELAS' ZMODL ZSIG0 ZMAT2 ZVAR0;
  1052. XXX4 ='ELAS' ZMODL ZSIG0 ZMAT1 ZVAR0;
  1053. 'SINON';
  1054. XXX3 ='ELAS' ZMODL ZSIG0 ZMAT2 ;
  1055. XXX4 ='ELAS' ZMODL ZSIG0 ZMAT1 ;
  1056. 'FINS';
  1057. DDEF0 = XXX4 - XXX3;
  1058. 'DETR' XXX3; 'DETR' XXX4;
  1059. DEPST0 = -1.* DDEF0;
  1060. 'FINS';
  1061. *
  1062. * - En cas de chargement thermiques
  1063. 'SI' ITHER ;
  1064. 'SI' ('EGA' ('TYPE' ZTEMP2) 'CHPOINT');
  1065. MCHTETA2 ='CHAN' 'CHAM' ZTEMP2 ZMODL 'STRESSES' 'TEMPERATURES' ;
  1066. 'SINON';
  1067. MCHTETA2 = ZTEMP2;
  1068. 'FINS';
  1069. ETREF = 'REDU' WTAB.'ETREF' ZMODL ; 'COMM' 'REDU au cas ou le modele soit parallele';
  1070. ETT =('EPTH' ZMODL ZMAT2 MCHTETA2 ) - ETREF ;
  1071. DTT = ETT '-' ETT0 ;
  1072. XDENDEF = XDENDEF '+' ETT ;
  1073.  
  1074. 'SI' WTAB.'POR1';
  1075. *------ Cas du milieu poreux avec chargement thermique ----------
  1076. * cas isotrope seulement pour le moment
  1077. * et on ne s'occupe pas du alpha-reference !!
  1078. MSRTT = PAS_EPTH PRECED ZMODL ZMAT2 MCHTETA2 ;
  1079. DMSRT0= MSRTT '-' MSRTT0 ;
  1080. 'FINS' ;
  1081.  
  1082. DEPST0 = DTT '+' DEPST0 ;
  1083. 'FINS';
  1084. *
  1085. * - Contraintes et forces equivalentes associees a ces deformations
  1086. 'SI' (DEPST0 'NEG' 0);
  1087. 'SI' ('OU' IENDOM IVIDOM ICERAM);
  1088. DSIGT0 = 'ELAS' ZMODL DEPST0 ZMAT2 ZVAR0;
  1089. 'SI' (XDENDEF 'NEG' 0) ;
  1090. XDENDEF = 'ELAS' ZMODL XDENDEF ZMAT2 ZVAR0;
  1091. 'FINSI' ;
  1092. 'SINON';
  1093. DSIGT0 = 'ELAS' ZMODL DEPST0 ZMAT2 ;
  1094. 'SI' (XDENDEF 'NEG' 0) ;
  1095. XDENDEF = 'ELAS' ZMODL XDENDEF ZMAT2 ;
  1096. 'FINSI' ;
  1097. 'FINS' ;
  1098.  
  1099. 'SI' (DMSRT0 'NEG' 0);
  1100. DSIGT0 = DSIGT0 '+' DMSRT0 ;
  1101. XDENDEF = XDENDEF '+' MSRTT ;
  1102. 'FINS';
  1103. DFTHE ='BSIG' ZMODL DSIGT0 ZMAT2 ;
  1104. 'SI' (XDENDEF 'NEG' 0) ;
  1105. FDENDEF = 'BSIG' ZMODL XDENDEF ZMAT2 ;
  1106. 'FINSI' ;
  1107. 'FINS';
  1108. *
  1109. * - Chargement en deformations imposees
  1110. 'SI' LOGDEF;
  1111. DDEFOR0 = 'REDU' (ZDEFOR2 - ZDEFOR1) ZMODL;
  1112. XDENDEF = 'REDU' ZDEFOR2 ZMODL ;
  1113. 'SI' ('OU' IENDOM IVIDOM ICERAM);
  1114. DSI1 = 'ELAS' ZMODL DDEFOR0 ZMAT2 ZVAR0;
  1115. XDENDEF = 'ELAS' ZMODL XDENDEF ZMAT2 ZVAR0;
  1116. 'SINON';
  1117. DSI1 = 'ELAS' ZMODL DDEFOR0 ZMAT2 ;
  1118. XDENDEF = 'ELAS' ZMODL XDENDEF ZMAT2 ;
  1119. 'FINS';
  1120. DFDEF = 'BSIG' ZMODL DSI1 ZMAT2 ;
  1121. FDENDEF = FDENDEF '+' ('BSIG' ZMODL XDENDEF ZMAT2) ;
  1122. 'FINS';
  1123. *
  1124. * - Chargement de pression suiveuse en grands deplacements
  1125. 'SI' (LOGPRE 'ET' IGRD) ;
  1126. MOP = WTAB.'MOD_PRE' ;
  1127. ZPEXTF = 'TIRE' WTAB.'CHARGEMENT' 'PRES' TI ;
  1128. 'SI' ('EXIS' WTAB 'MAT_PRE') ;
  1129. ZFPEXTF = 'BSIG' MOP ZPEXTF WTAB.'MAT_PRE' ;
  1130. 'SINON' ;
  1131. ZFPEXTF = 'BSIG' MOP ZPEXTF ;
  1132. 'FINS' ;
  1133. ZFSUIV = ZFSUIV '+' ZFPEXTF ;
  1134. 'FINS';
  1135. *
  1136. * - Procedure UTILISATEUR : y a-t-il des forces non conservatives
  1137. 'SI' WTAB.'PROCEDURE_CHARMECA';
  1138. * on ajoute l indice ADDI_MATRICE pour signaler a charmeca qu on
  1139. * souhaite aussi l operateur linearisé des Forces NL de charmeca
  1140. PRECED.'ADDI_MATRICE' = vrai;
  1141. TFP22 = CHARMECA PRECED TI ;
  1142. PRECED.'ADDI_MATRICE' = faux;
  1143. *
  1144. * FP22 = F^suiv_n+1
  1145. ADDISEC2 = 'EXIS' TFP22 'ADDI_SECOND' ;
  1146. 'SI' ADDISEC2 ;
  1147. FP22 = TFP22.'ADDI_SECOND' ;
  1148. FP022 = 'COPIER' FP22 ;
  1149. ZFSUIV = ZFSUIV '+' FP22 ;
  1150. 'FINS';
  1151. 'SI' ('EXIS' TFP22 'ADDI_MATRICE');
  1152. ZRAID = ZRAID 'ET' TFP22.'ADDI_MATRICE';
  1153. 'FINS';
  1154. 'FINS';
  1155. *
  1156. * A-t-on des C.L. unilaterales ?
  1157. XXX1 = 'EXTR' ZCLIM 'MAIL' 'UNIL';
  1158. IMPO12 = ('NBEL' XXX1) '>' 0 ;
  1159. *----------------------------------------------------------------------*
  1160. * Second membre RESIDU *
  1161. * *
  1162. * Calcul du premier residu : desequilibre entre les forces externes et *
  1163. * le calcul B*SIGMA. Le sigma qui sert est celui qui existerait si *
  1164. * le champ de deplacement ne changeait pas (ZDEP0). *
  1165. * *
  1166. * | forces exterieures sans reactions - forces interieures | *
  1167. * | F^ext_n+1 + DF^suiv + DF^ther + DF^defi - F^int_n | *
  1168. * | |-------------------------| | *
  1169. * | -1.*ZFPLO | *
  1170. * | (composantes de forces FX FY FZ ...) | *
  1171. * RESIDU = | | *
  1172. * | increment des relations imposees Du^imp | *
  1173. * | (composantes de depl. FLX) | *
  1174. * *
  1175. * A F^ext peuvent s'ajouter des termes supplementaires p.ex. en *
  1176. * dynamique ou en poreux. *
  1177. *----------------------------------------------------------------------*
  1178. * ZFCONSTA = [ F^ext_n+1 ; u^imp_n+1 ]
  1179. ZFEXT = 'ENLE' ZFCONSTA 'FLX';
  1180. ZFLX1 = 'EXCO' ZFCONSTA 'FLX' 'NOID' 'FLX' 'NATURE' 'DISCRET';
  1181. *
  1182. 'SI' WTAB.'FOR_MECA' ;
  1183. *
  1184. * Forces externes deja equilibrees au debut du pas par B*sigma
  1185. * ZF1 = F^int_n = B*sigma_n + K^cst*u_n
  1186. *HHO : Modifications appel a BSIGMA (ajout du champ de deplacements necessaire en HHO)
  1187. ZF1 = 'BSIG' ZMODL ZSIG0 ZMAT1 ZDEP0 ;
  1188. 'SI' IRCON;
  1189. ZF1 = ZF1 '+' ('REDU' MAI_CONS (RIG_CONS '*' ZDEP0));
  1190. 'FINS';
  1191. *
  1192. 'SI' IDYN ;
  1193. FFDYN = 'COPIER' ZF1;
  1194. 'FINS';
  1195. *
  1196. 'SI' ISOL ;
  1197. GRAP0= 'GRAD' MO_POR ZDEP0 MA_POR0 'CONS' ;
  1198. XXX1 = 'GRAD' MO_POR ZDEP0 MA_POR 'CONS' ;
  1199. XXXS =((1.- WTAB. 'TETA' )*GRAP0)+ (WTAB. 'TETA' * XXX1);
  1200. XXX2 = WTAB.'DT' '*' ('GNFL' MO_POR XXXS) ;
  1201. XXX3 = ZF1 ;
  1202. ZF1 = XXX3 - XXX2;
  1203. 'DETR' XXX3; 'DETR' XXX2 ;
  1204. 'FINS';
  1205. *
  1206. * [ FLXINI ; FREAP ] = [ u^imp_n ; -F^reac_n ]
  1207. FZDEP0 = ZDEP0 '*' ZCLIM;
  1208. FLXINI = 'EXCO' FZDEP0 'FLX' 'NOID' 'FLX' 'NATURE' 'DISCRET';
  1209. FREAP = 'ENLE' FZDEP0 'FLX';
  1210. *
  1211. * XXX1 = [ F^ext_n+1 ; u^imp_n+1 - u^imp_n ]
  1212. XXX1 = ZFCONSTA '-' FLXINI ;
  1213. *
  1214. * FEXT0 : chargement externe (sans reactions) au pas precedent
  1215. * FEXT0 = F^int_n - F^reac_n = F^ext_n
  1216. FEXT0 = ZF1 '+' FREAP;
  1217. *
  1218. * RESIDU : forces exterieures sans reactions (avec des termes
  1219. * supplementaires le cas echeant p.ex. en dynamique ou en poreux)
  1220. * - forces interieures et increment des relations imposees
  1221. *
  1222. * RESIDU = [DF^tot ; Du^imp]
  1223. * RESIDU = [F^ext_n+1 + F^suiv_n+1 + DF^ther + DF^defi - F^int_n ; Du^imp]
  1224. RESIDU = XXX1 '+' DFTHE '+' DFDEF '-' ZF1;
  1225. RESIDU = RESIDU '+' ZFSUIV ;
  1226. *
  1227. * ZFPLO = -1.* [- F^int_n + DF^ther + DF^defi ]
  1228. ZFPLO = ZF1 '-' DFTHE '-' DFDEF ;
  1229. *
  1230. * ZDFINI = [F^ext_n+1 + DF^ther + DF^defi - (F^int_n - F^reac_n) ; Du^imp]
  1231. ZDFINI = XXX1 '-' FEXT0 ;
  1232. 'SI' IPILOT ;
  1233. ZLX = 'EXCO' ZDEP0 'LX' 'NOID' 'LX' 'NATURE' 'DIFFUS';
  1234. DFEXT0F = 'ENLE' ZDFINI 'FLX';
  1235. DFEXT0L = 'EXCO' ZDFINI 'FLX' 'NOID' 'FLX';
  1236. 'SI' IMPO12 ;
  1237. DUUNIL = 'REDU' DFEXT0L XXX1 ;
  1238. XXX2 = 'DIFF' ('EXTR' ZCLIM 'MAIL') XXX1 ;
  1239. DUIMPO = 'REDU' DFEXT0L XXX2 ;
  1240. 'FINS';
  1241. 'FINSI' ;
  1242. 'SI' (ITHER 'OU' WTAB.'MATVAR');
  1243. ZDFINI = ZDFINI '+' DFTHE;
  1244. 'FINS';
  1245. 'SI' LOGDEF;
  1246. ZDFINI = ZDFINI '+' DFDEF;
  1247. 'FINS';
  1248. *
  1249. ZSDMBR = RESIDU '*' 1.D0;
  1250. *
  1251. 'FINSI' ;
  1252. * FOR_MECA
  1253. *
  1254. 'SI' WTAB.'NVSTNL' ;
  1255. ZVIFL = conti.'VITESSES_FLUIDE' ;
  1256. movifl = 'EXTR' WTAB.'MOD_NSL' 'DEPL' ;
  1257. UN = 'EXCO' ZVIFL movifl ;
  1258. ch_un = 'MANU' 'CHML' WTAB.'MOD_NSL' 'ADVE' un 'TYPE' 'CARACTERISTIQUES';
  1259. nugrad = 'ADVE' WTAB.'MOD_NSL' ch_un ;
  1260. 'SI' ('EXIS' PRECED 'OPTION_BDF2') ;
  1261. ZRNSL = WTAB.'RNSL' 'ET' nugrad 'ET' (WTAB.'NSMA' '*' 1.5) ;
  1262. unm = 'EXCO' conti.'VITESSES_FLUIDE_0' movifl ;
  1263. fmass = WTAB.'NSMA' '*' ((2. '*' un) '-' (0.5 '*' unm)) ;
  1264. 'SINON' ;
  1265. ZRNSL = WTAB.'RNSL' 'ET' nugrad 'ET' WTAB.'NSMA' ;
  1266. fmass = WTAB.'NSMA' * un ;
  1267. 'FINSI' ;
  1268. ftnsl = ZFCONSTA 'ET' fmass ;
  1269.  
  1270. 'SI' LOGADV ;
  1271. ch_ad = 'MANU' 'CHML' WTAB.'MOD_NSL' 'ADVE' UADV 'TYPE' 'CARACTERISTIQUES';
  1272. nuadv = 'ADVE' WTAB.'MOD_NSL' ch_ad ;
  1273. ZRNSL = ZRNSL 'ET' nuadv ;
  1274. fadv = nugrad * UADV ;
  1275. mess 'advec / inertie' (maxi (abs fadv)) (maxi (abs ftnsl)) (maxi (abs UADV)) ;
  1276. ftnsl = ftnsl '-' fadv ;
  1277. 'FINSI' ;
  1278.  
  1279. resnsl = ftnsl '-' (ZRNSL * ZVIFL) ;
  1280.  
  1281. * resolution
  1282. ZVIF1 = 'KRES' ZRNSL resnsl ;
  1283. ZVIFL = ZVIFL '+' ZVIF1 ;
  1284. 'FINSI' ;
  1285. *
  1286. *----------------------------------------------------------------------*
  1287. * 1ere Resolution - ZDEP1 = [ Du^0 ; LX_n+1 ] *
  1288. * *
  1289. * -> estimation a TI a partir de l'instant TEMPS0 *
  1290. * *
  1291. * K*Du^0 + At*LX_n+1 = F^ext_n+1 - Bsigma_n *
  1292. * A*Du^0 = d1 - A*u_n *
  1293. * |----------------| |-------------------| *
  1294. * ZRAID * ZDEP1 = RESIDU *
  1295. *----------------------------------------------------------------------*
  1296. 'SI' WTAB.'FOR_MECA' ;
  1297. FEXCI = 'VIDE' 'CHPOINT' / 'DISCRET' ;
  1298. 'SI' (WTAB.'ADHERENCE') ;
  1299. 'SI' ('EXIS' WTAB 'MAIL_BLOM') ;
  1300. FADRE = 'REDU' FADHE WTAB.'MAIL_BLOM' ;
  1301. 'SINON' ;
  1302. FADRE = FADHE ;
  1303. 'FINSI' ;
  1304. FEXCI = FEXCI 'ET' FADRE ;
  1305. 'FINSI' ;
  1306. *
  1307. * Modification du residu si automatique
  1308. 'SI' LOGPIL ;
  1309. RESIDU = RESIDU 'ET' (ETA0 '*' ZFPILIN);
  1310. 'FINSI' ;
  1311. *
  1312. * Force limite de frottement
  1313. 'SI' (WTAB.'CAFROTTE' 'ET' IMPO12);
  1314. 'SI' WTAB.'FROCABL' ;
  1315. zsigfT = 'REDU' zsigf ZMODLI;
  1316. FFROT = 'EXCC' ZCLIM ZDEPT MODCON WTAB.'MATCONTA' ZSIGFT;
  1317. RTRSF = 'VIDE' 'RIGIDITE' / 'RIGIDITE' ;
  1318. RFNS = 'VIDE' 'RIGIDITE' / 'RIGIDITE' ;
  1319. FEXCI = FEXCI '+' FFROT ;
  1320. 'FINSI';
  1321. 'SI' WTAB.'FROCOUL' ;
  1322. excfconv = excfconv + 1 ;
  1323. RTRSF*'RIGIDITE' RFNS*'RIGIDITE' FFROT = 'EXCF' ZCLIM ZDEPT MODCON WTAB.'MATCONTA' excfconv;
  1324. FEXCI = FEXCI '+' FFROT ;
  1325. MFROT = 'EXTR' FFROT 'MAIL' ;
  1326. 'FINSI';
  1327. 'FINS';
  1328. *
  1329. 'SI' ('NEG' ZRIBLO_M 'INCONNU') ;
  1330. 'SI' ((WTAB.'CAFROTTE' 'OU' WTAB.'ADHERENCE') 'ET' IMPO12);
  1331. ZDEP1 BID BID BID BID = 'RESO' ZRAID 'SOUC' RESIDU 'INIB'
  1332. ZRIBLO_M ZLISEA_M FEXCI RFNS RTRSF ;
  1333. 'SINON';
  1334. ZDEP1 BID BID BID = 'RESO' ZRAID 'SOUC' RESIDU 'INIB'
  1335. ZRIBLO_M ZLISEA_M ;
  1336. 'FINS';
  1337. ZRIBLO_M = 'MOT' 'INCONNU' ;
  1338. 'SINON';
  1339. 'SI' ((WTAB.'CAFROTTE' 'OU' WTAB.'ADHERENCE') 'ET' IMPO12);
  1340. ZDEP1 = 'RESO' ZRAID 'SOUC' RESIDU FEXCI RFNS RTRSF ;
  1341. 'SINON';
  1342. ZDEP1 = 'RESO' ZRAID 'SOUC' RESIDU;
  1343. 'FINS';
  1344. 'FINS';
  1345. *
  1346. 'SI' LOGPIL ;
  1347. ZDEPII = 'RESO' ZRAID ZFPILIN ;
  1348. ZDEPILO = 0. * ZDEPII ;
  1349. * Calcul de D_eta par appel a la procedure PILOINDI
  1350. D_ETA = PILOINDI PRECED ZDEP0 ZDEPILO ZDEP1 ZDEPII DTAU;
  1351. * Mise à jour
  1352. ZDEP1 = ZDEP1 '+' (D_ETA '*' ZDEPII) ;
  1353. ETA = ETA0 ;
  1354. 'FINSI' ;
  1355. *
  1356. * Pour initialiser le statut des CL unil lors du prochain appel a RESO
  1357. ZRAID_T = 'EXTR' ZRAID 'CONT';
  1358. 'SI' ('EXIS' ZRAID_T 'NITER');
  1359. * Conserver les conditions unilaterales pour VITEUNIL
  1360. 'SI' IDYN;
  1361. WTAB.'ZRAIDV' = ZRAID;
  1362. 'FINS';
  1363. ZRIBLO_M = ZRAID_T. 7 ;
  1364. ZLISEA_M = ZRAID_T. 6 ;
  1365. 'FINS';
  1366. *
  1367. * Pour eventuel calcul de augauto
  1368. zdep1d = zdep1 'ENLE' 'LX';
  1369. * et initialisation zdepl (sert pour xnum)
  1370. ZDEPL = 'EXCO' ZDEP1 'LX' 'NOID' 'LX' ;
  1371. *----------------------------------------------------------------------*
  1372. * Calcul d'une norme pour la convergence *
  1373. *----------------------------------------------------------------------*
  1374. XXX1 = ZFEXT '+' ZFSUIV ;
  1375. 'SI' LOGPIL ;
  1376. XXX1 = XXX1 '+' ((ETA '+' D_ETA) '*' ZFPILIN);
  1377. 'FINSI' ;
  1378. *
  1379. ZDEP1P50 = ZDEP1 ;
  1380. XDENO='XTY' ZDEP1P50 ( XXX1 -( RESIDU 'EXCO'
  1381. 'FLX' 'NOID' 'FLX' 'NATURE' 'DISCRET')) MLPRIM MLDUAL;
  1382. 'SI' ('VERI' xdeno) ; 'SINON'; xdeno = 1; 'FINSI';
  1383. MZDEP1M = 'MAXI' ZDEP1P50 'ABS' 'AVEC' MLDEPL;
  1384. MZFM = 'MAXI' (FDENDEF + DFTHE + DFDEF + ZF1) 'ABS' 'AVEC' MLDUAL;
  1385. XDENO1 = 'ABS' XDENO + (MZFM * MZDEP1M);
  1386. MZDEP1M = MZDEP1M + XPETIT;
  1387. XDENO=XDENO1/MZDEP1M;
  1388. XDENO = XDENO + MZFM;
  1389. XDENO = 'MAXI' XDENO ('MAXI' XXX1 'ABS' 'AVEC' MLDUAL);
  1390. 'SI' (XDENO < XPETIT); XDENO = 1.; 'FINSI';
  1391. XDENOM=XDENO;
  1392. 'SI' TSTMOM ;
  1393. ZDEP1P50 = ZDEP1 + XPETIT; ;
  1394. XDENOM = XDENO1/('MAXI' ZDEP1P50 'ABS' 'AVEC' MLROTA);
  1395. XDENOM = XDENOM + XPETIT ;
  1396. 'FINS' ;
  1397. *
  1398. ** mess 'xdeno propose' ' ' xdeno ' xdeno precedent' ' ' xdenoo;
  1399. ** XDENO = MAXI XDENOo XDENO;
  1400. ** XDENOo = XDENO;
  1401. ** XDENOM = MAXI XDENOMo XDENOM;
  1402. ** XDENOMo = XDENOM;
  1403. *
  1404. 'SI' WTAB.'CONV';
  1405. ZINCREMENT = ZDFINI '+' ZFSUIV ;
  1406. 'SINON' ;
  1407. * 'MESS' 'non convergence on garde xdeno xdenom ' xdeno xdenom;
  1408. 'SI' (XDENO < ZXDENO); XDENO = ZXDENO ; 'FINSI';
  1409. 'SI' (XDENOM < ZXDENOM); XDENOM = ZXDENOM ; 'FINSI';
  1410. 'FINS';
  1411. *
  1412. RESIDNOR = 'COPIER' RESIDU ;
  1413. *
  1414. 'SI' IPILOT;
  1415. 'SI' WTAB.'AUTODEUX' ;
  1416. XDENO = ZXDENO ;
  1417. XDENOM = ZXDENOM ;
  1418. 'FINS';
  1419. DFEXT = DFEXT0F '*' COEPI ;
  1420. ZFEXT = DFEXT '+' FEXT0 ;
  1421. XXX1 = RESIDU '*' COEPI;
  1422. XXX1 = XXX1 '+' ((1. '-' COEPI) '*' FREAP) ;
  1423. 'SI' IMPO12;
  1424. RESIDU = 'ENLE' XXX1 'FLX' ;
  1425. XJEUX = DUIMPO '*' COEPI ;
  1426. RESIDU = RESIDU '+' XJEUX '+' DUUNIL ;
  1427. 'SINON' ;
  1428. RESIDU = XXX1 ;
  1429. 'FINS';
  1430. 'FINS';
  1431. *----------------------------------------------------------------------*
  1432. * Corriger le residu a partir du pas precedent *
  1433. *----------------------------------------------------------------------*
  1434. * petite correction du residu pour esperer gagner du temps ***********
  1435. INIT = FAUX ;
  1436. 'SI'(('NEG' ZFNONL 'INCONNU') 'ET' WTAB.'INITIALISATION');
  1437. 'SI' IPILOT;
  1438. 'SI' ( WTAB.'AUTODEUX' 'ET' (COEPI 'NEG' 1.D0)) ;
  1439. 'MESS' 'Initialisation a partir du pas precedent 'COEPI;
  1440. INIT = VRAI;
  1441. XXX1 = RED1 * ZFNONL ;
  1442. XXX2 = XXX1 + RESIDU ;
  1443. 'DETR' RESIDU ;'DETR' XXX1 ;
  1444. RESIDU = XXX2;
  1445. 'FINS' ;
  1446. 'SINON';
  1447. * on fait la correction si le pas precedent a converge sans non convergence
  1448. 'SI' (WTAB.'CONV' 'ET' (ISOUSPPP 'EGA' 0));
  1449. * on fait la correction si le pas precedent etait non lineaire.
  1450. 'SI' (('MAXI' 'ABS' ZFNONL) > (ZPREC * XDENO)) ;
  1451. * on enleve le residu du pas precedent pour recuperer l'increment
  1452. * nominal du second membre e imposer f1.
  1453. * f2 et l'increment du second membre du pas precedent
  1454. ZINCREMENT = ZINCREMENT - WTAB.'RESIDU' ;
  1455. zdeps = WTAB.'ZDEP1' + zdep1;
  1456. FFNO= 'XTY' ZFNONL zdeps MNDUAL MNPRIM;
  1457. F1 = ZINCREMENT ;
  1458. F2 = INCRPREC ;
  1459. f12 = 'XTY' f1 zdeps MNDUAL MNPRIM;
  1460. f22 = 'XTY' f2 zdeps MNDUAL MNPRIM;
  1461. AMPL = f12/(f22 + XPETIT);
  1462. DTPREC=1E30; AMPLT = 0;
  1463. 'SI' ('NEG' WTAB.'DTPREC' 0);
  1464. DTPREC = WTAB.'DTPREC' ;
  1465. 'SI' (DTPREC > XPETIT);
  1466. AMPLT = WTAB.'DT_INIT' '/' DTPREC;
  1467. 'FINS';
  1468. 'FINS';
  1469. * Le chargement n'est il pas de fluage ou de thermique ?';
  1470. XDCOMP = ('XTY' ZDEP1 F1 MNPRIM MNDUAL) ;
  1471. 'SI' (('ABS' XDCOMP) < (ZPREC * XDENO * mzdep1m)
  1472. 'OU' (('ABS' FFNO) > (('ABS' F22) * 2.e2 )));
  1473. DTPREC = WTAB.'DTPREC' ;
  1474. AMPL=AMPLT;
  1475. * la decharge est-elle significative
  1476. 'SI'((F12/(f22 + XPETIT)) < -0.05);ampl=0.;'FINS';
  1477. ** 'MESS' 'F12 F22_' f12 f22;
  1478. 'MESS' 'Pas d increment de charge, initialisation calculee avec le temps' ;
  1479. 'SINON';
  1480. AMPL = F12 / (F22 + XPETIT);
  1481. 'FINS';
  1482. * changement de modele on n'initialise pas
  1483. 'SI' ('NEG' ZMODLI WTAB.'MO_TOT_PREC'); AMPL = 0; 'FINSI';
  1484. AMPL = MINI (prog AMPL AMPLT);
  1485. 'SI' ((AMPL > 0) 'ET' (AMPL < 2e1)) ;
  1486. 'MESS' 'Initialisation a partir de la solution precedente Coeff'AMPL;
  1487. XXX1 = AMPL * ZFNONL ;
  1488. XXX2 =RESIDU+ XXX1;
  1489. 'DETR' XXX1;'DETR' RESIDU;
  1490. RESIDU = XXX2;
  1491. INIT = VRAI;
  1492. 'FINS';
  1493. 'FINS';
  1494. 'FINS';
  1495. 'FINS';
  1496. 'FINS';
  1497. WTAB.'ZDEP1'=zdep1;
  1498. *
  1499. * Initialisation pour l'acceleration de convergence (on peut mettre
  1500. * n'importe quoi c'est pour ne pas faire de tests dans la boucle)
  1501. ACFP1 = 'COPIER' ZFEXT2 *0. ;
  1502. ACFP2 = ACFP1 ;
  1503. ACFP3 = ACFP1 ;
  1504. ACFEP1 = ACFP1 ;
  1505. ACFEP2 = ACFP1 ;
  1506. FCORF = 'COPIER' FREAP ;
  1507. CORREC = 0;
  1508. *
  1509. zdept = (zdept 'EXCO' 'LX' 'NOID' 'LX') + (ZDEPT * 0.);
  1510. *----------------------------------------------------------------------*
  1511. * Initialisation des messages pour le pas de temps courant *
  1512. *----------------------------------------------------------------------*
  1513. 'SI' IKLFFF ;
  1514. 'MESSAGE'
  1515. 'Iter'*13 'Nplas'*26 'Critere'*39 'Deps.max'*52 'Eps.max'*65 'Crit.flex'*78
  1516. ;
  1517. 'SINON';
  1518. 'MESSAGE'
  1519. 'Iter'*13 'Nplas'*26 'Fresidu'*39 'Deps.max'*52 'Eps.max'*65 'Mresidu'*78
  1520. ;
  1521. 'FINS';
  1522. 'FINSI' ;
  1523. *FOR_MECA
  1524. *
  1525. XCONVMIN = 1e20;
  1526. NONCONV = FAUX;
  1527. PASREINI = VRAI;
  1528. IPREM = VRAI ; ITACC = 0 ;
  1529. RECA_N = 0 ; IT_RECA = 0 ;
  1530. IKT = FAUX ;
  1531. URG = FAUX ; RED_URG = 0 ;
  1532. ZICONV = VRAI; MMCMAX = 0 ;
  1533. DPSMAXP = 0; DEPSTDM = 0. ;
  1534. zprecnc=1e-5; ITNV = -5 ; PASTEST = FAUX; TABCONV = 'TABL';
  1535. IT = 0 ; MMC = 0 ; XCONV = 0. ; DPSMAX = xpetit ; EPSM = 0. ;
  1536. RE_CNTRL = 'VIDE' 'RIGIDITE' / 'RIGIDITE'; DE_CNTRL = FAUX;
  1537. *######################################################################*
  1538. *----------------------------------------------------------------------*
  1539. * Boucle de convergence ETIQ (Configuration GEOM1) *
  1540. *----------------------------------------------------------------------*
  1541. * preparation pour la gestion du pas avant etiq
  1542. depstp zdeptp zsigfp fcorfp correcp = depst zdept zsigf fcorf 0;
  1543. * initialisation sous decoupage pour epsi
  1544. n_sous_dec = 1 ;
  1545. 'REPETER' ETIQ ;
  1546. * A tester mais ne fonctionne pas en l'etat
  1547. * nsoincr = n_sous_dec;
  1548. *
  1549. * Champ materiau final
  1550. ZMAT2F = ZMAT2 ;
  1551. *
  1552. * partir du bon etat (mis a jour si le sous increment est acceptable)
  1553. depst zdept zsigf fcorf correc = depstp zdeptp zsigfp fcorfp correcp;
  1554.  
  1555. 'SI' WTAB.'FOR_MECA' ;
  1556. zdep1 = (zdept 'EXCO' 'LX' 'NOID' 'LX') + (ZDEPT * 0.);
  1557. 'FINSI' ;
  1558. *
  1559. *
  1560. 'SI' autaug;
  1561. resmul = 2. * resmul;
  1562. 'SI' ((dpsmax < 1d-3) 'ET' (resmul > 1D-20) 'ET' ('NON' PASTEST));
  1563. resmul = resmul / (dpsmax + xpetit) * 1d-3;
  1564. 'FINSI';
  1565. 'SI' (resmul > 1.0); resmul = 1.0; 'SINON'; itacc = 4; 'FINSI';
  1566. 'FINSI';
  1567. *
  1568. INSTAB = FAUX;
  1569. nconvr = FAUX;
  1570. PASTEST = FAUX;
  1571. RECA_K = FAUX;
  1572. *
  1573. * IT est le compteur de ETIQ, ITACC doit etre =< 0 pour qu'on accelere
  1574. IT = IT + 1 ;
  1575. ITACC = ITACC - 1;
  1576. *
  1577. * Pour conserver les criteres de convergence
  1578. tabconv. it = 1;
  1579. *----------------------------------------------------------------------*
  1580. * Recalcul de la rigidite si necessaire *
  1581. *----------------------------------------------------------------------*
  1582. 'SI' IGRD ;
  1583. 'FORM' GEOM1;
  1584. *
  1585. 'SI' ((IT > 1) 'ET' (URG 'OU' (ITACC > 3) 'OU' INSTAB 'OU'
  1586. ('MULT' (IT - IT_RECA) ITRCLC) ));
  1587. HPP_EPS = FAUX;
  1588. URG = FAUX;
  1589. IT_RECA = IT;
  1590. INSTAB = FAUX;
  1591. PASREINI = FAUX;
  1592. * zdeptr = zdept - (zdep1/2);
  1593. zdeptr = zdept ;
  1594. GEOR ZMATTEMP = 'FORM' ZDEPTR ZMODLI ZMAT22 ;
  1595. * A cause de FORM, ZMATTEMP n'est plus parallele...
  1596. ZMATTEMP = 'REDU' ZMATTEMP ZMODL;
  1597. RECA_K = VRAI;
  1598. RECA_N = RECA_N + 1;
  1599. 'SI' (RECA_N > 20) ;
  1600. nonconv = vrai;
  1601. 'FINS';
  1602. *
  1603. * ----------------------------------------------------------------
  1604. * Recalcul de la rigidite a la fin du pas
  1605. txt_k ='CHAI' ' Recalcul de K (= K^el';
  1606. 'SOUC' 0;
  1607. 'SI' ('EGA' LAG_TOT 1) ;
  1608. 'FORM' GEOREF0;
  1609. HOOKRH = 'HOOK' ZMODL ('REDU' ZMAT2I ZMODL) ;
  1610. HOOKRH2 = 'PICA' HYPDEF HOOKRH ZMODL (ZDEP0 + ZDEPTR);
  1611. 'FORM' GEOR;
  1612. RITC = 'RIGI' HOOKRH2 ZMODL zmattemp 'NOER';
  1613. 'DETR' HOOKRH; 'DETR' HOOKRH2;
  1614. 'SINON';
  1615. RITC = 'RIGI' ZMODL ZMATTEMP 'NOER' ;
  1616. 'FINSI';
  1617. 'FORM' GEOM1;
  1618. ZMATTEMP = 0 ;
  1619. PASOK = 'SOUCI';
  1620. 'SI' PASOK;
  1621. 'MESS' 'rigi rate. on reessaye avec la configuration geom1';
  1622. RITC = 'RIGI' ZMODL ZMAT2F ;
  1623. 'FINSI';
  1624. *
  1625. * RITC peut contenir des CL
  1626. ZRI = 'EXTR' RITC 'RIGI' 'NOMU' ;
  1627. ZCL = 'EXTR' RITC 'RIGI' 'MULT' ;
  1628. 'SI' ('NEG' ('DIME' ZCL) 0) ;
  1629. ZCLIM = ZCLIM0 'ET' ZCL ;
  1630. 'SINON' ;
  1631. ZCLIM = ZCLIM0 ;
  1632. 'FINSI' ;
  1633. *
  1634. * ----------------------------------------------------------------
  1635. * Recalcul de la masse, option 'AUGMENTATION_AUTOMATIQUE'
  1636. 'SI' (AUTAUG 'ET' ('NON' IRAUGLU));
  1637. 'FORM' GEOR ;
  1638. RIG_AUG = 'MASSE' ZMODL ZMAT2F ;
  1639. RH = RITC;
  1640. *** mess 'actualisation rig_aug';
  1641. 'FINSI';
  1642. *
  1643. * ----------------------------------------------------------------
  1644. * Ajout de relations issues de la procedure DEFO_IMP
  1645. ZFCONT = 0 ;
  1646. IACTURES = FAUX ;
  1647. CDEPSLX = 'ENLE' (ZDEP0 '+' ZDEPT) 'LX' ;
  1648.  
  1649. ** DE_CNTRL = FAUX;
  1650. 'SI' (AUTAUG 'ET' DE_CNTRL);
  1651. IACTURES = VRAI ;
  1652. ZFCONT = RE_CNTRL '*' CDEPSLX 'EXCO' 'FLX' 'NOID' 'FLX' ;
  1653. ZCLIM = RE_CNTRL 'ET' ZCLIM ;
  1654. txt_k = 'CHAI' txt_k ' + K^ctr';
  1655. 'FINSI';
  1656. *
  1657. * ----------------------------------------------------------------
  1658. * Recalcul des conditions de contact-frottement ds la conf finale
  1659. 'SI' WTAB.'CONTACT';
  1660. 'FORM' GEOR ;
  1661. ** 'FORM' ZDEPT ;
  1662. *
  1663. ZFLX = 'EXCO' (CRR '*' ZDEPT) 'FLX' 'NOID' 'FLX' ;
  1664. *
  1665. MODCON = wtab.'MODCONTA';
  1666. 'SI' ('EXIS' WTAB 'MATCONTA');
  1667. CJEU CRR = 'RFCO' MODCON WTAB.'CONV' WTAB.'MATCONTA';
  1668. 'SINON' ;
  1669. CJEU CRR = 'RFCO' MODCON WTAB.'CONV' ;
  1670. 'FINS' ;
  1671. *
  1672. * Separer ce qui concerne l'adherence (FADH) et les jeux (FLX)
  1673. 'SI' ('NEG' CJEU 0) ;
  1674. FADHE = 'EXCO' CJEU 'FADH' 'NOID' 'FADH' ;
  1675. CDAP = 'EXCO' CJEU 'FLX' 'NOID' 'FLX' ;
  1676. *
  1677. * Ajout du glissement deja parcouru sur les conditions de frot.
  1678. 'SI' WTAB.'CAFROTTE';
  1679. 'SI' ('EGA' ('TYPE' MFROT) 'MAILLAGE');
  1680. CGLI = 'REDU' ZFLX MFROT ;
  1681. CDAP = CDAP ET (-1*CGLI) ;
  1682. 'FINSI';
  1683. 'FINSI';
  1684. 'FINSI';
  1685. *
  1686. 'SI' WTAB.'MODAL' ;
  1687. CRR = 'PJBA' CRR ZMODL ZMAT2F ;
  1688. MCRR = 'EXTR' CRR 'MAIL' 'MULT';
  1689. MCDAP = 'EXTR' CDAP 'MAIL' ;
  1690.  
  1691. 'REPETER' BCDA ('NBNO' MCDAP) ;
  1692. PBCDA = MCDAP 'POINT' &BCDA ;
  1693. PCRR ='POINT' MCRR &BCDA ;
  1694. CHCR ='MANU' 'CHPO' PCRR 1 'FLX' ('EXTR' CDAP 'FLX' PBCDA)
  1695. 'NATURE' 'DISCRETE' ;
  1696. 'SI' ('EGA' 1 &BCDA) ;
  1697. CCDA = CHCR ;
  1698. 'SINON' ;
  1699. CCDA = CCDA 'ET' CHCR ;
  1700. 'FINS' ;
  1701. 'FIN' BCDA ;
  1702.  
  1703. CDAP = CCDA ;
  1704. 'FINS' ;
  1705. *
  1706. 'SI' ('NEG' CRR 0) ;
  1707. IACTURES = VRAI ;
  1708. CHPZ = CRR '*' CDEPSLX 'EXCO' 'FLX' 'NOID' 'FLX' ;
  1709. 'SI' ('NEG' ZFCONT 0) ;
  1710. ZFCONT = ZFCONT '+' CHPZ ;
  1711. 'SINON' ;
  1712. ZFCONT = CHPZ ;
  1713. 'FINSI' ;
  1714. ZCLIM = CRR 'ET' ZCLIM ;
  1715. 'FINS';
  1716. *
  1717. 'SI' ('NEG' CJEU 0) ;
  1718. IACTURES = VRAI ;
  1719. ZFCONT = ZFCONT '+' CDAP ;
  1720. 'FINS';
  1721. 'FINS';
  1722. *
  1723. * Mise a jour des jeux dans le residu (partie force inchangee)
  1724. ZFCONSTA = ZFEXT2 ;
  1725. ZFLX1 = ZFLXB ;
  1726. 'SI' IACTURES ;
  1727. ZFCONSTA = ZFCONSTA '+' ZFCONT ;
  1728. RFORCE = 'ENLE' RESIDU 'FLX' ;
  1729. FLXPRE = 'EXCO' (ZCLIM '*' (ZDEP0 '+' ZDEPT)) 'FLX' 'NOID' 'FLX';
  1730. RJEUX = ZFCONT '+' ZFLXB '-' FLXPRE ;
  1731. RESIDU = RFORCE 'ET' RJEUX ;
  1732. ZFLX1 = ZFLX1 'ET' ZFCONT ;
  1733. 'FINS';
  1734. *
  1735. ZRAID = ZCLIM 'ET' ZRI;
  1736. * Mise a jour de FREAP
  1737. FREAP = (ZDEP0 'EXCO' 'LX' 'LX') * ZCLIM;
  1738.  
  1739.  
  1740.  
  1741. *
  1742. * ----------------------------------------------------------------
  1743. * Ajout de la partie dynamique
  1744. 'SI' IDYN;
  1745. * bp : en toute rigueur, il faudrait aussi recalculer la MASSE ...
  1746. * et ajouter l'amortissement le cas échéant ...
  1747. ZRAID = ZRAID 'ET' ZRMAS ;
  1748. 'FINS';
  1749. *
  1750. * ----------------------------------------------------------------
  1751. * Recalcul de la raideur geometrique (option K_SIGMA)
  1752. 'SI' IKSIA ;
  1753. 'FORM' GEOR;
  1754. ksigtc = 'KSIGM' ZSIGF ZMODL ZMAT2F ;
  1755. ZRAID = ZRAID 'ET' ksigtc ;
  1756. txt_k = 'CHAI' txt_k ' + K^sig';
  1757. 'FINS';
  1758. *
  1759. * ----------------------------------------------------------------
  1760. * Prise en compte d'eventuelles 'RIGIDITE_CONSTANTE'
  1761. 'SI' IRCON;
  1762. ZRAID = ZRAID 'ET' RIG_CONS;
  1763. txt_k = 'CHAI' txt_k ' + K^cst';
  1764. 'FINS';
  1765. *
  1766. * ----------------------------------------------------------------
  1767. * Traitement en de la rigidite augmentee
  1768. 'SI' IRAUG;
  1769. 'SI' iprem; zdep1d = ZDEP0 'ENLE' 'LX'; 'FINSI';
  1770. xkx = (xtmx zdep1d rh) + xpetit;
  1771. xmx = (xtmx zdep1d rig_aug) + xpetit;
  1772. 'SI' (('VERI' xkx) 'ET' ('VERI' xmx) 'NON');
  1773. zu1l = ZDEP0 'ENLE' 'LX';
  1774. xkx = (xtmx zu1l rh) + xpetit;
  1775. xmx = (xtmx zu1l rig_aug) + xpetit;
  1776. 'FINSI';
  1777. augauto = xkx / xmx; augk = 1.d-3;
  1778. augauto = abs augauto;
  1779. augk = abs augk;
  1780. augmult = augmult * 1.02;
  1781. 'SI' (augmult > 1D4 ); augmult=1D4 ; 'FINSI';
  1782.  
  1783. augauto = augauto * augmult ;
  1784. augk = augk * augmult;
  1785.  
  1786. 'SI' (iraug 'ET' autaug) ;
  1787. 'MESS' 'multiplicateur d augmentation masse' augauto ' raideur' augk ;
  1788. 'FINSI';
  1789. ZRAIDNA= ZRAID ;
  1790. 'SI' AUTAUG;
  1791. ZRAID = ZRAID 'ET' (RIG_AUG * augauto) 'ET' (RH * augk) ;
  1792. 'SINON';
  1793. ZRAID = ZRAID 'ET' RIG_AUG ;
  1794. *** 'MESS' ' augmentation residuelle ';
  1795. 'FINSI';
  1796. txt_k = 'CHAI' txt_k ' + K^aug';
  1797. *** ZNACCE = FAUX;
  1798. 'SINON';
  1799. augmult = augmult * 0.55 ;
  1800. si (augmult < 1d-3); augmult = 1d-3; finsi;
  1801. 'SI' iraug;
  1802. 'MESS' 'nouveau augauto augmult' ' ' augauto augmult;
  1803. 'FINSI';
  1804. 'SI' IRAUGLU;
  1805. zraid = zraid 'ET' rig_aug;
  1806. 'FINSI';
  1807. 'FINSI';
  1808. 'FORM' GEOM1;
  1809. *
  1810. XXX1 = 'EXTR' ZCLIM 'MAIL' 'UNIL';
  1811. IMPO12 = ('NBEL' XXX1) '>' 0 ;
  1812. 'SI' (IPILOT 'ET' IMPO12) ;
  1813. DUUNIL = 'REDU' DFEXT0L XXX1 ;
  1814. XXX2 = 'DIFF' ('EXTR' ZCLIM 'MAIL') XXX1 ;
  1815. DUIMPO = 'REDU' DFEXT0L XXX1 ;
  1816. 'FINS';
  1817. *
  1818. * ----------------------------------------------------------------
  1819. * Procedure UTILISATEUR
  1820. 'SI' WTAB.'PROCEDURE_CHARMECA';
  1821. PRECED.'ADDI_MATRICE' = vrai;
  1822. TFP22= CHARMECA PRECED TI ;
  1823. PRECED.'ADDI_MATRICE' = faux;
  1824. 'SI' (EXIS TFP22 'ADDI_MATRICE');
  1825. zraid = zraid 'ET' TFP22.'ADDI_MATRICE';
  1826. txt_k = 'CHAI' txt_k ' + K^cent';
  1827. 'FINS';
  1828. 'FINS';
  1829. *
  1830. 'MESS' ('CHAI' txt_k ' ) dans config deformee ') ;
  1831. *
  1832. * on impose le recalcul de K a la prochaine iteration si it=2 (pq ?)
  1833. 'SI' (IT 'EGA' 2) ;
  1834. ITACC = 4 ;
  1835. 'FINS' ;
  1836. 'FINS';
  1837. 'FINS';
  1838. *
  1839. 'SI' WTAB.'NVSTNL' ;
  1840. UN = 'EXCO' ZVIFL movifl ;
  1841. ch_un = 'MANU' 'CHML' WTAB.'MOD_NSL' 'ADVE' un 'TYPE' 'CARACTERISTIQUES';
  1842. nugrad = 'ADVE' WTAB.'MOD_NSL' ch_un ;
  1843. 'SI' ('EXIS' PRECED 'OPTION_BDF2') ;
  1844. ZRNSL = WTAB.'RNSL' 'ET' nugrad 'ET' (WTAB.'NSMA' '*' 1.5) ;
  1845. 'SINON' ;
  1846. ZRNSL = WTAB.'RNSL' 'ET' nugrad 'ET' WTAB.'NSMA' ;
  1847. 'FINSI' ;
  1848. ftnsl = ZFCONSTA 'ET' fmass ;
  1849.  
  1850. 'SI' LOGADV ;
  1851. ch_ad = 'MANU' 'CHML' WTAB.'MOD_NSL' 'ADVE' UADV 'TYPE' 'CARACTERISTIQUES';
  1852. NUADV = 'ADVE' WTAB.'MOD_NSL' ch_ad ;
  1853. ZRNSL = ZRNSL 'ET' NUADV ;
  1854. ftnsl = ftnsl '-' fadv ;
  1855. 'FINSI' ;
  1856.  
  1857. resnsl = ftnsl '-' (ZRNSL * ZVIFL) ;
  1858.  
  1859. 'FINSI' ;
  1860. *----------------------------------------------------------------------*
  1861. * Acceleration de convergence effective *
  1862. *----------------------------------------------------------------------*
  1863. 'SI' WTAB.'FOR_MECA' ;
  1864. CORRECA = CORREC ;
  1865. CORREC = 0;
  1866. ACFEP0 = RESIDU - FCORF;
  1867. ACFEP0 = ACFEP0 'ENLE' 'FLX';
  1868. ACFEPT = ACFEP0;
  1869. ACFP0 = ACFEP0 * ZJAC ;
  1870. ACFEP0 = ACFEP0 - CORRECA;
  1871. ACCEL = FAUX;
  1872. 'SI' (('NON' instab) 'ET' (ITACC '&lt;EG' 0) 'ET' (IT '>' 3) 'ET' ZNACCE );
  1873. ITACC = 2;
  1874. ACCEL = VRAI;
  1875. CORREC = 'ACT3' ACFEP2 ACFEP1 ACFEP0
  1876. ACFP3 ACFP2 ACFP1 ACFP0 ;
  1877. *
  1878. * verif que l'acceleration ne renvoie pas en arriere
  1879. 'SI' (wtab.'STABILITE' );
  1880.  
  1881. acc_ref = xty acfept zdep1d MNDUAL MNPRIM;
  1882. acc_ref = acc_ref + xpetit;
  1883.  
  1884. acc_dir = xty (acfept - correc) zdep1d MNDUAL MNPRIM;
  1885. acc_rap = acc_dir/acc_ref;
  1886. acc_lim = 256.; 'SI' (acc_ref < 0.); acc_lim = 1.; 'FINSI';
  1887. 'SI' (acc_rap '>' acc_lim) ;
  1888. * en cas d'acceleration trop grande, on la limite
  1889. 'MESS' 'Limitation acceleration ' ' 'acc_rap; CORREC = correc * (acc_lim/acc_rap);
  1890. 'FINSI';
  1891. 'SI' (acc_rap '&lt;EG' 0.) ;
  1892. * en cas d'acceleration retrograde, on n'accelere pas
  1893. 'MESS' 'Annulation acceleration: retrograde' ' ' acc_rap;
  1894. correc = 0.;
  1895. 'FINSI';
  1896. 'FINSI';
  1897. *
  1898. RESIDU = RESIDU '-' CORREC;
  1899. 'FINS';
  1900. *
  1901. ACFP3 = ACFP2 ;
  1902. ACFP2 = ACFP1 ;
  1903. ACFP1 = ACFP0 ;
  1904. ACFEP2 = ACFEP1 ;
  1905. ACFEP1 = ACFEP0 ;
  1906. 'FINSI' ;
  1907. *FOR_MECA
  1908. *
  1909. 'SI' (resmul < 0.99);
  1910. 'MESS' 'Reduction du chargement. Coefficient: ' resmul;
  1911. 'SI' (resmul < 1.d-40); 'ERREUR' 996; 'FINSI';
  1912. residu = ((residu -fcorf) * resmul) + fcorf;
  1913. ACFP0 = (RESIDU - FCORF) ;
  1914. 'FINSI';
  1915. 'SOUCI' 0;
  1916. *----------------------------------------------------------------------*
  1917. * Systeme a resoudre a l'iteration (i) *
  1918. * *
  1919. * | K At | | du | | Fsur + Fvol - Bsigma^(i) | *
  1920. * | | | | = | | *
  1921. * | A 0 | | LX(i) | | d1 - A*u^(i) | *
  1922. *----------------------------------------------------------------------*
  1923. 'SI' WTAB.'FOR_MECA' ;
  1924. FEXCI = 'VIDE' 'CHPOINT' / 'DISCRET' ;
  1925. 'SI' (WTAB.'ADHERENCE') ;
  1926. 'SI' ('EXIS' WTAB 'MAIL_BLOM') ;
  1927. FADRE = 'REDU' FADHE WTAB.'MAIL_BLOM' ;
  1928. 'SINON' ;
  1929. FADRE = FADHE ;
  1930. 'FINSI' ;
  1931. FEXCI = FEXCI 'ET' FADRE ;
  1932. 'FINSI' ;
  1933. *
  1934. * Force limite de frottement
  1935. 'SI' (WTAB.'CAFROTTE' 'ET' IMPO12);
  1936. 'SI' WTAB.'FROCABL' ;
  1937. zsigfT = 'REDU' zsigf ZMODLI;
  1938. FFROT = 'EXCC' ZCLIM ZDEPT MODCON WTAB.'MATCONTA' ZSIGFT;
  1939. RTRSF = 'VIDE' 'RIGIDITE' / 'RIGIDITE';
  1940. RFNS = 'VIDE' 'RIGIDITE' / 'RIGIDITE';
  1941. FEXCI = FEXCI '+' FFROT ;
  1942. 'FINSI';
  1943. 'SI' WTAB.'FROCOUL' ;
  1944. excfconv = excfconv + 1 ;
  1945. RTRSF*'RIGIDITE' RFNS*'RIGIDITE' FFROT = 'EXCF' ZCLIM ZDEPT MODCON WTAB.'MATCONTA' excfconv;
  1946. FEXCI = FEXCI '+' FFROT ;
  1947. MFROT = 'EXTR' FFROT 'MAIL' ;
  1948. 'FINSI';
  1949. 'FINS';
  1950. *
  1951. 'SI' ( ('NEG' ZRIBLO_M 'INCONNU') 'ET' (IPREM 'OU' RECA_K) ) ;
  1952. 'SI' ((WTAB.'CAFROTTE' 'OU' WTAB.'ADHERENCE') 'ET' IMPO12);
  1953. ZDEP1 BID BID BID BID = 'RESO' ZRAID 'SOUC' RESIDU 'INIB'
  1954. ZRIBLO_M ZLISEA_M
  1955. FEXCI RFNS RTRSF;
  1956. 'SINON';
  1957. ZDEP1 BID BID BID = 'RESO' ZRAID 'SOUC' RESIDU 'INIB'
  1958. ZRIBLO_M ZLISEA_M ;
  1959. 'FINS';
  1960. 'SINON';
  1961. 'SI' ((WTAB.'CAFROTTE' 'OU' WTAB.'ADHERENCE') 'ET' IMPO12);
  1962. ZDEP1 = 'RESO' ZRAID 'SOUC' RESIDU FEXCI RFNS RTRSF;
  1963. 'SINON';
  1964. ZDEP1 = 'RESO' ZRAID 'SOUC' RESIDU 'NOID';
  1965. 'SI' (WTAB.'MAN' 'ET' IPREM);
  1966. ORDRE = WTAB.'ORDRE' ;
  1967. ZDEP2 IOUT = CORMAN ZRAIDINI ZMODL ZMAT2F ORDRE
  1968. ZDEP0 ZSIG0 RESIDNOR WTAB ;
  1969. 'SI' (IOUT 'EGA' 1) ;
  1970. ZDEP1=ZDEP2;
  1971. 'FINS';
  1972. 'FINSI' ;
  1973. 'FINS';
  1974. 'FINS';
  1975. *
  1976. * verif stabilite matrice elastique
  1977. 'SI' (wtab.'STABILITE' 'ET' ('NON' HPP_EPS) 'ET' AUTAUG);
  1978. RESIDD = RESIDU - (ZRAID * (ZDEP1 'EXCO' 'LX' 'LX'));
  1979. acc_ref= xty (RESIDD ) zdep1 MNDUAL MNPRIM;
  1980. ** diag refactorise la matrise. On evite donc provisoirement de l'appeler
  1981. ** nbng = 'DIAG' ZRAID;
  1982. nbng = 0;
  1983. 'SI' ((acc_ref < 0.) 'OU' (nbng > 0));
  1984. 'MESS' 'Raideur negative' ' ' acc_ref ;
  1985. URG = VRAI;
  1986. HPP_EPS = FAUX;
  1987. iraug = vrai;
  1988. ** resmul = resmul * 0.25;
  1989. pastest = VRAI;
  1990. instab = vrai;
  1991. itacc = 4;
  1992. augmult = augmult * 1.4;
  1993. 'SI' FAUX;
  1994. * controler la deformation max
  1995. DEPST = 'EPSI' 'LINE' ZMODL (ZDEPT + ZDEP1) ZMAT2F;
  1996. DPSMAX = 'MAXI' 'ABS' DEPST 'AVEC' MLDEFOR;
  1997. RE_CNTRL2*'RIGIDITE' = DEFO_IMP DEPST (zdept + zdep1) ZRAID ;
  1998. DE_CNTRL = FAUX;
  1999. ZFCONT = RE_CNTRL2 '*' ( ZDEP1) 'EXCO' 'FLX' 'NOID' 'FLX' ;
  2000. cmul = MINI 1. ( 1d-3 / DPSMAX);
  2001. ZFEXT2 = ZFEXT2 + (ZFCONT * cmul);
  2002. ZFLXB = 'EXCO' ZFEXT2 'FLX' 'NOID' 'FLX' ;
  2003. ZFCONSTA = ZFEXT2 ;
  2004. RE_CNTRL = RE_CNTRL 'ET' RE_CNTRL2;
  2005. 'FINSI';
  2006. 'ITERER' ETIQ;
  2007. 'FINSI';
  2008. 'FINSI';
  2009. ** 'FINSI';
  2010. monsouc ='SOUCI';
  2011. 'SI' monsouc;
  2012. 'SI' ('NON' autaug);
  2013. 'MESS' 'Souci': (vale 'SOUCI');
  2014. HPP_EPS = FAUX;
  2015. pastest = VRAI;
  2016. instab = vrai;
  2017. ** IKSIA = FAUX;
  2018. URG=VRAI;
  2019. itacc = 4;
  2020. *** 'ITERER' ETIQ;
  2021. 'SINON';
  2022. 'MESS' 'Annulation du sous-increment: souci dans la resolution';
  2023. URG = VRAI;
  2024. iraug = vrai;
  2025. ** resmul = resmul * 0.25;
  2026. HPP_EPS = FAUX;
  2027. pastest = VRAI;
  2028. augmult = augmult * 1.4;
  2029. itacc = 4;
  2030. 'ITERER' ETIQ;
  2031. 'FINSI';
  2032. 'FINSI';
  2033. *
  2034. 'SI' LOGPIL ;
  2035. * Calcul de D_eta par appel a la procedure PILOINDI
  2036. D_ETA = PILOINDI PRECED ZDEP0 ZDEPT ZDEP1 ZDEPII DTAU ;
  2037. * Mise à jour
  2038. ZDEP1 = ZDEP1 '+' (D_ETA '*' ZDEPII) ;
  2039. ETA = ETA '+' D_ETA ;
  2040. 'FINSI' ;
  2041. *
  2042. * verif sens ----------------------------------------------------------
  2043. *
  2044. 'SI' (resmul < 1d-2) ; souci 0; 'FINSI';
  2045. zdep1s = zdep1 'ENLE' 'LX';
  2046. monsouc = souci;
  2047. pasunil = monsouc;
  2048. *
  2049. ZRAID_T = 'EXTR' ZRAID 'CONT';
  2050. 'SI' ('EXIS' ZRAID_T 'NITER');
  2051. * Conserver les conditions unilaterales pour VITEUNIL
  2052. 'SI' IDYN ; WTAB.'ZRAIDV' = ZRAID; 'FINS';
  2053. ZRIBLO_M = ZRAID_T.7;
  2054. ZLISEA_M = ZRAID_T.6;
  2055. WTAB.'MAIL_BLOM' = 'EXTR' ZRIBLO_M 'MAIL' 'MULT' ;
  2056. 'FINS';
  2057. 'FINSI' ;
  2058. *FOR_MECA
  2059.  
  2060. 'SI' WTAB.'NVSTNL' ;
  2061. ZVIF1 = 'KRES' ZRNSL Resnsl ;
  2062. ZVIFL = ZVIFL '+' ZVIF1 ;
  2063. 'FINSI' ;
  2064. *----------------------------------------------------------------------*
  2065. * Increment de deplacements et mult. de Lagrange *
  2066. * *
  2067. * ZDEP1 est un increment d'increment de deplacements, note du, et les *
  2068. * multiplicateurs de Lagrange sont resolus "en total". *
  2069. * La solution de la boucle ETIQ a l'iteration i a fourni : *
  2070. * ZDEP1 = [ du ; LX(i) ] *
  2071. * Le cumul est realise de maniere a avoir : *
  2072. * ZDEPT = [ Du^(i) ; LX(i) ] *
  2073. * avec Du^(i) = Du^(i-1) + du *
  2074. * Le champ de deplacements total s'ecrit : *
  2075. * ZDETOT = ZDEP0 '+' ZDEPT *
  2076. *----------------------------------------------------------------------*
  2077. 'SI' WTAB.'FOR_MECA' ;
  2078. 'SI' IPREM ;
  2079. ZDEPT = 'COPIER' ZDEP1 ;
  2080. ZDELA = 'COPIER' ZDEPT ;
  2081. 'SINON';
  2082. XXX1 = ZDEPT 'ENLEVER' 'LX' ;
  2083. ZDEPT = XXX1 '+' ZDEP1 ;
  2084. 'DETR' XXX1 ;
  2085. 'FINS' ;
  2086. *
  2087. 'SI' WTAB.'MODAL' ;
  2088. 'SI' ('EXIS' WTAB 'MODCONTA') ;
  2089. 'SI' IPREM ;
  2090. mamoco1 = 'EXTR' (ZDEPT 'ENLEVER' 'LX') 'MAIL' ;
  2091. 'SINON' ;
  2092. * mettre les point materiels dans zdept
  2093. zdeptu1 = 'REDU' zdept mamoco1 ;
  2094. ch_dco = 'RECO' zdeptu1 ZMODL ZMAT22 ;
  2095. ZDEPT = ('EXCO' zdept 'LX' 'LX') 'ET' zdeptu1 'ET' ch_dco ;
  2096. 'FINS';
  2097. 'FINS';
  2098. 'FINS' ;
  2099. *
  2100. * On part dans le decor : redemarrer a 0
  2101. 'SI' ((XCONV > 1E8) 'ET' PASREINI 'ET' ('NON' IPILOT)) ;
  2102. 'MESS' 'Reinitialisation du schema';
  2103. PASREINI=FAUX;
  2104. ITACC=4;
  2105. 'ITERER' ETIQ;
  2106. 'FINS';
  2107. *
  2108. * Pilotage automatique
  2109. 'SI' IPILOT;
  2110. 'SI' ACCEL ;
  2111. RED1 = ZAUTOREDU ;
  2112.  
  2113. * reduction du critere de pilotage red2 sert d'incateur si
  2114. * pas accelere tous les 2 pas
  2115.  
  2116. * MODIFICATION CB215821 : 18/06/2015
  2117. 'SI' (RED2 < (IT '/' 20));
  2118. ZAUTOREDU = ZAUTOREDU '*' 3.D0 ;
  2119. RED2 = RED2 '+' 1 ;
  2120. ITACC= 3;
  2121. 'FINS';
  2122.  
  2123. RED1 = RED1 '/' ZAUTOREDU ;
  2124. 'SI' (RED1 '<' 0.9) ; ZAUTORED1 = 10 ; 'FINS' ;
  2125.  
  2126. 'SI' (ZAUTOREDU '>' 1.D0 ) ;
  2127. 'MESS' 'On divise le critere de pilotage par 'ZAUTOREDU;
  2128. 'FINS' ;
  2129.  
  2130. OO = WTAB.'AUTOCRIT' '/' ZAUTOREDU ;
  2131. U1MA = AUTOPILO ZDEPT (COEPI'*'ZDELA) ZMODLI ZMAT22 WTAB;
  2132. AL1 = OO '/' U1MA ;
  2133. * al1 est le coefficient de normalisation
  2134. 'SI' (( al1 '>EG' 1.d0) 'ET' (coepi '>EG' 1.d0)) ;
  2135. * pour eviter d'aller au dela de alpha=1 on ignore le critere
  2136. AL1 = 1.d0 ;
  2137. 'FINS' ;
  2138. 'SI'((AL1 '>' 1.D0) 'ET' (COEPI '>' 0.D0));
  2139. 'SI' ( AL1 '>' (1.D0 '/' COEPI));
  2140. AL1 = 1.D0 '/' COEPI;
  2141. 'FINS';
  2142. 'FINS';
  2143. * normalisation
  2144. PASTEST = VRAI;
  2145. XXX1 = ZLX '*' ( 1.d0 '-' AL1);
  2146. XXX3 = ZDEPT '*' AL1 ;
  2147. ZDEPT = XXX3 '+' XXX1;
  2148. 'DETR' XXX3;
  2149. 'SINON';
  2150. * le calcul d'al1 est fait pour eviter la convergence
  2151. * s'il est externe a l'interval 0.5 --- 2.
  2152. OO = WTAB.'AUTOCRIT' '/' ZAUTOREDU ;
  2153. U1MA = AUTOPILO ZDEPT (COEPI'*'ZDELA) ZMODLI ZMAT22 WTAB;
  2154. AL1 = OO '/' U1MA ;
  2155. 'SI' (( AL1 '>EG' 1.d0) 'ET' (COEPI '>EG' 1.d0)) ;
  2156. AL1 = 1.d0 ;
  2157. 'FINS' ;
  2158. 'FINS' ;
  2159. 'FINS';
  2160. *
  2161. * garder les reactions pour le test de convergence
  2162. ZDEPLP = ZDEPL ;
  2163. ZDEPL = ZDEPT 'EXCO' 'LX' 'NOID' 'LX' ;
  2164. ZGRDUF = 'MOT' 'INCONNU';
  2165. *
  2166. *----------------------------------------------------------------------
  2167. * le nouveau champ est fixe on va tester l'equilibre(convergence)
  2168. * et calculer la force motrice pour l'iteration suivante
  2169. *----------------------------------------------------------------------
  2170. 'DETRUIRE' FCORF;
  2171. *
  2172. * Champ de deplacement total
  2173. ZDETOT = ZD0SLX '+' ZDEPT ;
  2174. *
  2175. * XXX1 = [ FCORF ; FCORU ] = [ -Freac_it ; u^imp_it ]
  2176. XXX1 = ZCLIM '*' ZDETOT;
  2177. FCORF = 'ENLE' XXX1 'FLX';
  2178. FCORU = 'EXCO' XXX1 'FLX' 'NOID' 'FLX';
  2179. 'DETR' XXX1 ;
  2180. *----------------------------------------------------------------------*
  2181. * Increment de deformations totales DEPST *
  2182. *----------------------------------------------------------------------*
  2183. * Update or total lagrangian ---------------------------------------
  2184. 'SI' IFEFP;
  2185. 'SI' IFEFPUL ;
  2186. * mess ' update lagrangian ZRIKTA';
  2187. GEOM2 = 'FORM' ZDEPT ;
  2188. ZRIKTA ZSIGF ZVARF ZDEIF = 'ECFEFP'
  2189. ZMODL ZDEI0 ZVAR0 ZDEPT ZMAT2F ZPREK NITMA 1 ;
  2190. 'SINON';
  2191. * mess ' total lagrangian ZRIKTA';
  2192. chp_z = ZDEPT + ZDEP0 ;
  2193. GEOM2 = 'FORM' chp_z ;
  2194. ZRIKTA ZSIGF ZVARF ZDEIF = 'ECFEFP'
  2195. ZMODL ZDEI0 ZVAR0 chp_z ZMAT2F ZPREK NITMA ;
  2196. chp_z = 1 ;
  2197. 'FINS';
  2198. FEQU2 = 'BSIG' ZMODL ZSIGF ZMAT2F ;
  2199. ZRAID = ZRIKTA 'ET' ZCLIM ;
  2200. 'FORM' GEOM1 ;
  2201. *
  2202. XXX1 = 'EXCO' (WTAB.'MOVA') ZVARF ;
  2203. EPSM = 'MAXI' XXX1 'AVEC' MLDEFOR ;
  2204. ACC = 'ABS' ( XXX1 - ACC0 ) ;
  2205. MMC = 'MASQUE' ACC 'SUPERIEUR' 1.D-10 'SOMME' ;
  2206. 'SI' (MMC '>' MMCMAX) ; MMCMAX = MMC ; 'FINS' ;
  2207. DPSMAX = 'MAXI' ACC 'AVEC' MLDEFOR ;
  2208. DEPST = 'CHAN' 'TYPE' ZDEIF 'DEFORMATIONS' ;
  2209. ZDEFF = ZDEF0 '+' DEPST ;
  2210. *
  2211. * cas standard -----------------------------------------------------
  2212. 'SINON';
  2213. 'SI' IGRD;
  2214. zdepm = zdept '*' 0.5;
  2215. 'FORM' GEOM1 ;
  2216. 'SI' WTAB.'ITCAR';
  2217. GEOM2M ZMATTEMP = 'FORM' ZDEPM ZMODLI ZMAT22 ;
  2218. * A cause de FORM, ZMATTEMP n'est plus parallele...
  2219. ZMATTEMP = 'REDU' ZMATTEMP ZMODL ;
  2220. 'SINON';
  2221. GEOM2M = 'FORM' ZDEPM ;
  2222. ZMATTEMP = ZMAT2F;
  2223. 'FINS';
  2224. GEOM2 = 'FORM' zdepm ;
  2225. 'FORM' GEOM1 ;
  2226. 'FINS';
  2227. *
  2228. * Pour les hypotheses de deformations non lineaires,
  2229. * on calcule l'increment de deformation sur le pas avec EPSI 'LINE'
  2230. * mais en se placant sur une configuration intermediaire (mi pas)
  2231. iDEFOK = VRAI ; PASOK = vrai;
  2232. 'SI' (('NEG' HYPDEF 'LINEAIRE') 'ET' IGRD);
  2233.  
  2234. 'SOUC' 0;
  2235. zdepti = zdept * (1.0 / n_sous_dec);
  2236. depst = 0;
  2237. 'REPE' b_eps n_sous_dec;
  2238. teps = (&b_eps * (1. / n_sous_dec)) - (1. / n_sous_dec /2.);
  2239. zdepa = zdept '*' teps;
  2240.  
  2241. geom2a = form zdepa;
  2242. 'SI' ('OU' IENDOM IVIDOM ICERAM);
  2243. HOOKENDO = 'HOOK' ZMODL ZMAT2F ZVAR0 ;
  2244. DEPSTat = 'EPSI' 'LINE' 'NOER' ZMODL ZDEPTi HOOKENDO ZMAT2F;
  2245. 'DETR' HOOKENDO;
  2246. 'SINON';
  2247. DEPSTat = 'EPSI' 'LINE' 'NOER' ZMODL ZDEPti ZMATTEMP ;
  2248. 'FINSI';
  2249. 'SI' ('SOUCI'); 'FORM' GEOM1; IDEFOK = faux; 'QUIT' b_eps; 'FINSI';
  2250. DEPSTa = 'PICA' HYPDEF DEPSTat ZMODL (-1. '*' ZDEPa) ;
  2251. FORM geom1;
  2252. depst = depst + depsta;
  2253. 'DETR' geom2a; 'DETR' zdepa; 'DETR' depsta; 'DETR' depstat;
  2254. 'FIN' b_eps;
  2255.  
  2256.  
  2257.  
  2258.  
  2259. PASOK = 'SOUCI' 'OU' ('NON' IDEFOK);
  2260.  
  2261. DPSMAX = 1 ;
  2262. 'SI' ('NON' PASOK);
  2263. 'FORM' GEOM1 ;
  2264.  
  2265. DPSMAX = 'MAXI' 'ABS' DEPST 'AVEC' MLDEFOR;
  2266. 'SI' (DPSMAX '>' 1000);
  2267. PASOK = VRAI;
  2268. si faux;
  2269. RE_CNTRL2*'RIGIDITE' = DEFO_IMP DEPST (zdept + zdep1) ZRAID;
  2270. DE_CNTRL = FAUX;
  2271. ZFCONT = RE_CNTRL2 '*' ( ZDEP1) 'EXCO' 'FLX' 'NOID' 'FLX' ;
  2272. cmul = MINI 1. ( 1d-3 / DPSMAX);
  2273. ZFEXT2 = ZFEXT2 + (ZFCONT * cmul);
  2274. ZFLXB = 'EXCO' ZFEXT2 'FLX' 'NOID' 'FLX' ;
  2275. ZFCONSTA = ZFEXT2 ;
  2276. RE_CNTRL = RE_CNTRL 'ET' RE_CNTRL2;
  2277. finsi;
  2278. HPP_EPS = FAUX;
  2279. pastest = VRAI;
  2280. urg = VRAI;
  2281. augmult = augmult * 1.4;
  2282. itacc = 4;
  2283. 'ITERER' ETIQ;
  2284. 'FINS';
  2285. 'FINS';
  2286.  
  2287. 'SI' PASOK;
  2288. MSGZ = 'CHAI' 'Deformation non lineaire non calculable' ;
  2289. pastest = vrai;
  2290. iDEFOK = FAUX ;
  2291. 'SI' autaug;
  2292. mess ('CHAI' MSGZ ' : annulation sous-increment');
  2293. iraug = vrai;
  2294. itacc = 4;
  2295. HPP_EPS = FAUX;
  2296. pastest = VRAI;
  2297. urg = VRAI;
  2298. augmult = augmult * 1.4;
  2299. ITACC=4;
  2300. 'ITERER' ETIQ;
  2301. dpsmax=10000;
  2302. 'FINSI';
  2303. MESS MSGZ ;
  2304. 'FINSI';
  2305.  
  2306. 'SI' PASOK;
  2307. 'MESS' 'deformation lineaire';
  2308. 'FINSI';
  2309. 'FINS';
  2310. *
  2311. ** on n'a pas pu calcule DEPST : on utilise les deformations lineaires
  2312. 'SI' PASOK;
  2313. 'SI' ('OU' IENDOM IVIDOM ICERAM);
  2314. HOOKENDO = 'HOOK' ZMODL ZMAT2F ZVAR0 ;
  2315. DEPST = 'EPSI' 'LINE' ZMODL ZDEPT HOOKENDO ZMAT2F;
  2316. 'DETR' HOOKENDO;
  2317. 'SINON';
  2318. DEPST = 'EPSI' 'LINE' ZMODL ZDEPT ZMAT2F;
  2319. 'FINSI';
  2320. DPSMAX = 'MAXI' 'ABS' DEPST 'AVEC' MLDEFOR;
  2321. 'SI' (('NEG' HYPDEF 'LINEAIRE') 'ET' IGRD);
  2322. 'MESS' 'Attention utilisation des deformations lineaires';
  2323. 'FINS';
  2324. 'FINS';
  2325. *
  2326. * Exprimer l'increment de deformations (DEPST) et les contraintes
  2327. * initiales (ZSIG0a) sur la configuration adequate
  2328. ZSIG0a = ZSIG0;
  2329. 'SI' (IGRD 'ET' iDEFOK);
  2330. 'FORM' GEOM1;
  2331. * Lagrangien total : sur GEOREF0
  2332. 'SI' ('EGA' LAG_TOT 1);
  2333. DEPSTa = 'PICA' HYPDEF DEPST ZMODL ZDEP0N ;
  2334. ZSIG0a = 'PICA' HYPDEF ZSIG0 ZMODL ZDEP0N ;
  2335. 'FORM' GEOREF0;
  2336. 'FINSI';
  2337. * Lagrangien fin pas : sur GEOM2
  2338. 'SI' ('EGA' LAG_TOT 2);
  2339. DEPSTa = 'PICA' HYPDEF DEPST ZMODL ZDEPT ;
  2340. ZSIG0a = 'PICA' HYPDEF ZSIG0 ZMODL ZDEPT ;
  2341. 'FORM' GEOM2;
  2342. 'FINS';
  2343. * Lagrangien mi pas : sur GEOM2M
  2344. 'SI' ('EGA' LAG_TOT 3);
  2345. DEPSTa = 'PICA' HYPDEF DEPST ZMODL ZDEPM ;
  2346. ZSIG0a = 'PICA' HYPDEF ZSIG0 ZMODL ZDEPM ;
  2347. 'FORM' GEOM2M;
  2348. 'FINS';
  2349. 'DETR' depst;
  2350. depst = depsta;
  2351. 'FINS';
  2352. *
  2353. 'SI' (knoconv '>' 1) ;
  2354. ZDT = 0. ;
  2355. 'SINON' ;
  2356. ZDT = WTAB.'DT' '*' COEPI ;
  2357. 'FINS' ;
  2358. * Verifier si les deformations obtenues ne sont pas trop
  2359. * importantes
  2360. 'SI' (IGRD 'ET' (('NON' HPP_EPS) 'OU' (IT > 0 )));
  2361. ref = 20;
  2362. n_sous_dec = MAXI 1 (n_sous_dec - 3) (((dpsmax/1d-2) + 1) ENTIER);
  2363. n_sous_dec = MINI 10 n_sous_dec ;
  2364. ** mess 'sous decoupage pour epsi : ' n_sous_dec;
  2365. coefmul = 1.d0 / (DPSMAX + XPETIT);
  2366. 'SI' (((coefmul < ref) 'OU' PASUNIL) 'ET' AUTAUG);
  2367. 'MESS' 'annulation du sous-increment: trop grand ' dpsmax;
  2368. si faux;
  2369. RE_CNTRL2*'RIGIDITE' = DEFO_IMP DEPST (zdept + zdep1) ZRAID ;
  2370. **pv DE_CNTRL = VRAI ;
  2371. ZFCONT = RE_CNTRL2 '*' ( ZDEP1) 'EXCO' 'FLX' 'NOID' 'FLX' ;
  2372. cmul = MINI 1. ( 1d-3 / DPSMAX);
  2373. ZFEXT2 = ZFEXT2 + (ZFCONT * cmul);
  2374. ZFLXB = 'EXCO' ZFEXT2 'FLX' 'NOID' 'FLX' ;
  2375. ZFCONSTA = ZFEXT2 ;
  2376. RE_CNTRL = RE_CNTRL 'ET' RE_CNTRL2;
  2377. finsi;
  2378. * il faut reduire d'urgence les deplacements
  2379. coefmul = 0.0;
  2380. 'SI' ((red_urg > 3) 'ET' ('NON' instab));
  2381. resmul = resmul * ('MINI' 0.25 (1d-2 / dpsmax));
  2382. 'FINSI';
  2383. augmult = augmult * 1.5;
  2384. urg = vrai;
  2385. * on garde zdep1d pour le calcul des coeff d'acceleration
  2386. zdep1d = zdep1 'ENLE' 'LX';
  2387. 'SI' AUTAUG;
  2388. iraug = vrai ;
  2389. *pv itacc = 3;
  2390. 'FINSI';
  2391. 'SI' ((dpsmaxp < 5e-4) 'OU' (dpsmaxp > 2e-2 )) ;
  2392. PASTEST=VRAI;
  2393. 'SINON';
  2394. 'SI' (red_urg > 4); nconvr = vrai; 'FINSI';
  2395. 'FINSI';
  2396. 'SI' (RED_URG > 2 ) ;
  2397. 'SI' ('NON' nonconv);
  2398. 'MESS' ' non convergence detectee 1' ;
  2399. ZMAXIT = IT+5;
  2400. 'FINS';
  2401. nonconv = vrai;
  2402. 'FINS';
  2403. itacc = 4;
  2404. 'SI' (resmul < 0.99);
  2405. resmul = resmul * ('MINI' 0.25 (1d-2 / dpsmax));
  2406. 'FINSI';
  2407. HPP_EPS = FAUX;
  2408. RED_URG = RED_URG + 1;
  2409. 'SI' ((dpsmaxp > 1.8d-2) 'ET' (augmult > 100));
  2410. zprecnc= 2.1d-2;
  2411. 'SI' (zmaxit > it) ;
  2412. zmaxit = it - 1;
  2413. 'FINSI';
  2414. 'FINSI';
  2415. 'SI' ((zmaxit > it) 'ET' (RED_URG > 4 )); ZMAXIT = IT-1; 'FINS';
  2416. 'ITERER' etiq;
  2417. 'SINON';
  2418. ** RE_CNTRL = 'VIDE' 'RIGIDITE' / 'RIGIDITE';
  2419. ** DE_CNTRL = FAUX;
  2420. 'FINS';
  2421. 'FINS';
  2422. *----------------------------------------------------------------------*
  2423. * Loi de comportement (operateur ELAS ou COMP) *
  2424. * sur la configuration definie par LAG_TOT *
  2425. *----------------------------------------------------------------------*
  2426. * Comportement elastique lineaire
  2427. 'SI' ('NON' IPLAVI) ;
  2428. *
  2429. * Increment de contraintes
  2430. DSIGT = 'ELAS' ZMODL DEPST ZMAT2F ;
  2431. *
  2432. * Soustraire les parties thermique et deformations imposees
  2433. 'SI' (ITHER 'OU' WTAB.'MATVAR') ;
  2434. XXX1 = DSIGT0 '*' COEPI;
  2435. XXX2 = DSIGT '-' XXX1;
  2436. DSIGT = XXX2 ;
  2437. 'FINS';
  2438.  
  2439. 'SI' LOGDEF ;
  2440. XXX1 = DSI1 '*' COEPI;
  2441. XXX2 = DSIGT '-' XXX1;
  2442. DSIGT = XXX2 ;
  2443. 'FINS';
  2444. ZSIGF = ZSIG0a + DSIGT ;
  2445. *
  2446. * Comportement non lineaire : integration avec l'operateur COMP
  2447. 'SINON';
  2448. *
  2449. * dans le cas des modeles endommageables de Lemaitre, on ecoule
  2450. * en tenant compte, dans les iterations internes, de la variation du
  2451. * materiau avec la temperature
  2452. *
  2453. * Soustraire les parties thermique et deformations imposees
  2454. 'SI' (ITHER 'OU' WTAB.'MATVAR') ;
  2455. XXX1 = DEPST0 '*' COEPI;
  2456. XXX2 = DEPST '-' XXX1 ;
  2457. DEPST = XXX2 ;
  2458. 'FINS' ;
  2459. *
  2460. 'SI' LOGDEF ;
  2461. XXX1 = DDEFOR0 '*' COEPI;
  2462. XXX2 = DEPST '-' XXX1 ;
  2463. DEPST = XXX2 ;
  2464. 'FINS' ;
  2465. *
  2466. ZMATT = ZMAT2F ;
  2467. 'SI' ('ET' ('ET' ITHER ('OU' IENDOM IVIDOM )) WTAB.'MATVAR') ;
  2468. * on recupere certain materiau avec les parametres fct de la temperature
  2469. * voir PAS_mate (il ne faut que la dependance thermique)
  2470. ZMATT = 'REDU' WTAB.'MA_COMP' ZMODL;
  2471. 'FINS';
  2472. *
  2473. * ...cas SSTE...
  2474. 'SI' ISSTE;
  2475. ZRIKTA ZSIGF ZVARF ZDEIF =
  2476. 'SSTE' ZMODL ZSIG0a ZVAR0 DEPST ZMATT
  2477. ZPREK WTAB.'NMAXSUBSTEPS' NITMA;
  2478. zvar0 = ('EXCO' zvar0 ('ENLE' com_var ('DIME' com_var)))
  2479. 'ET' ('EXCO' zvarf ssii);
  2480. zvar0 = 'CHANGER' 'TYPE' zvar0 'VARIABLES INTERNES';
  2481. *
  2482. * ...cas non SSTE...
  2483. 'SINON';
  2484. *
  2485. * Gradient des deplacements (total) a la fin du pas de temps
  2486. * par rapport a la configuration de reference GEOREF0
  2487. 'SI' IGRD;
  2488. 'SI' ('NEG' ZGRDU0 'INCONNU');
  2489. 'SI' ('NEG' LAG_TOT 1);
  2490. geot = 'FORM';
  2491. 'FORM' GEOREF0 ;
  2492. 'FINS' ;
  2493. ZGRDUF = 'GRAD' ZMODL ZMAT2F ZDETOT ;
  2494. D_GR_U = ZGRDUF '-' ZGRDU0 ;
  2495. 'SI' ('NEG' LAG_TOT 1); 'FORM' GEOT; 'FINS';
  2496. 'FINS';
  2497. 'FINS';
  2498. *
  2499. ZSIG01 = ZSIG0a;
  2500. ZVAR01 = ZVAR0;
  2501. ZDEI01 = ZDEI0;
  2502.  
  2503. 'SI' (nsoincr 'NEG' 1);
  2504. ZSIGM = ZSIG0a '/' 2;
  2505. 'FINS';
  2506. *
  2507. 'REPETER' sousinc nsoincr;
  2508. *
  2509. codeb = 'FLOT' (&sousinc '-' 1) '/' nsoincr ;
  2510. cofin = 'FLOT' &sousinc '/' nsoincr ;
  2511. *
  2512. tsodeb = WTAB.'TEMPS0' '+' (ZDT '*' codeb) ;
  2513. tsofin = WTAB.'TEMPS0' '+' (ZDT '*' cofin) ;
  2514. 'SI' (knoconv '>' 1) ;
  2515. tsodeb = tsofin ;
  2516. 'FINS' ;
  2517. che1 = 'ADET' 'NOUV' ZMODL CHASANST tsodeb 'TEMP' tsodeb;
  2518. che2 = 'ADET' 'NOUV' ZMODL CHASANST tsofin 'TEMP' tsofin;
  2519. *
  2520. 'SI' ('EXIS' WTAB 'MOD_LIA') ;
  2521. vite1 = 'CHAN' 'COMP' conti.'VITESSES' MLPRIM_LIA MVPRIM_LIA ;
  2522. che1 = che1
  2523. 'ET' ('CHAN' 'CHAM' ZMODL vite1 'STRESSES')
  2524. 'ET' ('CHAN' 'CHAM' ZMODL conti.'DEPLACEMENTS' 'STRESSES');
  2525. che2 = che2
  2526. 'ET' ('CHAN' 'CHAM' ZMODL ZDETOT 'STRESSES')
  2527. 'ET' ('CHAN' 'CHAM' ZMODL ZFCONST1 'STRESSES');
  2528. 'FINS';
  2529.  
  2530. 'SI' ither ;
  2531. TETDE = ZTET1 + (DTETD '*' codeb) ;
  2532. TETDF = ZTET1 + (DTETD '*' cofin) ;
  2533. che3 = 'CHAN' 'CHAM' ZMODL TETDE 'STRESSES' ;
  2534. che4 = 'CHAN' 'CHAM' ZMODL TETDF 'STRESSES' ;
  2535. 'FINS' ;
  2536.  
  2537. che1 = che1 'ET' che3 ;
  2538. che2 = che2 'ET' che4 ;
  2539. aa5 = DEPST '*' codeb;
  2540. che5 = ZDEF0 '+' aa5;
  2541. cc5 = DEPST '*' cofin;
  2542. che6 = ZDEF0 '+' CC5;
  2543. *On met les deformations en tete des champs pour COMP ('KTAN' 'PERT')
  2544. che1 = che5 'ET' che1 ;
  2545. che2 = che6 'ET' che2 ;
  2546. * pour les materiaux on garde toujours la valeur a la fin du pas
  2547. * et au debut du pas
  2548. che1 = che1 'ET' ZMAT1 ;
  2549. *
  2550. * Modele NON LINEAIRE UTILISATEUR + GRANDES DEFORMATIONS
  2551. * On ajoute les gradients de deplacements qui seront transformes en
  2552. * gradients de transformation avant appel a UMAT (dans WKUMA1.ESO)
  2553. 'SI' (( 'NEG' ZGRDU0 'INCONNU') 'ET' igrd) ;
  2554. gru1 = ZGRDU0 + (D_GR_U '*' codeb) ;
  2555. che1 = che1 'ET' gru1 ;
  2556. gru2 = ZGRDU0 + (D_GR_U '*' cofin) ;
  2557. che2 = che2 'ET' gru2 ;
  2558. 'FINS';
  2559. *
  2560. 'SI' LNLOC ;
  2561. ISTEP = 1 ;
  2562. chm_z = 'ADET' 'NOUV' ZMODL CHASANST tsodeb 'STEP' ISTEP ;
  2563. Che1a = che1 ;
  2564. che1 = che1 'ET' chm_z ;
  2565. chm_z = 'ADET' 'NOUV' ZMODL CHASANST tsofin 'STEP' ISTEP ;
  2566. Che2a = che2 ;
  2567. che2 = che2 'ET' chm_z ;
  2568. 'FINS';
  2569. *
  2570. che11 = che1 'ET' ZSIG01 'ET' ZVAR01 'ET' ZDEI01 ;
  2571. che22 = che2 'ET' ZMATT ;
  2572. 'SI' PARTLOCA;
  2573. cha11 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che11 ;
  2574. cha22 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che22 ;
  2575. souci 0 ;
  2576. cho22 = 'ASSIS' 'TOUS' 'COMP' MODRELOC cha11 cha22 ;
  2577. isoucomp = souci ;
  2578. cho2 = 'REDU' ZMODL ('ETG' cho22) ;
  2579. cha11 = 0 ; cha22 = 0 ; cho22 = 0 ;
  2580. 'SINON';
  2581. souci 0 ;
  2582. cho2 = 'COMP' ZMODL che11 che22 ;
  2583. isoucomp = souci ;
  2584. 'FINS';
  2585. *
  2586. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  2587. * DEBUT MODELE NON LOCAL
  2588. 'SI' LNLOC ;
  2589. ZVARF = 'CHANGER' ('EXCO' cho2 com_var 'NOID') 'TYPE' 'VARIABLES INTERNES' ;
  2590. * NLOC ne traitant pas des champs //, on reduit tout au modele initial.
  2591. 'SI' (PARALLEL 'ET' ('NON' PARTLOCA)) ;
  2592. ZVARF = 'REDU' ZVARF ZMODLI ;
  2593. 'FINS' ;
  2594.  
  2595. ********** boucle sous-iterations Helmholtz **************
  2596.  
  2597. * donnees pour acceleration Helmholtz
  2598. * periode de recours a l acceleration ACT3 pour Helmholtz
  2599. ZNACCEHL = 3;
  2600. * NBR min de sous-iterations d Helmholtz a faire avant d appeler ACT3 pour la 1ere fois
  2601. ITDEPHL = 4;
  2602.  
  2603. 'REPE' BIHL WTAB.'MAXSOUSITERATION' ;
  2604.  
  2605. 'SI' ('NEG' WTAB.'NON_LOCAL' 'HELM');
  2606. * cas non-local MOYE ou SB
  2607. 'SI' ('EGA' WTAB.'NON_LOCAL' 'SB') ;
  2608. MOD_SB = WTAB.'NLOC_MODL' ;
  2609. CONTP = 'PRIN' ('REDU' ZSIG01 MOD_SB) MOD_SB ;
  2610. ZVARF = ZVARF '+' CONTP '+' WTAB.'NLOC_SB_REGU' ;
  2611. 'FINS' ;
  2612. ZVARN = 'NLOC' ZVARF WTAB.'CONN' ;
  2613. ZVARN = 'CHANGER' ('EXCO' ZVARN com_var) 'TYPE' 'VARIABLES INTERNES' ;
  2614. LIERR1='PROG' 0 ;
  2615. LOGHLIN = VRAI ;
  2616. 'SINON' ;
  2617.  
  2618. * cas non-local HELM Sellier Millard
  2619. MOD_HELM = WTAB.'NLOC_MODL' ;
  2620. ZVARN = 'COPI' ZVARF ;
  2621. ZVARF1 = 'REDU' ZVARN MOD_HELM ;
  2622. LICOHELM = 'EXTR' ZVARF1 'COMP' ;
  2623. * (re) lecture du materiau si init ou evolutif
  2624. ZMATHL = PAS_MATE PRECED cho2 ;
  2625. * (re) calcul des matrices de rigidite de Helmholtz Sellier
  2626. PAS_HELM PRECED ZMATHL;
  2627. * variable pour la liste des residus de convergence sur Helmholtz
  2628. LIERR1='PROG';
  2629.  
  2630. * ----- debut boucle sur les formulations de Helmholtz ---
  2631. * on reste en lineaire si toutes les vari sont lineaires (sinon on aura des istep a 3)
  2632. LOGHLIN = VRAI ;
  2633. 'REPE' BH NHELM ;
  2634. LOGHLIN=LOGHLIN 'ET' TAHELM . &BH. 'LINEAIRE' ;
  2635. * nom de la variable a delocaliser
  2636. LEMO = TAHELM . &BH. 'NOM' ;
  2637. * recuperation de la variable sur le modele mecanique
  2638. ZVAUX = 'EXCO' ZVARF1 LEMO 'SCAL';
  2639. ZVAUX2 = 'CHAN' 'NOEUD' MOD_HELM ZVAUX;
  2640. * projection sur le modele de Helmholtz associe
  2641. ZVAUX2 = 'PROI' TAHELM . &BH . 'H_MODELE' ZVAUX2 'MINI';
  2642. * TRAC zvaux2 TAHELM . &BH . 'H_MODELE' titre ('CHAI' 'Vari a diffuser:'LEMO );
  2643. * la variable devient la source
  2644. FSOUR = 'SOUR' (TAHELM . &BH . 'H_MODELE') ZVAUX2 ;
  2645. * trac fsour (extr TAHELM . &BH . 'H_MODELE' 'MAIL') titre 'Source';
  2646. * trac (TAHELM . &BH . 'H_OPER') titre 'H_Oper';
  2647. * *** pour istep=1 *** 1ere etape non locale ***
  2648. 'SI' ('EGA' &BIHL 1) ;
  2649. * &BIHL est le compteur de sous iteration de Helmholtz
  2650. * il vaut forcement 1 si istep vaut 1, on fait la resolution classique
  2651. ZVNEW = 'RESO' (TAHELM . &BH . 'H_OPER') FSOUR ;
  2652. * sauvegarde pour possible acceleration ulterieure
  2653. TAHELM . &BH . 'XN' = ZVNEW ;
  2654. 'SINON';
  2655. * recuperation de la solution de la sous iteration de Helmholtz precedente
  2656. ZVNEW = TAHELM . &BH . 'XN' ;
  2657. 'FINSI';
  2658. * *** traitement non lineaire de la formulation ***
  2659. 'SI' ( 'NON' TAHELM . &BH . 'LINEAIRE' ) ;
  2660. * calcul du residu sur l evolution de la source
  2661. RESN = FSOUR - (TAHELM . &BH . 'H_OPER' '*' TAHELM . &BH . 'XN' );
  2662. * stockage du residu pour la mesure d erreur
  2663. RESN0='NOMC' RESN 'SCAL' ;
  2664. * normalisation des residus par l inverse de la capacite
  2665. * (calculee dans pas_helm.procedur)
  2666. ZERR0 = RESN0 '*' (TAHELM . &BH . 'INV_CAPA' );
  2667. * erreur residu
  2668. ZERR2 = 'MAXI' ( 'ABS' (ZERR0) );
  2669. * 'SI' ('EGA' &BIHL 1);
  2670. * TAHELM . &BH . 'RES1' = ZERR2;
  2671. * 'FINSI';
  2672. * carre de norme du residu
  2673. * ZERR1 = 'XTX' ZERR0;
  2674. * Residu
  2675. * ZERR2 = (ZERR1)**0.5;
  2676. * normalisation de la source par la capacite
  2677. ZERR3 = FSOUR '*' (TAHELM . &BH . 'INV_CAPA' )('MOTS' 'Q') ('MOTS' 'SCAL') ('MOTS' 'SCAL');
  2678. * calcul erreur
  2679. ZERR6 = 'MAXI' ( 'ABS' (ZERR3) );
  2680. * carre de la norme de la source
  2681. * ZERR5 = 'XTX' ZERR3;
  2682. * norme de la source
  2683. * ZERR6 = ZERR5 **0.5;
  2684. * valeur maximale du champs
  2685. * ZERR3='MAXI' (TAHELM . &BH . 'XN') 'ABS';
  2686. * mess 'variable Helmholtz', LEMO;
  2687. * mess 'residu', ZERR2;
  2688. * mess 'norme de la source ', ZERR5;
  2689. * mess 'norme Source', ZERR3;
  2690. * on teste si la valeur maximale est capable de normaliser
  2691. 'SI' (ZERR6 '>' 0.);
  2692. ZERR4 = ZERR2 '/' ZERR6;
  2693. 'SINON';
  2694. ZERR4 = ZERR2;
  2695. 'FINS';
  2696. * preparation acceleration
  2697. 'SI' ( 'EGA' &BIHL 1 ) ;
  2698. TAHELM . &BH . 'H_RESO' = TAHELM . &BH . 'H_OPER' ;
  2699. FWOR = FSOUR * 0. ;
  2700. TAHELM . &BH . 'ACFP1' = FWOR ;
  2701. TAHELM . &BH . 'ACFP2' = FWOR ;
  2702. TAHELM . &BH . 'RESNP' = FWOR ;
  2703. TAHELM . &BH . 'ACFP3' = FWOR ;
  2704. TAHELM . &BH . 'ACFEP1' = FWOR ;
  2705. TAHELM . &BH . 'ACFEP2' = FWOR ;
  2706. TAHELM . &BH . 'CORREC' = 0.;
  2707. 'SINON' ;
  2708. * acceleration de convergence act3
  2709. CORRECPHL = TAHELM . &BH . 'CORREC';
  2710. CORRECHL = 0;
  2711. ACFP0HL = RESN - TAHELM . &BH . 'RESNP' ;
  2712. ACFEP0HL = ACFP0HL - CORRECPHL ;
  2713. 'SI' ('MULT' &BIHL ZNACCEHL) ;
  2714. 'SI' (&BIHL > ITDEPHL);
  2715. CORRECHL = 'ACT3'
  2716. TAHELM . &BH . 'ACFEP2'
  2717. TAHELM . &BH . 'ACFEP1' ACFEP0HL
  2718. TAHELM . &BH . 'ACFP3'
  2719. TAHELM . &BH . 'ACFP2'
  2720. TAHELM . &BH . 'ACFP1' ACFP0HL ;
  2721. RESN = RESN - CORRECHL ;
  2722. 'FINS';
  2723. 'FINS';
  2724. TAHELM . &BH . 'RESNP' = RESN ;
  2725. TAHELM . &BH . 'ACFP3' = TAHELM . &BH . 'ACFP2' ;
  2726. TAHELM . &BH . 'ACFP2' = TAHELM . &BH . 'ACFP1' ;
  2727. TAHELM . &BH . 'ACFP1' = ACFP0HL ;
  2728. TAHELM . &BH . 'ACFEP2' = TAHELM . &BH . 'ACFEP1' ;
  2729. TAHELM . &BH . 'ACFEP1' = ACFEP0HL ;
  2730. * calcul correction delta-sigma
  2731. DZVN = 'RESO' TAHELM . &BH . 'H_RESO' RESN;
  2732. ZVNEW = ZVNEW + DZVN ;
  2733. TAHELM . &BH . 'XN' = ZVNEW ;
  2734. 'FINSI' ;
  2735. 'SINON';
  2736. * cette variable est traitee de facon lineaire
  2737. ZERR4=0.;
  2738. 'FINSI';
  2739. ZVNEW1 = 'NOMC' ZVNEW 'SCAL' ;
  2740. ZVNEW = 'CHAN' 'CHAM' ZVNEW1 MOD_HELM 'STRESSES' 'VARIABLES INTERNES';
  2741. * trac zvnew mod_helm titre 'ZVNEW';
  2742. DZVN = ZVNEW '-' ZVAUX ;
  2743. DZVN = 'NOMC' DZVN LEMO ;
  2744. ZVARN = ZVARN '+' DZVN ;
  2745.  
  2746. * mess 'Erreur retenue', ZERR4;
  2747. * stockage de l erreur normalisee dans la liste des erreurs
  2748. LIERR1 = LIERR1 'ET' ('PROG' ZERR4);
  2749.  
  2750. * cas particulier TWL2 Sellier theorie WL2
  2751. 'SI' ('EGA' LEMO 'TWL2') ;
  2752. ZSMAX = 'EXCO' ZVARF1 'SMAX' 'SMAX' ;
  2753. ZSMAN = 'MANU' 'CHML' MOD_HELM 'SMAX' ('MAXI' ZSMAX) 'TYPE' 'VARIABLES INTERNES' 'STRESSES' ;
  2754. ZVARN = ZVARN '+' ZSMAN '-' ZSMAX ;
  2755. 'FINS' ;
  2756. * ---- cas particulier SER Sellier renfort avec gradient
  2757. * calcul GRADIENT des contraintes sur les renforts (SERi)
  2758. * stockage dans la VARI HELM1(&BH,2) du gradient selon l absc curviligne
  2759. MOPRE = 'EXTR' LEMO 1 3 ;
  2760. 'SI' ('EGA' MOPRE 'SHR');
  2761. NUME='EXTR' LEMO 4;
  2762. NDIM1='VALE' 'DIME';
  2763. 'SI' ( 'EGA' NDIM1 3 );
  2764. * valeur SERi du pas precedent dans Helm(&bh,1)
  2765. LEMO1=CHAI 'H' &BH '1';
  2766. * valeur SERi du pas actuel
  2767. ZAUX1='EXCO' ZVARN LEMO 'SCAL';
  2768. ZVNEW1='NOMC' ZAUX1 'T' ;
  2769. ZVAUX1 = 'CHAN' 'NOEUD' MOD_HELM ZVNEW1;
  2770. ZVNEW2='PROI' TAHELM . &BH . 'H_MODELE' ZVAUX1 'MINI';
  2771. GRA1=GRAD TAHELM . &BH . 'H_MODELE' ZVNEW2;
  2772. GRA2=CHAN 'CHAM' 'NOEUD' GRA1 TAHELM . &BH . 'H_MODELE';
  2773. GRA3=PROI GRA2 MOD_HELM;
  2774. LMO3=EXTR GRA3 COMP;
  2775. NOM1=EXTR LMO3 1;
  2776. NOM2=EXTR LMO3 2;
  2777. NOM3=EXTR LMO3 3;
  2778. GRA11=EXCO NOM1 GRA3 'SCAL';
  2779. GRA22=EXCO NOM2 GRA3 'SCAL';
  2780. GRA33=EXCO NOM3 GRA3 'SCAL';
  2781. * recuperation vecteur direction de projection
  2782. LEMON1=CHAI 'VR' NUME '1';
  2783. NX1=EXCO LEMON1 ZMATHL 'SCAL';
  2784. NX1=CHAN 'NOEUD' NX1 MOD_HELM;
  2785. NX11=CHAN NX1 TYPE 'SCALAIRE';
  2786. LEMON2=CHAI 'VR' NUME '2';
  2787. NX2=EXCO LEMON2 ZMATHL 'SCAL';
  2788. NX2=CHAN 'NOEUD' NX2 MOD_HELM;
  2789. NX22=CHAN NX2 TYPE 'SCALAIRE';
  2790. LEMON3=CHAI 'VR' NUME '3';
  2791. NX3=EXCO LEMON3 ZMATHL 'SCAL';
  2792. NX3=CHAN 'NOEUD' NX3 MOD_HELM;
  2793. NX33=CHAN NX3 TYPE 'SCALAIRE';
  2794. * produit scalaire
  2795. ZVAUX3=(NX11*GRA11)+(NX22*GRA22)+(NX33*GRA33);
  2796. * stockage du gradient de epse dans YU2j
  2797. mot4=CHAI 'H' NUME '2';
  2798. zvaux5= 'NOMC' zvaux3 mot4;
  2799. zvaux6= 'CHAN' 'CHAM' zvaux5 MOD_HELM 'STRESSES' 'VARIABLES INTERNES';
  2800. * trac zvaux6 MOD_HELM titre '3D H2';
  2801. ZVAUX7 = 'EXCO' ZVARF1 MOT4 ;
  2802. DZVN = ZVAUX6 '-' ZVAUX7;
  2803. ZVARN = ZVARN '+' DZVN ;
  2804. 'SINON';
  2805. * valeur SERi du pas precedent dans Helm(&bh,1)
  2806. LEMO1=CHAI 'H' &BH '1';
  2807. ZAUX1='EXCO' ZVARF1 LEMO1 'SCAL';
  2808. ZVNEW1='NOMC' ZAUX1 'T' ;
  2809. ZVAUX1 = 'CHAN' 'NOEUD' MOD_HELM ZVNEW1;
  2810. ZVNEW2='PROI' TAHELM . &BH . 'H_MODELE' ZVAUX1 'MINI';
  2811. GRA1=GRAD TAHELM . &BH . 'H_MODELE' ZVNEW2;
  2812. GRA2=