Télécharger la plaquette Cast3M


fr en

Télécharger cube.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  cube.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. SAUT PAGE;                                                                      
  5. ************************************************************************        
  6. *     CUBE EN TRACTION UNIAXIALE                                                
  7. *                                                                               
  8. *     TEST ELEMENTAIRE DU GROUPE DE TRAVAIL 'STATIQUE NON LINEAIRE'             
  9. *     COMMISSION VPCS                                                           
  10. *                                                                               
  11. *     LE CALCUL MARCHE ET DONNE EXACTEMENT LES RESULTATS THEORIQUES             
  12. *     LE 05/05/90                                                               
  13. *                                                                               
  14. *     Un cube constitu{ de 2 {l{ments CUB8 est soumis @ une traction            
  15. *     uniaxiale (blocages des d{placements suivant Y et Z et d{place-           
  16. *     ment impos{ suivant X aux deux extr{mit{s (surfaces SU1 et SU2))          
  17. *                                                                               
  18. *     Les contraintes obtenues sont compar{es aux r{sultats th{oriques          
  19. *     et donnent les r{sultats pr{cis{s plus bas.                               
  20. *
  21. *     Ref: Hinton Ezatt .Fundamentals tests for two and three dimensio-
  22. *      -nal small strain, elastoplastic finite element analysis. April
  23. *      1987 NAFEMS
  24. **********************************************************************
  25. opti echo 0;
  26.  
  27. SAUT PAGE;                                                                      
  28. TITRE 'CUBE EN TRACTION UNIAXIALE';                                             
  29. OPTION DIME 3 ECHO 0;
  30.  
  31. graph = 'N';
  32.  
  33. TEMPS ;                                                                         
  34. DENS 1;                                                                         
  35. *-----------------------------------------------------------------------        
  36. ***** DEFINITION DE LA GEOMETRIE                                                
  37. *-----------------------------------------------------------------------        
  38. A = 1.;                                                                         
  39. A1 = -1.;                                                                       
  40. A2 = .5;                                                                        
  41. B = 1.;                                                                         
  42. C = 1.;                                                                         
  43. P1 = A 0 0; Q1 = A2 0 0; R1 = 0 0 0;                                            
  44. P2 = A B 0; Q2 = A2 B 0; R2 = 0 B 0;                                            
  45. P3 = A B C; Q3 = A2 B C; R3 = 0 B C;                                            
  46. P4 = A 0 C; Q4 = A2 0 C; R4 = 0 0 C;                                            
  47. TR = A1 0 0;                                                                    
  48. OPTION ELEM SEG2;                                                               
  49. LI1 = P1 DROIT 1 P2;                                                            
  50. LI2 = P2 DROIT 1 P3;                                                            
  51. LI3 = P3 DROIT 1 P4;                                                            
  52. LI4 = P4 DROIT 1 P1;                                                            
  53. CO1 = LI1 ET LI2 ET LI3 ET LI4;                                                 
  54. LJ1 = Q1 DROIT 1 Q2;                                                            
  55. LJ2 = Q2 DROIT 1 Q3;                                                            
  56. LJ3 = Q3 DROIT 1 Q4;                                                            
  57. LJ4 = Q4 DROIT 1 Q1;                                                            
  58. CO2 = LJ1 ET LJ2 ET LJ3 ET LJ4;                                                 
  59. LK1 = R1 DROIT 1 R2;                                                            
  60. LK2 = R2 DROIT 1 R3;                                                            
  61. LK3 = R3 DROIT 1 R4;                                                            
  62. LK4 = R4 DROIT 1 R1;                                                            
  63. CO3 = LK1 ET LK2 ET LK3 ET LK4;                                                 
  64. OPTION ELEM QUA4;                                                               
  65. SU1 = SURFACE CO1 PLANE;                                                        
  66. SU2 = SURFACE CO2 PLANE;                                                        
  67. SU3 = SURFACE CO3 PLANE;                                                        
  68. OPTION ELEM CUB8;                                                               
  69. VOL1 = SU1 VOLU 1 SU2;                                                          
  70. VOL2 = SU2 VOLU 1 SU3;                                                          
  71. VOL = VOL1 ET VOL2;                                                             
  72. ELIM .001 VOL;                                                                  
  73. *-----------------------------------------------------------------------        
  74. ***** TRACE                                                                     
  75. *-----------------------------------------------------------------------        
  76.  
  77. SI ( NEG graph 'N');
  78. OEIL = (30 20 20);                                                              
  79. TRAC OEIL VOL;                                                                  
  80. finsi;
  81.  
  82. *-----------------------------------------------------------------------        
  83. ***** CONDITIONS AUX LIMITES                                                    
  84. *-----------------------------------------------------------------------        
  85. CL1 = BLOQUE SU3 UX;                                                            
  86. CL2 = BLOQUE VOL UY UZ;                                                         
  87. CL3 = BLOQUE SU2 UX;                                                            
  88. CL4 = BLOQ SU1 UX;                                                              
  89. CL = CL1 ET CL2 ET CL3 ET CL4;                                                  
  90. *-----------------------------------------------------------------------        
  91. * VALEURS IMPOSEES DES DEPLACEMENTS                                             
  92. *-----------------------------------------------------------------------        
  93. EPS0 = .25E-4;                                                                  
  94. COEF1 = A * EPS0 * .5;                                                          
  95. COEF2 = A * EPS0;                                                               
  96. D1 = DEPI CL3 COEF1;                                                            
  97. D2 = DEPI CL4 COEF2;                                                            
  98. DEP = D1 ET D2;                                                
  99. *-----------------------------------------------------------------------        
  100. ***** CALCUL MODELE
  101. *-----------------------------------------------------------------------        
  102. MO = MODE VOL 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'PARFAIT' CUB8;
  103.  
  104. MAT1 = MATE MO YOUN 2.5E11 SIGY 5.0E6 NU .25;
  105.  
  106. *-----------------------------------------------------------------------        
  107. ***** CHAMPS RELATIFS A L'OPERATEUR PASAPAS                                     
  108. *-----------------------------------------------------------------------        
  109. LP1 = PROG 0. 1. 2. 3. 4. ;
  110. LP2 = PROG 0. 1. 2. 1. 0.;
  111. EV = EVOL MANU T LP1 F(T) LP2;  
  112. CHA1 = CHAR 'DIMP' DEP EV;  
  113. LIS1 = PROG 1. 2. 3. 4. ;
  114. *-----------------------------------------------------------------------        
  115. *****  CALCUL NON LINEAIRE                                                      
  116. *-----------------------------------------------------------------------        
  117. TAB1 = TABLE;
  118. TAB1.'BLOCAGES_MECANIQUES' = CL;
  119. TAB1.'CARACTERISTIQUES' = MAT1;
  120. TAB1.'TEMPS_CALCULES' = LIS1;
  121. TAB1.'MODELE' = MO;
  122. TAB1.'CHARGEMENT' = CHA1;                                                                                
  123. PASAPAS TAB1;
  124. *-----------------------------------------------------------------------        
  125. ***** EXTRACTION DES RESULTATS                                                  
  126. *-----------------------------------------------------------------------        
  127. SI1 = TAB1.CONTRAINTES.1;                                                         
  128. SI2 = TAB1.CONTRAINTES.2;                                                         
  129. SI3 = TAB1.CONTRAINTES.3;                                                         
  130. SI4 = TAB1.CONTRAINTES.4;                                                         
  131. SXX1 = EXTR SI1 SMXX 1 1 1;                                                     
  132. SYY1 = EXTR SI1 SMYY 1 1 1;                                                     
  133. SZZ1 = EXTR SI1 SMZZ 1 1 1;                                                     
  134. SXX2 = EXTR SI2 SMXX 1 1 1;                                                     
  135. SYY2 = EXTR SI2 SMYY 1 1 1;                                                     
  136. SZZ2 = EXTR SI2 SMZZ 1 1 1;                                                     
  137. SXX3 = EXTR SI3 SMXX 1 1 1;                                                     
  138. SYY3 = EXTR SI3 SMYY 1 1 1;                                                     
  139. SZZ3 = EXTR SI3 SMZZ 1 1 1;                                                     
  140. SXX4 = EXTR SI4 SMXX 1 1 1;                                                     
  141. SYY4 = EXTR SI4 SMYY 1 1 1;                                                     
  142. SZZ4 = EXTR SI4 SMZZ 1 1 1;                                                     
  143.  
  144. ***** RESULTATS THEORIQUES POUR CHAQUE PAS *****
  145.  
  146. sigxx1 = 7.5e6; sigxx2 = 11.66e6;  sigxx3 = 4.16e6;
  147. sigyy1 = 2.5e6; sigyy2 =  6.66e6;  sigyy3 = 4.16e6;
  148. sigzz1 = 2.5e6; sigzz2 =  6.66e6;  sigzz3 = 4.16e6;
  149.  
  150. sigxx4 = -3.33e6;
  151. sigyy4 =  1.66e6;
  152. sigzz4 =  1.66e6;
  153.  
  154.  
  155. *-----------------------------------------------------------------------        
  156. ***** AFFICHAGE DES RESULTATS                                                   
  157. *-----------------------------------------------------------------------        
  158. SAUT PAGE;                                                                      
  159. *////////// PARAMETRE D'EVOLUTION  T : 1 //////////';                           
  160. *-----------------------------------------------------------------------        
  161. MESS '                 T = 1 ';                                                 
  162. MESS ' RESULTATS THEORIQUES : SIGMAXX =  7.5 E 6' ;                             
  163. MESS '                        SIGMAYY =  2.5 E 6' ;                             
  164. MESS '                        SIGMAZZ =  2.5 E 6' ;                             
  165. MESS '                                          ' ;                             
  166. MESS ' RESULTATS CALCULES   : SIGMAXX = 'SXX1 ;                                 
  167. MESS '                        SIGMAYY = 'SYY1 ;                                 
  168. MESS '                        SIGMAZZ = 'SZZ1 ;                                 
  169. ergxx1 = abs((sigxx1 -sxx1)/sxx1);
  170. ergyy1 = abs((sigyy1 -syy1)/syy1);
  171. ergzz1 = abs((sigzz1 -szz1)/szz1);
  172.  
  173. MESS 'ergxx1 =' ergxx1;
  174. MESS 'ergyy1 =' ergyy1;
  175. MESS 'ergzz1 =' ergzz1;
  176. SAUT 1 LIGNE ;                                                                  
  177. *                                                                               
  178. *////////// PARAMETRE D'EVOLUTION  T : 2 //////////';                           
  179. *-----------------------------------------------------------------------        
  180. MESS '                 T = 2 ';                                                 
  181. MESS ' RESULTATS THEORIQUES : SIGMAXX = 11.66 E 6' ;                            
  182. MESS '                        SIGMAYY =  6.66 E 6' ;                            
  183. MESS '                        SIGMAZZ =  6.66 E 6' ;                            
  184. MESS '                                           ' ;                            
  185. MESS ' RESULTATS CALCULES   : SIGMAXX = 'SXX2 ;                                 
  186. MESS '                        SIGMAYY = 'SYY2 ;                                 
  187. MESS '                        SIGMAZZ = 'SZZ2 ;                                 
  188.  
  189. ergxx2 = abs((sigxx2 -sxx2)/sxx2);
  190. ergyy2 = abs((sigyy2 -syy2)/syy2);
  191. ergzz2 = abs((sigzz2 -szz2)/szz2);
  192.  
  193. MESS 'ergxx2 =' ergxx2;
  194. MESS 'ergyy2 =' ergyy2;
  195. MESS 'ergzz2 =' ergzz2;
  196.  
  197. SAUT 1 LIGNE ;                                                                  
  198. *                                                                               
  199. *////////// PARAMETRE D'EVOLUTION  T : 3 //////////';                           
  200. *-----------------------------------------------------------------------        
  201. MESS '                 T = 3 ' ;                                                
  202. MESS ' RESULTATS THEORIQUES : SIGMAXX =  4.16 E 6' ;                            
  203. MESS '                        SIGMAYY =  4.16 E 6' ;                            
  204. MESS '                        SIGMAZZ =  4.16 E 6' ;                            
  205. MESS '                                           ' ;                            
  206. MESS ' RESULTATS CALCULES   : SIGMAXX = 'SXX3 ;                                 
  207. MESS '                        SIGMAYY = 'SYY3 ;                                 
  208. MESS '                        SIGMAZZ = 'SZZ3 ;                                 
  209.  
  210. ergxx3 = abs((sigxx3 -sxx3)/sxx3);
  211. ergyy3 = abs((sigyy3 -syy3)/syy3);
  212. ergzz3 = abs((sigzz3 -szz3)/szz3);
  213.  
  214. MESS 'ergxx3 =' ergxx3;
  215. MESS 'ergyy3 =' ergyy3;
  216. MESS 'ergzz3 =' ergzz3;
  217.  
  218.  
  219. SAUT 1 LIGNE ;                                                                  
  220. *                                                                               
  221. *////////// PARAMETRE D'EVOLUTION  T : 4 //////////';                           
  222. *-----------------------------------------------------------------------        
  223. MESS '                 T = 4 ' ;                                                
  224. MESS ' RESULTATS THEORIQUES : SIGMAXX = -3.33 E 6' ;                            
  225. MESS '                        SIGMAYY =  1.66 E 6' ;                            
  226. MESS '                        SIGMAZZ =  1.66 E 6' ;                            
  227. MESS '                                           ' ;                            
  228. MESS ' RESULTATS CALCULES   : SIGMAXX = 'SXX4 ;                                 
  229. MESS '                        SIGMAYY = 'SYY4 ;                                 
  230. MESS '                        SIGMAZZ = 'SZZ4 ;                                 
  231. MESS '                                  ' ;                                     
  232.  
  233. ergxx4 = abs((sigxx4 -sxx4)/sxx4);
  234. ergyy4 = abs((sigyy4 -syy4)/syy4);
  235. ergzz4 = abs((sigzz4 -szz4)/szz4);
  236.  
  237. MESS 'ergxx4 =' ergxx4;
  238. MESS 'ergyy4 =' ergyy4;
  239. MESS 'ergzz4 =' ergzz4;
  240.  
  241.  
  242.  
  243.  
  244. ***** TEST DE FONCTIONNEMENT *****
  245.  
  246. lise = prog ergxx1 ergyy1 ergzz1 ergxx2 ergyy2 ergzz2
  247. ergxx3 ergyy3 ergzz3 ergxx4 ergyy4 ergzz4 ;
  248.  
  249. lima = maxi lise;
  250.  
  251. si ( lima < .05 );
  252.          erre 0;
  253. sinon;
  254.          erre 5;
  255. finsi;
  256.  
  257. lima = 100*lima;
  258. MESS 'L ECART MAXIMUM EST DE ' lima '%';
  259. temps;
  260. fin;
  261.  
  262.  
  263.  
  264.  
  265.  
  266.  
  267.  
  268.  
  269.  
  270.  
  271.  
  272.  
  273.  
  274.  

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