Télécharger multilayer.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : multilayer.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *******************************************************************
  5. * Cas test de non-régression CAST3M
  6. * Propriétés physiques de murs multicouches
  7. * auteur G.bonic 19/04/05
  8. *******************************************************************
  9.  
  10. * Procédure d'affichage
  11. 'DEBP' BOUYOU ;
  12. 'ARGU' RXT * 'TABLE' ;
  13. COMPT = PERSO.'COMPT';
  14. FREQ = PERSO.'FREQ' ;
  15. 'SI' (RXT.'TIC'.'NUPADT' 'EGA' ('ENTI'(COMPT * FREQ)));
  16. PERSO.'COMPT' = PERSO.'COMPT' + 1 ;
  17. VN = 'VECT' (RXT.'TIC'.'UN') 0.1 'UX' 'UY' 'BLANC';
  18. * TRAC VN RXT.'GEO'.'Mvtf' RXT.'TIC'.'TF' (CONT RXT.'GEO'.'Mvtf')
  19. * 'NCLK' ('TITR' RXT.'TIC'.'Tps') ;
  20. 'FINSI';
  21. FINPROC;
  22. ***********
  23. 'OPTI' 'DIME' 2;
  24. 'OPTI' 'ELEM' 'QUA8' ;
  25.  
  26. EPSI = 0.0001;
  27. P1 = -5. 0. ;
  28. P2 = 5. 0. ;
  29. P3 = 5.1 0. ;
  30. P4 = -5.1 0. ;
  31. Df = P1 'DROI' P2 50 ;
  32. Dwd = P2 'DROI' P3 1 COUL BLEU;
  33. Dwg = P4 'DROI' P1 1 COUL ROUGE;
  34. *trac (DF ET Dwd ET Dwg) 'NCLK';
  35. BASE = Dwg 'ET' (Df 'COUL' 'BLAN') 'ET' Dwd ;
  36. PLAF = BASE 'PLUS' (0. 10.) ;
  37. V1 = 'TRAN' 50 BASE (0. 10.) ;
  38. 'ELIM' EPSI (PLAF 'ET' V1);
  39. *trac V1 'NCLK';
  40. Wbas = 'TRAN' 1 BASE (0. -0.1) ;
  41. Whaut = 'TRAN' 1 PLAF (0. 0.1) ;
  42. Wcor = Wbas 'ET' Whaut 'COUL' 'VERT' ;
  43. Wd = V1 'ELEM' 'BLEU' ;
  44. Wg = V1 'ELEM' 'ROUG' ;
  45. Vflu = DIFF V1 Wd ;
  46. Vflu = DIFF Vflu Wg ;
  47. CAPDRO = Vflu 'POIN' 'PROC' (P1 'PLUS' ((49. * 10. / 50.) 5.));
  48. CAPGAU = Vflu 'POIN' 'PROC' (P1 'PLUS' ((1. * 10. / 50.) 5.));
  49. ELIM (Vflu 'ET' Wcor 'ET' WD 'ET' Wg 'ET' CAPDRO 'ET' CAPGAU) EPSI ;
  50. *trac (Vflu ET Wcor ET WD ET Wg) ;
  51. 'LIST' ('NBNO'(Vflu ET Wcor ET WD ET Wg)) ;
  52. Wex = CONT (Vflu ET Wcor ET WD ET Wg);
  53. *trac Wex 'NCLK';
  54. ****
  55. RXT = 'TABLE' ;
  56. RXT.'vtf' = Vflu ;
  57. *RXT.'breche' = AB ;
  58. *RXT.'axe' = AXE ;
  59. *RXT.'diru1' = (0 1) ;
  60. RXT.'pi' = ( Vflu 'POIN' 'PROC' (0.0 2.0)) ;
  61. RXT.'IMPR' = 2 ;
  62.  
  63. RXT.'THERMP' = VRAI ;
  64. RXT.'vtp' = (Wcor 'ET' Wg 'ET' Wd) ;
  65. RXT.'LAMBDA' = 3000.0 ;
  66. RXT.'ROCP' = (78.0 '*' 500.0) ;
  67. *RXT.'ROCP' = (78.0 '*' 500.0) ;
  68. RXT.'Tp0' = 500.0 ;
  69. RXT.'ECHAN' = 4000.0 ;
  70.  
  71. *RXT.'Lv' = 2.2E6 ;
  72.  
  73. RXT.'ECHEXT' = VRAI ;
  74. RXT.'parext' = Wex ;
  75. RXT.'HEXT' = 5000.0 ;
  76. RXT.'TPEXT' = 0.0 ;
  77.  
  78.  
  79. *RXT.'VAPEUR' = VRAI ;
  80. RXT.'AIR' = VRAI ;
  81. *RXT.'HE' = VRAI ;
  82. RXT.'TF0' = 500.1 ;
  83. RXT.'PT0' = 1.007E5 ;
  84. *RXT.'Yvap0' = 0.001 ;
  85. *RXT.'Yhe0' = 0.0 ;
  86.  
  87. *RXT.'MODTURB' = 'NUTURB' ;
  88. *RXT.'NUT' = 5.E-3 ;
  89. RXT.'MODTURB' = 'LMEL' ;
  90. *RXT.'LMEL' = 0.2 ;
  91. RXT.'LMEL' = 0.08 ;
  92.  
  93. RXT.'epsi' = 1.E-5 ;
  94. *RXT.'DISCR' = 'LINE' ;
  95. *RXT.'KPRE' = 'MSOMMET' ;
  96. RXT.'DISCR' = 'MACRO' ;
  97. RXT.'KPRE' = 'CENTREP1' ;
  98. RXT.'GRAPH' = FAUX ;
  99. RXT.'FRPREC' = 2 ;
  100.  
  101. PERSO=TABLE;
  102. PERSO.'COMPT' = 2 ;
  103. PERSO.'FREQ' = 1 ;
  104. rxt.'PERSO' = VRAI ;
  105. rxt.'PRCPERSO'= 'BOUYOU' ;
  106. rxt.'TABPERSO'= PERSO ;
  107.  
  108. RXT.'DT0' = 0.05 ;
  109. * Iteration fictive pour mettre en place les chpo à partir des flottants
  110. * initiaux
  111. EXECRXT 0 RXT ;
  112. * Mise en évidences des zones de la paroi où il faut imposer les nouveaux
  113. * champs
  114. cvtp = 'DOMA' RXT.'GEO'.'$vtp' 'CENTRE' ;
  115. cWg = 'CONT' Wg;
  116. cWd = 'CONT' Wd;
  117. cWh = 'CONT' Whaut;
  118. Mg = incl cvtp cWg 'BARY' ;
  119. Md = incl cvtp cWd 'BARY' ;
  120. Mh = incl cvtp cWh 'BARY' ;
  121. * Création des nouveaux CHPO
  122. l_g = 'MANU' 'CHPO' Mg 1 'SCAL' 15.0 'NATURE' 'DISCRET' ;
  123. l_d = 'MANU' 'CHPO' Md 1 'SCAL' 300.0 'NATURE' 'DISCRET' ;
  124. l_h = 'MANU' 'CHPO' Mh 1 'SCAL' 300.0 'NATURE' 'DISCRET' ;
  125. rxt.'TIC'.'LAMBDA' = 'KCHT' rxt.'GEO'.'$vtp' 'SCAL' 'CENTRE'
  126. 15.0 l_d l_g l_h;
  127. r_g = 'MANU' 'CHPO' Mg 1 'SCAL' (7800.0 '*' 500.0) 'NATURE' 'DISCRET' ;
  128. r_d = 'MANU' 'CHPO' Md 1 'SCAL' (78.0 '*' 500.0) 'NATURE' 'DISCRET' ;
  129. r_h = 'MANU' 'CHPO' Mh 1 'SCAL' (78.0 '*' 500.0) 'NATURE' 'DISCRET' ;
  130. rxt.'TIC'.'ROCP' = 'KCHT' rxt.'GEO'.'$vtp' 'SCAL' 'CENTRE'
  131. (7800.0 '*' 500.0) r_d r_g r_h;
  132.  
  133. EXECRXT 30 RXT;
  134. TOLER = 5.E-4 ;
  135. VALGr = 4.7680E+02 ;
  136. VALDR = 4.7679E+02 ;
  137. *VALD = 'REDU' rxt.tic.'TF' CAPDRO ;
  138. *VALG = 'REDU' rxt.tic.'TF' CAPGAU ;
  139.  
  140. VALD = 'EXTR' rxt.tic.'TF' 'SCAL' CAPDRO;
  141. VALG = 'EXTR' rxt.tic.'TF' 'SCAL' CAPGAU;
  142. 'MESS' VALGr VALDR;
  143. 'MESS' VALG VALD ;
  144. ERRE1 = 'ABS' (VALGr '-' VALG) '/' VALGr ;
  145. ERRE2 = 'ABS' (VALDr '-' VALD) '/' VALDr ;
  146. 'SI' (('MAXI' ('PROG' ERRE1 ERRE2)) '>EG' TOLER);
  147. 'MESSAGE' 'Il y a eu des erreurs' ;
  148. 'ERREUR' 5 ;
  149. 'FINSI';
  150. FIN;
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  

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