Télécharger la plaquette Cast3M


fr en

Télécharger poudre5.dgibi

Retour à la liste

Numérotation des lignes : 

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

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