Télécharger la plaquette Cast3M


fr en

Télécharger thpl4.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  thpl4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. complet = faux;
  5. *
  6. *************************************************
  7. *  pour calcul complet mettre complet à : vrai;
  8. ***************************************************
  9. * CAS TEST DU 92/03/20     PROVENANCE : TC1                                     
  10. * CAS TEST DU 92/03/19     PROVENANCE : PHIL                                    
  11. OPTI ECHO 0;                                                                    
  12. *                                                                               
  13. *********************************************************************           
  14. *   COMMISSION V.P.C.S.  : THERMO-PLASTICITE                        *           
  15. *                          THPL4.DGIBI                              *           
  16. *                                                                   *           
  17. *   SPHERE CREUSE EPAISSE EN THERMO-PLASTICITE                      *           
  18. *   CALCUL AXISYMETRIQUE                                            *           
  19. *                                                                   *           
  20. * FONCTIONS TESTEES : ELASTO-PLASTICITE PARFAITE SOUS               *           
  21. *       CHARGEMENT THERMIQUE.                                       *           
  22. *                                                                   *           
  23. *********************************************************************           
  24.  
  25. OPTION DIME 2 ;                                                                 
  26. OPTION MODE AXIS;                                                               
  27. OPTION TRACE BENS;                                                              
  28. GRAPH='N';                                                                      
  29. OPTION ELEM QUA8;                                                               
  30.  
  31. *                                                                               
  32. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*             
  33. *                  GEOMETRIE DU CYLINDRE                          *             
  34. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*             
  35. *                                                                               
  36. DENS 1. ; RA = 1. ; ZA = 0. ; PA = RA ZA;                                       
  37. dens 1.1; RM = 1.5; ZM = 0. ; PM1 = RM ZM;                                      
  38. DENS 1.11; RB = 2. ; ZB = 0. ; PB = RB ZB;                                      
  39. DENS 1. ; RC = 0. ; ZC = 1. ; PC = RC ZC;                                       
  40. DENS 1.1;                     PM2 = 0. 1.5;                                     
  41. DENS 1.11; RD = 0. ; ZD = 2. ; PD = RD ZD;                                      
  42. DENS 1. ; RO = 0. ; ZO = 0. ; PO = RO ZO ;                                      
  43. *                                                                               
  44.  
  45. DAM1 = PA DROI -7 PM1;   DM1B = PM1 D -4 PB; DAB = DAM1 et DM1B;                
  46. DDM2 = PD DROI -4 PM2;   DM2C = PM2 D -7 PC; DDC = DDM2 et DM2C;                
  47. *                                                                               
  48. N2 = 10;                                                                        
  49. SUR1 = DAB ROTA N2 90 PO;                                                       
  50. *                                                                               
  51. ELIM SUR1  0.0001;                                                              
  52. *                                                                               
  53. SI(NEG GRAPH 'N');                                                              
  54.   TITR ' SPHERE CREUSE EPAISSE : MAILLAGE';                                     
  55.   TRACE  SUR1;                                                                  
  56.  
  57. FINSI;                                                                          
  58.  
  59. list ( NBNO SUR1);                                                              
  60. list ( NBEL SUR1);                                                              
  61.  
  62. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  63. *                  CONDITIONS AUX LIMITES                     *                 
  64. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  65. *                                                                               
  66. CL1 = SYMT DEPL ROTA PA PB SUR1 0.0001;                                         
  67. CL2 = SYMT DEPL ROTA PC PD SUR1 0.0001;                                         
  68. CL = CL1 ET CL2 ;                                                               
  69. *                                                                               
  70. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  71. *          MODELES - MATERIAUX - CARACTERISTIQUES                  
  72. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  73. *                                                                               
  74. MO = MODE SUR1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'PARFAIT' QUA8;              
  75.  
  76. MAT = MATE MO YOUN 200e3 NU 0.3 ALPHA 1.E-5 SIGY 200 'TALP' 0. 'TREF' 0.;                          
  77. *                                                                               
  78. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*          
  79. *   LISTE DE VALEURS DE TEMPS OU L'ON VEUT DES RESULTATS : LIS1                 
  80. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*          
  81. *                                                                               
  82. *    LE NOMBRE DE PAS DE CALCUL (PAS DE CHARGEMENT), CORRESPOND AU              
  83. * NOMBRE DE VALEURS DE TEMPS OU L'ON VEUT DES RESULTATS.                        
  84. *    UN PAS ELASTIQUE (DEBUT PLASTIFICATION) JUSQU'A 100C.                     
  85. *    UN PAS PLASTIQUE TOUT LES 10C APRES.                                      
  86. *                                                                               
  87. *                                                                               
  88. T0 = 0.;    T1 = 1.;     T2 = 2.;    T3 = 3.;                                   
  89. T4 = 4.;    T5 = 5.;     T6 = 6.;    T7 = 7.;                                   
  90. T8 = 8.;    T9 = 9.;     T10 = 10.;  T11 = 11.;                                 
  91. T12 = 12.;  T13 = 13.;   T14 = 14.;  T15 = 15.;                                 
  92. T16 = 16.;  T17 = 17.;   T18 = 18.;                                             
  93. * 
  94. si complet;                                                                              
  95. LIS1 = PROG T0 T1 T2 T3 T4 T5 T6 T7 T8 T9 T10 T11 T12 T13                       
  96.               T14 T15 T16 T17 T18;
  97. sinon;
  98. LIS1 = PROG T0 T1 ;
  99. finsi;                                              
  100. *                                                                               
  101. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*           
  102. *   CHAMPS RELATIFS A LA PROCEDURE NONLIN                           *           
  103. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*           
  104. *                                                                               
  105. *  CHAMPS DE TEMPERATURE REPARTIE SUR SUR1 AUX DIFFERENTS INSTANTS              
  106. *     DE LA LISTE OU L'ON VEUT DES RESULTATS (CF. LIST1)                        
  107. *                                                                               
  108. *===================================================================*           
  109. *                                                                   *           
  110. * CALCUL D'UN CHPOINT CHP1 DE TEMPERATURE POUR UN GRADIENT DE       *           
  111. *   TEMPERATURE DE UN DEGRE CELSIUS, SOIT DT = 1C                  *           
  112. *                                                                   *           
  113. * DONNEES : - DT EST CROISSANT DE 0 @ DTMAX,                        *           
  114. *           - POUR A=1. < R < B=2., LA REPARTITION DE TEMPERATURE   *           
  115. *          EST : DT = T(A) - T(B) > 0.                              *           
  116. *                T(R) = DT/R * AB/(A-B)                             *           
  117. *                                                                   *           
  118. * CALCUL POUR DT = 1C, DE LA REPARTITION DE TEMPERATURE            *           
  119. *      SOIT      T(R) = 2. / R                                      *           
  120. *                                                                   *           
  121. *===================================================================*           
  122. *                                                                               
  123. * EXTRACTION DES CHPOINTS DES COORDONNEES R ET Z DU MAILLAGE SUR1.              
  124. *                                                                               
  125. CHPSUR1 = COOR 1 SUR1;                                                          
  126. CHPSUR2 = COOR 2 SUR1;                                                          
  127. *                                                                               
  128. * CREATION DU CHPOINT DE TEMPERATURE POUR DT = 1C                              
  129. *                                                                               
  130. CHPSUR = ( (CHPSUR1**2 ) + (CHPSUR2**2) )**(-0.5);                              
  131. CHP2 = 2. * CHPSUR;                                                             
  132. *                                                                               
  133. * TRANSFORMATION DU SOUS-TYPE "SCAL" EN SOUS-TYPE "T" POUR CHP2                 
  134. *    UTILISATION DE L'OPERATEUR VARI                                            
  135. *                                                                               
  136. LIS10 = PROG 0. 1000.;                                                          
  137. LIS20 = PROG 0. 1000.;                                                          
  138. EVOL2 = EVOL MANU SCAL LIS10 T LIS20;                                           
  139. *                                                                               
  140. CHP1 = VARI CHP2 EVOL2 T;                                                       
  141. *                                                                               
  142. * CHARGEMENT DE LA TABLE DES CHPOINTS DE TEMPERATURE                            
  143. *                                                                               
  144. *
  145. TEMP = TABLE;TEMPER = TABLE;
  146. TEMP . 0 = T0;TEMP . 1 = T1;
  147. TEMP . 2 = T2;TEMP . 3 = T3;
  148. TEMP . 4 = T4;TEMP . 5 = T5;
  149. TEMP . 6 = T6;TEMP . 7 = T8;
  150. TEMP . 8 = T8;TEMP . 9 = T9;
  151. TEMP . 10 = T10;TEMP . 11 = T11;
  152. TEMP . 12 = T12;TEMP . 13 = T13;
  153. TEMP . 14 = T14;TEMP . 15 = T15;
  154. TEMP . 16 = T16;TEMP . 17 = T17;
  155. TEMP . 18 = T18;
  156.  
  157. TEMPER . 0 =  MANU CHPO SUR1 1 T 0.;DT = 98.;
  158. TEMPER . 1 = (CHP1 * DT);DT = 110.;
  159. TEMPER . 2 = (CHP1 * DT);DT = 120.;
  160. TEMPER . 3 = (CHP1 * DT);DT = 130.;
  161. TEMPER . 4 = (CHP1 * DT);DT = 140.;
  162. TEMPER . 5 = (CHP1 * DT);DT = 150.;
  163. TEMPER . 6 = (CHP1 * DT);DT = 160.;
  164. TEMPER . 7 = (CHP1 * DT);DT = 170.;
  165. TEMPER . 8 = (CHP1 * DT);DT = 180.;
  166. TEMPER . 9 = (CHP1 * DT);DT = 188.;
  167. TEMPER . 10 = (CHP1 * DT);DT = 196.;
  168. TEMPER . 11 = (CHP1 * DT);DT = 210.;
  169. TEMPER . 12 = (CHP1 * DT);DT = 220.;
  170. TEMPER . 13 = (CHP1 * DT);DT = 230.;
  171. TEMPER . 14 = (CHP1 * DT);DT = 240.;
  172. TEMPER . 15 = (CHP1 * DT);DT = 250.;
  173. TEMPER . 16 = (CHP1 * DT);DT = 255.;
  174. TEMPER . 17 = (CHP1 * DT);DT = 264.83;
  175. TEMPER . 18 = (CHP1 * DT);
  176.  
  177. CHA2 = CHAR 'T' TEMP TEMPER;
  178.  
  179. *                                                                               
  180. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*          
  181. *   CREATION D'UN CHAMPS DE FORCE NUL : THERMO-ELASTICITE            *          
  182. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*          
  183. *                                                                               
  184. LFORCE = PROG 0. 0.;                                                            
  185. LLIS1  = PROG 0. 18.;                                                           
  186. *                                                                               
  187. FOR1 = FORCE FR 1. PA;                                                          
  188. CHPO1 = RESU FOR1;                                                              
  189. EVOL1 = EVOL MANU TEMPS LLIS1 FORCE LFORCE;                                     
  190. CHA1 = CHAR 'MECA' CHPO1 EVOL1;                                                        
  191. *                                                                               
  192. *                                                                               
  193. TAB1 = TABLE;  
  194. TAB1.'MODELE' = MO ;
  195. TAB1.'CARACTERISTIQUES' = MAT;
  196. TAB1.'CHARGEMENT' = (CHA1 ET CHA2);
  197. TAB1.'TEMPS_CALCULES' = LIS1;
  198. TAB1.'TEMPS_SAUVES' = LIS1;
  199. TAB1.'BLOCAGES_MECANIQUES' = CL;
  200.  
  201. *                                                                               
  202. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  203. *               PROCEDURES PASAPAS                               *                
  204. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  205. *                                                                               
  206. PASAPAS TAB1;                                                                                
  207.  
  208. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  209. *               EXTRACTION DES RESULTATS                       *                
  210. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  211. *                                                                               
  212. MESS ' **************************************************************';         
  213. MESS ' *   COMMISSION V.P.C.S.  : THERMO-PLASTICITE                 *';         
  214. MESS ' *               FICHIER  : SSNA06PLAPN.DGIBI                 *';         
  215. MESS ' *                                                            *';         
  216. MESS ' *   SPHERE CREUSE EPAISSE EN THERMO-PLASTICITE               *';         
  217. MESS ' *   CALCUL AXISYMETRIQUE                                     *';         
  218. MESS ' *                                                            *';         
  219. MESS ' * FONCTIONS TESTEES : ELASTO-PLASTICITE PARFAITE             *';         
  220. MESS ' *                     SOUS CHARGEMENT THERMIQUE.             *';         
  221. MESS ' *                                                            *';         
  222. MESS ' **************************************************************';         
  223. MESS ' *     CAS PLASTICITE SANS ECROUISSAGE ISOTROPE LINEAIRE      *';         
  224. MESS ' **************************************************************';         
  226. *                                                                               
  227. MESS ' Les rsultats thoriques sont tirs de l article';                       
  228. MESS ' "Thermo-elastic behaviour of strain-hardening thich-walled';             
  229. MESS ' sphere", D Durban, International journal of solids structures';          
  230. MESS ' vol 19, nx7, pp 643-652, 1983';                                          
  231. *===================================================================*           
  232. * resultats a 98 degres : indice T1                                             
  233. *===================================================================*           
  234. *                                                                               
  235. * on cherche pour r = 1.mm, 1.5mm, et 2.mm les valeurs suivantes :              
  236. *      smrr (contrainte radiale), smtt (contrainte tangentielle)                
  237. *          et ur (deplacement radial).                                          
  238. *    La temperature de 98 degres celsius correspond a la premiere               
  239. * plastification.                                                               
  240. *                                                                               
  241. NN1 = 1; NN2 = 2;
  242. NN3 = 3; NN4 = 4;
  243. NN5 = 5; NN6 = 6;
  244. NN7 = 7; NN8 = 8;
  245. NN9 = 9; NN10 = 10;
  246. NN11 = 11; NN12 = 12;
  247. NN13 = 13; NN14 = 14;
  248. NN15 = 15; NN16 = 16;
  249. NN17 = 17; NN18 = 18;                                                                                
  250.  
  251. MESS ' *------------------------------------------------*';                     
  252. MESS ' * Resultats a 98 degres : premiere plastifiction *';                     
  253. MESS ' *------------------------------------------------*';                     
  255. *                                                                               
  256. * Deplacements                                                                  
  257. *                         
  258. DEP98 = TAB1.DEPLACEMENTS.NN1;                                                       
  259. *                                                                               
  260. UR98 = TABLE; UR980 = TABLE;                                                    
  261. UR98.1 = EXTR DEP98 UR PA; UR980.1 = 1.26E-3;                                   
  262. UR98.2 = EXTR DEP98 UR PM1; UR980.2 = 2.08E-3;                                  
  263. UR98.3 = EXTR DEP98 UR PB; UR980.3 = 2.52E-3;                                   
  264. *                                                                               
  265. ERGUR98 = TABLE;                                                                
  266. I = 0;                                                                          
  267. REPETER BLOC1 3;                                                                
  268.   I = I + 1;                                                                    
  269.   SI (I &lt;EG 1) ; MESS ' r = 1.0 mm';FINSI;                                      
  270.   SI ((I &lt;EG 2) ET (I >EG 2)); MESS ' r = 1.5 mm';FINSI;                        
  271.   SI (I >EG 3); MESS ' r = 2.0 mm'; FINSI;                                      
  272. *                                                                               
  273. ERGUR98.I = 100 * ( ABS (( UR98.I - UR980.I) / (UR980.I+1e-6)));                
  274.   MESS ' Deplacement radial theorique :' UR980.I 'mm';                          
  275.   MESS ' Deplacement radial calcule   :' UR98.I 'mm';                           
  276.   MESS '    Soit un ecart de : ' ERGUR98.I '%';                                 
  277.   SAUT 1 LIGN;                                                                  
  278. FIN BLOC1;                                                                      
  279. *                                                                               
  280. * Contraintes                                                                   
  281. *                                                                               
  282. CHEL1 = TAB1.CONTRAINTES.NN1;                                                       
  283. CONT98 = CHAN CHPO MO CHEL1;                                                    
  284. *                                                                               
  285. SMRR98 = TABLE; SMRR980 = TABLE;                                                
  286. SMRR98.1 = EXTR CONT98 SMRR PA; SMRR980.1 = 0.;                                 
  287. SMRR98.2 = EXTR CONT98 SMRR PM1; SMRR980.2 = -38.52;                            
  288. SMRR98.3 = EXTR CONT98 SMRR PB; SMRR980.3 = 0.;                                 
  289. *                                                                               
  290. SMTT98 = TABLE; SMTT980 = TABLE;                                                
  291. SMTT98.1 = EXTR CONT98 SMTT PA; SMTT980.1 = -200.;                              
  292. SMTT98.2 = EXTR CONT98 SMTT PM1; SMTT980.2 = 5.93;                              
  293. SMTT98.3 = EXTR CONT98 SMTT PB; SMTT980.3 = 80.;                                
  294. *                                                                               
  295. ERGSRR98 = TABLE; ERGSTT98 = TABLE;                                             
  296. I = 0;                                                                          
  297. REPETER BLOC2 3;                                                                
  298.    I = I + 1;                                                                   
  299.    SI (I &lt;EG 1); MESS ' r = 1.0 mm';   FINSI;                                   
  300.    SI ((I &lt;EG 2) ET (I >EG 2)); MESS ' r = 1.5 mm';FINSI;                       
  301.    SI (I >EG 3); MESS ' r = 2.0 mm';   FINSI;                                   
  302. *                                                                               
  303.    ERGSRR98.I = 100 * ( ABS (( SMRR98.I - SMRR980.I ) /                         
  304.    ( SMRR980.I+1e-6 )));                                                        
  305.  
  306.    MESS '  Contrainte radiale theorique      :' SMRR980.I 'MPa';                
  307.    MESS '  Contrainte radiale calculee       :' SMRR98.I 'MPa';                 
  308.    MESS '        Soit un ecart de : ' ERGSRR98.I '%';                           
  309.    si ((( abs smrr980.i ) &lt;eg 10.) et ((abs smrr98.i) &lt;eg 10.));                
  310.        mess 'CET ECART EST NON-SIGNIFICATIF';                                   
  311.    finsi;                                                                       
  312.  
  313.    SAUT 1 LIGN;                                                                 
  314. *                                                                               
  315.    ERGSTT98.I = 100 * ( ABS (( SMTT98.I - SMTT980.I ) /                         
  316.      (  SMTT980.I+1e-6)));                                                      
  317.    MESS '  Contrainte tangentielle theorique :' SMTT980.I 'MPa';                
  318.    MESS '  Contrainte tangentielle calculee  :' SMTT98.I 'MPa';                 
  319.    MESS '        Soit un ecart de : ' ERGSTT98.I '%';                           
  320.    si ((( abs smtt980.i ) &lt;eg 10.) et ((abs smtt98.i) &lt;eg 10.));                
  321.         mess 'CET ECART EST NON-SIGNIFICATIF';                                  
  322.    finsi;                                                                       
  323.  
  324.    SAUT 1 LIGN;                                                                 
  325. FIN BLOC2;    
  326. si complet;                                                                  
  327. *                                                                               
  328. *===================================================================*           
  329. * resultats a 196 degres : indice T11                                           
  330. *===================================================================*           
  331. *                                                                               
  332. * On cherche - pour r = 1. mm et r = 2. mm smtt et smrr,                        
  333. *            - la valeur du rayon correspondant a la limite de                  
  334. *    plastification, et pour ce rayon, on veut smtt et smrr. En                 
  335. *    remarque on peut dire que smeq = 200.MPa = smtt + smrr.                    
  336. *            - la valeur du  rayon ou il y a un etat de contrainte              
  337. *    hydrostatique, et la valeur de cette contrainte.                           
  338. *                                                                               
  340. MESS ' *------------------------*';                                             
  341. MESS ' * Resultats a 196 degres *';                                             
  342. MESS ' *------------------------*';                                             
  344. *                                                                               
  345. * smtt et smrr  pour r = 1. mm et r = 2. mm                                     
  346. *                           
  347. CHEL1 = TAB1.CONTRAINTES.NN11; CHV1 = VMIS MO CHEL1;                                
  348. CONT196 = CHAN CHPO MO chel1;                                                   
  349. *                                                                               
  350. SMTT196 = TABLE; SMTT1960 = TABLE;                                              
  351. ERTT196  = TABLE;                                                               
  352. SMTT196.1  = EXTR CONT196 SMTT PA; SMTT1960.1 = -200.;                          
  353. SMTT196.2  = EXTR CONT196 SMTT PB; SMTT1960.2 = 153.82;                         
  354. *                                                                               
  355. SMRR196 = TABLE; SMRR1960 = TABLE;                                              
  356. SMRR196.1  = EXTR CONT196 SMRR PA; SMRR1960.1 = 0.;                             
  357. SMRR196.2  = EXTR CONT196 SMRR PB;                                              
  358. * SMRR1960.2 = 'NON-CONNUE';                                                    
  359. *                                                                               
  361. MESS ' r = 1.0 mm';                                                             
  362. MESS ' Contrainte radiale theorique      :' SMRR1960.1 'MPa';                   
  363. MESS ' Contrainte radiale calculee       :' SMRR196.1 'MPa';                    
  364. si ((abs SMRR196.1) &lt;eg 5.);                                                    
  365.      MESS 'ECART NON-SIGNIFICATIF';                                             
  366. finsi;                                                                          
  368. *                                                                               
  369. MESS ' Contrainte tangentielle theorique :' SMTT1960.1 'MPa';                   
  370. MESS ' Contrainte tangentielle calculee  :' SMTT196.1 'MPa';                    
  371.  ERTT196.1 = 100 * ( ABS (( SMTT196.1 - SMTT1960.1 ) /                          
  372.  (  SMTT1960.1+1e-6)));                                                         
  373. MESS '        Soit un ecart de : ' ERTT196.1 '%';                               
  375. *                                                                               
  376. MESS ' r = 2.0 mm';                                                             
  377. * MESS ' Contrainte radiale theorique :' SMRR1960.1 'MPa';                      
  378. MESS ' Contrainte radiale calculee  :' SMRR196.2 'MPa';                         
  380. *                                                                               
  381. MESS ' Contrainte tangentielle theorique :' SMTT1960.2 'MPa';                   
  382. MESS ' Contrainte tangentielle calculee  :' SMTT196.2 'MPa';                    
  383. ERTT196.2 = 100 * ( ABS (( SMTT196.2 - SMTT1960.2 ) /                           
  384.   (  SMTT1960.2+1e-6)));                                                        
  385. MESS '        Soit un ecart de : ' ERTT196.2 '%';                               
  387. *                                                                               
  388. * Recherche de la valeur du rayon de limite plastique                           
  389. *                                                                               
  390. *     On cherche cette valeur aux points de gauss pour les contraintes          
  391. *                                                                               
  392. *    Creation d'un champs CHEL2 par elements des valeurs des rayons             
  393. * (en spherique) aux points de gauss pour les contraintes. Alors le             
  394. * champs CHEL2 a la meme structure que le champs de contrainte CHEL1.           
  395. *                                                                               
  396.  
  397. chsur1 = coor 1 chel1;                                                          
  398. chsur2 = coor 2 chel1;                                                          
  399. CHEL2 = ((CHSUR1 ** 2) + (CHSUR2 ** 2)) ** 0.5;                                 
  400.  
  401.  
  402. *   Boucle sur les elements numerotes de 1 a 11 de la sous-zone 1, la           
  403. * valeur maximale etant inferieure a 1.5 mm. On affiche seulement le            
  404. * rayon maximum (RAYON1MAX).                                                    
  405. *                                                                               
  406. RAY1MAX = 1.;                                                                   
  407. IELEM = 0;                                                                      
  408. REPETER BLOC5 11;                                                               
  409.     IELEM = IELEM + 1;                                                          
  410. *                                                                               
  411. * Boucle sur les 9 points de gauss de l'element numerote IELEM                  
  412. *                                                                               
  413.     IPGAUS = 0;                                                                 
  414.     REPETER BLOC4 9;                                                            
  415.        IPGAUS = IPGAUS + 1;                                                     
  416.        SMTT1 = EXTR CHEL1 SMTT 1 IELEM IPGAUS;                                  
  417.        SMRR1 = EXTR CHEL1 SMRR 1 IELEM IPGAUS;                                  
  418.        SEQU1 = SMRR1 - SMTT1;                                                   
  419. *                                                                               
  420. *  test sur la contrainte equivalente : proche de 200 MPa ?                     
  421. *                                                                               
  422.        SI (( SEQU1 &lt;EG 205.) ET (SEQU1 >EG 195.));                              
  423. *                                                                               
  424. *    Dans ce cas, on cherche la valeur du rayon pour laquelle la contrainte     
  425. * equivalente est proche de la contrainte limite plastique.                     
  426. *                                                                               
  427.         RAYON1 = EXTRA CHEL2 SCAL 1 IELEM IPGAUS;                               
  428.         RAY1MAX = MAXI (PROG RAY1MAX RAYON1);                                   
  429.  
  430.         Si ( RAYON1 EGA RAY1MAX ) ;                                             
  431.                  SMTTR1 = SMTT1;                                                
  432.                  SMRRR1 = SMRR1;                                                
  433.                  SMVMR1 = extr chv1 'SCAL' 1 IELEM IPGAUS;                      
  434.                  xel = extr CHSUR1 'SCAL' 1 IELEM IPGAUS;                       
  435.                  yel = extr CHSUR2 'SCAL' 1 IELEM IPGAUS;                       
  436.  
  437.          finsi;                                                                 
  438.       FINSI;                                                                    
  439. *                                                                               
  440. * Fin boucle sur les points de gauss, et reinitialisation de IPGAUSS a 0.       
  441. *                                                                               
  442.    FIN BLOC4;                                                                   
  443. *                                                                               
  444. * Fin de boucle sur les elements numerotes de 1 a 11;                           
  445. *                                                                               
  446. FIN BLOC5;                                                                      
  447. RAYTH1 = 1.1328;                                                                
  448. SMTTR10 = -249.88;                                                              
  449. SMRRR10 = -49.88;                                                               
  450. erra1 = 100*(abs ( ( RAYTH1 - RAY1MAX)/RAYTH1));                                
  451. erttra1 = 100*(abs ( (SMTTR1 - smttr10)/smttr10));                              
  452. errra1 = 100*(abs ( (SMRRR1 -smrrr10)/smrrr10));                                
  453.  
  454.  
  455. * Affichage de la limite plastique de la premiere zone.                         
  456. *                                                                               
  457. MESS ' Rayon definisant la fin de la premiere zone plastique.';                 
  458. MESS ' ------------------------------------------------------';                 
  460. MESS ' Rayon1 calcul   = ' RAY1MAX ' mm';                                      
  461. MESS ' Rayon1 thorique = ' RAYTH1  ' mm';                                      
  462. MESS ' ECART           = ' erra1   '% ';                                        
  464. MESS ' SMTT POUR CE RAYON =' SMTTR1 ' MPA';                                     
  465. MESS ' SMTT thorique     =' SMTTR10 ' MPA';                                    
  466. MESS ' ECART             =' erttra1 ' MPA';                                     
  468. MESS ' SMRR POUR CE RAYON =' SMRRR1 ' MPA';                                     
  469. MESS ' SMRR thorique     =' SMRRR10 ' MPA';                                    
  470. mESS ' ECART             =' errra1 ' MPA';                                      
  472. MESS ' CONTRAINTE DE VON MISES =' SMVMR1 'MPA';                                 
  473. MESS 'LES COORDONNEES DU POINT DE GAUSS VALENT ';                               
  474. MESS ' X= ';list xel;                                                           
  475. MESS ' Y= ';list yel;                                                           
  476.  SAUT 1 LIGN;                                                                   
  477. *===================================================================*           
  478. * resultats a 260 dgres : indice T18                                            
  479. *===================================================================*           
  480. *                                                                               
  482. MESS ' *--------------------------------*';                                     
  483. MESS ' * Resultats a 264 degres celsius *';                                     
  484. MESS ' *--------------------------------*';                                     
  485. MESS ' Apparition de la deuxieme zone plastique (r = 2.)';                      
  486. MESS ' On cherche la limite de la premiere zone plastique.';                    
  487. MESS ' La valeur theorique est proche de 1.1968mm';                             
  489. *                                                                               
  490. * On cherche la limite de la premiere zone plastique.                           
  491. *                       
  492. CHEL1 = TAB1.CONTRAINTES.NN18;  CP1 = chan CHPO mo chel1;                           
  493. cVM1  = VMIS mo chel1;                                                          
  494.  
  495.  
  496. RAY2MAX = 1.;                                                                   
  497. IELEM = 0;                                                                      
  498.  
  499.  
  500. REPETER BLOC5 10;                                                               
  501.     IELEM = IELEM + 1;                                                          
  502.     IPGAUS = 0;                                                                 
  503.     REPETER BLOC4 9;                                                            
  504.       IPGAUS = IPGAUS + 1;                                                      
  505.       SMTT1 = EXTR CHEL1 SMTT 1 IELEM IPGAUS;                                   
  506.       SMRR1 = EXTR CHEL1 SMRR 1 IELEM IPGAUS;                                   
  507.       SEQU1 = SMRR1 - SMTT1;                                                    
  508. *                                                                               
  509.       SI (( SEQU1 &lt;EG 201.) ET (SEQU1 >EG 199.));                               
  510.         RAYON2 = EXTRA CHEL2 SCAL 1 IELEM IPGAUS;                               
  511.         RAY2MAX = MAXI (PROG RAY2MAX RAYON2);                                   
  512. *       MESS ' Pour r =' RAYON2 'mm , sigma equivalente =' SEQU1 'MPa';         
  513.  
  514.         SI (RAYON2 EGA RAY2MAX);                                                
  515.             SMTTR2 = SMTT1;                                                     
  516.             SMRRR2 = SMRR1;                                                     
  517.             SMVMR2 = extr CVM1 'SCAL' 1 IELEM IPGAUS;                           
  518.             xel = extr chsur1 'SCAL' 1 IELEM IPGAUS;                            
  519.             yel = extr chsur2 'SCAL' 1 IELEM IPGAUS;                            
  520.         FINSI;                                                                  
  521.  
  522.       FINSI;                                                                    
  523.    FIN BLOC4;                                                                   
  524. FIN BLOC5;                                                                      
  525. *                                                                               
  526. * Affichage de la limite plastique de la premiere zone.                         
  527. *                                                                               
  528. MESS ' Rayon definisant la fin de la premiere zone plastique.';                 
  529. MESS ' ------------------------------------------------------';                 
  531. MESS '   Rayon2 = ' RAY2MAX 'mm';                                               
  532. RAY2TH = 1.1968;                                                                
  533. erra2 = 100*(abs (ray2th -ray2max)/ray2th);                                     
  534. MESS '   RAYON2 thorique = ' RAY2TH ' mm';                                     
  535. MESS 'ECART =' erra2;                                                           
  536.  
  538. *                                                                               
  539. MESS 'AU POINT DE GAUSS,LES CONTRAINTES SONT';                                  
  540. MESS 'SMTT =' SMTTR2 'MPa ;SMRR =' SMRRR2 'MPa';                                
  541. MESS 'SMVM =' SMVMR2 'MPa';                                                     
  542. MESS 'LES COORDONNEES DU POINT VALENT';                                         
  543. MESS ' X= ' ;LIST xel;                                                          
  544. MESS ' Y= ' ;list yel;                                                          
  546.  
  548. TEMPS;                                                                          
  550. *                                                                               
  551.  
  552.  
  553. *                                                                               
  554. * TRACE DES ISO-CONTRAINTES                                                     
  555. *                                                                               
  556. CHELCONT = TABLE;                                                               
  557. CHELCONT.1 = TAB1.CONTRAINTES.NN1; 
  558. CHELCONT.10 = TAB1.CONTRAINTES.NN10;               
  559. CHELCONT.2 = TAB1.CONTRAINTES.NN2;   
  560. CHELCONT.11 = TAB1.CONTRAINTES.NN11;               
  561. CHELCONT.3 = TAB1.CONTRAINTES.NN3;   
  562. CHELCONT.12 = TAB1.CONTRAINTES.NN12;               
  563. CHELCONT.4 = TAB1.CONTRAINTES.NN4;   
  564. CHELCONT.13 = TAB1.CONTRAINTES.NN13;               
  565. CHELCONT.5 = TAB1.CONTRAINTES.NN5;   
  566. CHELCONT.14 = TAB1.CONTRAINTES.NN14;               
  567. CHELCONT.6 = TAB1.CONTRAINTES.NN6;   
  568. CHELCONT.15 = TAB1.CONTRAINTES.NN15;               
  569. CHELCONT.7 = TAB1.CONTRAINTES.NN7;   
  570. CHELCONT.16 = TAB1.CONTRAINTES.NN16;               
  571. CHELCONT.8 = TAB1.CONTRAINTES.NN8;   
  572. CHELCONT.17 = TAB1.CONTRAINTES.NN17;               
  573. CHELCONT.9 = TAB1.CONTRAINTES.NN9;   
  574. CHELCONT.18 = TAB1.CONTRAINTES.NN18;               
  575. *                                                                               
  576. * RECHERCHE DES LIEUX OU IL Y A PLASTIFICATION : TRACEE                         
  577. SI (NEG GRAPH 'N');                                                             
  578.  
  579.       CHEL2 = VMIS mo CHELCONT.1;                                               
  580.       TITR 'T= 98C : ISOCONTRAINTES DE VON-MISES';                             
  581.       TRAC CHEL2 mo SUR1;                                                       
  582. *                                                                               
  583.       CHEL2 = VMIS mo CHELCONT.11;                                              
  584.       TITR 'T= 196C : ISOCONTRAINTES DE VON-MISES';                            
  585.       TRAC CHEL2 mo SUR1;                                                       
  586. *                                                                               
  587.       CHEL2 = VMIS mo CHELCONT.18;                                              
  588.       TITR 'T= 265C : ISOCONTRAINTES DE VON-MISES';                            
  589.       TRAC CHEL2 mo SUR1;                                                       
  590. *                                                                               
  591.   LLIS10 = PROG 195. 200. 205.;                                                 
  592.       CHEL2 = VMIS mo CHELCONT.1;                                               
  593.       TITR 'T= 98C : ISOCONTRAINTES DE VON-MISES';                             
  594.       TRAC LLIS10 CHEL2 mo SUR1;                                                
  595. *                                                                               
  596.       CHEL2 = VMIS mo CHELCONT.11;                                              
  597.       TITR 'T= 196C : ISOCONTRAINTES DE VON-MISES';                            
  598.       TRAC LLIS10  CHEL2 mo SUR1;                                               
  599. *                                                                               
  600.       CHEL2 = VMIS mo CHELCONT.18;                                              
  601.       TITR 'T= 265C : ISOCONTRAINTES DE VON-MISES';                            
  602.       TRAC LLIS10 CHEL2 mo SUR1;                                                
  603. finsi;                                                                          
  605.    TEMPS    ;                                                                   
  607. finsi;                                                                               
  608.  
  609. *****************************************************                           
  610. ******   CODE DE FONCTIONNEMENT               *******                           
  611. *****************************************************                           
  612. * ON NE TIENT PAS COMPTE DES VALEURS NON-SIGNIFICATIVES  *                      
  613. * IL Y A ERREUR SI UN DES POURCENTAGES EST SUPERIEUR A 5%*                      
  614. SUERG1 = PROG ERGUR98.1 ERGUR98.2 ERGUR98.3 ERGSRR98.2                          
  615.              ERGSTT98.1 ERGSTT98.3;
  616. si complet;                                             
  617. SUERG2 = PROG ERTT196.1  ERTT196.2 ERRA1 ERTTRA1 ERRRA1                         
  618.               ERRA2;                                                            
  619. ERMA1 = MAXI SUERG1;                                                            
  620. ERMA2 = MAXI SUERG2;
  621. sinon;
  622. ERMA1 = MAXI SUERG1;                                                            
  623. ERMA2 = 0.;
  624. finsi;                                                            
  625. ERMA = MAXI ( PROG ERMA1 ERMA2);                                                
  626. SI  (ERMA &lt;EG 3.);                                                              
  627.          ERRE 0;                                                                
  628. SINON;                                                                          
  629.          ERRE 5;                                                                
  630. FINSI;                                                                          
  631. MESS 'L ECART MAXIMAL EST DE ' ERMA '%';                                        
  632. FIN;                                                                            
  633.  
  634.  
  635.  
  636.  
  637.  
  638.  
  639.  
  640.  
  641.  
  642.  
  643.  
  644.  

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