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 &!4 ;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 (tabuti1.'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<¯span~ ;
  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. >div rtqle="font: normal normal 1em/1.2em monospace; margin:0; padding:0; background:none; vertical-align:top;"> 'GIYZ' 1 1 igau ;
  • FIN BGIYZU2 ;
  •  
  • * Contraintes
  • igau = 0 ;
  • REPETER BSMXXU2 nbgau ;
  • igau = igau + 1 ;
  • smxxu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  • FIN BSMXXU2 ;
  •  
  • igau = 0 ;
  • REPETER BSMYYU2 nbgau ;
  • igau = igau + 1 ;
  • smyyu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  • FIN BSMYYU2 ;
  •  
  • igau = 0 ;
  • REPETER BSMZZU2 nbgau ;
  • igau = igau + 1 ;
  • smzzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  • FIN BSMZZU2 ;
  •  
  • igau = 0 ;
  • REPETER BSMXYU2 nbgau ;
  • igau = igau + 1 ;
  • smxyu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  • FIN BSMXYU2 ;
  •  
  • igau = 0 ;
  • REPETER BSMXZU2 nbgau ;
  • igau = igau + 1 ;
  • smxzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  • FIN BSMXZU2 ;
  •  
  • igau = 0 ;
  • REPETER BSMYZU2 nbgau ;
  • igau = igau + 1 ;
  • smyzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  • FIN BSMYZU2 ;
  •  
  • * Variables internes
  • igau = 0 ;
  • REPETER BEC0U2 nbgau ;
  • igau = igau + 1 ;
  • ec0u2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  • 'EC0 ' 1 1 igau ;
  • FIN BEC0U2 ;
  •  
  • igau = 0 ;
  • REPETER BESW0U2 nbgau ;
  • igau = igau + 1 ;
  • esw0u2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  • 'ESW0' 1 1 igau ;
  • FIN BESW0U2 ;
  •  
  • igau = 0 ;
  • REPETER BPU2 nbgau ;
  • igau = igau + 1 ;
  • pu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  • 'P ' 1 1 igau ;
  • FIN BPU2 ;
  •  
  • igau = 0 ;
  • REPETER BQTLDU2 nbgau ;
  • igau = igau + 1 ;
  • qtldu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  • 'QTLD' 1 1 igau ;
  • FIN BQTLDU2 ;
  •  
  • * Variables internes supplementaires
  • igau = 0 ;
  • REPETER BPSUPU2 nbgau ;
  • igau = igau + 1 ;
  • psupu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  • 'PSUP' 1 1 igau ;
  • FIN BPSUPU2 ;
  •  
  • igau = 0 ;
  • REPETER BQSUPU2 nbgau ;
  • igau = igau + 1 ;
  • qsupu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  • 'QSUP' 1 1 igau ;
  • FIN BQSUPU2 ;
  •  
  • * Comparaison des deplacements-----------------------------------------
  •  
  • icomp = 1 ;
  • ino = 0 ;
  • REPETER BUX nbno ;
  • ino = ino + 1 ;
  • SI (EGA uxu1.ino 0.0) ;
  • SI (NON (EGA uxu2.ino 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' uxu2.ino
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 10 ; errabs = uxu2.ino ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dux = ABS (uxu1.ino - uxu2.ino) ;
  • errrel = ABS (dux / uxu1.ino) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message 'uxu1.ino =' uxu1.ino 'uxu2.ino =' uxu2.ino ;
  • Message ;
  • FINSI ;
  • ************ierr = 11 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BUX ;
  •  
  • icomp = 2 ;
  • ino = 0 ;
  • REPETER BUY nbno ;
  • ino = ino + 1 ;
  • SI (EGA uyu1.ino 0.0) ;
  • SI (NON (EGA uyu2.ino 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' uyu2.ino
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 10 ; errabs = uyu2.ino ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • duy = ABS (uyu1.ino - uyu2.ino) ;
  • errrel = ABS (duy / uyu1.ino) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message 'uyu1.ino =' uyu1.ino 'uyu2.ino =' uyu2.ino ;
  • Message ;
  • FINSI ;
  • ************ierr = 11 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BUY ;
  •  
  • icomp = 3 ;
  • ino = 0 ;
  • REPETER BUZ nbno ;
  • ino = ino + 1 ;
  • SI (EGA uzu1.ino 0.0) ;
  • SI (NON (EGA uzu2.ino 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' uzu2.ino
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 10 ; errabs = uzu2.ino ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • duz = ABS (uzu1.ino - uzu2.ino) ;
  • errrel = ABS (duz / uzu1.ino) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message 'uzu1.ino =' uzu1.ino 'uzu2.ino =' uzu2.ino ;
  • Message ;
  • FINSI ;
  • ************ierr = 11 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BUZ ;
  •  
  • * Comparaison des deformations-----------------------------------------
  •  
  • icomp = 1 ;
  • igau = 0 ;
  • REPETER BEPXX nbgau ;
  • igau = igau + 1 ;
  • SI (EGA epxxu1.igau 0.0) ;
  • SI (NON (EGA epxxu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' epxxu2.igau
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 20 ; errabs = epxxu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • depxx = ABS (epxxu1.igau - epxxu2.igau) ;
  • errrel = ABS (depxx / epxxu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'epxxu1.igau =' epxxu1.igau
  • 'epxxu2.igau =' epxxu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEPXX ;
  •  
  • icomp = 2 ;
  • igau = 0 ;
  • REPETER BEPYY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA epyyu1.igau 0.0) ;
  • SI (NON (EGA epyyu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' epyyu2.igau
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 20 ; errabs = epyyu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • depyy = ABS (epyyu1.igau - epyyu2.igau) ;
  • errrel = ABS (depyy / epyyu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'epyyu1.igau =' epyyu1.igau
  • 'epyyu2.igau =' epyyu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEPYY ;
  •  
  • icomp = 3 ;
  • igau = 0 ;
  • REPETER BEPZZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA epzzu1.igau 0.0) ;
  • SI (NON (EGA epzzu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' epzzu2.igau
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 20 ; errabs = epzzu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • depzz = ABS (epzzu1.igau - epzzu2.igau) ;
  • errrel = ABS (depzz / epzzu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'epzzu1.igau =' epzzu1.igau
  • 'epzzu2.igau =' epzzu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEPZZ ;
  •  
  • icomp = 4 ;
  • igau = 0 ;
  • REPETER BGAXY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA gaxyu1.igau 0.0) ;
  • SI (NON (EGA gaxyu2.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' gaxyu2.igau
  • ************ 'Composante ' (EXTR LDEFO icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 20 ; errabs = gaxyu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgaxy = ABS (gaxyu1.igau - gaxyu2.igau) ;
  • errrel = ABS (dgaxy / gaxyu1.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LDEFO icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'gaxyu1.igau =' gaxyu1.igau
  • ************ 'gaxyu2.igau =' gaxyu2.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGAXY ;
  •  
  • icomp = 5 ;
  • igau = 0 ;
  • REPETER BGAXZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA gaxzu1.igau 0.0) ;
  • SI (NON (EGA gaxzu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' gaxzu2.igau
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 20 ; errabs = gaxzu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgaxz = ABS (gaxzu1.igau - gaxzu2.igau) ;
  • errrel = ABS (dgaxz / gaxzu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'gaxzu1.igau =' gaxzu1.igau
  • 'gaxzu2.igau =' gaxzu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGAXZ ;
  •  
  • icomp = 6 ;
  • igau = 0 ;
  • REPETER BGAYZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA gayzu1.igau 0.0) ;
  • SI (NON (EGA gayzu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' gayzu2.igau
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 20 ; errabs = gayzu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgayz = ABS (gayzu1.igau - gayzu2.igau) ;
  • errrel = ABS (dgayz / gayzu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'gayzu1.igau =' gayzu1.igau
  • 'gayzu2.igau =' gayzu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGAYZ ;
  •  
  • * Comparaison des deformations inelastiques----------------------------
  •  
  • icomp = 1 ;
  • igau = 0 ;
  • REPETER BEIXX nbgau ;
  • igau = igau + 1 ;
  • SI (EGA eixxu1.igau 0.0) ;
  • SI (NON (EGA eixxu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' eixxu2.igau
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 30 ; errabs = eixxu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • deixx = ABS (eixxu1.igau - eixxu2.igau) ;
  • errrel = ABS (deixx / eixxu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'eixxu1.igau =' eixxu1.igau
  • 'eixxu2.igau =' eixxu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEIXX ;
  •  
  • icomp = 2 ;
  • igau = 0 ;
  • REPETER BEIYY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA eiyyu1.igau 0.0) ;
  • SI (NON (EGA eiyyu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' eiyyu2.igau
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 30 ; errabs = eiyyu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • deiyy = ABS (eiyyu1.igau - eiyyu2.igau) ;
  • errrel = ABS (deiyy / eiyyu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'eiyyu1.igau =' eiyyu1.igau
  • 'eiyyu2.igau =' eiyyu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEIYY ;
  •  
  • icomp = 3 ;
  • igau = 0 ;
  • REPETER BEIZZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA eizzu1.igau 0.0) ;
  • SI (NON (EGA eizzu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' eizzu2.igau
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 30 ; errabs = eizzu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • deizz = ABS (eizzu1.igau - eizzu2.igau) ;
  • errrel = ABS (deizz / eizzu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'eizzu1.igau =' eizzu1.igau
  • 'eizzu2.igau =' eizzu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEIZZ ;
  •  
  • icomp = 4 ;
  • igau = 0 ;
  • REPETER BGIXY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA gixyu1.igau 0.0) ;
  • SI (NON (EGA gixyu2.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' gixyu2.igau
  • ************ 'Composante ' (EXTR LDEIN icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 30 ; errabs = gixyu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgixy = ABS (gixyu1.igau - gixyu2.igau) ;
  • errrel = ABS (dgixy / gixyu1.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LDEIN icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'gixyu1.igau =' gixyu1.igau
  • ************ 'gixyu2.igau =' gixyu2.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGIXY ;
  •  
  • icomp = 5 ;
  • igau = 0 ;
  • REPETER BGIXZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA gixzu1.igau 0.0) ;
  • SI (NON (EGA gixzu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' gixzu2.igau
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 30 ; errabs = gixzu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgixz = ABS (gixzu1.igau - gixzu2.igau) ;
  • errrel = ABS (dgixz / gixzu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'gixzu1.igau =' gixzu1.igau
  • 'gixzu2.igau =' gixzu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGIXZ ;
  •  
  • icomp = 6 ;
  • igau = 0 ;
  • REPETER BGIYZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA giyzu1.igau 0.0) ;
  • SI (NON (EGA giyzu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' giyzu2.igau
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 30 ; errabs = giyzu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgiyz = ABS (giyzu1.igau - giyzu2.igau) ;
  • errrel = ABS (dgiyz / giyzu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'giyzu1.igau =' giyzu1.igau
  • 'giyzu2.igau =' giyzu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGIYZ ;
  •  
  • * Comparaison des contraintes------------------------------------------
  •  
  • icomp = 1 ;
  • igau = 0 ;
  • REPETER BSMXX nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smxxu1.igau 0.0) ;
  • SI (NON (EGA smxxu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' smxxu2.igau
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 40 ; errabs = smxxu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmxx = ABS (smxxu1.igau - smxxu2.igau) ;
  • errrel = ABS (dsmxx / smxxu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'smxxu1.igau =' smxxu1.igau
  • 'smxxu2.igau =' smxxu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMXX ;
  •  
  • icomp = 2 ;
  • igau = 0 ;
  • REPETER BSMYY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smyyu1.igau 0.0) ;
  • SI (NON (EGA smyyu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' smyyu2.igau
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 40 ; errabs = smyyu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmyy = ABS (smyyu1.igau - smyyu2.igau) ;
  • errrel = ABS (dsmyy / smyyu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'smyyu1.igau =' smyyu1.igau
  • 'smyyu2.igau =' smyyu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMYY ;
  •  
  • icomp = 3 ;
  • igau = 0 ;
  • REPETER BSMZZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smzzu1.igau 0.0) ;
  • SI (NON (EGA smzzu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' smzzu2.igau
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 40 ; errabs = smzzu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmzz = ABS (smzzu1.igau - smzzu2.igau) ;
  • errrel = ABS (dsmzz / smzzu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'smzzu1.igau =' smzzu1.igau
  • 'smzzu2.igau =' smzzu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMZZ ;
  •  
  • icomp = 4 ;
  • igau = 0 ;
  • REPETER BSMXY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smxyu1.igau 0.0) ;
  • SI (NON (EGA smxyu2.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' smxyu2.igau
  • ************ 'Composante ' (EXTR LCONT icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 40 ; errabs = smxyu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmxy = ABS (smxyu1.igau - smxyu2.igau) ;
  • errrel = ABS (dsmxy / smxyu1.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LCONT icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'smxyu1.igau =' smxyu1.igau
  • ************ 'smxyu2.igau =' smxyu2.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMXY ;
  •  
  • icomp = 5 ;
  • igau = 0 ;
  • REPETER BSMXZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smxzu1.igau 0.0) ;
  • SI (NON (EGA smxzu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' smxzu2.igau
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 40 ; errabs = smxzu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmxz = ABS (smxzu1.igau - smxzu2.igau) ;
  • errrel = ABS (dsmxz / smxzu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'smxzu1.igau =' smxzu1.igau
  • 'smxzu2.igau =' smxzu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMXZ ;
  •  
  • icomp = 6 ;
  • igau = 0 ;
  • REPETER BSMYZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smyzu1.igau 0.0) ;
  • SI (NON (EGA smyzu2.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' smyzu2.igau
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 40 ; errabs = smyzu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmyz = ABS (smyzu1.igau - smyzu2.igau) ;
  • errrel = ABS (dsmyz / smyzu1.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'smyzu1.igau =' smyzu1.igau
  • 'smyzu2.igau =' smyzu2.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMYZ ;
  •  
  • * Comparaison des variables internes-----------------------------------
  •  
  • icomp = 1 ;
  • igau = 0 ;
  • REPETER BEC0 nbgau ;
  • igau = igau + 1 ;
  • SI (EGA ec0u1.igau 0.0) ;
  • SI (NON (EGA ec0u2.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' ec0u2.igau
  • ************ 'Composante ' (EXTR LVARI icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 50 ; errabs = ec0u2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dec0 = ABS (ec0u1.igau - ec0u2.igau) ;
  • errrel = ABS (dec0 / ec0u1.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LVARI icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'ec0u1.igau =' ec0u1.igau
  • ************ 'ec0u2.igau =' ec0u2.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEC0 ;
  •  
  • icomp = 2 ;
  • igau = 0 ;
  • REPETER BESW0 nbgau ;
  • igau = igau + 1 ;
  • SI (EGA esw0u1.igau 0.0) ;
  • SI (NON (EGA esw0u2.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' esw0u2.igau
  • ************ 'Composante ' (EXTR LVARI icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 50 ; errabs = esw0u2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • desw0 = ABS (esw0u1.igau - esw0u2.igau) ;
  • errrel = ABS (desw0 / esw0u1.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LVARI icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'esw0u1.igau =' esw0u1.igau
  • ************ 'esw0u2.igau =' esw0u2.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BESW0 ;
  •  
  • icomp = 3 ;
  • igau = 0 ;
  • REPETER BP nbgau ;
  • igau = igau + 1 ;
  • SI (EGA pu1.igau 0.0) ;
  • SI (NON (EGA pu2.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' pu2.igau
  • ************ 'Composante ' (EXTR LVARI icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 50 ; errabs = pu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dp = ABS (pu1.igau - pu2.igau) ;
  • errrel = ABS (dp / pu1.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LVARI icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'pu1.igau =' pu1.igau
  • ************ 'pu2.igau =' pu2.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BP ;
  •  
  • icomp = 4 ;
  • igau = 0 ;
  • REPETER BQTLD nbgau ;
  • igau = igau + 1 ;
  • SI (EGA qtldu1.igau 0.0) ;
  • SI (NON (EGA qtldu2.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' qtldu2.igau
  • ************ 'Composante ' (EXTR LVARI icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 50 ; errabs = qtldu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dqtld = ABS (qtldu1.igau - qtldu2.igau) ;
  • errrel = ABS (dqtld / qtldu1.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LVARI icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'qtldu1.igau =' qtldu1.igau
  • ************ 'qtldu2.igau =' qtldu2.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BQTLD ;
  •  
  • * Comparaison des variables internes supplementaires (modele no24)-----
  •  
  • icomp = 1 ;
  • igau = 0 ;
  • REPETER BPSUP nbgau ;
  • igau = igau + 1 ;
  • SI (EGA pu2.igau 0.0) ;
  • SI (NON (EGA psupu2.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' psupu2.igau
  • ************ 'Composante ' (EXTR LVSUP icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 60 ; errabs = psupu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dpsup = ABS (pu2.igau - psupu2.igau) ;
  • errrel = ABS (dpsup / pu2.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LVSUP icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'pu2.igau =' pu2.igau
  • ************ 'psupu2.igau =' psupu2.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 61 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BPSUP ;
  •  
  • icomp = 2 ;
  • igau = 0 ;
  • REPETER BQSUP nbgau ;
  • igau = igau + 1 ;
  • SI (EGA qtldu2.igau 0.0) ;
  • SI (NON (EGA qsupu2.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' qsupu2.igau
  • ************ 'Composante ' (EXTR LVSUP icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 60 ; errabs = qsupu2.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dqsup = ABS (qtldu2.igau - qsupu2.igau) ;
  • errrel = ABS (dqsup / qtldu2.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LVSUP icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'qtldu2.igau =' qtldu2.igau
  • ************ 'qsupu2.igau =' qsupu2.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 61 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BQSUP ;
  •  
  • FIN BECART ;
  •  
  • * FIN BOUCLE COMPARAISON DES RESULTATS________________________________
  •  
  • * MESSAGES D'ERREUR_____________________________________________________
  •  
  • Message ;
  • Message 'CPU avec modele externe 22 = ' CPUuti1 'centiemes de seconde' ;
  • Message 'CPU avec modele externe 24 = ' CPUuti2 'centiemes de seconde' ;
  • Message ;
  •  
  • ifin = netat-1 ;
  •  
  • gaxyu1ma = 0.0 ;
  • igau = 0 ;
  • REPETER BGAXYFU1 nbgau ;
  • igau = igau + 1 ;
  • gaxy = EXTR (tabuti1.'DEFORMATIONS'.ifin)
  • 'GAXY' 1 1 igau ;
  • SI (gaxy > gaxyu1ma) ; gaxyu1ma = gaxy ; FINSI ;
  • FIN BGAXYFU1 ;
  • gaxyu2ma = 0.0 ;
  • igau = 0 ;
  • REPETER BGAXYFU2 nbgau ;
  • igau = igau + 1 ;
  • gaxy = EXTR (tabuti2.'DEFORMATIONS'.ifin)
  • 'GAXY' 1 1 igau ;
  • SI (gaxy > gaxyu2ma) ; gaxyu2ma = gaxy ; FINSI ;
  • FIN BGAXYFU2 ;
  • Message 'Deformation totale finale (maxi sur points de Gauss)' ;
  • Message 'GAXY max avec modele externe 22 = ' gaxyu1ma ;
  • Message 'GAXY max avec modele externe 24 = ' gaxyu2ma ;
  • Message ;
  •  
  • gixyu1ma = 0.0 ;
  • igau = 0 ;
  • REPETER BGIXYFU1 nbgau ;
  • igau = igau + 1 ;
  • gixy = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ifin)
  • 'GIXY' 1 1 igau ;
  • SI (gixy > gixyu1ma) ; gixyu1ma = gixy ; FINSI ;
  • FIN BGIXYFU1 ;
  • gixyu2ma = 0.0 ;
  • igau = 0 ;
  • REPETER BGIXYFU2 nbgau ;
  • igau = igau + 1 ;
  • gixy = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ifin)
  • 'GIXY' 1 1 igau ;
  • SI (gixy > gixyu2ma) ; gixyu2ma = gixy ; FINSI ;
  • FIN BGIXYFU2 ;
  • Message 'Deformation de fluage finale (maxi sur points de Gauss)' ;
  • Message 'GIXY max avec modele externe 22 = ' gixyu1ma ;
  • Message 'GIXY max avec modele externe 24 = ' gixyu2ma ;
  • Message ;
  •  
  • smxyu1ma = 0.0 ;
  • igau = 0 ;
  • REPETER BSMXYFU1 nbgau ;
  • igau = igau + 1 ;
  • smxy = EXTR (tabuti1.'CONTRAINTES'.ifin)
  • 'SMXY' 1 1 igau ;
  • SI (smxy > smxyu1ma) ; smxyu1ma = smxy ; FINSI ;
  • FIN BSMXYFU1 ;
  • smxyu2ma = 0.0 ;
  • igau = 0 ;
  • REPETER BSMXYFU2 nbgau ;
  • igau = igau + 1 ;
  • smxy = EXTR (tabuti2.'CONTRAINTES'.ifin)
  • 'SMXY' 1 1 igau ;
  • SI (smxy > smxyu2ma) ; smxyu2ma = smxy ; FINSI ;
  • FIN BSMXYFU2 ;
  • Message 'Contrainte finale (maxi sur points de Gauss)' ;
  • Message 'SMXY max avec modele externe 22 = ' smxyu1ma ;
  • Message 'SMXY max avec modele externe 24 = ' smxyu2ma ;
  • Message ;
  •  
  • SI (ierr > 0) ;
  •  
  • SI (EGA ierr 1) ;
  • Message 'Instants de calcul non coincidants' ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 10) ;
  • Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 11) ;
  • Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 20) ;
  • Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 21) ;
  • Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 30) ;
  • Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 31) ;
  • Message 'Erreur relative' errrel 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 40) ;
  • Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 41) ;
  • Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 50) ;
  • Message 'Erreur absolue' errabs 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 51) ;
  • Message 'Erreur relative' errrel 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 60) ;
  • Message 'Erreur absolue' errabs 'Composante ' (EXTR LVSUP icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 61) ;
  • Message 'Erreur relative' errrel 'Composante ' (EXTR LVSUP icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SINON ;
  •  
  • Message 'Resultats concordants entre les deux modeles' ;
  • Message 'Aucune erreur relative superieure a ' critere
  • ' pour les composantes' ;
  • Message (EXTR LDEFO idefo) ' ' (EXTR LDEIN idein) ' '
  • (EXTR LCONT icont) ' '
  • (EXTR LVARI ivar1) ' ' (EXTR LVARI ivar2) ' '
  • (EXTR LVARI ivar3) ' ' (EXTR LVARI ivar4) ;
  • Message ;
  • Message 'Modele externe 24 : aucune erreur relative superieure a '
  • critere ;
  • Message 'entre les composantes '
  • (EXTR LVARI ivar3) ' ' (EXTR LVARI ivar4) ' et '
  • (EXTR LVSUP ivas1) ' ' (EXTR LVSUP ivas2) ;
  • ERRE 0 ;
  •  
  • FINSI ;
  •  
  • FIN ;
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  • © Cast3M 2003 - Tous droits réservés.
    Mentions légales