Télécharger la plaquette Cast3M


fr en

Télécharger ohno_cis1.dgibi

Retour à la liste

Numérotation des lignes : 

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

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