Télécharger chaboche2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : chaboche2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Chaboche2.dgibi: Jeux de données *
  5. * ------------------------------------- *
  6. * *
  7. complet = vrai;
  8. *
  9. *************************************************
  10. * pour calcul complet mettre complet à : vrai;
  11. ***************************************************
  12. **************************************************
  13. * *
  14. * TEST DE VALIDATION D'UNE LOI DE COMPORTEMENT *
  15. * -------------------------------------------- *
  16. * DE MATERIAU *
  17. * ----------- *
  18. * *
  19. * LOI DE COMPORTEMENT VISCOPLASTIQUE DE: *
  20. * CHABOCHE *
  21. * *
  22. * COMPARAISON DE CALCULS SUR DES ELEMENTS: *
  23. * - COQUE EPAISSE ( MFR=5 ) *
  24. * - COQUE MINCE ( MFR=3 ) *
  25. * - MASSIF 3D ( MFR=1 ) *
  26. * *
  27. * STRUCTURE : UNE PLAQUE CARREE *
  28. * - COTE : .1 M *
  29. * - EPAISSEUR : 5.E-4 M *
  30. * *
  31. * CHARGEMENT : DEPLACEMENTS IMPOSES *
  32. * ESSAI DE TRACTION *
  33. * *
  34. **************************************************
  35. opti echo 0 dime 3 elem qua8 ;
  36. *
  37. *
  38. * Maillage
  39. *
  40. epais0 = .0005 ;
  41. rayon0 = .05 ;
  42. l10 = .1 ;
  43. n10 = 1 ;
  44. p10 = rayon0 0. 0. ;
  45. p20 = rayon0 l10 0. ;
  46. p11 = rayon0 0. l10 ;
  47. p12 = rayon0 l10 l10 ;
  48. li1 = d p10 n10 p20 ;
  49. li2 = d p20 n10 p12 ;
  50. li3 = d p12 n10 p11 ;
  51. li4 = d p11 n10 p10 ;
  52. mail1='SURF' ( li1 'ET' li2 'ET' li3 'ET' li4 ) 'PLAN';
  53. *titr 'maillage' ;
  54. *trace mail1 ;
  55. *
  56. * Modele de calcul
  57. *
  58. mod00 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  59. 'VISCOPLASTIQUE' 'CHABOCHE' 'COQ8' ;
  60. mat00 = 'MATE' mod00 'YOUN' 7.34E10 'NU' 0.33
  61. 'RHO' 7.8E3 'N' 24
  62. 'KK' 10.E6 'K0' 116.E6 'ALFK' 1.5 'ALFR' 0.35
  63. 'ALF' 2.E6 'A1' 67.5E6 'C1' 1300 'BET1' 4807.E6
  64. 'R1' 4 'A2' 80.E6 'C2' 45 'BET2' 58480.E6
  65. 'R2' 4 'PHI' 1. 'B' 12 'GAMA' 2.E-7 'M' 2
  66. 'QMAX' 455.E6 'QSTA' 200.E6 'MU' 19 'ETA' 0.50 ;
  67. car00 = 'CARA' mod00 'EPAI' epais0 ;
  68. *
  69. * Conditions aux limites
  70. *
  71. cl10 = 'BLOQ' mail1 'UX' ;
  72. cl20 = 'BLOQ' li1 'UZ' ;
  73. cl30 = 'BLOQ' li3 'UZ' ;
  74. cl40 = 'BLOQ' li4 'UY' ;
  75. cl00 = cl10 'ET' cl20 'ET' cl30 'ET' cl40 ;
  76. *
  77. * Chargement
  78. *
  79. depmax = .0050 ;
  80. dep10 = 'DEPI' cl30 depmax ;
  81. ev00 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y
  82. ( 'PROG' 0. 1. ) ;
  83. cha00 = 'CHAR' 'DIMP' dep10 ev00 ;
  84. *vec1 = 'VECT' dep10 1.E-5 'FX' 'FY' 'FZ' 'ROUG' ;
  85. *titr ' chargement ' ;
  86. *trace vec1 mail1 ;
  87. *
  88. * Temps du calcul
  89. *
  90. dt00 = 2. ;
  91. si complet;
  92. tfin00 = 200. ;
  93. sinon;
  94. tfin00 = 25.;
  95. finsi;
  96. *
  97. *
  98. * Resolution par PASAPAS
  99. *---------------------------------------
  100. *
  101. ta10 = 'TABLE' ;
  102. ta10 .'MODELE' = mod00 ;
  103. ta10 .'CARACTERISTIQUES' = mat00 'ET' car00 ;
  104. ta10 .'BLOCAGES_MECANIQUES' = cl00 ;
  105. ta10 .'CHARGEMENT' = cha00 ;
  106. ta10 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt00 tfin00 ;
  107. *( 'PROG' 22. 'PAS' 2. 40. ) 'ET'
  108. *( 'PROG' 41. 'PAS' dt00 tfin00 ) ) ;
  109. *
  110. TMASAU=table;
  111. ta10 . 'MES_SAUVEGARDES'=TMASAU;
  112. TMASAU .'DEFTO'=VRAI;
  113. TMASAU .'DEFIN'=VRAI;
  114. PASAPAS ta10 ;
  115. *
  116. *
  117. * Post-traitement
  118. *-----------------------------------------
  119. *
  120. dim10 = 'DIME' ta10 .'TEMPS' ;
  121. depl120 = 'PROG' 0. ;
  122. *
  123. reac0 = ta10 .'REACTIONS'.( dim10 - 1 ) ;
  124. depl0 = ta10 .'DEPLACEMENTS'.( dim10 - 1 ) ;
  125. vec2 = 'VECT' reac0 1.E-2 'FX' 'FY' 'FZ' 'ROUG' ;
  126. *titr ' reactions ' ;
  127. *trace vec2 mail1 ;
  128. defo0 = 'DEFO' mail1 depl0 0. ;
  129. defo1 = 'DEFO' mail1 depl0 1. 'ROUG' ;
  130. titr ' deformees ' ;
  131. *trace ( defo0 'ET' defo1 ) ;
  132. *
  133. 'REPETER' bloc0 ( dim10 - 1 ) ;
  134. i10 = &bloc0 ;
  135. reac0 = ta10 .'REACTIONS'.i10 ;
  136. depl10 = ta10 .'DEPLACEMENTS'.i10 ;
  137. temp10 = ta10 .'TEMPS'.i10 ;
  138. depl11 = 'EXTR' depl10 'UY' p11 ;
  139. sigm0 = ta10 .'CONTRAINTES'.i10 ;
  140. sigm0 = 'RTENS' sigm0 mod00 ( 0 1 0 ) car00 ;
  141. sigm1 = 'CHAN' 'CHPO' sigm0 mod00 ;
  142. 'SI' ( i10 'EGA' ( dim10 - 1 ) ) ;
  143. titr ' Contrainte de cisaillement ' ;
  144. * trace ( 'EXCO' sigm1 'SMST' ) mail1 ;
  145. titr ' Contrainte suivant Y' ;
  146. * trace ( 'EXCO' sigm1 'SMSS' ) mail1 ;
  147. titr ' Contrainte suivant Z ' ;
  148. * trace ( 'EXCO' sigm1 'SMTT' ) mail1 ;
  149. 'FINSI' ;
  150. reac0 = 'REDU' reac0 li3 ;
  151. reac0 = 'EXCO' reac0 'FZ' 'SCAL' ;
  152. chp1 = 'MANU' 'CHPO' li3 1 'SCAL' 1. ;
  153. depl11 = 'XTY' reac0 chp1 ( 'MOTS' 'SCAL' )
  154. ( 'MOTS' 'SCAL' ) ;
  155. depl120 = depl120 'ET' ( 'PROG' depl11 ) ;
  156. 'FIN' bloc0 ;
  157. *
  158. *
  159. opti echo 0 dime 3 elem tri3 ;
  160. *
  161. *
  162. * Maillage
  163. *
  164. epais0 = .0005 ;
  165. rayon0 = .05 ;
  166. l10 = .1 ;
  167. n10 = 1 ;
  168. p10 = rayon0 0. 0. ;
  169. p20 = rayon0 l10 0. ;
  170. p11 = rayon0 0. l10 ;
  171. p12 = rayon0 l10 l10 ;
  172. li1 = d p10 n10 p20 ;
  173. li2 = d p20 n10 p12 ;
  174. li3 = d p12 n10 p11 ;
  175. li4 = d p11 n10 p10 ;
  176. mail1='SURF' ( li1 'ET' li2 'ET' li3 'ET' li4 ) 'PLAN';
  177. *titr 'maillage' ;
  178. *trace mail1 ;
  179. *
  180. * Modele de calcul
  181. *
  182. mod00 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  183. 'VISCOPLASTIQUE' 'CHABOCHE' 'DKT' ;
  184. mat00 = 'MATE' mod00 'YOUN' 7.34E10 'NU' 0.33
  185. 'RHO' 7.8E3 'N' 24
  186. 'KK' 10.E6 'K0' 116.E6 'ALFK' 1.5 'ALFR' 0.35
  187. 'ALF' 2.E6 'A1' 67.5E6 'C1' 1300 'BET1' 4807.E6
  188. 'R1' 4 'A2' 80.E6 'C2' 45 'BET2' 58480.E6
  189. 'R2' 4 'PHI' 1. 'B' 12 'GAMA' 2.E-7 'M' 2
  190. 'QMAX' 455.E6 'QSTA' 200.E6 'MU' 19 'ETA' 0.50 ;
  191. car00 = 'CARA' mod00 'EPAI' epais0 'CALF' ( 2. / 3. ) ;
  192. *
  193. * Conditions aux limites
  194. *
  195. cl10 = 'BLOQ' mail1 'UX' ;
  196. cl20 = 'BLOQ' li1 'UZ' ;
  197. cl30 = 'BLOQ' li3 'UZ' ;
  198. cl40 = 'BLOQ' li4 'UY' ;
  199. cl00 = cl10 'ET' cl20 'ET' cl30 'ET' cl40 ;
  200. *
  201. * Chargement
  202. *
  203. depmax = .0050 ;
  204. dep10 = 'DEPI' cl30 depmax ;
  205. ev00 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y
  206. ( 'PROG' 0. 1. ) ;
  207. cha00 = 'CHAR' 'DIMP' dep10 ev00 ;
  208. *vec1 = 'VECT' dep10 1.E-5 'FX' 'FY' 'FZ' 'ROUG' ;
  209. *titr ' chargement ' ;
  210. *trace vec1 mail1 ;
  211. *
  212. * Temps du calcul
  213. *
  214. *
  215. *
  216. * Resolution par PASAPAS
  217. *---------------------------------------
  218. *
  219. ta10 = 'TABLE' ;
  220. ta10 .'MODELE' = mod00 ;
  221. ta10 .'CARACTERISTIQUES' = mat00 'ET' car00 ;
  222. ta10 .'BLOCAGES_MECANIQUES' = cl00 ;
  223. ta10 .'CHARGEMENT' = cha00 ;
  224. ta10 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt00 tfin00 ;
  225. *( 'PROG' 22. 'PAS' 2. 40. ) 'ET'
  226. *( 'PROG' 41. 'PAS' dt00 tfin00 ) ) ;
  227. *
  228. TMASAU=table;
  229. ta10 . 'MES_SAUVEGARDES'=TMASAU;
  230. TMASAU .'DEFTO'=VRAI;
  231. TMASAU .'DEFIN'=VRAI;
  232. PASAPAS ta10 ;
  233. *
  234. *
  235. * Post-traitement
  236. *-----------------------------------------
  237. *
  238. dim10 = 'DIME' ta10 .'TEMPS' ;
  239. depl12 = 'PROG' 0. ;
  240. var2 = 'PROG' 0. ;
  241. *
  242. reac0 = ta10 .'REACTIONS'.( dim10 - 1 ) ;
  243. depl0 = ta10 .'DEPLACEMENTS'.( dim10 - 1 ) ;
  244. vec2 = 'VECT' reac0 1.E-2 'FX' 'FY' 'FZ' 'ROUG' ;
  245. *titr ' reactions ' ;
  246. *trace vec2 mail1 ;
  247. defo0 = 'DEFO' mail1 depl0 0. ;
  248. defo1 = 'DEFO' mail1 depl0 1. 'ROUG' ;
  249. titr ' deformees ' ;
  250. *trace ( defo0 'ET' defo1 ) ;
  251. *
  252. 'REPETER' bloc0 ( dim10 - 1 ) ;
  253. i10 = &bloc0 ;
  254. reac0 = ta10 .'REACTIONS'.i10 ;
  255. depl10 = ta10 .'DEPLACEMENTS'.i10 ;
  256. temp10 = ta10 .'TEMPS'.i10 ;
  257. sigm0 = ta10 .'CONTRAINTES'.i10 ;
  258. def0 = ta10 .'DEFORMATIONS_INELASTIQUES'.i10 ;
  259. var0 = ta10. 'VARIABLES_INTERNES'.i10 ;
  260. * var1 = 'EXTR' ( 'EXCO' var0 'GPN3' ) p11 ;
  261. eps0 = 'EPSI' depl10 mod00 car00 ;
  262. sigm0 = 'RTENS' sigm0 mod00 ( 0 1 0 ) ;
  263. eps0 = 'RTENS' eps0 mod00 ( 0 1 0 ) ;
  264. sigm1 = 'CHAN' 'CHPO' sigm0 mod00 ;
  265. eps1 = 'CHAN' 'CHPO' eps0 mod00 ;
  266. def0 = 'RTENS' def0 mod00 ( 0 1 0 ) ;
  267. def1 = 'CHAN' 'CHPO' def0 mod00 ;
  268. 'SI' ( i10 'EGA' ( dim10 - 1 ) ) ;
  269. titr ' Contrainte de cisaillement ' ;
  270. * trace ( 'EXCO' sigm1 'N12' ) mail1 ;
  271. titr ' Contrainte suivant Y' ;
  272. * trace ( 'EXCO' sigm1 'N11' ) mail1 ;
  273. titr ' Contrainte suivant Z ' ;
  274. * trace ( 'EXCO' sigm1 'N22' ) mail1 ;
  275. titr ' déformation plastique de cisaillement ' ;
  276. * trace ( 'EXCO' def1 'GAST' ) mail1 ;
  277. titr ' déformation plastique suivant Y' ;
  278. * trace ( 'EXCO' def1 'EPSS' ) mail1 ;
  279. titr ' déformation plastique Z ' ;
  280. * trace ( 'EXCO' def1 'EPTT' ) mail1 ;
  281. titr ' déformation de cisaillement ' ;
  282. * trace ( 'EXCO' eps1 'GAST' ) mail1 ;
  283. titr ' déformation suivant Y' ;
  284. * trace ( 'EXCO' eps1 'EPSS' ) mail1 ;
  285. titr ' déformation suivant Z ' ;
  286. * trace ( 'EXCO' eps1 'EPTT' ) mail1 ;
  287. 'FINSI' ;
  288. depl11 = 'EXTR' depl10 'UZ' p11 ;
  289. reac0 = 'REDU' reac0 li3 ;
  290. reac0 = 'EXCO' reac0 'FZ' 'SCAL' ;
  291. chp1 = 'MANU' 'CHPO' li3 1 'SCAL' 1. ;
  292. depl11 = 'XTY' reac0 chp1 ( 'MOTS' 'SCAL' )
  293. ( 'MOTS' 'SCAL' ) ;
  294. depl12 = depl12 'ET' ( 'PROG' depl11 ) ;
  295. * var2 = var2 'ET' ( 'PROG' var1 ) ;
  296. 'FIN' bloc0 ;
  297. *
  298. *
  299. opti echo 0 dime 3 elem cub8 ;
  300. *
  301. *
  302. * Maillage
  303. *
  304. rayon0 = .05 ;
  305. epais0 = .0005 ;
  306. l1 = epais0 + rayon0 ;
  307. l2 = .1 ;
  308. n1 = 1 ;
  309. n2 = 1 ;
  310. p1 = rayon0 0. 0. ;
  311. p2 = l1 0. 0. ;
  312. p3 = l1 l2 0. ;
  313. p4 = rayon0 l2 0. ;
  314. p11 = rayon0 0. l2 ;
  315. p21 = l1 0. l2 ;
  316. p31 = l1 l2 l2 ;
  317. p41 = rayon0 l2 l2 ;
  318. d1 = d p1 n1 p2 ;
  319. d2 = d p2 n2 p3 ;
  320. d3 = d p3 n1 p4 ;
  321. d4 = d p4 n2 p1 ;
  322. d11 = d p11 n1 p21 ;
  323. d21 = d p21 n2 p31 ;
  324. d31 = d p31 n1 p41 ;
  325. d41 = d p41 n2 p11 ;
  326. p01 = ( rayon0 + ( epais0 / 2. ) ) 0. l2 ;
  327. p02 = ( rayon0 + ( epais0 / 2. ) ) l2 l2 ;
  328. sur1 = 'DALL' d1 d2 d3 d4 ;
  329. sur2 = 'DALL' d11 d21 d31 d41 ;
  330. vol1 = sur1 'VOLU' n2 sur2 ;
  331. *titr ' maillage ' ;
  332. *trace vol1 ;
  333. *
  334. * Modele de calcul
  335. *
  336. mod0 = 'MODE' vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  337. 'VISCOPLASTIQUE' 'CHABOCHE' ;
  338. mat0 = 'MATE' mod0 'YOUN' 7.34E10 'NU' 0.33
  339. 'RHO' 7.8E3 'ALPHA' 1.E-5 'N' 24
  340. 'KK' 10.E6 'K0' 116.E6 'ALFK' 1.5 'ALFR' 0.35
  341. 'ALF' 2.E6 'A1' 67.5E6 'C1' 1300 'BET1' 4807.E6
  342. 'A2' 80.E6 'C2' 45 'BET2' 58480.E6 'R2' 4
  343. 'PHI' 1. 'R1' 4 'B' 12 'GAMA' 2.E-7 'M' 2
  344. 'QMAX' 455.E6 'QSTA' 200.E6 'MU' 19 'ETA' 0.50 ;
  345. *
  346. * Conditions aux limites
  347. *
  348. cl1 = 'BLOQ' ( d2 'ET' d21 ) 'UX' ;
  349. cl2 = 'BLOQ' sur1 'UZ' ;
  350. cl3 = 'BLOQ' sur2 'UZ' ;
  351. cl4 = 'BLOQ' ( d1 'ET' d11 ) 'UY' ;
  352. cl0 = cl1 'ET' cl2 'ET' cl3 'ET' cl4 ;
  353. *
  354. * Chargement
  355. *
  356. depmax = .0050 ;
  357. dep1 = 'DEPI' cl3 depmax ;
  358. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y
  359. ( 'PROG' 0. 1. ) ;
  360. cha0 = 'CHAR' 'DIMP' dep1 ev0 ;
  361.  
  362. *
  363. *
  364. * Resolution par PASAPAS
  365. *---------------------------------------
  366. *
  367. ta1 = 'TABLE' ;
  368. ta1 .'MODELE' = mod0 ;
  369. ta1 .'CARACTERISTIQUES' = mat0 ;
  370. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  371. ta1 .'CHARGEMENT' = cha0 ;
  372. ta1 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt00 tfin00 ;
  373. *( 'PROG' 22. 'PAS' 2. 40. ) 'ET'
  374. *( 'PROG' 41. 'PAS' dt00 tfin0 ) ) ;
  375. *
  376. TMASAU=table;
  377. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  378. TMASAU .'DEFTO'=VRAI;
  379. TMASAU .'DEFIN'=VRAI;
  380. PASAPAS ta1 ;
  381. *
  382. *
  383. * Post-traitement
  384. *-----------------------------------------
  385. *
  386. dim0 = 'DIME' ( ta1 .'TEMPS' ) ;
  387. char2 = 'PROG' 0. ;
  388. depl2 = 'PROG' 0. ;
  389. defi2 = 'PROG' 0. ;
  390. var10 = 'PROG' 0. ;
  391. *
  392. reac0 = ta1 .'REACTIONS'.( dim0 - 1 ) ;
  393. depl0 = ta1 .'DEPLACEMENTS'.( dim0 - 1 ) ;
  394. vec2 = 'VECT' reac0 1.E-2 'FX' 'FY' 'FZ' 'ROUG' ;
  395. *titr ' reactions ' ;
  396. *trace vec2 vol1 ;
  397. defo0 = 'DEFO' vol1 depl0 0. 'ROUG' ;
  398. defo1 = 'DEFO' vol1 depl0 1. ;
  399. titr ' deformees ' ;
  400. *trace ( defo0 'ET' defo1 ) ;
  401. reac1 = 'REDU' reac0 p1 ;
  402. reac2 = 'REDU' reac0 p2 ;
  403. vec01 = 'VECT' reac1 1.E-2 'FX' 'FY' 'FZ' 'ROUG' ;
  404. vec02 = 'VECT' reac2 1.E-2 'FX' 'FY' 'FZ' 'VERT' ;
  405. *titr ' reactions 2 ' ;
  406. *trace ( vec01 'ET' vec02 ) vol1 ;
  407. *
  408. 'REPETER' bloc0 ( dim0 - 1 ) ;
  409. i0 = &bloc0 ;
  410. reac0 = ta1 .'REACTIONS'.i0 ;
  411. depl0 = ta1 .'DEPLACEMENTS'.i0 ;
  412. sigm0 = ta1 .'CONTRAINTES'.i0 ;
  413. temp0 = ta1 .'TEMPS'.i0 ;
  414. def0 = ta1 .'DEFORMATIONS_INELASTIQUES'.i0 ;
  415. char1 = 'EXTR' depl0 'UZ' p21 ;
  416. reac0 = 'REDU' reac0 sur2 ;
  417. reac0 = 'EXCO' reac0 'FZ' 'SCAL' ;
  418. chp1 = 'MANU' 'CHPO' sur2 1 'SCAL' 1. ;
  419. depl11 = 'XTY' reac0 chp1 ( 'MOTS' 'SCAL' )
  420. ( 'MOTS' 'SCAL' ) ;
  421. depl2 = depl2 'ET' ( 'PROG' depl11 ) ;
  422. sigm0 = 'CHAN' 'CHPO' sigm0 mod0 ;
  423. def1 = 'CHAN' 'CHPO' def0 mod0 ;
  424. var0 = ta1 .'VARIABLES_INTERNES'.i0 ;
  425. * var1 = 'EXTR' ( 'EXCO' var0 'GPYZ' ) p21 ;
  426. 'SI' ( i0 'EGA' ( dim0 - 1 ) ) ;
  427. titr ' déformation plastique de cisaillement ' ;
  428. * trace ( 'EXCO' def1 'GAYZ' ) vol1 ;
  429. titr ' déformation plastique suivant Y' ;
  430. * trace ( 'EXCO' def1 'EPYY' ) vol1 ;
  431. titr ' déformation plastique Z ' ;
  432. * trace ( 'EXCO' def1 'EPZZ' ) vol1 ;
  433. titr ' déformation plastique suivant X' ;
  434. * trace ( 'EXCO' def1 'EPXX' ) vol1 ;
  435. titr ' Contrainte de cisaillement' ;
  436. * trace ( 'EXCO' sigm0 'SMYZ' ) vol1 ;
  437. titr ' Contrainte suivant Y' ;
  438. * trace ( 'EXCO' sigm0 'SMYY' ) vol1 ;
  439. titr ' Contrainte suivant Z' ;
  440. * trace ( 'EXCO' sigm0 'SMZZ' ) vol1 ;
  441. titr ' Contrainte suivant X' ;
  442. * trace ( 'EXCO' sigm0 'SMXX' ) vol1 ;
  443. 'FINSI' ;
  444. * cont0 = 'MAXI' ( 'EXCO' 'SMXX' sigm0 ) ;
  445. * cont1 = cont1 'ET' ( 'PROG' cont0 ) ;
  446. char2 = char2 'ET' ( 'PROG' char1 ) ;
  447. * var10 = var10 'ET' ( 'PROG' var1 ) ;
  448. 'FIN' bloc0 ;
  449. *
  450. *
  451. opti echo 1 dime 3 ;
  452. *
  453. *
  454. * Données du calcul coque
  455. *
  456. * Erreur
  457. *
  458. *
  459. err_e1 = 'PROG' 0. ;
  460. err_e2 = 'PROG' 0. ;
  461. dim10 = 'DIME' char2 ;
  462. *
  463. bool0 = faux ;
  464. 'REPETER' bloc0 ( dim10 - 1 ) ;
  465. i10 = &bloc0 ;
  466. depl1 = 'EXTR' depl2 i10 ;
  467. depl11 = 'EXTR' depl12 i10 ;
  468. depl111 = 'EXTR' depl120 i10 ;
  469. *
  470. 'SI' (depl1 'NEG' 0. ) ;
  471. err_e0 = 'ABS' ( ( depl1 - depl11 ) / depl1 ) ;
  472. err_e0 = err_e0 * 100. ;
  473. err_e00 = 'ABS' ( ( depl1 - depl111 ) / depl1 ) ;
  474. err_e00 = err_e00 * 100. ;
  475. 'SINON' ;
  476. err_e0 = 0. ;
  477. err_e00 = 0. ;
  478. 'FINSI' ;
  479. *
  480. 'SI' ( ( err_e0 '>' 5.E-2 ) 'OU'
  481. ( err_e00 '>' 5.E-2 ) ) ;
  482. bool0 = vrai ;
  483. 'FINSI' ;
  484. err_e1 = err_e1 'ET' ( 'PROG' err_e0 ) ;
  485. err_e2 = err_e2 'ET' ( 'PROG' err_e00 ) ;
  486. *
  487. 'FIN' bloc0 ;
  488. *
  489. 'SI' ( bool0 'EGA' vrai ) ;
  490. 'ERRE' 5 ;
  491. 'SINON' ;
  492. 'ERRE' 0 ;
  493. 'FINSI' ;
  494. *
  495. * Traces
  496. *
  497. evo1 = 'EVOL' 'MANU' 'depl' char2 'force' depl2 ;
  498. evo2 = 'EVOL' 'MANU' 'depl' char2 'force' depl12 ;
  499. evo3 = 'EVOL' 'MANU' 'chargement' char2 'erreur_%'
  500. err_e1 ;
  501. evo4 = 'EVOL' 'MANU' 'depl' char2 'chargement' depl120;
  502. evo5 = 'EVOL' 'MANU' 'chargement' char2 'erreur_%'
  503. err_e2 ;
  504. TAB = 'TABLE' ;
  505. TAB. 1 = 'MARQ PLUS ' ;
  506. TAB. 2 = 'MARQ CROI ' ;
  507. TAB. 3 = 'MARQ CARR ' ;
  508. TAB.'TITRE'= 'TABLE' ;
  509. TAB.'TITRE'. 1 = MOT 'MASSIF' ;
  510. TAB.'TITRE'. 2 = MOT 'COQUE' ;
  511. TAB.'TITRE'. 3 = MOT 'EPAISSE' ;
  512. *'DESS' ( evo1 'ET' evo2 'ET' evo4 )
  513. *'TITR' ' Deplacement de l extremite du cylindre'
  514. *'LEGE' TAB ;
  515. *
  516. TAB = 'TABLE' ;
  517. TAB. 1 = 'MARQ PLUS ' ;
  518. TAB. 2 = 'MARQ CROI ' ;
  519. TAB.'TITRE' = 'TABLE' ;
  520. TAB.'TITRE'. 1 = MOT 'MINCE' ;
  521. TAB.'TITRE'. 2 = MOT 'EPAISSE' ;
  522. *'DESS' ( evo3 'ET' evo5 )
  523. *'TITR' 'Erreur sur les deplacements '
  524. *'LEGE' TAB ;
  525. *
  526. *
  527. 'FIN' ;
  528.  
  529.  
  530.  
  531.  
  532.  
  533.  
  534.  
  535.  
  536.  
  537.  
  538.  
  539.  
  540.  
  541.  
  542.  
  543.  
  544.  
  545.  
  546.  

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