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

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