Télécharger la plaquette Cast3M


fr en

Télécharger pore2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  pore2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *           Test Pore2.dgibi: Jeux de données         *
  5. *           ---------------------------------         *
  6. *                                                     *
  7. complet = faux; 
  8. *************************************************
  9. *  pour calcul complet mettre complet à : vrai;
  10. ***************************************************
  11. *******************************************************            
  12. *                                                     *            
  13. *  TEST PORE2                                         *
  14. *                                                     *            
  15. *  CYLINDRE EPAIS EN MILIEU POREUX                    *
  16. *  AVEC EFFETS THERMIQUES                             *
  17. *                                                     *            
  18. *  REFERENCE : Benchmark INTERCLAY 1.4  Basic Model   *
  19. *                                                     *            
  20. *      Le milieu est elastoplastique,                 *
  21. *      modele Von-Mises associe                       *
  22. *                                                     *
  23. *      On decharge la contrainte a l'interieur du     *
  24. *      trou puis on impose un champ de temperature    *
  25. *                                                     *
  26. *      On compare avec les valeurs de reference       *
  27. *      obtenues aussi par CASTEM2000 et en accord     *
  28. *      avec les resultats INTERCLAY                   *
  29. *                                                     *            
  30. *******************************************************            
  31. *                                                                               
  32. OPTI ECHO 1;OPTI DEBU 1;
  33.  
  34. TITRE   ' TEST PORE2  ' ;
  35. GRAPH = 'N';                                                                    
  36. OPTIO DIME 2 MODEL AXIS   ELEM QUA8  ;
  37. *                                                                               
  38. *=====================================================*            
  39. *             geometrie : maillage                    *            
  40. *=====================================================*            
  41. *                                                                               
  42. * Dimension en metres                                                           
  43. *                                                                               
  44. DENS 0.5;                                                                       
  45. P1 = 2.5 0. ; P2 = 25. 0. ; Q1 = 2.5 0.5 ;Q2 = 25. 0.5;                        
  46. PP1 = 5 0. ; QQ1 = 5 0.5 ;                                                      
  47. LI1 = P1 D  P2 'DINI' 0.25 'DFIN' 2.5 ;
  48. LI2 =  P2 D Q2 ;
  49. LI3  = INVE ( Q1 D Q2 'DINI' 0.25 'DFIN' 2.5 );
  50. LI4  = Q1 D P1 ;
  51. LI  = LI1  ET LI2  ET LI3  ET LI4  ;
  52. CADRE = DALLER LI1  LI2  LI3  LI4  PLAN ;
  53. SI (NEG GRAPH 'N');                                                             
  54.    TRAC CADRE ;
  55. FINSI;                                                                          
  56. *                                                                               
  57. *=====================================================*            
  58. *      modele - materiau - caracteristiques en Pa     *
  59. *                 conditions aux limites              *            
  60. *=====================================================*            
  61. *                                                                               
  62.  MO=MODE CADRE POREUX ELASTIQUE PLASTIQUE
  63.                      DRUCKER_PRAGER ;
  64. *
  65.  MAT=MATE MO YOUNG 200.13E6 NU 0.25 MOB  7500E6 COB 1.
  66.  PERM 4.E-12   VISCO 4.
  67.  ALPH (5.E-5 / 3. ) TALP 0. TREF 0.  KF 3000E6   RHOF 0.
  68.  ALPF 40E-5 ALPM (19.E-5 / 3.)  RHO 0.
  69.  ALFA 0.         BETA 0.57735027 K 1.E6
  70.  GAMM 0.         DELT 0.57735027
  71.  ETA  0.         MU   0.57735027 KL 1.E6 H 0. ;
  72. *
  73. *  on bloque les faces superieure et inferieure
  74. *
  75. BL1=BLOQ UZ  LI1 ;                                                              
  76. BL2=BLOQ UZ  LI3 ;                                                              
  77. *
  78. BL = BL1 ET BL2 ;
  79. *                                                                               
  80. *=====================================================*            
  81. *         1-ere etape  : chargement mecanique         *
  82. *   Initialisation de la table pour appel a PASAPAS   *
  83. *=====================================================*            
  84. *                                                                               
  85. *   contraintes initiales
  86. *
  87. ZOZO = MANU CHML MO    SMRZ 0.
  88.    SMRR -5.E6 SMTT -5.E6 SMZZ -5.E6 MSR0 0.
  89.    TYPE CONTRAINTES  STRESSES ;
  90. ZUZU = MANU CHML MO    GARZ 0.
  91.    EPRR 0.    EPTT  0.   EPZZ 0.    PRES 2.5E6
  92.    TYPE DEFORMATIONS  STRESSES ;
  93. FF = BSIGMA MO ZOZO;
  94. EVT1= EVOL MANU 'TEMPS' (PROG 0       14 )
  95.                 'CHARGE' (PROG 1     1 );
  96.  CHA1 = CHAR 'MECA' FF  EVT1 ;
  97. *
  98. *    pression
  99. *
  100. FPI2 = PRESS MASS MO -1.E6   LI4 ;
  101. EVT2= EVOL MANU 'TEMPS' (PROG 0 PAS 1 14 )
  102.                 'CHARGE' (PROG 0 PAS 0.25 3.5) ;
  103. CHA2 = CHAR 'MECA' FPI2 EVT2 ;
  104. CHATT = CHA1 ET CHA2 ;
  105. *
  106. TAB1 = TABLE;
  107. TAB1 . 'CONTRAINTES' = TABLE;
  108. TAB1 . 'DEFORMATIONS' = TABLE;
  109. TAB1 . 'DEPLACEMENTS' = TABLE;
  110. TAB1.'CONTRAINTES' . 0 = ZOZO;
  111. TAB1.'DEFORMATIONS'. 0 = ZUZU;
  112. TAB1.'DEPLACEMENTS' . 0 = MANU CHPO CADRE  1 P 2.5E6 ;
  113. TAB1.'TEMPS_CALCULES' = PROG 0 6  10 14;
  114. TAB1.'MODELE' = MO;
  115. TAB1.'CARACTERISTIQUES' = MAT;
  116. TAB1.'BLOCAGES_MECANIQUES' = BL;
  117. TAB1.'CHARGEMENT' = CHATT;
  118. TAB1.'PRECISION'=5.E-7;
  119. *
  120. TMASAU=table;
  121. tab1 . 'MES_SAUVEGARDES'=TMASAU;
  122. TMASAU .'DEFTO'=VRAI;
  123. TMASAU .'DEFIN'=VRAI;
  124. TAB1.'MONO'=monop;tab1.'NB_EL_PARA'=0.9;
  125. TAB1.'HYPOTHESE_DEFORMATIONS'='LINEAIRE' ;
  126. PASAPAS TAB1 ;
  127. *
  128. *=====================================================*            
  129. *            2-eme etape  : temperature               *
  130. *   Initialisation de la table pour appel a PASAPAS   *
  131. *=====================================================*            
  132. *
  133. TAB2 = TABLE;
  134. TAB2 . 'VARIABLES_INTERNES' = TABLE;
  135. TAB2 . 'DEPLACEMENTS' = TABLE;
  136. TAB2 . 'CONTRAINTES' = TABLE;
  137. TAB2 . 'DEFORMATIONS' = TABLE;
  138. SS0 = PECHE TAB1 CONTRAINTES ;
  139. EE0 = PECHE TAB1 DEFORMATIONS ;
  140. TAB2.'CONTRAINTES' . 0 = SS0;
  141. TAB2.'DEFORMATIONS' . 0 = EE0;
  142. TAB2.'PRECISION'=5.E-7;
  143. TAB2.'VARIABLES_INTERNES' . 0 = 
  144.                      PECHE TAB1 VARIABLES_INTERNES ;
  145. TAB2.'DEPLACEMENTS' . 0 = PECHE TAB1 DEPLACEMENTS  ;
  146. *
  147. FF = BSIGMA MO SS0 ;
  148. EVT1= EVOL MANU 'TEMPS'  (PROG 0     100)
  149.                 'CHARGE' (PROG 1     1 );
  150. CHA1 = CHAR 'MECA' FF EVT1 ;
  151. *
  152. LEFAC = 1. / ( LOG ( 2.5/25.) ) ;
  153. RAD   = COOR 1 CADRE ;
  154. DT  = NOMC (  LEFAC * ( LOG (RAD / 25. ) ) ) T ;
  155. *
  156.  TABTER = TABLE ;TEMP = TABLE;
  157.  TEMP . 0 = 0.;
  158.  TEMP . 1 = 100.;
  159.  TABTER . 0 = MANU CHPO CADRE 1 T 0. ;
  160.  TABTER . 1 = (DT * 100.)  ;
  161.  CHA2 = CHAR 'T' TEMP TABTER;
  162. *
  163. CHATT = CHA1 ET CHA2;
  164. *
  165. si complet;
  166. tfin = 70.;
  167. sinon;
  168. tfin = 15.;
  169. finsi;
  170. TAB2.'TEMPS_CALCULES' = PROG 0  PAS 5.  tfin; ;
  171. TAB2.'MODELE' = MO;
  172. TAB2.'CARACTERISTIQUES' = MAT;
  173. TAB2.'BLOCAGES_MECANIQUES' = BL;
  174. TAB2.'CHARGEMENT' = CHATT;
  175. *
  176. TMASAU=table;
  177. tab2 . 'MES_SAUVEGARDES'=TMASAU;
  178. TMASAU .'DEFTO'=VRAI;
  179. TMASAU .'DEFIN'=VRAI;
  180. TAB2.'PROCESSEURS'= MOT 'AUTOMATIQUE';
  181. TAB2.'HYPOTHESE_DEFORMATIONS'='LINEAIRE' ;
  182. PASAPAS TAB2 ;
  183. *                                                                               
  184. *=====================================================*            
  185. *              affichage des resultats                *            
  186. *=====================================================*            
  187. *                                                                               
  188. MESS ' RESULTATS TEST PORE2 ';
  189. MESS ' -------------------- ';                                                  
  192. '   On cherche les contraintes et la pression ';
  194. '   Elles sont comparees a des valeurs de reference ';
  197. *                                                                               
  198. *----------- r{sultats --------------------------------      
  199. *                                                                               
  200. *    DISTRIBUTION DES CONTRAINTES ET PRESSION
  201. *    POUR T = 70.
  202. *
  203. si complet;
  204. TT = 70. ;
  205. sinon;
  206. TT = 15.;
  207. finsi;
  208. SIGAA = PECHE TAB2 CONTRAINTES  TT;
  209. EPSAA = PECHE TAB2 DEFORMATIONS  TT;
  210. SIGAB = CHAN 'CHPO' MO SIGAA ;
  211.  EPSAB = CHAN 'CHPO' MO EPSAA ;
  212. SIGARR= EVOL CHPO  SIGAB SMRR LI1 ;
  213. SIGATT= EVOL CHPO  SIGAB SMTT LI1 ;
  214. LAPP =  EVOL CHPO  EPSAB PRES LI1 ;
  215. SI (NEG GRAPH 'N');                                                             
  216.    DESS SIGARR ; DESS SIGATT ; DESS LAPP ;
  217. FINSI;                                                                          
  218. EVSIGAR = EXTR SIGARR ORDO 1;
  219. EVSIGAT = EXTR SIGATT ORDO 1;
  220. EVPP    = EXTR LAPP   ORDO 1;
  221. *
  222. *   VALEURS DE REFERENCE  calculees a la date du 1/2/93
  223. *
  224. si complet;
  225. *  CONTRAINTE RADIALE
  226. SOLRR = PROG
  227.  -1.59459E+06 -1.56051E+06 -1.78383E+06 -1.76850E+06 
  228.  -1.96931E+06 -1.97473E+06 -2.15711E+06 -2.17942E+06 
  229.  -2.34668E+06 -2.38270E+06 -2.53753E+06 -2.58467E+06 
  230.  -2.72927E+06 -2.78543E+06 -2.92159E+06 -2.98507E+06 
  231.  -3.11427E+06 -3.18371E+06 -3.30718E+06 -3.38150E+06
  232.  -3.50030E+06 -3.57863E+06 -3.69361E+06 -3.77540E+06 
  233.  -3.88744E+06;
  234. SOLRR = SOLRR ET (PROG
  235.  -3.97178E+06 -4.08114E+06 -4.16891E+06 -4.29345E+06 
  236.  -4.34768E+06 -4.48547E+06 -4.54414E+06 -4.62704E+06 
  237.  -4.68421E+06 -4.74775E+06 -4.79065E+06 -4.83863E+06 
  238.  -4.86963E+06 -4.90481E+06 -4.92612E+06 -4.95084E+06 
  239.  -4.96421E+06 -4.98036E+06 -4.98723E+06 -4.99637E+06
  240.  -4.99792E+06 -5.00142E+06 );
  241. *
  242. *  CONTRAINTE TANGENTIELLE
  243. SOLTT = PROG
  244.  -3.59051E+06 -3.55538E+06 -3.77813E+06 -3.76143E+06
  245.  -3.96142E+06 -3.96508E+06 -4.14632E+06 -4.16639E+06 
  246.  -4.33216E+06 -4.36546E+06 -4.51847E+06 -4.56247E+06 
  247.  -4.70501E+06 -4.75775E+06 -4.89181E+06 -4.95180E+06 
  248.  -5.07914E+06 -5.14532E+06 -5.26743E+06 -5.33906E+06
  249.  -5.45722E+06 -5.53370E+06 -5.64905E+06 -5.72971E+06 
  250.  -5.84300E+06;
  251. SOLTT = SOLTT ET ( PROG
  252.  -5.92761E+06 -6.03950E+06 -6.12667E+06 -6.24105E+06 
  253.  -6.30245E+06 -6.28413E+06 -6.13481E+06 -6.01293E+06 
  254.  -5.89978E+06 -5.80418E+06 -5.70545E+06 -5.62098E+06 
  255.  -5.53670E+06 -5.46384E+06 -5.39075E+06 -5.32680E+06 
  256.  -5.26302E+06 -5.20651E+06 -5.15031E+06 -5.09985E+06
  257.  -5.04987E+06 -5.00444E+06 );
  258. *
  259. *  PRESSION
  260. SOLPP = PROG
  261.   1.36747E+06  1.44034E+06  1.51321E+06  1.58674E+06  
  262.   1.66027E+06  1.73387E+06  1.80747E+06  1.88113E+06  
  263.   1.95479E+06  2.02829E+06  2.10179E+06  2.17490E+06  
  264.   2.24800E+06  2.32049E+06  2.39297E+06  2.46467E+06  
  265.   2.53636E+06  2.60709E+06  2.67783E+06  2.74773E+06
  266.   2.81762E+06  2.88621E+06  2.95480E+06  3.02386E+06  
  267.   3.09291E+06;
  268. SOLPP = SOLPP ET (PROG
  269.   3.15682E+06  3.22073E+06  3.29927E+06  3.37782E+06  
  270.   3.40824E+06  3.43867E+06  3.36939E+06  3.30011E+06  
  271.   3.24558E+06  3.19105E+06  3.13296E+06  3.07487E+06  
  272.   3.01765E+06  2.96042E+06  2.90299E+06  2.84556E+06  
  273.   2.78819E+06  2.73081E+06  2.67343E+06  2.61604E+06
  274.   2.55866E+06  2.50128E+06 );
  275. sinon;
  276. SOLRR = prog 
  277.  -1.56503E+06 -1.57213E+06 -1.75778E+06 -1.77784E+06 
  278.  -1.94803E+06 -1.98244E+06 -2.14000E+06 -2.18611E+06 
  279.  -2.33334E+06 -2.38903E+06 -2.52782E+06 -2.59133E+06 
  280.  -2.72321E+06 -2.79313E+06 -2.91938E+06 -2.99454E+06 
  281.  -3.11618E+06 -3.19563E+06 -3.31352E+06 -3.39647E+06
  282.  -3.51129E+06 -3.59714E+06 -3.70949E+06 -3.79760E+06 
  283.  -3.90787E+06;
  284. solrr = solrr et (prog 
  285.  -3.99816E+06 -4.12859E+06 -4.17682E+06 -4.32275E+06 
  286.  -4.37916E+06 -4.46316E+06 -4.52542E+06 -4.59310E+06 
  287.  -4.64335E+06 -4.69793E+06 -4.73815E+06 -4.78194E+06 
  288.  -4.81401E+06 -4.84898E+06 -4.87436E+06 -4.90211E+06 
  289.  -4.92201E+06 -4.94386E+06 -4.95928E+06 -4.97629E+06
  290.  -4.98804E+06 -5.00116E+06);
  291. soltt = prog 
  292.  -3.56189E+06 -3.56846E+06 -3.75410E+06 -3.77364E+06 
  293.  -3.94386E+06 -3.97775E+06 -4.13540E+06 -4.18103E+06 
  294.  -4.32842E+06 -4.38367E+06 -4.52270E+06 -4.58583E+06 
  295.  -4.71802E+06 -4.78765E+06 -4.91426E+06 -4.98918E+06 
  296.  -5.11126E+06 -5.19054E+06 -5.30890E+06 -5.39174E+06
  297.  -5.50708E+06 -5.59285E+06 -5.70571E+06 -5.79390E+06 
  298.  -5.90481E+06;
  299. soltt = soltt et ( prog
  300.  -5.99480E+06 -6.11303E+06 -6.16665E+06 -6.10728E+06 
  301.  -6.00175E+06 -5.91312E+06 -5.83171E+06 -5.76581E+06 
  302.  -5.69639E+06 -5.63972E+06 -5.58251E+06 -5.53566E+06 
  303.  -5.48786E+06 -5.44850E+06 -5.40854E+06 -5.37545E+06 
  304.  -5.34188E+06 -5.31389E+06 -5.28556E+06 -5.26174E+06
  305.  -5.23769E+06 -5.21735E+06);
  306.  
  307. solpp = prog 
  308.   3.52566E+05  4.45884E+05  5.39202E+05  6.32610E+05  
  309.   7.26018E+05  8.19101E+05  9.12184E+05  1.00502E+06  
  310.   1.09786E+06  1.19042E+06  1.28298E+06  1.37529E+06  
  311.   1.46759E+06  1.55963E+06  1.65167E+06  1.74345E+06  
  312.   1.83523E+06  1.92692E+06  2.01860E+06  2.10969E+06
  313.   2.20078E+06  2.29305E+06  2.38533E+06  2.47274E+06  
  314.   2.56016E+06;
  315. solpp = solpp et ( prog  
  316.   2.66334E+06  2.76651E+06  2.80560E+06  2.84469E+06    
  317.   2.82168E+06  2.79866E+06  2.78895E+06  2.77924E+06  
  318.   2.76631E+06  2.75339E+06  2.74124E+06  2.72909E+06  
  319.   2.71676E+06  2.70442E+06  2.69213E+06  2.67984E+06  
  320.   2.66755E+06  2.65525E+06  2.64295E+06  2.63066E+06
  321.   2.61836E+06  2.60606E+06);
  322.  
  323. finsi;
  324.  
  325. *
  326. *   ecarts
  327. *
  328.  
  329. ERSIGAR =(MAXI ABS ( EVSIGAR - SOLRR ))/
  330.          (MAXI ABS SOLRR);
  331. ERSIGAT =(MAXI ABS ( EVSIGAT - SOLTT ))/
  332.          (MAXI ABS SOLTT);
  333. ERPP    =( MAXI ABS (EVPP - SOLPP))/
  334.          (MAXI ABS SOLPP);
  335. *
  338. 'Ecart relatif en Contraintes radiales           :' ERSIGAR ;
  341. 'Ecart relatif en Contraintes circonferentielles :' ERSIGAT ;
  344. 'Ecart relatif en Pressions                      :' ERPP    ;
  345. *                                                                               
  346. *                                                                               
  347. *=====================================================*            
  348. *                  code fonctionnement                *            
  349. *=====================================================*            
  350. *                                                                               
  351. * L'ecart maximum entre valeurs de reference et 
  352. * calculee doit etre inferieur a 4.d-4
  353. *                                                                               
  354. SAUT 2 LIGNE;                                                                   
  355. SI (ERSIGAR <EG 0.001) ;
  356.   SI (ERSIGAT <EG 0.001) ;
  357.     SI (ERPP <EG 0.001) ;
  358.         ERRE 0;
  359.     SINON;
  360.         ERRE 5;
  361.     FINSI;
  362.   SINON;
  363.     ERRE 5;
  364.   FINSI;
  365. SINON;
  366.   ERRE 5;
  367. FINSI;
  368. *                                                                               
  370. *                                                                               
  371. FIN;                                                                            
  372.  
  373.  
  374.  
  375.  
  376.  
  377.  
  378.  
  379.  
  380.  
  381.  
  382.  
  383.  
  384.  
  385.  
  386.  
  387.  
  388.  
  389.  
  390.  
  391.  
  392.  
  393.  
  394.  
  395.  
  396.  
  397.  
  398.  
  399.  
  400.  
  401.  
  402.  

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