Télécharger la plaquette Cast3M


fr en

Télécharger endocp1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  endocp1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *           Test Endocp1.dgibi: Jeux de données       *
  5. *           -----------------------------------       *
  6. *                                                     *
  7. *
  8. complet = faux;
  9. *
  10. *************************************************
  11. *  pour calcul complet mettre complet à : vrai;
  12. *************************************************
  13. **
  14. *******************************************************
  15. *                                                     *
  16. * TEST;  MATERIAU ELASTO-PLASTIQUE ENDOMMAGEABLE      *
  17. *                                                     *
  18. * EPROUVETTE EN TRACTION AVEC DEPLACEMENTS IMPOSES.   *
  19. * MAILLAGE AXISYMETRIQUE.                             *
  20. *                                                     *
  21. * LES RESULTATS OBTENUS SONT COMPARES  A LA SOLUTION  *
  22. * ANALYTIQUE OBTENUE PAR ALGORITHME.                  *
  23. * ON TESTE - LA DEFORMATION PLASTIQUE CUMULEE         *
  24. *            (composante EPSE dans VARIABLES_INTERNES)*
  25. *          - D, LA VARIABLE D'ENDOMMAGEMENT           *
  26. *            dans VARIABLES_INTERNES).                *
  27. *          - SIGMAyy ( dans CONTRAINTES ).            *
  28. *                                                     *
  29. * Test effectue par Mr STRUB (cf rapport DMT/91 .558).*
  30. * Test modifie pour evolution 2/06/92.                *
  31. *******************************************************
  32.  
  33. GRAPH = 'N';
  34. OPTI ECHO 0 DIME 2 ELEM QUA8 MODE PLAN CONT ;
  35. *                                                                               
  36. *------------------------------------------------------                   
  37.    P1 = 0 0;          P2 = 5E-3 0;   
  38.    P3 = 5E-3 30E-3;   P4 = 0 30E-3;                  
  39. *                                                                               
  40.    L1 = P1 P2 DROIT 1  ;                                                        
  41.    L2 = P2 P3 DROIT 2 ;                                                         
  42.    L3 = P3 P4 DROIT 1  ;                                                        
  43.    L4 = P4 P1 DROIT 2 ;                                                         
  44. *                                                                               
  45.    EPROU = L1 L2 L3 L4 DALLER PLAN ;                                            
  46.  
  47. SI ( NEG GRAPH 'N');
  48.    trac eprou;
  49. finsi;
  50.    CL1 = BLOQ L1 UY ;                                                           
  51.    CL2 = BLOQ L3 UY ;                                                           
  52.    CL3 = BLOQ L4 UX ;                                                           
  53.    CL = CL1 ET CL2 ET CL3 ;                                                     
  54.    DEP = DEPI CL2 1E-3 ;                                                        
  55. *                                                                               
  56.    MO = MODE EPROU 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE'
  57.         'ENDOMMAGEABLE';
  58. *                                                                               
  59.   PROGEPS = PROG 0.    .47285E-2 .76066E-2 .89226E-2                                       
  60.    1.0533E-2 1.4011E-2 1.9096E-2 2.6081E-2 3.4654E-2                  
  61.    4.5828E-2 6.0813E-2 7.2835E-2 8.3328E-2 9.5526E-2                   
  62.    10.691E-2 11.808E-2 12.784E-2 13.891E-2 14.936E-2                  
  63.    15.700E-2 16.746E-2 17.720E-2 18.946E-2 19.860E-2                   
  64.    20.703E-2 21.607E-2 22.802E-2 23.685E-2 24.919E-2 
  65.    26.083E-2 ;                                              
  66.   PROGSIG = PROG 0.    350.E06   370.59E06 379.18E06                                          
  67.    387.37E06 400.52E06 414.22E06 427.94E06 440.69E06                  
  68.    453.63E06 467.26E06 476.30E06 483.24E06 490.46E06                  
  69.    496.56E06 502.06E06 506.53E06 511.28E06 515.51E06                  
  70.    518.45E06 522.29E06 525.71E06 529.81E06 532.73E06                  
  71.    535.33E06 538.03E06 541.47E06 543.92E06 547.23E06 
  72.    550.24E06 ;                                              
  73.    CTRAC = EVOL MANU EPS PROGEPS SIG PROGSIG ;                                  
  74. *                                                                               
  75. SI ( NEG GRAPH 'N');
  76.     DESS CTRAC;
  77. FINSI;
  78.   MAT = MATE MO YOUN 74020.E06 NU 0.3 RHO 7800. ALPH 0.
  79.         TRAC CTRAC EPSD 0.03 DC 0.23 EPSR 0.25 ;                      
  80. *  
  81.    LI1 = PROG 0. 30. ; LI2 = PROG 0. 30. ;                                      
  82.    EV = EVOL MANU T LI1 DEP(T) LI2 ;                                            
  83.    CHA = CHAR 'DIMP' DEP EV ;                                                                                                                   
  84. *                                                                               
  85. * LISTE DES PAS DE CALCUL :                                                     
  86. *                         0.13 MM INCREMENT ELASTIQUE                           
  87. *                         0.16 MM INCREMENT PLASTIQUE  
  88. si complet;                         
  89.  LIS = PROG 0.13 0.16 0.33 0.5 0.75 1. 1.25 1.5 1.75 2.                      
  90.             2.25 2.5 2.75 3. 3.25 3.5 3.75 4. 4.25 4.5                       
  91.             4.75 5. 5.25 5.5 5.75 6. 6.25 6.5 6.75 7.                        
  92.             7.25 7.5 7.75 8. ;                                               
  93. sinon;
  94.  LIS = prog  0.13 0.16  0.33 0.5 ;
  95. finsi;
  96.     TAB = TABLE ;
  97.     TAB.'BLOCAGES_MECANIQUES' = CL;
  98.     TAB.'MODELE' = MO;
  99.     TAB.'CARACTERISTIQUES' = MAT;
  100.     TAB.'TEMPS_CALCULES' = LIS;
  101.     TAB.'CHARGEMENT' = CHA;
  102.     TAB.'PRECISION' = 1E-08;                                                                                                              
  103. *                                                                               
  104.     PASAPAS TAB;
  105.  
  106. *******************************************************
  107. ************     TEST DE FONCTIONNEMENT  **************
  108. *******************************************************
  109.  
  110. dpc0 = extr TAB.VARIABLES_INTERNES.2 EPSE 1 1 1;
  111. dpc1 = extr TAB.VARIABLES_INTERNES.4 EPSE 1 1 1;
  112. dpc0t = 5.46e-4; dpc1t =1.12694E-02;
  113. si complet;
  114. dpc2 = extr TAB.VARIABLES_INTERNES.10 EPSE 1 1 1;
  115. dpc4 = extr TAB.VARIABLES_INTERNES.18 EPSE 1 1 1;
  116. dpc6 = extr TAB.VARIABLES_INTERNES.26 EPSE 1 1 1;
  117. dpc75 =extr TAB.VARIABLES_INTERNES.32 EPSE 1 1 1;
  118. dpc2t = 6.03e-2; dpc4t = 12.65e-2; dpc6t = 19.29e-2;
  119. dpc75t = 24.27e-2;
  120. finsi;
  121.  
  122.  
  123. D0 =  extr TAB.VARIABLES_INTERNES.2 'VHWD' 1 1 1;
  124. D1 =  extr TAB.VARIABLES_INTERNES.4 'VHWD' 1 1 1;
  125. D0T = 0.; D1T= 0.0; 
  126. si complet;
  127. D2 =  extr TAB.VARIABLES_INTERNES.10 'VHWD' 1 1 1;
  128. D4 =  extr TAB.VARIABLES_INTERNES.18 'VHWD' 1 1 1;
  129. D6 =  extr TAB.VARIABLES_INTERNES.26 'VHWD' 1 1 1;
  130. D75 =  extr TAB.VARIABLES_INTERNES.32 'VHWD' 1 1 1;
  131. D2T= 0.0257;  D4T= 0.0901; D6T= 0.1622; D75T= 0.2199;
  132. finsi;
  133.  
  134. SZZ0 = EXTR TAB.CONTRAINTES.2 SMYY 1 1 1;
  135. SZZ1 = EXTR TAB.CONTRAINTES.4 SMYY 1 1 1;
  136. SZZ0T = 354.3E6; SZZ1T =  4.07984E+08;
  137. si complet;
  138. SZZ2 = EXTR TAB.CONTRAINTES.10 SMYY 1 1 1;
  139. SZZ4 = EXTR TAB.CONTRAINTES.18 SMYY 1 1 1;
  140. SZZ6 = EXTR TAB.CONTRAINTES.26 SMYY 1 1 1;
  141. SZZ75 = EXTR TAB.CONTRAINTES.32 SMYY 1 1 1;
  142.  SZZ2T = 459.3E6; SZZ4T = 461.4E6; SZZ6T = 443.3E6;
  143. SZZ75T = 422.4E6;
  144. finsi;
  145.  
  146. EDPC0 = ABS((DPC0-DPC0T)/DPC0T);
  147. EDPC1 = ABS((DPC1-DPC1T)/DPC1T);
  148. si complet;
  149. EDPC2 = ABS((DPC2-DPC2T)/DPC2T);
  150. EDPC4 = ABS((DPC4-DPC4T)/DPC4T);
  151. EDPC6 = ABS((DPC6-DPC6T)/DPC6T);
  152. EDPC75 = ABS((DPC75-DPC75T)/DPC75T);
  153. finsi;
  154. ED0 = ABS((D0-D0T)/(D0T+0.001));
  155. ED1 = ABS((D1-D1T)/(D1T+0.001));
  156. si complet;
  157. ED2 = ABS((D2-D2T)/D2T);
  158. ED4 = ABS((D4-D4T)/D4T);
  159. ED6 = ABS((D6-D6T)/D6T);
  160. ED75 = ABS((D75-D75T)/D75T);
  161. finsi;
  162. ESZZ0 = ABS((SZZ0-SZZ0T)/SZZ0T);
  163. ESZZ1 = ABS((SZZ1-SZZ1T)/SZZ1T);
  164. si complet;
  165. ESZZ2 = ABS((SZZ2-SZZ2T)/SZZ2T);
  166. ESZZ4 = ABS((SZZ4-SZZ4T)/SZZ4T);
  167. ESZZ6 = ABS((SZZ6-SZZ6T)/SZZ6T);
  168. ESZZ75 = ABS((SZZ75-SZZ75T)/SZZ75T);
  169. finsi;
  170. si complet;
  171.  
  172. LERR = PROG EDPC0 EDPC1 EDPC2 EDPC4 EDPC6 EDPC75 ED0 
  173.        ED1 ED2 ED4 ED6 ED75
  174.        ESZZ0 ESZZ1 ESZZ2 ESZZ4 ESZZ6 ESZZ75;
  175. sinon;
  176. lerr =  PROG EDPC0 EDPC1 ED0 ED1 ESZZ0 ESZZ1;
  177. finsi;
  178. list LERR;
  179. ermax = maxi LERR;
  180. LIST ERMAX;
  181.  
  182. temps;
  183. SI ( ERMAX <EG 0.05 );
  184.    ERRE 0;
  185. SINON;
  186.    ERRE 5;
  187. FINSI;
  188.  
  189. FIN;
  190.  
  191.  
  192.  
  193.  
  194.  
  195.  
  196.  
  197.  
  198.  
  199.  
  200.  
  201.  
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  
  209.  
  210.  

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