Télécharger unpas.procedur

Retour à la liste

Numérotation des lignes :

  1. * UNPAS PROCEDUR FD218221 17/10/18 21:15:10 9590
  2. 'DEBPROC' UNPAS PRECED*'TABLE';
  3. *
  4. *----------------------------------------------------------------------*
  5. * PROCEDURE UNPAS *
  6. * *
  7. * calcul d'un increment de solution en grand deplacement plastique *
  8. * par la methode des residus *
  9. *----------------------------------------------------------------------*
  10. * En entree
  11. * PRECED la table passee à PASAPAS
  12. *
  13. * kich nota champ de materiau : etablit au debut de procedure UNPAS
  14. * estimation des caracteristiques fin de pas (TI ): ZMAT22 -->> ZMATXX
  15. * caracteristiques debut de pas (TEMPS0) : ZMAT11
  16. * puis dans la boucle de non-convergence
  17. * ZMAT1 initialise avec ZMAT11
  18. * materiau fin de pas de temps : ZMAT2 sorti de COMP,
  19. * range dans WTAB.'MAT1' en sortie unpas
  20. ************************************************************************
  21. * sortie STAB12 indice : *
  22. * *
  23. * DEPT increment de deplacement sur le pas *
  24. * SIGF contraintes a la fin du pas *
  25. * VARF variables internes a la fin du pas *
  26. * DFPF deformation inelastique a la fin du pas *
  27. * CONV logique valant vrai si pas de probleme de convergence *
  28. * DEFF deformations a la fin du pas si grandes deformations *
  29. *----------------------------------------------------------------------*
  30. *
  31.  
  32. *CB215821 : Recuperation de XPETIT (07/12/2016)
  33. XPETIT = 'VALE' 'PETI' ;
  34.  
  35. WTAB = PRECED.'WTABLE' ;
  36. LAG_TOT=VRAI;
  37. 'SI' ('EXIS' WTAB 'LAG_TOT'); LAG_TOT=WTAB.'LAG_TOT'; 'FINS';
  38. conti= PRECED.'CONTINUATION';
  39. estim= PRECED.'ESTIMATION' ;
  40. *
  41. * initialisation et reprise des valeurs des tables ****************
  42. *
  43. MXMYMZ = 'MOTS' 'MX' 'MY' 'MZ' 'MT' 'FP' 'FPQ' 'FTP'
  44. 'IMX' 'IMY' 'IMZ' 'IMT' 'IFP' 'IFPQ' 'IFTP' ;
  45. MXMFLX = 'MOTS' 'MX' 'MY' 'MZ' 'MT' 'FLX' 'FP' 'FPQ' 'FTP'
  46. 'IMX' 'IMY' 'IMZ' 'IMT' 'FLX' 'IFP' 'IFPQ' 'IFTP' ;
  47. MLPRIM = 'MOTS' 'UX' 'UY' 'UZ' 'UR' 'UT' 'RX' 'RY' 'RZ' 'RT'
  48. 'LX' 'P' 'PQ' 'TP' 'ALFA' 'BETA'
  49. 'IUX' 'IUY' 'IUZ' 'IUR' 'IUT' 'IRX' 'IRY' 'IRZ' 'IRT'
  50. 'IP' 'IPQ' 'ITP' 'IALF' 'IBET' ;
  51. MLDUAL = 'MOTS' 'FX' 'FY' 'FZ' 'FR' 'FT' 'MX' 'MY' 'MZ' 'MT'
  52. 'FLX' 'FP' 'FPQ' 'FTP' 'FALF' 'FBET'
  53. 'IFX' 'IFY' 'IFZ' 'IFR' 'IFT' 'IMX' 'IMY' 'IMZ' 'IMT'
  54. 'IFP' 'IFPQ' 'IFTP' 'IFAL' 'IFBE' ;
  55. MLDEPL = 'MOTS' 'UX' 'UY' 'UZ' 'UR' 'UT' 'IUX' 'IUY' 'IUZ' 'IUR' 'IUT'
  56. 'ALFA' 'BETA' 'IALF' 'IBET' ;
  57. MLROTA = 'MOTS' 'RX' 'RY' 'RZ' 'RT' 'P' 'PQ' 'TP'
  58. 'IRX' 'IRY' 'IRZ' 'IRT' 'IP' 'IPQ' 'ITP';
  59. MLDEFOR = 'MOTS' 'EPXX' 'EPYY' 'EPZZ' 'EPSS' 'EPTT' 'EPRR'
  60. 'GAXY' 'GAXZ' 'GAYZ' 'GAST' 'GASN' 'GATN'
  61. 'GARZ' 'GART' 'GAZT' 'RTSS' 'RTTT' 'RTST'
  62. 'RTZZ' 'RTXX' 'EPS ' 'GXY ' 'CX ' 'CY ' 'CZ '
  63. 'EPSE' ;
  64. MLDEFOR = 'MOTS' 'EPXX' 'EPYY' 'EPZZ' 'EPSS' 'EPTT' 'EPRR'
  65. 'GAXY' 'GAXZ' 'GAYZ' 'GAST' 'GASN' 'GATN'
  66. 'GARZ' 'GART' 'GAZT' 'GXY ' 'CX ' 'CY ' 'CZ '
  67. 'EPSE' ;
  68. MVPRIM = 'MOTS' 'VTX' 'VTY' 'VTZ' 'VTR' 'VTT' 'VWX' 'VWY' 'VWZ' 'VWT'
  69. 'VLX' 'VVP' 'VVPQ' 'VVTP' 'VALF' 'VBET'
  70. 'IVTX' 'IVTY' 'IVTZ' 'IVTR' 'IVTT' 'IVWX' 'IVWY' 'IVWZ' 'IVWT'
  71. 'IVVP' 'IVPQ' 'IVTP' 'IVAL' 'IVBE' ;
  72. MOCA = 'MOTS' 'VECT' 'VX ' 'VY ' 'VZ ' 'VXF ' 'VYF ' 'VZF '
  73. 'V1X ' 'V1Y ' 'V1Z ' 'V2X ' 'V2Y ' 'V2Z ' ;
  74.  
  75. * definition de variables locales
  76. *'SI' ( 'EXIS' PRECED 'ECRIT' ) ; list wtab; 'FINS';
  77. ICERAM = WTAB.'CERAMIQUE' ;
  78. IDYN = WTAB.'DYNAMIQUE';
  79. IELANL = WTAB.'NON_LINEAIRE';
  80. IENDOM = WTAB.'ENDOMMAGEMENT';
  81. IFEFP = WTAB.'FEFP_FORMULATION' ;
  82. IFEFPUL= WTAB.'UPDATE_LAGRANGIAN';
  83. IFTOL = 'NEG' WTAB.'FTOL' 'INCONNU' ;
  84. IGRD = WTAB.'GRANDS_DEPLACEMENTS';
  85. IKSIA = WTAB.'K_SIGMA';
  86. IKTAN = WTAB.'K_TANGENT' ;
  87. IMPLP = WTAB.'LIAISON_PERSISTANTE';
  88. IMTOL = 'NEG' WTAB.'MTOL' 'INCONNU';
  89. IPILOT = WTAB.'AUTOMATIQUE';
  90. IPLAST = WTAB.'PLASTIQUE';
  91. IPLAVI = WTAB.'IPLAVI';
  92. IPREDIC = WTAB.'PREDICTEUR';
  93. IRCON = WTAB.'RAIDCONST';
  94. IRAUG = WTAB.'RAIDAUGM';
  95. ISOL = WTAB.'CONSOLIDATION';
  96. ISSTE = WTAB.'SUBSTEPPING';
  97. ITHER = WTAB.'CHAR_THE' 'OU' WTAB.'FOR_THER';
  98. IVIEXT = WTAB.'VISCO_EXTERNE';
  99. IVIDOM = WTAB.'VISCODOMMAGE';
  100. IVISCO = WTAB.'VISCOPLASTIQUE';
  101. LNLOC = WTAB.'NLOC';
  102. LOGDEF = WTAB.'CHAR_DEFI';
  103. LOGPRE = WTAB.'CHAR_PRES' ;
  104. NITMA = WTAB.'NITERINTER_MAX';
  105. NSSTE = WTAB.'NMAXSUBSTEPS';
  106. POR1 = WTAB.'POR1' ;
  107. TI = WTAB.'T_FINAL';
  108. EKREAC = WTAB.'REAC_GRANDS';
  109. ZMAXIT = WTAB.'MAXITERATION' ;
  110. ZNACCE = 2 ;
  111. ZNCONS = WTAB.'NITER_KTANGENT' ;
  112. ZPREC = WTAB.'PRECISION' ;
  113. ZPREK = WTAB.'PRECISINTER' ;
  114. ZPRECD = WTAB.'PRECDECHARGE' ;
  115. ZPRECM = WTAB.'PRECFLEX' ;
  116. ZCLIM0 = WTAB.'BLOCAGES_MECANIQUES';
  117. *-- On initialise ZCLIM qui va contenir l'ensemble des C.L.
  118. ZCLIM = ZCLIM0 ;
  119. DT = WTAB.'DT' ;
  120. TEMPS0=WTAB.'TEMPS0';
  121. *--- doit on reactualiser la geometrie temp0 ne 0 et grand depl ? ---
  122. 'SI' WTAB.'RECALCUL' ;
  123. WTAB . 'RECARI' = VRAI ;
  124. WTAB . 'RECADET' = VRAI ;
  125. WTAB . 'REA_GEOM' = VRAI ;
  126. * on suppose que l'on est sur la bonne configuration
  127. GEOM1 = 'FORM';
  128. * 'FORM' GEOM1;
  129. 'SINON';
  130. GEOM1 = WTAB.'FOR0' ;
  131. WTAB . 'RECARI' = FAUX ;
  132. WTAB . 'RECADET' = FAUX ;
  133. 'FINS';
  134.  
  135. 'SI' ('OU' ('OU' IENDOM IVIDOM) ICERAM);
  136. WTAB.'RECARI'= VRAI;
  137. 'FINS';
  138.  
  139. 'SI' IRAUG;
  140. RIG_AUG = WTAB.'RIGIDITE_AUGMENTEE';
  141. 'FINS';
  142.  
  143. 'SI' IRCON;
  144. RIG_CONS = WTAB.'RIGIDITE_CONSTANTE';
  145. MAI_CONS ='EXTR' RIG_CONS 'MAIL' ;
  146. 'FINS';
  147.  
  148. *---------- chamelem etat estime pour la fin du pas de temps ------
  149. ZETAT2 = PAS_ETAT PRECED TI ;
  150. ZMAT22 = PAS_MATE PRECED ZETAT2;
  151. ZMATFI = ZMAT22 ;
  152.  
  153. *----------------------- nouveau chargement ----------------------
  154. 'SI' ('EXIS' (WTAB.'CHARGEMENT') 'MECA');
  155. ZFEXT2 = 'TIRE' WTAB.'CHARGEMENT' 'MECA' TI;
  156. TYP_2 = 'TYPE' ZFEXT2;
  157. 'SI' ('NEG' TYP_2 'CHPOINT ');
  158. 'MESS' '*** Erreur dans la definition du chargement ( MECA ) ***';
  159. 'ERRE' '*** Le type du champ n est pas CHPOINT ***';
  160. 'FINS';
  161. 'SINON';
  162. ZFEXT2= 'VIDE' 'CHPOINT' / 'DISCRET';
  163. 'FINS';
  164.  
  165. 'SI' ('EXIS' (WTAB.'CHARGEMENT') 'FORC');
  166. F2_FOR = 'TIRE' WTAB.'CHARGEMENT' 'FORC' TI;
  167. TYP_2 = 'TYPE' F2_FOR;
  168. 'SI' ('NEG' TYP_2 'CHPOINT ');
  169. 'MESS' '*** Erreur dans la definition du chargement ( FORC ) ***';
  170. 'ERRE' '*** Le type du champ n est pas CHPOINT ***';
  171. 'FINS';
  172. ZFEXT2 = ZFEXT2 '+' F2_FOR ;
  173. 'FINS';
  174.  
  175. 'SI' (LOGPRE 'ET' ('NON' IGRD)) ;
  176. MOP = WTAB.'MOD_PRE' ;
  177. ZPEXT ='TIRE' WTAB.'CHARGEMENT' 'PRES' TI ;
  178. 'SI' ('EXIS' WTAB 'MAT_PRE') ;
  179. ZFPEXT = 'BSIG' MOP ZPEXT ('REDU' ZMATFI MOP) ;
  180. 'SINON' ;
  181. ZFPEXT = 'BSIG' MOP ZPEXT ;
  182. 'FINS' ;
  183. ZFEXT2 = ZFEXT2 '+' ZFPEXT ;
  184. 'FINS' ;
  185.  
  186.  
  187. *--------------- Si il existe des deplacements imposes --------------
  188. 'SI' ('EXIS' (WTAB.'CHARGEMENT') 'DIMP');
  189. F2_DEP = 'TIRE' WTAB.'CHARGEMENT' 'DIMP' TI;
  190. ZFEXT2 = ZFEXT2 + F2_DEP;
  191. 'FINS';
  192. *---------------- si chargement deformation actualisation DEFOR ------
  193. 'SI' WTAB.'CHAR_DEFI' ;
  194. WTAB.'DEFOR2' = 'TIRE' WTAB.'CHARGEMENT' 'DEFI' TI;
  195. 'FINS';
  196.  
  197. *---------- dynamique : preparation du second membre --------------
  198. 'SI' ( WTAB.'DYNAMIQUE' ) ;
  199. 'SI' ('EGA' WTAB.'FREA1' 'INCONNU');
  200. 'SI' ('EXIS' WTAB.'CHARGEMENT' 'MECA');
  201. F1 = 'TIRE' WTAB.'CHARGEMENT' 'MECA' TEMP0;
  202. 'FINS';
  203. 'SI' ('EXIS' WTAB.'CHARGEMENT' 'FORC');
  204. F1F = 'TIRE' WTAB.'CHARGEMENT' 'FORC' TEMP0;
  205. 'SI' ('EGA' ('TYPE' F1) 'CHPOINT ') ;
  206. F1 = F1 + F1F ;
  207. 'SINON';
  208. F1 = F1F ;
  209. 'FINS' ;
  210. 'FINS';
  211.  
  212. 'SI' WTAB.'PROCEDURE_CHARMECA';
  213. TFF1 = CHARMECA PRECED TEMPS0 ;
  214. 'SI' ('EXIS' TFF1 'ADDI_SECOND');
  215. FF1=TFF1 .'ADDI_SECOND';
  216. 'FINS';
  217. 'SI' ('EGA' ('TYPE' F1) 'CHPOINT ');
  218. F1=F1 + FF1;
  219. 'SINON';
  220. F1 = FF1;
  221. 'FINS';
  222. 'FINS';
  223.  
  224. 'SI' LOGPRE ;
  225. MOP = WTAB.'MOD_PRE' ;
  226. ZPEXT0 = 'TIRE' WTAB.'CHARGEMENT' 'PRES' TEMPS0 ;
  227. 'SI' ('EXIS' WTAB 'MAT_PRE') ;
  228. FF1 = 'BSIG' MOP ZPEXT0 ('REDU' WTAB.'MAT1' MOP) ;
  229. 'SINON' ;
  230. FF1 = 'BSIG' MOP ZPEXT0 ;
  231. 'FINS' ;
  232. 'SI' ('EGA' ('TYPE' F1) 'CHPOINT ');
  233. F1=F1 + FF1;
  234. 'SINON';
  235. F1 = FF1;
  236. 'FINS';
  237. 'FINS';
  238. *
  239. LAF0 = 'BSIG' WTAB.'MO_TOT' conti.'CONTRAINTES' ZMAT22;
  240. 'SI' IRCON ;
  241. LAF0 = LAF0 'ET'
  242. ('REDU' MAI_CONS (RIG_CONS '*' conti.'DEPLACEMENTS'));
  243. 'FINS';
  244. 'SI' (('EXIS' WTAB.'CHARGEMENT' 'MECA') 'OU'
  245. ('EXIS' WTAB.'CHARGEMENT' 'FORC') 'OU'
  246. wtab.'PROCEDURE_CHARMECA');
  247. LAF1 = F1 - LAF0 ;
  248. 'SINON';
  249. LAF1 = -1 * LAF0;
  250. 'FINS';
  251. LAF2 = 'ENLEVER' (ZCLIM0 '*' conti.'DEPLACEMENTS') 'FLX';
  252. * forces exterieures + reactions - forces interieures au debut du calcul
  253. * c.a.d. (masse*acceleration initiale)+(amortissement*vitesse initiale)
  254. WTAB.'FREA1' = LAF1 - LAF2 ;
  255. 'FINS';
  256. 'SI' wtab.'LIAISON_PERSISTANTE' ;
  257. * forces d'acceleration au debut du pas
  258. 'SI' ('NEG' WTAB.'AMORTISSEMENT' 'INCONNU');
  259. FF4 = WTAB.'AMORTISSEMENT'* conti.'VITESSES';
  260. WTAB.'FMAN'= WTAB.'FREA1' - FF4 ;
  261. 'SINON' ;
  262. WTAB.'FMAN'= WTAB.'FREA1' ;
  263. 'FINS';
  264. 'FINS';
  265. *------------ il faut calculer la matrice de masse tout de suite -------------
  266. 'SI' ( ('NON' ('EXIS' WTAB MASSE )) 'OU'
  267. WTAB.'GRANDS_DEPLACEMENTS');
  268. WTAB.'MASSE' = 'MASS' WTAB.'MO_TOT' ('CHAN' 'MASSE'
  269. ZMAT22 WTAB.'MO_TOT' ) ;
  270. 'SI' WTAB.'MASSCONST';
  271. WTAB.'MASSE'=WTAB.'MASSE' 'ET' WTAB.'MASSE_CONSTANTE';
  272. 'FINS';
  273. 'FINS';
  274. FF = WTAB.'MASSE' *conti.'VITESSES';
  275. FF4 = 4. / DT * FF; 'DETR' FF;
  276. * partie du second membre qui ne depend que des informations du pas prec
  277. WTAB.'FREA1' = FF4 + WTAB.'FREA1';
  278. 'FINS';
  279. *--------- consolidation : preparation du second membre ---------
  280. 'SI' WTAB.'CONSOLIDATION' ;
  281. FF = BSIGMA WTAB.'MOD_POR' conti.'CONTRAINTES' ;
  282. FF4 = 'EXCO' WTAB.'MOT_POR' FF
  283. WTAB.'MOT_POR' 'NOID' 'NATURE' 'DISCRET' ;
  284. 'SI' WTAB.'DYNAMIQUE';
  285. WTAB.'FREA1' = WTAB.'FREA1' + FF4 ;
  286. 'SINON';
  287. WTAB.'FREA1' = FF4 ;
  288. 'FINS';
  289. ZFEXT2 = ZFEXT2 '+' FF4 ; 'DETR' FF ;
  290. * ---- traitement des flux si besoin ----
  291. 'SI' ( 'EXIS' WTAB.'CHARGEMENT' 'FLUX' ) ;
  292. FLUXT0= 'TIRE' WTAB.'CHARGEMENT' 'FLUX' TEMPS0;
  293. FLUXTI= 'TIRE' WTAB.'CHARGEMENT' 'FLUX' TI ;
  294. FACFLU = -1. * WTAB.'DT';
  295. FLUXT = ( FACFLU * (1 - WTAB.'TETA') * FLUXT0 )
  296. + ( FACFLU * WTAB.'TETA' *FLUXTI ) ;
  297. ZFEXT2 = ZFEXT2 '+' FLUXT ;
  298. 'DETR' FLUXT ; 'DETR' FLUXT0; 'DETR' FLUXTI;
  299. 'FINS' ;
  300. 'FINS';
  301. *----------------- calcul de la masse si frequentiel************
  302. 'SI' ( WTAB.'FREQUENTIEL' 'ET' ('NON' ('EXIS' WTAB 'MASSE' )));
  303. WTAB.'MASSE' = 'MASS' WTAB.'MO_TOT'('CHAN' 'MASSE'
  304. ZMAT22 WTAB.'MO_TOT' ) ;
  305. 'FINS';
  306.  
  307. *- Second membre
  308. ZFCONSTA = ZFEXT2 ;
  309.  
  310. ************************************************************************
  311. * Parallélisation du GIBIANE via les ASSISTANTS *
  312. ************************************************************************
  313. ZMODLI = WTAB.'MO_TOT' ;
  314. ZMODLP = WTAB.'MO_TOT_PREC';
  315. NBPART = WTAB.'NBPART' ;
  316.  
  317. PARALLEL = FAUX ;
  318. PARTLOCA = FAUX ;
  319. ZMODL = ZMODLI ;
  320. 'SI' ('EGA' WTAB.'PROCESSEURS' 'COMPORTEMENT') ;
  321. PARALLEL = VRAI ;
  322. PARTLOCA = VRAI ;
  323. MODRELOC = 'PART' 'ARLE' ZMODLI NBPART ;
  324. 'FINS' ;
  325. 'SI' ('EGA' WTAB.'PROCESSEURS' 'AUTOMATIQUE') ;
  326. PARALLEL = VRAI ;
  327. PARTLOCA = FAUX ;
  328. ZMODL = 'PART' 'ARLE' ZMODLI NBPART ;
  329. 'OPTI' 'PARA' VRAI ;
  330. 'FINS' ;
  331.  
  332. * Test sur un MODELE qui aurait changé
  333. 'SI' ('NEG' ZMODLI ZMODLP) ;
  334. CONTI.'CONTRAINTES'=
  335. ('REDU' CONTI. 'CONTRAINTES' ZMODLI) '+'
  336. ('ZERO' ZMODLI 'CONTRAINTES' ) ;
  337. ESTIM.'CONTRAINTES'=
  338. ('REDU' ESTIM. 'CONTRAINTES' ZMODLI) '+'
  339. ('ZERO' ZMODLI 'CONTRAINTES' ) ;
  340. CONTI.'DEFORMATIONS'=
  341. ('REDU' CONTI.'DEFORMATIONS' ZMODLI) '+'
  342. ('ZERO' ZMODLI 'DEFORMAT' ) ;
  343. ESTIM.'DEFORMATIONS'=
  344. ('REDU' ESTIM.'DEFORMATIONS' ZMODLI) '+'
  345. ('ZERO' ZMODLI 'DEFORMAT' ) ;
  346. 'SI' IPLAVI ;
  347. CONTI.'DEFORMATIONS_INELASTIQUES'=
  348. ('REDU' CONTI.'DEFORMATIONS_INELASTIQUES' ZMODLI) '+'
  349. ('ZERO' ZMODLI 'DEFINELA' ) ;
  350. ESTIM.'DEFORMATIONS_INELASTIQUES'=
  351. ('REDU' ESTIM.'DEFORMATIONS_INELASTIQUES' ZMODLI) '+'
  352. ('ZERO' ZMODLI 'DEFINELA' ) ;
  353. CONTI.'VARIABLES_INTERNES'=
  354. ('REDU' CONTI.'VARIABLES_INTERNES' ZMODLI) '+'
  355. ('ZERO' ZMODLI 'VARINTER' ) ;
  356. ESTIM.'VARIABLES_INTERNES'=
  357. ('REDU' ESTIM.'VARIABLES_INTERNES' ZMODLI) '+'
  358. ('ZERO' ZMODLI 'VARINTER' ) ;
  359. 'FINS';
  360.  
  361. WTAB.'ETAT1' = PAS_ETAT PRECED TEMPS0 ;
  362. WTAB.'MAT1' = PAS_MATE PRECED WTAB.'ETAT1';
  363.  
  364. 'SI' ITHER ;
  365. WTAB.'ETHER1'=('REDU' WTAB.'ETHER1' ZMODLI) '+'
  366. ('ZERO' ZMODLI 'DEFORMAT' ) ;
  367. 'SI' POR1;
  368. WTAB.'MSRTHER1'=
  369. ('REDU' WTAB.'MSRTHER1' ZMODLI) '+'
  370. ('ZERO' ZMODLI 'DEFORMAT' ) ;
  371. 'FINS';
  372. 'FINS';
  373.  
  374. 'SI' ('EXIS' WTAB 'DEFOR1' );
  375. WTAB.'DEFOR1' =
  376. ('REDU' WTAB.'DEFOR1' ZMODLI) '+'
  377. ('ZERO' ZMODLI 'DEFORMAT' ) ;
  378. 'FINS';
  379.  
  380. 'SI' ('EXIS' WTAB 'DEFOR2' );
  381. WTAB.'DEFOR2' =
  382. ('REDU' WTAB.'DEFOR2' ZMODLI) '+'
  383. ('ZERO' ZMODLI 'DEFORMAT' ) ;
  384. 'FINS';
  385. 'FINS';
  386.  
  387. 'SI' ITHER ;
  388. WTAB.'ETHER2' WTAB.'MSRTHER2' =PAS_EPTH PRECED WTAB.'MO_TOT'
  389. WTAB.'MAT1' WTAB.'TET2' ;
  390. 'FINS';
  391.  
  392.  
  393. *-----------materiau au debut du pas ------------------------------
  394. ZMAT11 = WTAB.'MAT1' ;
  395. 'SI' ('EGA' WTAB.'PROCESSEURS' 'AUTOMATIQUE') ;
  396. ZMAT1 ='REDU' ZMAT11 ZMODL ; COMM 'Version PARALLEL';
  397. 'SINO';
  398. ZMAT1 = ZMAT11;
  399. 'FINS';
  400. ZMAT1I = ZMAT1; COMM 'Version Initiale de ZMAT1 pour BONOCONV';
  401.  
  402.  
  403. ************************************************************************
  404. * Quelques initialisations *
  405. ************************************************************************
  406. STAB12 = 'TABL' ;
  407. STAB12.'ZU1' = CONTI.'DEPLACEMENTS' ;
  408. STAB12.'SIGF' = CONTI.'CONTRAINTES' ;
  409. STAB12.'DEFF' = CONTI.'DEFORMATIONS' ;
  410. STAB12.'FNONL' = WTAB.'FNONL' ;
  411. STAB12.'RESIDU' = WTAB.'RESIDU' ;
  412. STAB12.'XDENO' = WTAB.'XDENO' ;
  413. STAB12.'XDENOM' = WTAB.'XDENOM' ;
  414. 'SI' ( 'EXIS' WTAB 'ETAT1' );
  415. STAB12.'ETAT1' = WTAB.'ETAT1';
  416. 'FINS';
  417. STAB12.'ETAT2'= ZETAT2;
  418. 'SI' ('EXIS' WTAB 'DEFOR1' );
  419. STAB12.'DEFOR1' = WTAB.'DEFOR1' ;
  420. STAB12.'DEFOR2' = WTAB.'DEFOR2' ;
  421. 'FINS';
  422. 'SI' ( 'EXIS' WTAB 'FNONL');
  423. STAB12.'FNONL' = WTAB.'FNONL' ;
  424. 'FINS';
  425. 'SI' ('EXIS' WTAB 'TET1') ;
  426. STAB12.'TET1' = WTAB.'TET1' ;
  427. STAB12.'TET2' = WTAB.'TET2' ;
  428. 'FINS';
  429. * STAB12.'SUCCES' = VRAI ;
  430. 'SI' ('NEG' WTAB.'AUTOCOEF' 'INCONNU') ;
  431. STAB12.'AUTOCOEF' = WTAB.'AUTOCOEF' ;
  432. 'FINS' ;
  433. 'SI' ('NEG' WTAB.'AUTOREDU' 'INCONNU') ;
  434. STAB12.'AUTOREDU' = WTAB.'AUTOREDU' ;
  435. 'FINS' ;
  436. 'SI' ('NEG' WTAB.'SECOND_MEMBRE' 'INCONNU') ;
  437. STAB12.'SECOND_MEMBRE' = WTAB.'SECOND_MEMBRE' ;
  438. 'FINS' ;
  439. 'SI' ('NEG' WTAB.'LASTKTAN' 'INCONNU') ;
  440. STAB12.'LASTKTAN' = WTAB.'LASTKTAN' ;
  441. 'FINS' ;
  442. 'SI' ('NEG' WTAB.'AUTORED1' 'INCONNU') ;
  443. STAB12.'AUTORED1' = WTAB.'AUTORED1' ;
  444. 'FINS' ;
  445. 'SI' ('NEG' WTAB.'LISEA_M' 'INCONNU') ;
  446. STAB12.'LISEA_M' = WTAB.'LISEA_M' ;
  447. STAB12.'RIBLO_M' = WTAB.'RIBLO_M' ;
  448. 'FINS' ;
  449. 'SI' ('NEG' WTAB.'INCREMENT' 'INCONNU' );
  450. STAB12.'INCREMENT' = WTAB.'INCREMENT';
  451. INCRPREC = STAB12.'INCREMENT' ;
  452. 'FINS';
  453. STAB12.'FFROT' = WTAB.'FFROT' ;
  454. STAB12.'INITEMPS' = WTAB.'INITEMPS' ;
  455. STAB12.'DT' = WTAB.'DT' ;
  456. STAB12.'DTPREC' = WTAB.'DTPREC' ;
  457. 'SI' ITHER ;
  458. ETT0 ='REDU' WTAB.'ETHER1' ZMODL;
  459. STAB12.'TETA1' = WTAB.'TET1';
  460. STAB12.'TETA2' = WTAB.'TET2';
  461. 'SI' POR1;
  462. MSRTT0 = 'REDU' WTAB.'MSRTHER1' ZMODL;
  463. 'FINS';
  464. 'SINON' ;
  465. che3 = 'MANU' 'CHML' ZMODL 'T' 20. 'STRESSES' ;
  466. che4 = 'MANU' 'CHML' ZMODL 'T' 20. 'STRESSES' ;
  467. 'FINS' ;
  468.  
  469. ZMATI = 'REDU' ZMAT22 ZMODLI ;
  470. ZMAT = 'REDU' ZMATI ZMODL ;
  471.  
  472. * SP : initialisation DFGRAD en presence d'un modele MECANIQUE
  473. * (DFGRAD mis a INCONNU par defaut dans PAS_INIT)
  474. 'SI' ('NEG' WTAB.'MOD_MEC' 'INCONNU') ;
  475. 'SI' ('EGA' WTAB.'DFGRAD' 'INCONNU') ;
  476. STAB12.'DFGRAD' = 'GRAD' ZMODL ZMAT STAB12.'ZU1' ;
  477. 'SINO';
  478. STAB12.'DFGRAD' = 'REDU' WTAB.'DFGRAD' ZMODL;
  479. 'FINS' ;
  480. 'SINON' ;
  481. STAB12.'DFGRAD' = WTAB.'DFGRAD' ;
  482. 'FINS' ;
  483.  
  484. com_sig = 'EXTR' ZMODLI 'CONT';
  485. HPP_EPS = FAUX ;
  486. EPS_EPS = 'TEXT' ' ' ;
  487. EPS_NLIN = VRAI ;
  488.  
  489. 'SI' ('EGA' ('VALE' 'EPSI') 'LINEAIRE'); EPS_NLIN = FAUX; 'FINS';
  490. * Option a n'utiliser que par les utilisateurs avertis
  491. 'SI' ('EXIS' PRECED 'ACCELERATION') ;
  492. III = PRECED.'ACCELERATION' ;
  493. 'SI' ('EGA' ('TYPE' III) 'ENTIER') ; ZNACCE = III ; 'FINS' ;
  494. 'FINS' ;
  495. * Matrice tangente : non utilisee si IPLAVI a FAUX
  496. IKTAN = IKTAN 'ET' IPLAVI ;
  497. 'SI' (WTAB.'K_TANGENT' 'ET' ('NON' IPLAVI)) ;
  498. 'MESS' 'IPLAVI faux : pas de matrice tangente ->'
  499. ' on utilise la rigidite elastique' ;
  500. 'FINS' ;
  501. * Matrice tangente par perturbation :
  502. * Option non disponible si non local ou si IPLAVI a FAUX
  503. IPERT = WTAB.'K_TANGENT_PERT' 'ET' ('NON' LNLOC) 'ET' IPLAVI ;
  504. ZPERC1 = WTAB.'K_TANG_PERT_C1' ; ZPERC2 = WTAB.'K_TANG_PERT_C2' ;
  505. * Matrice tangente : partie symetrique utilisee
  506. 'SI' WTAB.'K_TANGENT_SYME' ;
  507. ZKTASYM = 'MOT' 'SYME' ;
  508. 'SINON' ;
  509. ZKTASYM = 'TEXTE' ' ' ;
  510. 'FINS';
  511. * Matrice tangente : pas d'acceleration en cas de modele FEFP ou SSTE
  512. 'SI' IKTAN ;
  513. 'SI' (IFEFP 'OU' ISSTE) ; ZNACCE = 999 ; 'FINS';
  514. 'FINS' ;
  515. *
  516. 'SI' IFTOL ;
  517. ZFTOL = 'ABS' WTAB.'FTOL' ;
  518. 'FINS';
  519. 'SI' IMTOL ;
  520. ZMTOL = 'ABS' WTAB.'MTOL' ;
  521. 'FINS';
  522. ITCAR = ( 'EXIS' ZMATI 'EPAI') 'OU' ('EXIS' ZMATI 'INRY') 'OU'
  523. ( 'EXIS' ZMATI 'MODS') 'OU'
  524. (( 'EXIS' ZMODLI 'ELEM' 'JOI1') 'ET' ('EXIS' ZMATI 'V1X ')
  525. 'ET' ('EXIS' ZMATI 'V1Y ')) ;
  526. *
  527. * CB215821 : Devrait t-on mutualiser ITCAR et WTAB.ITCAR? ==>
  528. * Ce ne sont pas tout a fait les memes tests aujourd'hui
  529. 'SI' (ITCAR 'EGA' FAUX);
  530. ZMAT2 = ZMATI ;
  531. ZMAT2R= ZMAT ;
  532. 'FINS';
  533.  
  534. 'SI' WTAB.'ITCAR';
  535. 'SI' ('EGA' WTAB.'CARA' ('MOT' 'INCONNU'));
  536. WTAB.'CARA' = ZMAT11 ;
  537. 'FINS';
  538. CARA1 = ZMAT1;
  539. 'FINS';
  540.  
  541. 'SI' ('OU' ('OU' IVISCO IVIDOM) IVIEXT); ZPREK = 5.E-7 ; 'FINS';
  542. 'SI' IENDOM; ZPREK = ZPREC ; 'FINS';
  543.  
  544. * on fait ici la séparation poreux .. pour l'avoir sur
  545. *les modèles partitionnes
  546. 'SI' POR1;
  547. MO_PORI = 'EXTR' ZMODLI 'FORM' 'POREUX';
  548. MO_POR = 'EXTR' ZMODL 'FORM' 'POREUX';
  549.  
  550. MA_POR = 'REDU' ZMAT22 MO_POR ;
  551. ** kich ma_por0 intervient si ISOL
  552. ** initialement MA_POR0= 'REDU' MO_POR (STAB12.'MAT1');
  553. MA_POR0 = 'REDU' ZMAT11 MO_POR ;
  554.  
  555. MAI_POR = 'EXTR' MO_POR 'MAILLAGE' ;
  556. MAI_PORI= 'EXTR' MO_PORI 'MAILLAGE' ;
  557. 'FINS';
  558.  
  559. * recuperation de certains champs, si nbpart>1 zmodl est partitionné
  560. * sinon c'est le modele initial
  561. * *
  562. DEFT0 = 'REDU' conti.'DEFORMATIONS' ZMODL ;
  563. ZSIGF = 'REDU' STAB12.'SIGF' ZMODL ;
  564. 'SI' IPLAVI ;
  565. ZDEIF = 'REDU' CONTI.'DEFORMATIONS_INELASTIQUES' ZMODL ;
  566. ZVARF = 'REDU' CONTI.'VARIABLES_INTERNES' ZMODL ;
  567. com_var = 'EXTR' ZMODLI 'VARI' ;
  568. com_dei = 'EXTR' ZMODLI 'DEIN' ;
  569. lnom = com_var ;
  570. 'SI' ISOL ;
  571. com_maa = 'EXTR' ZMODLI 'MATE' ;
  572. 'FINS' ;
  573. 'FINS' ;
  574. *
  575. * teste t'on les moments ?
  576. TSTMOM = ITCAR 'OU' ('EGA' ('VALE' 'MODE') 'PLANGENE') ;
  577. * teste t'on les POREUX ?
  578. 'SI' POR1 ; TSTMOM=VRAI ; 'FINS';
  579. *
  580. IKLFFF=VRAI;
  581. 'SI'TSTMOM; 'SI' IFTOL; 'SI' IMTOL;
  582. IKLFFF=FAUX;
  583. 'FINS'; 'FINS'; 'FINS';
  584. *
  585. 'SI'('NON' TSTMOM); 'SI' IFTOL;
  586. IKLFFF=FAUX;
  587. 'FINS';'FINS';
  588. *
  589. GEOREF0 = WTAB.'FOR0' ;
  590. WTAB.'CONV'=VRAI;
  591. WTAB . 'ISOUSPAS' = 0;
  592. NSOUSPAS = WTAB . 'MAXSOUSPAS';
  593. ZCCONV = VRAI ;
  594. KNOCONV = 0 ;
  595.  
  596.  
  597. ************************************************************************
  598. ****** boucle de non convergence
  599. ************************************************************************
  600. 'REPETER' BONOCONV NSOUSPAS ;
  601.  
  602. KNOCONV = KNOCONV+1 ;
  603. STAB12.'CONV' = FAUX ;
  604. DT_INIT = STAB12.'DT' ;
  605. DTINI = STAB12.'DT' ;
  606.  
  607. ZSIG0 = ZSIGF ;
  608. 'SI' IPLAVI ;
  609. ZEPS0 = ZDEIF ;
  610. ZVAR0 = ZVARF ;
  611. 'FINS';
  612.  
  613. ZU1 = STAB12.'ZU1' ;
  614. 'SI' LAG_TOT;
  615. mzu1 = zu1 '*' -1;
  616. 'FINS';
  617. GR_U_DEB = STAB12.'DFGRAD';
  618. 'SI' ITHER ;
  619. TETA1 = STAB12.'TETA1' ;
  620. TETA2 = STAB12.'TETA2' ;
  621. DTETD = TETA2 '-' TETA1;
  622. 'FINS' ;
  623.  
  624. * materiau au debut du pas en cas de non convergence etat1=etat2
  625. 'SI' (knoconv > 1) ;
  626. ZMAT1 = 'REDU' (PAS_MATE PRECED STAB12.'ETAT1') ZMODL ;
  627. 'SINON';
  628. ZMAT1 = ZMAT1I ;
  629. 'FINS' ;
  630. *------ caracteristiques initiales en cas de grands deplacements ------
  631. 'SI' WTAB.'ITCAR';
  632. MECAR1 = 'EXCO' MOCA ZMAT1 'NOID' ;
  633. MECAR2 = 'EXCO' MOCA CARA1 'NOID' ;
  634. MECAR1 = 'CHAN' 'TYPE' MECAR1 'CARACTERISTIQUES';
  635. MECAR2 = 'CHAN' 'TYPE' MECAR2 'CARACTERISTIQUES';
  636. ZMAT1 = ZMAT1 '-' MECAR1 '+' MECAR2 ;
  637. 'FINS';
  638.  
  639. *----------- Calcul du champ de materiau a la fin du pas ------------
  640. 'SI' (knoconv '>' 1) ;
  641. MCHC = PAS_ETAT PRECED TI ;
  642. MMMM = PAS_MATE PRECED MCHC;
  643. 'SINON' ;
  644. MMMM = ZMAT22 ;
  645. 'FINS' ;
  646. MMMM = 'REDU' MMMM ZMODL;
  647.  
  648. *----- Caracteristiques initiales en cas de grands deplacements -----
  649. 'SI' WTAB.'ITCAR';
  650. MECAR1 = 'EXCO' MOCA MMMM 'NOID' ;
  651. MECAR2 = 'EXCO' MOCA CARA1 'NOID' ;
  652. MECAR1 = 'CHAN' 'TYPE' MECAR1 'CARACTERISTIQUES';
  653. MECAR2 = 'CHAN' 'TYPE' MECAR2 'CARACTERISTIQUES';
  654. MMMM = MMMM '-' MECAR1 '+' MECAR2 ;
  655. * (fdp) on reporte les variations du materiau sur le pas
  656. * dans les autres instances du champ materiau
  657. * (ZMAT et ZMATI semblent suffirent)
  658. MECAR1 = 'EXCO' MOCA ZMAT 'NOID' ;
  659. MECAR2 = 'EXCO' MOCA MMMM 'NOID' ;
  660. MECAR1 = 'CHAN' 'TYPE' MECAR1 'CARACTERISTIQUES';
  661. MECAR2 = 'CHAN' 'TYPE' MECAR2 'CARACTERISTIQUES';
  662. ZMAT = ZMAT '-' MECAR1 '+' MECAR2 ;
  663. ZMATI = 'REDU' ZMAT ZMODLI ;
  664. 'FINS';
  665.  
  666. *------------ Calcul de la rigidite a la fin du pas ----------------
  667. * Test ci-dessous est inutile, WTAB.'BLOCAGES_MECANIQUES' et ZCLIM0
  668. * correspondent aux memes rigidites ... (cf. ligne 117)
  669. AA1 ='EXTR' WTAB.'BLOCAGES_MECANIQUES' 'MAIL' ;
  670. AA2 ='EXTR' ZCLIM0 'MAIL' ;
  671. AA3 ='DIFF' AA1 AA2 ;
  672. AA4 = NBNO AA3 ;
  673. 'SI' ( (WTAB.'RECARI' 'ET' WTAB.'CONV') 'OU'
  674. ('NON' ('EXIS' WTAB 'RRRR') ) 'OU' ('NEG' AA4 0));
  675. 'SI' ('OU' ('OU' IENDOM IVIDOM) ICERAM);
  676. HOOKENDO ='HOOK' ZMODL MMMM ZVAR0;
  677. RH ='RIGI' ZMODL HOOKENDO MMMM ;
  678. 'DETR' HOOKENDO;
  679. 'SINON';
  680. RH = 'RIGI' ZMODL MMMM ;
  681. 'FINS';
  682. * RH peut contenir des CL
  683. ZCL = 'EXTR' RH 'RIGI' 'MULT' ;
  684. 'SI' ('NEG' ('DIME' ZCL) 0) ;
  685. ZCLIM = ZCLIM0 'ET' ZCL ;
  686. 'FINSI' ;
  687. RRRR = RH 'ET' ZCLIM0 ;
  688. * Prise en compte d'eventuelles RIGIDITE_CONSTANTE
  689. 'SI' IRCON;
  690. RRRR = RRRR 'ET' RIG_CONS;
  691. 'FINS';
  692. *
  693. 'SI' IRAUG;
  694. RRRR = RRRR 'ET' RIG_AUG ;
  695. 'FINS';
  696.  
  697. 'SI' ('EGA' ('DIME' ZCL) 0) ;
  698. * Stockage de la rigidite pour eviter de la recalculer
  699. WTAB.'RRRR'=RRRR;
  700. 'FINSI';
  701. *
  702. 'SINON';
  703. RRRR=WTAB.'RRRR';
  704. 'FINS';
  705. ZRAID=RRRR;
  706.  
  707. *------------ consolidation ou dynamique faut-il recalculer l'operateur?
  708. * -----: preparation du pas de temps ------------
  709. 'SI' ( WTAB.'CONSOLIDATION' 'OU' WTAB.'DYNAMIQUE');
  710. DT = TI '-' TEMP0;
  711. 'SI' ( '>' (DELTAN '*' 0.9999) DT) ;
  712. WTAB . 'RECAOP' = VRAI;
  713. 'FINS';
  714. 'SI' ( '<' (DELTAN '*' 1.0001) DT) ;
  715. WTAB . 'RECAOP' = VRAI;
  716. 'FINS';
  717. 'SI' WTAB.'MATVAR';
  718. WTAB . 'RECAOP' = VRAI;
  719. 'FINS';
  720. DELTAN=WTAB.'DT';
  721. 'FINS';
  722. *---------------------- Formation de l operateur -----------------------
  723. 'SI' ( WTAB.'DYNAMIQUE' 'OU' WTAB.'CONSOLIDATION');
  724. 'SI' ('NEG' WTAB.'OPERATEUR' 'INCONNU');
  725. ZRAID = WTAB.'OPERATEUR';
  726. 'SI' (WTAB . 'RECAOP') ;
  727. ZRAID = RRRR ;
  728. 'FINS';
  729. 'FINS';
  730. 'FINS';
  731.  
  732. *------------ operateur frequentiel -----------------------
  733. *------------ operateur amortissement en frequentiel
  734. 'SI' WTAB.'FREQUENTIEL' ;
  735. RRR2 = 'AMOR' ZMODL ZMAT ;
  736. RR2 = 'CHAN' 'INCO' RRR2
  737. ('MOTS' 'ALFA' 'BETA') ('MOTS' 'IALF' 'IBET')
  738. ('MOTS' 'FALF' 'FBET') ('MOTS' 'FALF' 'FBET') 'QUEL' ;
  739. RR3 = 'CHAN' 'INCO' RRR2
  740. ('MOTS' 'ALFA' 'BETA') ('MOTS' 'ALFA' 'BETA')
  741. ('MOTS' 'FALF' 'FBET') ('MOTS' 'IFAL' 'IFBE') 'QUEL' ;
  742. RRR2 = RR2 'ET' RR3 ;
  743.  
  744. RR1 = ZRAID ;
  745. OMEGI= 2.* PI * TI ;
  746. RRR1 = OMEGI * OMEGI * (-1.) * WTAB.'MASSE' ;
  747. RR1 = ZRAID 'ET' RRR1 ;
  748. RR4 = 'CHAN' 'INCO' (RR1 '*' (-1.D0))
  749. ('MOTS' 'ALFA' 'BETA') ('MOTS' 'IALF' 'IBET')
  750. ('MOTS' 'FALF' 'FBET') ('MOTS' 'IFAL' 'IFBE') 'QUEL';
  751. ZRAID= RR1 'ET' RR4 ;
  752. RR5 = OMEGI '*' RRR2 ;
  753. ZRAID= ZRAID 'ET' RR5 ;
  754. 'FINS' ;
  755.  
  756. *--------------- et la perméabilité ----------------------------------
  757. 'SI' (WTAB.'CONSOLIDATION') ;
  758. 'SI' (WTAB.'GRANDS_DEPLACEMENTS' 'OU' WTAB.'MATVAR') ;
  759. WTAB .'PERMEABILITE'= 'PERM' WTAB.'MOD_POR' MMMM ;
  760. WTAB .'RECAOP' = VRAI ;
  761. 'FINS';
  762. 'FINS';
  763.  
  764. *------------- Cas de la consolidation ou de la dynamique -------------
  765. *------------- il faut recalculer l'operateur d'iteration -------------
  766. 'SI' (WTAB . 'RECAOP') ;
  767. 'SI' ( WTAB.'DYNAMIQUE');
  768. ZRAID = 4.D0 '/'( DT ** 2) '*' WTAB.'MASSE' 'ET' ZRAID;
  769. 'SI' ('NEG' WTAB.'AMORTISSEMENT' 'INCONNU');
  770. ZRAID = WTAB.'AMORTISSEMENT' '*' (2.D0 '/' DT) 'ET' ZRAID;
  771. 'FINS';
  772. 'FINS' ;
  773.  
  774. 'SI' ( WTAB.'CONSOLIDATION');
  775. ZRAID =-1.* DT* WTAB.'TETA'* WTAB.'PERMEABILITE' 'ET' ZRAID ;
  776. 'FINS' ;
  777. WTAB . 'OPERATEUR'= ZRAID ;
  778. WTAB . 'RECAOP' = FAUX ;
  779. 'FINS';
  780.  
  781. *-------------- traitement des contacts frottements automatiques -------
  782. CDEP = STAB12.'ZU1' ;
  783. CDEPSLX ='ENLE' CDEP 'LX' ;
  784.  
  785. BFCONT = FAUX ;
  786. BCLIM2 = FAUX ;
  787. 'SI' WTAB.'CONTACT';
  788. MODCON= WTAB.'MODCONTA';
  789. WTAB.'CONT_REEL'=MODCON;
  790. 'SI' ('EXIS' WTAB 'MAT_JEU') ;
  791. CDAP CRR RFROT='RFCO' MODCON WTAB.'CONV' WTAB.'MAT_JEU' ;
  792. 'SINON' ;
  793. CDAP CRR RFROT='RFCO' MODCON WTAB.'CONV' ;
  794. 'FINS';
  795.  
  796. 'SI' (WTAB . 'MODAL') ;
  797. CRR = 'PJBA' CRR ZMODL MMMM ;
  798. MCRR = 'EXTR' CRR 'MAIL' 'MULT' ;
  799. MCDAP= 'EXTR' CDAP 'MAIL' ;
  800. 'REPETER' BCDA ('NBNO' MCDAP) ;
  801. PBCDA ='POIN' MCDAP &BCDA ;
  802. PCRR ='POIN' MCRR &BCDA ;
  803. CHCR ='MANU' 'CHPO' PCRR 1 'FLX' ('EXTR' CDAP 'FLX' PBCDA)
  804. 'NATURE' 'DISCRETE' ;
  805. 'SI' ('EGA' 1 &BCDA) ;
  806. CCDA = CHCR ;
  807. 'SINON' ;
  808. CCDA = CHCR 'ET' CCDA ;
  809. 'FINS' ;
  810. 'FIN' BCDA ;
  811. CDAP = CCDA ;
  812. 'FINS' ;
  813.  
  814. * ATTENTION : mettre les conditions de frottement en premier pour
  815. * les numeroter en dernier
  816. ZCLIM2 = 'VIDE' 'RIGIDITE';
  817. ZFCONT ='VIDE' 'CHPOINT' / 'DISCRET';
  818. 'SI' ('NEG' CRR 0);
  819. BCLIM2 = VRAI ;
  820. BFCONT = VRAI ;
  821. ZCLIM2 = CRR 'ET' ZCLIM2 ;
  822. CCOR = CRR '*' CDEPSLX ;
  823. ZFCONT = ZFCONT '+' CCOR ;
  824. 'FINS';
  825.  
  826. 'SI' ( 'NEG' 0 CDAP);
  827. BFCONT = VRAI ;
  828. ZFCONT = ZFCONT '+' CDAP ;
  829. 'FINS';
  830.  
  831. 'SI' ('NEG' 0 RFROT);
  832. BCLIM2 = VRAI ;
  833. BFCONT = VRAI ;
  834. ZCLIM2 = RFROT 'ET' ZCLIM2 ;
  835. CCOR = RFROT '*' CDEPSLX ;
  836. ZFCONT = ZFCONT '+' CCOR ;
  837. 'FINS';
  838. *
  839. * Mise a jour de ZCLIM et ZFCONSTA si necessaire
  840. 'SI' BCLIM2;
  841. ZCLIM = ZCLIM2 'ET' ZCLIM ;
  842. ZRAID = ZCLIM2 'ET' ZRAID;
  843. 'FINS';
  844. *
  845. 'SI' BFCONT;
  846. ZFCONSTA = ZFEXT2 '+' ZFCONT ;
  847. 'FINS';
  848. 'FINS';
  849. *
  850. ZFCONST1 = ZFCONSTA;
  851.  
  852. * --------------- pilotage automatique ********************************
  853. ISNPB = FAUX ;
  854. AL1 = 1. ;COEPI = 1.d0; COEINC=0.d0;COEPI0=1.d0;DAL1=100.D0;
  855. CORPREC = 1. ;
  856. * CORPREC = 10.;
  857. 'SI' (('EGA' ipredic 'HPP') 'ET' WTAB.'CONV');
  858. EPS_EPS = 'TEXTE' 'LINEAIRE';
  859. HPP_EPS = VRAI;
  860. 'FINS';
  861. *
  862. 'SI' IPILOT ;
  863. 'SI' ( WTAB.'AUTODEUX' ) ;
  864. COEPI = 'ABS' ( STAB12.'AUTOCOEF');
  865. COEPI = COEPI / (1.-COEPI);
  866. 'SI' (COEPI > 1.D0) ; COEPI=1.D0;'FINS';
  867. COEPI0=COEPI;STAB12.'AUTOCOEF'=COEPI;
  868. 'SINON';
  869. STAB12.'AUTOCOEF' = 1.D0;
  870. 'FINS' ;
  871. RED1 = 1. ;
  872. RED2 = 0 ;
  873. 'SI' ('NEG' WTAB.'AUTORED1' 'INCONNU') ;
  874. 'SI' (STAB12.'AUTORED1' > 0);
  875. STAB12.'AUTORED1' = STAB12.'AUTORED1' - 1;
  876. 'SI' (STAB12.'AUTORED1' 'EGA' 0) ;
  877. COEPI = 3 * COEPI ;
  878. STAB12.'AUTOREDU' = STAB12.'AUTOREDU' / 3.;
  879. 'SI' (STAB12.'AUTOREDU' > 1.1 );
  880. * on travaille encore avec un critere reduit
  881. STAB12.'AUTORED1' = 4 ;
  882. RED1 = 3. ;
  883. 'FINS';
  884. 'MESS' 'On multiplie le critere de pilotage par 3';
  885. 'FINS';
  886. 'FINS';
  887. 'FINS';
  888. 'SI' (COEPI > 1d0);
  889. RED1 = RED1 / COEPI ;COEPI =1d0;
  890. 'FINS';
  891. 'SI' ( 'NEG' WTAB.'NBPLAS' 'INCONNU') ;
  892. 'SI' ( WTAB.'NBPLAS' 'EGA' 0) ;
  893. 'SI' (COEPI < 0.) ;COEPI = COEPI * -2.;'FINS';
  894. 'FINS';
  895. 'FINS';
  896. COEPI = 'ABS' COEPI ;
  897. COEPI0 = COEPI;
  898. * sans pilotage
  899. 'SINON';
  900. STAB12.'AUTOCOEF'= 1.D0;
  901. 'FINS';
  902.  
  903. ************************************************************************
  904. *------------- quelques initialisations pour la boucle ETIQ ********
  905. ************************************************************************
  906. URG = FAUX; RED_URG = 0 ; IT = 0 ; c_zdepr = faux;ITACC = 0;
  907. ZICONV = VRAI; MMC = 0 ; MMCMAX = 0 ; EPSM = 0.; DPSMAX = 0. ;
  908. DEPSTDM = 0. ; DEKREAC1 = 0. ;XCONVNOR = 0. ; ITNORM1 = 0 ;
  909. DEPSTREF = 100. * (WTAB . 'MAXDEFOR') ;ZDEPL=0;
  910. GR_U_K = GR_U_DEB ;
  911. DITNORM1 = 0 ;NBCYCLE1 = 0 ;zdept = 0; zdeptq = 0; zdeptp = 0;
  912. zprecnc=1e-5;FTHE = 0. ; FDEF = 0. ; ITNV = -5 ;
  913. TABCONV = 'TABL';
  914. 'SI' IFEFPUL;
  915. XUPDA = 1;
  916. 'FINS';
  917.  
  918. DEPST0=0;
  919. *********** en cas de materiaux variables ****************************
  920. 'SI' (WTAB.'MATVAR');
  921. 'SI' ('OU' ('OU' IENDOM IVIDOM) ICERAM);
  922. XXX3 = 'ELAS' ZMODL ZSIG0 ZMAT ZVAR0;
  923. XXX4 = 'ELAS' ZMODL ZSIG0 ZMAT1 ZVAR0;
  924. 'SINON';
  925. XXX3 = 'ELAS' ZMODL ZSIG0 ZMAT ;
  926. XXX4 = 'ELAS' ZMODL ZSIG0 ZMAT1 ;
  927. 'FINS';
  928. DDEF0 = XXX4 - XXX3;
  929. 'DETR' XXX3;
  930. 'DETR' XXX4;
  931. DEPST0=-1.* DDEF0;
  932. 'FINS';
  933.  
  934. ***** En cas de chargement thermiques *********************************
  935. 'SI' ITHER ;
  936. 'SI' (WTAB.'MATVAR' 'ET' IPILOT);
  937. 'MESS' 'Le pilotage n est pas possible avec un'
  938. ' materiau qui depend de la temperature' ;
  939. 'ERREUR' 19 ;
  940. 'FINS' ;
  941.  
  942. ETT MSRTT = PAS_EPTH PRECED ZMODL ZMAT TETA2 ;
  943. DTT = ETT '-' ETT0 ;
  944. *--------- Cas du milieu poreux avec chargement thermique ----------
  945. * cas isotrope seulement pour le moment
  946. * et on ne s'occupe pas du alpha-reference !!
  947. 'SI' POR1 ;
  948. DMSRT0 = MSRTT '-' MSRTT0 ;
  949. 'FINS' ;
  950.  
  951. 'SI' ('EGA' DEPST0 0);
  952. DEPST0 = DTT;
  953. 'SINON';
  954. DEPST0 = DTT '+' DEPST0 ;
  955. 'FINS' ;
  956. 'FINS';
  957.  
  958. * calcul de dsig0 et fthe en cas de chargement si necessaire ***********
  959. 'SI' ( ITHER 'OU' WTAB.'MATVAR');
  960. 'SI' ('OU' ('OU' IENDOM IVIDOM) ICERAM);
  961. DSIGT0 = 'ELAS' ZMODL DEPST0 ZMAT ZVAR0;
  962. 'SINON';
  963. DSIGT0 = 'ELAS' ZMODL DEPST0 ZMAT ;
  964. 'FINS';
  965.  
  966. 'SI' (POR1 'ET' ITHER);
  967. DSIGT0 = DSIGT0 '+' DMSRT0;
  968. 'FINS';
  969. FTHE = 'BSIGMA' ZMODL DSIGT0 ZMAT ;
  970. 'FINS';
  971.  
  972. *-------------- deformations imposes **********************************
  973. 'SI' LOGDEF;
  974. DDEFOR0 ='REDU' (STAB12.'DEFOR2' - STAB12.'DEFOR1') ZMODL;
  975. 'SI' ('OU' ('OU' IENDOM IVIDOM) ICERAM);
  976. DSI1 = 'ELAS' ZMODL DDEFOR0 ZMAT ZVAR0;
  977. 'SINON';
  978. DSI1 = 'ELAS' ZMODL DDEFOR0 ZMAT ;
  979. 'FINS';
  980. FDEF = 'BSIGMA' ZMODL DSI1 ZMAT ;
  981. 'FINS';
  982.  
  983. *-------------- pression imposee et grands deplacements ?***************
  984. 'SI' (LOGPRE 'ET' IGRD) ;
  985. MOP = WTAB.'MOD_PRE' ;
  986. ZPEXT0 = 'TIRE' WTAB.'CHARGEMENT' 'PRES' TEMPS0 ;
  987. ZPEXTF = 'TIRE' WTAB.'CHARGEMENT' 'PRES' TI ;
  988. *
  989. 'SI' ('EXIS' WTAB 'MAT_PRE') ;
  990. ZFPEXT0 = 'BSIG' MOP ZPEXT0 ('REDU' ZMAT11 MOP) ;
  991. ZFPEXTF = 'BSIG' MOP ZPEXTF ('REDU' ZMAT22 MOP) ;
  992. 'SINON' ;
  993. ZFPEXT0 = 'BSIG' MOP ZPEXT0 ;
  994. ZFPEXTF = 'BSIG' MOP ZPEXTF ;
  995. 'FINS' ;
  996. ZFP0F = 'COPIER' ZFPEXTF ; COEFP=1.D0;
  997. 'FINS';
  998.  
  999. * --------- ktangent et fefp******************************************
  1000. 'SI' IKTAN ;
  1001. 'SI' IFEFP ;
  1002. IKT_SAUV = VRAI ;
  1003. 'SI' ('NEG' WTAB.'LASTKTAN' 'INCONNU') ;
  1004. 'MESS' 'FEFP: Start with LASTKTAN' ;
  1005. ZRIKTA = STAB12.'LASTKTAN' ;
  1006. ZRAID = ZCLIM 'ET' ZRIKTA ;
  1007. 'SINON' ;
  1008. 'MESS' 'FEFP: Previous KTAN not available' ;
  1009. ZRAID = ZRAID 'ET' ('KSIGMA' ZMODL ZSIG0 ZMAT) ;
  1010. 'FINS' ;
  1011. 'SINON' ;
  1012. 'SI' ('NEG' WTAB.'LASTKTAN' 'INCONNU') ;
  1013. IKT_SAUV = VRAI ;
  1014. 'SI' IPERT ;
  1015. 'MESS' 'Matrice tangente par perturbation - '
  1016. 'Demarrage avec KTAN = LASTKTAN' ;
  1017. 'SINON' ;
  1018. 'MESS' 'Matrice tangente "coherente" - '
  1019. 'Demarrage avec KTAN = LASTKTAN' ;
  1020. 'FINS' ;
  1021. ZRIKTA = STAB12.'LASTKTAN' ;
  1022. ZRAID = ZCLIM 'ET' ZRIKTA ;
  1023. 'SINON' ;
  1024. 'SI' IPERT ;
  1025. IKT_SAUV='NEG' WTAB.'K_TANGENT_ITER0'
  1026. 'MAT_ELASTIQUE';
  1027. 'MESS' 'Matrice tangente par perturbation - '
  1028. 'Demarrage avec KTAN = rigidite elastique' ;
  1029. 'SINON' ;
  1030. IKT_SAUV = ('NEG' WTAB.'K_TANGENT_ITER0' 'MAT_ELASTIQUE')
  1031. 'ET' ('NEG' WTAB.'K_TANGENT_ITER0' 'MAT_TANGENTE') ;
  1032. 'SI' ('EGA' WTAB.'K_TANGENT_ITER0'
  1033. 'MAT_ELASTIQUE') ;
  1034. 'MESS' 'Matrice tangente "coherente" - '
  1035. 'Demarrage avec KTAN = rigidite elastique' ;
  1036. 'SINON' ;
  1037. 'MESS' 'Matrice tangente "coherente" - '
  1038. 'Demarrage avec KTAN (DTTAN = 0.)' ;
  1039. DTTAN = 0. ;
  1040. ZRIKTA = 'KTAN' ZMODL ZSIG0 ZVAR0 ZMAT
  1041. 'PREC' ZPREK 'DT ' DTTAN ZKTASYM ;
  1042. ZRAID = ZCLIM 'ET' ZRIKTA ;
  1043. 'FINS' ;
  1044. 'FINS' ;
  1045. 'FINS' ;
  1046. 'FINS' ;
  1047. 'FINS' ;
  1048.  
  1049. *-------- en grands deplacements option K_SIGMA ***********************
  1050. *
  1051. 'SI' (IGRD 'ET' ('NON' HPP_EPS));
  1052. 'SI' (IKSIA 'ET' ('NON' IFEFP)) ;
  1053. KSI1 ='KSIGMA' ZMODL ZSIG0 ZMAT;
  1054. ZRAIDINI= ZRAID ;
  1055. ZRAID = ZRAID 'ET' KSI1 ;
  1056. 'FINS' ;
  1057. 'FINS' ;
  1058.  
  1059. * Y a-t-il des forces non conservatives ( forces suiveuses)? ***********
  1060. ADDISEC0 = FAUX;
  1061. ADDISEC2 = FAUX;
  1062. 'SI' WTAB.'PROCEDURE_CHARMECA';
  1063. * on ajoute l indice ADDI_MATRICE pour signaler a charmeca qu on
  1064. * souhaite aussi l operateur linearisé des Forces NL de charmeca
  1065. PRECED . 'ADDI_MATRICE' = vrai;
  1066. TFP22 = CHARMECA PRECED WTAB.'T_FINAL';
  1067. PRECED . 'ADDI_MATRICE' = faux;
  1068. * FP22 = F^suiv_n+1
  1069. DMZPRES = 0.D0 ;
  1070. 'SI' ('EXIS' TFP22 'ADDI_SECOND') ;
  1071. FP22 = TFP22.'ADDI_SECOND' ;
  1072. FP022 = 'COPIER' FP22 ;
  1073. MZPRES = 'MAXI' 'ABS' FP22 ;
  1074. DMZPRES = MZPRES ;
  1075. ADDISEC2= VRAI ;
  1076. 'FINS';
  1077. 'SI' ('EXIS' TFP22 'ADDI_MATRICE');
  1078. ZRAID = ZRAID 'ET' TFP22.'ADDI_MATRICE';
  1079. 'FINS';
  1080.  
  1081. * FP0 = F^suiv_n
  1082. TFP0 = CHARMECA PRECED TEMPS0 ;
  1083. 'SI'('EXIS' TFP0 'ADDI_SECOND') ;
  1084. FP0 = TFP0.'ADDI_SECOND' ;
  1085. MZPRES0 = 'MAXI' 'ABS' FP0 ;
  1086. DMZPRES = DMZPRES '-' MZPRES0 ;
  1087. ADDISEC0= VRAI ;
  1088. 'FINS';
  1089. * 'SI' ('EXIS' TFP0 'ADDI_MATRICE');
  1090. * ZRAID = ZRAID 'ET' TFP0.'ADDI_MATRICE';
  1091. * 'FINS';
  1092. COEFP = 1.D0 ;
  1093. 'FINS';
  1094.  
  1095.  
  1096. * -----------calcul de la partie constante du second membre **********
  1097. * en consolidation
  1098. * ZFP1 est cense contenir : - B0*SIG0 et
  1099. * DT*(1-TETA)*FI0 + DT*H*P
  1100. *
  1101. * dans ZFCONSTA on met le second membre de u ***********************
  1102.  
  1103. * en dynamique ********************************************************
  1104. * ZFP1 est cense contenir : F0 + 4/DT*M*V0 - B0*SIG0
  1105. 'SI' IDYN ;
  1106. UNSURH = 1.D0 '/' STAB12.'DT' ;
  1107. ZFP1 = WTAB.'FREA1' ;
  1108. ZDYFEXT = ZFCONSTA 'ENLEVER' 'FLX';
  1109. ZFCONSTA= ZFCONSTA '+' ZFP1 ;
  1110. 'FINS';
  1111. ZFEXT = ZFCONSTA 'ENLEVER' 'FLX';
  1112.  
  1113. *---------deplacement (ou jeu) e imposer e la fin du pas *************
  1114. * on separe les efforts ZFEXT (=F^ext_n+1) ***************************
  1115. * et deplacement (ou jeu) ZFLX1 (u^imp_n+1) a imposer a la fin du pas
  1116. ZFLX1 = 'EXCO' ZFCONSTA 'FLX' 'NOID' 'FLX' 'NATURE' 'DISCRET';
  1117. 'SI' ('NEG' STAB12.'FFROT' 'INCONNU'); FFROT = STAB12.'FFROT';
  1118. 'SINON' ; FFROT=ZFEXT * 0; 'FINS';
  1119. FFROTP = FFROT;
  1120.  
  1121. * calcul des forces externes deja equilibrees au debut du pas ********
  1122. * par B*SIGMA : ZF1 = F^int_n = B*sigma_n + K^cst*u_n
  1123. ZF1 = 'BSIGMA' ZMODL ZSIG0 ZMAT ;
  1124. 'SI' IRCON;
  1125. ZF1 = ZF1 '+' ('REDU' MAI_CONS (RIG_CONS '*' ZU1));
  1126. 'FINS';
  1127. 'SI' IDYN ; FFDYN = 'COPIER' ZF1; 'FINS';
  1128. 'SI' ISOL ;
  1129. GRAP0= 'GRAD' MO_POR ZU1 MA_POR0 'CONS' ;
  1130. XXX1 = 'GRAD' MO_POR ZU1 MA_POR 'CONS' ;
  1131. XXXS =((1.- WTAB. 'TETA' )*GRAP0)+ (WTAB. 'TETA' * XXX1);
  1132. XXX2 = STAB12.'DT' * ('GNFL' MO_POR XXXS) ;
  1133. XXX3 = ZF1 ;
  1134. ZF1 = XXX3 - XXX2;'DETR' XXX3;
  1135. 'DETR' XXX2 ;
  1136. 'FINS';
  1137.  
  1138. * initialisation des variables forces et deplacement*******************
  1139. * zzd est le deplacement au pas precedent (=u_n) et ZLX=lambda_n
  1140. ZZD ='ENLE' ZU1 'LX';
  1141. ZLX = ZU1 'EXCO' 'LX' 'NOID' 'LX' 'NATURE' 'DIFFUS';
  1142. * --------- flxini est la partie des FLX deja realisee au debut du pas
  1143. FLXINI= ZZD '*' ZCLIM;
  1144. * FREAP : -1*reactions du pas precedent (=F^reac_n) transportees sur les nou
  1145. FREAP = ZLX '*' ZCLIM;
  1146. FEXT0 = ZF1 '+' FREAP;
  1147. * FEXT0 est le chargement externe (sans reactions vu par la
  1148. * structure le pas d'avant) = F^ext_n = F^int_n - F^reac_n
  1149. *
  1150. * on va calculer le premier residu c'est a dire le desequilibre *******
  1151. * entre les forces externes et le calcul B*SIGMA.
  1152. * le sigma qui sert est celui qui existerait si le champ de
  1153. * deplacement ne changeait pas (ZU1). faire attention aux FLX
  1154. * En pilotage on reprend ce residu que l'on multiplie par COEPI
  1155. * XXX1 = [F^ext_n+1 ; Du^imp]
  1156. * DFEXT0 = increment des forces et des FLX a imposer (le residu
  1157. * du pas precedent) = [DF^ext ; Du^imp]
  1158. XXX1 = ZFCONSTA '-' FLXINI ;
  1159. DFEXT0 = XXX1 '-' FEXT0 ;
  1160.  
  1161. * si pression suiveuse dfext0 contient en plus l'increment des forces
  1162. * de pression du uniquement a la reactualisation de la geometrie (sans
  1163. * augmentation du module)
  1164. * mais comme F^int_n equilibre deja F^ext_n + F^suiv_n +... et qu'on est
  1165. * toujours sur config_n, on doit avoir : DFEXT0= [DF^ext ; Du^imp]
  1166. * avec DF^ext qui ne contient pas de forces suiveuses (...a vérifier)
  1167. 'SI' ADDISEC0;
  1168. DFEXT0 = DFEXT0 '+' FP0 ;
  1169. 'FINS';
  1170. 'SI' (LOGPRE 'ET' IGRD) ;
  1171. DFEXT0 = DFEXT0 '+' ZFPEXT0 ;
  1172. 'FINS' ;
  1173. DFEXT0F = DFEXT0 'ENLEVER' 'FLX';
  1174. DFEXT0L = DFEXT0 'EXCO' 'FLX' 'NOID' 'FLX';
  1175.  
  1176. * DFEXT0 = [DFEXT0F ; DFEXT0L]
  1177. * = [F^ext_n+1 - (F^int_n - F^reac_n - F^suiv_n) ; Du^imp]
  1178. * RESIDU = forces exterieures sans reactions
  1179. * (avec des termes supplementaires le cas echeant p.ex. en
  1180. * dynamique ou en poreux) - forces interieures
  1181. * et increment des relations imposees
  1182. * ---> la resolution fournira dU et dR
  1183.  
  1184. * RESIDU = [F^ext_n+1 + F^ther_n+1 + F^defi_n+1 - F^int_n ; Du^imp]
  1185. RESIDU = XXX1 '+'FTHE '+'FDEF '-'ZF1;
  1186. ZDFINI ='COPIER' DFEXT0 ;
  1187. ZFPLO = ZF1 '-' FTHE '-' FDEF ;
  1188. 'SI' (ITHER 'OU' WTAB.'MATVAR');
  1189. ZDFINI=ZDFINI '+' FTHE;
  1190. 'FINS';
  1191. 'SI' LOGDEF;
  1192. ZDFINI= ZDFINI '+' FDEF;
  1193. 'FINS';
  1194. 'SI' (LOGPRE 'ET' IGRD) ;
  1195. RESIDU = RESIDU '+' ZFPEXTF ;
  1196. ZDFINI = ZDFINI '+' ZFPEXTF '-' ZFPEXT0 ;
  1197. 'FINS';
  1198.  
  1199. 'SI' ADDISEC0 ;
  1200. ZDFINI = ZDFINI '-' FP0 ;
  1201. 'FINS';
  1202.  
  1203. 'SI' ADDISEC2 ;
  1204. * mess ' fp22 ' ; list resu fp22;
  1205. RESIDU = RESIDU '+' FP22 ;
  1206. * on tient compte de l'increment des forces suiveuses en direction
  1207. * et en module
  1208. ZDFINI = ZDFINI '+' FP22 ;
  1209. 'FINS';
  1210.  
  1211. * ici, on a :
  1212. * RESIDU = [F^ext_n+1 + F^ther_n+1 + F^defi_n+1 + F^suiv_n+1
  1213. * - F^int_n ; Du^imp]
  1214. * = [ DF^tot ; Du^imp]
  1215. * ZDFINI = [F^ext_n+1 + F^ther_n+1 + F^defi_n+1 + F^suiv_n+1
  1216. * - (F^int_n - F^reac_n) ; Du^imp]
  1217. *
  1218. IMPO12= FAUX;
  1219. XXX1 ='EXTR' ZCLIM 'MAIL' 'UNIL';
  1220. 'SI' (('NBEL' XXX1) '>' 0);
  1221. IMPO12 = VRAI;
  1222. DIMPO12='REDU' DFEXT0L XXX1 ;
  1223. DIMPOV = DFEXT0L '-' DIMPO12;
  1224. 'FINS';
  1225. stab12.'SECOND_MEMBRE' = RESIDU '*' 1.D0;
  1226.  
  1227.  
  1228. ************************************************************************
  1229. *** 1ERE RESOLUTION ***
  1230. ************************************************************************
  1231. 'SI' ('EXIS' STAB12 'RIBLO_M' ) ;
  1232. 'SI' ( WTAB. 'CAFROTTE' );
  1233. ZDEP1 BID BID BID BID = 'RESO' ZRAID RESIDU 'INIB'
  1234. STAB12.'RIBLO_M'
  1235. STAB12.'LISEA_M' FFROT;
  1236. 'SINON';
  1237. ZDEP1 BID BID BID = 'RESO' ZRAID RESIDU 'INIB'
  1238. STAB12.'RIBLO_M'
  1239. STAB12.'LISEA_M' ;
  1240. 'FINS';
  1241. 'OUBLIER' STAB12.'RIBLO_M';
  1242. 'SINON';
  1243. 'SI' (WTAB. 'CAFROTTE');
  1244. ZDEP1 = 'RESO' ZRAID RESIDU FFROT;
  1245. 'SINON';
  1246. ZDEP1 = 'RESO' ZRAID RESIDU;
  1247. 'FINS';
  1248. 'FINS';
  1249.  
  1250. ZRAID_T = 'EXTR' ZRAID 'CONT';
  1251.  
  1252. 'SI' IDYN;
  1253. STAB12.'ZRAIDV'= ZRAID;
  1254. 'FINS';
  1255.  
  1256. * 'SI' ('EXIS' ZRAID_T 'NITER');
  1257. * STAB12.'RIBLO_M' = ZRAID_T. 7 ;
  1258. * STAB12.'LISEA_M' = ZRAID_T. 6 ;
  1259. * 'FINS';
  1260. * RENORMALISATION de du (ajout bp le 28/11/2012) ****************
  1261. 'SI' (exis WTAB 'RENORMALISATION');
  1262. 'SI' (WTAB . 'RENORMALISATION');
  1263. * juste apres le resou,
  1264. * on calcule et on limite ZDEP1 de l iteration
  1265. * en limitant DEPS1 par MAXDEFOR
  1266. coefmul = 1.;
  1267. DEPS1 = EPSI 'LINE' ZMODL ZDEP1 ZMAT;
  1268. DPS1MAX = MAXI 'ABS' DEPS1 'AVEC' MLDEFOR ;
  1269. 'SI' (DPS1MAX > WTAB . 'MAXDEFOR');
  1270. coefmul = WTAB . 'MAXDEFOR' / DPS1MAX;
  1271. ZDEP1 = coefmul * ZDEP1;
  1272. mess 'MAXDEFOR dépassé : Renormalisation initiale par ' coefmul;
  1273. 'FINS';
  1274. 'FINS';
  1275. 'FINS';
  1276. * On sauve le deplacement initial pour la convergence forcee ***********
  1277. 'SI' ('EXIS' WTAB 'DEPI');
  1278. zdeptini = WTAB.'DEPI';
  1279. 'SINON';
  1280. zdeptini = zdep1;
  1281. 'FINS';
  1282. *
  1283. * calcul d'une norme pour la convergence**************************
  1284. *
  1285. XXX1= ZFEXT;
  1286. 'SI' ADDISEC2;
  1287. XXX1=ZFEXT + FP22;
  1288. 'FINS';
  1289. 'SI' (LOGPRE 'ET' IGRD) ;
  1290. XXX1 = ZFEXT + ZFPEXTF ;
  1291. 'FINS';
  1292.  
  1293. ZDEP1P50 = ZDEP1 ;
  1294. XDENO='XTY' ZDEP1P50 ( XXX1 -( RESIDU 'EXCO'
  1295. 'FLX' 'NOID' 'FLX' 'NATURE' 'DISCRET')) MLPRIM MLDUAL;
  1296. MZDEP1M = 'MAXI' ZDEP1P50 'ABS' 'AVEC' MLDEPL;
  1297. MZFM = 'MAXI' (FTHE + FDEF + ZF1) 'ABS' 'AVEC' MLDUAL;
  1298. XDENO1 = 'ABS' XDENO + (MZFM * MZDEP1M);
  1299. MZDEP1M = MZDEP1M + XPETIT;
  1300. XDENO=XDENO1/MZDEP1M;
  1301. XDENO = XDENO + MZFM;
  1302. XDENO = XDENO + XPETIT ;
  1303. XDENOM=XDENO;
  1304. 'SI' TSTMOM ;
  1305. ZDEP1P50 = ZDEP1 + XPETIT; ;
  1306. XDENOM = XDENO1/('MAXI' ZDEP1P50 'ABS' 'AVEC' MLROTA);
  1307. XDENOM = XDENOM + XPETIT ;
  1308. 'FINS' ;
  1309. 'SI' IPILOT ;
  1310. 'SI' ( WTAB.'AUTODEUX' ) ;
  1311. XDENO = STAB12.'XDENO';
  1312. XDENOM = STAB12.'XDENOM';
  1313. 'FINS';
  1314. 'FINS';
  1315. 'SI' ('NON' WTAB.'CONV');
  1316. XDENO = STAB12.'XDENO';
  1317. XDENOM = STAB12.'XDENOM';
  1318. 'FINS';
  1319.  
  1320. IAFAIR=FAUX;
  1321. 'SI' WTAB.'CONV';
  1322. STAB12.'INCREMENT' = 'COPIER' ZDFINI ; 'FINS';
  1323. RESIDNOR = 'COPIER' RESIDU ;
  1324. 'SI' IPILOT;
  1325. XXX3=DFEXT0F * ( 1-COEPI) ;
  1326. XXX2 = ZFEXT - XXX3;'DETR' XXX3;
  1327. 'DETR' ZFEXT; ;ZFEXT = XXX2;
  1328. XXX1 = RESIDU * COEPI; 'DETR' RESIDU;
  1329. XXX2 = 1.D0 -COEPI * FREAP;RESIDU = XXX1 + XXX2;
  1330. 'SI' IMPO12;
  1331. RESIDU = 1.D0 - COEPI * DIMPO12 + RESIDU;
  1332. 'FINS';
  1333. IAFAIR=VRAI;
  1334. 'FINS';
  1335. *
  1336. * petite correction du residu pour esperer gagner du temps ************
  1337. *
  1338. INIT = FAUX ;
  1339. 'SI'(('NEG' STAB12.'FNONL' 'INCONNU') 'ET'
  1340. (WTAB.'INITIALISATION'));
  1341. 'SI' IPILOT;
  1342. 'SI' ( WTAB.'AUTODEUX' 'ET' (COEPI 'NEG' 1.D0)) ;
  1343. 'MESS' 'Initialisation a partir du pas precedent ' COEPI;
  1344. IAFAIR=VRAI; INIT = VRAI;
  1345. XXX1= RED1 * STAB12.'FNONL';
  1346. XXX2= XXX1 + RESIDU ;
  1347. 'DETR' RESIDU ;'DETR' XXX1 ;RESIDU =XXX2;
  1348. 'FINS' ;
  1349. 'SINON';
  1350. * on fait la correction si le pas precedent a converge
  1351. 'SI' WTAB.'CONV';
  1352. * on fait la correction si le pas precedent etait non lineaire.
  1353. 'SI' (('MAXI' 'ABS' STAB12.'FNONL') > (ZPREC * XDENO)) ;
  1354. * on enleve le residu du pas precedent pour recuperer l'increment
  1355. * nominal du second membre e imposer f1.
  1356. * f2 et l'increment du second membre du pas precedent
  1357. STAB12.'INCREMENT'=STAB12.'INCREMENT' - STAB12.'RESIDU';
  1358. zdeps = WTAB.'ZDEP1' + zdep1;
  1359. FFNO= 'XTY' STAB12.'FNONL' zdeps MLDUAL MLPRIM;
  1360. F1 = STAB12.'INCREMENT' ;
  1361. F2 = INCRPREC ;
  1362. f12 = 'XTY' f1 zdeps MLDUAL MLPRIM;
  1363. f22 = 'XTY' f2 zdeps MLDUAL MLPRIM;
  1364. AMPL = f12/(f22 + XPETIT);
  1365. DTPREC=1E30; AMPLT = 0;
  1366. 'SI' ('NEG' STAB12.'DTPREC' 'INCONNU');
  1367. DTPREC = STAB12.'DTPREC';
  1368. 'SI' (DTPREC > XPETIT);
  1369. AMPLT=WTAB.'DT_INIT' /DTPREC;
  1370. 'FINS';
  1371. 'FINS';
  1372. * Le chargement n'est il pas de fluage ou de thermique ?';
  1373. XDCOMP = ('XTY' ZDEP1 F1 MLPRIM MLDUAL) ;
  1374. 'SI' (('ABS' XDCOMP) < (ZPREC * XDENO * mzdep1m)
  1375. 'OU' (('ABS' FFNO) > (('ABS' F22) * 2.e2 )));
  1376. DTPREC = STAB12.'DTPREC';
  1377. AMPL=AMPLT;
  1378. * la decharge est-elle significative
  1379. 'SI'((F12/(f22 + XPETIT)) < -0.05);ampl=0.;'FINS';
  1380. STAB12.'INITEMPS'=VRAI;LOGTEMP=FAUX;;
  1381. 'MESS'
  1382. 'Pas d increment de charge, initialisation calculee avec le temps';
  1383. 'SINON';
  1384. AMPL = F12 / (F22 + XPETIT);
  1385. LOGTEMP=STAB12.'INITEMPS';
  1386. STAB12.'INITEMPS'=FAUX;
  1387. 'FINS';
  1388. AMPL = MINI (prog AMPL AMPLT);
  1389. 'SI' ((AMPL > 0) 'ET' (AMPL < 2e1)) ;
  1390. 'MESS' 'Initialisation a partir de la solution precedente Coeff'
  1391. AMPL;
  1392. XXX1 = AMPL * STAB12.'FNONL';
  1393. XXX2 =RESIDU+ XXX1;
  1394. 'DETR' XXX1;'DETR' RESIDU;
  1395. RESIDU = XXX2;
  1396. IAFAIR=VRAI; INIT = VRAI;
  1397. 'FINS';
  1398. 'FINS';
  1399. 'FINS';
  1400. 'FINS';
  1401. 'FINS';
  1402. WTAB.'ZDEP1'=zdep1;
  1403. *
  1404. * initialisation en plastique *****************************************
  1405. *
  1406. 'SI' IPLAVI ;
  1407. 'SI' ('NEG' WTAB.'MOVA' 'RIEN') ;
  1408. ACC0 = 'EXCO' (WTAB.'MOVA') ZVAR0 ;
  1409. 'FINS' ;
  1410. 'FINS' ;
  1411. *
  1412. * debut des iterations internes boucle etiquette ********************
  1413. *
  1414. IKT = FAUX ;
  1415. IPREM = VRAI ; RECA_K = FAUX; RECA_N = 0;
  1416. DEPSTP = 'ZERO' ZMODL 'DEFORMATIONS' ;
  1417. DEPSTK = DEPSTP ;
  1418. *
  1419. * initialisation acceleration de convergence ********************
  1420. *
  1421. iafair = vrai;
  1422. PASTEST = FAUX;
  1423. ZITAC= 0 ;
  1424. * on peut mettre n'importe quoi c'est pour
  1425. * ne pas faire de tests dans la boucle
  1426. ACFP1 = 'COPIER' ZFEXT2 *0. ;
  1427. ACFP2 = ACFP1 ;
  1428. ZDEPLD = ACFP1 ;
  1429. ACFP3 = ACFP1 ;
  1430. ACFEP1 = ACFP1 ;
  1431. ACFEP2 = ACFP1 ;
  1432. FCORF = 'COPIER' FREAP;
  1433. * initialisation du meilleur critere
  1434. XCONVMIN = 1e20;
  1435. DPSMREF = 0 ;
  1436. XCONV = 0. ;
  1437. XCONVP = 1. ;
  1438. NSOINCR = 1 ;
  1439. 'SI' ('NEG' WTAB.'SOUS_INCREMENT' 'INCONNU') ;
  1440. NSOINCR = WTAB.'SOUS_INCREMENT' ;
  1441. 'FINS';
  1442. NONCONV = FAUX;
  1443. CORREC=0;
  1444. PASREINI=VRAI;
  1445. coefmt=1.;
  1446. *
  1447. * initialisation des messages pour l'iteration en cours ****************
  1448. *
  1449. 'SI' IKLFFF ;
  1450. 'MESSAGE'
  1451. ' Iter Nplas Critere Deps.max Eps.max Crit.flex'
  1452. ;
  1453. 'SINON';
  1454. 'MESSAGE'
  1455. ' Iter Nplas Fresidu Deps.max Eps.max Mresidu '
  1456. ;
  1457. 'FINS';
  1458.  
  1459.  
  1460.  
  1461. *=======================================================================
  1462. *======= DEBUT DE LA BOUCLE DE CONVERGENCE =======
  1463. *=======================================================================
  1464. 'REPETER' ETIQ ;
  1465. *IT est le compteur de ETIQ, ITACC doit etre =< 0 pour qu'on accelere
  1466. IT= IT + 1 ;
  1467. ITACC = ITACC - 1;
  1468. ZITAC = ZITAC + 1 ;
  1469. *
  1470. *---------------------------------------------------------------------
  1471. * La force motrice de l'iteration est fixee: RESIDU
  1472. * on va calculer un nouveau champ de deplacement
  1473. *
  1474. * calcul de l'increment de l'increment de deplacement zdep1
  1475. * par resolution lineaire
  1476. *-----------------------------------------------------------------------
  1477. *
  1478. * petits travaux pour acceleration de convergence
  1479. *
  1480. CORRECP = CORREC;
  1481. CORREC = 0;
  1482. PASTEST=FAUX;
  1483. ACFP0 = (RESIDU - FCORF) 'ENLE' FLX ;
  1484. ACFEP0 = ACFP0;
  1485. ACFEP0 = ACFEP0 - CORRECP ;
  1486. 'SI' IGRD ;
  1487. 'FORM' GEOM1;
  1488. * si (it>2 et (Depsi>0.1 ou URG ou itacc>4))
  1489. * ou si on l'impose via nrecalc (FAUX ici)...
  1490. * rem : 0.1 est la valeur par defaut de EKREAC (= . 'REAC_GRANDS')
  1491. 'SI' ((&etiq '>EG' 2) 'ET' ((DEKREAC1 > EKREAC)
  1492. 'OU' URG 'OU'(ITACC >EG 4)) 'OU' (FAUX ));
  1493. URG=FAUX;
  1494. sigttcca= zsigf ;
  1495. geor='FORM' zdept;
  1496. zdepr = zdept;
  1497. c_zdepr = faux;
  1498. RECA_K = VRAI;
  1499. RECA_N = RECA_N + 1;
  1500.  
  1501. 'SI' (RECA_N > 20) ;
  1502. nonconv = vrai;
  1503. 'FINS';
  1504.  
  1505. txt_k ='CHAI' ' Recalcul de K (= K^el';
  1506. ritc ='RIGI' zmodl zmat ;
  1507. * RH peut contenir des CL
  1508. ZRI = 'EXTR' RITC 'RIGI' 'NOMU' ;
  1509. ZCL = 'EXTR' RITC 'RIGI' 'MULT' ;
  1510. 'SI' ('NEG' ('DIME' ZCL) 0) ;
  1511. ZCLIM = ZCLIM0 'ET' ZCL ;
  1512. 'SINON' ;
  1513. ZCLIM = ZCLIM0 ;
  1514. 'FINSI' ;
  1515.  
  1516. BFCONT = FAUX ;
  1517. BCLIM2 = FAUX ;
  1518. 'SI' WTAB.'CONTACT';
  1519. * recalcul raideur de contact et jeu
  1520. MODCON = wtab.'MODCONTA';
  1521. WTAB.'CONT_REEL'=MODCON;
  1522. 'SI' ('EXIS' WTAB 'MAT_JEU');
  1523. CDAP CRR RFROT='RFCO' MODCON WTAB.'CONV' WTAB.'MAT_JEU';
  1524. 'SINON' ;
  1525. CDAP CRR RFROT='RFCO' MODCON WTAB.'CONV' ;
  1526. 'FINS' ;
  1527.  
  1528. 'SI' (WTAB . 'MODAL') ;
  1529. CRR = 'PJBA' CRR ZMODL MMMM ;
  1530. MCRR = 'EXTR' CRR 'MAIL' 'MULT' ;
  1531. MCDAP = 'EXTR' CDAP 'MAIL' ;
  1532.  
  1533. 'REPETER' BCDA ('NBNO' MCDAP) ;
  1534. PBCDA = MCDAP 'POINT' &BCDA ;
  1535. PCRR ='POINT' MCRR &BCDA ;
  1536. CHCR ='MANU' 'CHPO' PCRR 1 'FLX' ('EXTR' CDAP 'FLX' PBCDA)
  1537. 'NATURE' 'DISCRETE' ;
  1538. 'SI' ('EGA' 1 &BCDA) ;
  1539. CCDA = CHCR ;
  1540. 'SINON' ;
  1541. CCDA = CCDA 'ET' CHCR ;
  1542. 'FINS' ;
  1543. 'FIN' BCDA ;
  1544.  
  1545. CDAP = CCDA ;
  1546. 'FINS' ;
  1547.  
  1548. * ATTENTION : mettre les conditions de frottement en premier pour
  1549. * les numeroter en dernier
  1550. ZCLIM2 ='VIDE' 'RIGIDITE' ;
  1551. ZFCONT ='VIDE' 'CHPOINT' / 'DISCRET';
  1552. 'SI' ('NEG' CRR 0);
  1553. BCLIM2 = VRAI ;
  1554. BFCONT = VRAI ;
  1555. ZCLIM2 = CRR 'ET' ZCLIM2 ;
  1556. CCOR = CRR '*' ((zdept '+' zu1) 'ENLE' 'LX');
  1557. ZFCONT = ZFCONT '+' CCOR ;
  1558. 'FINS';
  1559.  
  1560. 'SI' ( 'NEG' 0 CDAP);
  1561. BFCONT = VRAI ;
  1562. ZFCONT = ZFCONT '+' CDAP ;
  1563. 'FINS';
  1564.  
  1565. 'SI' ('NEG' 0 RFROT);
  1566. BCLIM2 = VRAI ;
  1567. BFCONT = VRAI ;
  1568. ZCLIM2 = RFROT 'ET' ZCLIM2 ;
  1569. CCOR = RFROT '*' CDEPSLX ;
  1570. ZFCONT = ZFCONT '+' CCOR ;
  1571. 'FINS';
  1572.  
  1573. * Mise a jour de ZCLIM, ZFCONSTA et RESIDU si necessaire
  1574. 'SI' BCLIM2;
  1575. * Calcul des reactions totales AVANT
  1576. freacav= ZCLIM '*' zdetot 'ENLE' 'FLX';
  1577. ZCLIM = ZCLIM2 'ET' ZCLIM ;
  1578.  
  1579. * Calcul des reactions totales APRES
  1580. freacap= ZCLIM '*' zdetot 'ENLE' 'FLX';
  1581. RESIDU = RESIDU '+' freacav '-' freacap;
  1582. 'FINS';
  1583.  
  1584. 'FINS';
  1585.  
  1586. 'SI' BFCONT;
  1587. ZFCONSTA = ZFEXT2 '+' ZFCONT;
  1588. 'SINO';
  1589. ZFCONSTA = ZFEXT2 ;
  1590. 'FINS';
  1591. ZFCONST1 = ZFCONSTA;
  1592.  
  1593. 'SI' IDYN;
  1594. ZFCONSTA = ZFCONSTA '+' ZFP1 ;
  1595. 'FINS';
  1596.  
  1597. ZFLX1 = 'EXCO' ZFCONSTA 'FLX' 'NOID' 'FLX' 'NATURE' 'DISCRET';
  1598. ZRAID = ZCLIM 'ET' ZRI;
  1599. 'SI' iksia ;
  1600. ksigtc= 'KSIGM' sigttcca zmodl zmat;
  1601. ZRAID= ZRAID 'ET' ksigtc;
  1602. 'FINS';
  1603. 'SI' IRCON;
  1604. ZRAID = ZRAID 'ET' RIG_CONS;
  1605. txt_k = 'CHAI' txt_k ' + K^cst';
  1606. 'FINS';
  1607. 'SI' IRAUG;
  1608. ZRAID = ZRAID 'ET' RIG_AUG;
  1609. txt_k = 'CHAI' txt_k ' + K^aug';
  1610. 'FINS';
  1611.  
  1612. 'SI' WTAB.'PROCEDURE_CHARMECA';
  1613. PRECED . 'ADDI_MATRICE' = vrai;
  1614. TFP22= CHARMECA PRECED WTAB.'T_FINAL';
  1615. PRECED . 'ADDI_MATRICE' = faux;
  1616. 'SI' (EXIS TFP22 'ADDI_MATRICE');
  1617. zraid = zraid 'ET' TFP22.'ADDI_MATRICE';
  1618. txt_k = 'CHAI' txt_k ' + K^cent';
  1619. 'FINS';
  1620. 'FINS';
  1621.  
  1622. 'SI' IDYN;
  1623. * bp : en toute rigueur, il faudrait aussi recalculer la MASSE ...
  1624. * et ajouter l'amortissement le cas échéant ...
  1625. MMA = WTAB.'MASSE' '*' ( 4.D0 '/' WTAB.'DT' '/' WTAB.'DT');
  1626. ZRAID= ZRAID 'ET' MMA;
  1627. 'FINS';
  1628. 'MESS' ( 'CHAI' txt_k ' ) dans config deformee ' DEKREAC1 );
  1629. 'FORM' GEOM1;
  1630. 'DETR' GEOR;
  1631. * on impose le recalcul de K a la prochaine iteration si it=2
  1632. 'SI' (IT 'EGA' 2) ;
  1633. ITACC = 5 ;
  1634. 'SINO' ;
  1635. ITACC = 3 ;
  1636. 'FINS' ;
  1637. GR_U_K = GR_U_FIN ;
  1638. zdeptm = zdept ;
  1639. 'FINS';
  1640. 'FINS';
  1641.  
  1642. *
  1643. * acceleration de convergence effective ------------------------------
  1644. *
  1645. 'SI' (('NON' (NONCONV 'ET' ZICONV)) 'ET' ('MULT' IT ZNACCE)
  1646. 'ET' (ITACC '&lt;EG' 0) 'ET' (IT '>' 3) );
  1647. CORREC = 'ACT3' ACFEP2 ACFEP1 ACFEP0
  1648. ACFP3 ACFP2 ACFP1 ACFP0 ;
  1649. RESIDU = RESIDU '-' CORREC;
  1650. 'FINS';
  1651.  
  1652. ACFP3 = ACFP2 ;
  1653. ACFP2 = ACFP1 ;
  1654. ACFP1 = ACFP0 ;
  1655. ACFEP2 = ACFEP1 ;
  1656. ACFEP1 = ACFEP0 ;
  1657. *
  1658. * Resolution ---------------------------------------------------------
  1659. * on obtient ZDEP1 = [ du ; Dlx ]
  1660. 'SI' IAFAIR;
  1661. 'SI' ( ('EXIS' STAB12 'RIBLO_M') 'ET' (IPREM 'OU' RECA_K) ) ;
  1662. 'SI' ( WTAB.'CAFROTTE' );
  1663. ZDEP1 BID BID BID BID = 'RESO' ZRAID RESIDU 'INIB'
  1664. STAB12.'RIBLO_M' STAB12.'LISEA_M' FFROT;
  1665. 'SINON';
  1666. ZDEP1 BID BID BID = 'RESO' ZRAID RESIDU 'INIB'
  1667. STAB12.'RIBLO_M' STAB12.'LISEA_M' ;
  1668. 'FINS';
  1669.  
  1670. 'SINON';
  1671. 'SI' ( WTAB.'CAFROTTE' );
  1672. ZDEP1 = 'RESO' ZRAID RESIDU FFROT ;
  1673. 'SINON';
  1674. ZDEP1 = 'RESO' ZRAID RESIDU;
  1675. 'SI' (WTAB.'MAN' 'ET' IPREM);
  1676. ORDRE = WTAB.'ORDRE' ;
  1677. ZDEP2 IOUT = CORMAN ZRAIDINI ZMODL ZMAT ORDRE
  1678. ZU1 ZSIG0 RESIDNOR WTAB ;
  1679. 'SI' (IOUT 'EGA' 1) ;
  1680. ZDEP1=ZDEP2;
  1681. 'FINS';
  1682. 'FINS';
  1683. 'FINS';
  1684. 'FINS';
  1685.  
  1686.  
  1687. ZRAID_T = 'EXTR' ZRAID 'CONT';
  1688. 'SI' ('EXIS' ZRAID_T 'NITER');
  1689. 'SI' IDYN ; STAB12.'ZRAIDV'= ZRAID; 'FINS';
  1690. STAB12.'RIBLO_M' = ZRAID_T.7;
  1691. STAB12.'LISEA_M' = ZRAID_T.6;
  1692. 'FINS';
  1693. 'FINS';
  1694. RECA_K = FAUX;
  1695. * RENORMALISATION de du (ajout bp le 28/11/2012) ***********************
  1696. 'SI' ('EXIS' WTAB 'RENORMALISATION');
  1697. 'SI' (WTAB . 'RENORMALISATION');
  1698. * juste apres le resou,
  1699. * on calcule et on limite ZDEP1 de l iteration
  1700. * en limitant DEPS1 par MAXDEFOR
  1701. coefmul = 1.;
  1702. DEPS1 = EPSI 'LINE' ZMODL ZDEP1 ZMAT;
  1703. DPS1MAX = MAXI 'ABS' DEPS1 'AVEC' MLDEFOR ;
  1704. 'SI' (DPS1MAX > WTAB . 'MAXDEFOR');
  1705. coefmul = WTAB . 'MAXDEFOR' / DPS1MAX;
  1706. 'SI' IPREM;
  1707. ZDEP1 = coefmul * ZDEP1 ;
  1708. 'SINO';
  1709. ZDLX0 = ZDEPTP 'EXCO' 'LX' 'NOID' 'LX';
  1710. ZDEP1 = (coefmul * ZDEP1) + ((1.-coefmul) * ZDLX0);
  1711. 'DETR' ZDLX0;
  1712. 'FINS';
  1713. 'MESS' 'MAXDEFOR dépassé : Renormalisation de l iteration par '
  1714. coefmul;
  1715. 'FINS';
  1716. 'FINS';
  1717. 'FINS';
  1718. 'REPETER' bzdept2 7;
  1719. *
  1720. * verif sens ----------------------------------------------------------
  1721. *
  1722. 'SI' WTAB.'STABILITE';
  1723. sens = xty ((residu enle 'FLX') +
  1724. (-1 * residu exco 'FLX' 'NOID' 'FLX')) zdep1 MLDUAL MLPRIM;
  1725. nbneg = diagn zraid;
  1726. 'SI' (( nbneg > 0));
  1727. ZDEP1 = (zdep1 enle lx) * -0.2 et (zdep1 exco 'LX' 'NOID' 'LX') ;
  1728. 'MESS' 'instabilite detectee' nbneg sens;
  1729. HPP_EPS = FAUX;
  1730. 'FINS';
  1731. 'FINS';
  1732. *
  1733. * 1ere iteration ------------------------------------------------------
  1734. 'SI' IPREM ;
  1735. ZDEPT = 'COPIER' ZDEP1 ;
  1736. ZDEPTM = ZDEPT ;
  1737. ZDEPTP = ZDEPT ;
  1738. ZDELA = 'COPIER' ZDEPT ;
  1739. 'SI' (WTAB . 'MODAL') ;
  1740. 'SI' ('EXIS' WTAB 'MODCONTA') ;
  1741. mamoco1 = 'EXTR' (ZDEPT 'ENLEVER' 'LX') 'MAIL' ;
  1742. 'FINS';
  1743. 'FINS' ;
  1744. * iprem est faux: on est apres la premiere operation---------------
  1745. *
  1746. 'SINON';
  1747. * zdept est l'increment de deplacement total avec les lagrangiens
  1748. * de la solution complete
  1749. XXX1 = ZDEPT 'ENLEVER' 'LX' ;
  1750. zdeptp = ZDEPT ;
  1751. zsigfp = zsigf ;
  1752. fcorfp = fcorf ;
  1753. zclimp = zclim ;
  1754. * on cumule les deplacements mais pas les lx
  1755. * Du^(i) = Du^(i-1) + du ; Dlx^(i) = 0 + Dlx
  1756.  
  1757. ZDEPT = XXX1 '+' ZDEP1 ;
  1758. 'DETR' XXX1 ;
  1759. 'FINS' ;
  1760. *
  1761. 'SI' (WTAB . 'MODAL') ;
  1762. 'SI' ('EXIS' WTAB 'MODCONTA') ;
  1763. * mettre les point materiels dans zdept
  1764. zdeptu1 = 'REDU' zdept mamoco1 ;
  1765. ch_dco = 'RECO' zdeptu1 ZMODL ZMATFI ;
  1766. ZDEPT = ('EXCO' zdept 'LX' 'LX') 'ET' zdeptu1 'ET' ch_dco ;
  1767. 'FINS';
  1768. 'FINS' ;
  1769. * boucle de reduction de zdept si necessaire
  1770.  
  1771. ACC_Q = FAUX;
  1772. ACC_R = VRAI;
  1773. 'REPETER' bzdept 40;
  1774. *
  1775. *--- CAS 1 ------------------------------------------------------------*
  1776. * Limitation du maximum d'increment de deformation (DPSMAX) : *
  1777. * En non-convergence, on limite l'increment de deformation pour faci- *
  1778. * -liter l'integration du comportement. *
  1779. * En convergence normale, on limite l'increment de deformation max. *
  1780. * entre 2 iterations pour tenter d'eviter de partir dans le decors *
  1781. * en grands deplacement (Jacobien negatifs, etc...) *
  1782. * nul si mouvement corps rigide (contact...) : *
  1783. 'SI' ((IT '>' 5 ) 'ET' (DPSMAX '>' 0.) 'ET' (&BZDEPT 'EGA' 1)) ;
  1784. coefmul = WTAB.maxdefor / (DPSMAX + XPETIT) ;
  1785. * cas ou le residu augmente...
  1786. 'SI' (XCONV '>' TABCONV.(IT '-' 2) ) ;
  1787. * on n'a pas convergé ET DPSMAX trop grand : MAXDEFOR/zprec < DPSMAX
  1788. * ou trop petit : DPSMAX < 1E-4*zprec
  1789. 'SI' ((NONCONV 'ET' ((coefmul '<' zprec) 'OU' (dpsmax*1D4 < zprec))));
  1790. coefmul = coefmul '*' 1D-1;
  1791. ZDEPT = zdept '*' coefmul ;
  1792. 'SI' ((RED_URG '>' 1) 'ET' NONCONV);
  1793. zflx1 = flxini '+' ((zflx1 '-' flxini) '*' 0.5 );
  1794. 'FINS';
  1795. RED_URG = RED_URG '+' 1;
  1796. ZNACCE = 999;
  1797. 'MESS' ' Renormalisation de l increment de deplacement'
  1798. dpsmax coefmul;
  1799. 'SI' (DEPSTDM '>' DEPSTREF) ;
  1800. DEPSTREF = 2. * DEPSTDM ;
  1801. 'FINS' ;
  1802. * Est-ce qu'on cycle des renormalisations en Non-convergence ? *
  1803. * On regarde le nombre d'iterations entre 2 renormalisations *
  1804. * & la valeur de XCONV. *
  1805. 'SI' NONCONV ;
  1806. 'SI' (ITNORM1 'EGA' 0) ;
  1807. ITNORM1 = IT ;
  1808. 'SINO' ;
  1809. 'SI' (DITNORM1 'EGA' 0) ;
  1810. XCONVNOR = 2.D0 '*' XCONV ;
  1811. DITNORM1 = IT '-' ITNORM1 ;
  1812. ITNORM1 = IT ;
  1813. 'SINO' ;
  1814. 'SI' (XCONV '&lt;EG' XCONVNOR) ;
  1815. XCONVNOR = 2.D0 '*' XCONV ;
  1816. 'SI' ((IT - ITNORM1) 'EGA' DITNORM1) ;
  1817. NBCYCLE1 = NBCYCLE1 '+' 1 ;
  1818. 'SI' (NBCYCLE1 'EGA' 2) ;
  1819. 'MESS'
  1820. ' Renormalisation cyclique en non-convergence detectee ' ;
  1821. 'FINS' ;
  1822. 'FINS' ;
  1823. DITNORM1 = IT '-' ITNORM1 ;
  1824. ITNORM1 = IT ;
  1825. 'FINS' ;
  1826. 'FINS' ;
  1827. 'FINS' ;
  1828. 'FINS' ;
  1829. * 'SI' (coefmtt < 1.);
  1830. **pv zflx1 = flxini + ((zflx1 - flxini) * coefmul);
  1831. * coefmt=coefmtt;
  1832. * 'FINS';
  1833. 'FINS' ;
  1834. 'FINS' ;
  1835. * fin du cas ou le residu augmente...
  1836. 'FINS';
  1837.  
  1838. *
  1839. *--- CAS 2 ------------------------------------------------------------*
  1840. * si on part dans les decors on redemarre a 0
  1841. 'SI' ((XCONV > 1E8) 'ET' PASREINI 'ET' ('NON' IPILOT)) ;
  1842. 'MESS' 'Reinitialisation du schema';
  1843. XXX1 = ZDEPT 'ENLEVER' 'LX' ;
  1844. ZDEPT = ZDEPT '-' XXX1 ;
  1845. PASREINI=FAUX;
  1846. ITACC=3;
  1847. 'FINS';
  1848.  
  1849. *
  1850. *--- CAS 3 ------------------------------------------------------------*
  1851. 'SI' (('MULT' IT ZNACCE) 'ET' (ITACC '&lt;EG' 0)) ;
  1852.  
  1853. * pilotage automatique
  1854. 'SI' IPILOT;
  1855. RED1 = STAB12.'AUTOREDU' ;
  1856.  
  1857. * reduction du critere de pilotage red2 sert d'incateur si
  1858. * pas accelere tous les 2 pas
  1859.  
  1860. * MODIFICATION CB215821 : 18/06/2015
  1861. 'SI' (RED2 < (IT '/' 20));
  1862. STAB12.'AUTOREDU' = STAB12.'AUTOREDU' '*' 3.D0 ;
  1863. RED2 = RED2 '+' 1 ;
  1864. ITACC= 3;
  1865. 'FINS';
  1866.  
  1867. RED1 = RED1 '/' (STAB12.'AUTOREDU') ;
  1868. 'SI' (RED1 '<' 0.9) ; STAB12.AUTORED1 = 10 ; 'FINS' ;
  1869.  
  1870. 'SI' ( STAB12.'AUTOREDU' '>' 1.D0 ) ;
  1871. 'MESS' 'On divise le critere de pilotage par '
  1872. STAB12.'AUTOREDU';
  1873. 'FINS' ;
  1874.  
  1875. * test si snap back et si refus de l'acceleration
  1876. 'SI' IPLAVI ;
  1877. XXX2 = 'EXCO' ZDEPT 'LX' 'NOID' 'LX' ;
  1878. XXX1 = 2.D0 '*' XXX2;
  1879. XXX3 = ZDEPT '-' XXX1;
  1880. 'SINON' ;
  1881. XXX1 = COEPI '*' ZDELA;
  1882. XXX3 = ZDEPT '-' XXX1 ;
  1883. XXX2 = 'EXCO' XXX3 'LX' 'NOID' 'LX';
  1884. XXX1 = XXX2 '*' 2.D0 ;
  1885. XXX4 = COEPI '*' ZDELA;
  1886. XXX3 = ZDEPT '-' XXX4 '-' XXX1;
  1887. 'FINS' ;
  1888. SRE = 'XTY' ZDFINI XXX3 MLDUAL MLPRIM ;
  1889.  
  1890. 'DETR' XXX2 ;'DETR' XXX1; 'DETR' XXX3;
  1891. XXX1 ='EXCO' ZDEPT 'LX' 'NOID' 'LX' ;
  1892. XXX2 = XXX1 '*' 2.D0 ;
  1893. XXX3 = ZDEPT '-' XXX2;
  1894. SRT = 'XTY' ZDFINI XXX3 MLDUAL MLPRIM;
  1895. 'DETR' XXX1 ; 'DETR' XXX2; 'DETR' XXX3;
  1896.  
  1897. ISNPB = FAUX;
  1898. 'SI' (SRT '<' 0) ;
  1899. ISNPB = VRAI ;
  1900. 'FINS';
  1901.  
  1902. OO = WTAB.AUTOCRIT ;
  1903. OO = OO '/' STAB12.'AUTOREDU' ;
  1904. U1MA = AUTOPILO ZDEPT (COEPI'*'ZDELA) ZMODLI ZMATFI WTAB;
  1905. AL1 = OO '/' U1MA ;
  1906. * al1 est le coefficient de normalisation
  1907. 'SI' (( al1 '>EG' 1.d0) 'ET' (coepi '>EG' 1.d0)) ;
  1908. * pour eviter d'aller au dela de alpha=1 on ignore le critere
  1909. AL1 = 1.d0 ;
  1910. 'FINS' ;
  1911. 'SI'((AL1 '>' 1.D0) 'ET' (COEPI '>' 0.D0));
  1912. 'SI' ( AL1 '>' (1.D0 '/' COEPI));
  1913. AL1 = 1.D0 '/' COEPI;
  1914. 'FINS';
  1915. 'FINS';
  1916. * normalisation
  1917.  
  1918. PASTEST=VRAI;
  1919. XXX1 = ZLX '*' ( 1.d0 '-' AL1);
  1920. XXX3 = ZDEPT '*' AL1 ;
  1921. ZDEPT = XXX3 '+' XXX1;
  1922. 'DETR' XXX3;
  1923. 'FINS';
  1924. 'SINON';
  1925. 'SI' IPILOT ;
  1926. **
  1927. * le calcul d'al1 est fait pour eviter la convergence
  1928. * s'il est externe a l'interval 0.5 --- 2.
  1929. OO = WTAB.'AUTOCRIT' ;
  1930. OO = OO '/' STAB12.'AUTOREDU' ;
  1931. U1MA = AUTOPILO ZDEPT (COEPI'*'ZDELA) ZMODLI ZMATFI WTAB;
  1932. AL1 = OO '/' U1MA ;
  1933. 'SI' (( AL1 '>EG' 1.d0) 'ET' (COEPI '>EG' 1.d0)) ;
  1934. AL1 = 1.d0 ;
  1935. 'FINS' ;
  1936. 'FINS' ;
  1937. 'FINS';
  1938. *
  1939. * garder les reactions pour le test de convergence
  1940. *
  1941. ZDEPLP = ZDEPL ;
  1942. ZDEPL = ZDEPT 'EXCO' 'LX' 'NOID' 'LX' ;
  1943.  
  1944. *----------------------------------------------------------------------
  1945. * le nouveau champ est fixe on va tester l'equilibre(convergence)
  1946. * et calculer la force motrice pour l'iteration suivante
  1947. *----------------------------------------------------------------------
  1948. *
  1949. * calcul de fcorf = lambda * m force de reaction -----------------
  1950. * fcoru = m * u depimp (flx)
  1951.  
  1952. 'DETRUIRE' FCORF;
  1953. * calcul du deplacement total u
  1954. ZDETOT = ZZD '+' ZDEPT ;
  1955. XXX1 = ZCLIM '*' ZDETOT;
  1956. * -1*reactions a l'iteration it
  1957. FCORF = 'ENLEVER' XXX1 'FLX';
  1958. * valeur des relations imposees a l'iteration it
  1959. FCORU = 'EXCO' XXX1 'FLX' 'NOID' 'FLX';
  1960. 'DETR' XXX1 ;
  1961. *
  1962. * dans le cas des modeles endommageables de Lemaitre, on ecoule
  1963. * en tenant compte, dans les iterations internes, de la variation du
  1964. * materiau avec la temperature
  1965. *
  1966. ZMATT = ZMAT ;
  1967. 'SI' ('ET' ('ET' ITHER ('OU' IENDOM IVIDOM ) ) WTAB.'MATVAR') ;
  1968. * on recupere certain materiau avec les parametres fct de la temperatue
  1969. * voir PAS_mate (il ne faut que la dependance thermique)
  1970. ZMATT = 'REDU' WTAB.'MA_COMP' ZMODL;
  1971. 'FINS';
  1972.  
  1973. 'SI' IFEFP;
  1974. * Update or total lagrangian -----------------------------------------
  1975. 'SI' IFEFPUL;
  1976. * mess ' update lagrangian ZRIKTA';
  1977. GEOM2 = 'FORM' ZDEPT ;
  1978. ZRIKTA ZSIGF ZVARF ZDEIF = 'ECFEFP'
  1979. ZMODL ZEPS0 ZVAR0 ZDEPT ZMAT ZPREK NITMA XUPDA;
  1980. 'SINON';
  1981. * mess ' total lagrangian ZRIKTA';
  1982. chp_z = ZDEPT + ZU1 ;
  1983. GEOM2 = 'FORM' chp_z ;
  1984. ZRIKTA ZSIGF ZVARF ZDEIF = 'ECFEFP'
  1985. ZMODL ZEPS0 ZVAR0 chp_z ZMAT ZPREK NITMA ;
  1986. chp_z = 1 ;
  1987. 'FINS';
  1988. FEQU2 = 'BSIGMA' ZMODL ZSIGF ZMAT ;
  1989. ZRAID = ZRIKTA 'ET' ZCLIM ;
  1990. 'FORM' GEOM1 ;
  1991. **
  1992. XXX1 = 'EXCO' (WTAB.'MOVA') ZVARF ;
  1993. EPSM = 'MAXI' XXX1 'AVEC' MLDEFOR ;
  1994. ACC = 'ABS' ( XXX1 - ACC0 ) ;
  1995. *'DETR' XXX1 ;
  1996. MMC = 'MASQUE' ACC 'SUPERIEUR' 1.D-10 'SOMME' ;
  1997. 'SI' (MMC '>' MMCMAX) ; MMCMAX = MMC ; 'FINS' ;
  1998. DPSMAXP = DPSMAX ;
  1999. DPSMAX = 'MAXI' ACC 'AVEC' MLDEFOR ;
  2000. DEPST = 'CHAN' 'TYPE' ZDEIF 'DEFORMATIONS' ;
  2001.  
  2002. GR_U_FIN= 'MOT' 'INCONNU';
  2003.  
  2004. 'SINON';
  2005. * cas standard --------------------------------------------------------
  2006. ZMAT05=ZMAT;
  2007. 'SI' IGRD;
  2008. 'FORM' GEOM1 ;
  2009. 'FINS';
  2010. *
  2011. * calcul de l'increment de deformation elast totale DEPST--------
  2012. *
  2013. 'SI' (EPS_NLIN 'ET' IGRD);
  2014. * (fdp) s'il y a des caracteristiques a mettre a jour, on le fait aussi
  2015. * sur les config intermediaires
  2016. zdepm = zdept '*' 0.5;
  2017. mzdepm = zdepm '*' -1 ;
  2018. 'SI' ITCAR;
  2019. GEOMIL ZMATTEMP = 'FORM' ZDEPM ZMODL ZMAT05 ;
  2020. DEPST = 'EPSI' 'LINE' 'NOER' ZMODL ZDEPT ZMATTEMP ;
  2021. 'DETR' ZMATTEMP ;
  2022. 'SINON';
  2023. GEOMIL = 'FORM' ZDEPM ;
  2024. DEPST = 'EPSI' 'LINE' 'NOER' ZMODL ZDEPT ZMAT05 ;
  2025. 'FINS';
  2026. PASOK = FAUX;
  2027.  
  2028. 'SI' ((NBPART '<' 1) 'OU' PARTLOCA 'OU' ('NON' PARALLEL));
  2029. 'SI' ('EGA' ('TYPE' DEPST) 'ENTIER');
  2030. PASOK = VRAI;
  2031. 'FINS';
  2032.  
  2033. 'SINON';
  2034. 'REPETER' BTEST NBPART;
  2035. 'SI' ('EGA' ('TYPE' (DEPST . &BTEST) ) 'ENTIER');
  2036. ** 'MESS' 'epsi trop grand ';
  2037. PASOK = VRAI;
  2038. 'FINS';
  2039. 'FIN' BTEST;
  2040. 'FINS';
  2041.  
  2042. 'SI' ('NON' PASOK);
  2043. depst = 'PICA' depst ZMODL mzdepm;
  2044. 'FORM' geom1;
  2045. 'SI' ('MAXI' 'ABS' DEPST '>' 2);
  2046. PASOK = VRAI;
  2047. 'FINS';
  2048. 'FINS';
  2049.  
  2050. 'SINON';
  2051. PASOK = vrai;
  2052. 'FINS';
  2053.  
  2054. ** on n'a pas calculé DEPST : on utilise les deformations lineaires
  2055. 'SI' PASOK;
  2056. * mess 'deformation lineaire';
  2057. 'SI' IGRD;
  2058. 'FORM' GEOM1 ;
  2059. 'FINS';
  2060. DEPST='EPSI' ZMODL ZDEPT ZMAT05 ;
  2061. 'FINS';
  2062.  
  2063. * calcul des increments de deformation et contrainte sur la configuration initia
  2064. * en lagrangien total
  2065. 'SI' IGRD;
  2066. 'SI' LAG_TOT;
  2067. * passage de depst et sig0 sur for0
  2068. depst = 'PICA' depst ZMODL mzu1 ;
  2069. zsig0 = 'PICA' zsig0 ZMODL mzu1 ;
  2070. form wtab.'FOR0';
  2071. 'FINS';
  2072. 'FINS';
  2073. * ICI depst et zsig0 sont sur la configuration initiale en lagrangien total et
  2074. * il faut soustraire les parties thermique et deformations imposees
  2075. 'SI' IPLAVI;
  2076. 'SI' ('OU' ('OU' IENDOM IVIDOM) ICERAM);
  2077. HOOKENDO = 'HOOK' ZMODL ZMAT05 ZVAR0 ;
  2078. 'SI' IGRD;
  2079. DEPST = 'EPSI' 'II' ZMODL ZDEPT HOOKENDO ZMAT05;
  2080. 'SINON';
  2081. DEPST = 'EPSI' 'LINE' ZMODL ZDEPT HOOKENDO ZMAT05;
  2082. 'FINS';
  2083. 'FINS';
  2084.  
  2085. * si thermoplastique on enleve e alpha *dt et d'autres termes
  2086. * si le materiau depend de la temperature
  2087. 'SI' (ITHER 'OU' WTAB.'MATVAR') ;
  2088. XXX1 = DEPST0 '*' COEPI;
  2089. XXX2 = DEPST '-' XXX1 ;
  2090. * 'DETR' DEPST ; 'DETR' XXX1;
  2091. DEPST = XXX2 ;
  2092. 'FINS' ;
  2093.  
  2094. 'SI' LOGDEF ;
  2095. XXX1 = DDEFOR0 '*' COEPI;
  2096. XXX2 = DEPST '-' XXX1 ;
  2097. * 'DETR' DEPST ; 'DETR' XXX1;
  2098. DEPST = XXX2 ;
  2099. 'FINS' ;
  2100.  
  2101. 'SINON';
  2102. * en elasticite, on a directement l'increment de contrainte
  2103. DSIGT= 'ELAS' ZMODL DEPST ZMAT05 ;
  2104. 'SI' (('NON' ITHER) 'ET' WTAB.'MATVAR') ;
  2105. XXXXX3 = 'ELAS' ZMODL ZSIG0 ZMAT05 ;
  2106. XXXXX4 = 'ELAS' ZMODL ZSIG0 ZMAT1 ;
  2107. XXDEFO = XXXXX4 '-' XXXXX3 ;
  2108. DSIGTV = 'ELAS' ZMODL XXDEFO ZMAT05 ;
  2109. DSIGT = DSIGT '+' DSIGTV ;
  2110. * 'DETR' XXXXX3 ; 'DETR' XXXXX4 ;
  2111. * 'DETR' XXDEFO ; 'DETR' DSIGTV ;
  2112. 'FINS' ;
  2113. * si thermoplastique on enleve e alpha *dt et d'autres termes
  2114. * si le materiau depend de la temperature
  2115. 'SI' ITHER ;
  2116. XXX1 = DSIGT0 '*' COEPI;
  2117. XXX2 = DSIGT '-' XXX1;
  2118. 'DETR' DSIGT ; 'DETR' XXX1;
  2119. DSIGT = XXX2 ;
  2120. 'FINS';
  2121.  
  2122. 'SI' LOGDEF ;
  2123. XXX1 = DSI1 '*' COEPI;
  2124. XXX2 = DSIGT '-' XXX1;
  2125. * 'DETR' DSIGT ; 'DETR' XXX1;
  2126. DSIGT = XXX2 ;
  2127. 'FINS';
  2128. 'FINS';
  2129. *
  2130. * Calcul du gradient des deplacements (total) a la fin du pas de ------
  2131. * temps par rapport a la configuration de reference GEOREF0 (initiale)
  2132. * rappel : On a : ZDETOT = ZZD + ZDEPT ;
  2133. 'SI' ('NEG' GR_U_DEB 'INCONNU');
  2134. 'SI' ('NON' LAG_TOT); 'FORM' WTAB.FOR0; 'FINS' ;
  2135. GR_U_FIN = 'GRAD' ZMODL ZMAT ZDETOT ;
  2136. D_GR_U = GR_U_FIN '-' GR_U_DEB ;
  2137. 'SI' ('NON' LAG_TOT); 'FORM' GEOM1 ; 'FINS';
  2138. 'SINON';
  2139. GR_U_FIN = 'MOT' 'INCONNU';
  2140. 'FINS';
  2141. *
  2142. 'SI' ZCCONV ;
  2143. STAB12.'DT' = DTINI '*' COEPI ;
  2144. 'SINON' ;
  2145. STAB12.'DT' = 0. ;
  2146. 'FINS' ;
  2147. *
  2148. * si on est plastique ou viscoplastique on ecoule pour avoir-----------
  2149. * le nouveau champ de contraintes
  2150. *
  2151. MMC = 0 ;
  2152. 'SI' IPLAVI ;
  2153.  
  2154. *...cas SSTE...
  2155. 'SI' ISSTE;
  2156. ZRIKTA ZSIGF ZVARF ZDEIF =
  2157. 'SSTE' ZMODL ZSIG0 ZVAR0 DEPST ZMATT
  2158. ZPREK NSSTE NITMA;
  2159. zvar0 = ('EXCO' zvar0 ('ENLE' lnom ('DIME' lnom))) 'ET'
  2160. ('EXCO' zvarf ssii);
  2161. zvar0 = 'CHANGER' 'TYPE' zvar0 'VARIABLES INTERNES';
  2162.  
  2163. *...cas non SSTE...
  2164. 'SINON';
  2165. ZSIG01=ZSIG0;
  2166. ZVAR01=ZVAR0;
  2167. ZEPS01=ZEPS0;
  2168. ** defineto = DEPIN0 ; TABCONT='TABLE';
  2169. CHASANST= WTAB.'CHARGEMENT' 'ENLE' 'T';
  2170.  
  2171. 'SI' (nsoincr 'NEG' 1);
  2172. ZSIGM=ZSIG0 '/' 2;
  2173. 'FINS';
  2174.  
  2175. 'REPETER' sousinc nsoincr;
  2176.  
  2177. tsodeb='FLOT' stab12 .'DT'/nsoincr*( &sousinc - 1)+conti.'TEMPS';
  2178. tsofin='FLOT' stab12 .'DT'/nsoincr* &sousinc +conti.'TEMPS';
  2179. 'SI' ('NON' ZCCONV) ;
  2180. tsodeb = tsofin ;
  2181. 'FINS' ;
  2182.  
  2183. che1= 'ADET' 'NOUV' ZMODL CHASANST tsodeb 'TEMP' tsodeb;
  2184. che2= 'ADET' 'NOUV' ZMODL CHASANST tsofin 'TEMP' tsofin;
  2185.  
  2186. 'SI' ('EXIS' WTAB 'MOD_LIA') ;
  2187. vite1 = 'CHAN' 'COMP' conti.'VITESSES' MLPRIM MVPRIM ;
  2188. che1 = che1
  2189. 'ET' ('CHAN' 'CHAM' ZMODL vite1 'STRESSES')
  2190. 'ET' ('CHAN' 'CHAM' ZMODL conti.'DEPLACEMENTS' 'STRESSES');
  2191. che2 = che2
  2192. 'ET' ('CHAN' 'CHAM' ZMODL ZDETOT 'STRESSES')
  2193. 'ET' ('CHAN' 'CHAM' ZMODL ZFCONST1 'STRESSES');
  2194. 'FINS';
  2195.  
  2196. 'SI' ither ;
  2197. TETDE=STAB12.'TET1' + (DTETD *('FLOT' (&sousinc - 1)/nsoincr));
  2198. TETDF=STAB12.'TET1' + (DTETD *('FLOT' &sousinc /nsoincr));
  2199. che3 = 'CHAN' 'CHAM' ZMODL TETDE 'STRESSES' ;
  2200. che4 = 'CHAN' 'CHAM' ZMODL TETDF 'STRESSES' ;
  2201. 'FINS' ;
  2202.  
  2203. che1 = che1 'ET' che3 ;
  2204. che2 = che2 'ET' che4 ;
  2205. aa5 = DEPST '*' ( (&sousinc '-' 1.) '/' nsoincr);
  2206. che5 = DEFT0 '+' aa5;
  2207. cc5 = DEPST '*' ('FLOT' &sousinc '/' nsoincr);
  2208. che6 = DEFT0 '+' CC5;
  2209. *On met les deformations en tete des champs pour COMP ('KTAN' 'PERT')
  2210. che1 = che5 'ET' che1 ;
  2211. che2 = che6 'ET' che2 ;
  2212. * pour les materiaux on garde toujours la valeur a la fin du pas
  2213. * et au debut du pas
  2214. * 'SI' ('NEG' ZMAT1 'INCONNU') ;
  2215. che1 = che1 'ET' ZMAT1 ;
  2216. * 'FINS' ;
  2217. *
  2218. * Modele NON LINEAIRE UTILISATEUR + GRANDES DEFORMATIONS
  2219. * On ajoute les gradients de deplacements qui seront transformes en
  2220. * gradients de transformation avant appel a UMAT (dans WKUMA1.ESO)
  2221. 'SI' ( 'NEG' GR_U_DEB 'INCONNU') ;
  2222. gru1 = GR_U_DEB + ( D_GR_U * ('FLOT' (&sousinc-1) / nsoincr));
  2223. che1 = che1 'ET' gru1 ;
  2224. gru2 = GR_U_DEB + ( D_GR_U * ('FLOT' &sousinc / nsoincr));
  2225. che2 = che2 'ET' gru2 ;
  2226. 'FINS';
  2227.  
  2228. 'SI' LNLOC ;
  2229. CHE1A = che1 ;
  2230. CHE2A = che2 ;
  2231.  
  2232. ISTEP = 1 ;
  2233. chm_z = 'ADET' 'NOUV' ZMODL CHASANST tsodeb 'STEP' ISTEP ;
  2234. che11 = CHE1A 'ET' chm_z ;
  2235. che11 = che11 'ET' ZSIG01 'ET' ZVAR01 'ET' ZEPS01 ;
  2236.  
  2237. chm_z = 'ADET' 'NOUV' ZMODL CHASANST tsofin 'STEP' ISTEP ;
  2238. che22 = CHE2A 'ET' chm_z ;
  2239. che22 = che22 'ET' ZMATT ;
  2240.  
  2241. 'SI' PARTLOCA ;
  2242. cha11 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che11 ;
  2243. cha22 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che22 ;
  2244. cho22 = 'ASSIS' 'TOUS' 'COMP' MODRELOC cha11 cha22 ;
  2245. cho2 = 'REDU' ZMODL ('ETG' cho22) ;
  2246. cha11 = 0 ; cha22 = 0 ; cho22 = 0 ;
  2247. 'SINON';
  2248. cho2 = 'COMP' ZMODL che11 che22 ;
  2249. 'FINS' ;
  2250.  
  2251. ZVARF = 'CHANGER' ('EXCO' cho2 com_var 'NOID')
  2252. 'TYPE' 'VARIABLES INTERNES' ;
  2253. * NLOC ne traitant pas des champs //, on reduit tout au modele initial.
  2254. 'SI' (PARALLEL 'ET' ('NON' PARTLOCA)) ;
  2255. ZVARF = 'REDU' ZVARF ZMODLI ;
  2256. 'FINS' ;
  2257. * cas non-local MOYE ou SB
  2258. 'SI' ('NEG' WTAB.'NON_LOCAL' 'HELM');
  2259. 'SI' ('EGA' WTAB.'NON_LOCAL' 'SB') ;
  2260. MOD_SB = WTAB.'NLOC_MODL' ;
  2261. CONTP = 'PRIN' ('REDU' ZSIG01 MOD_SB) MOD_SB ;
  2262. ZVARF = ZVARF '+' CONTP '+' WTAB.'NLOC_SB_REGU' ;
  2263. 'FINS' ;
  2264. ZVARN = 'NLOC' ZVARF WTAB.'CONN' ;
  2265. ZVARN = 'CHANGER' ('EXCO' ZVARN com_var)
  2266. 'TYPE' 'VARIABLES INTERNES' ;
  2267. 'SINON' ;
  2268. * cas non-local HELM
  2269. ERREUR ' Formulation non-locale indisponible' ;
  2270. 'FINSI' ;
  2271. * On rend ZVARN // si besoin :
  2272. 'SI' (PARALLEL 'ET' ('NON' PARTLOCA)) ;
  2273. ZVARN = 'REDU' ZVARN ZMODL ;
  2274. 'FINS' ;
  2275.  
  2276. ISTEP = 2 ;
  2277. chm_z = 'ADET' 'NOUV' ZMODL CHASANST tsodeb 'STEP' ISTEP ;
  2278. che11 = CHE1A 'ET' chm_z ;
  2279. che11 = che11 'ET' ZSIG01 'ET' ZVARN 'ET' ZEPS01 ;
  2280.  
  2281. chm_z = 'ADET' 'NOUV' ZMODL CHASANST tsofin 'STEP' ISTEP ;
  2282. che22 = CHE2A 'ET' chm_z ;
  2283. che22 = che22 'ET' ZMATT ;
  2284.  
  2285. 'SI' PARTLOCA ;
  2286. cha11 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che11 ;
  2287. cha22 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che22 ;
  2288. cho22 = 'ASSIS' 'TOUS' 'COMP' MODRELOC cha11 cha22 ;
  2289. cho2 = 'REDU' ZMODL ('ETG' cho22) ;
  2290. cha11 = 0 ; cha22 = 0 ; cho22 = 0 ;
  2291. 'SINON';
  2292. cho2 = 'COMP' ZMODL che11 che22 ;
  2293. 'FINS' ;
  2294.  
  2295. * Un peu de menage
  2296. CHE1A = 1 ;
  2297. CHE2A = 1 ;
  2298. ZVARN = 1 ;
  2299. ZVARF = 1 ;
  2300. chm_z = 1 ;
  2301.  
  2302. 'SINON' ;
  2303. che11 = che1 'ET' ZSIG01 'ET' ZVAR01 'ET' ZEPS01 ;
  2304. che22 = che2 'ET' ZMATT ;
  2305. * si ( 'EXIST' PRECED 'ECRIT' ) ;
  2306. * 'FINS';
  2307. 'SI' PARTLOCA;
  2308. cha11 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che11 ;
  2309. cha22 = 'ASSIS' 'TOUS' 'REDU' MODRELOC che22 ;
  2310. cho22 = 'ASSIS' 'TOUS' 'COMP' MODRELOC cha11 cha22 ;
  2311. cho2 = 'REDU' ZMODL ('ETG' cho22) ;
  2312. cha11 = 0 ; cha22 = 0 ; cho22 = 0 ;
  2313. 'SINON';
  2314. * 'SI' ( 'EXIS' PRECED 'ECRIT' );
  2315. * 'FINS';
  2316. cho2 = 'COMP' ZMODL che11 che22 ;
  2317. 'FINS';
  2318. 'FINS' ;
  2319. ZSIGF = 'CHANGER' ('EXCO' cho2 com_sig 'NOID')
  2320. 'TYPE' 'CONTRAINTES' ;
  2321. ZVARF = 'CHANGER' ('EXCO' cho2 com_var 'NOID')
  2322. 'TYPE' 'VARIABLES INTERNES' ;
  2323. ZDEIF = 'CHANGER' ('EXCO' cho2 com_dei 'NOID')
  2324. 'TYPE' 'DEFORMATIONS INELASTIQUES' ;
  2325.  
  2326. 'SI' ('EXIS' WTAB 'MOD_LIA') ;
  2327. ZFLIA = 'EXCO' cho2 'FLIA' 'NOID';
  2328. 'FINS' ;
  2329. 'SI' (nsoincr 'NEG' 1);
  2330. 'SI' (&sousinc 'EGA' nsoincr);
  2331. aaa1=zsigf/2;
  2332. aaa2=aaa1+ zsigm;
  2333. 'DETR' aaa1; 'DETR' zsigm;
  2334. zsigm=aaa2;
  2335. 'SINON';
  2336. aaa2=zsigf+ zsigm;
  2337. 'DETR' zsigm;
  2338. zsigm=aaa2;
  2339. 'FINS';
  2340. 'FINS';
  2341. zsig01 = zsigf ;
  2342. ZVAR01 = ZVARF ;
  2343. ZEPS01 = ZDEIF ;
  2344. ** defineTO=ZDEIF+defineTO;
  2345. 'FIN' sousinc;
  2346. *
  2347. * Matrice tangente par perturbation evaluee pour la derniere iteration calculee
  2348. * A voir : cas grand deplacement ZSIGF et ZMAT, cas poreux et thermique ZSIGF
  2349. 'SI' (IKTAN 'ET' IPERT) ;
  2350. Z1COMP = che11 ;
  2351. Z2COMP = che22 'ET' ZSIGF ;
  2352. 'FINS' ;
  2353. *
  2354. 'SI' (nsoincr 'NEG' 1);
  2355. * pour tenir compte de ce que le travail de la correction est 1/2 FU,
  2356. * on la multiplie par 2
  2357. CONT =(ZSIGF + ZSIG0) ;
  2358. ZSIGM = ZSIGM*(2. /nsoincr) ;
  2359. ZSIGM = ZSIGM - CONT ;
  2360. BZSIGM='BSIGMA' ZMODL ZSIGM ZMAT;
  2361. ** ZDEPSPL=defineTO - ZEPS0 ;
  2362. 'FINS';
  2363. * on enleve toutes traces de champs inutiles
  2364. zsig01=1;ZVAR01=1;ZEPS01=1;
  2365. ** defineTO=1; tabcont =1;
  2366.  
  2367. * le fin d'en dessous est le fin de "si isste ... sinon ... finsi"
  2368. * on est encore dans "si iplavi ... finsi"
  2369. 'FINS';
  2370. *...fin du si ISSTE sinon ...
  2371. *
  2372. * cas particulier poreux et thermique
  2373. *
  2374. 'SI' ITHER;
  2375. 'SI' POR1 ;
  2376. ZSIGF = ZSIGF - ( COEPI * DMSRT0 ) ;
  2377. 'FINS';
  2378. 'FINS';
  2379. *
  2380. * max de epse pendant l'iteration
  2381. * nbre de points qui ont une evolution non lineaire
  2382. *
  2383. 'SI' ('EGA' WTAB.'MOVA' 'RIEN') ;
  2384. ** EPSM = 0.; DPSMAX=0.;
  2385. MMC=0;
  2386. 'SINON';
  2387. XXX1 = 'EXCO' (WTAB.'MOVA') ZVARF ;
  2388. ** EPSM = 'MAXI' XXX1 ;
  2389. ACC = 'ABS' ( XXX1 - ACC0 ) ;
  2390. ** 'DETR' XXX1 ;
  2391. MMC = 'MASQUE' ACC 'SUPERIEUR' 1.D-10 'SOMME' ;
  2392. 'SI' (MMC > MMCMAX) ; MMCMAX = MMC ; 'FINS' ;
  2393. ** DPSMAX = 'MAXI' ACC ;
  2394. 'FINS';
  2395. *
  2396. * dans le cas SANS plasticite,viscoplas,endommagement
  2397. * on calcule directement le nouveau champ de contrainte
  2398. *
  2399. 'SINON' ;
  2400. ZSIGF = ZSIG0 + DSIGT ;
  2401. 'FINS';
  2402. * si il y a lieu transport sur configuration debut de pas
  2403. 'SI' IGRD ;
  2404. 'SI' LAG_TOT;
  2405. zsig0 = pica zsig0 zmodl zu1;
  2406. zsigf = pica zsigf zmodl zu1;
  2407. depst = pica depst zmodl zu1;
  2408. form geom1;
  2409. 'FINS';
  2410. 'FINS';
  2411. *
  2412. DPSMAXP = DPSMAX ;
  2413. DPSMAX = 'MAXI' 'ABS' DEPST 'AVEC' MLDEFOR ;
  2414. EPSM = 'MAXI' 'ABS' (ZDETOT 'EPSI' ZMAT ZMODL) 'AVEC' MLDEFOR ;
  2415. *
  2416. *--- CAS 4 -------------------------------------------------------------
  2417. * on determine les forces resultantes de l'increment------------------
  2418. *
  2419. * en grand deplacements on transforme pi en sigma
  2420. * on met a jour les coordonnees
  2421. *
  2422. ZMAT2=ZMAT;
  2423.  
  2424. 'SI' (IGRD 'ET' ('NON' HPP_EPS) );
  2425.  
  2426. coefmul = 1.d0 / (DPSMAX + XPETIT);
  2427. * coefmulr = -1;
  2428. ref = 75;
  2429. 'SI' ((&bzdept > 1) 'ET' ACC_Q 'ET' ACC_R );
  2430. ref = 200 ; 'FINS';
  2431. 'SI' ((coefmul < ref) 'ET' WTAB.'LINESEARCH');
  2432. * il faut reduire d'urgence les deplacements ;
  2433. zdeptlx = zdept 'EXCO' 'LX' 'LX';
  2434. ACC_Q = VRAI;
  2435. znacce=999; zmaxit = it;
  2436. * coefmulr = (2*dpsmax+1) ** 0.5 - 1 ;
  2437. * coefmulr = 5e1 / coefmulr;
  2438. * coefmuls = 8e-3 * coefmul;
  2439. * coefmul = coefmulr;
  2440. ** si (coefmuls < coefmulr); coefmul = coefmuls; finsi;
  2441. * si (coefmul > 0.3) ; coefmul = 0.3; finsi;
  2442. coefmul = 0.5;
  2443. 'SI' ((&bzdept 'EGA' 1) 'ET' (RED_URG < 5));
  2444. coefmul = coefmul * -1;
  2445. urg = vrai;
  2446. 'FINS';
  2447. zdeptq = zdeptp;
  2448. 'SI' (&bzdept > 1);zdeptq=zdeptini exco 'LX' 'NOID' 'LX'; finsi;
  2449. ZDEPT = zdeptq +( coefmul * (ZDEPT -zdeptq)) ;
  2450. * zdept = (zdept enle 'LX') + zdeptlx;
  2451.  
  2452. 'SI' (RED_URG > 2 ) ;
  2453. 'SI' ('NON' nonconv); 'MESS' ' non convergence detectee 1' ;
  2454. ZMAXIT = IT+5;
  2455. 'FINS';
  2456. nonconv = vrai;
  2457. 'FINS';
  2458. 'SI' (&bzdept 'EGA' 1);
  2459. 'MESS' 'reduction d urgence des deplacements';
  2460. hpp_eps=faux; RED_URG = RED_URG + 1;
  2461. 'FINS';
  2462. SI (RED_URG > 6 ); ZMAXIT = IT; 'FINS';
  2463.  
  2464. 'ITERER' bzdept;
  2465. 'FINS';
  2466. 'FINS';
  2467.  
  2468. 'FINS';
  2469.  
  2470. * au dessus le finsi du si fefp sinon(Update or total lagrangian) ----
  2471. *
  2472. * linesearch -----------------------------------------------------------
  2473. *
  2474. 'SI' ('NON' WTAB.'LINESEARCH'); 'QUITTER' bzdept; 'FINS';
  2475. 'SI' (it < 2); 'QUITTER' bzdept; 'FINS';
  2476. 'SI' ('MULT' IT ZNACCE); 'QUITTER' bzdept; 'FINS';
  2477. 'SI' (HPP_EPS); 'QUITTER' bzdept; 'FINS';
  2478. 'SI' ACC_Q ; 'QUITTER' bzdept; 'FINS';
  2479. ACC_Q = VRAI;
  2480. ACC_R = FAUX;
  2481.  
  2482. r1 = residc '*' -1. 'ENLE' 'FLX';
  2483.  
  2484. zdeptdif = (zdept '-' zdeptp) ;
  2485. krr = (zclim '*' zdept) '-' (zclimp '*' zdeptp) 'ENLE' 'FLX';
  2486. GEOTEMP='FORM' zdeptp;
  2487. k1 = bsigma (sigma zdeptdif zmodl zmat) zmodl zmat;
  2488. k1 = k1 '-' krr ;
  2489. 'FORM' geom1 ;
  2490. 'DETR' GEOTEMP;
  2491.  
  2492. GEOTEMP='FORM' zdept ;
  2493. k2 = bsigma (sigma zdeptdif zmodl zmat) zmodl zmat;
  2494. k2 = k2 '-' krr;
  2495.  
  2496. 'FORM' geom1 ;
  2497. 'DETR' GEOTEMP;
  2498.  
  2499. alpha = k2 '-' k1 ;
  2500. beta = k1 ;
  2501.  
  2502. * nouvelle acceleration de convergence pour les grands deplacements
  2503. a0 = xty beta r1 mldual mldual ;
  2504. a1 =(xty beta beta mldual mldual) '+'(xty alpha r1 mldual mldual);
  2505. a2 =(3./2.) * (xty alpha beta mldual mldual);
  2506. a3 =(1./2.) * (xty alpha alpha mldual mldual);
  2507.  
  2508. xr1 = racpol a0 a1 a2 a3;
  2509. lambd = extr xr1 1;
  2510.  
  2511. 'SI' WTAB.'STABILITE' ;
  2512. 'SI' (lambd < 0);lambd = 0.2;'FINS';
  2513. 'FINS';
  2514.  
  2515. 'SI' ((dime xr1) > 1 );
  2516. lambd2 = extr xr1 2;
  2517. lambd3 = extr xr1 3;
  2518.  
  2519. 'SI' WTAB.'STABILITE' ;
  2520. 'SI' (lambd2 < 0); lambd2 = 0.2; 'FINS';
  2521. 'SI' (lambd3 < 0); lambd3 = 0.2; 'FINS';
  2522. 'FINS';
  2523.  
  2524. 'SI' ((lambd2 '-' 1.5 'ABS') '<' (lambd '-' 1.5 'ABS'));
  2525. lambd = lambd2;
  2526. 'FINS';
  2527.  
  2528. 'SI' ((lambd3 '-' 1.5 'ABS') '<' (lambd '-' 1.5 'ABS'));
  2529. lambd = lambd3;
  2530. 'FINS';
  2531. 'FINS';
  2532.  
  2533. lambdav=lambd;
  2534. 'SI' (lambd > 4 ); lambd = 4 ; 'FINS';
  2535. 'SI' (('ABS' lambd < 1e-4) 'ET' (itacc 'NEG' 3) 'ET'
  2536. (dpsmax > zprecnc));
  2537. lambd = 0.1 ; urg = vrai; znacce = 999; itacc=4;
  2538. 'FINS';
  2539. 'SI' ((lambd < 0) 'ET' WTAB.'STABILITE'); lambd = 0.3 ; 'FINS';
  2540. 'SI' ('ABS' lambd < 1e-3 ); lambd = 1e-3; 'FINS';
  2541. 'SI' ( lambd < -3 ); lambd =-3. ; 'FINS';
  2542. 'SI' ((lambd < 0) 'ET' (itacc 'NEG' 3)) ; znacce = 999;'FINS';
  2543. *
  2544. * verif sens
  2545. *
  2546. **'SI' WTAB.'STABILITE' ;
  2547. ** 'SI' (lambd < 0); lambd = 0.3; finsi;
  2548. ** znacce = 999;
  2549. **'FINS';
  2550. * zdeptm = (zdep1 'ENLE' 'LX') * (lambd - 1);
  2551. zdeptmm = zdep1 'EXTRAI' 'MAIL';
  2552. zdeptm = (zdep1 - ((zdeptp 'EXCO' 'LX' 'LX') 'REDU' zdeptmm))
  2553. * (lambd - 1);
  2554. zdept = zdept + zdeptm ;
  2555.  
  2556. * pour l(autre acceleration de convergence
  2557. correc = correc - (ZRAID * zdeptm);
  2558.  
  2559. fin bzdept;
  2560. * fin de boucle de reduction de zdept
  2561. *---------------------------------------------------------------------
  2562.  
  2563. 'SI' ('NON' IFEFP);
  2564. 'SI' (IGRD 'ET' ('NON' HPP_EPS));
  2565. * grands deplacements -----------------------------------------------
  2566.  
  2567. ZPI = ZSIGF ;
  2568. zdepf = zdept ;
  2569. ZSIGF = 'PICA' ZMODL ZPI ZDEPF ;
  2570. 'DETR' ZPI ;
  2571.  
  2572. 'SI' ITCAR ;
  2573. GEOM2 ZMAT2 ='FORM' ZDEPF ZMODLI ZMATI;
  2574. ZMAT2 ='REDU' ZMAT2 ZMODL ;
  2575. 'SINON';
  2576. GEOM2 = 'FORM' ZDEPF ;
  2577. 'FINS';
  2578. FEQU2 = 'BSIGMA' ZMODL ZSIGF ZMAT2 ;
  2579. 'SI' IRCON;
  2580. FEQU2 = FEQU2 '+'('REDU' MAI_CONS (RIG_CONS '*' ZDETOT));
  2581. 'FINS';
  2582. 'SI' LOGPRE ;
  2583. ZDFINI = ZDFINI - ZFPEXTF ;
  2584. 'DETR' ZFPEXTF ;
  2585. ZPEXTF = 'TIRE' WTAB.'CHARGEMENT' 'PRES' TI ;
  2586. 'SI' ('EXIS' WTAB 'MAT_PRE') ;
  2587. ZFPEXTF = 'BSIG' MOP ZPEXTF ('REDU' ZMATFI MOP) ;
  2588. 'SINON' ;
  2589. ZFPEXTF = 'BSIG' MOP ZPEXTF ;
  2590. 'FINS' ;
  2591. ZDFINI = ZDFINI + ZFPEXTF ;
  2592. 'FINS' ;
  2593.  
  2594. 'SI' ADDISEC2;
  2595. ZDFINI = ZDFINI '-' FP22 ;
  2596. TFP22 = CHARMECA PRECED WTAB.'T_FINAL' ;
  2597. ADDISEC2 = FAUX;
  2598. 'SI' ('EXIS' TFP22 'ADDI_SECOND') ;
  2599. FP22 = TFP22.'ADDI_SECOND' ;
  2600. ZDFINI = ZDFINI + FP22 ;
  2601. ADDISEC2 = VRAI ;
  2602. 'FINS';
  2603. 'SI' IPILOT ;
  2604. COEFP =(COEPI*DMZPRES+ MZPRES0)/(MZPRES+XPETIT);
  2605. 'FINS' ;
  2606. 'FINS' ;
  2607. 'FORM' GEOM1 ;
  2608. *
  2609. * dans les autres que IGRD cas on calcule fequ2(desequilibre)---------
  2610. *
  2611. 'SINON';
  2612. FEQU2 = 'BSIGMA' ZMODL ZSIGF ZMAT ;
  2613. 'SI' IRCON;
  2614. FEQU2 = FEQU2 '+' ('REDU' MAI_CONS (RIG_CONS '*' ZDETOT));
  2615. 'FINS';
  2616. 'SI' ADDISEC2;
  2617. ZDFINI = ZDFINI '-' FP22 ;
  2618. TFP22 = CHARMECA PRECED WTAB.'T_FINAL' ;
  2619. ADDISEC2 = FAUX;
  2620. 'SI' ('EXIS' TFP22 'ADDI_SECOND');
  2621. FP22 = TFP22.'ADDI_SECOND';
  2622. ZDFINI = ZDFINI '+' FP22 ;
  2623. ADDISEC2 = VRAI ;
  2624. 'FINS';
  2625. 'SI' IPILOT ;
  2626. COEFP =(COEPI*DMZPRES+ MZPRES0)/(MZPRES+XPETIT);
  2627. 'FINS' ;
  2628. 'FINS' ;
  2629. 'FINS';
  2630. 'FINS';
  2631. * le finsi au dessus est le finsi de si iefefp
  2632. 'SI' ('EXIS' WTAB 'MOD_LIA') ;
  2633. FLIAI = 0.d0 ;
  2634. 'SI' ('EGA' ('TYPE' ZFLIA) 'MCHAML ') ;
  2635. NZLIA = 'EXTR' ZFLIA 'NBZO' ;
  2636. 'SI' (NZLIA > 0) ;
  2637. 'REPETER' BZLIA NZLIA ;
  2638. * un point support par zone - 2010 kich
  2639. FLIAI= FLIAI +
  2640. ('EXCO' ('EXTR' ZFLIA 'FLIA' &BZLIA 1 1) MLDUAL 'NOID') ;
  2641. 'FIN' BZLIA ;
  2642. 'FINS' ;
  2643. FEQU2 = FEQU2 + FLIAI ;
  2644. 'FINS' ;
  2645. 'FINS' ;
  2646. *
  2647. * calcul des forces correctrices de frottement-------------------------
  2648. *
  2649. 'SI' (WTAB.'CAFROTTE');
  2650. * apres non convergence, on maintient les forces de frottement
  2651. MAINT=FAUX;
  2652. 'SI' ('NON' WTAB.'CONV' 'ET' ('NEG' WTAB.'FFROT' 'INCONNU'));
  2653. MAINT=VRAI;
  2654. 'FINS';
  2655. zsigfT = 'REDU' zsigf ZMODLI;
  2656. MODFRO = 'EXTR' WTAB.'CONT_REEL' 'COMP' 'FROTTANT';
  2657. MATFRO = 'REDU' WTAB.'MAT_FRO' WTAB.'CONT_REEL';
  2658. 'SI' MAINT;
  2659. FFROT1 = 'EXCF' ZRAID ZDEPT MODFRO MATFRO WTAB.'FFROT' ZSIGFt ;
  2660. 'SINON';
  2661. FFROT1 = 'EXCF' ZRAID ZDEPT MODFRO MATFRO ZSIGFt ;
  2662. 'FINS';
  2663. FFROT = ((1*FFROT1) + FFROTP)/2;
  2664. 'DETR' FFROTP;
  2665. FFROTP = FFROT1;
  2666. 'FINS';
  2667.  
  2668.  
  2669. 'SI' ISOL ;
  2670. 'SI' IPLAVI ;
  2671. ZMATPF = 'CHANGER' ('EXCO' cho2 com_maa 'NOID')
  2672. 'TYPE' 'CARACTERISTIQUES';
  2673. MA_POR = 'REDU' ZMATPF MO_POR ;
  2674. 'FINS';
  2675. XXX1 = 'GRAD' MO_POR ZDETOT MA_POR 'CONS' ;
  2676. XXXS = (1. - WTAB.'TETA' )*GRAP0 + (WTAB.'TETA'*XXX1) ;
  2677. XXX2 = STAB12.'DT' * ('GNFL' MO_POR XXXS) ;
  2678. XXX3 = FEQU2 ;
  2679. FEQU2 = XXX3 - XXX2;'DETR' XXX3; 'DETR' XXX2 ;
  2680. * 'DETR' XXX1; 'DETR' XXXS ;
  2681. 'FINS' ;
  2682.  
  2683. * si on a fait de la sous incrementation ( a cause du dynamique
  2684. 'SI' ( nsoincr > 1);
  2685. FEQU2 = FEQU2 + BZSIGM ;
  2686. 'FINS';
  2687.  
  2688. *
  2689. * forces qui varient en dynamique------------------------------------
  2690. * b*sigma + 4/h2*m*dx
  2691. * + 2/h *c*dx
  2692. *
  2693. 'SI' IDYN;
  2694. 'DETR' FFDYN;
  2695. FFDYN = 'COPIER' FEQU2 ;
  2696. XXX1 = WTAB.'MASSE' * ZDEPT;
  2697. XXX3 = 4. * UNSURH * UNSURH * XXX1 ;
  2698. 'DETR' XXX1;
  2699. 'SI' ('EXIS' WTAB 'MOD_LIA');
  2700. 'SI' (NZLIA > 0) ;
  2701. XVIT2 = 0.d0;
  2702. 'REPETER' BZLIA NZLIA ;
  2703. XVIT2 = XVIT2 +
  2704. ('EXCO' ('EXTR' ZFLIA 'FLIA' &BZLIA 1 1) MVPRIM MLPRIM 'NOID') ;
  2705. 'FIN' BZLIA ;
  2706. 'FINS' ;
  2707. 'FINS' ;
  2708. 'SI' ( 'NEG' WTAB.'AMORTISSEMENT' 'INCONNU') ;
  2709. XXX1 = WTAB.'AMORTISSEMENT' * ZDEPT;
  2710. XXX2 = 2. * UNSURH * XXX1 ;
  2711. XXX4 = XXX3 + XXX2 ;
  2712. 'DETR' XXX3 ;
  2713. 'DETR' XXX1 ;
  2714. 'DETR' XXX2 ;
  2715. XXX3 = XXX4 ;
  2716. 'FINS' ;
  2717. XXX4 = FEQU2 + XXX3;
  2718. 'DETR' XXX3 ;
  2719. 'DETR' FEQU2 ;
  2720. FEQU2 = XXX4 ;
  2721. *
  2722. * forces correctrices en cas de liaison persistante :
  2723. * on veut avoir (forces inertielles + forces visqueuses) compatibles
  2724. * avec accelerations et vitesses relatives nulles aux points de contact
  2725. * (pendant le contact). On modifie fequ2 ---> residu et l'iteration sui
  2726. * fournira les bonnes reactions
  2727. *
  2728. 'SI' IMPLP;
  2729. VADD XXX3 = VITETFOR ZRAID_T WTAB FCORU ;
  2730. XXX4 = FEQU2 - XXX3;
  2731. 'DETR' XXX3 ;
  2732. 'DETR' FEQU2 ;
  2733. FEQU2 = XXX4;
  2734. 'FINS' ;
  2735. 'FINS';
  2736. *
  2737. * preparation des tests de convergence.-------------------------------
  2738. * on calcule la plus grande composante de f + freac(sans flx)
  2739. * fequi est de plus le second membre que l'on veut equilibrer
  2740. * ne pas oublier d'y ajouter les forces suiveuses
  2741. * test de convergence: equilibre de la structure
  2742. * il se fait sur residc, on calcule auusi le nouveau second
  2743. * membre qui est residc moins les forces d'appuis, il faut aussi
  2744. * remettre a jour les FLX qui travaillent eux en incremental
  2745. * pour les depl imposes cela revient a imposer un increment nul et
  2746. * pour les rela unilaterales cela revient a mettre a jour le jeu
  2747. *
  2748. RESIDUP= RESIDU ;
  2749. XXX1 = ZFEXT ;
  2750. 'SI' ADDISEC2 ;
  2751. XXX2= XXX1 '+' (FP22 '*' COEFP) ;
  2752. XXX1= XXX2 ;
  2753. 'FINS';
  2754.  
  2755. 'SI' (LOGPRE 'ET' IGRD) ;
  2756. XXX2 = XXX1 '+' (ZFPEXTF '*' COEFP) ;
  2757. XXX1 = XXX2 ;
  2758. 'FINS' ;
  2759.  
  2760. * forces exterieures(+ autres termes p.ex. en dynamique ou en poreux)
  2761. * sans reactions
  2762. * - forces interieures(+ autres termes p.ex. en dynamique ou en poreux)
  2763.  
  2764. RESIDU = XXX1 - FEQU2;
  2765. * forces exterieures(+ autres termes p.ex. en dynamique ou en poreux)
  2766. * + reactions
  2767. FEQUI = XXX1 - FCORF;
  2768.  
  2769. * forces exterieures(+ autres termes p.ex. en dynamique ou en poreux)
  2770. * + reactions - forces interieures(+ autres termes p.ex. en dynamique
  2771. * ou en poreux)
  2772. RESIDC= RESIDU - FCORF;
  2773.  
  2774. * relations unilaterales et autres(Reinitialisation)
  2775. *'SI' IMPO12 ;
  2776. XXX1 = RESIDU +ZFLX1;'DETR' RESIDU;
  2777. RESIDU = XXX1 - FCORU;'DETR' XXX1;
  2778. 'SI' (IPILOT 'ET' IMPO12);
  2779. XXX3 =COEPI-1.D0*DIMPOV;
  2780. XXX1= RESIDU + XXX3;'DETR' RESIDU;RESIDU=XXX1;
  2781. 'DETR' XXX3;
  2782. 'FINS';
  2783. *'FINS';
  2784. XXX1 = FEQUI - ZFPLO;
  2785. ZFAU1 = 'ABS' XXX1 ;'DETR' XXX1;
  2786. XAUXF = 'MAXI' ZFAU1 'ABS' 'SANS' MXMFLX;
  2787. 'SI' TSTMOM ;
  2788. XAUXM = 'MAXI' ZFAU1 'ABS' 'AVEC' MXMYMZ;
  2789. 'FINS';
  2790. 'DETR' ZFAU1;
  2791. **'DETR' ZDETOT ;
  2792. COEINC = 0.D0;
  2793. *
  2794. * Calcul du nouveau alpha pour le pilotage -----------------------------
  2795. *
  2796. 'SI' (IPILOT 'ET' PASTEST) ;
  2797. * seulement si on a accelere
  2798. *
  2799. 'SI' (WTAB.'VISCOPLASTIQUE' 'OU' WTAB.'VISCODOMMAGE'
  2800. 'OU' ITHER 'OU' LOGDEF );
  2801. 'MESS' 'ALPHA calcule avec la norme de l increment' ;
  2802. COEPI0 = COEPI;
  2803. COEPI = 'MINI' ('PROG' 1. (STAB12.'AUTOCOEF' * al1) );
  2804. * les acceleration auront lieu toutes les 4 pas
  2805. COEINC =COEPI0 - COEPI ;
  2806. * ITACC=3 ;
  2807. 'SINON' ;
  2808. * les multiplicateurs sont affectes d'un signe -
  2809. XXX1 = 'EXCO' ZDEPT 'LX' 'NOID' 'LX';
  2810. XXX2=2.D0 * XXX1;'DETR' XXX1;
  2811. XXX4=ZDEPT - XXX2;'DETR' XXX2;
  2812. XXX3=DFEXT0L *( COEPI * ( 1 - AL1 ));
  2813. XXX2 = XXX3 + RESIDC;'DETR' XXX3;XXX5=RESIDNOR* 1.D0;
  2814. 'SI' IMPO12 ;
  2815. XXX1= COEPI * ( 1 - AL1 ) * DIMPO12;
  2816. XXX3= XXX2 - XXX1;'DETR' XXX1;'DETR' XXX2;
  2817. XXX2 = XXX3;XXX1 = XXX5 - DIMPO12;
  2818. 'DETR' XXX5; XXX5=XXX1;
  2819. 'FINS';
  2820. 'SI' ADDISEC2;
  2821. XXX7 = XXX5 '+' FP22 ;
  2822. CORFR = XXX7 '-' FP022;
  2823. 'FINS';
  2824.  
  2825. XX3 ='XTY' RESIDNOR XXX4 MLDUAL MLPRIM;
  2826. XX1 ='XTY' XXX2 XXX4 MLDUAL MLPRIM;
  2827. XX2 ='XTY' XXX5 XXX4 MLDUAL MLPRIM;
  2828. COEINC = XX1 '/' XX2;
  2829. 'DETR' XXX4;'DETR' XXX2;'DETR' XXX5;
  2830. COEPI0 = COEPI;
  2831. COEPI = 'MINI' ('PROG' (COEPI0 - COEINC) 1.d0 ) ;
  2832. 'SI' ((COEPI0 'EGA' 1d0) 'ET' (AL1 'EGA' 1D0));
  2833. COEPI = 1d0;
  2834. 'FINS';
  2835. COEINC = COEPI0 - COEPI;
  2836. 'SI' ((COEPI 'EGA' 1d0) 'ET' (COEPI0 'NEG' 1d0));ITACC=3;
  2837. 'FINS';
  2838. 'FINS' ;
  2839. * mise a jour de ZFEXT en fonction de coepi on ne met a jour
  2840. * de maniere definitive que les forces, les modifs sur les FLX
  2841. * pour reimposer les conditions de deplacements ne
  2842. * sont faites que pour le premier residu attention au cas des jeux
  2843. * sortis de IMPO qu'il ne faut pas toucher
  2844. XXX3 = DFEXT0F '*' COEINC ;
  2845. XXX2 = ZFEXT '-' XXX3;'DETR' XXX3;
  2846. 'DETR' ZFEXT; ;ZFEXT = XXX2;
  2847. 'SI' ADDISEC2;
  2848. MZPRES2=COEPI '*' DMZPRES '+' MZPRES0;
  2849. COEFP =MZPRES2 '/'(MZPRES '+' XPETIT);
  2850. XXX3 = FP22 '*' COEFP;
  2851. XXX2 = ZFEXT '+' XXX3;
  2852. 'DETR' XXX3;
  2853. 'FINS';
  2854. XXX1= XXX2 - FEQU2;
  2855. 'DETR' RESIDU;XXX3=COEPI * DFEXT0L; XXX4= XXX3 + FLXINI;
  2856. 'DETR' XXX3; XXX3 = XXX4 - FCORU;'DETR' XXX4;
  2857. RESIDU = XXX1 + XXX3;'DETR' XXX3;
  2858. 'SI' IMPO12;
  2859. XXX3= 1.D0 - COEPI * DIMPO12;
  2860. XXX2= RESIDU + XXX3; 'DETR' RESIDU;
  2861. RESIDU = XXX2; 'DETR' XXX3;
  2862. 'FINS';
  2863. 'SI' (ITHER 'OU' WTAB.'MATVAR') ;
  2864. XXX3 = FTHE * COEINC;
  2865. XXX2= RESIDU + XXX3; 'DETR' RESIDU;
  2866. RESIDU = XXX2; 'DETR' XXX3;
  2867. 'FINS';
  2868. <