Télécharger oscicyl2.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************************
  2. * NOM : OSCICYL2
  3. * DESCRIPTION : Cylindre oscillant dans une cavité fermée
  4. *
  5. * Deux paramètres : le nombre de Reynolds (convection/viscosité)
  6. * le nombre de Keulegan-Carpenter (convection/masse,trainee)
  7. *
  8. *
  9. * Bibliographie : thèse de Marion Duclercq
  10. * avec Daniel BROC
  11. *
  12. * LANGAGE : GIBIANE-CAST3M
  13. * AUTEUR : Stephane GOUNAND et Dihya MEZI (CEA/DES/ISAS/DM2S/SEMT/LTA)
  14. * mel : stephane.gounand@cea.fr
  15. **********************************************************************
  16. * VERSION : v1, 28/01/2022, version initiale
  17. * HISTORIQUE : v1, 28/01/2022, création
  18. * HISTORIQUE :
  19. * HISTORIQUE :
  20. ************************************************************************
  21. *
  22. *
  23. 'OPTION' 'DIME' 2 'ELEM' 'QUA8' ;
  24. 'OPTI' 'ISOV' suli ;
  25.  
  26. graph = faux ;
  27. *
  28. ****Definitions des parametres des longueurs
  29.  
  30. Long = 2. ;
  31. Haut = 1. ;
  32. Rayon = 20.E-2;
  33.  
  34. * Parametres pour le maillage
  35.  
  36. nx=20;
  37. ny=10;
  38. nc=5;
  39. *nx = nx * 2 ;
  40. *ny = ny * 2 ;
  41. *nc = nc * 2 ;
  42.  
  43. ********** Creation du contour exterieur
  44. pA= 0. 0. ; pB=Long 0. ; pC= Long Haut ; pD= 0. Haut ;
  45. bas = 'DROI' nx pA pB;
  46. dro= 'DROI' ny pB pC;
  47. hau= 'DROI' nx pC pD;
  48. gau = 'DROI' ny pD pA;
  49. **** Creation du cercle
  50. pcen= (Long * 0.5) (Haut * 0.5) ;
  51. pE = ((0.5 * Long) '-' Rayon) (0.5 * Haut);
  52. pF = ((0.5 * Long) ) ((0.5 * Haut) '+'Rayon );
  53. pG= ((0.5 * Long)'+' Rayon ) ((0.5 * Haut) );
  54. pH=((0.5 * Long) ) ((0.5 * Haut) '-' Rayon );
  55. CE1 = 'CERCLE' nc pE pcen pF;
  56. CE2= 'CERCLE' nc CE1 pcen pG;
  57. CE3= 'CERCLE' nc CE2 pcen pH;
  58. CE= 'CERCLE' nc CE3 pcen pE;
  59. * Contour du maillage exterieur
  60. cmt1= bas 'ET' dro 'ET' hau 'ET' gau ;
  61. *** Contour total avec cercle
  62. cmt= bas 'ET' dro 'ET' hau 'ET' gau 'DIFF' CE;
  63. mt = 'SURFACE' cmt;
  64. 'ELIMINATION' mt 1.E-9;
  65. * Changer le maillage en un maillage quadratique fluide
  66. _mt='CHAN' mt 'QUAF';
  67. *'TRAC' _mt nclk ;
  68.  
  69.  
  70. **** Conditions de blocage
  71. *** Vitesse
  72. mblox = 'BLOQUE' 'UX' cmt1 ;
  73. mbloy = 'BLOQUE' 'UY' cmt1 ;
  74. mblocx = 'BLOQUE' 'UX' CE ;
  75. mblocy = 'BLOQUE' 'UY' CE ;
  76.  
  77. **** Physique
  78. *** nombre de Reynolds de 100
  79. *** nombre de Keulegan-Carpenter de 2
  80. Re = 100. ;
  81. nre= 1 '/' Re ;
  82. KC = 2. ;
  83.  
  84. **** Numérique
  85. *** Pas de temps (cf. thèse Marion Duclercq)
  86. 'SI' ('>' KC 2.) ;
  87. dt = 0.05 '/' KC ;
  88. 'SINO' ;
  89. dt = 0.025 ;
  90. 'FINS' ;
  91. *dt = dt * 4. ;
  92.  
  93.  
  94. mns1 = mode _mt 'NAVIER_STOKES' 'NLIN' 'QUAF' 'LINE' 'CHPO' 'GAP7' ;
  95. cns1 = mate mns1 REYN Re FREQ (dt**(-1));
  96.  
  97. *** Pression
  98. mapcen = extr mns1 'CENT' ;
  99. pp2 = 'POIN' mapcen 'PROC' (0.5 0.5) ;
  100. mblop = 'BLOQ' 'LX1' pp2 ;
  101.  
  102.  
  103. mtbid = changer _mt line ;
  104. mbid = mode mtbid mecanique elastique ;
  105. cbid = mate mbid youn 1.e10 nu 0.3 rho 2.;
  106.  
  107. **** Matrices fixes du probleme
  108. *** Rigidité (objet RIGIDITE)
  109. * mrigi = 'NLIN' gdisc _mt A B mot1 mot2 ;
  110. nrigi = rigi (mns1 et mbid) (cns1 et cbid) ;
  111. rigbid = rigi (mbid) (cbid) ;
  112. *** Contrainte div U
  113. ndiv = rigi 'DIVE' (mns1 et mbid) (cns1 et cbid) ;
  114. *** Masse pondérée par le pas de temps
  115. nmass = mass (mns1 et mbid) (cns1 et cbid) ;
  116.  
  117.  
  118.  
  119. **** Conditions initiales
  120. U0 = 'MANUEL' 'CHPO' _mt 2 'UX' 0. 'UY' 0. ;
  121. UNM = U0 ;
  122. UN = U0 ;
  123. solu= UN;
  124. EPS= 1.0D-9;
  125. tps= 0. ;
  126.  
  127. *
  128. *
  129. abs_t = prog 0. pas Dt (105*Dt) ;
  130. ord_t = 'SIN' (360. '*' abs_t) ;
  131. ev_t = evol vert manu ' t' abs_t ord_t ;
  132. abcs_CE = ((coor 1 CE) - (coor 1 pcen))/Rayon ; abcs_CE = abs abcs_CE ;
  133. fbl_CE = 'DEPI' mblocx ('NOMC' 'UX' abcs_CE) ;
  134. char_t = 'CHAR' 'DIMP' fbl_CE ev_t ;
  135.  
  136. *
  137.  
  138. t_fltub = table 'PASAPAS' ;
  139. t_fltub . modele = mns1 ;
  140. t_fltub . caracteristiques = cns1 ;
  141. t_fltub . blocages_mecaniques = mblox 'ET'
  142. mbloy 'ET' mblop 'ET' mblocx 'ET' mblocy ;
  143. t_fltub . chargement = char_t ;
  144. t_fltub . temps_calcules = abs_t ;
  145. *t_fltub . temps_calcules = prog 0. Dt ;
  146.  
  147. pasapas t_fltub ;
  148.  
  149.  
  150. nomvit = 'EXTR' mns1 'DEPL' ;
  151.  
  152. si (graph) ;
  153. repeter BTRA ((dime abs_t) - 1) ;
  154. solu = t_fltub . vitesses_fluide . &btra ;
  155. unpip = 'EXCO' solu nomvit ;
  156. *** synthèse
  157. vun = 'VECTEUR' unpip (* (4. '/' nx) 1.2) nomvit 'JAUN' ;
  158. pnp = t_fltub . pression_fluide . &btra ;
  159. 'TRAC' pnp vun _mt (contour mt) 16 'TITR' ' PRESSION AUX NOEUDS' ;
  160.  
  161. fin btra ;
  162.  
  163. 'TRAC' pnp vun _mt 16 'TITR' ' PRESSION AUX NOEUDS' ;
  164. finsi ;
  165. *opti donn 5 ;
  166. *
  167. * End of dgibi file oscillations
  168. *
  169. 'FIN' ;
  170.  
  171.  
  172.  
  173.  

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