Télécharger la plaquette Cast3M


fr en

Télécharger isotro_cis.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  isotro_cis.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **************************************************
  5. *                                                *
  6. *  TEST DE VALIDATION D'UNE LOI DE COMPORTEMENT  *
  7. *  --------------------------------------------  *
  8. *              DE MATERIAU                       *
  9. *              -----------                       *
  10. *                                                *
  11. *  LOI DE COMPORTEMENT ELASTO-PLASTIQUE          *
  12. *             ISOTROPE                           *
  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. graph = 'N' ;
  28.  
  29. opti debu 1 ;
  30.  
  31. opti echo 1 dime 3 elem cub8 ;
  32. *  
  33. *  
  34. *  Maillage 
  35. *  
  36. rayon0 = .05 ; 
  37. epais0 = .0005 ; 
  38. l1 = epais0 + rayon0 ; 
  39. l2 = .1 ; 
  40. n1 = 1 ; 
  41. n2 = 1 ; 
  42. p1 = rayon0 0. 0. ; 
  43. p2 = l1 0. 0. ; 
  44. p3 = l1 l2 0. ; 
  45. p4 = rayon0 l2 0. ; 
  46. p11 = rayon0 0. l2 ;
  47. p21 = l1 0. l2 ;
  48. p31 = l1 l2 l2 ;
  49. p41 = rayon0 l2 l2 ;
  50. d1 = d p1 n1 p2 ; 
  51. d2 = d p2 n2 p3 ; 
  52. d3 = d p3 n1 p4 ; 
  53. d4 = d p4 n2 p1 ; 
  54. d11 = d p11 n1 p21 ;
  55. d21 = d p21 n2 p31 ;
  56. d31 = d p31 n1 p41 ;
  57. d41 = d p41 n2 p11 ; 
  58. *p01 = ( rayon0 + ( epais0 / 2. ) ) 0. l2 ;
  59. *p02 = ( rayon0 + ( epais0 / 2. ) ) l2 l2 ;
  60. sur1 = 'DALL' d1 d2 d3 d4 ; 
  61. sur2 = 'DALL' d11 d21 d31 d41 ; 
  62. vol1 = sur1 'VOLU' n2 sur2 ; 
  63. *titr ' maillage ' ; 
  64. *trace vol1 ; 
  65. *  
  66. *  Deplacement maximal
  67. *
  68. depmax = .01 ;
  69. *
  70. *  Courbe d'ecrouissage
  71. *
  72. beta0 = 1250. ;
  73. sig0 = 410.E6 ;
  74. xm = .5 ;
  75. youn0 = 7.34E10 ;
  76. epsp0 = 'PROG' 0. ;
  77. deb0 = 'PROG' 0. ;
  78. epsmax0 = depmax / l2 ;
  79. pas0 = epsmax0 / 50 ;
  80. epsp1 = 'PROG' 0. 'PAS' pas0 epsmax0 ;
  81. dim00 = 'DIME' epsp1 ;
  82. un = 'PROG' dim00 '*' 1. ;
  83. Y0 = sig0 ; 
  84. epsp2 = beta0 * epsp1 ;
  85. epsp3 = epsp2 + un ;
  86. epsp31 = log epsp3 ;
  87. epsp32 = epsp31 * xm ;
  88. epsp4 = exp epsp32 ;
  89. Y1 = Y0 * epsp4 ;
  90. *
  91. epsp1 = epsp0 'ET' epsp1 ;
  92. Y1 = deb0 'ET' Y1 ;
  93. epsp2 = epsp1 + ( Y1 / youn0 ) ;
  94. trac00 = 'EVOL' 'MANU' 'defo' epsp1 'sigm' Y1 ;
  95. * 
  96. *    Modele de calcul 
  97. * 
  98. trac0 = 'EVOL' 'MANU' 'defo' epsp2 'sigm' Y1 ; 
  99. y1    = y1 enle 1 ;
  100. epsp2 = (epsp2 enle 1) - (y1 / youn0) ;
  101. ecro0 = evol vert manu ep epsp2 sig y1 ;
  102. si (neg graph 'N') ;
  103.   'DESS' (trac0 et ecro0) titr ' Courbes de traction et d ecrouissage (vert)' ; 
  104. fins ;
  105.  
  106. mod0 = 'MODE' vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'   
  107.        'PLASTIQUE' 'ISOTROPE' ; 
  108. mat0 = 'MATE' mod0 'YOUN' youn0 'NU' 0.33   
  109.        'ECRO' ecro0 ;
  110. *  
  111. *    Conditions aux limites 
  112. *  
  113. cl1 = 'BLOQ' sur1 'UX' 'RY' ; 
  114. cl2 = 'BLOQ' sur2 'UX' 'RY' ;  
  115. cl3 = 'BLOQ' vol1 'UY' 'UZ' ; 
  116. cl0 = cl1 'ET' cl2 'ET' cl3 ;  
  117. *  
  118. *    Chargement 
  119. *  
  120. chp1 = 'COOR' 3 vol1 ;
  121. chp1 = chp1 * depmax / l2 ;
  122. chp1 = 'NOMC' chp1 'UY' ;
  123. chp2 = 'MANU' 'CHPO' vol1 1 'UZ' 0. ;
  124. *trace chp1 vol1 ;
  125. chp0 = chp1 + chp2 ; 
  126. dep1 = 'DEPI' cl3 chp0 ;  
  127. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y ( 'PROG' 0. 1. ) ; 
  128. cha0 = 'CHAR' 'DIMP' dep1 ev0 ; 
  129. *vec1 = 'VECT' dep1 1.E-5 'FX' 'FY' 'FZ' 'ROUG' ; 
  130. *titr ' chargement ' ; 
  131. *trace vec1 vol1 ; 
  132. *  
  133. *  Temps du calcul  
  134. *  
  135. dt0 = 50. ; 
  136. tfin0 = 1000. ;  
  137. *   
  138. *   
  139. *  Resolution par PASAPAS  
  140. *---------------------------------------
  141. *    
  142. ta1 = 'TABLE' ; 
  143. ta1 .'MODELE' = mod0 ; 
  144. ta1 .'CARACTERISTIQUES' = mat0 ; 
  145. ta1 .'BLOCAGES_MECANIQUES' = cl0 ; 
  146. ta1 .'CHARGEMENT' = cha0 ; 
  147. *ta1 .'TEMPS_CALCULES' = ( ( 'PROG' 0. 'PAS' dt0 80. ) 'ET'
  148. *( 'PROG' 85. 'PAS' 5. 90. ) 'ET' 
  149. *( 'PROG' 92.5 'PAS' 2.5 tfin0 ) ) ; 
  150. ta1 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt0 tfin0 ; 
  151. * 
  152. TMASAU=table;
  153. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  154. TMASAU .'DEFTO'=VRAI;
  155. TMASAU .'DEFIN'=VRAI; 
  156. ta1 .'HYPOTHESE_DEFORMATIONS'='LINEAIRE' ;
  157. PASAPAS ta1 ; 
  158. *  
  159. *  
  160. *  Post-traitement 
  161. *-----------------------------------------
  162. * 
  163. dim0 = 'DIME' ( ta1 .'TEMPS' ) ; 
  164. char2 = 'PROG' 0. ; 
  165. depl2 = 'PROG' 0. ; 
  166. defi2 = 'PROG' 0. ; 
  167. cont1 = 'PROG' 0. ; 
  168. *  
  169. reac0 = ta1 .'REACTIONS'.( dim0 - 1 ) ;  
  170. depl0 = ta1 .'DEPLACEMENTS'.( dim0 - 1 ) ;  
  171. vec2 = 'VECT' reac0 1.E-2 'FX' 'FY' 'FZ' 'ROUG' ; 
  172. *titr ' reactions ' ; 
  173. *trace vec2 vol1 ;
  174. defo0 = 'DEFO' vol1 depl0 0. 'ROUG' ; 
  175. defo1 = 'DEFO' vol1 depl0 1. ; 
  176. *titr ' deformees ' ; 
  177. *trace ( defo0 'ET' defo1 ) ; 
  178. *  
  179. 'REPETER' bloc0 ( dim0 - 1 ) ; 
  180.   i0 = &bloc0 ; 
  181.   reac0 = ta1 .'REACTIONS'.i0 ; 
  182.   depl0 = ta1 .'DEPLACEMENTS'.i0 ; 
  183.   sigm0 = ta1 .'CONTRAINTES'.i0 ; 
  184.   def0 = ta1 .'DEFORMATIONS_INELASTIQUES'.i0 ; 
  185.   eps0 = 'EPSI' mod0 depl0 'LINE' ;
  186.   temp0 = ta1 .'TEMPS'.i0 ;   
  187.   reac0 = 'REDU' reac0 sur2 ; 
  188.   reac0 = 'EXCO' reac0 'FY' 'SCAL' ; 
  189.   chp1 = 'MANU' 'CHPO' sur2 1 'SCAL' 1. ; 
  190.   reac1 = 'XTY' reac0 chp1 ( 'MOTS' 'SCAL' ) ( 'MOTS' 'SCAL' ) ; 
  191.   depl1 = 'EXTR' depl0 'UY' p21 ; 
  192.   depl2 = depl2 'ET' ( 'PROG' reac1 ) ; 
  193.   char2 = char2 'ET' ( 'PROG' depl1 ) ; 
  194. 'FIN' bloc0 ; 
  195. * 
  196. *   
  197. opti echo 1 dime 3 elem tri3 ; 
  198. *  
  199. *  
  200. *  Maillage 
  201. *  
  202. epais0 = .0005 ; 
  203. rayon0 = .05 ; 
  204. l10 = .1 ; 
  205. n10 = 1 ; 
  206. p10 = rayon0 0. 0. ; 
  207. p20 = rayon0 l10 0. ; 
  208. p11 = rayon0 0. l10 ; 
  209. p12 = rayon0 l10 l10 ; 
  210. li1 = d p10 n10 p20 ; 
  211. li2 = d p20 n10 p12 ; 
  212. li3 = d p12 n10 p11 ; 
  213. li4 = d p11 n10 p10 ; 
  214. mail1 = 'SURF' ( li1 'ET' li2 'ET' li3 'ET' li4 ) 'PLAN' ; 
  215. *titr 'maillage' ; 
  216. *trace mail1 ; 
  217. *  
  218. depmax = .01 ;
  219. *
  220. *  Courbe d'ecrouissage
  221. *
  222. beta0 = 1250. ;
  223. sig0 = 410.E6 ;
  224. xm = .5 ;
  225. youn0 = 7.34E10 ;
  226. epsp0 = 'PROG' 0. ;
  227. deb0 = 'PROG' 0. ;
  228. epsmax0 = depmax / l10 ;
  229. pas0 = epsmax0 / 50 ;
  230. epsp1 = 'PROG' 0. 'PAS' pas0 epsmax0 ;
  231. dim00 = 'DIME' epsp1 ;
  232. un = 'PROG' dim00 '*' 1. ;
  233. Y0 = sig0 ; 
  234. epsp2 = beta0 * epsp1 ;
  235. epsp3 = epsp2 + un ;
  236. epsp31 = log epsp3 ;
  237. epsp32 = epsp31 * xm ;
  238. epsp4 = exp epsp32 ;
  239. Y1 = Y0 * epsp4 ;
  240. *
  241. epsp1 = epsp0 'ET' epsp1 ;
  242. Y1 = deb0 'ET' Y1 ;
  243. epsp2 = epsp1 + ( Y1 / youn0 ) ;
  244. trac00 = 'EVOL' 'MANU' 'defo' epsp1 'sigm' Y1 ; 
  245. *  
  246. *    Modele de calcul 
  247. * 
  248. trac0 = 'EVOL' 'MANU' 'defo' epsp2 'sigm' Y1 ; 
  249. *'DESS' trac0 ; 
  250. mod00 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'   
  251.        'PLASTIQUE' 'ISOTROPE' 'COQ3' ; 
  252. mat00 = 'MATE' mod00 'YOUN' youn0 'NU' 0.33   
  253.         'ECRO' ecro0 ;
  254. car00 = 'CARA' mod00 'EPAI' epais0 ;   
  255. *  
  256. *    Conditions aux limites 
  257. *  
  258. cl10 = 'BLOQ' li1 'UX' 'RY' ; 
  259. cl20 = 'BLOQ' li3 'UX' 'RY' ; 
  260. cl30 = 'BLOQ' mail1 'UY' 'UZ' ; 
  261. cl00 = cl10 'ET' cl20 'ET' cl30 ;  
  262. *  
  263. *    Chargement 
  264. *  
  265. chp1 = 'COOR' 3 mail1 ;
  266. chp1 = chp1 * depmax / l10 ;
  267. chp1 = 'NOMC' chp1 'UY' ;
  268. *trace chp1 mail1 ;
  269. chp2 = 'MANU' 'CHPO' mail1 1 'UZ' 0. ;
  270. chp0 = chp1 + chp2 ; 
  271. dep10 = 'DEPI' cl30 chp0 ;  
  272. ev00 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y ( 'PROG' 0. 1. ) ; 
  273. cha00 = 'CHAR' 'DIMP' dep10 ev00 ; 
  274. *vec1 = 'VECT' dep10 1.E-5 'FX' 'FY' 'FZ' 'ROUG' ; 
  275. *titr ' chargement ' ; 
  276. *trace vec1 mail1 ; 
  277. *  
  278. *  Temps du calcul  
  279. *  
  280. dt00 = 50. ; 
  281. tfin00 = 1000. ;  
  282. *   
  283. *   
  284. *  Resolution par PASAPAS  
  285. *---------------------------------------
  286. *    
  287. ta10 = 'TABLE' ; 
  288. ta10 .'MODELE' = mod00 ; 
  289. ta10 .'CARACTERISTIQUES' = mat00 'ET' car00 ; 
  290. ta10 .'BLOCAGES_MECANIQUES' = cl00 ; 
  291. ta10 .'CHARGEMENT' = cha00 ; 
  292. *ta10 .'TEMPS_CALCULES' = ( ( 'PROG' 0. 'PAS' dt00 80. ) 'ET'
  293. *( 'PROG' 85. 'PAS' 5. 90. )  'ET' 
  294. *( 'PROG' 92.5 'PAS' 2.5 tfin00 ) ) ; 
  295. ta10 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt00 tfin00 ; 
  296. *ta10 .'TEMPS_CALCULES' = ( ( 'PROG' 0. 'PAS' 20. 60. ) 'ET' 
  297. *( 'PROG' 62. 'PAS' 2. 70. ) 'ET' 
  298. *( 'PROG' 71. 'PAS' dt00 79. ) 'ET' 
  299. *( 'PROG' 79.5 'PAS' .5 tfin00 ) ) ; 
  300. *
  301. TMASAU=table;
  302. ta10 . 'MES_SAUVEGARDES'=TMASAU;
  303. TMASAU .'DEFTO'=VRAI;
  304. TMASAU .'DEFIN'=VRAI;  
  305. ta10 .'HYPOTHESE_DEFORMATIONS'='LINEAIRE' ;
  306. PASAPAS ta10 ; 
  307. *  
  308. *  
  309. *  Post-traitement 
  310. *-----------------------------------------
  311. * 
  312. dim10 = 'DIME' ta10 .'TEMPS' ; 
  313. depl12 = 'PROG' 0. ; 
  314. *  
  315. reac0 = ta10 .'REACTIONS'.( dim10 - 1 ) ;  
  316. depl0 = ta10 .'DEPLACEMENTS'.( dim10 - 1 ) ; 
  317. vec2 = 'VECT' reac0 1.E-2 'FX' 'FY' 'FZ' 'ROUG' ; 
  318. *titr ' reactions ' ; 
  319. *trace vec2 mail1 ;   
  320. defo0 = 'DEFO' mail1 depl0 0. ;
  321. defo1 = 'DEFO' mail1 depl0 1. 'ROUG' ;
  322. *titr ' deformees ' ;  
  323. *trace ( defo0 'ET' defo1 ) ;  
  324. *   
  325. 'REPETER' bloc0 ( dim10 - 1 ) ; 
  326.   i10 = &bloc0 ; 
  327.   reac0 = ta10 .'REACTIONS'.i10 ; 
  328.   depl10 = ta10 .'DEPLACEMENTS'.i10 ; 
  329.   temp10 = ta10 .'TEMPS'.i10 ; 
  330. *  depl11 = 'EXTR' depl10 'UY' p11 ; 
  331.   reac1 = 'REDU' reac0 li3 ; 
  332. *  list reac1 ; 
  333.   eps0 = 'EPSI' depl10 mod00 car00 'LINE' ; 
  334.   sigm0 = ta10 .'CONTRAINTES'.i10 ; 
  335.   def0 = ta10 .'DEFORMATIONS_INELASTIQUES'.i10 ; 
  336.   eps0 = 'RTENS' eps0 mod00 ( 0 1 0 ) ; 
  337.   sigm0 = 'RTENS' sigm0 mod00 ( 0 1 0 ) ; 
  338.   def0 = 'RTENS' def0 mod00 ( 0 1 0 ) ;
  339.   reac1 = 'EXCO' reac1 'FY' 'SCAL' ; 
  340.   chp1 = 'MANU' 'CHPO' li3 1 'SCAL' 1 ; 
  341.   depl11 = 'XTY' reac1 chp1 ( 'MOTS' 'SCAL' ) ( 'MOTS' 'SCAL' ) ; 
  342.   depl12 = depl12 'ET' ( 'PROG' depl11 ) ; 
  343. 'FIN' bloc0 ; 
  344. *  
  345. *
  346. opti echo 1 dime 3 elem qua8 ; 
  347. *  
  348. *  
  349. *  Maillage 
  350. *  
  351. epais0 = .0005 ; 
  352. rayon0 = .05 ; 
  353. l10 = .1 ; 
  354. n10 = 1 ; 
  355. p10 = rayon0 0. 0. ; 
  356. p20 = rayon0 l10 0. ; 
  357. p11 = rayon0 0. l10 ; 
  358. p12 = rayon0 l10 l10 ; 
  359. li1 = d p10 n10 p20 ; 
  360. li2 = d p20 n10 p12 ; 
  361. li3 = d p12 n10 p11 ; 
  362. li4 = d p11 n10 p10 ; 
  363. mail1 = 'SURF' ( li1 'ET' li2 'ET' li3 'ET' li4 ) 'PLAN' ; 
  364. *titr 'maillage' ; 
  365. *trace mail1 ; 
  366. *  
  367. depmax = .01 ;
  368. *
  369. *  Courbe d'ecrouissage
  370. *
  371. beta0 = 1250. ;
  372. sig0 = 410.E6 ;
  373. xm = .5 ;
  374. youn0 = 7.34E10 ;
  375. epsp0 = 'PROG' 0. ;
  376. deb0 = 'PROG' 0. ;
  377. epsmax0 = depmax / l10 ;
  378. pas0 = epsmax0 / 50 ;
  379. epsp1 = 'PROG' 0. 'PAS' pas0 epsmax0 ;
  380. dim00 = 'DIME' epsp1 ;
  381. un = 'PROG' dim00 '*' 1. ;
  382. Y0 = sig0 ; 
  383. epsp2 = beta0 * epsp1 ;
  384. epsp3 = epsp2 + un ;
  385. epsp31 = log epsp3 ;
  386. epsp32 = epsp31 * xm ;
  387. epsp4 = exp epsp32 ;
  388. Y1 = Y0 * epsp4 ;
  389. *
  390. epsp1 = epsp0 'ET' epsp1 ;
  391. Y1 = deb0 'ET' Y1 ;
  392. epsp2 = epsp1 + ( Y1 / youn0 ) ;
  393. trac00 = 'EVOL' 'MANU' 'defo' epsp1 'sigm' Y1 ; 
  394. *  
  395. *    Modele de calcul 
  396. * 
  397. trac0 = 'EVOL' 'MANU' 'defo' epsp2 'sigm' Y1 ; 
  398. *'DESS' trac0 ; 
  399. mod00 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'   
  400.        'PLASTIQUE' 'ISOTROPE' 'COQ8' ; 
  401. mat00 = 'MATE' mod00 'YOUN' youn0 'NU' 0.33   
  402.         'ECRO' ecro0 ;
  403. car00 = 'CARA' mod00 'EPAI' epais0 ;   
  404. *  
  405. *    Conditions aux limites 
  406. *  
  407. cl10 = 'BLOQ' li1 'RY' ; 
  408. cl20 = 'BLOQ' li3 'RY' ; 
  409. cl30 = 'BLOQ' mail1 'UX' 'UY' 'UZ' ; 
  410. cl00 = cl10 'ET' cl20 'ET' cl30 ;  
  411. *  
  412. *    Chargement 
  413. *  
  414. chp1 = 'COOR' 3 mail1 ;
  415. chp1 = chp1 * depmax / l10 ;
  416. chp1 = 'NOMC' chp1 'UY' ;
  417. *trace chp1 mail1 ;
  418. chp2 = 'MANU' 'CHPO' mail1 1 'UZ' 0. ;
  419. chp3 = 'MANU' 'CHPO' mail1 1 'UX' 0. ; 
  420. chp0 = chp1 + chp2 + chp3 ; 
  421. dep10 = 'DEPI' cl30 chp0 ;  
  422. ev00 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y ( 'PROG' 0. 1. ) ; 
  423. cha00 = 'CHAR' 'DIMP' dep10 ev00 ; 
  424. *vec1 = 'VECT' dep10 1.E-5 'FX' 'FY' 'FZ' 'ROUG' ; 
  425. *titr ' chargement ' ; 
  426. *trace vec1 mail1 ; 
  427. *  
  428. *  Temps du calcul  
  429. *  
  430. dt00 = 50. ; 
  431. tfin00 = 1000. ;  
  432. *   
  433. *   
  434. *  Resolution par PASAPAS  
  435. *---------------------------------------
  436. *    
  437. ta10 = 'TABLE' ; 
  438. ta10 .'MODELE' = mod00 ; 
  439. ta10 .'CARACTERISTIQUES' = mat00 'ET' car00 ; 
  440. ta10 .'BLOCAGES_MECANIQUES' = cl00 ; 
  441. ta10 .'CHARGEMENT' = cha00 ; 
  442. *ta10 .'TEMPS_CALCULES' = ( ( 'PROG' 0. 'PAS' dt00 80. ) 'ET' 
  443. *( 'PROG' 85. 'PAS' 5. 90. ) 'ET' 
  444. *( 'PROG' 92.5 'PAS' 2.5 tfin00 ) ) ; 
  445. ta10 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt00 tfin00 ; 
  446. *ta10 .'TEMPS_CALCULES' = ( ( 'PROG' 0. 'PAS' 20. 60. ) 'ET' 
  447. *( 'PROG' 62. 'PAS' 2. 70. ) 'ET' 
  448. *( 'PROG' 71. 'PAS' dt00 79. ) 'ET' 
  449. *( 'PROG' 79.5 'PAS' .5 tfin00 ) ) ; 
  450. *
  451. TMASAU=table;
  452. ta10 . 'MES_SAUVEGARDES'=TMASAU;
  453. TMASAU .'DEFTO'=VRAI;
  454. TMASAU .'DEFIN'=VRAI;  
  455. ta10 .'HYPOTHESE_DEFORMATIONS'='LINEAIRE' ;
  456. PASAPAS ta10 ; 
  457. *  
  458. *  
  459. *  Post-traitement 
  460. *-----------------------------------------
  461. * 
  462. dim10 = 'DIME' ta10 .'TEMPS' ; 
  463. depl120 = 'PROG' 0. ; 
  464. *  
  465. reac0 = ta10 .'REACTIONS'.( dim10 - 1 ) ;  
  466. depl0 = ta10 .'DEPLACEMENTS'.( dim10 - 1 ) ; 
  467. vec2 = 'VECT' reac0 1.E-2 'FX' 'FY' 'FZ' 'ROUG' ; 
  468. *titr ' reactions ' ; 
  469. *trace vec2 mail1 ;   
  470. defo0 = 'DEFO' mail1 depl0 0. ;
  471. defo1 = 'DEFO' mail1 depl0 1. 'ROUG' ;
  472. *titr ' deformees ' ;  
  473. *trace ( defo0 'ET' defo1 ) ;  
  474. *   
  475. 'REPETER' bloc0 ( dim10 - 1 ) ; 
  476.   i10 = &bloc0 ; 
  477.   reac0 = ta10 .'REACTIONS'.i10 ;  
  478.   depl10 = ta10 .'DEPLACEMENTS'.i10 ; 
  479.   temp10 = ta10 .'TEMPS'.i10 ; 
  480. *  depl11 = 'EXTR' depl10 'UY' p11 ; 
  481.   reac0 = 'REDU' reac0 li3 ; 
  482.   reac0 = 'EXCO' reac0 'FY' 'SCAL' ; 
  483.   chp1 = 'MANU' 'CHPO' li3 1 'SCAL' 1. ; 
  484.   depl11 = 'XTY' reac0 chp1 ( 'MOTS' 'SCAL' ) ( 'MOTS' 'SCAL' ) ; 
  485.   depl120 = depl120 'ET' ( 'PROG' depl11 ) ; 
  486.   eps0 = 'EPSI' depl10 mod00 car00 'LINE' ;
  487.   sigm0 = ta10 .'CONTRAINTES'.i10 ;
  488.   def0 = ta10 .'DEFORMATIONS_INELASTIQUES'.i10 ; 
  489. *  sigm0 = 'RTENS' sigm0 mod00 ( 0 1 0 ) ;
  490. *  sigm1 = 'CHAN' 'CHPO' sigm0 mod00 ;
  491. *  titr ' Contrainte de cisaillement ' ;
  492. *  trace ( 'EXCO' sigm1 'SMST' ) mail1 ;
  493. *  titr ' Contrainte suivant Y' ;
  494. *  trace ( 'EXCO' sigm1 'SMSS' ) mail1 ;
  495. *  titr ' Contrainte suivant Z ' ;
  496. *  trace ( 'EXCO' sigm1 'SMTT' ) mail1 ; 
  497. 'FIN' bloc0 ; 
  498. *  
  499. *  
  500. *  
  501. opti echo 1 dime 3 ;  
  502. *  
  503. *  Donnees du calcul coque  
  504. *  
  505. *  Erreur   
  506. *  
  507. * 
  508. err_e1 = 'PROG' 0. ; 
  509. err_e2 = 'PROG' 0. ; 
  510. dim10 = 'DIME' char2 ; 
  511. * 
  512. bool0 = faux ; 
  513. 'REPETER' bloc0 ( dim10 - 1 ) ;
  514.   i10 = &bloc0 ;
  515.   depl1 = 'EXTR' depl2 i10 ; 
  516.   depl11 = 'EXTR' depl12 i10 ; 
  517.   depl111 = 'EXTR' depl120 i10 ; 
  518. * 
  519. 'SI' (depl1 'NEG' 0. ) ; 
  520.   err_e0 = 'ABS' ( ( depl1 - depl11 ) / depl1 ) ; 
  521.   err_e0 = err_e0 * 100. ; 
  522.   err_e00 = 'ABS' ( ( depl1 - depl111 ) / depl1 ) ; 
  523.   err_e00 = err_e00 * 100. ; 
  524. 'SINON' ; 
  525.   err_e0 = 0. ; 
  526.   err_e00 = 0. ; 
  527. 'FINSI' ; 
  528. *  
  529. err_e1 = err_e1 'ET' ( 'PROG' err_e0 ) ; 
  530. err_e2 = err_e2 'ET' ( 'PROG' err_e00 ) ; 
  531. 'SI' ( ( err_e0 '>' 1.E-6 ) 'OU' ( err_e00 '>' 1.E-6 ) ) ; 
  532.   bool0 = vrai ; 
  533. 'FINSI' ; 
  534. *  
  535. 'FIN' bloc0 ;  
  536. *  
  537. 'SI' ( bool0 'EGA' vrai ) ; 
  538.     'ERRE' 5 ; 
  539. 'SINON' ; 
  540.     'ERRE' 0 ; 
  541. 'FINSI' ; 
  542. *
  543. *  Traces 
  544. *  
  545. evo1 = 'EVOL' 'MANU' 'depl' char2 'force' depl2 ; 
  546. evo2 = 'EVOL' 'MANU' 'depl' char2 'force' depl12 ; 
  547. evo3 = 'EVOL' 'MANU' 'chargement' char2 'erreur_%' err_e1 ; 
  548. evo4 = 'EVOL' 'MANU' 'depl' char2 'chargement' depl120 ; 
  549. evo5 = 'EVOL' 'MANU' 'chargement' char2 'erreur_%' err_e2 ; 
  550. TAB = 'TABLE' ; 
  551. TAB. 1 = 'MARQ PLUS ' ;
  552. TAB. 2 = 'MARQ CROI ' ;
  553. TAB. 3 = 'MARQ CARR ' ;
  554. TAB.'TITRE' = 'TABLE' ;
  555. TAB.'TITRE'. 1 = MOT 'MASSIF' ;
  556. TAB.'TITRE'. 2 = MOT 'COQUE' ;
  557. TAB.'TITRE'. 3 = MOT 'EPAISSE' ;
  558. *'DESS' ( evo1 'ET' evo2 'ET' evo4 )
  559. *'TITR' ' Deplacement de l extremite du cylindre'
  560. *'LEGE' TAB ;
  561. * 
  562. TAB = 'TABLE' ;
  563. TAB. 1 = 'MARQ PLUS ' ;
  564. TAB. 2 = 'MARQ CROI ' ;
  565. TAB.'TITRE' = 'TABLE' ;
  566. TAB.'TITRE'. 1 = MOT 'MINCE' ;
  567. TAB.'TITRE'. 2 = MOT 'EPAISSE' ;
  568. *'DESS' ( evo3 'ET' evo5 )
  569. *'TITR' 'Erreur sur les deplacements '
  570. *'LEGE' TAB ;
  571. *  
  572. 'FIN' ;  
  573.  
  574.  
  575.  
  576.  

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