Télécharger creep01_cisXY.dgibi

Retour à la liste

Numérotation des lignes :

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

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