Télécharger la plaquette Cast3M


fr en

Télécharger diffu4.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : diffu4.dgibi
  2.  
  3. ************************************************************************
  4. ************************************************************************
  5.  
  6. ************************************************************************
  7. * Ce Cas-Test Vérifie le bon fonctionnement du modele de DIFFUSION dans
  8. * les éléments finis suivants : COQ3,COQ4,COQ6,COQ8,MASSIFS(3D)
  9. *    - MODE
  10. *    - MATE
  11. *    - COND
  12. *    - BLOQ
  13. *    - CAPA
  14. *
  15. * Le test consiste a comparer le CHPOINT solution obtenue entre la
  16. * THERMIQUE et la DIFFUSION
  17. ************************************************************************
  18.  
  19. OPTI 'DIME' 3;
  20. OPTI 'ELEM' 'CUB8';
  21. OPTI 'TRAC' 'X';
  22.  
  23. PRECIS = VALE 'PREC';
  24.  
  25. BLIST  = FAUX ; 'COMM' 'Mettre VRAI pour avoir les LIST des OBJETS  ';
  26. BTRAC  = FAUX ; 'COMM' 'Mettre VRAI pour avoir les TRAC des CHPOINTS';
  27.  
  28. ************************************************************************
  29. *                               MAILLAGE                               *
  30. ************************************************************************
  31.  
  32. * Nom des composantes PRIMALES et DUALES
  33. MP = 'B';
  34. MD = 'Z';
  35.  
  36. * POINTS *
  37. P1 = 0. 0. 0. ;
  38. P2 = 1. 0. 0. ;
  39.  
  40. * LIGNES *
  41. NBE1 = 10 ;
  42. L1   = DROI NBE1 P1 P2;
  43. L2   = L1 'PLUS' (0. 1. 0.);
  44.  
  45. * SURFACES *
  46. S1   = 'REGL' NBE1 L1 L2   ;
  47. S2   = S1 'PLUS' (0. 0. 1.);
  48. S3   = CHAN S2 'TRI3';
  49. P3   = POIN S3 'PROC' (0. 0. 1.);
  50. P4   = POIN S3 'PROC' (1. 1. 1.);
  51. P5   = POIN S1 'PROC' (1. 1. 0.);
  52.  
  53. S11  = CHAN S1 'QUA8';
  54. S31  = CHAN S3 'TRI6';
  55.  
  56. * VOLUMES *
  57. V1   = S1 VOLU NBE1 S2;
  58.  
  59. * COULEURS *
  60. S1   = S1  'COUL' 'BLEU' ;
  61. S11  = S11 'COUL' 'BLEU' ;
  62. S2   = S2  'COUL' 'ROUG' ;
  63. S3   = S3  'COUL' 'JAUN' ;
  64. S31  = S31 'COUL' 'JAUN' ;
  65.  
  66. * TRACE *
  67. SI BTRAC;
  68.   TRAC 'CACH' 'QUAL'  V1         ;
  69.   TRAC 'CACH' 'QUAL' (S1  ET S3 );
  70.   TRAC 'CACH' 'QUAL' (S11 ET S31);
  71. FINS;
  72.  
  73. ************************************************************************
  74. *                                MODELE                                *
  75. ************************************************************************
  76. MODTHE0 ='MODE'  V1          'THERMIQUE' 'CONDUCTION'              ;
  77. MODTHE1 ='MODE' (S1 ET S3)   'THERMIQUE' 'CONDUCTION' 'COQ3' 'COQ4';
  78. MODTHE2 ='MODE' (S11 ET S31) 'THERMIQUE' 'CONDUCTION' 'COQ6' 'COQ8';
  79.  
  80. MODDIF0 ='MODE'  V1          
  81.                    'DIFFUSION' 'FICK'               'INCO' MP MD;
  82. MODDIF1 ='MODE' (S1 ET S3)
  83.                    'DIFFUSION' 'FICK' 'COQ3' 'COQ4' 'INCO' MP MD;
  84. MODDIF2 ='MODE' (S11 ET S31)
  85.                    'DIFFUSION' 'FICK' 'COQ6' 'COQ8' 'INCO' MP MD;
  86. SI BLIST;
  87.   LIST MODTHE0;
  88.   LIST MODDIF0;
  89.   LIST MODTHE1;
  90.   LIST MODDIF1;
  91.   LIST MODTHE2;
  92.   LIST MODDIF2;
  93. FINS;
  94.  
  95. ************************************************************************
  96. *                               MATERIAU                               *
  97. ************************************************************************
  98. MATTHE0 = MATE MODTHE0 'K'  2.D0 'RHO'  1.D0 'C' 1.D0             ;
  99. MATTHE1 = MATE MODTHE1 'K'  2.D0 'RHO'  1.D0 'C' 1.D0 'EPAI' 1.D-5;
  100. MATTHE2 = MATE MODTHE2 'K'  2.D0 'RHO'  1.D0 'C' 1.D0 'EPAI' 1.D-5;
  101. MATDIF0 = MATE MODDIF0 'KD' 2.D0 'CDIF' 1.D0                      ;
  102. MATDIF1 = MATE MODDIF1 'KD' 2.D0 'CDIF' 1.D0          'EPAI' 1.D-5;
  103. MATDIF2 = MATE MODDIF2 'KD' 2.D0 'CDIF' 1.D0          'EPAI' 1.D-5;
  104.  
  105. SI BLIST;
  106.   LIST MATTHE0;
  107.   LIST MATDIF0;
  108.   LIST MATTHE1;
  109.   LIST MATDIF1;
  110.   LIST MATTHE2;
  111.   LIST MATDIF2;
  112. FINS;
  113.  
  114. ************************************************************************
  115. *                             CONDUCTIVITE                             *
  116. ************************************************************************
  117. CONDTH0 = COND MODTHE0 MATTHE0;
  118. CONDDI0 = COND MODDIF0 MATDIF0;
  119. CONDTH1 = COND MODTHE1 MATTHE1;
  120. CONDDI1 = COND MODDIF1 MATDIF1;
  121. CONDTH2 = COND MODTHE2 MATTHE2;
  122. CONDDI2 = COND MODDIF2 MATDIF2;
  123.  
  124. SI BLIST;
  125.   LIST CONDTH0;
  126.   LIST CONDDI0;
  127.   LIST CONDTH1;
  128.   LIST CONDDI1;
  129.   LIST CONDTH2;
  130.   LIST CONDDI2;
  131. FINS;
  132.  
  133.  
  134. ************************************************************************
  135. *                             BLOCAGES                                 *
  136. ************************************************************************
  137. BLTH1 = BLOQ 'T' P1;
  138. DDTH1 = DEPI BLTH1 0.D0;
  139. BLTH2 = BLOQ 'T' P5;
  140. DDTH2 = DEPI BLTH2 1.D0;
  141. BLTH3 = BLOQ 'T' P3;
  142. DDTH3 = DEPI BLTH3 0.D0;
  143. BLTH4 = BLOQ 'T' P4;
  144. DDTH4 = DEPI BLTH4 1.D0;
  145.  
  146. BLDI1 = BLOQ MP P1;
  147. DDDI1 = DEPI BLDI1 0.D0;
  148. BLDI2 = BLOQ MP P5;
  149. DDDI2 = DEPI BLDI2 1.D0;
  150. BLDI3 = BLOQ MP P3;
  151. DDDI3 = DEPI BLDI3 0.D0;
  152. BLDI4 = BLOQ MP P4;
  153. DDDI4 = DEPI BLDI4 1.D0;
  154.  
  155. ************************************************************************
  156. *                             RESOLUTION                               *
  157. ************************************************************************
  158. RIGTH0 = CONDTH0 ET BLTH1 ET BLTH2 ET BLTH3 ET BLTH4;
  159. SECTH0 = DDTH1 ET DDTH2 ET DDTH3 ET DDTH4;
  160. T0 = RESO RIGTH0 SECTH0;
  161.  
  162. RIGTH1 = CONDTH1 ET BLTH1 ET BLTH2 ET BLTH3 ET BLTH4;
  163. SECTH1 = DDTH1 ET DDTH2 ET DDTH3 ET DDTH4;
  164. T1 = RESO RIGTH1 SECTH1;
  165.  
  166. RIGTH2 = CONDTH2 ET BLTH1 ET BLTH2 ET BLTH3 ET BLTH4;
  167. SECTH2 = DDTH1 ET DDTH2 ET DDTH3 ET DDTH4;
  168. T2 = RESO RIGTH2 SECTH2;
  169.  
  170. RIGDI0 = CONDDI0 ET BLDI1 ET BLDI2 ET BLDI3 ET BLDI4;
  171. SECDI0 = DDDI1 ET DDDI2 ET DDDI3 ET DDDI4;
  172. D0 = RESO RIGDI0 SECDI0;
  173.  
  174. RIGDI1 = CONDDI1 ET BLDI1 ET BLDI2 ET BLDI3 ET BLDI4;
  175. SECDI1 = DDDI1 ET DDDI2 ET DDDI3 ET DDDI4;
  176. D1 = RESO RIGDI1 SECDI1;
  177.  
  178. RIGDI2 = CONDDI2 ET BLDI1 ET BLDI2 ET BLDI3 ET BLDI4;
  179. SECDI2 = DDDI1 ET DDDI2 ET DDDI3 ET DDDI4;
  180. D2 = RESO RIGDI2 SECDI2;
  181.  
  182. SI BTRAC;
  183.   TRAC (T0 ET D0)  V1         ;
  184.   TRAC (T1 ET D1) (S1  ET S3 );
  185.   TRAC (T2 ET D2) (S11 ET S31);
  186. FINS;
  187.  
  188.  
  189. ************************************************************************
  190. *                               CAPACITE                               *
  191. ************************************************************************
  192. CAPATH0 = CAPA MODTHE0 MATTHE0;
  193. CAPATH1 = CAPA MODTHE1 MATTHE1;
  194. CAPATH2 = CAPA MODTHE2 MATTHE2;
  195. CAPADI0 = CAPA MODDIF0 MATDIF0;
  196. CAPADI1 = CAPA MODDIF1 MATDIF1;
  197. CAPADI2 = CAPA MODDIF2 MATDIF2;
  198.  
  199. SI BLIST;
  200.   LIST CAPATH0;
  201.   LIST CAPADI0;
  202.   LIST CAPATH1;
  203.   LIST CAPADI1;
  204.   LIST CAPATH2;
  205.   LIST CAPADI2;
  206. FINS;
  207.  
  208.  
  209. ************************************************************************
  210. *                            COMPARAISONS                              *
  211. ************************************************************************
  212. LMOT0a = MOTS 'LX' 'B';
  213. LMOT0b = MOTS 'LX' 'T';
  214. D0b    = CHAN 'COMP' LMOT0a LMOT0b D0;
  215. T0c    = EXCO 'T' T0 ;
  216. D0c    = EXCO 'T' D0b;
  217. RESU0  = RESU (ABS (T0c - D0c));
  218. MAX0   = MAXI RESU0;
  219.  
  220. LMOT1a = MOTS 'LX' 'BIN ' 'B   ' 'BSU ';
  221. LMOT1b = MOTS 'LX' 'TINF' 'T   ' 'TSUP';
  222. D1b    = CHAN 'COMP' LMOT1a LMOT1b D1;
  223. T1c    = EXCO (MOTS 'TINF' 'T   ' 'TSUP') T1 ;
  224. D1c    = EXCO (MOTS 'TINF' 'T   ' 'TSUP') D1b;
  225. RESU1  = RESU (ABS (T1c - D1c));
  226. MAX1   = MAXI RESU1;
  227.  
  228. LMOT2a = MOTS 'LX' 'BIN ' 'B   ' 'BSU ';
  229. LMOT2b = MOTS 'LX' 'TINF' 'T   ' 'TSUP';
  230. D2b    = CHAN 'COMP' LMOT2a LMOT2b D2;
  231. T2c    = EXCO (MOTS 'TINF' 'T   ' 'TSUP') T2 ;
  232. D2c    = EXCO (MOTS 'TINF' 'T   ' 'TSUP') D2b;
  233. RESU2  = RESU (ABS (T2c - D2c));
  234. MAX2   = MAXI RESU2;
  235.  
  236.  
  237. ************************************************************************
  238. *                            VERIFICATION                              *
  239. ************************************************************************
  240. MESS 'MAXI DE LA RESULTANTE MAX0:' MAX0;
  241. PRECIS0 = PRECIS * (MAXI T0);
  242. SI (MAX0 >EG PRECIS0);
  243.   ERRE 5;
  244. FINS;
  245.  
  246. MESS 'MAXI DE LA RESULTANTE MAX1:' MAX1;
  247. PRECIS1 = PRECIS * (MAXI T1);
  248. SI (MAX1 >EG PRECIS1);
  249.   ERRE 5;
  250. FINS;
  251.  
  252. MESS 'MAXI DE LA RESULTANTE MAX2:' MAX2;
  253. PRECIS2 = PRECIS * (MAXI T2);
  254. SI (MAX2 >EG PRECIS2);
  255.   ERRE 5;
  256. FINS;
  257.  
  258.  
  259. FIN;
  260.  
  261.  
  262.  
  263.  

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