Télécharger umat05.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : umat05.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 'NON_LINEAIRE' 'UTILISATEUR', integrateur specifique UMAT
  13. *
  14. * Test No5 : Modele PLASTIQUE ISOTROPE
  15. * Reprise du test isotro_cis.dgibi
  16. *
  17. *=======================================================================
  18.  
  19.  
  20. opti echo 0 ;
  21. * ajout de option epsilon lineaire pour la precision des test!
  22. OPTION epsilon lineaire;
  23.  
  24.  
  25. *=======================================================================
  26. *
  27. * PROCEDURES DEPENDANTES
  28. *
  29. *=======================================================================
  30.  
  31.  
  32.  
  33.  
  34. *=======================================================================
  35. *
  36. * JEU DE DONNEES PRINCIPAL
  37. *
  38. *=======================================================================
  39.  
  40. *=======================================================================
  41. * Resolution avec le modele 'PLASTIQUE' 'ISOTROPE' de CAST3M
  42. * Reprise du jeu de donnees isotro_cis.dgibi
  43. * Element CUB8
  44.  
  45. **************************************************
  46. * *
  47. * TEST DE VALIDATION D'UNE LOI DE COMPORTEMENT *
  48. * -------------------------------------------- *
  49. * DE MATERIAU *
  50. * ----------- *
  51. * *
  52. * LOI DE COMPORTEMENT ELASTO-PLASTIQUE *
  53. * ISOTROPE *
  54. * *
  55. * COMPARAISON DE CALCULS SUR DES ELEMENTS: *
  56. * - COQUE EPAISSE ( MFR=5 ) *
  57. * - COQUE MINCE ( MFR=3 ) *
  58. * - MASSIF 3D ( MFR=1 ) *
  59. * *
  60. * STRUCTURE : UNE PLAQUE CARREE *
  61. * - COTE : .1 M *
  62. * - EPAISSEUR : 5.E-4 M *
  63. * *
  64. * CHARGEMENT : DEPLACEMENTS IMPOSES *
  65. * CISAILLEMENT PUR *
  66. * *
  67. **************************************************
  68. opti echo 0 dime 3 elem cub8 ;
  69. *
  70. *
  71. * Maillage
  72. *
  73. rayon0 = .05 ;
  74. epais0 = .0005 ;
  75. l1 = epais0 + rayon0 ;
  76. l2 = .1 ;
  77. n1 = 1 ;
  78. n2 = 1 ;
  79. p1 = rayon0 0. 0. ;
  80. p2 = l1 0. 0. ;
  81. p3 = l1 l2 0. ;
  82. p4 = rayon0 l2 0. ;
  83. p11 = rayon0 0. l2 ;
  84. p21 = l1 0. l2 ;
  85. p31 = l1 l2 l2 ;
  86. p41 = rayon0 l2 l2 ;
  87. d1 = d p1 n1 p2 ;
  88. d2 = d p2 n2 p3 ;
  89. d3 = d p3 n1 p4 ;
  90. d4 = d p4 n2 p1 ;
  91. d11 = d p11 n1 p21 ;
  92. d21 = d p21 n2 p31 ;
  93. d31 = d p31 n1 p41 ;
  94. d41 = d p41 n2 p11 ;
  95. *p01 = ( rayon0 + ( epais0 / 2. ) ) 0. l2 ;
  96. *p02 = ( rayon0 + ( epais0 / 2. ) ) l2 l2 ;
  97. sur1 = 'DALL' d1 d2 d3 d4 ;
  98. sur2 = 'DALL' d11 d21 d31 d41 ;
  99. vol1 = sur1 'VOLU' n2 sur2 ;
  100. *titr ' maillage ' ;
  101. *trace vol1 ;
  102. *
  103. * Deplacement maximal
  104. *
  105. depmax = .01 ;
  106. *
  107. * Courbe d'ecrouissage
  108. *
  109. beta0 = 1250. ;
  110. sig0 = 410.E6 ;
  111. xm = .5 ;
  112. youn0 = 7.34E10 ;
  113. epsp0 = 'PROG' 0. ;
  114. deb0 = 'PROG' 0. ;
  115. epsmax0 = depmax / l2 ;
  116. pas0 = epsmax0 / 50 ;
  117. epsp1 = 'PROG' 0. 'PAS' pas0 epsmax0 ;
  118. dim00 = 'DIME' epsp1 ;
  119. un = 'PROG' dim00 '*' 1. ;
  120. Y0 = sig0 ;
  121. epsp2 = beta0 * epsp1 ;
  122. epsp3 = epsp2 + un ;
  123. epsp31 = log epsp3 ;
  124. epsp32 = epsp31 * xm ;
  125. epsp4 = exp epsp32 ;
  126. Y1 = Y0 * epsp4 ;
  127. *
  128. epsp1 = epsp0 'ET' epsp1 ;
  129. Y1 = deb0 'ET' Y1 ;
  130. epsp2 = epsp1 + ( Y1 / youn0 ) ;
  131. trac00 = 'EVOL' 'MANU' 'defo' epsp1 'sigm' Y1 ;
  132. *
  133. * Modele de calcul
  134. *
  135. trac0 = 'EVOL' 'MANU' 'defo' epsp2 'sigm' Y1 ;
  136. lsm1 = y1 enle 1 ;
  137. lep1 = (epsp2 enle 1) - (lsm1 / youn0) ;
  138. evec0 = evol vert manu eps lep1 sig lsm1 ;
  139. *dess (trac0 et evec0) titr ' Courbes de traction et d ecrouissage (vert)' ;
  140. mod0 = 'MODE' vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  141. 'PLASTIQUE' 'ISOTROPE' ;
  142. mat0 = 'MATE' mod0 'YOUN' youn0 'NU' 0.33
  143. 'ECRO' evec0 ;
  144. *
  145. * Conditions aux limites
  146. *
  147. cl1 = 'BLOQ' sur1 'UX' 'RY' ;
  148. cl2 = 'BLOQ' sur2 'UX' 'RY' ;
  149. cl3 = 'BLOQ' vol1 'UY' 'UZ' ;
  150. cl0 = cl1 'ET' cl2 'ET' cl3 ;
  151. *
  152. * Chargement
  153. *
  154. chp1 = 'COOR' 3 vol1 ;
  155. chp1 = chp1 * depmax / l2 ;
  156. chp1 = 'NOMC' chp1 'UY' ;
  157. chp2 = 'MANU' 'CHPO' vol1 1 'UZ' 0. ;
  158. *trace chp1 vol1 ;
  159. chp0 = chp1 + chp2 ;
  160. dep1 = 'DEPI' cl3 chp0 ;
  161. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y ( 'PROG' 0. 1. ) ;
  162. cha0 = 'CHAR' 'DIMP' dep1 ev0 ;
  163. *vec1 = 'VECT' dep1 1.E-5 'FX' 'FY' 'FZ' 'ROUG' ;
  164. *titr ' chargement ' ;
  165. *trace vec1 vol1 ;
  166. *
  167. * Temps du calcul
  168. *
  169. dt0 = 50. ;
  170. tfin0 = 1000. ;
  171. *
  172. *
  173. * Resolution par PASAPAS
  174. *---------------------------------------
  175. *
  176. ta1 = 'TABLE' ;
  177. ta1 .'MODELE' = mod0 ;
  178. ta1 .'CARACTERISTIQUES' = mat0 ;
  179. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  180. ta1 .'CHARGEMENT' = cha0 ;
  181. *ta1 .'TEMPS_CALCULES' = ( ( 'PROG' 0. 'PAS' dt0 80. ) 'ET'
  182. *( 'PROG' 85. 'PAS' 5. 90. ) 'ET'
  183. *( 'PROG' 92.5 'PAS' 2.5 tfin0 ) ) ;
  184. ta1 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt0 tfin0 ;
  185. *
  186. * MODIF
  187. * MODIF_________________________________________________________________
  188. TEMPS 'ZERO' ;
  189. * MODIF_________________________________________________________________
  190. * MODIF
  191. TMASAU=table;
  192. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  193. TMASAU .'DEFTO'=VRAI;
  194. TMASAU .'DEFIN'=VRAI;
  195. PASAPAS ta1 ;
  196. * MODIF
  197. * MODIF_________________________________________________________________
  198. TABTPS = TEMP 'NOEC';
  199. CPUint = TABTPS.'TEMPS_CPU'.'INITIAL';
  200. * MODIF_________________________________________________________________
  201. * MODIF
  202.  
  203.  
  204. *=======================================================================
  205. * Resolution avec le modele 'NON_LINEAIRE' 'UTILISATEUR'
  206.  
  207. * Definition du modele et affectation des proprietes de materiau--------
  208.  
  209. LCM05_0 = MOTS 'YOUN' 'NU ' ;
  210.  
  211. LCM05_1 = MOTS 'E001' 'S001' 'E002' 'S002' 'E003' 'S003' 'E004' 'S004'
  212. 'E005' 'S005' 'E006' 'S006' 'E007' 'S007' 'E008' 'S008'
  213. 'E009' 'S009' 'E010' 'S010' ;
  214.  
  215. LCM05_2 = MOTS 'E011' 'S011' 'E012' 'S012' 'E013' 'S013' 'E014' 'S014'
  216. 'E015' 'S015' 'E016' 'S016' 'E017' 'S017' 'E018' 'S018'
  217. 'E019' 'S019' 'E020' 'S020' ;
  218.  
  219. LCM05_3 = MOTS 'E021' 'S021' 'E022' 'S022' 'E023' 'S023' 'E024' 'S024'
  220. 'E025' 'S025' 'E026' 'S026' 'E027' 'S027' 'E028' 'S028'
  221. 'E029' 'S029' 'E030' 'S030' ;
  222.  
  223. LCM05_4 = MOTS 'E031' 'S031' 'E032' 'S032' 'E033' 'S033' 'E034' 'S034'
  224. 'E035' 'S035' 'E036' 'S036' 'E037' 'S037' 'E038' 'S038'
  225. 'E039' 'S039' 'E040' 'S040' ;
  226.  
  227. LCM05_5 = MOTS 'E041' 'S041' 'E042' 'S042' 'E043' 'S043' 'E044' 'S044'
  228. 'E045' 'S045' 'E046' 'S046' 'E047' 'S047' 'E048' 'S048'
  229. 'E049' 'S049' 'E050' 'S050' 'E051' 'S051' 'E052' 'S052' ;
  230.  
  231. LCMAT05 = LCM05_0 et LCM05_1 et LCM05_2 et LCM05_3 et
  232. LCM05_4 et LCM05_5 ;
  233.  
  234. LCVAR05 = MOTS 'EPSE' 'ENXX' 'ENYY' 'ENZZ' 'GNXY' 'GNXZ' 'GNYZ' ;
  235.  
  236. mo_util = MODE vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  237. 'NON_LINEAIRE' 'UTILISATEUR' 'NUME_LOI' 5
  238. 'C_MATERIAU' LCMAT05 'C_VARINTER' LCVAR05 ;
  239.  
  240. ma_u0 = MATE mo_util 'YOUN' youn0 'NU ' 0.33 ;
  241.  
  242. ma_u11 = MATE mo_util 'E001' (EXTR epsp2 1) 'S001' (EXTR Y1 1)
  243. 'E002' (EXTR epsp2 2) 'S002' (EXTR Y1 2)
  244. 'E003' (EXTR epsp2 3) 'S003' (EXTR Y1 3)
  245. 'E004' (EXTR epsp2 4) 'S004' (EXTR Y1 4)
  246. 'E005' (EXTR epsp2 5) 'S005' (EXTR Y1 5) ;
  247.  
  248. ma_u12 = MATE mo_util 'E006' (EXTR epsp2 6) 'S006' (EXTR Y1 6)
  249. 'E007' (EXTR epsp2 7) 'S007' (EXTR Y1 7)
  250. 'E008' (EXTR epsp2 8) 'S008' (EXTR Y1 8)
  251. 'E009' (EXTR epsp2 9) 'S009' (EXTR Y1 9)
  252. 'E010' (EXTR epsp2 10) 'S010' (EXTR Y1 10) ;
  253.  
  254. ma_u21 = MATE mo_util 'E011' (EXTR epsp2 11) 'S011' (EXTR Y1 11)
  255. 'E012' (EXTR epsp2 12) 'S012' (EXTR Y1 12)
  256. 'E013' (EXTR epsp2 13) 'S013' (EXTR Y1 13)
  257. 'E014' (EXTR epsp2 14) 'S014' (EXTR Y1 14)
  258. 'E015' (EXTR epsp2 15) 'S015' (EXTR Y1 15) ;
  259.  
  260. ma_u22 = MATE mo_util 'E016' (EXTR epsp2 16) 'S016' (EXTR Y1 16)
  261. 'E017' (EXTR epsp2 17) 'S017' (EXTR Y1 17)
  262. 'E018' (EXTR epsp2 18) 'S018' (EXTR Y1 18)
  263. 'E019' (EXTR epsp2 19) 'S019' (EXTR Y1 19)
  264. 'E020' (EXTR epsp2 20) 'S020' (EXTR Y1 20) ;
  265.  
  266. ma_u31 = MATE mo_util 'E021' (EXTR epsp2 21) 'S021' (EXTR Y1 21)
  267. 'E022' (EXTR epsp2 22) 'S022' (EXTR Y1 22)
  268. 'E023' (EXTR epsp2 23) 'S023' (EXTR Y1 23)
  269. 'E024' (EXTR epsp2 24) 'S024' (EXTR Y1 24)
  270. 'E025' (EXTR epsp2 25) 'S025' (EXTR Y1 25) ;
  271.  
  272. ma_u32 = MATE mo_util 'E026' (EXTR epsp2 26) 'S026' (EXTR Y1 26)
  273. 'E027' (EXTR epsp2 27) 'S027' (EXTR Y1 27)
  274. 'E028' (EXTR epsp2 28) 'S028' (EXTR Y1 28)
  275. 'E029' (EXTR epsp2 29) 'S029' (EXTR Y1 29)
  276. 'E030' (EXTR epsp2 30) 'S030' (EXTR Y1 30) ;
  277.  
  278. ma_u41 = MATE mo_util 'E031' (EXTR epsp2 31) 'S031' (EXTR Y1 31)
  279. 'E032' (EXTR epsp2 32) 'S032' (EXTR Y1 32)
  280. 'E033' (EXTR epsp2 33) 'S033' (EXTR Y1 33)
  281. 'E034' (EXTR epsp2 34) 'S034' (EXTR Y1 34)
  282. 'E035' (EXTR epsp2 35) 'S035' (EXTR Y1 35) ;
  283.  
  284. ma_u42 = MATE mo_util 'E036' (EXTR epsp2 36) 'S036' (EXTR Y1 36)
  285. 'E037' (EXTR epsp2 37) 'S037' (EXTR Y1 37)
  286. 'E038' (EXTR epsp2 38) 'S038' (EXTR Y1 38)
  287. 'E039' (EXTR epsp2 39) 'S039' (EXTR Y1 39)
  288. 'E040' (EXTR epsp2 40) 'S040' (EXTR Y1 40) ;
  289.  
  290. ma_u51 = MATE mo_util 'E041' (EXTR epsp2 41) 'S041' (EXTR Y1 41)
  291. 'E042' (EXTR epsp2 42) 'S042' (EXTR Y1 42)
  292. 'E043' (EXTR epsp2 43) 'S043' (EXTR Y1 43)
  293. 'E044' (EXTR epsp2 44) 'S044' (EXTR Y1 44)
  294. 'E045' (EXTR epsp2 45) 'S045' (EXTR Y1 45) ;
  295.  
  296. ma_u52 = MATE mo_util 'E046' (EXTR epsp2 46) 'S046' (EXTR Y1 46)
  297. 'E047' (EXTR epsp2 47) 'S047' (EXTR Y1 47)
  298. 'E048' (EXTR epsp2 48) 'S048' (EXTR Y1 48)
  299. 'E049' (EXTR epsp2 49) 'S049' (EXTR Y1 49)
  300. 'E050' (EXTR epsp2 50) 'S050' (EXTR Y1 50) ;
  301.  
  302. ma_u61 = MATE mo_util 'E051' (EXTR epsp2 51) 'S051' (EXTR Y1 51)
  303. 'E052' (EXTR epsp2 52) 'S052' (EXTR Y1 52) ;
  304.  
  305. ma_util = ma_u0 et ma_u11 et ma_u12 et ma_u21 et ma_u22 et
  306. ma_u31 et ma_u32 et ma_u41 et ma_u42 et
  307. ma_u51 et ma_u52 et ma_u61 ;
  308.  
  309. * Resolution par PASAPAS------------------------------------------------
  310.  
  311. TABU = TABLE ;
  312. TABU.'MODELE' = mo_util ;
  313. TABU.'CARACTERISTIQUES' = ma_util ;
  314. TABU.'BLOCAGES_MECANIQUES' = cl0 ;
  315. TABU.'CHARGEMENT' = cha0 ;
  316. TABU.'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt0 tfin0 ;
  317. TMASAU=table;
  318. tabu . 'MES_SAUVEGARDES'=TMASAU;
  319. TMASAU .'DEFTO'=VRAI;
  320. TMASAU .'DEFIN'=VRAI;
  321. TEMPS 'ZERO' ;
  322. PASAPAS TABU ;
  323. TABTPS = TEMP 'NOEC';
  324. CPUext = TABTPS.'TEMPS_CPU'.'INITIAL';
  325.  
  326.  
  327. *=======================================================================
  328. * Comparaison des resultats obtenus par les deux modeles
  329. *
  330. * N.B. Comme on simule un essai de cisaillement YZ :
  331. * - Critere strict (errrel < 1.d-12) sur GAYZ, GIYZ et SMYZ
  332. * - Pour les autres composantes on note les depassements du critere
  333. * si impression demandee
  334.  
  335. ***opti echo 1 ;
  336.  
  337. uxc = TABLE ; uyc = TABLE ; uzc = TABLE ;
  338. epxxc = TABLE ; epyyc = TABLE ; epzzc = TABLE ;
  339. gaxyc = TABLE ; gaxzc = TABLE ; gayzc = TABLE ;
  340. eixxc = TABLE ; eiyyc = TABLE ; eizzc = TABLE ;
  341. gixyc = TABLE ; gixzc = TABLE ; giyzc = TABLE ;
  342. smxxc = TABLE ; smyyc = TABLE ; smzzc = TABLE ;
  343. smxyc = TABLE ; smxzc = TABLE ; smyzc = TABLE ;
  344. epsec = TABLE ;
  345.  
  346. uxu = TABLE ; uyu = TABLE ; uzu = TABLE ;
  347. epxxu = TABLE ; epyyu = TABLE ; epzzu = TABLE ;
  348. gaxyu = TABLE ; gaxzu = TABLE ; gayzu = TABLE ;
  349. enxxu = TABLE ; enyyu = TABLE ; enzzu = TABLE ;
  350. gnxyu = TABLE ; gnxzu = TABLE ; gnyzu = TABLE ;
  351. smxxu = TABLE ; smyyu = TABLE ; smzzu = TABLE ;
  352. smxyu = TABLE ; smxzu = TABLE ; smyzu = TABLE ;
  353. epseu = TABLE ;
  354.  
  355. LDEPL = MOTS 'UX ' 'UY ' 'UZ ' ;
  356. LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  357. LDEIN = MOTS 'EIXX' 'EIYY' 'EIZZ' 'GIXY' 'GIXZ' 'GIYZ' ;
  358. LDENN = MOTS 'ENXX' 'ENYY' 'ENZZ' 'GNXY' 'GNXZ' 'GNYZ' ;
  359. LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  360. LVARI = MOTS 'EPSE' ;
  361.  
  362. idefo = 6 ;
  363. idein = 6 ;
  364. icont = 6 ;
  365.  
  366. ierr = 0 ;
  367. critere = 1.d-12 ;
  368. IMESS = 0 ;
  369. ***IMESS = 1 ;
  370.  
  371. nbno = 8 ;
  372. nbgau = 8 ;
  373.  
  374. netat = DIME (ta1.'TEMPS') ;
  375. ietat = -1 ;
  376.  
  377. * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  378.  
  379. REPETER BECART netat ;
  380.  
  381. ietat = ietat + 1 ;
  382.  
  383. temps1 = ta1.'TEMPS'.ietat ;
  384. temps2 = TABU.'TEMPS'.ietat ;
  385. SI ((ABS (temps1-temps2)) > 1.0E-9) ;
  386. ierr = 1 ;
  387. QUITTER BECART ;
  388. FINSI ;
  389.  
  390. * Extraction des resultats de calcul : Modele interne CAST3M-----------
  391.  
  392. * Deplacements
  393. uxc.1 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p1 ;
  394. uxc.2 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p2 ;
  395. uxc.3 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p3 ;
  396. uxc.4 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p4 ;
  397. uxc.5 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p11 ;
  398. uxc.6 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p21 ;
  399. uxc.7 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p31 ;
  400. uxc.8 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p41 ;
  401.  
  402. uyc.1 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p1 ;
  403. uyc.2 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p2 ;
  404. uyc.3 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p3 ;
  405. uyc.4 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p4 ;
  406. uyc.5 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p11 ;
  407. uyc.6 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p21 ;
  408. uyc.7 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p31 ;
  409. uyc.8 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p41 ;
  410.  
  411. uzc.1 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p1 ;
  412. uzc.2 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p2 ;
  413. uzc.3 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p3 ;
  414. uzc.4 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p4 ;
  415. uzc.5 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p11 ;
  416. uzc.6 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p21 ;
  417. uzc.7 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p31 ;
  418. uzc.8 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p41 ;
  419.  
  420. * Deformations totales
  421. igau = 0 ;
  422. REPETER BEPXXC nbgau ;
  423. igau = igau + 1 ;
  424. epxxc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  425. FIN BEPXXC ;
  426.  
  427. igau = 0 ;
  428. REPETER BEPYYC nbgau ;
  429. igau = igau + 1 ;
  430. epyyc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  431. FIN BEPYYC ;
  432.  
  433. igau = 0 ;
  434. REPETER BEPZZC nbgau ;
  435. igau = igau + 1 ;
  436. epzzc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  437. FIN BEPZZC ;
  438.  
  439. igau = 0 ;
  440. REPETER BGAXYC nbgau ;
  441. igau = igau + 1 ;
  442. gaxyc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  443. FIN BGAXYC ;
  444.  
  445. igau = 0 ;
  446. REPETER BGAXZC nbgau ;
  447. igau = igau + 1 ;
  448. gaxzc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  449. FIN BGAXZC ;
  450.  
  451. igau = 0 ;
  452. REPETER BGAYZC nbgau ;
  453. igau = igau + 1 ;
  454. gayzc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  455. FIN BGAYZC ;
  456.  
  457. * Deformations inelastiques
  458. igau = 0 ;
  459. REPETER BEIXXC nbgau ;
  460. igau = igau + 1 ;
  461. eixxc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  462. 'EIXX' 1 1 igau ;
  463. FIN BEIXXC ;
  464.  
  465. igau = 0 ;
  466. REPETER BEIYYC nbgau ;
  467. igau = igau + 1 ;
  468. eiyyc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  469. 'EIYY' 1 1 igau ;
  470. FIN BEIYYC ;
  471.  
  472. igau = 0 ;
  473. REPETER BEIZZC nbgau ;
  474. igau = igau + 1 ;
  475. eizzc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  476. 'EIZZ' 1 1 igau ;
  477. FIN BEIZZC ;
  478.  
  479. igau = 0 ;
  480. REPETER BGIXYC nbgau ;
  481. igau = igau + 1 ;
  482. gixyc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  483. 'GIXY' 1 1 igau ;
  484. FIN BGIXYC ;
  485.  
  486. igau = 0 ;
  487. REPETER BGIXZC nbgau ;
  488. igau = igau + 1 ;
  489. gixzc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  490. 'GIXZ' 1 1 igau ;
  491. FIN BGIXZC ;
  492.  
  493. igau = 0 ;
  494. REPETER BGIYZC nbgau ;
  495. igau = igau + 1 ;
  496. giyzc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  497. 'GIYZ' 1 1 igau ;
  498. FIN BGIYZC ;
  499.  
  500. * Contraintes
  501. igau = 0 ;
  502. REPETER BSMXXC nbgau ;
  503. igau = igau + 1 ;
  504. smxxc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  505. FIN BSMXXC ;
  506.  
  507. igau = 0 ;
  508. REPETER BSMYYC nbgau ;
  509. igau = igau + 1 ;
  510. smyyc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  511. FIN BSMYYC ;
  512.  
  513. igau = 0 ;
  514. REPETER BSMZZC nbgau ;
  515. igau = igau + 1 ;
  516. smzzc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  517. FIN BSMZZC ;
  518.  
  519. igau = 0 ;
  520. REPETER BSMXYC nbgau ;
  521. igau = igau + 1 ;
  522. smxyc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  523. FIN BSMXYC ;
  524.  
  525. igau = 0 ;
  526. REPETER BSMXZC nbgau ;
  527. igau = igau + 1 ;
  528. smxzc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  529. FIN BSMXZC ;
  530.  
  531. igau = 0 ;
  532. REPETER BSMYZC nbgau ;
  533. igau = igau + 1 ;
  534. smyzc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  535. FIN BSMYZC ;
  536.  
  537. * Variables internes
  538. igau = 0 ;
  539. REPETER BEPSEC nbgau ;
  540. igau = igau + 1 ;
  541. epsec.igau = EXTR (ta1.'VARIABLES_INTERNES'.ietat)
  542. 'EPSE' 1 1 igau ;
  543. FIN BEPSEC ;
  544.  
  545. * Extraction des resultats de calcul : Modele externe utilisateur------
  546.  
  547. uxu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p1 ;
  548. uxu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p2 ;
  549. uxu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p3 ;
  550. uxu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p4 ;
  551. uxu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p11 ;
  552. uxu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p21 ;
  553. uxu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p31 ;
  554. uxu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p41 ;
  555.  
  556. uyu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p1 ;
  557. uyu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p2 ;
  558. uyu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p3 ;
  559. uyu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p4 ;
  560. uyu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p11 ;
  561. uyu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p21 ;
  562. uyu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p31 ;
  563. uyu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p41 ;
  564.  
  565. uzu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p1 ;
  566. uzu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p2 ;
  567. uzu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p3 ;
  568. uzu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p4 ;
  569. uzu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p11 ;
  570. uzu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p21 ;
  571. uzu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p31 ;
  572. uzu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p41 ;
  573.  
  574. * Deformations totales
  575. igau = 0 ;
  576. REPETER BEPXXU nbgau ;
  577. igau = igau + 1 ;
  578. epxxu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  579. FIN BEPXXU ;
  580.  
  581. igau = 0 ;
  582. REPETER BEPYYU nbgau ;
  583. igau = igau + 1 ;
  584. epyyu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  585. FIN BEPYYU ;
  586.  
  587. igau = 0 ;
  588. REPETER BEPZZU nbgau ;
  589. igau = igau + 1 ;
  590. epzzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  591. FIN BEPZZU ;
  592.  
  593. igau = 0 ;
  594. REPETER BGAXYU nbgau ;
  595. igau = igau + 1 ;
  596. gaxyu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  597. FIN BGAXYU ;
  598.  
  599. igau = 0 ;
  600. REPETER BGAXZU nbgau ;
  601. igau = igau + 1 ;
  602. gaxzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  603. FIN BGAXZU ;
  604.  
  605. igau = 0 ;
  606. REPETER BGAYZU nbgau ;
  607. igau = igau + 1 ;
  608. gayzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  609. FIN BGAYZU ;
  610.  
  611. * Deformations inelastiques (variables internes supplementaires)
  612. igau = 0 ;
  613. REPETER BENXXU nbgau ;
  614. igau = igau + 1 ;
  615. enxxu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  616. 'ENXX' 1 1 igau ;
  617. FIN BENXXU ;
  618.  
  619. igau = 0 ;
  620. REPETER BENYYU nbgau ;
  621. igau = igau + 1 ;
  622. enyyu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  623. 'ENYY' 1 1 igau ;
  624. FIN BENYYU ;
  625.  
  626. igau = 0 ;
  627. REPETER BENZZU nbgau ;
  628. igau = igau + 1 ;
  629. enzzu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  630. 'ENZZ' 1 1 igau ;
  631. FIN BENZZU ;
  632.  
  633. igau = 0 ;
  634. REPETER BGNXYU nbgau ;
  635. igau = igau + 1 ;
  636. gnxyu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  637. 'GNXY' 1 1 igau ;
  638. FIN BGNXYU ;
  639.  
  640. igau = 0 ;
  641. REPETER BGNXZU nbgau ;
  642. igau = igau + 1 ;
  643. gnxzu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  644. 'GNXZ' 1 1 igau ;
  645. FIN BGNXZU ;
  646.  
  647. igau = 0 ;
  648. REPETER BGNYZU nbgau ;
  649. igau = igau + 1 ;
  650. gnyzu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  651. 'GNYZ' 1 1 igau ;
  652. FIN BGNYZU ;
  653.  
  654. * Contraintes
  655. igau = 0 ;
  656. REPETER BSMXXU nbgau ;
  657. igau = igau + 1 ;
  658. smxxu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  659. FIN BSMXXU ;
  660.  
  661. igau = 0 ;
  662. REPETER BSMYYU nbgau ;
  663. igau = igau + 1 ;
  664. smyyu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  665. FIN BSMYYU ;
  666.  
  667. igau = 0 ;
  668. REPETER BSMZZU nbgau ;
  669. igau = igau + 1 ;
  670. smzzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  671. FIN BSMZZU ;
  672.  
  673. igau = 0 ;
  674. REPETER BSMXYU nbgau ;
  675. igau = igau + 1 ;
  676. smxyu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  677. FIN BSMXYU ;
  678.  
  679. igau = 0 ;
  680. REPETER BSMXZU nbgau ;
  681. igau = igau + 1 ;
  682. smxzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  683. FIN BSMXZU ;
  684.  
  685. igau = 0 ;
  686. REPETER BSMYZU nbgau ;
  687. igau = igau + 1 ;
  688. smyzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  689. FIN BSMYZU ;
  690.  
  691. * Variables internes
  692. igau = 0 ;
  693. REPETER BEPSEU nbgau ;
  694. igau = igau + 1 ;
  695. epseu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  696. 'EPSE' 1 1 igau ;
  697. FIN BEPSEU ;
  698.  
  699. * Comparaison des deplacements-----------------------------------------
  700.  
  701. icomp = 1 ;
  702. ino = 0 ;
  703. REPETER BUX nbno ;
  704. ino = ino + 1 ;
  705. SI (EGA uxc.ino 0.0) ;
  706. SI (NON (EGA uxu.ino 0.0)) ;
  707. SI (EGA IMESS 1) ;
  708. Message 'Erreur absolue' uxu.ino
  709. 'Composante ' (EXTR LDEPL icomp)
  710. ' Noeud' ino 'Etat' ietat ;
  711. Message ;
  712. FINSI ;
  713. ************ierr = 10 ; errabs = uxu.ino ; QUITTER BECART ;
  714. FINSI ;
  715. SINON ;
  716. dux = ABS (uxc.ino - uxu.ino) ;
  717. errrel = ABS (dux / uxc.ino) ;
  718. SI (errrel > critere) ;
  719. SI (EGA IMESS 1) ;
  720. Message 'Erreur relative' errrel
  721. 'Composante ' (EXTR LDEPL icomp)
  722. ' Noeud' ino 'Etat' ietat ;
  723. Message 'uxc.ino =' uxc.ino 'uxu.ino =' uxu.ino ;
  724. Message ;
  725. FINSI ;
  726. ************ierr = 11 ; QUITTER BECART ;
  727. FINSI ;
  728. FINSI ;
  729. FIN BUX ;
  730.  
  731. icomp = 2 ;
  732. ino = 0 ;
  733. REPETER BUY nbno ;
  734. ino = ino + 1 ;
  735. SI (EGA uyc.ino 0.0) ;
  736. SI (NON (EGA uyu.ino 0.0)) ;
  737. SI (EGA IMESS 1) ;
  738. Message 'Erreur absolue' uyu.ino
  739. 'Composante ' (EXTR LDEPL icomp)
  740. ' Noeud' ino 'Etat' ietat ;
  741. Message ;
  742. FINSI ;
  743. ************ierr = 10 ; errabs = uyu.ino ; QUITTER BECART ;
  744. FINSI ;
  745. SINON ;
  746. duy = ABS (uyc.ino - uyu.ino) ;
  747. errrel = ABS (duy / uyc.ino) ;
  748. SI (errrel > critere) ;
  749. SI (EGA IMESS 1) ;
  750. Message 'Erreur relative' errrel
  751. 'Composante ' (EXTR LDEPL icomp)
  752. ' Noeud' ino 'Etat' ietat ;
  753. Message 'uyc.ino =' uyc.ino 'uyu.ino =' uyu.ino ;
  754. Message ;
  755. FINSI ;
  756. ************ierr = 11 ; QUITTER BECART ;
  757. FINSI ;
  758. FINSI ;
  759. FIN BUY ;
  760.  
  761. icomp = 3 ;
  762. ino = 0 ;
  763. REPETER BUZ nbno ;
  764. ino = ino + 1 ;
  765. SI (EGA uzc.ino 0.0) ;
  766. SI (NON (EGA uzu.ino 0.0)) ;
  767. SI (EGA IMESS 1) ;
  768. Message 'Erreur absolue' uzu.ino
  769. 'Composante ' (EXTR LDEPL icomp)
  770. ' Noeud' ino 'Etat' ietat ;
  771. Message ;
  772. FINSI ;
  773. ************ierr = 10 ; errabs = uzu.ino ; QUITTER BECART ;
  774. FINSI ;
  775. SINON ;
  776. duz = ABS (uzc.ino - uzu.ino) ;
  777. errrel = ABS (duz / uzc.ino) ;
  778. SI (errrel > critere) ;
  779. SI (EGA IMESS 1) ;
  780. Message 'Erreur relative' errrel
  781. 'Composante ' (EXTR LDEPL icomp)
  782. ' Noeud' ino 'Etat' ietat ;
  783. Message 'uzc.ino =' uzc.ino 'uzu.ino =' uzu.ino ;
  784. Message ;
  785. FINSI ;
  786. ************ierr = 11 ; QUITTER BECART ;
  787. FINSI ;
  788. FINSI ;
  789. FIN BUZ ;
  790.  
  791. * Comparaison des deformations-----------------------------------------
  792.  
  793. icomp = 1 ;
  794. igau = 0 ;
  795. REPETER BEPXX nbgau ;
  796. igau = igau + 1 ;
  797. SI (EGA epxxc.igau 0.0) ;
  798. SI (NON (EGA epxxu.igau 0.0)) ;
  799. SI (EGA IMESS 1) ;
  800. Message 'Erreur absolue' epxxu.igau
  801. 'Composante ' (EXTR LDEFO icomp)
  802. ' Point Gauss' igau 'Etat' ietat ;
  803. Message ;
  804. FINSI ;
  805. ************ierr = 20 ; errabs = epxxu.igau ; QUITTER BECART ;
  806. FINSI ;
  807. SINON ;
  808. depxx = ABS (epxxc.igau - epxxu.igau) ;
  809. errrel = ABS (depxx / epxxc.igau) ;
  810. SI (errrel > critere) ;
  811. SI (EGA IMESS 1) ;
  812. Message 'Erreur relative' errrel
  813. 'Composante ' (EXTR LDEFO icomp)
  814. ' Point Gauss' igau 'Etat' ietat ;
  815. Message 'epxxc.igau =' epxxc.igau 'epxxu.igau =' epxxu.igau;
  816. Message ;
  817. FINSI ;
  818. ************ierr = 21 ; QUITTER BECART ;
  819. FINSI ;
  820. FINSI ;
  821. FIN BEPXX ;
  822.  
  823. icomp = 2 ;
  824. igau = 0 ;
  825. REPETER BEPYY nbgau ;
  826. igau = igau + 1 ;
  827. SI (EGA epyyc.igau 0.0) ;
  828. SI (NON (EGA epyyu.igau 0.0)) ;
  829. SI (EGA IMESS 1) ;
  830. Message 'Erreur absolue' epyyu.igau
  831. 'Composante ' (EXTR LDEFO icomp)
  832. ' Point Gauss' igau 'Etat' ietat ;
  833. Message ;
  834. FINSI ;
  835. ************ierr = 20 ; errabs = epyyu.igau ; QUITTER BECART ;
  836. FINSI ;
  837. SINON ;
  838. depyy = ABS (epyyc.igau - epyyu.igau) ;
  839. errrel = ABS (depyy / epyyc.igau) ;
  840. SI (errrel > critere) ;
  841. SI (EGA IMESS 1) ;
  842. Message 'Erreur relative' errrel
  843. 'Composante ' (EXTR LDEFO icomp)
  844. ' Point Gauss' igau 'Etat' ietat ;
  845. Message 'epyyc.igau =' epyyc.igau 'epyyu.igau =' epyyu.igau;
  846. Message ;
  847. FINSI ;
  848. ************ierr = 21 ; QUITTER BECART ;
  849. FINSI ;
  850. FINSI ;
  851. FIN BEPYY ;
  852.  
  853. icomp = 3 ;
  854. igau = 0 ;
  855. REPETER BEPZZ nbgau ;
  856. igau = igau + 1 ;
  857. SI (EGA epzzc.igau 0.0) ;
  858. SI (NON (EGA epzzu.igau 0.0)) ;
  859. SI (EGA IMESS 1) ;
  860. Message 'Erreur absolue' epzzu.igau
  861. 'Composante ' (EXTR LDEFO icomp)
  862. ' Point Gauss' igau 'Etat' ietat ;
  863. Message ;
  864. FINSI ;
  865. ************ierr = 20 ; errabs = epzzu.igau ; QUITTER BECART ;
  866. FINSI ;
  867. SINON ;
  868. depzz = ABS (epzzc.igau - epzzu.igau) ;
  869. errrel = ABS (depzz / epzzc.igau) ;
  870. SI (errrel > critere) ;
  871. SI (EGA IMESS 1) ;
  872. Message 'Erreur relative' errrel
  873. 'Composante ' (EXTR LDEFO icomp)
  874. ' Point Gauss' igau 'Etat' ietat ;
  875. Message 'epzzc.igau =' epzzc.igau 'epzzu.igau =' epzzu.igau;
  876. Message ;
  877. FINSI ;
  878. ************ierr = 21 ; QUITTER BECART ;
  879. FINSI ;
  880. FINSI ;
  881. FIN BEPZZ ;
  882.  
  883. icomp = 4 ;
  884. igau = 0 ;
  885. REPETER BGAXY nbgau ;
  886. igau = igau + 1 ;
  887. SI (EGA gaxyc.igau 0.0) ;
  888. SI (NON (EGA gaxyu.igau 0.0)) ;
  889. SI (EGA IMESS 1) ;
  890. Message 'Erreur absolue' gaxyu.igau
  891. 'Composante ' (EXTR LDEFO icomp)
  892. ' Point Gauss' igau 'Etat' ietat ;
  893. Message ;
  894. FINSI ;
  895. ************ierr = 20 ; errabs = gaxyu.igau ; QUITTER BECART ;
  896. FINSI ;
  897. SINON ;
  898. dgaxy = ABS (gaxyc.igau - gaxyu.igau) ;
  899. errrel = ABS (dgaxy / gaxyc.igau) ;
  900. SI (errrel > critere) ;
  901. SI (EGA IMESS 1) ;
  902. Message 'Erreur relative' errrel
  903. 'Composante ' (EXTR LDEFO icomp)
  904. ' Point Gauss' igau 'Etat' ietat ;
  905. Message 'gaxyc.igau =' gaxyc.igau 'gaxyu.igau =' gaxyu.igau;
  906. Message ;
  907. FINSI ;
  908. ************ierr = 21 ; QUITTER BECART ;
  909. FINSI ;
  910. FINSI ;
  911. FIN BGAXY ;
  912.  
  913. icomp = 5 ;
  914. igau = 0 ;
  915. REPETER BGAXZ nbgau ;
  916. igau = igau + 1 ;
  917. SI (EGA gaxzc.igau 0.0) ;
  918. SI (NON (EGA gaxzu.igau 0.0)) ;
  919. SI (EGA IMESS 1) ;
  920. Message 'Erreur absolue' gaxzu.igau
  921. 'Composante ' (EXTR LDEFO icomp)
  922. ' Point Gauss' igau 'Etat' ietat ;
  923. Message ;
  924. FINSI ;
  925. ************ierr = 20 ; errabs = gaxzu.igau ; QUITTER BECART ;
  926. FINSI ;
  927. SINON ;
  928. dgaxz = ABS (gaxzc.igau - gaxzu.igau) ;
  929. errrel = ABS (dgaxz / gaxzc.igau) ;
  930. SI (errrel > critere) ;
  931. SI (EGA IMESS 1) ;
  932. Message 'Erreur relative' errrel
  933. 'Composante ' (EXTR LDEFO icomp)
  934. ' Point Gauss' igau 'Etat' ietat ;
  935. Message 'gaxzc.igau =' gaxzc.igau 'gaxzu.igau =' gaxzu.igau;
  936. Message ;
  937. FINSI ;
  938. ************ierr = 21 ; QUITTER BECART ;
  939. FINSI ;
  940. FINSI ;
  941. FIN BGAXZ ;
  942.  
  943. icomp = 6 ;
  944. igau = 0 ;
  945. REPETER BGAYZ nbgau ;
  946. igau = igau + 1 ;
  947. SI (EGA gayzc.igau 0.0) ;
  948. SI (NON (EGA gayzu.igau 0.0)) ;
  949. ************SI (EGA IMESS 1) ;
  950. ************Message 'Erreur absolue' gayzu.igau
  951. ************ 'Composante ' (EXTR LDEFO icomp)
  952. ************ ' Point Gauss' igau 'Etat' ietat ;
  953. ************Message ;
  954. ************FINSI ;
  955. ierr = 20 ; errabs = gayzu.igau ; QUITTER BECART ;
  956. FINSI ;
  957. SINON ;
  958. dgayz = ABS (gayzc.igau - gayzu.igau) ;
  959. errrel = ABS (dgayz / gayzc.igau) ;
  960. SI (errrel > critere) ;
  961. ************SI (EGA IMESS 1) ;
  962. ************Message 'Erreur relative' errrel
  963. ************ 'Composante ' (EXTR LDEFO icomp)
  964. ************ ' Point Gauss' igau 'Etat' ietat ;
  965. ************Message 'gayzc.igau =' gayzc.igau 'gayzu.igau =' gayzu.igau;
  966. ************Message ;
  967. ************FINSI ;
  968. ierr = 21 ; QUITTER BECART ;
  969. FINSI ;
  970. FINSI ;
  971. FIN BGAYZ ;
  972.  
  973. * Comparaison des deformations inelastiques----------------------------
  974.  
  975. icomp = 1 ;
  976. igau = 0 ;
  977. REPETER BEIXX nbgau ;
  978. igau = igau + 1 ;
  979. SI (EGA eixxc.igau 0.0) ;
  980. SI (NON (EGA enxxu.igau 0.0)) ;
  981. SI (EGA IMESS 1) ;
  982. Message 'Erreur absolue' enxxu.igau
  983. 'Composante ' (EXTR LDEIN icomp)
  984. ' Point Gauss' igau 'Etat' ietat ;
  985. Message ;
  986. FINSI ;
  987. ************ierr = 30 ; errabs = enxxu.igau ; QUITTER BECART ;
  988. FINSI ;
  989. SINON ;
  990. deixx = ABS (eixxc.igau - enxxu.igau) ;
  991. errrel = ABS (deixx / eixxc.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 'eixxc.igau =' eixxc.igau 'enxxu.igau =' enxxu.igau;
  998. Message ;
  999. FINSI ;
  1000. ************ierr = 31 ; QUITTER BECART ;
  1001. FINSI ;
  1002. FINSI ;
  1003. FIN BEIXX ;
  1004.  
  1005. icomp = 2 ;
  1006. igau = 0 ;
  1007. REPETER BEIYY nbgau ;
  1008. igau = igau + 1 ;
  1009. SI (EGA eiyyc.igau 0.0) ;
  1010. SI (NON (EGA enyyu.igau 0.0)) ;
  1011. SI (EGA IMESS 1) ;
  1012. Message 'Erreur absolue' enyyu.igau
  1013. 'Composante ' (EXTR LDEIN icomp)
  1014. ' Point Gauss' igau 'Etat' ietat ;
  1015. Message ;
  1016. FINSI ;
  1017. ************ierr = 30 ; errabs = enyyu.igau ; QUITTER BECART ;
  1018. FINSI ;
  1019. SINON ;
  1020. deiyy = ABS (eiyyc.igau - enyyu.igau) ;
  1021. errrel = ABS (deiyy / eiyyc.igau) ;
  1022. SI (errrel > critere) ;
  1023. SI (EGA IMESS 1) ;
  1024. Message 'Erreur relative' errrel
  1025. 'Composante ' (EXTR LDEIN icomp)
  1026. ' Point Gauss' igau 'Etat' ietat ;
  1027. Message 'eiyyc.igau =' eiyyc.igau 'enyyu.igau =' enyyu.igau;
  1028. Message ;
  1029. FINSI ;
  1030. ************ierr = 31 ; QUITTER BECART ;
  1031. FINSI ;
  1032. FINSI ;
  1033. FIN BEIYY ;
  1034.  
  1035. icomp = 3 ;
  1036. igau = 0 ;
  1037. REPETER BEIZZ nbgau ;
  1038. igau = igau + 1 ;
  1039. SI (EGA eizzc.igau 0.0) ;
  1040. SI (NON (EGA enzzu.igau 0.0)) ;
  1041. SI (EGA IMESS 1) ;
  1042. Message 'Erreur absolue' enzzu.igau
  1043. 'Composante ' (EXTR LDEIN icomp)
  1044. ' Point Gauss' igau 'Etat' ietat ;
  1045. Message ;
  1046. FINSI ;
  1047. ************ierr = 30 ; errabs = enzzu.igau ; QUITTER BECART ;
  1048. FINSI ;
  1049. SINON ;
  1050. deizz = ABS (eizzc.igau - enzzu.igau) ;
  1051. errrel = ABS (deizz / eizzc.igau) ;
  1052. SI (errrel > critere) ;
  1053. SI (EGA IMESS 1) ;
  1054. Message 'Erreur relative' errrel
  1055. 'Composante ' (EXTR LDEIN icomp)
  1056. ' Point Gauss' igau 'Etat' ietat ;
  1057. Message 'eizzc.igau =' eizzc.igau 'enzzu.igau =' enzzu.igau;
  1058. Message ;
  1059. FINSI ;
  1060. ************ierr = 31 ; QUITTER BECART ;
  1061. FINSI ;
  1062. FINSI ;
  1063. FIN BEIZZ ;
  1064.  
  1065. icomp = 4 ;
  1066. igau = 0 ;
  1067. REPETER BGIXY nbgau ;
  1068. igau = igau + 1 ;
  1069. SI (EGA gixyc.igau 0.0) ;
  1070. SI (NON (EGA gnxyu.igau 0.0)) ;
  1071. SI (EGA IMESS 1) ;
  1072. Message 'Erreur absolue' gnxyu.igau
  1073. 'Composante ' (EXTR LDEIN icomp)
  1074. ' Point Gauss' igau 'Etat' ietat ;
  1075. Message ;
  1076. FINSI ;
  1077. ************ierr = 30 ; errabs = gnxyu.igau ; QUITTER BECART ;
  1078. FINSI ;
  1079. SINON ;
  1080. dgixy = ABS (gixyc.igau - gnxyu.igau) ;
  1081. errrel = ABS (dgixy / gixyc.igau) ;
  1082. SI (errrel > critere) ;
  1083. SI (EGA IMESS 1) ;
  1084. Message 'Erreur relative' errrel
  1085. 'Composante ' (EXTR LDEIN icomp)
  1086. ' Point Gauss' igau 'Etat' ietat ;
  1087. Message 'gixyc.igau =' gixyc.igau 'gnxyu.igau =' gnxyu.igau;
  1088. Message ;
  1089. FINSI ;
  1090. ************ierr = 31 ; QUITTER BECART ;
  1091. FINSI ;
  1092. FINSI ;
  1093. FIN BGIXY ;
  1094.  
  1095. icomp = 5 ;
  1096. igau = 0 ;
  1097. REPETER BGIXZ nbgau ;
  1098. igau = igau + 1 ;
  1099. SI (EGA gixzc.igau 0.0) ;
  1100. SI (NON (EGA gnxzu.igau 0.0)) ;
  1101. SI (EGA IMESS 1) ;
  1102. Message 'Erreur absolue' gnxzu.igau
  1103. 'Composante ' (EXTR LDEIN icomp)
  1104. ' Point Gauss' igau 'Etat' ietat ;
  1105. Message ;
  1106. FINSI ;
  1107. ************ierr = 30 ; errabs = gnxzu.igau ; QUITTER BECART ;
  1108. FINSI ;
  1109. SINON ;
  1110. dgixz = ABS (gixzc.igau - gnxzu.igau) ;
  1111. errrel = ABS (dgixz / gixzc.igau) ;
  1112. SI (errrel > critere) ;
  1113. SI (EGA IMESS 1) ;
  1114. Message 'Erreur relative' errrel
  1115. 'Composante ' (EXTR LDEIN icomp)
  1116. ' Point Gauss' igau 'Etat' ietat ;
  1117. Message 'gixzc.igau =' gixzc.igau 'gnxzu.igau =' gnxzu.igau;
  1118. Message ;
  1119. FINSI ;
  1120. ************ierr = 31 ; QUITTER BECART ;
  1121. FINSI ;
  1122. FINSI ;
  1123. FIN BGIXZ ;
  1124.  
  1125. icomp = 6 ;
  1126. igau = 0 ;
  1127. REPETER BGIYZ nbgau ;
  1128. igau = igau + 1 ;
  1129. SI (EGA giyzc.igau 0.0) ;
  1130. SI (NON (EGA gnyzu.igau 0.0)) ;
  1131. ************SI (EGA IMESS 1) ;
  1132. ************Message 'Erreur absolue' gnyzu.igau
  1133. ************ 'Composante ' (EXTR LDEIN icomp)
  1134. ************ ' Point Gauss' igau 'Etat' ietat ;
  1135. ************Message ;
  1136. ************FINSI ;
  1137. ierr = 30 ; errabs = gnyzu.igau ; QUITTER BECART ;
  1138. FINSI ;
  1139. SINON ;
  1140. dgiyz = ABS (giyzc.igau - gnyzu.igau) ;
  1141. errrel = ABS (dgiyz / giyzc.igau) ;
  1142. SI (errrel > critere) ;
  1143. ************SI (EGA IMESS 1) ;
  1144. ************Message 'Erreur relative' errrel
  1145. ************ 'Composante ' (EXTR LDEIN icomp)
  1146. ************ ' Point Gauss' igau 'Etat' ietat ;
  1147. ************Message 'giyzc.igau =' giyzc.igau 'gnyzu.igau =' gnyzu.igau;
  1148. ************Message ;
  1149. ************FINSI ;
  1150. ierr = 31 ; QUITTER BECART ;
  1151. FINSI ;
  1152. FINSI ;
  1153. FIN BGIYZ ;
  1154.  
  1155. * Comparaison des contraintes------------------------------------------
  1156.  
  1157. icomp = 1 ;
  1158. igau = 0 ;
  1159. REPETER BSMXX nbgau ;
  1160. igau = igau + 1 ;
  1161. SI (EGA smxxc.igau 0.0) ;
  1162. SI (NON (EGA smxxu.igau 0.0)) ;
  1163. SI (EGA IMESS 1) ;
  1164. Message 'Erreur absolue' smxxu.igau
  1165. 'Composante ' (EXTR LCONT icomp)
  1166. ' Point Gauss' igau 'Etat' ietat ;
  1167. Message ;
  1168. FINSI ;
  1169. ************ierr = 40 ; errabs = smxxu.igau ; QUITTER BECART ;
  1170. FINSI ;
  1171. SINON ;
  1172. dsmxx = ABS (smxxc.igau - smxxu.igau) ;
  1173. errrel = ABS (dsmxx / smxxc.igau) ;
  1174. SI (errrel > critere) ;
  1175. SI (EGA IMESS 1) ;
  1176. Message 'Erreur relative' errrel
  1177. 'Composante ' (EXTR LCONT icomp)
  1178. ' Point Gauss' igau 'Etat' ietat ;
  1179. Message 'smxxc.igau =' smxxc.igau 'smxxu.igau =' smxxu.igau;
  1180. Message ;
  1181. FINSI ;
  1182. ************ierr = 41 ; QUITTER BECART ;
  1183. FINSI ;
  1184. FINSI ;
  1185. FIN BSMXX ;
  1186.  
  1187. icomp = 2 ;
  1188. igau = 0 ;
  1189. REPETER BSMYY nbgau ;
  1190. igau = igau + 1 ;
  1191. SI (EGA smyyc.igau 0.0) ;
  1192. SI (NON (EGA smyyu.igau 0.0)) ;
  1193. SI (EGA IMESS 1) ;
  1194. Message 'Erreur absolue' smyyu.igau
  1195. 'Composante ' (EXTR LCONT icomp)
  1196. ' Point Gauss' igau 'Etat' ietat ;
  1197. Message ;
  1198. FINSI ;
  1199. ************ierr = 40 ; errabs = smyyu.igau ; QUITTER BECART ;
  1200. FINSI ;
  1201. SINON ;
  1202. dsmyy = ABS (smyyc.igau - smyyu.igau) ;
  1203. errrel = ABS (dsmyy / smyyc.igau) ;
  1204. SI (errrel > critere) ;
  1205. SI (EGA IMESS 1) ;
  1206. Message 'Erreur relative' errrel
  1207. 'Composante ' (EXTR LCONT icomp)
  1208. ' Point Gauss' igau 'Etat' ietat ;
  1209. Message 'smyyc.igau =' smyyc.igau 'smyyu.igau =' smyyu.igau;
  1210. Message ;
  1211. FINSI ;
  1212. ************ierr = 41 ; QUITTER BECART ;
  1213. FINSI ;
  1214. FINSI ;
  1215. FIN BSMYY ;
  1216.  
  1217. icomp = 3 ;
  1218. igau = 0 ;
  1219. REPETER BSMZZ nbgau ;
  1220. igau = igau + 1 ;
  1221. SI (EGA smzzc.igau 0.0) ;
  1222. SI (NON (EGA smzzu.igau 0.0)) ;
  1223. SI (EGA IMESS 1) ;
  1224. Message 'Erreur absolue' smzzu.igau
  1225. 'Composante ' (EXTR LCONT icomp)
  1226. ' Point Gauss' igau 'Etat' ietat ;
  1227. Message ;
  1228. FINSI ;
  1229. ************ierr = 40 ; errabs = smzzu.igau ; QUITTER BECART ;
  1230. FINSI ;
  1231. SINON ;
  1232. dsmzz = ABS (smzzc.igau - smzzu.igau) ;
  1233. errrel = ABS (dsmzz / smzzc.igau) ;
  1234. SI (errrel > critere) ;
  1235. SI (EGA IMESS 1) ;
  1236. Message 'Erreur relative' errrel
  1237. 'Composante ' (EXTR LCONT icomp)
  1238. ' Point Gauss' igau 'Etat' ietat ;
  1239. Message 'smzzc.igau =' smzzc.igau 'smzzu.igau =' smzzu.igau;
  1240. Message ;
  1241. FINSI ;
  1242. ************ierr = 41 ; QUITTER BECART ;
  1243. FINSI ;
  1244. FINSI ;
  1245. FIN BSMZZ ;
  1246.  
  1247. icomp = 4 ;
  1248. igau = 0 ;
  1249. REPETER BSMXY nbgau ;
  1250. igau = igau + 1 ;
  1251. SI (EGA smxyc.igau 0.0) ;
  1252. SI (NON (EGA smxyu.igau 0.0)) ;
  1253. SI (EGA IMESS 1) ;
  1254. Message 'Erreur absolue' smxyu.igau
  1255. 'Composante ' (EXTR LCONT icomp)
  1256. ' Point Gauss' igau 'Etat' ietat ;
  1257. Message ;
  1258. FINSI ;
  1259. ************ierr = 40 ; errabs = smxyu.igau ; QUITTER BECART ;
  1260. FINSI ;
  1261. SINON ;
  1262. dsmxy = ABS (smxyc.igau - smxyu.igau) ;
  1263. errrel = ABS (dsmxy / smxyc.igau) ;
  1264. SI (errrel > critere) ;
  1265. SI (EGA IMESS 1) ;
  1266. Message 'Erreur relative' errrel
  1267. 'Composante ' (EXTR LCONT icomp)
  1268. ' Point Gauss' igau 'Etat' ietat ;
  1269. Message 'smxyc.igau =' smxyc.igau 'smxyu.igau =' smxyu.igau;
  1270. Message ;
  1271. FINSI ;
  1272. ************ierr = 41 ; QUITTER BECART ;
  1273. FINSI ;
  1274. FINSI ;
  1275. FIN BSMXY ;
  1276.  
  1277. icomp = 5 ;
  1278. igau = 0 ;
  1279. REPETER BSMXZ nbgau ;
  1280. igau = igau + 1 ;
  1281. SI (EGA smxzc.igau 0.0) ;
  1282. SI (NON (EGA smxzu.igau 0.0)) ;
  1283. SI (EGA IMESS 1) ;
  1284. Message 'Erreur absolue' smxzu.igau
  1285. 'Composante ' (EXTR LCONT icomp)
  1286. ' Point Gauss' igau 'Etat' ietat ;
  1287. Message ;
  1288. FINSI ;
  1289. ************ierr = 40 ; errabs = smxzu.igau ; QUITTER BECART ;
  1290. FINSI ;
  1291. SINON ;
  1292. dsmxz = ABS (smxzc.igau - smxzu.igau) ;
  1293. errrel = ABS (dsmxz / smxzc.igau) ;
  1294. SI (errrel > critere) ;
  1295. SI (EGA IMESS 1) ;
  1296. Message 'Erreur relative' errrel
  1297. 'Composante ' (EXTR LCONT icomp)
  1298. ' Point Gauss' igau 'Etat' ietat ;
  1299. Message 'smxzc.igau =' smxzc.igau 'smxzu.igau =' smxzu.igau;
  1300. Message ;
  1301. FINSI ;
  1302. ************ierr = 41 ; QUITTER BECART ;
  1303. FINSI ;
  1304. FINSI ;
  1305. FIN BSMXZ ;
  1306.  
  1307. icomp = 6 ;
  1308. igau = 0 ;
  1309. REPETER BSMYZ nbgau ;
  1310. igau = igau + 1 ;
  1311. SI (EGA smyzc.igau 0.0) ;
  1312. SI (NON (EGA smyzu.igau 0.0)) ;
  1313. ************SI (EGA IMESS 1) ;
  1314. ************Message 'Erreur absolue' smyzu.igau
  1315. ************ 'Composante ' (EXTR LCONT icomp)
  1316. ************ ' Point Gauss' igau 'Etat' ietat ;
  1317. ************Message ;
  1318. ************FINSI ;
  1319. ierr = 40 ; errabs = smyzu.igau ; QUITTER BECART ;
  1320. FINSI ;
  1321. SINON ;
  1322. dsmyz = ABS (smyzc.igau - smyzu.igau) ;
  1323. errrel = ABS (dsmyz / smyzc.igau) ;
  1324. SI (errrel > critere) ;
  1325. ************SI (EGA IMESS 1) ;
  1326. ************Message 'Erreur relative' errrel
  1327. ************ 'Composante ' (EXTR LCONT icomp)
  1328. ************ ' Point Gauss' igau 'Etat' ietat ;
  1329. ************Message 'smyzc.igau =' smyzc.igau 'smyzu.igau =' smyzu.igau;
  1330. ************Message ;
  1331. ************FINSI ;
  1332. ierr = 41 ; QUITTER BECART ;
  1333. FINSI ;
  1334. FINSI ;
  1335. FIN BSMYZ ;
  1336.  
  1337. * Comparaison des variables internes-----------------------------------
  1338.  
  1339. icomp = 1 ;
  1340. igau = 0 ;
  1341. REPETER BEPSE nbgau ;
  1342. igau = igau + 1 ;
  1343. SI (EGA epsec.igau 0.0) ;
  1344. SI (NON (EGA epseu.igau 0.0)) ;
  1345. SI (EGA IMESS 1) ;
  1346. Message 'Erreur absolue' epseu.igau
  1347. 'Composante ' (EXTR LVARI icomp)
  1348. ' Point Gauss' igau 'Etat' ietat ;
  1349. Message ;
  1350. FINSI ;
  1351. ************ierr = 50 ; errabs = epseu.igau ; QUITTER BECART ;
  1352. FINSI ;
  1353. SINON ;
  1354. depse = ABS (epsec.igau - epseu.igau) ;
  1355. errrel = ABS (depse / epsec.igau) ;
  1356. SI (errrel > critere) ;
  1357. SI (EGA IMESS 1) ;
  1358. Message 'Erreur relative' errrel
  1359. 'Composante ' (EXTR LVARI icomp)
  1360. ' Point Gauss' igau 'Etat' ietat ;
  1361. Message 'epsec.igau =' epsec.igau 'epseu.igau =' epseu.igau;
  1362. Message ;
  1363. FINSI ;
  1364. ************ierr = 51 ; QUITTER BECART ;
  1365. FINSI ;
  1366. FINSI ;
  1367. FIN BEPSE ;
  1368.  
  1369. FIN BECART ;
  1370.  
  1371. * FIN BOUCLE COMPARAISON DES RESULTATS________________________________
  1372.  
  1373. * MESSAGES D'ERREUR_____________________________________________________
  1374.  
  1375. Message ;
  1376. Message 'CPU avec modele interne = ' CPUint 'centiemes de seconde' ;
  1377. Message 'CPU avec modele externe = ' CPUext 'centiemes de seconde' ;
  1378.  
  1379. ifin = netat-1 ;
  1380. giyzcmax = 0.0 ;
  1381. igau = 0 ;
  1382. REPETER BGIYZFC nbgau ;
  1383. igau = igau + 1 ;
  1384. giyz = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ifin)
  1385. 'GIYZ' 1 1 igau ;
  1386. SI (giyz > giyzcmax) ; giyzcmax = giyz ; FINSI ;
  1387. FIN BGIYZFC ;
  1388. giyzumax = 0.0 ;
  1389. igau = 0 ;
  1390. REPETER BGIYZFU nbgau ;
  1391. igau = igau + 1 ;
  1392. giyz = EXTR (TABU.'VARIABLES_INTERNES'.ifin)
  1393. 'GNYZ' 1 1 igau ;
  1394. SI (giyz > giyzumax) ; giyzumax = giyz ; FINSI ;
  1395. FIN BGIYZFU ;
  1396. Message ;
  1397. Message 'Deformation de fluage finale (maxi sur points de Gauss)' ;
  1398. Message 'GIYZ max avec modele interne = ' giyzcmax ;
  1399. Message 'GIYZ max avec modele externe = ' giyzumax ;
  1400. Message ;
  1401.  
  1402. SI (ierr > 0) ;
  1403.  
  1404. SI (EGA ierr 1) ;
  1405. Message 'Instants de calcul non coincidants' ;
  1406. ERRE 5 ;
  1407. FINSI ;
  1408.  
  1409. SI (EGA ierr 10) ;
  1410. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  1411. ' Noeud' ino 'Etat' ietat ;
  1412. ERRE 5 ;
  1413. FINSI ;
  1414.  
  1415. SI (EGA ierr 11) ;
  1416. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  1417. ' Noeud' ino 'Etat' ietat ;
  1418. ERRE 5 ;
  1419. FINSI ;
  1420.  
  1421. SI (EGA ierr 20) ;
  1422. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  1423. ' Point Gauss' igau 'Etat' ietat ;
  1424. ERRE 5 ;
  1425. FINSI ;
  1426.  
  1427. SI (EGA ierr 21) ;
  1428. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  1429. ' Point Gauss' igau 'Etat' ietat ;
  1430. ERRE 5 ;
  1431. FINSI ;
  1432.  
  1433. SI (EGA ierr 30) ;
  1434. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEIN icomp)
  1435. ' Point Gauss' igau 'Etat' ietat ;
  1436. ERRE 5 ;
  1437. FINSI ;
  1438.  
  1439. SI (EGA ierr 31) ;
  1440. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEIN icomp)
  1441. ' Point Gauss' igau 'Etat' ietat ;
  1442. ERRE 5 ;
  1443. FINSI ;
  1444.  
  1445. SI (EGA ierr 40) ;
  1446. Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  1447. ' Point Gauss' igau 'Etat' ietat ;
  1448. ERRE 5 ;
  1449. FINSI ;
  1450.  
  1451. SI (EGA ierr 41) ;
  1452. Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  1453. ' Point Gauss' igau 'Etat' ietat ;
  1454. ERRE 5 ;
  1455. FINSI ;
  1456.  
  1457. SI (EGA ierr 50) ;
  1458. Message 'Erreur absolue' errabs 'Composante ' (EXTR LVARI icomp)
  1459. ' Point Gauss' igau 'Etat' ietat ;
  1460. ERRE 5 ;
  1461. FINSI ;
  1462.  
  1463. SI (EGA ierr 51) ;
  1464. Message 'Erreur relative' errrel 'Composante ' (EXTR LVARI icomp)
  1465. ' Point Gauss' igau 'Etat' ietat ;
  1466. ERRE 5 ;
  1467. FINSI ;
  1468.  
  1469. SINON ;
  1470.  
  1471. Message 'Resultats concordants entre les deux modeles' ;
  1472. Message 'Aucune erreur relative superieure a ' critere
  1473. ' pour les composantes' ;
  1474. Message (EXTR LDEFO idefo) ' ' (EXTR LDEIN idein)
  1475. ' ' (EXTR LCONT icont) ;
  1476. ERRE 0 ;
  1477.  
  1478. FINSI ;
  1479.  
  1480. FIN ;
  1481.  
  1482.  
  1483.  
  1484.  
  1485.  
  1486.  
  1487.  
  1488.  
  1489.  
  1490.  

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