Télécharger creep01_cisYZ.dgibi

Retour à la liste

Numérotation des lignes :

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

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