Télécharger la plaquette Cast3M


fr en

Télécharger thpl3.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  thpl3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 92/03/20     PROVENANCE : TC1                                     
  5. * CAS TEST DU 92/03/19     PROVENANCE : PHIL                                    
  6. opti echo 0;                                                                    
  7. MESS ' ***************************************************************';        
  8. MESS ' * COMMISSION V.P.C.S.  : THERMO-PLASTICITE                    *';        
  9. MESS ' *                                                             *';        
  10. MESS ' * CYLINDRE SOUMIS A UNE REPARTITION LOGARITHMIQUE DE          *';        
  11. MESS ' * TEMPERATURE, CALCUL EN AXISYMETRIE.                         *';        
  12. MESS ' *                                                             *';        
  13. MESS ' * FONCTIONS TESTEES : ELASTO-PLASTICITE PARFAITE SOUS         *';        
  14. MESS ' *                     CHARGEMENT THERMIQUE.                   *';        
  15. MESS ' *                                                             *';        
  16. MESS ' ***************************************************************';        
  18. MESS ' LES RESULTATS DE REFERENCE SONT LA MOYENNE DES CALCULS';                 
  19. MESS ' EFFECTUES PAR FRAMASOFT ET LE CETIM';                                    
  20.  
  21. OPTION TRACE BENS;                                                              
  22. OPTION DIME 2 ;                                                                 
  23. OPTION MODE AXIS;                                                               
  24. OPTION ELEM QUA8;                                                               
  25.  
  26. GRAPH='N';                                                                      
  27. *                                                                               
  28. *++++++++++++++++++++++++++++++++*                                              
  29. *    GEOMETRIE DU CYLINDRE       *                                              
  30. *++++++++++++++++++++++++++++++++*                                              
  31. *                                                                               
  32. DENS 1. ; RA = 1. ; ZA = 0.  ; PA = RA ZA;                                      
  33. DENS 1.2; RB = 2. ; ZB = 0.  ; PB = RB ZB;                                      
  34. dens 1.1; pm = 1.497 0;                                                         
  35.  
  36.  
  37. DENS 1. ; RC = 1. ; ZC = 0.1 ; PC = RC ZC;                                      
  38. DENS 1.2; RD = 2. ; ZD = 0.1 ; PD = RD ZD;                                      
  39. DENS 1.1; pz = 1.497 0.1;                                                       
  40. *                                                                               
  41. DAm = PA DROIT -7 pm;                                                           
  42. DmB = pm droit -3  PB;                                                          
  43. DAB = DAm et DmB;                                                               
  44.  
  45. DBD = PB DROIT -1 PD;                                                           
  46. DDz = PD DROIT -3 Pz;                                                           
  47. DCA = PC DROIT -1 PA;                                                           
  48. DzC = Pz DROIT -7 PC;                                                           
  49. DDC = DDz et DzC;                                                               
  50.  
  51. *                                                                               
  52. SUR1 = DALL DAB DBD DDC DCA PLAN ;                                              
  53.  
  54.  
  55. ELIM SUR1 0.0001;                                                               
  56. *                                                                               
  57. NOSUR1 = NBNO SUR1;                                                             
  58. ELSUR1 = NBEL SUR1;                                                             
  59. *                                                                               
  60. SI (NEG GRAPH 'N');                                                             
  61.   TITR ' SSNA04 : MAILLAGE PAR QUA8 (NOEUDS)';                                  
  62.   TRACE NOEUD SUR1;                                                             
  63.   TITR ' SSNA04 : MAILLAGE PAR QUA8 (ELEMENTS)';                                
  64.   TRACE ELEM SUR1;                                                              
  65. FINSI;                                                                          
  66.  
  67. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  68. *                  CONDITIONS AUX LIMITES                     *                 
  69. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  70. *                                                                               
  71. CL1 = BLOQUE DAB UZ ;                                                           
  72. CL2 = BLOQUE DDC UZ ;                                                           
  73. CL = CL1 ET CL2 ;                                                               
  74. *                                                                               
  75. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  76. *         MODELES - MATERIAUX - CARACTERISTIQUES                 
  77. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  78. *                                                                               
  79. MO = MODE SUR1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'PARFAIT' QUA8;              
  80.  
  81. MAT = MATE MO YOUN 200000 NU 0.3 ALPHA 12E-6 SIGY 300 'TALP' 0. 'TREF' 0.;                          
  82. *                                                                               
  83. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  84. *        CHAMPS RELATIFS A LA PROCEDURE PASAPAS                  *                
  85. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  86. *                                                                               
  87. *==============================================================*                
  88. *       CHAMPS DE TEMPERATURE REPARTIE SUR SUR1                *                
  89. *==============================================================*                
  90. *                                                                               
  91. * 10 PAS DE CHARGEMENT ENTRE 0C ET 200C                                       
  92. *                                                                               
  93. T0 = 0.; T1 = 1. ; T2 = 2.; T3 = 3. ;                                           
  94. T4 = 4.; T5 = 5. ; T6 = 6.; T7 = 7. ;                                           
  95. T8 = 8.; T9 = 9. ; T10 = 10.;                                                   
  96. *                                                                               
  97. * CREATION DU CHPOINT DE TEMPERATURE POUR TE = 1C (TEMP. EXTE.)                
  98. *                                                                               
  99. CHPSUR = COOR 1 SUR1;                                                           
  100. *                                                                               
  101. CONST1 = 1. / (LOG 2.);                                                         
  102. CHP2 = (LOG CHPSUR) * CONST1;                                                   
  103.  
  104. * TRANSFORMATION DU SOUS-TYPE "SCAL" EN SOUS-TYPE "T" POUR CHP2                 
  105. *    UTILISATION DE L'OPERATEUR VARI                                            
  106. *                                                                               
  107. LIS10 = PROG 0. 200.;                                                           
  108. LIS20 = PROG 0. 200.;                                                           
  109. EVOL2 = EVOL MANU SCAL LIS10 T LIS20;                                           
  110. *                                                                               
  111. CHP1 = VARI CHP2 EVOL2 T;                                                       
  112.  
  113. * 
  114. TEMP = TABLE;TEMPER = TABLE;
  115. TEMP . 0 = T0;TEMP . 1 = T1;
  116. TEMP . 2 = T2;TEMP . 3 = T3;
  117. TEMP . 4 = T4;TEMP . 5 = T5;
  118. TEMP . 6 = T6;TEMP . 7 = T8;
  119. TEMP . 8 = T8;TEMP . 9 = T9;
  120. TEMP . 10 = T10;
  121.  
  122. TEMPER . 0 = MANU CHPO SUR1 1 T 0.;
  123. TEMPER . 1 = (CHP1 * 20.);
  124. TEMPER . 2 = (CHP1 * 40.);
  125. TEMPER . 3 = (CHP1 * 60.);
  126. TEMPER . 4 = (CHP1 * 80.);
  127. TEMPER . 5 = (CHP1 * 100.);
  128. TEMPER . 6 = (CHP1 * 120.);
  129. TEMPER . 7 = (CHP1 * 140.);
  130. TEMPER . 8 = (CHP1 * 160.);
  131. TEMPER . 9 = (CHP1 * 180.);
  132. TEMPER . 10 = (CHP1 * 200.);
  133.  
  134. CHA2 = CHAR 'T' TEMP TEMPER;
  135.  
  136. *                                                                               
  137. *================================================================               
  138. *   LISTE DE VALEURS DE TEMPS OU L'ON VEUT DES RESULTATS : LIS1                 
  139. *================================================================               
  140. *                                                                               
  141. LIS1 = PROG T0 T1 T2 T3 T4 T5 T6 T7 T8 T9 T10;                                  
  142. *                                                                               
  143. *================================================================               
  144. *   CREATION D'UN CHAMPS DE FORCE NUL : THERMO-ELASTICITE       *               
  145. *================================================================               
  146. *                                                                               
  147. LIS2 = PROG T0 T10;                                                             
  148. LFORCE = PROG 0. 0. ;                                                           
  149. FOR1 = FORCE FZ 1. PA;                                                          
  150. CHPO1 = RESU FOR1;                                                              
  151. EVOL1 = EVOL MANU TEMPS LIS2 FORCE LFORCE;                                      
  152. CHA1 = CHAR 'MECA' CHPO1 EVOL1;                                                        
  153.  
  154. *================================================================               
  155. *           REMPLISSAGE DE LA TABLE POUR PASAPAS                 *               
  156. *================================================================               
  157. *                                                                               
  158. TAB1 = TABLE;  
  159. TAB1.'MODELE' = MO ;
  160. TAB1.'CARACTERISTIQUES' = MAT;
  161. TAB1.'CHARGEMENT' = (CHA1 ET CHA2);
  162. TAB1.'TEMPS_CALCULES' = LIS1;
  163. TAB1.'TEMPS_SAUVES' = LIS1;
  164. TAB1.'BLOCAGES_MECANIQUES' = CL;
  165.  
  166. *                                                                               
  167. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  168. *               PROCEDURE PASAPAS                               *                
  169. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  170. *                                                                               
  171. PASAPAS TAB1;                                               
  172. *                                                                               
  173. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  174. *               EXTRACTION DES RESULTATS                       *                
  175. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  176.  
  177. NN = 10;
  178. CHEL1 = TAB1.CONTRAINTES.NN;                                                      
  179.  
  180.  
  181.  
  182. CHPO1 = CHAN CHPO MO CHEL1;                                                     
  183.  
  184. CHEL2 = VMIS MO CHEL1;                                                          
  185. *                                                                               
  186. SI (NEG GRAPH 'N');                                                             
  187.    TITR ' TRACE DES CONTRAINTES DE VON-MISES';                                  
  188.    TRAC CHEL2 MO  SUR1;                                                         
  189.    LLIS2 = PROG 0. 290. 295. 300. 305. 310. 600.;                               
  190.    TITR 'LIEU OU IL Y A PLASTIFICATION';                                        
  191.    TRAC LLIS2 CHEL2 MO SUR1;                                                    
  192. FINSI;                                                                          
  193. *                                                                               
  194. SMRRA = EXTR CHPO1 SMRR PA; SMRRA1 = 0.1894;                                    
  195. ERGRRA = ABS (100 * (SMRRA - SMRRA1) / SMRRA1) ;                                
  196. SMRRB = EXTR CHPO1 SMRR PB; SMRRB1 = 0.152 ;                                    
  197. ERGRRB = ABS (100 * (SMRRB - SMRRB1) / SMRRB1) ;                                
  198. SMRRM = extr CHPO1 SMRR Pm; SMRRM1 = 47.09;                                     
  199. ERGRRM = ABS (100 * (SMRRM - SMRRM1) / SMRRM1) ;                                
  200.  
  201.  
  202. SMZZA = EXTR CHPO1 SMZZ PA; SMZZA1 = 114.2 ;                                    
  203. ERGZZA = ABS (100 * (SMZZA - SMZZA1) / SMZZA1) ;                                
  204. SMZZB = EXTR CHPO1 SMZZ PB; SMZZB1 = -345.9;                                    
  205. ERGZZB = ABS (100 * (SMZZB - SMZZB1) / SMZZB1) ;                                
  206.  SMZZm = EXTR CHPO1 SMZZ Pm; SMZZm1 = -285.1 ;                                  
  207. ERGZZm = ABS (100 * (SMZZm - SMZZm1) / SMZZm1);                                 
  208.  
  209. SMTTA = EXTR CHPO1 SMTT PA; SMTTA1 = 340.8 ;                                    
  210. ERGTTA = ABS (100 * (SMTTA - SMTTA1) / SMTTA1) ;                                
  211. SMTTB = EXTR CHPO1 SMTT PB; SMTTB1 = -191. ;                                    
  212. ERGTTB = ABS (100 * (SMTTB - SMTTB1) / SMTTB1) ;                                
  213. SMTTm = EXTR CHPO1 SMTT Pm; SMTTm1 = -51.52;                                    
  214. ERGTTm = ABS (100 * (SMTTm - SMTTm1) / SMTTm1) ;                                
  215. opti echo 0;                                                                    
  216.  
  218. MESS ' RESULTATS : TEST SSNA04/90 ';                                            
  219. MESS ' -------------------------- ';                                            
  220. MESS '    rA = 1; rM = 1.497; rB =2;';                                          
  222.  
  223. MESS ' Resultats sur les contraintes interpolees (MPa)';                        
  224. MESS ' -----------------------------------------------';                        
  226. *                                                                               
  227. MESS ' CONT.  POINT    REFERENCE        CALCULEE         ECART%';               
  228. MESS ' SMRR     A       ' SMRRA1 SMRRA  'NON SIGNIFICATIF';                     
  229. MESS ' SMRR     B       ' SMRRB1 SMRRB  'NON SIGNIFICATIF';                     
  230. MESS ' SMRR     M       ' SMRRM1 SMRRM  ERGRRM  ;                               
  232.  
  233. MESS ' SMTT     A       ' SMTTA1 SMTTA  ERGTTA  ;                               
  234. MESS ' SMTT     B       ' SMTTB1 SMTTB  ERGTTB  ;                               
  235. MESS ' SMTT     M       ' SMTTm1 SMTTm  ERGTTm  ;                               
  237. *                                                                               
  238. MESS ' SMZZ     A       ' SMZZA1 SMZZA  ERGZZA  ;                               
  239. MESS ' SMZZ     B       ' SMZZB1 SMZZB  ERGZZB  ;                               
  240. MESS ' SMZZ     M       ' SMZZm1 SMZZm  ERGZZm  ;                               
  241.  SAUT 1 LIGN;                                                                   
  242.  
  244. *                                                                               
  245. DEP1 = TAB1.DEPLACEMENTS.NN;                                                       
  246. *                                                                               
  247. SI (NEG GRAPH 'N');                                                             
  248.   TITR ' DEFORMEE : CALCUL AXISYMETRIQUE (SSNA04)';                             
  249.   DEF0 = DEFO SUR1 DEP1 0. ;                                                    
  250.   DEF1 = DEFO SUR1 DEP1 20. ROUGE ;                                             
  251.   TRAC (DEF0 ET DEF1) ;                                                         
  252. FINSI;                                                                          
  253. *                                                                               
  254. URA = EXTR DEP1 UR PA;                                                          
  255. URB = EXTR DEP1 UR PB;                                                          
  256. URM = EXTR DEP1 UR Pm;                                                          
  257. *                                                                               
  258. URA1 =  1.808E-03;                                                              
  259. URB1 =  3.859E-03;                                                              
  260. *                                                                               
  261. ERGURA = ABS (100 * (URA - URA1) / URA1) ;                                      
  262. ERGURB = ABS (100 * (URB - URB1) / URB1) ;                                      
  263. *                                                                               
  264. MESS ' Resultats sur les deplacements radiaux ';                                
  265. MESS ' -------------------------------------- ';                                
  267. *                                                                               
  268. MESS ' Deplacement radial theorique pour r = 1 mm :' URA1 'mm';                 
  269. MESS ' Deplacement radial calcul   pour r = 1 mm :' URA 'mm';                  
  270. MESS '                Soit une erreur relative de :' ERGURA '%';                
  272.  
  273. MESS ' Deplacement radial calcul pour r = 1.497 mm : ' URM 'mm';               
  275. *                                                                               
  276. MESS ' Deplacement radial thorique pour r = 2 mm :' URB1 'mm';                 
  277. MESS ' Deplacement radial calcul   pour r = 2 mm :' URB 'mm';                  
  278. MESS '                Soit une erreur relative de :' ERGURB '%';                
  279. *                                                                               
  281. *                                                                               
  282. MESS ' Nombre de noeuds  : ' NOSUR1;                                            
  283. MESS ' Nombre d elements : ' ELSUR1;                                            
  284. SAUT 1 LIGN ;                                                                   
  285. *                                                                               
  286. TEMPS;                                                                          
  287. *                                                                               
  288. SAUT 1 LIGN ;                                                                   
  289. *                                                                               
  290. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  291. *              CODE FONCTIONNEMENT                             *                
  292. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  293. *                                                                               
  294. * Maximum des erreurs relatives                                                 
  295. *                                                                               
  296. ERGMAX = MAXI (PROG ERGURA ERGURB ERGTTA ERGTTB ERGZZA ERGZZB                   
  297.          ERGRRM ERGZZM ERGTTM);                                                 
  298. * Il y a erreur si un des pourcentage erg*** est superieur a 5%                 
  299. *                                                                               
  300. SI (ERGMAX <EG 5.);                                                             
  301.      ERRE  0 ;                                                                  
  302. SINON;                                                                          
  303.      ERRE  5 ;                                                                  
  304. FINSI;                                                                          
  305. *                                                                               
  306. FIN;                                                                            
  307. *                                                                               
  308.  
  309.  
  310.  
  311.  
  312.  
  313.  
  314.  
  315.  
  316.  

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