Télécharger creep04_cisXY.dgibi

Retour à la liste

Numérotation des lignes :

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

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