Télécharger pore3.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : pore3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Pore3.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. complet =faux;
  8. *
  9. *************************************************
  10. * pour calcul complet mettre complet à : vrai;
  11. *************************************************
  12. *******************************************************
  13. * *
  14. * TEST PORE3 *
  15. * *
  16. * CONSOLIDATION UNIDIMENSIONNELLE *
  17. * *
  18. * REFERENCE : Probleme de Terzaghi *
  19. * *
  20. * Le milieu est poroélastique *
  21. * On charge à la surface par une pression *
  22. * On compare avec la solution analytique *
  23. * *
  24. *******************************************************
  25. *
  26. OPTI ECHO 1 ;
  27. * ajout de option epsilon lineaire pour la precision des test!
  28. OPTION epsilon lineaire;
  29.  
  30. TITRE ' TEST PORE3 ' ;
  31. GRAPH = 'N';
  32. OPTIO DIME 2 ELEM QUA8 MODE PLAN DEFO ;
  33. MOTP = MOTS 'P' ;
  34. LREP = PROG 1.E6 ;
  35. OPTI NORM MOTP LREP ;
  36. *
  37. *=====================================================*
  38. * geometrie : maillage *
  39. *=====================================================*
  40. *
  41. DENS 1. ;
  42. *
  43. H = 10. ;H2 = H * H ;
  44. P1 = 0. 0. ; P2 = 1. 0. ; P3 = 1. 10. ; P4 = 0. 10. ;
  45. LI1= P1 D 1 P2 ;
  46. LI2= P2 D P3 DINI 1. DFIN 0.1 ;
  47. LI3= P3 D 1 P4 ;
  48. LI4= P4 D P1 DINI 0.1 DFIN 1. ;
  49. LI = LI1 ET LI2 ET LI3 ET LI4 ;
  50. CADRE = DALLER LI1 LI2 LI3 LI4 PLAN ;
  51. SI (NEG GRAPH 'N');
  52. TRAC CADRE ;
  53. FINSI;
  54. *
  55. *=====================================================*
  56. * modele - materiau *
  57. * rigidite - conditions aux limites *
  58. *=====================================================*
  59. *
  60. MO=MODE CADRE POREUX;
  61. *
  62. LAM = 1.D10 ;
  63. MU = 2.5D9 ;
  64. PK = 1.D-9 ;
  65. *
  66. MAT=MATE MO YOUNG 7000D6 NU 0.4 MOB 1.D18
  67. COB 1. PERM PK VISCO 1.
  68. ALPH 0. KF 1. RHOF 0. ALPF 0. ALPM 0. RHO 0. ;
  69. BL1=BLOQ UY LI1 ;
  70. BL2=BLOQ UX (LI2 ET LI4) ;
  71. BLP = BLOQ P LI3 ;
  72. *
  73. * on impose une pression sur le dessus
  74. *
  75. P0 = 4.D3 ;
  76. FPP = PRESSION MASS MO P0 LI3;
  77. *
  78. *=====================================================*
  79. * consolidation *
  80. *=====================================================*
  81. *
  82. PER = PERM MO MAT;
  83. ev = evol manu t ( prog 0 1.D-5 1000) F ( PROG 0 1 1);
  84. CHA = CHAR 'MECA' FPP EV ;
  85. si complet;
  86. LETEMP= PROG 0 PAS 1.D-2 1.D-1 PAS 1.D-1 1.D0
  87. PAS 5.D-1 1.D1;
  88. sinon;
  89. letemp = prog 0 PAS 1.D-2 1.D-1 PAS 1.D-1 5.d-1;
  90. finsi;
  91.  
  92. NBT = DIME LETEMP ;
  93. TFINAL = EXTR LETEMP NBT;
  94.  
  95. *
  96. *======================================================
  97. * remplissage de la table
  98. *======================================================
  99. *
  100. TAB1 = TABLE ;
  101. TAB1.'MODELE' = MO;
  102. TAB1.'CARACTERISTIQUES' = MAT;
  103. TAB1.'BLOCAGES_MECANIQUES' = BL1 ET BL2 ET BLP;
  104. TAB1.'TEMPS_CALCULES' = LETEMP;
  105. TAB1.'CHARGEMENT' = CHA;
  106. TAB1.'CONSOLIDATION' = VRAI;
  107. TAB1.'TETA' = 0.55 ;
  108. TAB1.'PERMEABILITE' = PER ;
  109. *
  110. * appel a PASAPAS
  111. *
  112. PASAPAS TAB1;
  113. *
  114. *=====================================================*
  115. * affichage des resultats *
  116. *=====================================================*
  117. *
  118. MESS ' RESULTATS TEST PORE3 ';
  119. MESS ' -------------------- ';
  120. ' On cherche le deplacement vertical et la pression ';
  121. ' on compare aux solutions analytiques ';
  122. *
  123. TTA = INDEX ( TAB1.DEPLACEMENTS);
  124. NT = DIME TTA ;
  125. LIST NT;
  126. IT = 0 ;
  127. REPETER BOPI NT;
  128. IT = IT + 1 ;
  129. LET = TTA . IT;
  130. DD= TAB1.DEPLACEMENTS.LET;
  131. SI ( EGA IT 1 ) ;
  132. PRT = PROG LET ;
  133. PRD = PROG ( EXTR DD UY P3);
  134. EV = EVOL CHPO DD P LI2 ;
  135. SINON;
  136. PRT = PRT ET (PROG LET) ;
  137. PRD = PRD ET (PROG ( EXTR DD UY P3)) ;
  138. EV = EV ET (EVOL CHPO DD P LI2) ;
  139. FINSI;
  140. FIN BOPI ;
  141. FV = EVOL ROSE MANU 'T' PRT 'UY' PRD;
  142. *
  143. DD= PECHE TAB1 DEPLACEMENTS ;
  144. EV1= EVOL ROSE CHPO DD P LI2 ;
  145. EVD1= EVOL ROSE CHPO DD UY LI2 ;
  146. *
  147. *=====================================================*
  148. * determination de la solution analytique
  149. *=====================================================*
  150. *
  151. X = COOR 2 LI2;
  152. XX = ( MANU CHPO LI2 1 SCAL 1 ) - (X / H) ;
  153. XX = NOMC XX P ;
  154. FACT = (LAM + ( 2. * MU )) * ( PK / H2 ) ;
  155. FACU = -1. * ( (LAM + ( 2. * MU )) / ( P0 * H ) ) ;
  156. PI = 3.141592654D0;
  157. PI2 = PI / 2.D0 ;
  158. DEGPI = 180. / PI ;
  159. *
  160. ITT = 0 ;
  161. REPE BOUT NBT;
  162. ITT = ITT + 1 ;
  163. T= EXTR LETEMP ITT ;
  164. PP = MANU CHPO LI2 1 P 0. ;
  165. UU = X / H ;
  166. UU = NOMC UU P ;
  167. TT = FACT * T;
  168. N = 0;
  169. si complet ;
  170. nfoi = 50;
  171. sinon;
  172. nfoi = 10;
  173. finsi;
  174. REPE BOU nfoi;
  175. M = PI2 * ( (2. * N ) + 1. ) ;
  176. M2 = M*M ;
  177. EX = EXP ( ( -1. * M2) * TT ) ;
  178. SSI = SIN ( (M * DEGPI) * XX ) ;
  179. CCO = COS ( (M * DEGPI) * XX ) ;
  180. PP = PP + ( ( ( 2. / M ) * EX) * SSI) ;
  181. UU = UU - ( ( ( 2. / M2) * EX) * CCO) ;
  182. N=N + 1 ;
  183. FIN BOU ;
  184. PP = PP * P0 ;
  185. UU = UU / FACU ;
  186. SI ( EGA ITT 1 ) ;
  187. PGT = PROG T ;
  188. PGU = PROG ( EXTR UU 'P' P3) ;
  189. SINON ;
  190. PGT = PGT ET (PROG T) ;
  191. PGU = PGU ET (PROG ( EXTR UU 'P' P3)) ;
  192. FINSI;
  193. *
  194. SI ( EGA ITT NBT ) ;
  195. EV2= EVOL TURQ CHPO PP 'P' LI2;
  196. SI (NEG GRAPH 'N');
  197. TITRE ' COMPARAISON DES PROFILS DE PRESSION ' ;
  198. DESS (EV1 ET EV2 );
  199. FINSI;
  200. EVD2= EVOL TURQ CHPO UU 'P' LI2;
  201. SI (NEG GRAPH 'N');
  202. TITRE
  203. ' COMPARAISON DES PROFILS DE DEPLACEMENT VERTICAL';
  204. DESS (EVD1 ET EVD2 );
  205. FINSI;
  206. FINSI;
  207.  
  208. FIN BOUT;
  209. FVT= EVOL TURQ MANU 'T' PGT 'UY' PGU;
  210. SI (NEG GRAPH 'N');
  211. DESS ( FV ET FVT) ;
  212. FINSI;
  213. *
  214. * ecarts
  215. *
  216. EVY1 = EXTR EV1 ORDO ;
  217. EVY2 = EXTR EV2 ORDO ;
  218. EVDY1 = EXTR EVD1 ORDO ;
  219. EVDY2 = EXTR EVD2 ORDO ;
  220. ERPP =(MAXI ABS ( EVY1 - EVY2 ))/(MAXI ABS EVY2 );
  221. ERDUY =( MAXI ABS (EVDY1 - EVDY2 ))/(MAXI ABS EVDY2);
  222. *
  223. *
  224. 'Ecart relatif en pression :' ERPP ;
  225. 'Ecart relatif en profil de deplacement vertical :' ERDUY ;
  226. *
  227. *
  228. *=====================================================*
  229. * code fonctionnement *
  230. *=====================================================*
  231. *
  232. * L'ecart maximum entre valeurs de reference et
  233. * calculee doit etre inferieur a 6.d-3
  234. *
  235. SAUT 2 LIGNE;
  236. SI (ERDUY <EG 0.1) ;
  237. SI (ERPP <EG 0.1) ;
  238. ERRE 0;
  239. SINON;
  240. ERRE 5;
  241. FINSI;
  242. SINON;
  243. ERRE 5;
  244. FINSI;
  245. *
  246. *
  247. FIN ;
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  
  255.  
  256.  
  257.  
  258.  
  259.  
  260.  
  261.  
  262.  
  263.  
  264.  
  265.  

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