Télécharger creep03_cisXY.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : creep03_cisXY.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *=======================================================================
  5. *
  6. * CAS TEST DE VALIDATION DE LA PRISE EN COMPTE D'UNE LOI DE
  7. *
  8. * COMPORTEMENT MECANIQUE NON LINEAIRE EXTERNE
  9. *
  10. *=======================================================================
  11. *
  12. * Modele 'VISCO_EXTERNE' 'GENERAL', integre par CCREEP (schema CCONST)
  13. *
  14. * Test No3 : Modele de FLUAGE POLYNOMIAL
  15. * CISAILLEMENT XY
  16. * Modele externalise non instanciable (ILOI=21) compare
  17. * au modele externalise non instanciable avec variables
  18. * internes supplementaires (ILOI=23)
  19. *
  20. *=======================================================================
  21.  
  22.  
  23. opti echo 0 ;
  24.  
  25.  
  26. *=======================================================================
  27. *
  28. * PROCEDURES DEPENDANTES
  29. *
  30. *=======================================================================
  31.  
  32.  
  33.  
  34.  
  35. *=======================================================================
  36. *
  37. * JEU DE DONNEES PRINCIPAL
  38. *
  39. *=======================================================================
  40.  
  41. *-----------------------------------------------------------------------
  42. * 0 - Maillage
  43. * PLAQUE CARREE COTE : 1 m EPAISSEUR : 1 cm
  44.  
  45. opti dime 3 elem cub8 mode tridim ;
  46.  
  47. l1 = 1.0 ;
  48. l2 = 1.0 ;
  49. ep0 = 0.01 ;
  50.  
  51. p01 = 0.0 0.0 0.0 ;
  52. p02 = ep0 0.0 0.0 ;
  53. p03 = ep0 l2 0.0 ;
  54. p04 = 0.0 l2 0.0 ;
  55.  
  56. p10 = 0.0 0.0 l1 ;
  57. p20 = ep0 0.0 l1 ;
  58. p30 = ep0 l2 l1 ;
  59. p40 = 0.0 l2 l1 ;
  60.  
  61. li01 = 'DROIT' 1 p01 p02 ;
  62. li02 = 'DROIT' 1 p02 p03 ;
  63. li03 = 'DROIT' 1 p03 p04 ;
  64. li04 = 'DROIT' 1 p04 p01 ;
  65.  
  66. li10 = 'DROIT' 1 p10 p20 ;
  67. li20 = 'DROIT' 1 p20 p30 ;
  68. li30 = 'DROIT' 1 p30 p40 ;
  69. li40 = 'DROIT' 1 p40 p10 ;
  70.  
  71. sur01 = 'DALL' li01 li02 li03 li04 ;
  72. sur10 = 'DALL' li10 li20 li30 li40 ;
  73. mail1 = sur01 'VOLU' 1 sur10 ;
  74.  
  75. li2 = 'DROIT' 1 p02 p20 ;
  76. li3 = 'DROIT' 1 p20 p10 ;
  77. li4 = 'DROIT' 1 p10 p01 ;
  78. sbloq = 'DALL' li01 li2 li3 li4 ;
  79.  
  80. li5 = 'DROIT' 1 p04 p03 ;
  81. li6 = 'DROIT' 1 p03 p30 ;
  82. li7 = 'DROIT' 1 p30 p40 ;
  83. li8 = 'DROIT' 1 p40 p04 ;
  84. schar = 'DALL' li5 li6 li7 li8 ;
  85.  
  86. ELIM (mail1 ET sbloq ET schar) 1.e-5 ;
  87.  
  88. *-----------------------------------------------------------------------
  89. * 1 - Calcul avec modele 'VISCO_EXTERNE' 'GENERAL' no21
  90. * Modele 'FLUAGE' 'POLYNOMIAL' externalise dans CREEP
  91. * Modele non instanciable au niveau des iterations internes
  92.  
  93. * 1.1 Definition du modele utilisateur externalise dans CREEP
  94.  
  95. mo_uti1 = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  96. 'VISCO_EXTERNE' 'GENERAL' 'NUME_LOI' 21 ;
  97.  
  98. * 1.2 Affectation des proprietes de materiau
  99.  
  100. densiu = 0.95 ;
  101. dgrain = 8.0E-6 ;
  102. tuo2 = 1500.0 ;
  103. tauxfisk = 1.5E19 ;
  104.  
  105. xyoun = (2.2693E11 - (1.5399E7*tuo2)) - (9.35973E3 * (tuo2**2.0)) ;
  106. xcis = (8.583E10 - (5.157E6*tuo2)) - (3.747E3 * (tuo2**2.0)) ;
  107. xnu = (xyoun/(2.0*xcis)) - 1.0 ;
  108. xrho = 10950.0 ;
  109. xalph = ( (1.1833E-5 - (5.013E-9*tuo2)) + (3.756E-12 * (tuo2**2.0)) )
  110. - (6.125E-17 * (tuo2**3.0)) ;
  111.  
  112. ma_uti1 = MATE mo_uti1 'YOUN' xyoun 'NU ' xnu 'RHO ' xrho
  113. 'ALPH' xalph ;
  114.  
  115. * 1.3 Definition des conditions aux limites et chargements
  116.  
  117. cl1 = BLOQ mail1 'UY' 'UZ' ;
  118. cl2 = BLOQ sbloq 'UX' ;
  119. cl0 = cl1 ET cl2 ;
  120.  
  121. depmax = 5.0E8 ;
  122. dep0 = FORCE schar 'FX' depmax 'FY' 0.0 'FZ' 0.0 ;
  123. ev0 = EVOL MANU temps (PROG 0.0 1.E-2 1000.0)
  124. y (PROG 0.0 1.0 1.0 ) ;
  125. cha0 = CHAR MECA dep0 ev0 ;
  126.  
  127. * 1.4 Definition des instants de calcul
  128.  
  129. pro0 = PROG 0.0 PAS 1.E-6 3.0E-5 ;
  130. pro1 = PROG 0.0 PAS 1.E-6 3.0E-5 ;
  131.  
  132. * 1.5 Resolution par PASAPAS
  133.  
  134. tabuti1 = TABLE ;
  135. tabuti1.'MODELE' = mo_uti1 ;
  136. tabuti1.'CARACTERISTIQUES' = ma_uti1 ;
  137. tabuti1.'CHARGEMENT' = cha0 ;
  138. tabuti1.'BLOCAGES_MECANIQUES' = cl0 ;
  139. tabuti1.'TEMPS_CALCULES' = pro0 ;
  140. tabuti1.'TEMPS_SAUVES' = pro1 ;
  141. TMASAU=table;
  142. tabuti1 . 'MES_SAUVEGARDES'=TMASAU;
  143. TMASAU .'DEFTO'=VRAI;
  144. TMASAU .'DEFIN'=VRAI;
  145. TEMPS 'ZERO' ;
  146. PASAPAS tabuti1 ;
  147. TABTPS = TEMP 'NOEC';
  148. CPUuti1 = TABTPS.'TEMPS_CPU'.'INITIAL' ;
  149.  
  150. *-----------------------------------------------------------------------
  151. * 2 - Calcul avec modele 'VISCO_EXTERNE' 'GENERAL' no23
  152. * Modele 'FLUAGE' 'POLYNOMIAL' externalise dans CREEP
  153. * Modele non instanciable au niveau des iterations internes,
  154. * comprenant deux variables internes supplementaires
  155.  
  156. * 2.1 Definition du modele utilisateur externalise dans CREEP
  157.  
  158. LCVAR23 = MOTS 'PSUP' 'QSUP' ;
  159.  
  160. mo_uti2 = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  161. 'VISCO_EXTERNE' 'GENERAL'
  162. 'NUME_LOI' 23 'C_VARINTER' LCVAR23 ;
  163.  
  164. * 2.2 Affectation des proprietes de materiau
  165.  
  166. ma_uti2 = MATE mo_uti2 'YOUN' xyoun 'NU ' xnu 'RHO ' xrho
  167. 'ALPH' xalph ;
  168.  
  169. * 2.3 Definition des conditions aux limites et chargements : idem 1.3
  170.  
  171. * 2.4 Definition des instants de calcul : idem 1.4
  172.  
  173. * 2.5 Resolution par PASAPAS
  174.  
  175. tabuti2 = TABLE ;
  176. tabuti2.'MODELE' = mo_uti2 ;
  177. tabuti2.'CARACTERISTIQUES' = ma_uti2 ;
  178. tabuti2.'CHARGEMENT' = cha0 ;
  179. tabuti2.'BLOCAGES_MECANIQUES' = cl0 ;
  180. tabuti2.'TEMPS_CALCULES' = pro0 ;
  181. tabuti2.'TEMPS_SAUVES' = pro1 ;
  182. TMASAU=table;
  183. tabuti2 . 'MES_SAUVEGARDES'=TMASAU;
  184. TMASAU .'DEFTO'=VRAI;
  185. TMASAU .'DEFIN'=VRAI;
  186. TEMPS 'ZERO' ;
  187. PASAPAS tabuti2 ;
  188. TABTPS = TEMP 'NOEC';
  189. CPUuti2 = TABTPS.'TEMPS_CPU'.'INITIAL' ;
  190.  
  191. *-----------------------------------------------------------------------
  192. * 3 - Post-traitement : comparaisons entre les deux modeles
  193. *
  194. * N.B. Comme on simule un essai de cisaillement XY :
  195. * - Critere strict (errrel < 1.d-8) :
  196. * * sur GAXY, GIXY et SMXY,
  197. * * sur les variables internes EC0, ESW0, P, QTILD,
  198. * * pour le modele externe no23 : sur les variables internes
  199. * supplementaires PSUP et QSUP (comparees a P et QTILD).
  200. * - Pour les autres composantes on note les depassements du
  201. * critere si impression demandee
  202.  
  203. uxu1 = TABLE ; uyu1 = TABLE ; uzu1 = TABLE ;
  204. epxxu1 = TABLE ; epyyu1 = TABLE ; epzzu1 = TABLE ;
  205. gaxyu1 = TABLE ; gaxzu1 = TABLE ; gayzu1 = TABLE ;
  206. eixxu1 = TABLE ; eiyyu1 = TABLE ; eizzu1 = TABLE ;
  207. gixyu1 = TABLE ; gixzu1 = TABLE ; giyzu1 = TABLE ;
  208. smxxu1 = TABLE ; smyyu1 = TABLE ; smzzu1 = TABLE ;
  209. smxyu1 = TABLE ; smxzu1 = TABLE ; smyzu1 = TABLE ;
  210. ec0u1 = TABLE ; esw0u1 = TABLE ; pu1 = TABLE ; qtldu1 = TABLE ;
  211.  
  212. uxu2 = TABLE ; uyu2 = TABLE ; uzu2 = TABLE ;
  213. epxxu2 = TABLE ; epyyu2 = TABLE ; epzzu2 = TABLE ;
  214. gaxyu2 = TABLE ; gaxzu2 = TABLE ; gayzu2 = TABLE ;
  215. eixxu2 = TABLE ; eiyyu2 = TABLE ; eizzu2 = TABLE ;
  216. gixyu2 = TABLE ; gixzu2 = TABLE ; giyzu2 = TABLE ;
  217. smxxu2 = TABLE ; smyyu2 = TABLE ; smzzu2 = TABLE ;
  218. smxyu2 = TABLE ; smxzu2 = TABLE ; smyzu2 = TABLE ;
  219. ec0u2 = TABLE ; esw0u2 = TABLE ; pu2 = TABLE ; qtldu2 = TABLE ;
  220. psupu2 = TABLE ; qsupu2 = TABLE ;
  221.  
  222. LDEPL = MOTS 'UX ' 'UY ' 'UZ ' ;
  223. LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  224. LDEIN = MOTS 'EIXX' 'EIYY' 'EIZZ' 'GIXY' 'GIXZ' 'GIYZ' ;
  225. LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  226. LVARI = MOTS 'EC0 ' 'ESW0' 'P ' 'QTLD' ;
  227. LVSUP = MOTS 'PSUP' 'QSUP' ;
  228.  
  229. idefo = 4 ;
  230. idein = 4 ;
  231. icont = 4 ;
  232. ivar1 = 1 ;
  233. ivar2 = 2 ;
  234. ivar3 = 3 ;
  235. ivar4 = 4 ;
  236. ivas1 = 1 ;
  237. ivas2 = 2 ;
  238.  
  239. ierr = 0 ;
  240. critere = 1.d-8 ;
  241. IMESS = 0 ;
  242. ***IMESS = 1 ;
  243.  
  244. nbno = 8 ;
  245. nbgau = 8 ;
  246.  
  247. netat = DIME (tabuti1.'TEMPS') ;
  248. ietat = -1 ;
  249.  
  250.  
  251. * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  252.  
  253. REPETER BECART netat ;
  254.  
  255. ietat = ietat + 1 ;
  256.  
  257. temps1 = tabuti1.'TEMPS'.ietat ;
  258. temps2 = tabuti2.'TEMPS'.ietat ;
  259. SI ((ABS (temps1-temps2)) > 1.0E-9) ;
  260. ierr = 1 ;
  261. QUITTER BECART ;
  262. FINSI ;
  263.  
  264. * Extraction des resultats de calcul : Modele utilisateur no21---------
  265.  
  266. * Deplacements
  267. uxu1.1 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  268. uxu1.2 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  269. uxu1.3 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  270. uxu1.4 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  271. uxu1.5 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  272. uxu1.6 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  273. uxu1.7 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  274. uxu1.8 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  275.  
  276. uyu1.1 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  277. uyu1.2 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  278. uyu1.3 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  279. uyu1.4 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  280. uyu1.5 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  281. uyu1.6 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  282. uyu1.7 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  283. uyu1.8 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  284.  
  285. uzu1.1 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  286. uzu1.2 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  287. uzu1.3 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  288. uzu1.4 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  289. uzu1.5 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  290. uzu1.6 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  291. uzu1.7 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  292. uzu1.8 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  293.  
  294. * Deformations totales
  295. igau = 0 ;
  296. REPETER BEPXXU1 nbgau ;
  297. igau = igau + 1 ;
  298. epxxu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau;
  299. FIN BEPXXU1 ;
  300.  
  301. igau = 0 ;
  302. REPETER BEPYYU1 nbgau ;
  303. igau = igau + 1 ;
  304. epyyu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau;
  305. FIN BEPYYU1 ;
  306.  
  307. igau = 0 ;
  308. REPETER BEPZZU1 nbgau ;
  309. igau = igau + 1 ;
  310. epzzu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau;
  311. FIN BEPZZU1 ;
  312.  
  313. igau = 0 ;
  314. REPETER BGAXYU1 nbgau ;
  315. igau = igau + 1 ;
  316. gaxyu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau;
  317. FIN BGAXYU1 ;
  318.  
  319. igau = 0 ;
  320. REPETER BGAXZU1 nbgau ;
  321. igau = igau + 1 ;
  322. gaxzu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau;
  323. FIN BGAXZU1 ;
  324.  
  325. igau = 0 ;
  326. REPETER BGAYZU1 nbgau ;
  327. igau = igau + 1 ;
  328. gayzu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau;
  329. FIN BGAYZU1 ;
  330.  
  331. * Deformations inelastiques
  332. igau = 0 ;
  333. REPETER BEIXXU1 nbgau ;
  334. igau = igau + 1 ;
  335. eixxu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  336. 'EIXX' 1 1 igau ;
  337. FIN BEIXXU1 ;
  338.  
  339. igau = 0 ;
  340. REPETER BEIYYU1 nbgau ;
  341. igau = igau + 1 ;
  342. eiyyu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  343. 'EIYY' 1 1 igau ;
  344. FIN BEIYYU1 ;
  345.  
  346. igau = 0 ;
  347. REPETER BEIZZU1 nbgau ;
  348. igau = igau + 1 ;
  349. eizzu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  350. 'EIZZ' 1 1 igau ;
  351. FIN BEIZZU1 ;
  352.  
  353. igau = 0 ;
  354. REPETER BGIXYU1 nbgau ;
  355. igau = igau + 1 ;
  356. gixyu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  357. 'GIXY' 1 1 igau ;
  358. FIN BGIXYU1 ;
  359.  
  360. igau = 0 ;
  361. REPETER BGIXZU1 nbgau ;
  362. igau = igau + 1 ;
  363. gixzu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  364. 'GIXZ' 1 1 igau ;
  365. FIN BGIXZU1 ;
  366.  
  367. igau = 0 ;
  368. REPETER BGIYZU1 nbgau ;
  369. igau = igau + 1 ;
  370. giyzu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  371. 'GIYZ' 1 1 igau ;
  372. FIN BGIYZU1 ;
  373.  
  374. * Contraintes
  375. igau = 0 ;
  376. REPETER BSMXXU1 nbgau ;
  377. igau = igau + 1 ;
  378. smxxu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  379. FIN BSMXXU1 ;
  380.  
  381. igau = 0 ;
  382. REPETER BSMYYU1 nbgau ;
  383. igau = igau + 1 ;
  384. smyyu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  385. FIN BSMYYU1 ;
  386.  
  387. igau = 0 ;
  388. REPETER BSMZZU1 nbgau ;
  389. igau = igau + 1 ;
  390. smzzu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  391. FIN BSMZZU1 ;
  392.  
  393. igau = 0 ;
  394. REPETER BSMXYU1 nbgau ;
  395. igau = igau + 1 ;
  396. smxyu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  397. FIN BSMXYU1 ;
  398.  
  399. igau = 0 ;
  400. REPETER BSMXZU1 nbgau ;
  401. igau = igau + 1 ;
  402. smxzu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  403. FIN BSMXZU1 ;
  404.  
  405. igau = 0 ;
  406. REPETER BSMYZU1 nbgau ;
  407. igau = igau + 1 ;
  408. smyzu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  409. FIN BSMYZU1 ;
  410.  
  411. * Variables internes
  412. igau = 0 ;
  413. REPETER BEC0U1 nbgau ;
  414. igau = igau + 1 ;
  415. ec0u1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  416. 'EC0 ' 1 1 igau ;
  417. FIN BEC0U1 ;
  418.  
  419. igau = 0 ;
  420. REPETER BESW0U1 nbgau ;
  421. igau = igau + 1 ;
  422. esw0u1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  423. 'ESW0' 1 1 igau ;
  424. FIN BESW0U1 ;
  425.  
  426. igau = 0 ;
  427. REPETER BPU1 nbgau ;
  428. igau = igau + 1 ;
  429. pu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  430. 'P ' 1 1 igau ;
  431. FIN BPU1 ;
  432.  
  433. igau = 0 ;
  434. REPETER BQTLDU1 nbgau ;
  435. igau = igau + 1 ;
  436. qtldu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  437. 'QTLD' 1 1 igau ;
  438. FIN BQTLDU1 ;
  439.  
  440. * Extraction des resultats de calcul : Modele utilisateur no23---------
  441.  
  442. * Deplacements
  443. uxu2.1 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  444. uxu2.2 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  445. uxu2.3 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  446. uxu2.4 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  447. uxu2.5 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  448. uxu2.6 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  449. uxu2.7 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  450. uxu2.8 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  451.  
  452. uyu2.1 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  453. uyu2.2 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  454. uyu2.3 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  455. uyu2.4 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  456. uyu2.5 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  457. uyu2.6 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  458. uyu2.7 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  459. uyu2.8 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  460.  
  461. uzu2.1 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  462. uzu2.2 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  463. uzu2.3 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  464. uzu2.4 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  465. uzu2.5 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  466. uzu2.6 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  467. uzu2.7 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  468. uzu2.8 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  469.  
  470. * Deformations totales
  471. igau = 0 ;
  472. REPETER BEPXXU2 nbgau ;
  473. igau = igau + 1 ;
  474. epxxu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau;
  475. FIN BEPXXU2 ;
  476.  
  477. igau = 0 ;
  478. REPETER BEPYYU2 nbgau ;
  479. igau = igau + 1 ;
  480. epyyu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau;
  481. FIN BEPYYU2 ;
  482.  
  483. igau = 0 ;
  484. REPETER BEPZZU2 nbgau ;
  485. igau = igau + 1 ;
  486. epzzu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau;
  487. FIN BEPZZU2 ;
  488.  
  489. igau = 0 ;
  490. REPETER BGAXYU2 nbgau ;
  491. igau = igau + 1 ;
  492. gaxyu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau;
  493. FIN BGAXYU2 ;
  494.  
  495. igau = 0 ;
  496. REPETER BGAXZU2 nbgau ;
  497. igau = igau + 1 ;
  498. gaxzu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau;
  499. FIN BGAXZU2 ;
  500.  
  501. igau = 0 ;
  502. REPETER BGAYZU2 nbgau ;
  503. igau = igau + 1 ;
  504. gayzu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau;
  505. FIN BGAYZU2 ;
  506.  
  507. * Deformations inelastiques
  508. igau = 0 ;
  509. REPETER BEIXXU2 nbgau ;
  510. igau = igau + 1 ;
  511. eixxu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  512. 'EIXX' 1 1 igau ;
  513. FIN BEIXXU2 ;
  514.  
  515. igau = 0 ;
  516. REPETER BEIYYU2 nbgau ;
  517. igau = igau + 1 ;
  518. eiyyu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  519. 'EIYY' 1 1 igau ;
  520. FIN BEIYYU2 ;
  521.  
  522. igau = 0 ;
  523. REPETER BEIZZU2 nbgau ;
  524. igau = igau + 1 ;
  525. eizzu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  526. 'EIZZ' 1 1 igau ;
  527. FIN BEIZZU2 ;
  528.  
  529. igau = 0 ;
  530. REPETER BGIXYU2 nbgau ;
  531. igau = igau + 1 ;
  532. gixyu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  533. 'GIXY' 1 1 igau ;
  534. FIN BGIXYU2 ;
  535.  
  536. igau = 0 ;
  537. REPETER BGIXZU2 nbgau ;
  538. igau = igau + 1 ;
  539. gixzu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  540. 'GIXZ' 1 1 igau ;
  541. FIN BGIXZU2 ;
  542.  
  543. igau = 0 ;
  544. REPETER BGIYZU2 nbgau ;
  545. igau = igau + 1 ;
  546. giyzu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  547. 'GIYZ' 1 1 igau ;
  548. FIN BGIYZU2 ;
  549.  
  550. * Contraintes
  551. igau = 0 ;
  552. REPETER BSMXXU2 nbgau ;
  553. igau = igau + 1 ;
  554. smxxu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  555. FIN BSMXXU2 ;
  556.  
  557. igau = 0 ;
  558. REPETER BSMYYU2 nbgau ;
  559. igau = igau + 1 ;
  560. smyyu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  561. FIN BSMYYU2 ;
  562.  
  563. igau = 0 ;
  564. REPETER BSMZZU2 nbgau ;
  565. igau = igau + 1 ;
  566. smzzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  567. FIN BSMZZU2 ;
  568.  
  569. igau = 0 ;
  570. REPETER BSMXYU2 nbgau ;
  571. igau = igau + 1 ;
  572. smxyu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  573. FIN BSMXYU2 ;
  574.  
  575. igau = 0 ;
  576. REPETER BSMXZU2 nbgau ;
  577. igau = igau + 1 ;
  578. smxzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  579. FIN BSMXZU2 ;
  580.  
  581. igau = 0 ;
  582. REPETER BSMYZU2 nbgau ;
  583. igau = igau + 1 ;
  584. smyzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  585. FIN BSMYZU2 ;
  586.  
  587. * Variables internes
  588. igau = 0 ;
  589. REPETER BEC0U2 nbgau ;
  590. igau = igau + 1 ;
  591. ec0u2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  592. 'EC0 ' 1 1 igau ;
  593. FIN BEC0U2 ;
  594.  
  595. igau = 0 ;
  596. REPETER BESW0U2 nbgau ;
  597. igau = igau + 1 ;
  598. esw0u2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  599. 'ESW0' 1 1 igau ;
  600. FIN BESW0U2 ;
  601.  
  602. igau = 0 ;
  603. REPETER BPU2 nbgau ;
  604. igau = igau + 1 ;
  605. pu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  606. 'P ' 1 1 igau ;
  607. FIN BPU2 ;
  608.  
  609. igau = 0 ;
  610. REPETER BQTLDU2 nbgau ;
  611. igau = igau + 1 ;
  612. qtldu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  613. 'QTLD' 1 1 igau ;
  614. FIN BQTLDU2 ;
  615.  
  616. * Variables internes supplementaires
  617. igau = 0 ;
  618. REPETER BPSUPU2 nbgau ;
  619. igau = igau + 1 ;
  620. psupu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  621. 'PSUP' 1 1 igau ;
  622. FIN BPSUPU2 ;
  623.  
  624. igau = 0 ;
  625. REPETER BQSUPU2 nbgau ;
  626. igau = igau + 1 ;
  627. qsupu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  628. 'QSUP' 1 1 igau ;
  629. FIN BQSUPU2 ;
  630.  
  631. * Comparaison des deplacements-----------------------------------------
  632.  
  633. icomp = 1 ;
  634. ino = 0 ;
  635. REPETER BUX nbno ;
  636. ino = ino + 1 ;
  637. SI (EGA uxu1.ino 0.0) ;
  638. SI (NON (EGA uxu2.ino 0.0)) ;
  639. SI (EGA IMESS 1) ;
  640. Message 'Erreur absolue' uxu2.ino
  641. 'Composante ' (EXTR LDEPL icomp)
  642. ' Noeud' ino 'Etat' ietat ;
  643. Message ;
  644. FINSI ;
  645. ************ierr = 10 ; errabs = uxu2.ino ; QUITTER BECART ;
  646. FINSI ;
  647. SINON ;
  648. dux = ABS (uxu1.ino - uxu2.ino) ;
  649. errrel = ABS (dux / uxu1.ino) ;
  650. SI (errrel > critere) ;
  651. SI (EGA IMESS 1) ;
  652. Message 'Erreur relative' errrel
  653. 'Composante ' (EXTR LDEPL icomp)
  654. ' Noeud' ino 'Etat' ietat ;
  655. Message 'uxu1.ino =' uxu1.ino 'uxu2.ino =' uxu2.ino ;
  656. Message ;
  657. FINSI ;
  658. ************ierr = 11 ; QUITTER BECART ;
  659. FINSI ;
  660. FINSI ;
  661. FIN BUX ;
  662.  
  663. icomp = 2 ;
  664. ino = 0 ;
  665. REPETER BUY nbno ;
  666. ino = ino + 1 ;
  667. SI (EGA uyu1.ino 0.0) ;
  668. SI (NON (EGA uyu2.ino 0.0)) ;
  669. SI (EGA IMESS 1) ;
  670. Message 'Erreur absolue' uyu2.ino
  671. 'Composante ' (EXTR LDEPL icomp)
  672. ' Noeud' ino 'Etat' ietat ;
  673. Message ;
  674. FINSI ;
  675. ************ierr = 10 ; errabs = uyu2.ino ; QUITTER BECART ;
  676. FINSI ;
  677. SINON ;
  678. duy = ABS (uyu1.ino - uyu2.ino) ;
  679. errrel = ABS (duy / uyu1.ino) ;
  680. SI (errrel > critere) ;
  681. SI (EGA IMESS 1) ;
  682. Message 'Erreur relative' errrel
  683. 'Composante ' (EXTR LDEPL icomp)
  684. ' Noeud' ino 'Etat' ietat ;
  685. Message 'uyu1.ino =' uyu1.ino 'uyu2.ino =' uyu2.ino ;
  686. Message ;
  687. FINSI ;
  688. ************ierr = 11 ; QUITTER BECART ;
  689. FINSI ;
  690. FINSI ;
  691. FIN BUY ;
  692.  
  693. icomp = 3 ;
  694. ino = 0 ;
  695. REPETER BUZ nbno ;
  696. ino = ino + 1 ;
  697. SI (EGA uzu1.ino 0.0) ;
  698. SI (NON (EGA uzu2.ino 0.0)) ;
  699. SI (EGA IMESS 1) ;
  700. Message 'Erreur absolue' uzu2.ino
  701. 'Composante ' (EXTR LDEPL icomp)
  702. ' Noeud' ino 'Etat' ietat ;
  703. Message ;
  704. FINSI ;
  705. ************ierr = 10 ; errabs = uzu2.ino ; QUITTER BECART ;
  706. FINSI ;
  707. SINON ;
  708. duz = ABS (uzu1.ino - uzu2.ino) ;
  709. errrel = ABS (duz / uzu1.ino) ;
  710. SI (errrel > critere) ;
  711. SI (EGA IMESS 1) ;
  712. Message 'Erreur relative' errrel
  713. 'Composante ' (EXTR LDEPL icomp)
  714. ' Noeud' ino 'Etat' ietat ;
  715. Message 'uzu1.ino =' uzu1.ino 'uzu2.ino =' uzu2.ino ;
  716. Message ;
  717. FINSI ;
  718. ************ierr = 11 ; QUITTER BECART ;
  719. FINSI ;
  720. FINSI ;
  721. FIN BUZ ;
  722.  
  723. * Comparaison des deformations-----------------------------------------
  724.  
  725. icomp = 1 ;
  726. igau = 0 ;
  727. REPETER BEPXX nbgau ;
  728. igau = igau + 1 ;
  729. SI (EGA epxxu1.igau 0.0) ;
  730. SI (NON (EGA epxxu2.igau 0.0)) ;
  731. SI (EGA IMESS 1) ;
  732. Message 'Erreur absolue' epxxu2.igau
  733. 'Composante ' (EXTR LDEFO icomp)
  734. ' Point Gauss' igau 'Etat' ietat ;
  735. Message ;
  736. FINSI ;
  737. ************ierr = 20 ; errabs = epxxu2.igau ; QUITTER BECART ;
  738. FINSI ;
  739. SINON ;
  740. depxx = ABS (epxxu1.igau - epxxu2.igau) ;
  741. errrel = ABS (depxx / epxxu1.igau) ;
  742. SI (errrel > critere) ;
  743. SI (EGA IMESS 1) ;
  744. Message 'Erreur relative' errrel
  745. 'Composante ' (EXTR LDEFO icomp)
  746. ' Point Gauss' igau 'Etat' ietat ;
  747. Message 'epxxu1.igau =' epxxu1.igau
  748. 'epxxu2.igau =' epxxu2.igau;
  749. Message ;
  750. FINSI ;
  751. ************ierr = 21 ; QUITTER BECART ;
  752. FINSI ;
  753. FINSI ;
  754. FIN BEPXX ;
  755.  
  756. icomp = 2 ;
  757. igau = 0 ;
  758. REPETER BEPYY nbgau ;
  759. igau = igau + 1 ;
  760. SI (EGA epyyu1.igau 0.0) ;
  761. SI (NON (EGA epyyu2.igau 0.0)) ;
  762. SI (EGA IMESS 1) ;
  763. Message 'Erreur absolue' epyyu2.igau
  764. 'Composante ' (EXTR LDEFO icomp)
  765. ' Point Gauss' igau 'Etat' ietat ;
  766. Message ;
  767. FINSI ;
  768. ************ierr = 20 ; errabs = epyyu2.igau ; QUITTER BECART ;
  769. FINSI ;
  770. SINON ;
  771. depyy = ABS (epyyu1.igau - epyyu2.igau) ;
  772. errrel = ABS (depyy / epyyu1.igau) ;
  773. SI (errrel > critere) ;
  774. SI (EGA IMESS 1) ;
  775. Message 'Erreur relative' errrel
  776. 'Composante ' (EXTR LDEFO icomp)
  777. ' Point Gauss' igau 'Etat' ietat ;
  778. Message 'epyyu1.igau =' epyyu1.igau
  779. 'epyyu2.igau =' epyyu2.igau;
  780. Message ;
  781. FINSI ;
  782. ************ierr = 21 ; QUITTER BECART ;
  783. FINSI ;
  784. FINSI ;
  785. FIN BEPYY ;
  786.  
  787. icomp = 3 ;
  788. igau = 0 ;
  789. REPETER BEPZZ nbgau ;
  790. igau = igau + 1 ;
  791. SI (EGA epzzu1.igau 0.0) ;
  792. SI (NON (EGA epzzu2.igau 0.0)) ;
  793. SI (EGA IMESS 1) ;
  794. Message 'Erreur absolue' epzzu2.igau
  795. 'Composante ' (EXTR LDEFO icomp)
  796. ' Point Gauss' igau 'Etat' ietat ;
  797. Message ;
  798. FINSI ;
  799. ************ierr = 20 ; errabs = epzzu2.igau ; QUITTER BECART ;
  800. FINSI ;
  801. SINON ;
  802. depzz = ABS (epzzu1.igau - epzzu2.igau) ;
  803. errrel = ABS (depzz / epzzu1.igau) ;
  804. SI (errrel > critere) ;
  805. SI (EGA IMESS 1) ;
  806. Message 'Erreur relative' errrel
  807. 'Composante ' (EXTR LDEFO icomp)
  808. ' Point Gauss' igau 'Etat' ietat ;
  809. Message 'epzzu1.igau =' epzzu1.igau
  810. 'epzzu2.igau =' epzzu2.igau;
  811. Message ;
  812. FINSI ;
  813. ************ierr = 21 ; QUITTER BECART ;
  814. FINSI ;
  815. FINSI ;
  816. FIN BEPZZ ;
  817.  
  818. icomp = 4 ;
  819. igau = 0 ;
  820. REPETER BGAXY nbgau ;
  821. igau = igau + 1 ;
  822. SI (EGA gaxyu1.igau 0.0) ;
  823. SI (NON (EGA gaxyu2.igau 0.0)) ;
  824. ************SI (EGA IMESS 1) ;
  825. ************Message 'Erreur absolue' gaxyu2.igau
  826. ************ 'Composante ' (EXTR LDEFO icomp)
  827. ************ ' Point Gauss' igau 'Etat' ietat ;
  828. ************Message ;
  829. ************FINSI ;
  830. ierr = 20 ; errabs = gaxyu2.igau ; QUITTER BECART ;
  831. FINSI ;
  832. SINON ;
  833. dgaxy = ABS (gaxyu1.igau - gaxyu2.igau) ;
  834. errrel = ABS (dgaxy / gaxyu1.igau) ;
  835. SI (errrel > critere) ;
  836. ************SI (EGA IMESS 1) ;
  837. ************Message 'Erreur relative' errrel
  838. ************ 'Composante ' (EXTR LDEFO icomp)
  839. ************ ' Point Gauss' igau 'Etat' ietat ;
  840. ************Message 'gaxyu1.igau =' gaxyu1.igau
  841. ************ 'gaxyu2.igau =' gaxyu2.igau;
  842. ************Message ;
  843. ************FINSI ;
  844. ierr = 21 ; QUITTER BECART ;
  845. FINSI ;
  846. FINSI ;
  847. FIN BGAXY ;
  848.  
  849. icomp = 5 ;
  850. igau = 0 ;
  851. REPETER BGAXZ nbgau ;
  852. igau = igau + 1 ;
  853. SI (EGA gaxzu1.igau 0.0) ;
  854. SI (NON (EGA gaxzu2.igau 0.0)) ;
  855. SI (EGA IMESS 1) ;
  856. Message 'Erreur absolue' gaxzu2.igau
  857. 'Composante ' (EXTR LDEFO icomp)
  858. ' Point Gauss' igau 'Etat' ietat ;
  859. Message ;
  860. FINSI ;
  861. ************ierr = 20 ; errabs = gaxzu2.igau ; QUITTER BECART ;
  862. FINSI ;
  863. SINON ;
  864. dgaxz = ABS (gaxzu1.igau - gaxzu2.igau) ;
  865. errrel = ABS (dgaxz / gaxzu1.igau) ;
  866. SI (errrel > critere) ;
  867. SI (EGA IMESS 1) ;
  868. Message 'Erreur relative' errrel
  869. 'Composante ' (EXTR LDEFO icomp)
  870. ' Point Gauss' igau 'Etat' ietat ;
  871. Message 'gaxzu1.igau =' gaxzu1.igau
  872. 'gaxzu2.igau =' gaxzu2.igau;
  873. Message ;
  874. FINSI ;
  875. ************ierr = 21 ; QUITTER BECART ;
  876. FINSI ;
  877. FINSI ;
  878. FIN BGAXZ ;
  879.  
  880. icomp = 6 ;
  881. igau = 0 ;
  882. REPETER BGAYZ nbgau ;
  883. igau = igau + 1 ;
  884. SI (EGA gayzu1.igau 0.0) ;
  885. SI (NON (EGA gayzu2.igau 0.0)) ;
  886. SI (EGA IMESS 1) ;
  887. Message 'Erreur absolue' gayzu2.igau
  888. 'Composante ' (EXTR LDEFO icomp)
  889. ' Point Gauss' igau 'Etat' ietat ;
  890. Message ;
  891. FINSI ;
  892. ************ierr = 20 ; errabs = gayzu2.igau ; QUITTER BECART ;
  893. FINSI ;
  894. SINON ;
  895. dgayz = ABS (gayzu1.igau - gayzu2.igau) ;
  896. errrel = ABS (dgayz / gayzu1.igau) ;
  897. SI (errrel > critere) ;
  898. SI (EGA IMESS 1) ;
  899. Message 'Erreur relative' errrel
  900. 'Composante ' (EXTR LDEFO icomp)
  901. ' Point Gauss' igau 'Etat' ietat ;
  902. Message 'gayzu1.igau =' gayzu1.igau
  903. 'gayzu2.igau =' gayzu2.igau;
  904. Message ;
  905. FINSI ;
  906. ************ierr = 21 ; QUITTER BECART ;
  907. FINSI ;
  908. FINSI ;
  909. FIN BGAYZ ;
  910.  
  911. * Comparaison des deformations inelastiques----------------------------
  912.  
  913. icomp = 1 ;
  914. igau = 0 ;
  915. REPETER BEIXX nbgau ;
  916. igau = igau + 1 ;
  917. SI (EGA eixxu1.igau 0.0) ;
  918. SI (NON (EGA eixxu2.igau 0.0)) ;
  919. SI (EGA IMESS 1) ;
  920. Message 'Erreur absolue' eixxu2.igau
  921. 'Composante ' (EXTR LDEIN icomp)
  922. ' Point Gauss' igau 'Etat' ietat ;
  923. Message ;
  924. FINSI ;
  925. ************ierr = 30 ; errabs = eixxu2.igau ; QUITTER BECART ;
  926. FINSI ;
  927. SINON ;
  928. deixx = ABS (eixxu1.igau - eixxu2.igau) ;
  929. errrel = ABS (deixx / eixxu1.igau) ;
  930. SI (errrel > critere) ;
  931. SI (EGA IMESS 1) ;
  932. Message 'Erreur relative' errrel
  933. 'Composante ' (EXTR LDEIN icomp)
  934. ' Point Gauss' igau 'Etat' ietat ;
  935. Message 'eixxu1.igau =' eixxu1.igau
  936. 'eixxu2.igau =' eixxu2.igau;
  937. Message ;
  938. FINSI ;
  939. ************ierr = 31 ; QUITTER BECART ;
  940. FINSI ;
  941. FINSI ;
  942. FIN BEIXX ;
  943.  
  944. icomp = 2 ;
  945. igau = 0 ;
  946. REPETER BEIYY nbgau ;
  947. igau = igau + 1 ;
  948. SI (EGA eiyyu1.igau 0.0) ;
  949. SI (NON (EGA eiyyu2.igau 0.0)) ;
  950. SI (EGA IMESS 1) ;
  951. Message 'Erreur absolue' eiyyu2.igau
  952. 'Composante ' (EXTR LDEIN icomp)
  953. ' Point Gauss' igau 'Etat' ietat ;
  954. Message ;
  955. FINSI ;
  956. ************ierr = 30 ; errabs = eiyyu2.igau ; QUITTER BECART ;
  957. FINSI ;
  958. SINON ;
  959. deiyy = ABS (eiyyu1.igau - eiyyu2.igau) ;
  960. errrel = ABS (deiyy / eiyyu1.igau) ;
  961. SI (errrel > critere) ;
  962. SI (EGA IMESS 1) ;
  963. Message 'Erreur relative' errrel
  964. 'Composante ' (EXTR LDEIN icomp)
  965. ' Point Gauss' igau 'Etat' ietat ;
  966. Message 'eiyyu1.igau =' eiyyu1.igau
  967. 'eiyyu2.igau =' eiyyu2.igau;
  968. Message ;
  969. FINSI ;
  970. ************ierr = 31 ; QUITTER BECART ;
  971. FINSI ;
  972. FINSI ;
  973. FIN BEIYY ;
  974.  
  975. icomp = 3 ;
  976. igau = 0 ;
  977. REPETER BEIZZ nbgau ;
  978. igau = igau + 1 ;
  979. SI (EGA eizzu1.igau 0.0) ;
  980. SI (NON (EGA eizzu2.igau 0.0)) ;
  981. SI (EGA IMESS 1) ;
  982. Message 'Erreur absolue' eizzu2.igau
  983. 'Composante ' (EXTR LDEIN icomp)
  984. ' Point Gauss' igau 'Etat' ietat ;
  985. Message ;
  986. FINSI ;
  987. ************ierr = 30 ; errabs = eizzu2.igau ; QUITTER BECART ;
  988. FINSI ;
  989. SINON ;
  990. deizz = ABS (eizzu1.igau - eizzu2.igau) ;
  991. errrel = ABS (deizz / eizzu1.igau) ;
  992. SI (errrel > critere) ;
  993. SI (EGA IMESS 1) ;
  994. Message 'Erreur relative' errrel
  995. 'Composante ' (EXTR LDEIN icomp)
  996. ' Point Gauss' igau 'Etat' ietat ;
  997. Message 'eizzu1.igau =' eizzu1.igau
  998. 'eizzu2.igau =' eizzu2.igau;
  999. Message ;
  1000. FINSI ;
  1001. ************ierr = 31 ; QUITTER BECART ;
  1002. FINSI ;
  1003. FINSI ;
  1004. FIN BEIZZ ;
  1005.  
  1006. icomp = 4 ;
  1007. igau = 0 ;
  1008. REPETER BGIXY nbgau ;
  1009. igau = igau + 1 ;
  1010. SI (EGA gixyu1.igau 0.0) ;
  1011. SI (NON (EGA gixyu2.igau 0.0)) ;
  1012. ************SI (EGA IMESS 1) ;
  1013. ************Message 'Erreur absolue' gixyu2.igau
  1014. ************ 'Composante ' (EXTR LDEIN icomp)
  1015. ************ ' Point Gauss' igau 'Etat' ietat ;
  1016. ************Message ;
  1017. ************FINSI ;
  1018. ierr = 30 ; errabs = gixyu2.igau ; QUITTER BECART ;
  1019. FINSI ;
  1020. SINON ;
  1021. dgixy = ABS (gixyu1.igau - gixyu2.igau) ;
  1022. errrel = ABS (dgixy / gixyu1.igau) ;
  1023. SI (errrel > critere) ;
  1024. ************SI (EGA IMESS 1) ;
  1025. ************Message 'Erreur relative' errrel
  1026. ************ 'Composante ' (EXTR LDEIN icomp)
  1027. ************ ' Point Gauss' igau 'Etat' ietat ;
  1028. ************Message 'gixyu1.igau =' gixyu1.igau
  1029. ************ 'gixyu2.igau =' gixyu2.igau;
  1030. ************Message ;
  1031. ************FINSI ;
  1032. ierr = 31 ; QUITTER BECART ;
  1033. FINSI ;
  1034. FINSI ;
  1035. FIN BGIXY ;
  1036.  
  1037. icomp = 5 ;
  1038. igau = 0 ;
  1039. REPETER BGIXZ nbgau ;
  1040. igau = igau + 1 ;
  1041. SI (EGA gixzu1.igau 0.0) ;
  1042. SI (NON (EGA gixzu2.igau 0.0)) ;
  1043. SI (EGA IMESS 1) ;
  1044. Message 'Erreur absolue' gixzu2.igau
  1045. 'Composante ' (EXTR LDEIN icomp)
  1046. ' Point Gauss' igau 'Etat' ietat ;
  1047. Message ;
  1048. FINSI ;
  1049. ************ierr = 30 ; errabs = gixzu2.igau ; QUITTER BECART ;
  1050. FINSI ;
  1051. SINON ;
  1052. dgixz = ABS (gixzu1.igau - gixzu2.igau) ;
  1053. errrel = ABS (dgixz / gixzu1.igau) ;
  1054. SI (errrel > critere) ;
  1055. SI (EGA IMESS 1) ;
  1056. Message 'Erreur relative' errrel
  1057. 'Composante ' (EXTR LDEIN icomp)
  1058. ' Point Gauss' igau 'Etat' ietat ;
  1059. Message 'gixzu1.igau =' gixzu1.igau
  1060. 'gixzu2.igau =' gixzu2.igau;
  1061. Message ;
  1062. FINSI ;
  1063. ************ierr = 31 ; QUITTER BECART ;
  1064. FINSI ;
  1065. FINSI ;
  1066. FIN BGIXZ ;
  1067.  
  1068. icomp = 6 ;
  1069. igau = 0 ;
  1070. REPETER BGIYZ nbgau ;
  1071. igau = igau + 1 ;
  1072. SI (EGA giyzu1.igau 0.0) ;
  1073. SI (NON (EGA giyzu2.igau 0.0)) ;
  1074. SI (EGA IMESS 1) ;
  1075. Message 'Erreur absolue' giyzu2.igau
  1076. 'Composante ' (EXTR LDEIN icomp)
  1077. ' Point Gauss' igau 'Etat' ietat ;
  1078. Message ;
  1079. FINSI ;
  1080. ************ierr = 30 ; errabs = giyzu2.igau ; QUITTER BECART ;
  1081. FINSI ;
  1082. SINON ;
  1083. dgiyz = ABS (giyzu1.igau - giyzu2.igau) ;
  1084. errrel = ABS (dgiyz / giyzu1.igau) ;
  1085. SI (errrel > critere) ;
  1086. SI (EGA IMESS 1) ;
  1087. Message 'Erreur relative' errrel
  1088. 'Composante ' (EXTR LDEIN icomp)
  1089. ' Point Gauss' igau 'Etat' ietat ;
  1090. Message 'giyzu1.igau =' giyzu1.igau
  1091. 'giyzu2.igau =' giyzu2.igau;
  1092. Message ;
  1093. FINSI ;
  1094. ************ierr = 31 ; QUITTER BECART ;
  1095. FINSI ;
  1096. FINSI ;
  1097. FIN BGIYZ ;
  1098.  
  1099. * Comparaison des contraintes------------------------------------------
  1100.  
  1101. icomp = 1 ;
  1102. igau = 0 ;
  1103. REPETER BSMXX nbgau ;
  1104. igau = igau + 1 ;
  1105. SI (EGA smxxu1.igau 0.0) ;
  1106. SI (NON (EGA smxxu2.igau 0.0)) ;
  1107. SI (EGA IMESS 1) ;
  1108. Message 'Erreur absolue' smxxu2.igau
  1109. 'Composante ' (EXTR LCONT icomp)
  1110. ' Point Gauss' igau 'Etat' ietat ;
  1111. Message ;
  1112. FINSI ;
  1113. ************ierr = 40 ; errabs = smxxu2.igau ; QUITTER BECART ;
  1114. FINSI ;
  1115. SINON ;
  1116. dsmxx = ABS (smxxu1.igau - smxxu2.igau) ;
  1117. errrel = ABS (dsmxx / smxxu1.igau) ;
  1118. SI (errrel > critere) ;
  1119. SI (EGA IMESS 1) ;
  1120. Message 'Erreur relative' errrel
  1121. 'Composante ' (EXTR LCONT icomp)
  1122. ' Point Gauss' igau 'Etat' ietat ;
  1123. Message 'smxxu1.igau =' smxxu1.igau
  1124. 'smxxu2.igau =' smxxu2.igau;
  1125. Message ;
  1126. FINSI ;
  1127. ************ierr = 41 ; QUITTER BECART ;
  1128. FINSI ;
  1129. FINSI ;
  1130. FIN BSMXX ;
  1131.  
  1132. icomp = 2 ;
  1133. igau = 0 ;
  1134. REPETER BSMYY nbgau ;
  1135. igau = igau + 1 ;
  1136. SI (EGA smyyu1.igau 0.0) ;
  1137. SI (NON (EGA smyyu2.igau 0.0)) ;
  1138. SI (EGA IMESS 1) ;
  1139. Message 'Erreur absolue' smyyu2.igau
  1140. 'Composante ' (EXTR LCONT icomp)
  1141. ' Point Gauss' igau 'Etat' ietat ;
  1142. Message ;
  1143. FINSI ;
  1144. ************ierr = 40 ; errabs = smyyu2.igau ; QUITTER BECART ;
  1145. FINSI ;
  1146. SINON ;
  1147. dsmyy = ABS (smyyu1.igau - smyyu2.igau) ;
  1148. errrel = ABS (dsmyy / smyyu1.igau) ;
  1149. SI (errrel > critere) ;
  1150. SI (EGA IMESS 1) ;
  1151. Message 'Erreur relative' errrel
  1152. 'Composante ' (EXTR LCONT icomp)
  1153. ' Point Gauss' igau 'Etat' ietat ;
  1154. Message 'smyyu1.igau =' smyyu1.igau
  1155. 'smyyu2.igau =' smyyu2.igau;
  1156. Message ;
  1157. FINSI ;
  1158. ************ierr = 41 ; QUITTER BECART ;
  1159. FINSI ;
  1160. FINSI ;
  1161. FIN BSMYY ;
  1162.  
  1163. icomp = 3 ;
  1164. igau = 0 ;
  1165. REPETER BSMZZ nbgau ;
  1166. igau = igau + 1 ;
  1167. SI (EGA smzzu1.igau 0.0) ;
  1168. SI (NON (EGA smzzu2.igau 0.0)) ;
  1169. SI (EGA IMESS 1) ;
  1170. Message 'Erreur absolue' smzzu2.igau
  1171. 'Composante ' (EXTR LCONT icomp)
  1172. ' Point Gauss' igau 'Etat' ietat ;
  1173. Message ;
  1174. FINSI ;
  1175. ************ierr = 40 ; errabs = smzzu2.igau ; QUITTER BECART ;
  1176. FINSI ;
  1177. SINON ;
  1178. dsmzz = ABS (smzzu1.igau - smzzu2.igau) ;
  1179. errrel = ABS (dsmzz / smzzu1.igau) ;
  1180. SI (errrel > critere) ;
  1181. SI (EGA IMESS 1) ;
  1182. Message 'Erreur relative' errrel
  1183. 'Composante ' (EXTR LCONT icomp)
  1184. ' Point Gauss' igau 'Etat' ietat ;
  1185. Message 'smzzu1.igau =' smzzu1.igau
  1186. 'smzzu2.igau =' smzzu2.igau;
  1187. Message ;
  1188. FINSI ;
  1189. ************ierr = 41 ; QUITTER BECART ;
  1190. FINSI ;
  1191. FINSI ;
  1192. FIN BSMZZ ;
  1193.  
  1194. icomp = 4 ;
  1195. igau = 0 ;
  1196. REPETER BSMXY nbgau ;
  1197. igau = igau + 1 ;
  1198. SI (EGA smxyu1.igau 0.0) ;
  1199. SI (NON (EGA smxyu2.igau 0.0)) ;
  1200. ************SI (EGA IMESS 1) ;
  1201. ************Message 'Erreur absolue' smxyu2.igau
  1202. ************ 'Composante ' (EXTR LCONT icomp)
  1203. ************ ' Point Gauss' igau 'Etat' ietat ;
  1204. ************Message ;
  1205. ************FINSI ;
  1206. ierr = 40 ; errabs = smxyu2.igau ; QUITTER BECART ;
  1207. FINSI ;
  1208. SINON ;
  1209. dsmxy = ABS (smxyu1.igau - smxyu2.igau) ;
  1210. errrel = ABS (dsmxy / smxyu1.igau) ;
  1211. SI (errrel > critere) ;
  1212. ************SI (EGA IMESS 1) ;
  1213. ************Message 'Erreur relative' errrel
  1214. ************ 'Composante ' (EXTR LCONT icomp)
  1215. ************ ' Point Gauss' igau 'Etat' ietat ;
  1216. ************Message 'smxyu1.igau =' smxyu1.igau
  1217. ************ 'smxyu2.igau =' smxyu2.igau;
  1218. ************Message ;
  1219. ************FINSI ;
  1220. ierr = 41 ; QUITTER BECART ;
  1221. FINSI ;
  1222. FINSI ;
  1223. FIN BSMXY ;
  1224.  
  1225. icomp = 5 ;
  1226. igau = 0 ;
  1227. REPETER BSMXZ nbgau ;
  1228. igau = igau + 1 ;
  1229. SI (EGA smxzu1.igau 0.0) ;
  1230. SI (NON (EGA smxzu2.igau 0.0)) ;
  1231. SI (EGA IMESS 1) ;
  1232. Message 'Erreur absolue' smxzu2.igau
  1233. 'Composante ' (EXTR LCONT icomp)
  1234. ' Point Gauss' igau 'Etat' ietat ;
  1235. Message ;
  1236. FINSI ;
  1237. ************ierr = 40 ; errabs = smxzu2.igau ; QUITTER BECART ;
  1238. FINSI ;
  1239. SINON ;
  1240. dsmxz = ABS (smxzu1.igau - smxzu2.igau) ;
  1241. errrel = ABS (dsmxz / smxzu1.igau) ;
  1242. SI (errrel > critere) ;
  1243. SI (EGA IMESS 1) ;
  1244. Message 'Erreur relative' errrel
  1245. 'Composante ' (EXTR LCONT icomp)
  1246. ' Point Gauss' igau 'Etat' ietat ;
  1247. Message 'smxzu1.igau =' smxzu1.igau
  1248. 'smxzu2.igau =' smxzu2.igau;
  1249. Message ;
  1250. FINSI ;
  1251. ************ierr = 41 ; QUITTER BECART ;
  1252. FINSI ;
  1253. FINSI ;
  1254. FIN BSMXZ ;
  1255.  
  1256. icomp = 6 ;
  1257. igau = 0 ;
  1258. REPETER BSMYZ nbgau ;
  1259. igau = igau + 1 ;
  1260. SI (EGA smyzu1.igau 0.0) ;
  1261. SI (NON (EGA smyzu2.igau 0.0)) ;
  1262. SI (EGA IMESS 1) ;
  1263. Message 'Erreur absolue' smyzu2.igau
  1264. 'Composante ' (EXTR LCONT icomp)
  1265. ' Point Gauss' igau 'Etat' ietat ;
  1266. Message ;
  1267. FINSI ;
  1268. ************ierr = 40 ; errabs = smyzu2.igau ; QUITTER BECART ;
  1269. FINSI ;
  1270. SINON ;
  1271. dsmyz = ABS (smyzu1.igau - smyzu2.igau) ;
  1272. errrel = ABS (dsmyz / smyzu1.igau) ;
  1273. SI (errrel > critere) ;
  1274. SI (EGA IMESS 1) ;
  1275. Message 'Erreur relative' errrel
  1276. 'Composante ' (EXTR LCONT icomp)
  1277. ' Point Gauss' igau 'Etat' ietat ;
  1278. Message 'smyzu1.igau =' smyzu1.igau
  1279. 'smyzu2.igau =' smyzu2.igau;
  1280. Message ;
  1281. FINSI ;
  1282. ************ierr = 41 ; QUITTER BECART ;
  1283. FINSI ;
  1284. FINSI ;
  1285. FIN BSMYZ ;
  1286.  
  1287. * Comparaison des variables internes-----------------------------------
  1288.  
  1289. icomp = 1 ;
  1290. igau = 0 ;
  1291. REPETER BEC0 nbgau ;
  1292. igau = igau + 1 ;
  1293. SI (EGA ec0u1.igau 0.0) ;
  1294. SI (NON (EGA ec0u2.igau 0.0)) ;
  1295. ************SI (EGA IMESS 1) ;
  1296. ************Message 'Erreur absolue' ec0u2.igau
  1297. ************ 'Composante ' (EXTR LVARI icomp)
  1298. ************ ' Point Gauss' igau 'Etat' ietat ;
  1299. ************Message ;
  1300. ************FINSI ;
  1301. ierr = 50 ; errabs = ec0u2.igau ; QUITTER BECART ;
  1302. FINSI ;
  1303. SINON ;
  1304. dec0 = ABS (ec0u1.igau - ec0u2.igau) ;
  1305. errrel = ABS (dec0 / ec0u1.igau) ;
  1306. SI (errrel > critere) ;
  1307. ************SI (EGA IMESS 1) ;
  1308. ************Message 'Erreur relative' errrel
  1309. ************ 'Composante ' (EXTR LVARI icomp)
  1310. ************ ' Point Gauss' igau 'Etat' ietat ;
  1311. ************Message 'ec0u1.igau =' ec0u1.igau
  1312. ************ 'ec0u2.igau =' ec0u2.igau;
  1313. ************Message ;
  1314. ************FINSI ;
  1315. ierr = 51 ; QUITTER BECART ;
  1316. FINSI ;
  1317. FINSI ;
  1318. FIN BEC0 ;
  1319.  
  1320. icomp = 2 ;
  1321. igau = 0 ;
  1322. REPETER BESW0 nbgau ;
  1323. igau = igau + 1 ;
  1324. SI (EGA esw0u1.igau 0.0) ;
  1325. SI (NON (EGA esw0u2.igau 0.0)) ;
  1326. ************SI (EGA IMESS 1) ;
  1327. ************Message 'Erreur absolue' esw0u2.igau
  1328. ************ 'Composante ' (EXTR LVARI icomp)
  1329. ************ ' Point Gauss' igau 'Etat' ietat ;
  1330. ************Message ;
  1331. ************FINSI ;
  1332. ierr = 50 ; errabs = esw0u2.igau ; QUITTER BECART ;
  1333. FINSI ;
  1334. SINON ;
  1335. desw0 = ABS (esw0u1.igau - esw0u2.igau) ;
  1336. errrel = ABS (desw0 / esw0u1.igau) ;
  1337. SI (errrel > critere) ;
  1338. ************SI (EGA IMESS 1) ;
  1339. ************Message 'Erreur relative' errrel
  1340. ************ 'Composante ' (EXTR LVARI icomp)
  1341. ************ ' Point Gauss' igau 'Etat' ietat ;
  1342. ************Message 'esw0u1.igau =' esw0u1.igau
  1343. ************ 'esw0u2.igau =' esw0u2.igau;
  1344. ************Message ;
  1345. ************FINSI ;
  1346. ierr = 51 ; QUITTER BECART ;
  1347. FINSI ;
  1348. FINSI ;
  1349. FIN BESW0 ;
  1350.  
  1351. icomp = 3 ;
  1352. igau = 0 ;
  1353. REPETER BP nbgau ;
  1354. igau = igau + 1 ;
  1355. SI (EGA pu1.igau 0.0) ;
  1356. SI (NON (EGA pu2.igau 0.0)) ;
  1357. ************SI (EGA IMESS 1) ;
  1358. ************Message 'Erreur absolue' pu2.igau
  1359. ************ 'Composante ' (EXTR LVARI icomp)
  1360. ************ ' Point Gauss' igau 'Etat' ietat ;
  1361. ************Message ;
  1362. ************FINSI ;
  1363. ierr = 50 ; errabs = pu2.igau ; QUITTER BECART ;
  1364. FINSI ;
  1365. SINON ;
  1366. dp = ABS (pu1.igau - pu2.igau) ;
  1367. errrel = ABS (dp / pu1.igau) ;
  1368. SI (errrel > critere) ;
  1369. ************SI (EGA IMESS 1) ;
  1370. ************Message 'Erreur relative' errrel
  1371. ************ 'Composante ' (EXTR LVARI icomp)
  1372. ************ ' Point Gauss' igau 'Etat' ietat ;
  1373. ************Message 'pu1.igau =' pu1.igau
  1374. ************ 'pu2.igau =' pu2.igau;
  1375. ************Message ;
  1376. ************FINSI ;
  1377. ierr = 51 ; QUITTER BECART ;
  1378. FINSI ;
  1379. FINSI ;
  1380. FIN BP ;
  1381.  
  1382. icomp = 4 ;
  1383. igau = 0 ;
  1384. REPETER BQTLD nbgau ;
  1385. igau = igau + 1 ;
  1386. SI (EGA qtldu1.igau 0.0) ;
  1387. SI (NON (EGA qtldu2.igau 0.0)) ;
  1388. ************SI (EGA IMESS 1) ;
  1389. ************Message 'Erreur absolue' qtldu2.igau
  1390. ************ 'Composante ' (EXTR LVARI icomp)
  1391. ************ ' Point Gauss' igau 'Etat' ietat ;
  1392. ************Message ;
  1393. ************FINSI ;
  1394. ierr = 50 ; errabs = qtldu2.igau ; QUITTER BECART ;
  1395. FINSI ;
  1396. SINON ;
  1397. dqtld = ABS (qtldu1.igau - qtldu2.igau) ;
  1398. errrel = ABS (dqtld / qtldu1.igau) ;
  1399. SI (errrel > critere) ;
  1400. ************SI (EGA IMESS 1) ;
  1401. ************Message 'Erreur relative' errrel
  1402. ************ 'Composante ' (EXTR LVARI icomp)
  1403. ************ ' Point Gauss' igau 'Etat' ietat ;
  1404. ************Message 'qtldu1.igau =' qtldu1.igau
  1405. ************ 'qtldu2.igau =' qtldu2.igau;
  1406. ************Message ;
  1407. ************FINSI ;
  1408. ierr = 51 ; QUITTER BECART ;
  1409. FINSI ;
  1410. FINSI ;
  1411. FIN BQTLD ;
  1412.  
  1413. * Comparaison des variables internes supplementaires (modele no23)-----
  1414.  
  1415. icomp = 1 ;
  1416. igau = 0 ;
  1417. REPETER BPSUP nbgau ;
  1418. igau = igau + 1 ;
  1419. SI (EGA pu2.igau 0.0) ;
  1420. SI (NON (EGA psupu2.igau 0.0)) ;
  1421. ************SI (EGA IMESS 1) ;
  1422. ************Message 'Erreur absolue' psupu2.igau
  1423. ************ 'Composante ' (EXTR LVSUP icomp)
  1424. ************ ' Point Gauss' igau 'Etat' ietat ;
  1425. ************Message ;
  1426. ************FINSI ;
  1427. ierr = 60 ; errabs = psupu2.igau ; QUITTER BECART ;
  1428. FINSI ;
  1429. SINON ;
  1430. dpsup = ABS (pu2.igau - psupu2.igau) ;
  1431. errrel = ABS (dpsup / pu2.igau) ;
  1432. SI (errrel > critere) ;
  1433. ************SI (EGA IMESS 1) ;
  1434. ************Message 'Erreur relative' errrel
  1435. ************ 'Composante ' (EXTR LVSUP icomp)
  1436. ************ ' Point Gauss' igau 'Etat' ietat ;
  1437. ************Message 'pu2.igau =' pu2.igau
  1438. ************ 'psupu2.igau =' psupu2.igau;
  1439. ************Message ;
  1440. ************FINSI ;
  1441. ierr = 61 ; QUITTER BECART ;
  1442. FINSI ;
  1443. FINSI ;
  1444. FIN BPSUP ;
  1445.  
  1446. icomp = 2 ;
  1447. igau = 0 ;
  1448. REPETER BQSUP nbgau ;
  1449. igau = igau + 1 ;
  1450. SI (EGA qtldu2.igau 0.0) ;
  1451. SI (NON (EGA qsupu2.igau 0.0)) ;
  1452. ************SI (EGA IMESS 1) ;
  1453. ************Message 'Erreur absolue' qsupu2.igau
  1454. ************ 'Composante ' (EXTR LVSUP icomp)
  1455. ************ ' Point Gauss' igau 'Etat' ietat ;
  1456. ************Message ;
  1457. ************FINSI ;
  1458. ierr = 60 ; errabs = qsupu2.igau ; QUITTER BECART ;
  1459. FINSI ;
  1460. SINON ;
  1461. dqsup = ABS (qtldu2.igau - qsupu2.igau) ;
  1462. errrel = ABS (dqsup / qtldu2.igau) ;
  1463. SI (errrel > critere) ;
  1464. ************SI (EGA IMESS 1) ;
  1465. ************Message 'Erreur relative' errrel
  1466. ************ 'Composante ' (EXTR LVSUP icomp)
  1467. ************ ' Point Gauss' igau 'Etat' ietat ;
  1468. ************Message 'qtldu2.igau =' qtldu2.igau
  1469. ************ 'qsupu2.igau =' qsupu2.igau;
  1470. ************Message ;
  1471. ************FINSI ;
  1472. ierr = 61 ; QUITTER BECART ;
  1473. FINSI ;
  1474. FINSI ;
  1475. FIN BQSUP ;
  1476.  
  1477. FIN BECART ;
  1478.  
  1479. * FIN BOUCLE COMPARAISON DES RESULTATS________________________________
  1480.  
  1481. * MESSAGES D'ERREUR_____________________________________________________
  1482.  
  1483. Message ;
  1484. Message 'CPU avec modele externe 21 = ' CPUuti1 'centiemes de seconde' ;
  1485. Message 'CPU avec modele externe 23 = ' CPUuti2 'centiemes de seconde' ;
  1486. Message ;
  1487.  
  1488. ifin = netat-1 ;
  1489.  
  1490. gaxyu1ma = 0.0 ;
  1491. igau = 0 ;
  1492. REPETER BGAXYFU1 nbgau ;
  1493. igau = igau + 1 ;
  1494. gaxy = EXTR (tabuti1.'DEFORMATIONS'.ifin)
  1495. 'GAXY' 1 1 igau ;
  1496. SI (gaxy > gaxyu1ma) ; gaxyu1ma = gaxy ; FINSI ;
  1497. FIN BGAXYFU1 ;
  1498. gaxyu2ma = 0.0 ;
  1499. igau = 0 ;
  1500. REPETER BGAXYFU2 nbgau ;
  1501. igau = igau + 1 ;
  1502. gaxy = EXTR (tabuti2.'DEFORMATIONS'.ifin)
  1503. 'GAXY' 1 1 igau ;
  1504. SI (gaxy > gaxyu2ma) ; gaxyu2ma = gaxy ; FINSI ;
  1505. FIN BGAXYFU2 ;
  1506. Message 'Deformation totale finale (maxi sur points de Gauss)' ;
  1507. Message 'GAXY max avec modele externe 21 = ' gaxyu1ma ;
  1508. Message 'GAXY max avec modele externe 23 = ' gaxyu2ma ;
  1509. Message ;
  1510.  
  1511. gixyu1ma = 0.0 ;
  1512. igau = 0 ;
  1513. REPETER BGIXYFU1 nbgau ;
  1514. igau = igau + 1 ;
  1515. gixy = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ifin)
  1516. 'GIXY' 1 1 igau ;
  1517. SI (gixy > gixyu1ma) ; gixyu1ma = gixy ; FINSI ;
  1518. FIN BGIXYFU1 ;
  1519. gixyu2ma = 0.0 ;
  1520. igau = 0 ;
  1521. REPETER BGIXYFU2 nbgau ;
  1522. igau = igau + 1 ;
  1523. gixy = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ifin)
  1524. 'GIXY' 1 1 igau ;
  1525. SI (gixy > gixyu2ma) ; gixyu2ma = gixy ; FINSI ;
  1526. FIN BGIXYFU2 ;
  1527. Message 'Deformation de fluage finale (maxi sur points de Gauss)' ;
  1528. Message 'GIXY max avec modele externe 21 = ' gixyu1ma ;
  1529. Message 'GIXY max avec modele externe 23 = ' gixyu2ma ;
  1530. Message ;
  1531.  
  1532. smxyu1ma = 0.0 ;
  1533. igau = 0 ;
  1534. REPETER BSMXYFU1 nbgau ;
  1535. igau = igau + 1 ;
  1536. smxy = EXTR (tabuti1.'CONTRAINTES'.ifin)
  1537. 'SMXY' 1 1 igau ;
  1538. SI (smxy > smxyu1ma) ; smxyu1ma = smxy ; FINSI ;
  1539. FIN BSMXYFU1 ;
  1540. smxyu2ma = 0.0 ;
  1541. igau = 0 ;
  1542. REPETER BSMXYFU2 nbgau ;
  1543. igau = igau + 1 ;
  1544. smxy = EXTR (tabuti2.'CONTRAINTES'.ifin)
  1545. 'SMXY' 1 1 igau ;
  1546. SI (smxy > smxyu2ma) ; smxyu2ma = smxy ; FINSI ;
  1547. FIN BSMXYFU2 ;
  1548. Message 'Contrainte finale (maxi sur points de Gauss)' ;
  1549. Message 'SMXY max avec modele externe 21 = ' smxyu1ma ;
  1550. Message 'SMXY max avec modele externe 23 = ' smxyu2ma ;
  1551. Message ;
  1552.  
  1553. SI (ierr > 0) ;
  1554.  
  1555. SI (EGA ierr 1) ;
  1556. Message 'Instants de calcul non coincidants' ;
  1557. ERRE 5 ;
  1558. FINSI ;
  1559.  
  1560. SI (EGA ierr 10) ;
  1561. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  1562. ' Noeud' ino 'Etat' ietat ;
  1563. ERRE 5 ;
  1564. FINSI ;
  1565.  
  1566. SI (EGA ierr 11) ;
  1567. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  1568. ' Noeud' ino 'Etat' ietat ;
  1569. ERRE 5 ;
  1570. FINSI ;
  1571.  
  1572. SI (EGA ierr 20) ;
  1573. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  1574. ' Point Gauss' igau 'Etat' ietat ;
  1575. ERRE 5 ;
  1576. FINSI ;
  1577.  
  1578. SI (EGA ierr 21) ;
  1579. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  1580. ' Point Gauss' igau 'Etat' ietat ;
  1581. ERRE 5 ;
  1582. FINSI ;
  1583.  
  1584. SI (EGA ierr 30) ;
  1585. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEIN icomp)
  1586. ' Point Gauss' igau 'Etat' ietat ;
  1587. ERRE 5 ;
  1588. FINSI ;
  1589.  
  1590. SI (EGA ierr 31) ;
  1591. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEIN icomp)
  1592. ' Point Gauss' igau 'Etat' ietat ;
  1593. ERRE 5 ;
  1594. FINSI ;
  1595.  
  1596. SI (EGA ierr 40) ;
  1597. Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  1598. ' Point Gauss' igau 'Etat' ietat ;
  1599. ERRE 5 ;
  1600. FINSI ;
  1601.  
  1602. SI (EGA ierr 41) ;
  1603. Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  1604. ' Point Gauss' igau 'Etat' ietat ;
  1605. ERRE 5 ;
  1606. FINSI ;
  1607.  
  1608. SI (EGA ierr 50) ;
  1609. Message 'Erreur absolue' errabs 'Composante ' (EXTR LVARI icomp)
  1610. ' Point Gauss' igau 'Etat' ietat ;
  1611. ERRE 5 ;
  1612. FINSI ;
  1613.  
  1614. SI (EGA ierr 51) ;
  1615. Message 'Erreur relative' errrel 'Composante ' (EXTR LVARI icomp)
  1616. ' Point Gauss' igau 'Etat' ietat ;
  1617. ERRE 5 ;
  1618. FINSI ;
  1619.  
  1620. SI (EGA ierr 60) ;
  1621. Message 'Erreur absolue' errabs 'Composante ' (EXTR LVSUP icomp)
  1622. ' Point Gauss' igau 'Etat' ietat ;
  1623. ERRE 5 ;
  1624. FINSI ;
  1625.  
  1626. SI (EGA ierr 61) ;
  1627. Message 'Erreur relative' errrel 'Composante ' (EXTR LVSUP icomp)
  1628. ' Point Gauss' igau 'Etat' ietat ;
  1629. ERRE 5 ;
  1630. FINSI ;
  1631.  
  1632. SINON ;
  1633.  
  1634. Message 'Resultats concordants entre les deux modeles' ;
  1635. Message 'Aucune erreur relative superieure a ' critere
  1636. ' pour les composantes' ;
  1637. Message (EXTR LDEFO idefo) ' ' (EXTR LDEIN idein) ' '
  1638. (EXTR LCONT icont) ' '
  1639. (EXTR LVARI ivar1) ' ' (EXTR LVARI ivar2) ' '
  1640. (EXTR LVARI ivar3) ' ' (EXTR LVARI ivar4) ;
  1641. Message ;
  1642. Message 'Modele externe 23 : aucune erreur relative superieure a '
  1643. critere ;
  1644. Message 'entre les composantes '
  1645. (EXTR LVARI ivar3) ' ' (EXTR LVARI ivar4) ' et '
  1646. (EXTR LVSUP ivas1) ' ' (EXTR LVSUP ivas2) ;
  1647. ERRE 0 ;
  1648.  
  1649. FINSI ;
  1650.  
  1651. FIN ;
  1652.  
  1653.  
  1654.  
  1655.  
  1656.  
  1657.  
  1658.  
  1659.  
  1660.  
  1661.  

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