Télécharger la plaquette Cast3M


fr en

Télécharger onera2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  onera2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *        Test onera2.dgibi: Jeux de données        *
  5. *        ----------------------------------        *
  6. *                                                  *
  7. complet = faux;
  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. *          ONERA (CHABOCHE unifie)               *
  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' 'ONERA' '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' 'ONERA' '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' 'ONERA' ; 
  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.  
  547.  
  548.  

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