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. modtot = moddif0 'ET' moddif1 'ET' moddif2 ;
  96. lincp = 'EXTR' modtot 'DEPL' ;
  97. lincd = 'EXTR' modtot 'FORC' ;
  98. opti inco lincp lincd ;
  99.  
  100. ************************************************************************
  101. *                               MATERIAU                               *
  102. ************************************************************************
  103. MATTHE0 = MATE MODTHE0 'K'  2.D0 'RHO'  1.D0 'C' 1.D0             ;
  104. MATTHE1 = MATE MODTHE1 'K'  2.D0 'RHO'  1.D0 'C' 1.D0 'EPAI' 1.D-5;
  105. MATTHE2 = MATE MODTHE2 'K'  2.D0 'RHO'  1.D0 'C' 1.D0 'EPAI' 1.D-5;
  106. MATDIF0 = MATE MODDIF0 'KD' 2.D0 'CDIF' 1.D0                      ;
  107. MATDIF1 = MATE MODDIF1 'KD' 2.D0 'CDIF' 1.D0          'EPAI' 1.D-5;
  108. MATDIF2 = MATE MODDIF2 'KD' 2.D0 'CDIF' 1.D0          'EPAI' 1.D-5;
  109.  
  110. SI BLIST;
  111.   LIST MATTHE0;
  112.   LIST MATDIF0;
  113.   LIST MATTHE1;
  114.   LIST MATDIF1;
  115.   LIST MATTHE2;
  116.   LIST MATDIF2;
  117. FINS;
  118.  
  119. ************************************************************************
  120. *                             CONDUCTIVITE                             *
  121. ************************************************************************
  122. CONDTH0 = COND MODTHE0 MATTHE0;
  123. CONDDI0 = COND MODDIF0 MATDIF0;
  124. CONDTH1 = COND MODTHE1 MATTHE1;
  125. CONDDI1 = COND MODDIF1 MATDIF1;
  126. CONDTH2 = COND MODTHE2 MATTHE2;
  127. CONDDI2 = COND MODDIF2 MATDIF2;
  128.  
  129. SI BLIST;
  130.   LIST CONDTH0;
  131.   LIST CONDDI0;
  132.   LIST CONDTH1;
  133.   LIST CONDDI1;
  134.   LIST CONDTH2;
  135.   LIST CONDDI2;
  136. FINS;
  137.  
  138.  
  139. ************************************************************************
  140. *                             BLOCAGES                                 *
  141. ************************************************************************
  142. BLTH1 = BLOQ 'T' P1;
  143. DDTH1 = DEPI BLTH1 0.D0;
  144. BLTH2 = BLOQ 'T' P5;
  145. DDTH2 = DEPI BLTH2 1.D0;
  146. BLTH3 = BLOQ 'T' P3;
  147. DDTH3 = DEPI BLTH3 0.D0;
  148. BLTH4 = BLOQ 'T' P4;
  149. DDTH4 = DEPI BLTH4 1.D0;
  150.  
  151. BLDI1 = BLOQ MP P1;
  152. DDDI1 = DEPI BLDI1 0.D0;
  153. BLDI2 = BLOQ MP P5;
  154. DDDI2 = DEPI BLDI2 1.D0;
  155. BLDI3 = BLOQ MP P3;
  156. DDDI3 = DEPI BLDI3 0.D0;
  157. BLDI4 = BLOQ MP P4;
  158. DDDI4 = DEPI BLDI4 1.D0;
  159.  
  160. ************************************************************************
  161. *                             RESOLUTION                               *
  162. ************************************************************************
  163. RIGTH0 = CONDTH0 ET BLTH1 ET BLTH2 ET BLTH3 ET BLTH4;
  164. SECTH0 = DDTH1 ET DDTH2 ET DDTH3 ET DDTH4;
  165. T0 = RESO RIGTH0 SECTH0;
  166.  
  167. RIGTH1 = CONDTH1 ET BLTH1 ET BLTH2 ET BLTH3 ET BLTH4;
  168. SECTH1 = DDTH1 ET DDTH2 ET DDTH3 ET DDTH4;
  169. T1 = RESO RIGTH1 SECTH1;
  170.  
  171. RIGTH2 = CONDTH2 ET BLTH1 ET BLTH2 ET BLTH3 ET BLTH4;
  172. SECTH2 = DDTH1 ET DDTH2 ET DDTH3 ET DDTH4;
  173. T2 = RESO RIGTH2 SECTH2;
  174.  
  175. RIGDI0 = CONDDI0 ET BLDI1 ET BLDI2 ET BLDI3 ET BLDI4;
  176. SECDI0 = DDDI1 ET DDDI2 ET DDDI3 ET DDDI4;
  177. D0 = RESO RIGDI0 SECDI0;
  178.  
  179. RIGDI1 = CONDDI1 ET BLDI1 ET BLDI2 ET BLDI3 ET BLDI4;
  180. SECDI1 = DDDI1 ET DDDI2 ET DDDI3 ET DDDI4;
  181. D1 = RESO RIGDI1 SECDI1;
  182.  
  183. RIGDI2 = CONDDI2 ET BLDI1 ET BLDI2 ET BLDI3 ET BLDI4;
  184. SECDI2 = DDDI1 ET DDDI2 ET DDDI3 ET DDDI4;
  185. D2 = RESO RIGDI2 SECDI2;
  186.  
  187. SI BTRAC;
  188.   TRAC (T0 ET D0)  V1         ;
  189.   TRAC (T1 ET D1) (S1  ET S3 );
  190.   TRAC (T2 ET D2) (S11 ET S31);
  191. FINS;
  192.  
  193.  
  194. ************************************************************************
  195. *                               CAPACITE                               *
  196. ************************************************************************
  197. CAPATH0 = CAPA MODTHE0 MATTHE0;
  198. CAPATH1 = CAPA MODTHE1 MATTHE1;
  199. CAPATH2 = CAPA MODTHE2 MATTHE2;
  200. CAPADI0 = CAPA MODDIF0 MATDIF0;
  201. CAPADI1 = CAPA MODDIF1 MATDIF1;
  202. CAPADI2 = CAPA MODDIF2 MATDIF2;
  203.  
  204. SI BLIST;
  205.   LIST CAPATH0;
  206.   LIST CAPADI0;
  207.   LIST CAPATH1;
  208.   LIST CAPADI1;
  209.   LIST CAPATH2;
  210.   LIST CAPADI2;
  211. FINS;
  212.  
  213.  
  214. ************************************************************************
  215. *                            COMPARAISONS                              *
  216. ************************************************************************
  217. LMOT0a = MOTS 'LX' 'B';
  218. LMOT0b = MOTS 'LX' 'T';
  219. D0b    = CHAN 'COMP' LMOT0a LMOT0b D0;
  220. T0c    = EXCO 'T' T0 ;
  221. D0c    = EXCO 'T' D0b;
  222. RESU0  = RESU (ABS (T0c - D0c));
  223. MAX0   = MAXI RESU0;
  224.  
  225. LMOT1a = MOTS 'LX' 'BIN ' 'B   ' 'BSU ';
  226. LMOT1b = MOTS 'LX' 'TINF' 'T   ' 'TSUP';
  227. D1b    = CHAN 'COMP' LMOT1a LMOT1b D1;
  228. T1c    = EXCO (MOTS 'TINF' 'T   ' 'TSUP') T1 ;
  229. D1c    = EXCO (MOTS 'TINF' 'T   ' 'TSUP') D1b;
  230. RESU1  = RESU (ABS (T1c - D1c));
  231. MAX1   = MAXI RESU1;
  232.  
  233. LMOT2a = MOTS 'LX' 'BIN ' 'B   ' 'BSU ';
  234. LMOT2b = MOTS 'LX' 'TINF' 'T   ' 'TSUP';
  235. D2b    = CHAN 'COMP' LMOT2a LMOT2b D2;
  236. T2c    = EXCO (MOTS 'TINF' 'T   ' 'TSUP') T2 ;
  237. D2c    = EXCO (MOTS 'TINF' 'T   ' 'TSUP') D2b;
  238. RESU2  = RESU (ABS (T2c - D2c));
  239. MAX2   = MAXI RESU2;
  240.  
  241.  
  242. ************************************************************************
  243. *                            VERIFICATION                              *
  244. ************************************************************************
  245. MESS 'MAXI DE LA RESULTANTE MAX0:' MAX0;
  246. PRECIS0 = PRECIS * (MAXI T0);
  247. SI (MAX0 >EG PRECIS0);
  248.   ERRE 5;
  249. FINS;
  250.  
  251. MESS 'MAXI DE LA RESULTANTE MAX1:' MAX1;
  252. PRECIS1 = PRECIS * (MAXI T1);
  253. SI (MAX1 >EG PRECIS1);
  254.   ERRE 5;
  255. FINS;
  256.  
  257. MESS 'MAXI DE LA RESULTANTE MAX2:' MAX2;
  258. PRECIS2 = PRECIS * (MAXI T2);
  259. SI (MAX2 >EG PRECIS2);
  260.   ERRE 5;
  261. FINS;
  262.  
  263.  
  264. FIN;
  265.  
  266.  
  267.  
  268.  
  269.  
  270.  

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