Télécharger mistral_axi.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : mistral_axi.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * Test mistral_axi.dgibi: Jeux de donnees *
  6. * --------------------------------------- *
  7. * *
  8. opti echo 0 ;
  9. opti dime 2 mode axis elem qua8 ;
  10. *
  11. *******************************************************
  12. * TEST DE VALIDATION *
  13. * ------------------ *
  14. * *
  15. * MODELE MISTRAL *
  16. * ELASTICITE ET VISCOPLASTICITE *
  17. * *
  18. * MAILLAGE: *
  19. * TUBE CYLINDRIQUE (AXISYMETRIQUE) *
  20. * *
  21. * CHARGEMENT: *
  22. * PRESSION IMPOSEE : RAMPE RAPIDE ET PLATEAU *
  23. *******************************************************
  24. *
  25. * Geometrie
  26. *
  27. * ------------------------------------------
  28. RI = 49.5E-3 ; RE = 50.5E-3 ; H = 0.010 ;
  29. * ------------------------------------------
  30. RM = (RI+RE)/2. ;
  31. RE2MRI2 = (RE**2)-(RI**2) ;
  32. COEFI = (RI**2)/RE2MRI2 ;
  33. E = RE-RI ;
  34. VECR = 1. 0. ;
  35. *
  36. PIB = RI 0 ; PEB = RE 0 ;
  37. PIH = RI H ; PEH = RE H ;
  38. *
  39. LB = droi 2 PIB PEB ;
  40. LE = droi 1 PEB PEH ;
  41. LH = droi 2 PEH PIH ;
  42. LI = droi 1 PIH PIB ;
  43. *
  44. * Maillage
  45. *
  46. EPROU = dall LB LE LH LI ;
  47. * trac EPROU ;
  48. *
  49. POMB = RM 0. ; POMH = RM H ;
  50. PMB = EPROU poin proc POMB ;
  51. PMH = EPROU poin proc POMH ;
  52. *
  53. * Modele et materiau
  54. *
  55. MODG = mode EPROU mecanique elastique orthotrope
  56. viscoplastique mistral ;
  57. *
  58. * ------------------------------------------------------------
  59. SENSIP1 = 1 ; SENSIP2 = -3 ; ICBASE = 0 ;
  60. fichier = '/u/castem/divers/mimataxi_par' ;
  61. PDILT E1 E2 E3 NU12 NU23 NU13 MU12 MU23 MU13
  62. PNBRE PCOHI PECOU PEDIR PRVCE PECRX PDVDI PCROI PINCR
  63. = @mistpar fichier SENSIP1 SENSIP2 ;
  64. * ------------------------------------------------------------
  65. MAZY = mate MODG 'YG1 ' E1 'YG2 ' E2 'YG3 ' E3
  66. 'NU12' NU12 'NU23' NU23 'NU13' NU13
  67. 'G12 ' MU12
  68. 'ALP1' 0. 'ALP2' 0. 'ALP3' 0.
  69. 'DILT' PDILT 'NBRE' PNBRE 'COHI' PCOHI 'ACOU' PECOU 'EDIR' PEDIR
  70. 'RVCE' PRVCE 'ECRX' PECRX 'DVDI' PDVDI 'CROI' PCROI 'INCR' PINCR
  71. 'SIP1' SENSIP1 'SIP2' SENSIP2 'IBAS' ICBASE
  72. 'DIRECTION' VECR ;
  73. *
  74. * Conditions aux limites
  75. *
  76. CLB = bloq uz LB ;
  77. CLH = rela ense uz LH ;
  78. CL = CLB et CLH ;
  79. *
  80. * Chargement
  81. *
  82. * -------------------------------------------------------
  83. T1 = 0.001 ; T2 = 2000.001 ;
  84. DT1 = 25. ;
  85. TEMPSCAL = prog 0. T1 pas DT1 100.001
  86. 150.001 200.001 300.001
  87. 500.001 700.001 1000.001 1400.001 T2 ;
  88. SIGTT1 = 400.E6 ;
  89. TT0 = 623.15 ; TT1 = 623.15 ; TT2 = 623.15 ;
  90. PHIT0 = 2.E25 ;
  91. * -------------------------------------------------------
  92. TEMPS = prog 0. T1 T2 ;
  93. *
  94. PRI = E/RI ;
  95. PRESI = pres mass MODG (PRI) LI ;
  96. SZPRI = PRI*COEFI ;
  97. PRHPI = (-1.)*SZPRI ;
  98. PRESHPI = pres mass MODG (PRHPI) LH ;
  99. EVPI = evol manu TEMPS (prog 0. SIGTT1 SIGTT1) ;
  100. CHAPI = char meca (PRESI et PRESHPI) EVPI ;
  101. *
  102. TT = manu chpo EPROU 1 'T' 1. ;
  103. EVTT = evol manu TEMPS (prog TT0 TT1 TT2) ;
  104. CHTT = char 'T' TT EVTT ;
  105. *
  106. PHI = manu chpo EPROU 1 'FI' 1. ;
  107. EVFI = evol manu TEMPS (prog 0. 0. 0.) ;
  108. CHFI = char 'FI' PHI EVFI ;
  109. *
  110. CHA = CHAPI et CHTT et CHFI ;
  111. *
  112. * Valeurs initiales
  113. *
  114. VINT0 = zero MODG 'VARINTER' ;
  115. FIT0 = manu chml MODG 'FIT ' PHIT0 type 'SCALAIRE' 'STRESSES' ;
  116. VINT0 = VINT0 + FIT0 ;
  117. *
  118. * Calcul
  119. *
  120. TAB = TABLE ;
  121. TAB.'VARIABLES_INTERNES' = TABLE ;
  122. TAB.'BLOCAGES_MECANIQUES' = CL ;
  123. TAB.'CARACTERISTIQUES' = MAZY ;
  124. TAB.'MODELE' = MODG ;
  125. TAB.'CHARGEMENT' = CHA ;
  126. TAB.'VARIABLES_INTERNES' . 0 = VINT0 ;
  127. TAB.'TEMPS_CALCULES' = TEMPSCAL ;
  128. PASAPAS TAB ;
  129. *
  130. * Traitement des resultats
  131. *
  132. TEM = TAB.'TEMPS' ;
  133. SIG = TAB.'CONTRAINTES' ;
  134. DEP = tab.'DEPLACEMENTS' ;
  135. VI = tab.'VARIABLES_INTERNES' ;
  136. NT = dime TEMPSCAL-1 ;
  137. * ---------------
  138. ERMAX = 0.01 ;
  139. * ---------------
  140. T = TEM.NT ;
  141. *
  142. SIGT = SIG.NT ;
  143. SIGT_PO = chang chpo MODG SIGT ;
  144. SIGRR = extr SIGT_PO SMRR PMB ;
  145. SIGTT = extr SIGT_PO SMTT PMB ;
  146. SIGZZ = extr SIGT_PO SMZZ PMB ;
  147. SIGRZ = extr SIGT_PO SMRZ PMB ;
  148. *
  149. DEPT = DEP.NT ;
  150. EPSIT = epsi DEPT MODG ;
  151. EPSIT_PO = chang chpo MODG EPSIT ;
  152. EPSRR = extr EPSIT_PO EPRR PMB ;
  153. EPSTT = extr EPSIT_PO EPTT PMB ;
  154. EPSZZ = extr EPSIT_PO EPZZ PMB ;
  155. EPSRZ = (extr EPSIT_PO GARZ PMB)/2. ;
  156. *
  157. VIT = VI.NT ;
  158. VIT_PO = chang chpo MODG VIT ;
  159. EPSFRR = extr VIT_PO EP11 PMB ;
  160. EPSFTT = extr VIT_PO EP12 PMB ;
  161. EPSFZZ = extr VIT_PO EP13 PMB ;
  162. EPSFRZ = extr VIT_PO EP15 PMB ;
  163. *
  164. mess ;
  165. mess 'T : ' T ;
  166. mess ;
  167. mess 'SIGRR SIGTT SIGZZ : ' SIGRR SIGTT SIGZZ ;
  168. mess 'SIGRZ : ' SIGRZ ;
  169. mess 'EPSRR EPSTT EPSZZ : ' EPSRR EPSTT EPSZZ ;
  170. mess 'EPSRZ : ' EPSRZ ;
  171. mess 'EPSFRR EPSFTT EPSFZZ : ' EPSFRR EPSFTT EPSFZZ ;
  172. mess 'EPSFRZ : ' EPSFRZ ;
  173. mess ;
  174. *
  175. * SIGT_RR = exco SIGT 'SMRR' ;
  176. * trac SIGT_RR MODG ;
  177. * SIGT_TT = exco SIGT 'SMTT' ;
  178. * trac SIGT_TT MODG ;
  179. * SIGT_ZZ = exco SIGT 'SMZZ' ;
  180. * trac SIGT_ZZ MODG ;
  181. * SIGT_RZ = exco SIGT 'SMRZ' ;
  182. * trac SIGT_RZ MODG ;
  183. *
  184. * ------------------------------------------------------------
  185. YOUNGR = 1.1E11 ; YOUNGZ = 0.9E11 ; YOUNGT = 1.E11 ;
  186. POISRZ = 0.3 ; POISZT = 0.3 ; POISRT = 0.3 ;
  187. MURZ = 4.E10 ;
  188. * ------------------------------------------------------------
  189. SIGRR1 = (-1.)*SIGTT1*E/(2.*RM) ;
  190. SIGZZ1 = SIGTT1*SZPRI ;
  191. EPSERR1 = (SIGRR1/YOUNGR)-((POISRT/YOUNGR)*SIGTT1)
  192. -((POISRZ/YOUNGR)*SIGZZ1) ;
  193. EPSETT1 = (SIGTT1/YOUNGT)-((POISZT/YOUNGZ)*SIGZZ1)
  194. -((POISRT/YOUNGR)*SIGRR1) ;
  195. EPSEZZ1 = (SIGZZ1/YOUNGZ)-((POISRZ/YOUNGR)*SIGRR1)
  196. -((POISZT/YOUNGZ)*SIGTT1) ;
  197. * --------------------------
  198. EPSFRR1 = -0.00199068 ;
  199. EPSFTT1 = 0.001931867 ;
  200. EPSFZZ1 = 0.0000588108 ;
  201. * --------------------------
  202. EPSRR1 = EPSERR1 + EPSFRR1 ;
  203. EPSTT1 = EPSETT1 + EPSFTT1 ;
  204. EPSZZ1 = EPSEZZ1 + EPSFZZ1 ;
  205. *
  206. ERSIGRR = (SIGRR/SIGRR1) - 1. ;
  207. ERSIGTT = (SIGTT/SIGTT1) - 1. ;
  208. ERSIGZZ = (SIGZZ/SIGZZ1) - 1. ;
  209. P_ER = prog ERSIGRR ERSIGTT ERSIGZZ ;
  210. ERSIGMA = maxi P_ER 'ABS' ;
  211. *
  212. EREPSRR = (EPSRR/EPSRR1) - 1. ;
  213. EREPSTT = (EPSTT/EPSTT1) - 1. ;
  214. EREPSZZ = (EPSZZ/EPSZZ1) - 1. ;
  215. P_ER = prog EREPSRR EREPSTT EREPSZZ ;
  216. EREPSMA = maxi P_ER 'ABS' ;
  217. *
  218. EREPSFRR = (EPSFRR/EPSFRR1) - 1. ;
  219. EREPSFTT = (EPSFTT/EPSFTT1) - 1. ;
  220. EREPSFZZ = (EPSFZZ/EPSFZZ1) - 1. ;
  221. P_ER = prog EREPSFRR EREPSFTT EREPSFZZ ;
  222. EREPSFMA = maxi P_ER 'ABS' ;
  223. *
  224. si (ERSIGMA > ERMAX) ;
  225. mess ;
  226. mess 'ABS(erreur relative) sur SIGRR, SIGTT ou SIGZZ > ' ERMAX ;
  227. mess ERSIGRR ERSIGTT ERSIGZZ ;
  228. finsi ;
  229. si (EREPSMA > ERMAX) ;
  230. mess ;
  231. mess 'ABS(erreur relative) sur EPSRR, EPSTT ou EPSZZ > ' ERMAX ;
  232. mess EREPSRR EREPSTT EREPSZZ ;
  233. finsi ;
  234. si (EREPSFMA > ERMAX) ;
  235. mess ;
  236. mess 'ABS(erreur relative) sur EPSFRR, EPSFTT ou EPSFZZ > ' ERMAX ;
  237. mess EREPSFRR EREPSFTT EREPSFZZ ;
  238. mess ;
  239. finsi ;
  240. P_ER = prog ERSIGMA EREPSMA EREPSFMA ;
  241. ERMA = maxi P_ER ;
  242. si (ERMA > ERMAX) ;
  243. ERRE 5 ;
  244. finsi ;
  245.  
  246. fin ;
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  

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