Télécharger la plaquette Cast3M


fr en

Télécharger ther4m.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  ther4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. GRAPH='N';                                                                      
  5. SAUT PAGE ;                                                                     
  6. *-------------------------------------------------------------------*           
  7. *                                                                   *           
  8. *   TEST DES CL DE TEMPERATURE POUR PY13 CU20 TE10 PR15 TRIDIM:     *           
  9. *   -----------------------------------------------------------     *           
  10. *                                                                   *           
  11. *               TEMPERATURE IMPOSEE + CONVECTION                    *           
  12. *                                   + FLUX                          *           
  13. *                                   + SOURCE                        *           
  14. *                                                                   *           
  15. * Cet exemple permet de tester les conditions aux limites de temp{- *           
  16. * rature pour des {lements massifs 3D quadratiques.                 *           
  17. *                                                                   *           
  18. * Un parall{l{pip}de rectangle constitu{ d'{l{ments PY13, CU20, TE10*           
  19. * et PR15 est soumis @ une temp{rature impos{e @ sa face sup{rieure *           
  20. * et successivement @ une convection forc{e, @ un flux impos{ et @  *           
  21. * une source volumique de chaleur.                                  *           
  22. *                                                                   *           
  23. * Les r{sultats sont pr{sent{s sous forme d'isovaleurs qui doivent  *           
  24. * etre parall}les aux faces sup{rieure et inf{rieure.               *           
  25. *                                                                   *           
  26. *                      21 DECEMBRE 1988                             *           
  27. *-------------------------------------------------------------------*           
  28. *                                                                               
  29. *---  DISPOSITION DES NOEUDS DES ELEMENTS:                                      
  30. *                                                                               
  31. *            A4__AF4___F4___DF4___D4                                            
  32. *            /         /          /|                                            
  33. *          AB4       EF4       CD4 D34                                          
  34. *          /         /          /  |   IN:CENTRE DE                             
  35. *        B4___BE4___E4___CE4___C4  D3  (AN BN EN FN AN+1 BN+1 EN+1 FN+1)        
  36. *         |  6PY13  |   1CU20  |  /                                             
  37. *        B34       E34       C34 CD3   JN:CENTRE DE                             
  38. *         |         |          |/      (DN CN EN FN DN+1 CN+1 EN+1 FN+1)        
  39. *        B3___BE3___E3___CE3___C3                                               
  40. *                                                                               
  41. *            A3__AF3___F3___DF3___D3                                            
  42. *            /         /          /|                                            
  43. *          AB3       EF3  DE3  CD3 D23                                          
  44. *          /         /          /  |                                            
  45. *        B3___BE3___E3___CE3___C3  D2                                           
  46. *         |  1CU20  | 4PY13    |  /                                             
  47. *        B23       E23   +   C23 CD2                                            
  48. *         |         |   4TE10  |/                                               
  49. *        B2___BE2___E2___CE2___C2                                               
  50. *                                                                               
  51. *            A2__AF2___F2___DF2___D2                                            
  52. *            /         /          /|                                            
  53. *          AB2       EF2  DE2  CD2 D12                                          
  54. *          /         /          /  |                                            
  55. *        B2___BE2___E2___CE2___C2  D1                                           
  56. *         |  6PY13  |   2PR15  |  /                                             
  57. *        B12       E12       C12 CD1                                            
  58. *         |         |          |/                                               
  59. *        B1___BE1___E1___CE1___C1                                               
  60. *                                                                               
  61. *                                                                               
  62. *------------------------------------------------------------------------       
  63. *---  OPTIONS GENERALES DE CALCUL:                                              
  64. *                                                                               
  65. TITRE 'THER4 - C.L. TEMPERATURE - ELEMENTS QUADRATIQUES' ;                      
  66. OPTION  ECHO 0 ;                                                                
  67. OPTION DIME 3 ;                                                                 
  68. TEMPS ;                                                                         
  69. *                                                                               
  70. *---  CREATION DE LA GEOMETRIE: POINTS SUPPORTS DES ELEMENTS                    
  71. *                                                                               
  72.   A1=0.  0.   0.;  B1=1.  0.   0.;  C1=1.  4.   0.;  D1=0. 4.   0.;             
  73.   A2=0.  0.   2.;  B2=1.  0.   2.;  C2=1.  4.   2.;  D2=0. 4.   2.;             
  74.   A3=0.  0.   4.;  B3=1.  0.   4.;  C3=1.  4.   4.;  D3=0. 4.   4.;             
  75.   A4=0.  0.   6.;  B4=1.  0.   6.;  C4=1.  4.   6.;  D4=0. 4.   6.;             
  76.   E1=1.  2.   0.;  F1=0.  2.   0.;  E2=1.  2.   2.;  F2=0. 2.   2.;             
  77.   E3=1.  2.   4.;  F3=0.  2.   4.;  E4=1.  2.   6.;  F4=0. 2.   6.;             
  78.   I1= .5 1.   1.;  J2= .5 3.   3.;  I3= .5 1.   5.;                             
  79.  A12=0.  0.   1.; B12=1.  0.   1.; C12=1.  4.   1.; D12=0. 4.   1.;             
  80.  A23=0.  0.   3.; B23=1.  0.   3.; C23=1.  4.   3.; D23=0  4.   3.;             
  81.  A34=0.  0.   5.; B34=1.  0.   5.; C34=1.  4.   5.; D34=0. 4.   5.;             
  82.  E12=1.  2.   1.; F12=0.  2.   1.; E23=1.  2.   3.; F23=0. 2.   3.;             
  83.  E34=1.  2.   5.; F34=0.  2.   5.;                                              
  84.  AB1= .5 0.   0.; AB2= .5 0.   2.; AB3= .5 0.   4.; AB4= .5 0.   6.;            
  85.  BE1=1.  1.   0.; BE2=1.  1.   2.; BE3=1.  1.   4.; BE4=1.  1.   6.;            
  86.  CE1=1.  3.   0.; CE2=1.  3.   2.; CE3=1.  3.   4.; CE4=1.  3.   6.;            
  87.  CD1= .5 4.   0.; CD2= .5 4.   2.; CD3= .5 4.   4.; CD4= .5 4.   6.;            
  88.  DF1=0.  3.   0.; DF2=0.  3.   2.; DF3=0.  3.   4.; DF4=0.  3.   6.;            
  89.  AF1=0.  1.   0.; AF2=0.  1.   2.; AF3=0.  1.   4.; AF4=0.  1.   6.;            
  90.  EF1= .5 2.   0.; EF2= .5 2.   2.; EF3= .5 2.   4.; EF4= .5 2.   6.;            
  91.  MEC1= 1. 3. 1. ; MFD1= 0. 3. 1. ;                                              
  92. *                                                                               
  93. *---  MILIEUX DES SEGMENTS (XN IN) ET (XN JN)                                   
  94. *                                                                               
  95. AI1= .25  .5  .5; BI1= .75  .5  .5; EI1= .75 1.5  .5; FI1= .25 1.5  .5;         
  96. EJ2= .75 2.5 2.5; FJ2= .25 2.5 2.5; DJ2= .25 3.5 2.5; CJ2= .75 3.5 2.5;         
  97. AI3= .25  .5 4.5; BI3= .75  .5 4.5; EI3= .75 1.5 4.5; FI3= .25 1.5 4.5;         
  98. *                                                                               
  99. *---  MILIEUX DES SEGMENTS (XN+1 IN) ET (XN+1 JN)                               
  100. *                                                                               
  101. A2I1=.25  .5 1.5; B2I1=.75  .5 1.5; E2I1=.75 1.5 1.5; F2I1=.25 1.5 1.5;         
  102. E3J2=.75 2.5 3.5; F3J2=.25 2.5 3.5; D3J2=.25 3.5 3.5; C3J2=.75 3.5 3.5;         
  103. A4I3=.25  .5 5.5; B4I3=.75  .5 5.5; E4I3=.75 1.5 5.5; F4I3=.25 1.5 5.5;         
  104. DE2 =.5  3.  2. ; DE3 =.5  3.  4. ;                                             
  105. *                                                                               
  106. *---  CREATION DES PYRAMIDES A 13 NOEUDS:                                       
  107. *                                                                               
  108. PP01 = MANU PY13                                                                
  109.  A1   AB1  B1   BE1  E1   EF1  F1   AF1  AI1  BI1  EI1  FI1  I1 ;               
  110. PP02 = MANU PY13                                                                
  111.  E1   EF1  F1   F12  F2   EF2  E2   E12  EI1  FI1  F2I1 E2I1 I1 ;               
  112. PP03 = MANU PY13                                                                
  113.  E2   EF2  F2   AF2  A2   AB2  B2   BE2  E2I1 F2I1 A2I1 B2I1 I1 ;               
  114. PP04 = MANU PY13                                                                
  115.  A2   AB2  B2   B12  B1   AB1  A1   A12  A2I1 B2I1 BI1  AI1  I1 ;               
  116. PP05 = MANU PY13                                                                
  117.  B1   BE1  E1   E12  E2   BE2  B2   B12  BI1  EI1  E2I1 B2I1 I1 ;               
  118. PP06 = MANU PY13                                                                
  119.  A1   AF1  F1   F12  F2   AF2  A2   A12  AI1  FI1  F2I1 A2I1 I1 ;               
  120. PP07 = MANU PY13                                                                
  121.  E2   EF2  F2   F23  F3   EF3  E3   E23  EJ2  FJ2  F3J2 E3J2 J2 ;               
  122. PP08 = MANU PY13                                                                
  123.  D3   CD3  C3   C23  C2   CD2  D2   D23  D3J2 C3J2 CJ2  DJ2  J2 ;               
  124. PP09 = MANU PY13                                                                
  125.  C2   CE2  E2   E23  E3   CE3  C3   C23  CJ2  EJ2  E3J2 C3J2 J2 ;               
  126. PP10 = MANU PY13                                                                
  127.  D2   DF2  F2   F23  F3   DF3  D3   D23  DJ2  FJ2  F3J2 D3J2 J2 ;               
  128. PP11 = MANU PY13                                                                
  129.  A3   AB3  B3   BE3  E3   EF3  F3   AF3  AI3  BI3  EI3  FI3  I3 ;               
  130. PP12 = MANU PY13                                                                
  131.  E3   EF3  F3   F34  F4   EF4  E4   E34  EI3  FI3  F4I3 E4I3 I3 ;               
  132. PP13 = MANU PY13                                                                
  133.  E4   EF4  F4   AF4  A4   AB4  B4   BE4  E4I3 F4I3 A4I3 B4I3 I3 ;               
  134. PP14 = MANU PY13                                                                
  135.  A4   AB4  B4   B34  B3   AB3  A3   A34  A4I3 B4I3 BI3  AI3  I3 ;               
  136. PP15 = MANU PY13                                                                
  137.  B3   BE3  E3   E34  E4   BE4  B4   B34  BI3  EI3  E4I3 B4I3 I3 ;               
  138. PP16 = MANU PY13                                                                
  139.  A3   AF3  F3   F34  F4   AF4  A4   A34  AI3  FI3  F4I3 A4I3 I3 ;               
  140. *                                                                               
  141. *---  VOLUME COMPOSE DES 16 PYRAMIDES A 13 NOEUDS:                              
  142. *                                                                               
  143. VOL1 = PP01 ET PP02 ET PP03 ET PP04 ET PP05 ET PP06 ET PP07 ET PP08             
  144.     ET PP09 ET PP10 ET PP11 ET PP12 ET PP13 ET PP14 ET PP15 ET PP16 ;           
  145. *                                                                               
  146. *---  CREATION DES TETRAEDRES A 10 NOEUDS:                                      
  147. *                                                                               
  148. TT01 = MANU TE10                                                                
  149.  D2   CD2  C2   CE2  E2   DE2  DJ2  CJ2  EJ2  J2 ;                              
  150. TT02 = MANU TE10                                                                
  151.  D2   DF2  F2   EF2  E2   DE2  DJ2  FJ2  EJ2  J2 ;                              
  152. TT03 = MANU TE10                                                                
  153.  D3   CD3  C3   CE3  E3   DE3  D3J2 C3J2 E3J2 J2 ;                              
  154. TT04 = MANU TE10                                                                
  155.  D3   DF3  F3   EF3  E3   DE3  D3J2 F3J2 E3J2 J2 ;                              
  156. *                                                                               
  157. *---  VOLUME COMPOSE DES 4 TETRAEDRES A 10 NOEUDS:                              
  158. *                                                                               
  159. VOL2 = TT01 ET TT02 ET TT03 ET TT04 ;                                           
  160. *                                                                               
  161. *---  CREATION DES PRISMES A 15 NOEUDS:                                         
  162. *                                                                               
  163. PM01 = MANU PR15 E2 E12 E1 MEC1 C2 CE2 EF2 EF1 CD2                              
  164.                  F2 F12 F1 MFD1 D2 DF2 ;                                        
  165. PM02 = MANU PR15 E1 CE1 C1 C12 C2 MEC1 EF1 CD1 CD2                              
  166.                  F1 DF1 D1 D12 D2 MFD1 ;                                        
  167. *                                                                               
  168. *---  VOLUME COMPOSE DES 2 PRISMES A 15 NOEUDS:                                 
  169. *                                                                               
  170. VOL4 = PM01 ET PM02 ;                                                           
  171. *                                                                               
  172. *---  CREATION DES CUBES A 20 NOEUDS:                                           
  173. *                                                                               
  174. CB02 = MANU CU20                                                                
  175. E2 BE2 B2 AB2 A2 AF2 F2 EF2 E23 B23 A23 F23 E3 BE3 B3 AB3                       
  176. A3 AF3 F3 EF3 ;                                                                 
  177. CB03 = MANU CU20                                                                
  178. E3 CE3 C3 CD3 D3 DF3 F3 EF3 E34 C34 D34 F34 E4 CE4 C4 CD4                       
  179. D4 DF4 F4 EF4 ;                                                                 
  180. *                                                                               
  181. *--- VOLUME COMPOSE DES 2 CUBES A 20 NOEUDS:                                    
  182. *                                                                               
  183. VOL3 = CB02 ET CB03 ;                                                           
  184. *                                                                               
  185. *--- REUNION DES VOLUMES ELEMENTAIRES:                                          
  186. *                                                                               
  187. VOLTOT = VOL1 ET VOL2 ET VOL3 ET VOL4 ;                                         
  188. *                                                                               
  189. VOLSOU = PP11 ET PP12 ET PP13 ET PP14 ET PP15 ET PP16 ET CB03 ;                 
  190. *                                                                               
  191. *---  SURFACE INFERIEURE:                                                       
  192. *                                                                               
  193. S_INF1 = MANU QUA8 A1 AB1 B1 BE1 E1 EF1 F1 AF1 ;                                
  194. S_INF2 = MANU QUA8 C1 CD1 D1 DF1 F1 EF1 E1 CE1 ;                                
  195. S_INF = S_INF1 ET S_INF2 ;                                                      
  196. *                                                                               
  197. *---  SURFACE SUPERIEURE:                                                       
  198. *                                                                               
  199. S_SUP1 = MANU QUA8 A4 AB4 B4 BE4 E4 EF4 F4 AF4 ;                                
  200. S_SUP2 = MANU QUA8 C4 CD4 D4 DF4 F4 EF4 E4 CE4 ;                                
  201. S_SUP = S_SUP1 ET S_SUP2 ;                                                      
  202. OEIL1 = 15 12 10 ;                                                              
  203. SI(NEG GRAPH 'N');                                                              
  204.      TRAC OEIL1 VOLTOT ;                                                        
  205. FINSI;                                                                          
  206. *                                                                               
  207. *---  DONNEES DU PROBLEME DE THERMIQUE:                                         
  208. *                                                                               
  209. *---  MODELISATION:                                                             
  210. *                                                                               
  211. MOD1 = MODE VOLTOT THERMIQUE ISOTROPE ;                                         
  212. MOD2 = MODE S_SUP thermique  CONVECTION ;                                                  
  213. *                                                                               
  214. *---  DONNEE DES CARACTERISTIQUES DU MATERIAU:                                  
  215. *                                                                               
  216. MAT1 = MATE MOD1 'K' 100. ;                                                     
  217. MAT2 = MATE MOD2 'H' 15. ;                                                      
  218. *                                                                               
  219. *---  CREATION DES MATRICES DE CONDUCTIVITE:                                    
  220. *                                                                               
  221. CND1 = CONDUCTIVITE MOD1 MAT1 ;                                                 
  222. CND2 = CONDUCTIVITE MOD2 MAT2 ;                                                 
  223. *                                                                               
  224. *---  TEMPERATURES IMPOSEES: BLOQUE + DEPI                                      
  225. *                                                                               
  226. BB1 = BLOQUE S_INF 'T' ;                                                        
  227. EE1 = DEPI BB1 10. ;                                                            
  228. *                                                                               
  229. *---  FLUX EQUIVALENTS A LA CONDITION DE CONVECTION:                            
  230. *                                                                               
  231. CNV1 = CONVECTION MOD2 MAT2 'T' 200. ;  
  232. list cnv1;
  233.  
  234. *                                                                               
  235. *---  FLUX IMPOSE:                                                              
  236. *                                                                               
  237. FLU1 = FLUX MOD1 20. S_SUP ;                                                    
  238. *                                                                               
  239. *---  SOURCE VOLUMIQUE:                                                         
  240. *                                                                               
  241. S1 = SOURCE MOD1 100. VOLSOU ;                                                  
  242. *                                                                               
  243. *---  ASSEMBLAGE DES TERMES DE CONDUCTIVITE:                                    
  244. *                                                                               
  245. CCC1 = CND1 ET CND2 ET BB1 ;CCC2 = CND1 ET BB1 ;                                
  246. *                                                                               
  247. *---  ASSEMBLAGE DES TERMES DE FLUX EQUIVALENTS:                                
  248. *                                                                               
  249. FFF1 = EE1 ET CNV1 ;FFF2 = EE1 ET FLU1 ;FFF3 = EE1 ET S1 ;                      
  250. *                                                                               
  251. *---  1ERE RESOLUTION: TEMPERATURE IMPOSEE + CONVECTION                         
  252. *                                                                               
  253. CHTER1 = RESOUDRE CCC1 FFF1 ;                                                   
  254. *                                                                               
  255. *---  2EME RESOLUTION: TEMPERATURE IMPOSEE + FLUX                               
  256. *                                                                               
  257. CHTER2 = RESOUDRE CCC2 FFF2 ;                                                   
  258. *                                                                               
  259. *---  3EME RESOLUTION: TEMPERATURE IMPOSEE + SOURCE                             
  260. *                                                                               
  261. CHTER3 = RESOUDRE CCC2 FFF3 ;                                                   
  262. *                                                                               
  263. * REMARQUE : LE CHAMP DE TEMPERATURE EST SENSIBLE A LA REPARTITION              
  264. * DES POINTS D'INTEGRATION DANS LA STRUCTURE.                                   
  265. *                                                                               
  266. *---  POST-TRAITEMENT: TRACE DU MAILLAGE ET DES CHAMPS RESULTATS                
  267. *                                                                               
  268. OEIL = 100. 100. 100. ;                                                         
  269. TITRE 'ELEMENTS PY13 CU20 TE10 PR15 : MAILLAGE' ;                               
  270. SI(NEG GRAPH 'N');                                                              
  271.     TRACER OEIL VOLTOT ;                                                        
  272. FINSI;                                                                          
  273. OEIL = 100. 0. 0. ;                                                             
  274. TITRE 'ELEMENTS PY13 CU20 TE10 PR15 : CONVECTION' ;                             
  275. SI(NEG GRAPH 'N');                                                              
  276.     TRACER OEIL VOLTOT CHTER1 ;                                                 
  277. FINSI;                                                                          
  278. TITRE 'ELEMENTS PY13 CU20 TE10 PR15 : FLUX' ;                                   
  279. SI(NEG GRAPH 'N');                                                              
  280.     TRACER OEIL VOLTOT CHTER2 ;                                                 
  281. FINSI;                                                                          
  282. TITRE 'ELEMENTS PY13 CU20 TE10 PR15 : SOURCE' ;                                 
  283. SI(NEG GRAPH 'N');                                                              
  284.     TRACER OEIL VOLTOT CHTER3 ;                                                 
  285. FINSI;                                                                          
  286. *                                                                               
  287. T1 = EXTR CHTER1 T C4 ;                                                         
  288. T2 = EXTR CHTER2 T D4 ;                                                         
  289. T3 = EXTR CHTER3 T D4 ;                                                         
  290. SAUT PAGE ;                                                                     
  291. SAUT 2 LIGN ;                                                                   
  292. MESS ' TEMPERATURE MAXIMUM              DE REFERENCE        CALCULEE ';         
  293. MESS 'CONVECTION                             102.14           ' T1 ;            
  294. MESS 'FLUX                                   11.243           ' T2 ;            
  295. MESS 'SOURCE                                 20.432           ' T3 ;            
  296. TEMPS ;                                                                         
  297. *      CODE DE FONCTIONNEMENT SUR LA TEMP. MAX. EN SOURCE IMPOSEE               
  298. TREF=20.4;RESI=ABS((T3-TREF)/TREF);                                             
  299. SI(RESI <EG 5E-2);                                                              
  300.   ERRE  0;                                                                      
  301. SINO;                                                                           
  302.   ERRE  5;                                                                      
  303. FINSI;                                                                          
  304. FIN;                                                                            
  305.  
  306.  
  307.  
  308.  
  309.  
  310.  
  311.  
  312.  

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