Télécharger poudre6.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : poudre6.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. ******************************************************************
  6. * *
  7. * Cylindre de poudre *
  8. * *
  9. * Programme test numero 6 de la loi poudre_A - le 14/03/1997 *
  10. * *
  11. * Laurent Sanchez *
  12. * *
  13. ******************************************************************
  14. * loi a 3 parametres log A,Vitesse de deplacement imposee
  15. * Un cylindre de poudre seul en TA6V est densifie.
  16. * Simulation d'un essai en matrice fermee
  17. * VITESSE DE DEPLACEMENT IMPOSEE EGALE A 10-3MM/S - RHO0 = 0.7
  18. * La densite initiale de la poudre est uniforme.
  19. * Calcul avec Gdes Rotations, Gdes deformations, Pression
  20. * Suiveuse, log A & les vrais coefficients thermiques, New E
  21. * ATTENTION : PRG fonctionnant avec k2001
  22. * utilation de POUDRE_A : loi a 3 parametres
  23. * Le maillage est regulier en Q8.
  24. *
  25. OPTION DIME 2 ELEM QUA8 MODE AXIS ;
  26. *
  27. *
  28. ******************************************************************
  29. * *
  30. * Maillage *
  31. * *
  32. ******************************************************************
  33. *
  34. *
  35. * Parametres du maillage en mm
  36. *
  37. rayon_1 = 2.0 ;
  38. haute_1 = 2.0 ;
  39. *
  40. *
  41. * Definition des coordonnees des points
  42. *
  43. xa00 = 0.0 ;
  44. ya00 = 0.0 ;
  45. *
  46. xa01 = rayon_1 ;
  47. ya01 = ya00 ;
  48. *
  49. xa02 = xa01 ;
  50. ya02 = ya01 + haute_1 ;
  51. *
  52. xa03 = 0.0 ;
  53. ya03 = ya02 ;
  54. *
  55. *
  56. * Definition des points
  57. *
  58. a00 = xa00 ya00 ;
  59. a01 = xa01 ya01 ;
  60. a02 = xa02 ya02 ;
  61. a03 = xa03 ya03 ;
  62. *
  63. *
  64. * Definition du nombre d'elements
  65. *
  66. nba00a01 = 1 ;
  67. nba01a02 = 1 ;
  68. nba02a03 = 1 ;
  69. nba03a00 = 1 ;
  70. *
  71. *
  72. * Definition des lignes
  73. *
  74. la00a01 = DROI nba00a01 a00 a01 ;
  75. la01a02 = DROI nba01a02 a01 a02 ;
  76. la02a03 = DROI nba02a03 a02 a03 ;
  77. la03a00 = DROI nba03a00 a03 a00 ;
  78. *
  79. *
  80. * Definition des surfaces
  81. *
  82. elmat1 = SURF ( la00a01 ET la01a02 ET la02a03 ET la03a00 )
  83. PLAN COUL bleu ;
  84. *
  85. piece_1 = elmat1 ;
  86. tot = piece_1 ;
  87. *
  88. *
  89. ******************************************************************
  90. * *
  91. * Contour exterieur et symetrie *
  92. * *
  93. ******************************************************************
  94. *
  95. *
  96. l_ext1 = la01a02 ;
  97. l_ext2 = la02a03 ;
  98. l_ext = l_ext1 ET l_ext2 ;
  99. *
  100. l_symZ = la00a01 ;
  101. l_symR = la03a00 ;
  102. *
  103. *TRAC l_ext ;
  104. *TRAC l_sym ;
  105. TITR 'Maillage regulier en Q8 d un quart de piece: tot' ;
  106. *TRAC tot ;
  107. *
  108. *
  109. *****************************************************************
  110. * *
  111. * Conditions de l'essai *
  112. * *
  113. *****************************************************************
  114. *
  115. *
  116. T0tempe = 800. ;
  117. T1tempe = 800. ;
  118. T2tempe = 800. ;
  119. T3tempe = 800. ;
  120. *
  121. T0temps = 0. ;
  122. T1temps = 600. ;
  123. T2temps = T1temps + 600. ;
  124. T3temps = T2temps + 600. ;
  125. *
  126. *
  127. *****************************************************************
  128. * *
  129. * Champ-point de temperature *
  130. * *
  131. *****************************************************************
  132. *
  133. *
  134. T0CHTEM = MANU 'CHPO' tot 1 'T' T0tempe ;
  135. *
  136. clth0001 = BLOQ 'T' l_ext ;
  137. p_bloth1 = DEPI clth0001 1. ;
  138. p_temps1 = PROG T0temps T1temps T2temps T3temps ;
  139. p_theta1 = PROG T0tempe T1tempe T2tempe T3tempe ;
  140. ev_ther1 = EVOL MANU TEMPS p_temps1 'T' p_theta1 ;
  141. *DESS ev_ther1 TITR 'Evolution du chargement en temperature (C)' ;
  142. ch_ther1 = CHAR 'TIMP' p_bloth1 ev_ther1 ;
  143. *
  144. *
  145. *****************************************************************
  146. * *
  147. * Fichier Materiau Mecanique - Poudre TA6V *
  148. * *
  149. *****************************************************************
  150. *
  151. *
  152. p1model = MODE piece_1 MECANIQUE ELASTIQUE VISCOPLASTIQUE
  153. POUDRE_A CONS DEUX ;
  154. *
  155. V1YOUN = 76.0e+03 ;
  156. V1NU = 0.33 ;
  157. V1ALPHA = 1.06e-5 ;
  158. *
  159. V1logA = -9.28 ;
  160. V1N = 3.04 ;
  161. V1QSRT = 0. ;
  162. *
  163. V1RHOR = 0.7 ;
  164. V1F0 = 0.1627 ;
  165. V1F1 = 0. ;
  166. V1F2 = 0.8743 ;
  167. V1F3 = 0. ;
  168. V1F4 = 2. ;
  169. V1F5 = 0.66 ;
  170. V1C0 = 1.9325 ;
  171. V1C1 = 0. ;
  172. V1C2 = 0.6013 ;
  173. V1C3 = 0. ;
  174. V1C4 = 0. ;
  175. V1C5 = 1000. ;
  176. *
  177. *
  178. p1mater = MATER p1model YOUN V1YOUN NU V1NU ALPHA V1ALPHA
  179. A V1logA N V1N QSRT V1QSRT F0 V1F0 F1 V1F1 F2 V1F2 F3 V1F3
  180. F4 V1F4 F5 V1F5 C0 V1C0 C1 V1C1 C2 V1C2 C3 V1C3 C4 V1C4
  181. C5 V1C5 RHOR V1RHOR ;
  182. *
  183. *
  184. *****************************************************************
  185. * *
  186. * Fichier Materiau Thermique - Poudre TA6V *
  187. * *
  188. *****************************************************************
  189. *
  190. *
  191. p2model = MODE piece_1 THERMIQUE ISOTROPE CONS UN ;
  192. *
  193. V2C = 845. ;
  194. V2K = 18.0e-3 ;
  195. V2RHO = 4.0e-6 ;
  196. *
  197. *
  198. p2mater = MATER p2model 'C' V2C 'K' V2K 'RHO' V2RHO ;
  199. *
  200. *
  201. *****************************************************************
  202. * *
  203. * Chargement *
  204. * *
  205. *****************************************************************
  206. *
  207. *
  208. p0bloca = BLOQ UZ l_symZ ;
  209. p1bloca = BLOQ UR ( l_symR ET l_ext1 ) ;
  210. p2bloca = BLOQ UZ l_ext2 ;
  211. *p_bloca = p0bloca ET p1bloca ;
  212. p_bloca = p0bloca ET p1bloca ET p2bloca ;
  213. *
  214. *T0pres = 0.1 ;
  215. *T1pres = 15. ;
  216. *T2pres = 15. ;
  217. *T3pres = 0.1 ;
  218. *
  219. li_temps = PROG T0temps T1temps T2temps T3temps ;
  220. *
  221. *li_pres = PROG T0pres T1pres T2pres T3pres ;
  222. *ev_pres = EVOL MANU TEMPS li_temps PRES li_pres ;
  223. *DESS ev_pres TITR 'Evolution du chargement en pression (MPa)' ;
  224. *p_force1 = MANU 'CHPO' l_ext 1 'P' 1.0 ;
  225. *p_charg1 = CHAR 'PSUI' p_force1 ev_pres ;
  226. *
  227. T0depla = 0.0 ;
  228. T1depla = -0.6 ;
  229. T2depla = -0.6 ;
  230. T3depla = 0.6 ;
  231. *
  232. li_depl1 = PROG T0depla T1depla T2depla T3depla ;
  233. ev_depl1 = EVOL MANU 'TEMPS' li_temps 'DEPLACEMENT' li_depl1 ;
  234. ch_depl1 = DEPIMP p2bloca 1.0 ;
  235. p_charg1 = CHAR DIMP ch_depl1 ev_depl1 ;
  236. *
  237. *
  238. *****************************************************************
  239. * *
  240. * Calcul *
  241. * *
  242. *****************************************************************
  243. *
  244. *
  245. p_model = p1model ET p2model ;
  246. p_mater = p1mater ET p2mater ;
  247. p_charg = p_charg1 ET ch_ther1 ;
  248. *
  249. *
  250. TAB1 = TABLE ;
  251. TAB1 . 'TEMPERATURES' = TABLE ;
  252. TAB1 . 'TEMPERATURES' . 0 = T0CHTEM ;
  253. TAB1 . 'CHARGEMENT' = p_charg ;
  254. TAB1 . 'MODELE' = p_model ;
  255. TAB1 . 'CARACTERISTIQUES' = p_mater ;
  256. TAB1 . 'BLOCAGES_MECANIQUES' = p_bloca ;
  257. TAB1 . 'BLOCAGES_THERMIQUES' = clth0001 ;
  258. *
  259. TAB1 . 'TEMPS_CALCULES' = PROG 0. 50. PAS 50. 300. ;
  260. *
  261. TAB1 . 'TEMPS_SAUVES' = PROG 0. 300. ;
  262. *
  263. TAB1 . 'MAXITERATION' = 169 ;
  264. TAB1 . 'MAXISOUSPAS' = 2000 ;
  265. TAB1 . 'MOVA' = 'MOT' 'VHOR' ;
  266. TAB1 . 'PRECISION' = 5.0E-4 ;
  267. TAB1 . 'TALPHA_REFERENCE' = 25. ;
  268. TAB1 . 'PROCEDURE_THERMIQUE' = 'MOT' 'LINEAIRE' ;
  269. TAB1 . 'GRANDS_DEPLACEMENTS' = VRAI ;
  270. TAB1 . 'CELSIUS' = VRAI ;
  271. TAB1.'LAGRANGIEN' = 'MOT' 'REACTUALISE';
  272. *
  273. PASAPAS TAB1 ;
  274. *
  275. *
  276. *****************************************************************
  277. * *
  278. * Test de validation *
  279. * *
  280. *****************************************************************
  281. *
  282. *
  283. ind1 = INDEX (TAB1 . 'VARIABLES_INTERNES') ;
  284. nb1 = DIME ind1 ;
  285. *
  286. *
  287. * Densite finale numerique
  288. *
  289. vm1 = TAB1 . 'VARIABLES_INTERNES' .(ind1.nb1) ;
  290. vm2 = EXCO VHOR vm1 ;
  291. vm3 = REDU vm2 tot ;
  292. rho1 = MAXI vm3 ;
  293. *
  294. *
  295. * Densite finale analytique
  296. *
  297. rho2 = 0.85 ;
  298. *
  299. *
  300. * Contrainte radiale finale numerique
  301. *
  302. ind1 = INDEX (TAB1 . 'CONTRAINTES' ) ;
  303. nb1 = DIME ind1 ;
  304. *
  305. *
  306. * Deplacement numerique du point C
  307. *
  308. ind2 = INDEX (TAB1 . 'DEPLACEMENTS' ) ;
  309. nb2 = DIME ind2 ;
  310. uz1 = EXTR TAB1 .'DEPLACEMENTS'.(ind2.nb2) UZ a02 ;
  311. *
  312. *
  313. * Deplacement analytique du point C
  314. *
  315. uz2 = haute_1 * (V1RHOR - rho2) ;
  316. *
  317. *
  318. * Erreurs sur le deplacement et la densite :
  319. *
  320. erho1 = ABS ((rho1 - rho2 ) / rho2);
  321. euz1 = ABS ((uz1 - uz2) / uz2) ;
  322. *
  323. *
  324. * Ecriture des erreurs a l'ecran :
  325. *
  326. OPTI ECHO 0 ;
  327. *
  328. SAUTER 3 LIGNES ;
  329. *
  330. MESS 'Calcul en matrice fermee :' ;
  331. MESS 'Erreur sur le deplacement vertical : ' euz1 ;
  332. MESS 'Erreur sur la densite finale : ' erho1 ;
  333. SAUTER 3 LIGNES ;
  334. *
  335. *
  336. SI ( (euz1 + erho1 ) >EG 0.05 ) ;
  337. MESS '--------- RESULTATS INCORRECTS -------------' ;
  338. SAUTER 3 LIGNES ;
  339. ERREUR 5;
  340. SINON ;
  341. MESS '--------- RESULTATS CORRECTS -------------' ;
  342. SAUTER 3 LIGNES ;
  343. FINSI ;
  344. *
  345. FIN ;
  346.  
  347.  
  348.  
  349.  
  350.  
  351.  
  352.  
  353.  
  354.  
  355.  

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