Télécharger la plaquette Cast3M


fr en

Télécharger Jpor1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  jpor1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. GRAPH = 'N' ;
  5. SAUT PAGE ;
  6. ************************************************************************
  7. *
  8. *
  9. *
  10. *       TEST SUR UN JOINT POREUX BIDIMENSIONNEL (ELEMENT JOP3)
  11. *
  12. *       UN JOINT POREUX 2D EST BLOQUE EN BAS ET CHARGE EN HAUT
  13. *       AVEC UNE FORCE TANGENTIELLE ET UNE FORCE DE COMPRESSION.
  14. *       UN DEBIT NUL EST IMPOSE PARTOUT. 
  15. *       LES DEPLACEMENTS DE LA FACE 'AB' ET LES PRESSIONS NODALES 
  16. *       SONT COMPARES AUX MEMES VALEURS OBTENUES PAR LA SOLUTION
  17. *       ANALYTIQUE.
  18. *
  19. *
  20. *
  21. *                   A ._______________._______________. B
  22. *
  23. *                   E ._______________________________. F
  24. *
  25. *                   C ._______________._______________. D
  26. *                 ////////////////////////////////////////
  27. *
  28. *
  29. *
  30. ************************************************************************
  31. *
  33. OPTI ECHO 0 ;
  34. *
  35. TITRE 'JOINT POREUX BIDIMENSIONNEL' ;
  36. *
  37. *------------------------------------------------------
  38. *
  39. *        DEFINITION DE LA GEOMETRIE 
  40. *
  41. *------------------------------------------------------
  42. *
  43. OPTI ELEM SEG3 ;
  44. *
  45. A = -1.00 0.00 ;
  46. B =  1.00 0.00 ;
  47. *
  48. AB =  A DROI 1 B ;
  49. *
  50. C = -1.00 0.00 ;
  51. D =  1.00 0.00 ;
  52. *
  53. CD =  C DROI 1 D ;
  54. *
  55. OPTION ELEM SEG2 ;
  56. *
  57. E = -1.00 0.00 ;
  58. F =  1.00 0.00 ;
  59. *
  60. EF =  E DROI 1 F ;
  61. *
  62. *------------------------------------------------------
  63. *
  64. *        DEFINITION DU MAILLAGE 
  65. *
  66. *------------------------------------------------------
  67. *
  68. OPTI ELEM RAP3 ;
  69. *
  70. JOI0 = RACC 0.001 AB CD EF ;
  71. *
  72. *------------------------------------------------------
  73. *
  74. *        DEFINITION DU MODELE ET DU MATERIAU
  75. *
  76. *------------------------------------------------------
  77. *
  78. MOD0 = MODE JOI0 POREUX ;
  79. *
  80. MAT0 = MATE MOD0 KS 1.E10 KN 1.E10 COB 1. MOB 1.E10
  81.                  PERT 1. PERH 1. PERB 1. VISC 1. ;
  82. *
  83. *------------------------------------------------------
  84. *
  85. *        DEFINITION DES CONDITIONS AUX LIMITES
  86. *
  87. *------------------------------------------------------
  88. *
  89. BLO1 = BLOQ UX CD ;
  90. BLO2 = BLOQ UY CD ;
  91. BLO0 = BLO1 ET BLO2 ;
  92. *
  93. FOR1 = FORC FX  1.E10 AB ;
  94. FOR2 = FORC FY -1.E10 AB ;
  95. FOR0 = FOR1 ET FOR2 ;
  96. *
  97. *------------------------------------------------------
  98. *
  99. *        MATRICE DE PERMEABILITE
  100. *
  101. *------------------------------------------------------
  102. *
  103. PER0 = PERM MOD0 MAT0 ;
  104. *
  105. *------------------------------------------------------
  106. *
  107. *        DEFINITION DU CHARGEMENT
  108. *
  109. *------------------------------------------------------
  110. *
  111. L50 = PROG  0.0  1.0 ;
  112. L51 = PROG  0.0  1.0 ;
  113. EV1 = EVOL MANU X L50 Y L51 ;
  114. *
  115. CHA0 = CHAR 'MECA' FOR0 EV1 ;
  116. LIS0 = PROG 0.0 PAS 0.1 1.0 ; 
  117. *
  118. *------------------------------------------------------
  119. *
  120. *        RESOLUTION
  121. *
  122. *------------------------------------------------------
  123. *
  124. TAB1 = TABLE ;
  125. *
  126. TAB1.'MODELE'              = MOD0 ;
  127. TAB1.'CARACTERISTIQUES'    = MAT0 ;
  128. TAB1.'BLOCAGES_MECANIQUES' = BLO0 ;
  129. TAB1.'TEMPS_CALCULES'      = LIS0 ;
  130. TAB1.'CHARGEMENT'          = CHA0 ;
  131. TAB1.'PERMEABILITE'        = PER0 ;
  132. TAB1.'CONSOLIDATION'       = VRAI ;
  133. *
  134. PASAPAS TAB1 ;
  135. *
  136. *------------------------------------------------------
  137. *
  138. *        POST TRAITEMENT
  139. *
  140. *------------------------------------------------------
  141. *
  142. pre_e1 = prog ; pre_s1 = prog ; pre_i1 = prog ;
  143. pre_e2 = prog ; pre_s2 = prog ; pre_i2 = prog ;
  144. def_s1 = prog ; def_s2 = prog ; 
  145. tem0   = prog ;
  146. *
  147. xks = extr mat0 'KS ' 1 1 1 ; xkn = extr mat0 'KN ' 1 1 1 ;
  148. xco = extr mat0 'COB' 1 1 1 ; xmo = extr mat0 'MOB' 1 1 1 ;
  149. *
  150. ndime = dime ( index(tab1.deplacements) ) ;
  151. *
  152. i=0 ;
  153. *
  154. dif1=0 ; dif2=0 ; dif3=0 ; dif4=0 ; 
  155. dif5=0 ; dif6=0 ; dif7=0 ; dif8=0 ;
  156. *
  157. repeter bou1 (ndime-1) ;
  158. *
  159. i=i+1 ;
  160. *
  161. dep0 = tab1.deplacements.i ;
  162. eps0 = epsi mod0 dep0 ;
  163. *
  164. pre1 =  extr dep0 'P   ' a  ;
  165. pre2 =  extr dep0 'P   ' b  ;
  166. pre3 =  extr dep0 'P   ' e  ;
  167. pre4 =  extr dep0 'P   ' f  ;
  168. pre5 =  extr dep0 'P   ' c  ;
  169. pre6 =  extr dep0 'P   ' d  ;
  170. val1 =  extr eps0 'DRSN' 1 1 1 ;
  171. val2 =  extr eps0 'DRSN' 1 1 2 ;
  172. val3 =  extr eps0 'DRSN' 1 1 3 ;
  173. pal1 =  extr eps0 'DRN ' 1 1 1 ;
  174. pal2 =  extr eps0 'DRN ' 1 1 2 ;
  175. pal3 =  extr eps0 'DRN ' 1 1 3 ;
  176. *--------------------------------------------------------------------
  177. *   l'opérateur INTG n'est pas encore disponible pour les éléments
  178. *   joints poreux: intégration à la main
  179. *--------------------------------------------------------------------
  180. P555=0.55555556 ;
  181. P888=0.88888889 ;
  182. *
  183. def1 = ( (P555 * val1) + (P888 * val2) + (P555 * val3) ) / 2. ;
  184. def2 = ( (P555 * pal1) + (P888 * pal2) + (P555 * pal3) ) / 2. ;
  185. *
  186. *
  187. sig1 = (-1.E10 / (ndime-1)) * i / 2. ;
  188. tau1 = ( 1.E10 / (ndime-1)) * i / 2. ;
  189. pref = -1. * sig1 * (xco * xmo) / (((xco**2)*xmo) + xkn) ;
  190. drex = tau1 / xks ;
  191. drey = (sig1 + (xco * pref)) / xkn ;
  192. *
  193. pre_s1 = inse pre_s1 i pre1 ;
  194. pre_s2 = inse pre_s2 i pre2 ;
  195. pre_e1 = inse pre_e1 i pre3 ;
  196. pre_e2 = inse pre_e2 i pre4 ;
  197. pre_i1 = inse pre_i1 i pre5 ;
  198. pre_i2 = inse pre_i2 i pre6 ;
  199. def_s1 = inse def_s1 i def1 ;
  200. def_s2 = inse def_s2 i def2 ;
  201. *
  202. tem0 = inse tem0 i i ;
  203. *
  204. dif1 = dif1 + ( abs ((pre1-pref) / pref) ) ;
  205. dif2 = dif2 + ( abs ((pre2-pref) / pref) ) ;
  206. dif3 = dif3 + ( abs ((pre3-pref) / pref) ) ;
  207. dif4 = dif4 + ( abs ((pre4-pref) / pref) ) ;
  208. dif5 = dif5 + ( abs ((pre5-pref) / pref) ) ;
  209. dif6 = dif6 + ( abs ((pre6-pref) / pref) ) ;
  210. dif7 = dif7 + ( abs ((def1-drex) / drex) ) ;
  211. dif8 = dif8 + ( abs ((def2-drey) / drey) ) ;
  212. *
  213. fin bou1 ;
  214. *
  215. si ( NEG graph 'N' ) ;
  216. *
  217. gra1 = evol manu t0 tem0 pr_su1 pre_s1 ;
  218. gra2 = evol manu t0 tem0 pr_su2 pre_s2 ;
  219. gra3 = evol manu t0 tem0 pr_el1 pre_e1 ;
  220. gra4 = evol manu t0 tem0 pr_el2 pre_e2 ;
  221. gra5 = evol manu t0 tem0 pr_in1 pre_i1 ;
  222. gra6 = evol manu t0 tem0 pr_in2 pre_i2 ;
  223. gra7 = evol manu t0 tem0 de_tan def_s1 ;
  224. gra8 = evol manu t0 tem0 de_nor def_s2 ;
  225. *
  226. dess gra7 ;
  227. dess gra8 ;
  228. dess ( gra1 et gra2 et gra3 et gra4 et gra5 et gra6 ) ;
  229. *
  230. finsi ;
  231. *
  232. *------------------------------------------------------
  233. *
  234. *        CODE DE FONCTIONNEMENT
  235. *
  236. *------------------------------------------------------
  237. *
  238. dif0 = ( dif1 + dif2 + dif3 + dif4 +
  239.          dif5 + dif6 + dif7 + dif8 ) / ( 8. * (ndime-1) ) ;
  240. si ( dif0 < 5e-2 ) ;
  241. erre 0 ;
  242. sinon ;
  243. erre 5 ;
  244. finsi ;
  245. *
  246. fin ;
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  
  255.  
  256.  
  257.  

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