Télécharger creep02_cisXY.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : creep02_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 No2 : Modele de FLUAGE POLYNOMIAL
  15. * CISAILLEMENT XY
  16. * Le modele externalise est instanciable au niveau des
  17. * iterations internes (uniquement les coef. du polynome
  18. * de la loi de fluage)
  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 le modele 'FLUAGE' 'POLYNOMIAL' interne CAST3M
  90.  
  91. * 1.1 Definition du modele
  92.  
  93. mo_flupo = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  94. 'FLUAGE' 'POLYNOMIAL' ;
  95.  
  96. * 1.2 Affectation des proprietes de materiau
  97.  
  98. densiu = 0.95 ;
  99. dgrain = 8.0E-6 ;
  100. tuo2 = 1500.0 ;
  101. tauxfisk = 1.5E19 ;
  102.  
  103. xyoun = (2.2693E11 - (1.5399E7*tuo2)) - (9.35973E3 * (tuo2**2.0)) ;
  104. xcis = (8.583E10 - (5.157E6*tuo2)) - (3.747E3 * (tuo2**2.0)) ;
  105. xnu = (xyoun/(2.0*xcis)) - 1.0 ;
  106. xrho = 10950.0 ;
  107. xalph = ( (1.1833E-5 - (5.013E-9*tuo2)) + (3.756E-12 * (tuo2**2.0)) )
  108. - (6.125E-17 * (tuo2**3.0)) ;
  109.  
  110. xsmax = xyoun / 1000.0 ;
  111. xaf0 = 0.0 ;
  112. xaf1 = (EXP(0.183*(100.0-densiu)))
  113. * ( (0.2031E-13+(0.67861E-32*tauxfisk)) / (dgrain**2.0) )
  114. * (EXP(-45294.4*(tuo2**-1.0))) ;
  115. xaf2 = 1.0 ;
  116. xaf3 = (EXP(0.241*(100.0-densiu)))
  117. * (0.1524E-25+(0.50825E-44*tauxfisk))
  118. * (EXP(-66431.8*(tuo2**-1.0))) ;
  119. xaf4 = 4.5 ;
  120. xaf5 = (1.447E-34*tauxfisk) * (EXP(-3624.0*(tuo2**-1.0))) ;
  121. xaf6 = 1.0 ;
  122.  
  123. ma_flupo = MATE mo_flupo 'YOUN' xyoun 'NU ' xnu 'RHO ' xrho
  124. 'ALPH' xalph
  125. 'SMAX' xsmax 'AF0 ' xaf0 'AF1 ' xaf1
  126. 'AF2 ' xaf2 'AF3 ' xaf3 'AF4 ' xaf4
  127. 'AF5 ' xaf5 'AF6 ' xaf6 ;
  128.  
  129. * 1.3 Definition des conditions aux limites et chargements
  130.  
  131. cl1 = BLOQ mail1 'UY' 'UZ' ;
  132. cl2 = BLOQ sbloq 'UX' ;
  133. cl0 = cl1 ET cl2 ;
  134.  
  135. depmax = 5.0E8 ;
  136. dep0 = FORCE schar 'FX' depmax 'FY' 0.0 'FZ' 0.0 ;
  137. ev0 = EVOL MANU temps (PROG 0.0 1.E-2 1000.0)
  138. y (PROG 0.0 1.0 1.0 ) ;
  139. cha0 = CHAR MECA dep0 ev0 ;
  140.  
  141. * 1.4 Definition des instants de calcul
  142.  
  143. pro0 = PROG 0.0 PAS 1.E-6 3.0E-5 ;
  144. pro1 = PROG 0.0 PAS 1.E-6 3.0E-5 ;
  145.  
  146. * 1.5 Resolution par PASAPAS
  147.  
  148. tabflup = TABLE ;
  149. tabflup.'MODELE' = mo_flupo ;
  150. tabflup.'CARACTERISTIQUES' = ma_flupo ;
  151. tabflup.'CHARGEMENT' = cha0 ;
  152. tabflup.'BLOCAGES_MECANIQUES' = cl0 ;
  153. tabflup.'TEMPS_CALCULES' = pro0 ;
  154. tabflup.'TEMPS_SAUVES' = pro1 ;
  155. TMASAU=table;
  156. tabflup . 'MES_SAUVEGARDES'=TMASAU;
  157. TMASAU .'DEFTO'=VRAI;
  158. TMASAU .'DEFIN'=VRAI;
  159. TEMPS 'ZERO' ;
  160. PASAPAS tabflup ;
  161. TABTPS = TEMP 'NOEC';
  162. CPUint = TABTPS.'TEMPS_CPU'.'INITIAL' ;
  163.  
  164. *-----------------------------------------------------------------------
  165. * 2 - Calcul avec modele 'VISCO_EXTERNE' 'GENERAL' no22
  166. * Modele 'FLUAGE' 'POLYNOMIAL' externalise dans CREEP,
  167. * instanciable au niveau des iterations internes
  168.  
  169. * 2.1 Definition du modele utilisateur externalise dans CREEP
  170.  
  171. LCPAR22 = MOTS 'TFIS' 'TUO2' 'FACF' 'DSIU' 'DGRA' ;
  172.  
  173. mo_util = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  174. 'VISCO_EXTERNE' 'GENERAL'
  175. 'NUME_LOI' 22 'PARA_LOI' LCPAR22 ;
  176.  
  177. * 2.2 Affectation des proprietes de materiau
  178.  
  179. ma_util = MATE mo_util 'YOUN' xyoun 'NU ' xnu 'RHO ' xrho
  180. 'ALPH' xalph ;
  181.  
  182. * 2.3 Definition des conditions aux limites et chargements
  183. * IDEM 1.3 pour les CL et le chargement mecaniques
  184. * Chargements supplementaires correspondant aux parametres externes
  185.  
  186. chTFIS = MANU CHPO mail1 'TFIS' tauxfisk ;
  187. chTUO2 = MANU CHPO mail1 'TUO2' tuo2 ;
  188. chFACF = MANU CHPO mail1 'FACF' 1.0 ;
  189. chDSIU = MANU CHPO mail1 'DSIU' densiu ;
  190. chDGRA = MANU CHPO mail1 'DGRA' dgrain ;
  191.  
  192. chTFIS = CHAN ATTRIBUT chTFIS 'NATURE' 'DISCRET' ;
  193. chTUO2 = CHAN ATTRIBUT chTUO2 'NATURE' 'DISCRET' ;
  194. chFACF = CHAN ATTRIBUT chFACF 'NATURE' 'DISCRET' ;
  195. chDSIU = CHAN ATTRIBUT chDSIU 'NATURE' 'DISCRET' ;
  196. chDGRA = CHAN ATTRIBUT chDGRA 'NATURE' 'DISCRET' ;
  197.  
  198. Heavysid = EVOL MANU 'temps' (PROG 0.0 1000.0)
  199. 'signal' (PROG 1.0 1.0) ;
  200.  
  201. charTFIS = CHAR 'TFIS' chTFIS Heavysid ;
  202. charTUO2 = CHAR 'TUO2' chTUO2 Heavysid ;
  203. charFACF = CHAR 'FACF' chFACF Heavysid ;
  204. charDSIU = CHAR 'DSIU' chDSIU Heavysid ;
  205. charDGRA = CHAR 'DGRA' chDGRA Heavysid ;
  206.  
  207. charutil = cha0 ET
  208. charTFIS ET charTUO2 ET charFACF ET charDSIU ET charDGRA ;
  209.  
  210. * 2.4 Definition des instants de calcul : idem 1.4
  211.  
  212. * 2.5 Resolution par PASAPAS
  213.  
  214. tabutil = TABLE ;
  215. tabutil.'MODELE' = mo_util ;
  216. tabutil.'CARACTERISTIQUES' = ma_util ;
  217. tabutil.'CHARGEMENT' = charutil ;
  218. tabutil.'BLOCAGES_MECANIQUES' = cl0 ;
  219. tabutil.'TEMPS_CALCULES' = pro0 ;
  220. tabutil.'TEMPS_SAUVES' = pro1 ;
  221. TMASAU=table;
  222. tabutil . 'MES_SAUVEGARDES'=TMASAU;
  223. TMASAU .'DEFTO'=VRAI;
  224. TMASAU .'DEFIN'=VRAI;
  225. TEMPS 'ZERO' ;
  226. PASAPAS tabutil ;
  227. TABTPS = TEMP 'NOEC';
  228. CPUext = TABTPS.'TEMPS_CPU'.'INITIAL' ;
  229.  
  230. *-----------------------------------------------------------------------
  231. * 3 - Post-traitement : comparaisons entre les deux modeles
  232. *
  233. * N.B. Comme on simule un essai de cisaillement XY :
  234. * - Critere strict (errrel < 1.d-8) sur GAXY, GIXY et SMXY
  235. * - Pour les autres composantes on note les depassements du
  236. * critere si impression demandee
  237.  
  238. uxc = TABLE ; uyc = TABLE ; uzc = TABLE ;
  239. epxxc = TABLE ; epyyc = TABLE ; epzzc = TABLE ;
  240. gaxyc = TABLE ; gaxzc = TABLE ; gayzc = TABLE ;
  241. eixxc = TABLE ; eiyyc = TABLE ; eizzc = TABLE ;
  242. gixyc = TABLE ; gixzc = TABLE ; giyzc = TABLE ;
  243. smxxc = TABLE ; smyyc = TABLE ; smzzc = TABLE ;
  244. smxyc = TABLE ; smxzc = TABLE ; smyzc = TABLE ;
  245.  
  246. uxu = TABLE ; uyu = TABLE ; uzu = TABLE ;
  247. epxxu = TABLE ; epyyu = TABLE ; epzzu = TABLE ;
  248. gaxyu = TABLE ; gaxzu = TABLE ; gayzu = TABLE ;
  249. eixxu = TABLE ; eiyyu = TABLE ; eizzu = TABLE ;
  250. gixyu = TABLE ; gixzu = TABLE ; giyzu = TABLE ;
  251. smxxu = TABLE ; smyyu = TABLE ; smzzu = TABLE ;
  252. smxyu = TABLE ; smxzu = TABLE ; smyzu = TABLE ;
  253.  
  254. LDEPL = MOTS 'UX ' 'UY ' 'UZ ' ;
  255. LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  256. LDEIN = MOTS 'EIXX' 'EIYY' 'EIZZ' 'GIXY' 'GIXZ' 'GIYZ' ;
  257. LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  258.  
  259. idefo = 4 ;
  260. idein = 4 ;
  261. icont = 4 ;
  262.  
  263. ierr = 0 ;
  264. *critere = 1.d-8 ; kich
  265. critere = 2.d-8 ;
  266. IMESS = 0 ;
  267. ***IMESS = 1 ;
  268.  
  269. nbno = 8 ;
  270. nbgau = 8 ;
  271.  
  272. netat = DIME (tabflup.'TEMPS') ;
  273. ietat = -1 ;
  274.  
  275.  
  276. * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  277.  
  278. REPETER BECART netat ;
  279.  
  280. ietat = ietat + 1 ;
  281.  
  282. temps1 = tabflup.'TEMPS'.ietat ;
  283. temps2 = tabutil.'TEMPS'.ietat ;
  284. SI ((ABS (temps1-temps2)) > 1.0E-9) ;
  285. ierr = 1 ;
  286. QUITTER BECART ;
  287. FINSI ;
  288.  
  289. * Extraction des resultats de calcul : Modele interne CAST3M-----------
  290.  
  291. * Deplacements
  292. uxc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  293. uxc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  294. uxc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  295. uxc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  296. uxc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  297. uxc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  298. uxc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  299. uxc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  300.  
  301. uyc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  302. uyc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  303. uyc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  304. uyc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  305. uyc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  306. uyc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  307. uyc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  308. uyc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  309.  
  310. uzc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  311. uzc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  312. uzc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  313. uzc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  314. uzc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  315. uzc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  316. uzc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  317. uzc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  318.  
  319. * Deformations totales
  320. igau = 0 ;
  321. REPETER BEPXXC nbgau ;
  322. igau = igau + 1 ;
  323. epxxc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  324. FIN BEPXXC ;
  325.  
  326. igau = 0 ;
  327. REPETER BEPYYC nbgau ;
  328. igau = igau + 1 ;
  329. epyyc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  330. FIN BEPYYC ;
  331.  
  332. igau = 0 ;
  333. REPETER BEPZZC nbgau ;
  334. igau = igau + 1 ;
  335. epzzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  336. FIN BEPZZC ;
  337.  
  338. igau = 0 ;
  339. REPETER BGAXYC nbgau ;
  340. igau = igau + 1 ;
  341. gaxyc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  342. FIN BGAXYC ;
  343.  
  344. igau = 0 ;
  345. REPETER BGAXZC nbgau ;
  346. igau = igau + 1 ;
  347. gaxzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  348. FIN BGAXZC ;
  349.  
  350. igau = 0 ;
  351. REPETER BGAYZC nbgau ;
  352. igau = igau + 1 ;
  353. gayzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  354. FIN BGAYZC ;
  355.  
  356. * Deformations inelastiques
  357. igau = 0 ;
  358. REPETER BEIXXC nbgau ;
  359. igau = igau + 1 ;
  360. eixxc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  361. 'EIXX' 1 1 igau ;
  362. FIN BEIXXC ;
  363.  
  364. igau = 0 ;
  365. REPETER BEIYYC nbgau ;
  366. igau = igau + 1 ;
  367. eiyyc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  368. 'EIYY' 1 1 igau ;
  369. FIN BEIYYC ;
  370.  
  371. igau = 0 ;
  372. REPETER BEIZZC nbgau ;
  373. igau = igau + 1 ;
  374. eizzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  375. 'EIZZ' 1 1 igau ;
  376. FIN BEIZZC ;
  377.  
  378. igau = 0 ;
  379. REPETER BGIXYC nbgau ;
  380. igau = igau + 1 ;
  381. gixyc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  382. 'GIXY' 1 1 igau ;
  383. FIN BGIXYC ;
  384.  
  385. igau = 0 ;
  386. REPETER BGIXZC nbgau ;
  387. igau = igau + 1 ;
  388. gixzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  389. 'GIXZ' 1 1 igau ;
  390. FIN BGIXZC ;
  391.  
  392. igau = 0 ;
  393. REPETER BGIYZC nbgau ;
  394. igau = igau + 1 ;
  395. giyzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  396. 'GIYZ' 1 1 igau ;
  397. FIN BGIYZC ;
  398.  
  399. * Contraintes
  400. igau = 0 ;
  401. REPETER BSMXXC nbgau ;
  402. igau = igau + 1 ;
  403. smxxc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  404. FIN BSMXXC ;
  405.  
  406. igau = 0 ;
  407. REPETER BSMYYC nbgau ;
  408. igau = igau + 1 ;
  409. smyyc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  410. FIN BSMYYC ;
  411.  
  412. igau = 0 ;
  413. REPETER BSMZZC nbgau ;
  414. igau = igau + 1 ;
  415. smzzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  416. FIN BSMZZC ;
  417.  
  418. igau = 0 ;
  419. REPETER BSMXYC nbgau ;
  420. igau = igau + 1 ;
  421. smxyc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  422. FIN BSMXYC ;
  423.  
  424. igau = 0 ;
  425. REPETER BSMXZC nbgau ;
  426. igau = igau + 1 ;
  427. smxzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  428. FIN BSMXZC ;
  429.  
  430. igau = 0 ;
  431. REPETER BSMYZC nbgau ;
  432. igau = igau + 1 ;
  433. smyzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  434. FIN BSMYZC ;
  435.  
  436. * Extraction des resultats de calcul : Modele externe utilisateur------
  437.  
  438. * Deplacements
  439. uxu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  440. uxu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  441. uxu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  442. uxu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  443. uxu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  444. uxu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  445. uxu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  446. uxu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  447.  
  448. uyu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  449. uyu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  450. uyu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  451. uyu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  452. uyu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  453. uyu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  454. uyu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  455. uyu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  456.  
  457. uzu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  458. uzu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  459. uzu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  460. uzu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  461. uzu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  462. uzu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  463. uzu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  464. uzu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  465.  
  466. * Deformations totales
  467. igau = 0 ;
  468. REPETER BEPXXU nbgau ;
  469. igau = igau + 1 ;
  470. epxxu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  471. FIN BEPXXU ;
  472.  
  473. igau = 0 ;
  474. REPETER BEPYYU nbgau ;
  475. igau = igau + 1 ;
  476. epyyu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  477. FIN BEPYYU ;
  478.  
  479. igau = 0 ;
  480. REPETER BEPZZU nbgau ;
  481. igau = igau + 1 ;
  482. epzzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  483. FIN BEPZZU ;
  484.  
  485. igau = 0 ;
  486. REPETER BGAXYU nbgau ;
  487. igau = igau + 1 ;
  488. gaxyu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  489. FIN BGAXYU ;
  490.  
  491. igau = 0 ;
  492. REPETER BGAXZU nbgau ;
  493. igau = igau + 1 ;
  494. gaxzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  495. FIN BGAXZU ;
  496.  
  497. igau = 0 ;
  498. REPETER BGAYZU nbgau ;
  499. igau = igau + 1 ;
  500. gayzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  501. FIN BGAYZU ;
  502.  
  503. * Deformations inelastiques
  504. igau = 0 ;
  505. REPETER BEIXXU nbgau ;
  506. igau = igau + 1 ;
  507. eixxu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  508. 'EIXX' 1 1 igau ;
  509. FIN BEIXXU ;
  510.  
  511. igau = 0 ;
  512. REPETER BEIYYU nbgau ;
  513. igau = igau + 1 ;
  514. eiyyu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  515. 'EIYY' 1 1 igau ;
  516. FIN BEIYYU ;
  517.  
  518. igau = 0 ;
  519. REPETER BEIZZU nbgau ;
  520. igau = igau + 1 ;
  521. eizzu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  522. 'EIZZ' 1 1 igau ;
  523. FIN BEIZZU ;
  524.  
  525. igau = 0 ;
  526. REPETER BGIXYU nbgau ;
  527. igau = igau + 1 ;
  528. gixyu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  529. 'GIXY' 1 1 igau ;
  530. FIN BGIXYU ;
  531.  
  532. igau = 0 ;
  533. REPETER BGIXZU nbgau ;
  534. igau = igau + 1 ;
  535. gixzu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  536. 'GIXZ' 1 1 igau ;
  537. FIN BGIXZU ;
  538.  
  539. igau = 0 ;
  540. REPETER BGIYZU nbgau ;
  541. igau = igau + 1 ;
  542. giyzu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  543. 'GIYZ' 1 1 igau ;
  544. FIN BGIYZU ;
  545.  
  546. * Contraintes
  547. igau = 0 ;
  548. REPETER BSMXXU nbgau ;
  549. igau = igau + 1 ;
  550. smxxu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  551. FIN BSMXXU ;
  552.  
  553. igau = 0 ;
  554. REPETER BSMYYU nbgau ;
  555. igau = igau + 1 ;
  556. smyyu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  557. FIN BSMYYU ;
  558.  
  559. igau = 0 ;
  560. REPETER BSMZZU nbgau ;
  561. igau = igau + 1 ;
  562. smzzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  563. FIN BSMZZU ;
  564.  
  565. igau = 0 ;
  566. REPETER BSMXYU nbgau ;
  567. igau = igau + 1 ;
  568. smxyu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  569. FIN BSMXYU ;
  570.  
  571. igau = 0 ;
  572. REPETER BSMXZU nbgau ;
  573. igau = igau + 1 ;
  574. smxzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  575. FIN BSMXZU ;
  576.  
  577. igau = 0 ;
  578. REPETER BSMYZU nbgau ;
  579. igau = igau + 1 ;
  580. smyzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  581. FIN BSMYZU ;
  582.  
  583. * Comparaison des deplacements-----------------------------------------
  584.  
  585. icomp = 1 ;
  586. ino = 0 ;
  587. REPETER BUX nbno ;
  588. ino = ino + 1 ;
  589. SI (EGA uxc.ino 0.0) ;
  590. SI (NON (EGA uxu.ino 0.0)) ;
  591. SI (EGA IMESS 1) ;
  592. Message 'Erreur absolue' uxu.ino
  593. 'Composante ' (EXTR LDEPL icomp)
  594. ' Noeud' ino 'Etat' ietat ;
  595. Message ;
  596. FINSI ;
  597. ************ierr = 10 ; errabs = uxu.ino ; QUITTER BECART ;
  598. FINSI ;
  599. SINON ;
  600. dux = ABS (uxc.ino - uxu.ino) ;
  601. errrel = ABS (dux / uxc.ino) ;
  602. SI (errrel > critere) ;
  603. SI (EGA IMESS 1) ;
  604. Message 'Erreur relative' errrel
  605. 'Composante ' (EXTR LDEPL icomp)
  606. ' Noeud' ino 'Etat' ietat ;
  607. Message 'uxc.ino =' uxc.ino 'uxu.ino =' uxu.ino ;
  608. Message ;
  609. FINSI ;
  610. ************ierr = 11 ; QUITTER BECART ;
  611. FINSI ;
  612. FINSI ;
  613. FIN BUX ;
  614.  
  615. icomp = 2 ;
  616. ino = 0 ;
  617. REPETER BUY nbno ;
  618. ino = ino + 1 ;
  619. SI (EGA uyc.ino 0.0) ;
  620. SI (NON (EGA uyu.ino 0.0)) ;
  621. SI (EGA IMESS 1) ;
  622. Message 'Erreur absolue' uyu.ino
  623. 'Composante ' (EXTR LDEPL icomp)
  624. ' Noeud' ino 'Etat' ietat ;
  625. Message ;
  626. FINSI ;
  627. ************ierr = 10 ; errabs = uyu.ino ; QUITTER BECART ;
  628. FINSI ;
  629. SINON ;
  630. duy = ABS (uyc.ino - uyu.ino) ;
  631. errrel = ABS (duy / uyc.ino) ;
  632. SI (errrel > critere) ;
  633. SI (EGA IMESS 1) ;
  634. Message 'Erreur relative' errrel
  635. 'Composante ' (EXTR LDEPL icomp)
  636. ' Noeud' ino 'Etat' ietat ;
  637. Message 'uyc.ino =' uyc.ino 'uyu.ino =' uyu.ino ;
  638. Message ;
  639. FINSI ;
  640. ************ierr = 11 ; QUITTER BECART ;
  641. FINSI ;
  642. FINSI ;
  643. FIN BUY ;
  644.  
  645. icomp = 3 ;
  646. ino = 0 ;
  647. REPETER BUZ nbno ;
  648. ino = ino + 1 ;
  649. SI (EGA uzc.ino 0.0) ;
  650. SI (NON (EGA uzu.ino 0.0)) ;
  651. SI (EGA IMESS 1) ;
  652. Message 'Erreur absolue' uzu.ino
  653. 'Composante ' (EXTR LDEPL icomp)
  654. ' Noeud' ino 'Etat' ietat ;
  655. Message ;
  656. FINSI ;
  657. ************ierr = 10 ; errabs = uzu.ino ; QUITTER BECART ;
  658. FINSI ;
  659. SINON ;
  660. duz = ABS (uzc.ino - uzu.ino) ;
  661. errrel = ABS (duz / uzc.ino) ;
  662. SI (errrel > critere) ;
  663. SI (EGA IMESS 1) ;
  664. Message 'Erreur relative' errrel
  665. 'Composante ' (EXTR LDEPL icomp)
  666. ' Noeud' ino 'Etat' ietat ;
  667. Message 'uzc.ino =' uzc.ino 'uzu.ino =' uzu.ino ;
  668. Message ;
  669. FINSI ;
  670. ************ierr = 11 ; QUITTER BECART ;
  671. FINSI ;
  672. FINSI ;
  673. FIN BUZ ;
  674.  
  675. * Comparaison des deformations-----------------------------------------
  676.  
  677. icomp = 1 ;
  678. igau = 0 ;
  679. REPETER BEPXX nbgau ;
  680. igau = igau + 1 ;
  681. SI (EGA epxxc.igau 0.0) ;
  682. SI (NON (EGA epxxu.igau 0.0)) ;
  683. SI (EGA IMESS 1) ;
  684. Message 'Erreur absolue' epxxu.igau
  685. 'Composante ' (EXTR LDEFO icomp)
  686. ' Point Gauss' igau 'Etat' ietat ;
  687. Message ;
  688. FINSI ;
  689. ************ierr = 20 ; errabs = epxxu.igau ; QUITTER BECART ;
  690. FINSI ;
  691. SINON ;
  692. depxx = ABS (epxxc.igau - epxxu.igau) ;
  693. errrel = ABS (depxx / epxxc.igau) ;
  694. SI (errrel > critere) ;
  695. SI (EGA IMESS 1) ;
  696. Message 'Erreur relative' errrel
  697. 'Composante ' (EXTR LDEFO icomp)
  698. ' Point Gauss' igau 'Etat' ietat ;
  699. Message 'epxxc.igau =' epxxc.igau 'epxxu.igau =' epxxu.igau;
  700. Message ;
  701. FINSI ;
  702. ************ierr = 21 ; QUITTER BECART ;
  703. FINSI ;
  704. FINSI ;
  705. FIN BEPXX ;
  706.  
  707. icomp = 2 ;
  708. igau = 0 ;
  709. REPETER BEPYY nbgau ;
  710. igau = igau + 1 ;
  711. SI (EGA epyyc.igau 0.0) ;
  712. SI (NON (EGA epyyu.igau 0.0)) ;
  713. SI (EGA IMESS 1) ;
  714. Message 'Erreur absolue' epyyu.igau
  715. 'Composante ' (EXTR LDEFO icomp)
  716. ' Point Gauss' igau 'Etat' ietat ;
  717. Message ;
  718. FINSI ;
  719. ************ierr = 20 ; errabs = epyyu.igau ; QUITTER BECART ;
  720. FINSI ;
  721. SINON ;
  722. depyy = ABS (epyyc.igau - epyyu.igau) ;
  723. errrel = ABS (depyy / epyyc.igau) ;
  724. SI (errrel > critere) ;
  725. SI (EGA IMESS 1) ;
  726. Message 'Erreur relative' errrel
  727. 'Composante ' (EXTR LDEFO icomp)
  728. ' Point Gauss' igau 'Etat' ietat ;
  729. Message 'epyyc.igau =' epyyc.igau 'epyyu.igau =' epyyu.igau;
  730. Message ;
  731. FINSI ;
  732. ************ierr = 21 ; QUITTER BECART ;
  733. FINSI ;
  734. FINSI ;
  735. FIN BEPYY ;
  736.  
  737. icomp = 3 ;
  738. igau = 0 ;
  739. REPETER BEPZZ nbgau ;
  740. igau = igau + 1 ;
  741. SI (EGA epzzc.igau 0.0) ;
  742. SI (NON (EGA epzzu.igau 0.0)) ;
  743. SI (EGA IMESS 1) ;
  744. Message 'Erreur absolue' epzzu.igau
  745. 'Composante ' (EXTR LDEFO icomp)
  746. ' Point Gauss' igau 'Etat' ietat ;
  747. Message ;
  748. FINSI ;
  749. ************ierr = 20 ; errabs = epzzu.igau ; QUITTER BECART ;
  750. FINSI ;
  751. SINON ;
  752. depzz = ABS (epzzc.igau - epzzu.igau) ;
  753. errrel = ABS (depzz / epzzc.igau) ;
  754. SI (errrel > critere) ;
  755. SI (EGA IMESS 1) ;
  756. Message 'Erreur relative' errrel
  757. 'Composante ' (EXTR LDEFO icomp)
  758. ' Point Gauss' igau 'Etat' ietat ;
  759. Message 'epzzc.igau =' epzzc.igau 'epzzu.igau =' epzzu.igau;
  760. Message ;
  761. FINSI ;
  762. ************ierr = 21 ; QUITTER BECART ;
  763. FINSI ;
  764. FINSI ;
  765. FIN BEPZZ ;
  766.  
  767. icomp = 4 ;
  768. igau = 0 ;
  769. REPETER BGAXY nbgau ;
  770. igau = igau + 1 ;
  771. SI (EGA gaxyc.igau 0.0) ;
  772. SI (NON (EGA gaxyu.igau 0.0)) ;
  773. ************SI (EGA IMESS 1) ;
  774. ************Message 'Erreur absolue' gaxyu.igau
  775. ************ 'Composante ' (EXTR LDEFO icomp)
  776. ************ ' Point Gauss' igau 'Etat' ietat ;
  777. ************Message ;
  778. ************FINSI ;
  779. ierr = 20 ; errabs = gaxyu.igau ; QUITTER BECART ;
  780. FINSI ;
  781. SINON ;
  782. dgaxy = ABS (gaxyc.igau - gaxyu.igau) ;
  783. errrel = ABS (dgaxy / gaxyc.igau) ;
  784. SI (errrel > critere) ;
  785. ************SI (EGA IMESS 1) ;
  786. ************Message 'Erreur relative' errrel
  787. ************ 'Composante ' (EXTR LDEFO icomp)
  788. ************ ' Point Gauss' igau 'Etat' ietat ;
  789. ************Message 'gaxyc.igau =' gaxyc.igau 'gaxyu.igau =' gaxyu.igau;
  790. ************Message ;
  791. ************FINSI ;
  792. ierr = 21 ; QUITTER BECART ;
  793. FINSI ;
  794. FINSI ;
  795. FIN BGAXY ;
  796.  
  797. icomp = 5 ;
  798. igau = 0 ;
  799. REPETER BGAXZ nbgau ;
  800. igau = igau + 1 ;
  801. SI (EGA gaxzc.igau 0.0) ;
  802. SI (NON (EGA gaxzu.igau 0.0)) ;
  803. SI (EGA IMESS 1) ;
  804. Message 'Erreur absolue' gaxzu.igau
  805. 'Composante ' (EXTR LDEFO icomp)
  806. ' Point Gauss' igau 'Etat' ietat ;
  807. Message ;
  808. FINSI ;
  809. ************ierr = 20 ; errabs = gaxzu.igau ; QUITTER BECART ;
  810. FINSI ;
  811. SINON ;
  812. dgaxz = ABS (gaxzc.igau - gaxzu.igau) ;
  813. errrel = ABS (dgaxz / gaxzc.igau) ;
  814. SI (errrel > critere) ;
  815. SI (EGA IMESS 1) ;
  816. Message 'Erreur relative' errrel
  817. 'Composante ' (EXTR LDEFO icomp)
  818. ' Point Gauss' igau 'Etat' ietat ;
  819. Message 'gaxzc.igau =' gaxzc.igau 'gaxzu.igau =' gaxzu.igau;
  820. Message ;
  821. FINSI ;
  822. ************ierr = 21 ; QUITTER BECART ;
  823. FINSI ;
  824. FINSI ;
  825. FIN BGAXZ ;
  826.  
  827. icomp = 6 ;
  828. igau = 0 ;
  829. REPETER BGAYZ nbgau ;
  830. igau = igau + 1 ;
  831. SI (EGA gayzc.igau 0.0) ;
  832. SI (NON (EGA gayzu.igau 0.0)) ;
  833. SI (EGA IMESS 1) ;
  834. Message 'Erreur absolue' gayzu.igau
  835. 'Composante ' (EXTR LDEFO icomp)
  836. ' Point Gauss' igau 'Etat' ietat ;
  837. Message ;
  838. FINSI ;
  839. ************ierr = 20 ; errabs = gayzu.igau ; QUITTER BECART ;
  840. FINSI ;
  841. SINON ;
  842. dgayz = ABS (gayzc.igau - gayzu.igau) ;
  843. errrel = ABS (dgayz / gayzc.igau) ;
  844. SI (errrel > critere) ;
  845. SI (EGA IMESS 1) ;
  846. Message 'Erreur relative' errrel
  847. 'Composante ' (EXTR LDEFO icomp)
  848. ' Point Gauss' igau 'Etat' ietat ;
  849. Message 'gayzc.igau =' gayzc.igau 'gayzu.igau =' gayzu.igau;
  850. Message ;
  851. FINSI ;
  852. ************ierr = 21 ; QUITTER BECART ;
  853. FINSI ;
  854. FINSI ;
  855. FIN BGAYZ ;
  856.  
  857. * Comparaison des deformations inelastiques----------------------------
  858.  
  859. icomp = 1 ;
  860. igau = 0 ;
  861. REPETER BEIXX nbgau ;
  862. igau = igau + 1 ;
  863. SI (EGA eixxc.igau 0.0) ;
  864. SI (NON (EGA eixxu.igau 0.0)) ;
  865. SI (EGA IMESS 1) ;
  866. Message 'Erreur absolue' eixxu.igau
  867. 'Composante ' (EXTR LDEIN icomp)
  868. ' Point Gauss' igau 'Etat' ietat ;
  869. Message ;
  870. FINSI ;
  871. ************ierr = 30 ; errabs = eixxu.igau ; QUITTER BECART ;
  872. FINSI ;
  873. SINON ;
  874. deixx = ABS (eixxc.igau - eixxu.igau) ;
  875. errrel = ABS (deixx / eixxc.agau<+sPan>) ;
  876. SI (errrel > critere) ;
  877. SI (EGA IMESS 1) ;
  878. Message 'Erreur relative' errrel
  879. 'Composante ' (EXTR LDEIN icomp)
  880. ' Point Gauss' igau 'Etat' ietat ;
  881. Message 'eixxc.igau =' eixxc.igau 'eixxu.igau =' eixxu.igau;
  882. Message ;
  883. FINSI ;
  884. ************ierr = 31 ; QUITTER BECART ;
  885. FINSI ;
  886. FINSI ;
  887. FIN BEIXX ;
  888.  
  889. icomp = 2 ;
  890. igau = 0 ;
  891. REPETER BEIYY nbgau ;
  892. igau = igau + 1 ;
  893. SI (EGA eiyyc.igau 0.0) ;
  894. SI (NON (EGA eiyyu.igau 0.0)) ;
  895. SI (EGA IMESS 1) ;
  896. Message 'Erreur absolue' eiyyu.igau
  897. 'Composante ' (EXTR LDEIN icomp)
  898. ' Point Gauss' igau 'Etat' ietat ;
  899. Message ;
  900. FINSI ;
  901. ************ierr = 30 ; errabs = eiyyu.igau ; QUITTER BECART ;
  902. FINSI ;
  903. SINON ;
  904. deiyy = ABS (eiyyc.igau - eiyyu.igau) ;
  905. errrel = ABS (deiyy / eiyyc.igau) ;
  906. SI (errrel > critere) ;
  907. SI (EGA IMESS 1) ;
  908. Message 'Erreur relative' errrel
  909. 'Composante ' (EXTR LDEIN icomp)
  910. ' Point Gauss' igau 'Etat' ietat ;
  911. Message 'eiyyc.igau =' eiyyc.igau 'eiyyu.igau =' eiyyu.igau;
  912. Message ;
  913. FINSI ;
  914. ************ierr = 31 ; QUITTER BECART ;
  915. FINSI ;
  916. FINSI ;
  917. FIN BEIYY ;
  918.  
  919. icomp = 3 ;
  920. igau = 0 ;
  921. REPETER BEIZZ nbgau ;
  922. igau = igau + 1 ;
  923. SI (EGA eizzc.igau 0.0) ;
  924. SI (NON (EGA eizzu.igau 0.0)) ;
  925. SI (EGA IMESS 1) ;
  926. Message 'Erreur absolue' eizzu.igau
  927. 'Composante ' (EXTR LDEIN icomp)
  928. ' Point Gauss' igau 'Etat' ietat ;
  929. Message ;
  930. FINSI ;
  931. ************ierr = 30 ; errabs = eizzu.igau ; QUITTER BECART ;
  932. FINSI ;
  933. SINON ;
  934. deizz = ABS (eizzc.igau - eizzu.igau) ;
  935. errrel = ABS (deizz / eizzc.igau) ;
  936. SI (errrel > critere) ;
  937. SI (EGA IMESS 1) ;
  938. Message 'Erreur relative' errrel
  939. 'Composante ' (EXTR LDEIN icomp)
  940. ' Point Gauss' igau 'Etat' ietat ;
  941. Message 'eizzc.igau =' eizzc.igau 'eizzu.igau =' eizzu.igau;
  942. Message ;
  943. FINSI ;
  944. ************ierr = 31 ; QUITTER BECART ;
  945. FINSI ;
  946. FINSI ;
  947. FIN BEIZZ ;
  948.  
  949. icomp = 4 ;
  950. igau = 0 ;
  951. REPETER BGIXY nbgau ;
  952. igau = igau + 1 ;
  953. SI (EGA gixyc.igau 0.0) ;
  954. SI (NON (EGA gixyu.igau 0.0)) ;
  955. ************SI (EGA IMESS 1) ;
  956. ************Message 'Erreur absolue' gixyu.igau
  957. ************ 'Composante ' (EXTR LDEIN icomp)
  958. ************ ' Point Gauss' igau 'Etat' ietat ;
  959. ************Message ;
  960. ************FINSI ;
  961. ierr = 30 ; errabs = gixyu.igau ; QUITTER BECART ;
  962. FINSI ;
  963. SINON ;
  964. dgixy = ABS (gixyc.igau - gixyu.igau) ;
  965. errrel = ABS (dgixy / gixyc.igau) ;
  966. SI (errrel > critere) ;
  967. ************SI (EGA IMESS 1) ;
  968. ************Message 'Erreur relative' errrel
  969. ************ 'Composante ' (EXTR LDEIN icomp)
  970. ************ ' Point Gauss' igau 'Etat' ietat ;
  971. ************Message 'gixyc.igau =' gixyc.igau 'gixyu.igau =' gixyu.igau;
  972. ************Message ;
  973. ************FINSI ;
  974. ierr = 31 ; QUITTER BECART ;
  975. FINSI ;
  976. FINSI ;
  977. FIN BGIXY ;
  978.  
  979. icomp = 5 ;
  980. igau = 0 ;
  981. REPETER BGIXZ nbgau ;
  982. igau = igau + 1 ;
  983. SI (EGA gixzc.igau 0.0) ;
  984. SI (NON (EGA gixzu.igau 0.0)) ;
  985. SI (EGA IMESS 1) ;
  986. Message 'Erreur absolue' gixzu.igau
  987. 'Composante ' (EXTR LDEIN icomp)
  988. ' Point Gauss' igau 'Etat' ietat ;
  989. Message ;
  990. FINSI ;
  991. ************ierr = 30 ; errabs = gixzu.igau ; QUITTER BECART ;
  992. FINSI ;
  993. SINON ;
  994. dgixz = ABS (gixzc.igau - gixzu.igau) ;
  995. errrel = ABS (dgixz / gixzc.igau) ;
  996. SI (errrel > critere) ;
  997. SI (EGA IMESS 1) ;
  998. Message 'Erreur relative' errrel
  999. 'Composante ' (EXTR LDEIN icomp)
  1000. ' Point Gauss' igau 'Etat' ietat ;
  1001. Message 'gixzc.igau =' gixzc.igau 'gixzu.igau =' gixzu.igau;
  1002. Message ;
  1003. FINSI ;
  1004. ************ierr = 31 ; QUITTER BECART ;
  1005. FINSI ;
  1006. FINSI ;
  1007. FIN BGIXZ ;
  1008.  
  1009. icomp = 6 ;
  1010. igau = 0 ;
  1011. REPETER BGIYZ nbgau ;
  1012. igau = igau + 1 ;
  1013. SI (EGA giyzc.igau 0.0) ;
  1014. SI (NON (EGA giyzu.igau 0.0)) ;
  1015. SI (EGA IMESS 1) ;
  1016. Message 'Erreur absolue' giyzu.igau
  1017. 'Composante ' (EXTR LDEIN icomp)
  1018. ' Point Gauss' igau 'Etat' ietat ;
  1019. Message ;
  1020. FINSI ;
  1021. ************ierr = 30 ; errabs = giyzu.igau ; QUITTER BECART ;
  1022. FINSI ;
  1023. SINON ;
  1024. dgiyz = ABS (giyzc.igau - giyzu.igau) ;
  1025. errrel = ABS (dgiyz / giyzc.igau) ;
  1026. SI (errrel > critere) ;
  1027. SI (EGA IMESS 1) ;
  1028. Message 'Erreur relative' errrel
  1029. 'Composante ' (EXTR LDEIN icomp)
  1030. ' Point Gauss' igau 'Etat' ietat ;
  1031. Message 'giyzc.igau =' giyzc.igau 'giyzu.igau =' giyzu.igau;
  1032. Message ;
  1033. FINSI ;
  1034. ************ierr = 31 ; QUITTER BECART ;
  1035. FINSI ;
  1036. FINSI ;
  1037. FIN BGIYZ ;
  1038.  
  1039. * Comparaison des contraintes------------------------------------------
  1040.  
  1041. icomp = 1 ;
  1042. igau = 0 ;
  1043. REPETER BSMXX nbgau ;
  1044. igau = igau + 1 ;
  1045. SI (EGA smxxc.igau 0.0) ;
  1046. SI (NON (EGA smxxu.igau 0.0)) ;
  1047. SI (EGA IMESS 1) ;
  1048. Message 'Erreur absolue' smxxu.igau
  1049. 'Composante ' (EXTR LCONT icomp)
  1050. ' Point Gauss' igau 'Etat' ietat ;
  1051. Message ;
  1052. FINSI ;
  1053. ************ierr = 40 ; errabs = smxxu.igau ; QUITTER BECART ;
  1054. FINSI ;
  1055. SINON ;
  1056. dsmxx = ABS (smxxc.igau - smxxu.igau) ;
  1057. errrel = ABS (dsmxx / smxxc.igau) ;
  1058. SI (errrel > critere) ;
  1059. SI (EGA IMESS 1) ;
  1060. Message 'Erreur relative' errrel
  1061. 'Composante ' (EXTR LCONT icomp)
  1062. ' Point Gauss' igau 'Etat' ietat ;
  1063. Message 'smxxc.igau =' smxxc.igau 'smxxu.igau =' smxxu.igau;
  1064. Message ;
  1065. FINSI ;
  1066. ************ierr = 41 ; QUITTER BECART ;
  1067. FINSI ;
  1068. FINSI ;
  1069. FIN BSMXX ;
  1070.  
  1071. icomp = 2 ;
  1072. igau = 0 ;
  1073. REPETER BSMYY nbgau ;
  1074. igau = igau + 1 ;
  1075. SI (EGA smyyc.igau 0.0) ;
  1076. SI (NON (EGA smyyu.igau 0.0)) ;
  1077. SI (EGA IMESS 1) ;
  1078. Message 'Erreur absolue' smyyu.igau
  1079. 'Composante ' (EXTR LCONT icomp)
  1080. ' Point Gauss' igau 'Etat' ietat ;
  1081. Message ;
  1082. FINSI ;
  1083. ************ierr = 40 ; errabs = smyyu.igau ; QUITTER BECART ;
  1084. FINSI ;
  1085. SINON ;
  1086. dsmyy = ABS (smyyc.igau - smyyu.igau) ;
  1087. errrel = ABS (dsmyy / smyyc.igau) ;
  1088. SI (errrel > critere) ;
  1089. SI (EGA IMESS 1) ;
  1090. Message 'Erreur relative' errrel
  1091. 'Composante ' (EXTR LCONT icomp)
  1092. ' Point Gauss' igau 'Etat' ietat ;
  1093. Message 'smyyc.igau =' smyyc.igau 'smyyu.igau =' smyyu.igau;
  1094. Message ;
  1095. FINSI ;
  1096. ************ierr = 41 ; QUITTER BECART ;
  1097. FINSI ;
  1098. FINSI ;
  1099. FIN BSMYY ;
  1100.  
  1101. icomp = 3 ;
  1102. igau = 0 ;
  1103. REPETER BSMZZ nbgau ;
  1104. igau = igau + 1 ;
  1105. SI (EGA smzzc.igau 0.0) ;
  1106. SI (NON (EGA smzzu.igau 0.0)) ;
  1107. SI (EGA IMESS 1) ;
  1108. Message 'Erreur absolue' smzzu.igau
  1109. 'Composante ' (EXTR LCONT icomp)
  1110. ' Point Gauss' igau 'Etat' ietat ;
  1111. Message ;
  1112. FINSI ;
  1113. ************ierr = 40 ; errabs = smzzu.igau ; QUITTER BECART ;
  1114. FINSI ;
  1115. SINON ;
  1116. dsmzz = ABS (smzzc.igau - smzzu.igau) ;
  1117. errrel = ABS (dsmzz / smzzc.igau) ;
  1118. SI (errrel > critere) ;
  1119. SI (EGA IMESS 1) ;
  1120. Message 'Erreur relative' errrel
  1121. 'Composante ' (EXTR LCONT icomp)
  1122. ' Point Gauss' igau 'Etat' ietat ;
  1123. Message 'smzzc.igau =' smzzc.igau 'smzzu.igau =' smzzu.igau;
  1124. Message ;
  1125. FINSI ;
  1126. ************ierr = 41 ; QUITTER BECART ;
  1127. FINSI ;
  1128. FINSI ;
  1129. FIN BSMZZ ;
  1130.  
  1131. icomp = 4 ;
  1132. igau = 0 ;
  1133. REPETER BSMXY nbgau ;
  1134. igau = igau + 1 ;
  1135. SI (EGA smxyc.igau 0.0) ;
  1136. SI (NON (EGA smxyu.igau 0.0)) ;
  1137. ************SI (EGA IMESS 1) ;
  1138. ************Message 'Erreur absolue' smxyu.igau
  1139. ************ 'Composante ' (EXTR LCONT icomp)
  1140. ************ ' Point Gauss' igau 'Etat' ietat ;
  1141. ************Message ;
  1142. ************FINSI ;
  1143. ierr = 40 ; errabs = smxyu.igau ; QUITTER BECART ;
  1144. FINSI ;
  1145. SINON ;
  1146. dsmxy = ABS (smxyc.igau - smxyu.igau) ;
  1147. errrel = ABS (dsmxy / smxyc.igau) ;
  1148. SI (errrel > critere) ;
  1149. ************SI (EGA IMESS 1) ;
  1150. ************Message 'Erreur relative' errrel
  1151. ************ 'Composante ' (EXTR LCONT icomp)
  1152. ************ ' Point Gauss' igau 'Etat' ietat ;
  1153. ************Message 'smxyc.igau =' smxyc.igau 'smxyu.igau =' smxyu.igau;
  1154. ************Message ;
  1155. ************FINSI ;
  1156. ierr = 41 ; QUITTER BECART ;
  1157. FINSI ;
  1158. FINSI ;
  1159. FIN BSMXY ;
  1160.  
  1161. icomp = 5 ;
  1162. igau = 0 ;
  1163. REPETER BSMXZ nbgau ;
  1164. igau = igau + 1 ;
  1165. SI (EGA smxzc.igau 0.0) ;
  1166. SI (NON (EGA smxzu.igau 0.0)) ;
  1167. SI (EGA IMESS 1) ;
  1168. Message 'Erreur absolue' smxzu.igau
  1169. 'Composante ' (EXTR LCONT icomp)
  1170. ' Point Gauss' igau 'Etat' ietat ;
  1171. Message ;
  1172. FINSI ;
  1173. ************ierr = 40 ; errabs = smxzu.igau ; QUITTER BECART ;
  1174. FINSI ;
  1175. SINON ;
  1176. dsmxz = ABS (smxzc.igau - smxzu.igau) ;
  1177. errrel = ABS (dsmxz / smxzc.igau) ;
  1178. SI (errrel > critere) ;
  1179. SI (EGA IMESS 1) ;
  1180. Message 'Erreur relative' errrel
  1181. 'Composante ' (EXTR LCONT icomp)
  1182. ' Point Gauss' igau 'Etat' ietat ;
  1183. Message 'smxzc.igau =' smxzc.igau 'smxzu.igau =' smxzu.igau;
  1184. Message ;
  1185. FINSI ;
  1186. ************ierr = 41 ; QUITTER BECART ;
  1187. FINSI ;
  1188. FINSI ;
  1189. FIN BSMXZ ;
  1190.  
  1191. icomp = 6 ;
  1192. igau = 0 ;
  1193. REPETER BSMYZ nbgau ;
  1194. igau = igau + 1 ;
  1195. SI (EGA smyzc.igau 0.0) ;
  1196. SI (NON (EGA smyzu.igau 0.0)) ;
  1197. SI (EGA IMESS 1) ;
  1198. Message 'Erreur absolue' smyzu.igau
  1199. 'Composante ' (EXTR LCONT icomp)
  1200. ' Point Gauss' igau 'Etat' ietat ;
  1201. Message ;
  1202. FINSI ;
  1203. ************ierr = 40 ; errabs = smyzu.igau ; QUITTER BECART ;
  1204. FINSI ;
  1205. SINON ;
  1206. dsmyz = ABS (smyzc.igau - smyzu.igau) ;
  1207. errrel = ABS (dsmyz / smyzc.igau) ;
  1208. SI (errrel > critere) ;
  1209. SI (EGA IMESS 1) ;
  1210. Message 'Erreur relative' errrel
  1211. 'Composante ' (EXTR LCONT icomp)
  1212. ' Point Gauss' igau 'Etat' ietat ;
  1213. Message 'smyzc.igau =' smyzc.igau 'smyzu.igau =' smyzu.igau;
  1214. Message ;
  1215. FINSI ;
  1216. ************ierr = 41 ; QUITTER BECART ;
  1217. FINSI ;
  1218. FINSI ;
  1219. FIN BSMYZ ;
  1220.  
  1221. FIN BECART ;
  1222.  
  1223. * FIN BOUCLE COMPARAISON DES RESULTATS________________________________
  1224.  
  1225. * MESSAGES D'ERREUR_____________________________________________________
  1226.  
  1227. Message ;
  1228. Message 'CPU avec modele interne = ' CPUint 'centiemes de seconde' ;
  1229. Message 'CPU avec modele externe = ' CPUext 'centiemes de seconde' ;
  1230. Message ;
  1231.  
  1232. ifin = netat-1 ;
  1233.  
  1234. gaxycmax = 0.0 ;
  1235. igau = 0 ;
  1236. REPETER BGAXYFC nbgau ;
  1237. igau = igau + 1 ;
  1238. gaxy = EXTR (tabflup.'DEFORMATIONS'.ifin)
  1239. 'GAXY' 1 1 igau ;
  1240. SI (gaxy > gaxycmax) ; gaxycmax = gaxy ; FINSI ;
  1241. FIN BGAXYFC ;
  1242. gaxyumax = 0.0 ;
  1243. igau = 0 ;
  1244. REPETER BGAXYFU nbgau ;
  1245. igau = igau + 1 ;
  1246. gaxy = EXTR (tabutil.'DEFORMATIONS'.ifin)
  1247. 'GAXY' 1 1 igau ;
  1248. SI (gaxy > gaxyumax) ; gaxyumax = gaxy ; FINSI ;
  1249. FIN BGAXYFU ;
  1250. Message 'Deformation totale finale (maxi sur points de Gauss)' ;
  1251. Message 'GAXY max avec modele interne = ' gaxycmax ;
  1252. Message 'GAXY max avec modele externe = ' gaxyumax ;
  1253. Message ;
  1254.  
  1255. gixycmax = 0.0 ;
  1256. igau = 0 ;
  1257. REPETER BGIXYFC nbgau ;
  1258. igau = igau + 1 ;
  1259. gixy = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ifin)
  1260. 'GIXY' 1 1 igau ;
  1261. SI (gixy > gixycmax) ; gixycmax = gixy ; FINSI ;
  1262. FIN BGIXYFC ;
  1263. gixyumax = 0.0 ;
  1264. igau = 0 ;
  1265. REPETER BGIXYFU nbgau ;
  1266. igau = igau + 1 ;
  1267. gixy = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ifin)
  1268. 'GIXY' 1 1 igau ;
  1269. SI (gixy > gixyumax) ; gixyumax = gixy ; FINSI ;
  1270. FIN BGIXYFU ;
  1271. Message 'Deformation de fluage finale (maxi sur points de Gauss)' ;
  1272. Message 'GIXY max avec modele interne = ' gixycmax ;
  1273. Message 'GIXY max avec modele externe = ' gixyumax ;
  1274. Message ;
  1275.  
  1276. smxycmax = 0.0 ;
  1277. igau = 0 ;
  1278. REPETER BSMXYFC nbgau ;
  1279. igau = igau + 1 ;
  1280. smxy = EXTR (tabflup.'CONTRAINTES'.ifin)
  1281. 'SMXY' 1 1 igau ;
  1282. SI (smxy > smxycmax) ; smxycmax = smxy ; FINSI ;
  1283. FIN BSMXYFC ;
  1284. smxyumax = 0.0 ;
  1285. igau = 0 ;
  1286. REPETER BSMXYFU nbgau ;
  1287. igau = igau + 1 ;
  1288. smxy = EXTR (tabutil.'CONTRAINTES'.ifin)
  1289. 'SMXY' 1 1 igau ;
  1290. SI (smxy > smxyumax) ; smxyumax = smxy ; FINSI ;
  1291. FIN BSMXYFU ;
  1292. Message 'Contrainte finale (maxi sur points de Gauss)' ;
  1293. Message 'SMXY max avec modele interne = ' smxycmax ;
  1294. Message 'SMXY max avec modele externe = ' smxyumax ;
  1295. Message ;
  1296.  
  1297. SI (ierr > 0) ;
  1298.  
  1299. SI (EGA ierr 1) ;
  1300. Message 'Instants de calcul non coincidants' ;
  1301. ERRE 5 ;
  1302. FINSI ;
  1303.  
  1304. SI (EGA ierr 10) ;
  1305. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  1306. ' Noeud' ino 'Etat' ietat ;
  1307. ERRE 5 ;
  1308. FINSI ;
  1309.  
  1310. SI (EGA ierr 11) ;
  1311. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  1312. ' Noeud' ino 'Etat' ietat ;
  1313. ERRE 5 ;
  1314. FINSI ;
  1315.  
  1316. SI (EGA ierr 20) ;
  1317. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  1318. ' Point Gauss' igau 'Etat' ietat ;
  1319. ERRE 5 ;
  1320. FINSI ;
  1321.  
  1322. SI (EGA ierr 21) ;
  1323. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  1324. ' Point Gauss' igau 'Etat' ietat ;
  1325. ERRE 5 ;
  1326. FINSI ;
  1327.  
  1328. SI (EGA ierr 30) ;
  1329. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEIN icomp)
  1330. ' Point Gauss' igau 'Etat' ietat ;
  1331. ERRE 5 ;
  1332. FINSI ;
  1333.  
  1334. SI (EGA ierr 31) ;
  1335. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEIN icomp)
  1336. ' Point Gauss' igau 'Etat' ietat ;
  1337. ERRE 5 ;
  1338. FINSI ;
  1339.  
  1340. SI (EGA ierr 40) ;
  1341. Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  1342. ' Point Gauss' igau 'Etat' ietat ;
  1343. ERRE 5 ;
  1344. FINSI ;
  1345.  
  1346. SI (EGA ierr 41) ;
  1347. Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  1348. ' Point Gauss' igau 'Etat' ietat ;
  1349. ERRE 5 ;
  1350. FINSI ;
  1351.  
  1352. SI (EGA ierr 50) ;
  1353. Message 'Erreur absolue' errabs 'Composante ' (EXTR LVARI icomp)
  1354. ' Point Gauss' igau 'Etat' ietat ;
  1355. ERRE 5 ;
  1356. FINSI ;
  1357.  
  1358. SI (EGA ierr 51) ;
  1359. Message 'Erreur relative' errrel 'Composante ' (EXTR LVARI icomp)
  1360. ' Point Gauss' igau 'Etat' ietat ;
  1361. ERRE 5 ;
  1362. FINSI ;
  1363.  
  1364. SINON ;
  1365.  
  1366. Message 'Resultats concordants entre les deux modeles' ;
  1367. Message 'Aucune erreur relative superieure a ' critere
  1368. ' pour les composantes' ;
  1369. Message (EXTR LDEFO idefo) ' ' (EXTR LDEIN idein)
  1370. ' ' (EXTR LCONT icont) ;
  1371. ERRE 0 ;
  1372.  
  1373. FINSI ;
  1374.  
  1375. FIN ;
  1376.  
  1377.  
  1378.  
  1379.  
  1380.  
  1381.  
  1382.  
  1383.  
  1384.  
  1385.  
  1386.  
  1387.  

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