Télécharger la plaquette Cast3M


fr en

Télécharger thpl5.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  thpl5.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. MESS '*********************************************************';               
  14. MESS '* COMMISSION V.P.C.S.  : THERMO-PLASTICITE              *';               
  15. MESS '*                                                       *';               
  16. MESS '*CYLINDRE SOUMIS A UNE REPARTITION LOGARITHMIQUE        *';               
  17. MESS '* DE TEMPERATURE.                                       *';               
  18. MESS '*CALCUL EN DEFORMATIONS PLANES                          *';               
  19. MESS '*                                                       *';               
  20. MESS '*FONCTIONS TESTEES : ELASTO-PLASTICITE AVEC ECROUISSAGE *';               
  21. MESS '*  SOUS CHARGEMENT THERMIQUE.                           *';               
  22. MESS '*                                                       *';               
  23. MESS '*********************************************************';               
  24. *                                                                               
  26. OPTI MODE PLAN DEFO;                                                            
  27. OPTI ELEM QUA8;                                                                 
  28. OPTI TRAC BENS;                                                                 
  29. GRAPH='N';                                                                    
  30.  
  31. *++++++++++++++++++++++++++++++++*                                              
  32. *    GEOMETRIE DU CYLINDRE       *                                              
  33. *++++++++++++++++++++++++++++++++*                                              
  34. *                                                                               
  35. DENS 1. ; XA = 1. ; YA = 0. ; PA = XA YA;                                       
  36. DENS 1.2; XB = 2. ; YB = 0. ; PB = XB YB;                                       
  37. DENS 1. ; XC = 0. ; YC = 1. ; PC = XC YC;                                       
  38. DENS 1.2; XD = 0. ; YD = 2. ; PD = XD YD;                                       
  39. DENS 1. ; XO = 0. ; YO = 0. ; PO = XO YO ;                                      
  40. DENS 1.1; pz = 0. 1.497; pm = 1.497 0. ;                                        
  41.  
  42.  
  43. DAm = PA droi -7 pm;                                                            
  44. DmB = Pm DROI -3 PB;                                                            
  45. DDz = PD DROI -3 Pz;                                                            
  46. DzC = Pz DROI -7 pc;                                                            
  47. DAB = DAm et DmB;                                                               
  48. DDC = DDz et DzC;                                                               
  49. * 
  50. si complet;                                                                              
  51. N2 = 20;   
  52. sinon;
  53. n2 = 10;
  54. finsi;                                                                     
  55. SUR1 = DAB ROTA N2 90 PO;                                                       
  56.  
  57. *                                                                               
  58. NOSUR1 = NBNO SUR1;                                                             
  59. ELSUR1 = NBEL SUR1;                                                             
  60. *                                                                               
  61. ELIM (SUR1 ET DAB ET DDC) 0.0001;                                               
  62. *                                                                               
  63. SI(NEG GRAPH 'N');                                                              
  64.   TITR ' SSNP05 : MAILLAGE PAR QUA8 (NOEUDS)';                                  
  65.   TRACE NOEUD SUR1;                                                             
  66.   TITR ' SSNP05 : MAILLAGE PAR QUA8 (ELEMENTS)';                                
  67.   TRACE ELEM SUR1;                                                              
  68. FINSI;                                                                          
  69. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  70. *                  CONDITIONS AUX LIMITES                     *                 
  71. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  72. *                                                                               
  73. CL1 = SYMT DEPL ROTA PA PB SUR1 0.0001;                                         
  74. CL2 = SYMT DEPL ROTA PC PD SUR1 0.0001;                                         
  75. CL = CL1 ET CL2 ;                                                               
  76. *                                                                               
  77. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  78. *         MODELES - MATERIAUX - CARACTERISTIQUES                  
  79. *+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                 
  80. *                                                                               
  81. MO = MODE sur1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'PARFAIT' QUA8;              
  82.  
  83. MAT = MATE MO YOUN 200000 NU 0.3 ALPHA 12E-6 SIGY 300 'TALP' 0. 'TREF' 0.;                          
  84. *                                                                               
  85. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  86. *         CHAMPS RELATIFS A LA PROCEDURE                       *                
  87. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  88. *                                                                               
  89. *  CHAMPS DE TEMPERATURE REPARTIE SUR SUR1 AUX DIFFERENTS                       
  90. *     INSTANTS DE CHARGEMENT                                                    
  91. *        10 PAS DE CHARGEMENT ENTRE TE = 0°C ET TE = 200°C                      
  92. *         (TE = TEMPERATURE EXTERIEURE).                                        
  93. *                                                                               
  94. *==================================================================*            
  95. *  TABLE DE TRAVAIL : CHPOINT DE TEMPERATURE A T0 = 0 ET T1 = 1.   *            
  96. *==================================================================*            
  97. *                                                                               
  98. T0 = 0.; T1 = 1. ; T2 = 2.; T3 = 3.; T4 = 4.; T5 = 5.;                          
  99. T6 = 6.; T7 = 7. ; T8 = 8.; T9 = 9.; T10 = 10.;                                 
  100.  
  101. *                                                                               
  102. * EXTRACTIONS DES CHPOINTS DES COORDONNEES X ET Y DES NOEUDS DE SUR1            
  103. *                                                                               
  104. CHPSURX = COOR 1 SUR1;                                                          
  105. CHPSURY = COOR 2 SUR1;                                                          
  106. CHPSUR = ((CHPSURX ** 2) + (CHPSURY ** 2)) ** 0.5;                              
  107. *                                                                               
  108. * CREATION DU CHPOINT DE TEMPERATURE POUR TE = 1°C                              
  109. *                                                                               
  110. CONST1 = 1. / (LOG 2.);                                                         
  111. CHP2 = (LOG CHPSUR) * CONST1;                                                   
  112. *                                                                               
  113. * TRANSFORMATION DU SOUS-TYPE "SCAL" EN SOUS-TYPE "T" POUR CHP2                 
  114. *    UTILISATION DE L'OPERATEUR VARI                                            
  115. *                                                                               
  116. LIS10 = PROG 0. 200.;                                                           
  117. LIS20 = PROG 0. 200.;                                                           
  118. EVOL2 = EVOL MANU SCAL LIS10 T LIS20;                                           
  119. *                                                                               
  120. CHP1 = VARI CHP2 EVOL2 T;                                                       
  121. *                                                                               
  122. * CHARGEMENT DE LA TABLE DES CHPOINTS DE TEMPERATURE                            
  123. *
  124. TEMP = TABLE;WCH1 = TABLE;
  125. TEMP . 0 = T0;TEMP . 1 = T1;
  126. TEMP . 2 = T2;TEMP . 3 = T3;
  127. TEMP . 4 = T4;TEMP . 5 = T5;
  128. TEMP . 6 = T6;TEMP . 7 = T8;
  129. TEMP . 8 = T8;TEMP . 9 = T9;
  130. TEMP . 10 = T10;
  131.  
  132. WCH1.0 = MANU CHPO SUR1 1 T 0.  ;                                                                                
  133. WCH1.1  = CHP1 * 120.;                                                         
  134. WCH1.2  = CHP1 * 130.;                                                         
  135. WCH1.3  = CHP1 * 140.;                                                         
  136. WCH1.4  = CHP1 * 150.;                                                         
  137. WCH1.5  = CHP1 * 160.;                                                         
  138. WCH1.6  = CHP1 * 170.;                                                         
  139. WCH1.7  = CHP1 * 180.;                                                         
  140. WCH1.8  = CHP1 * 190.;                                                         
  141. WCH1.9  = CHP1 * 195.;                                                         
  142. WCH1.10 = CHP1 * 200.; 
  143.  
  144. CHA2 = CHAR 'T' TEMP WCH1;                                                        
  145. *                                                                               
  146. *                                                                               
  147. *                                                                               
  148. * LISTE DE VALEURS DE TEMPS OU L'ON VEUT DES RESULTATS : LIS1                   
  149. * 
  150. si complet;                                                                              
  151. LIS1 = PROG T0 T1 T2 T3 T4 T5 T6 T7 T8 T9 T10;   
  152. sinon;
  153. LIS1 =    PROG T0 T1 T2;
  154. finsi;
  155.  
  156. *                                                                               
  157. *   CREATION D'UN CHAMPS DE FORCE NUL : THERMO-ELASTICITE                       
  158. *
  159. si complet;                                                                               
  160. LFORCE = PROG 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. ;                                
  161. sinon;
  162. LFORCE = PROG 0. 0. 0.;
  163. finsi;
  164.  
  165. FOR1 = FORCE FX 1. PA;                                                          
  166. CHPO1 = RESU FOR1;                                                              
  167. EVOL1 = EVOL MANU TEMPS LIS1 FORCE LFORCE;                                      
  168. CHA1 = CHAR 'MECA' CHPO1 EVOL1;                                                        
  169. *                                                                               
  170. * 
  171. TAB1 = TABLE;  
  172. TAB1.'MODELE' = MO ;
  173. TAB1.'CARACTERISTIQUES' = MAT;
  174. TAB1.'CHARGEMENT' = (CHA1 ET CHA2);
  175. TAB1.'TEMPS_CALCULES' = LIS1;
  176. TAB1.'TEMPS_SAUVES' = LIS1;
  177. TAB1.'BLOCAGES_MECANIQUES' = CL;
  178. TAB1.'PROCESSEURS' = 'AUTOMATIQUE';                                                        
  179. *                                                                               
  180. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  181. *               PROCEDURE PASAPAS                              *                
  182. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  183. *                                                                               
  184. PASAPAS TAB1;                                           
  185. *                                                                               
  186. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  187. *         EXTRACTION et AFFICHAGE des RESULTATS                *                
  188. *++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*                
  189. *                                                                               
  190. MESS ' RESULTATS DU TEST SSNP05';                                               
  191. MESS ' ------------------------';                                               
  192. MESS ' Les resultats de reference sont la moyenne des calculs';                 
  193. MESS ' effectues par Framasoft et le Cetim';                                    
  195.  
  196. *                                                                               
  197. MESS ' Contraintes interpolees (MPa)';                                          
  198. MESS ' -----------------------------';                                          
  200. *                                                                               
  201. * Interpolations des contraintes aux noeuds                                     
  202. *
  203. si complet;
  204. NN10 = 10; 
  205. sinon;
  206. nn10 = 2;
  207. finsi;                                                                              
  208. CHEL1 = TAB1.CONTRAINTES.NN10;                                                      
  209. CHPO1 = CHAN CHPO MO CHEL1;                                                     
  210.  
  211. * Extraction des contraintes                                                    
  212. * 
  213. si complet ;   
  214. SMRR10 = 0.1731;     SMRR30 = 0.3806;      SMRR20 = 47.04;        
  215. SMTT10 = 340.6;   SMTT30 = -191.75 ;       SMTT20 = -50.63; 
  216. SMZZ10 = 114.1;  SMZZ30 = -345.95;   SMZZ20 = -284.  ;                                                     
  217. sinon;
  218. SMRR10 = .36175;    SMRR30 = -2.9683 ; SMRR20 = 35.943;
  219. SMTT10 =272.68;    SMTT30 = -169.49; SMTT20 =  -21.166;
  220. SMZZ10 = 81.913;  SMZZ30 = -348.10;   SMZZ20 = -177.17;  
  221. finsi;
  222. SMRR1 = EXTR CHPO1 SMXX PA;                                 
  223. SMRR3 = EXTR CHPO1 SMXX PB;                                  
  224. SMRR2 = EXTR CHPO1 SMXX pm;                                     
  225. *                                                                               
  226. SMTT1 = EXTR CHPO1 SMYY PA;                                 
  227. SMTT3 = EXTR CHPO1 SMYY PB;                                
  228. SMTT2 = EXTR CHPO1 SMYY pm;                                    
  229. * ;                                                                             
  230. SMZZ1 = EXTR CHPO1 SMZZ PA;                                    
  231. SMZZ3 = EXTR CHPO1 SMZZ PB;                                    
  232. SMZZ2 = EXTR CHPO1 SMZZ pm;                                     
  233. *                                                                               
  234. * Calcul des ecarts avec resultats de reference                                 
  235. *                                                                               
  236. ERGRR1 = 100 * (ABS ((SMRR1 - SMRR10) / SMRR10));                               
  237. ERGRR2 = 100 * (ABS ((SMRR2 - SMRR20) / SMRR20));                               
  238. ERGRR3 = 100 * (ABS ((SMRR3 - SMRR30) / SMRR30));                               
  239. *                                                                               
  240. ERGTT1 = 100 * (ABS ((SMTT1 - SMTT10) / SMTT10));                               
  241. ERGTT2 = 100 * (ABS ((SMTT2 - SMTT20) / SMTT20));                               
  242. ERGTT3 = 100 * (ABS ((SMTT3 - SMTT30) / SMTT30));                               
  243. *                                                                               
  244. ERGZZ1 = 100 * (ABS ((SMZZ1 - SMZZ10) / SMZZ10));                               
  245. ERGZZ2 = 100 * (ABS ((SMZZ2 - SMZZ20) / SMZZ20));                               
  246. ERGZZ3 = 100 * (ABS ((SMZZ3 - SMZZ30) / SMZZ30));                               
  247. *                                                                               
  248. MESS ' Resultats pour r = 1. mm ';                                              
  249. MESS ' ------------------------ ';                                              
  251. *                                                                               
  252. MESS ' Contrainte radiale de reference (smrr)  :' SMRR10 'MPa';                 
  253. MESS ' Contrainte radiale calculee  (smrr)     :' SMRR1 'MPa';                  
  254. MESS '            Soit un ecart de :' ERGRR1 '%';                               
  255. MESS '            ECART NON-SIGNIFICATIF';                                      
  257. *                                                                               
  258. MESS ' Contrainte tangentielle de reference     :' SMTT10 'MPa';                
  259. MESS ' Contrainte tangentielle calculee  (smtt) :' SMTT1 'MPa';                 
  260. MESS '            Soit un ecart de :' ERGTT1 '%';                               
  262. *                                                                               
  263. MESS ' Contrainte axiale de reference           :' SMZZ10 'MPa';                
  264. MESS ' Contrainte axiale calculee  (smzz)       :' SMZZ1 'MPa';                 
  265. MESS '            Soit un ecart de :' ERGZZ1 '%';                               
  267. *                                                                               
  268. MESS ' Resultats pour r = 1.497 mm ';                                           
  269. MESS ' ---------------------------- ';                                          
  271. *                                                                               
  272. MESS ' Contrainte radiale de reference          :' SMRR20 'MPa';                
  273. MESS ' Contrainte radiale calculee  (smrr)      :' SMRR2 'MPa';                 
  274. MESS '            Soit un ecart de :' ERGRR2 '%';                               
  276. *                                                                               
  277. MESS ' Contrainte tangentielle de reference     :' SMTT20 'MPa';                
  278. MESS ' Contrainte tangentielle calculee  (smtt) :' SMTT2 'MPa';                 
  279. MESS '            Soit un ecart de :' ERGTT2 '%';                               
  281. *                                                                               
  282. MESS ' Contrainte axiale de reference           :' SMZZ20 'MPa';                
  283. MESS ' Contrainte axiale calculee  (smzz)       :' SMZZ2 'MPa';                 
  284. MESS '            Soit un ecart de :' ERGZZ2 '%';                               
  286.  
  287. *                                                                               
  288. MESS ' Resultats pour r = 2. mm ';                                              
  289. MESS ' ------------------------ ';                                              
  291. *                                                                               
  292. MESS ' Contrainte radiale de reference          :' SMRR30 'MPa';                
  293. MESS ' Contrainte radiale calculee  (smrr)      :' SMRR3 'MPa';                 
  294. MESS '            Soit un ecart de :' ERGRR3 '%';                               
  295. MESS '            ECART NON-SIGNIFICATIF';                                      
  297. *                                                                               
  298. MESS ' Contrainte tangentielle de reference     :' SMTT30 'MPa';                
  299. MESS ' Contrainte tangentielle calculee  (smtt) :' SMTT3 'MPa';                 
  300. MESS '            Soit un ecart de :' ERGTT3 '%';                               
  302. *                                                                               
  303. MESS ' Contrainte axiale de reference           :' SMZZ30 'MPa';                
  304. MESS ' Contrainte axiale calculee  (smzz)       :' SMZZ3 'MPa';                 
  305. MESS '            Soit un ecart de :' ERGZZ3 '%';                               
  307. *                                                                               
  309. MESS ' Deplacements radiaux';                                                   
  310. MESS ' --------------------';                                                   
  312. *                                                                               
  313. * Extractions des deplacements radiaux                                          
  314. *                                                                               
  315. DEP1 = TAB1.DEPLACEMENTS.NN10;                                                       
  316. *                                                                               
  317. UR1 = EXTR DEP1 UX PA;                                                          
  318. UR2 = EXTR DEP1 UX PB;                                                          
  319. UR3 = EXTR DEP1 UX pm; 
  320. mess ur1 ur2 ur3;                                                         
  321. *
  322. si complet;                                                                               
  323. UR10 = 0.1804E-2;                                                               
  324. UR20 = 0.386E-2;
  325. sinon;
  326. ur10 =  1.24001E-03  ;ur20 = 2.48285E-03;
  327. finsi;                                                              
  328. *UR30 = NON CALCULE                                                             
  329. *                                                                               
  330. * Calcul des ecarts avec solution V.P.C.S.                                      
  331. *                                                                               
  332. ERGUR1 = 100 * ( ABS (( UR1 - UR10) / UR10 ));                                  
  333. ERGUR2 = 100 * ( ABS (( UR2 - UR20) / UR20 ));                                  
  334. *ERGUR3 = 100 * ( ABS (( UR3 - UR30) / UR30 ));                                 
  335. *                                                                               
  336. MESS ' Deplacement de reference radial pour r = 1. mm :' UR10 'mm';             
  337. MESS ' Deplacement calcul  radial  pour r = 1. mm    :' UR1 'mm';              
  338. MESS '    Soit un ecart de :' ERGUR1 '%';                                       
  339. SAUT 1 LIGN;;                                                                   
  340. *                                                                               
  341. *MESS ' Deplacement de reference radial pour r = 1.497 mm :' UR30 'mm';         
  342. MESS ' Deplacement calcul      radial pour r = 1.497 mm :' UR3 'mm';           
  343. *MESS '    Soit un ecart de :' ERGUR3 '%';                                      
  345. *                                                                               
  346. MESS ' Deplacement de reference radial pour r = 2. mm :' UR20 'mm';             
  347. MESS ' Deplacement calcul      radial pour r = 2. mm :' UR2 'mm';              
  348. MESS '    Soit un ecart de :' ERGUR2 '%';                                       
  350. *                                                                               
  351. MESS ' Nombre de noeuds  : ' NOSUR1;                                            
  352. MESS ' Nombre d elements : ' ELSUR1;                                            
  354. *                                                                               
  355. * Tracer de la deformee.                                                        
  356. *                                                                               
  357. SI (NEG GRAPH 'N');                                                             
  358.   DEF0 = DEFO SUR1 DEP1 0. ;                                                    
  359.   DEF1 = DEFO SUR1 DEP1 10. ROUGE ;                                             
  360.   TITR ' SSNP05 : DEFORMEE (DEFO. PLANE)';                                      
  361.   TRAC (DEF0 ET DEF1);                                                          
  362. FINSI;                                                                          
  363. *                                                                               
  364. * tracer de la contrainte equivalente de Von-Mises                              
  365. *                                                                               
  366. SI (NEG GRAPH 'N');                                                             
  367.    CHEL2 = VMIS mo CHEL1;                                                       
  368.    TITR ' TRACE DES CONTRAINTES DE VON-MISES';                                  
  369.    TRAC CHEL2 mo SUR1;                                                          
  370.    LLIS2 = PROG 0. 290. 295. 300. 305. 310. 600.;                               
  371.    TITR 'LIEU OU IL Y A PLASTIFICATION';                                        
  372.    TRAC LLIS2 CHEL2 mo SUR1;                                                    
  373. FINSI;                                                                          
  374. *                                                                               
  376. TEMPS;                                                                          
  378. *                                                                               
  379. *===================================================================*           
  380. *                  Code fonctionnement                              *           
  381. *===================================================================*           
  382. *                                                                               
  383. * Maximum des erreurs relatives                                                 
  384. *                                                                               
  385. ERG1 = MAXI (PROG ergtt1 ergzz1 ergrr2 ergtt2 ergzz2 ergtt3 ergzz3);            
  386. *                                                                               
  387. * remarque : on ne tient pas compte des ecarts sur la contrainte                
  388. *           radiale pour r = 1.mm et r = 2.mm                                   
  389. *                                                                               
  390. *                                                                               
  391. ERG2 = MAXI (PROG ERGur1 ergur2);                                               
  392. ERGMAX = MAXI (PROG ERG1 ERG2);                                                 
  393. *                                                                               
  394. * Il y a erreur si un des pourcentage erg*** est suprieur a 5%.                 
  395. *                                                                               
  396. SI (ERGMAX <EG 5.);                                                             
  397.      ERRE  0 ;                                                                  
  398. SINON;                                                                          
  399.      ERRE  5 ;                                                                  
  400. FINSI;                                                                          
  401. *                                                                               
  402. FIN;                                                                            
  403. *                                                                               
  404.  
  405.  
  406.  
  407.  
  408.  
  409.  
  410.  
  411.  
  412.  
  413.  
  414.  
  415.  
  416.  
  417.  
  418.  
  419.  

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