Télécharger creep01_cisXZ.dgibi

Retour à la liste

Numérotation des lignes :

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

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