Télécharger fsi7.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : fsi7.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *************************************************************
  5. *
  6. * Cas test FSI7
  7. *
  8. * Calcul d'une masse ajoutee en mode de Fourier
  9. * (lame fluide)
  10. *
  11. *
  12. *
  13. * D. Combescure Décembre 2006
  14. *
  15. *************************************************************
  16. *
  17. *FLAGVISU = VRAI;
  18. *FLAGLONG = VRAI;
  19. *
  20. FLAGVISU = FAUX;
  21. FLAGLONG = FAUX;
  22. *
  23. opti dime 2 elem qua4 mode four 1;
  24. *
  25. *
  26. HH1 = 9.D-3;
  27. R1 = 7.85D-3;
  28. h1 = 0.3D-3;
  29. *
  30. HH2 = 20.D-3;
  31. R2 = 10.D-3;
  32. h2 = 0.6D-3;
  33. *
  34. SI FLAGLONG;
  35. nz12 = 10;
  36. nz1 = 10;
  37. nx1 = 5;
  38. nx2 = 5;
  39. SINON;
  40. nz12 = 1;
  41. nz1 = 1;
  42. nx1 = 1;
  43. nx2 = 1;
  44. FINSI;
  45. *
  46. p1 = 0. 0.;
  47. p2 = (R1 - h1) 0.;
  48. p3 = (R1 + h1) 0.;
  49. surbas1 = droi nx1 p2 p3;
  50. p4 = (R2 - h2) 0.;
  51. p5 = (R2 + h2) 0.;
  52. surbas2 = droi nx2 p4 p5;
  53. *
  54. p3s = p3 plus (0. 0.);
  55. p4s = p4 plus (0. 0.);
  56. sursol = droi 1 p3s p4s;
  57. *
  58. vz1= 0. HH1;
  59. p2h = p2 plus (0.5*vz1);
  60. p3h = p3 plus (0.5*vz1);
  61. p4h = p4 plus (0.5*vz1);
  62. p5h = p5 plus (0.5*vz1);
  63. p2b = p2 plus (-0.5*vz1);
  64. p3b = p3 plus (-0.5*vz1);
  65. p4b = p4 plus (-0.5*vz1);
  66. p5b = p5 plus (-0.5*vz1);
  67. *
  68. vz12= 0. (HH2 - HH1);
  69. p4h2 = p4h plus (0.5*vz12);
  70. p5h2 = p5h plus (0.5*vz12);
  71. p4b2 = p4b plus (-0.5*vz12);
  72. p5b2 = p5b plus (-0.5*vz12);
  73. *
  74. volint1 = (tran surbas1 nz1 (0.5*vz1)) et
  75. (tran surbas1 nz1 ((-0.5)*vz1)) ;
  76. volint2 = ((tran surbas2 nz1 (0.5*vz1)) et
  77. (tran (surbas2 plus (0.5*vz1)) nz12 (0.5*vz12)))
  78. et ((tran surbas2 nz1 ((-0.5)*vz1)) et
  79. (tran (surbas2 plus ((-0.5)*vz1)) nz12 ((-0.5)*vz12)));
  80. volsol = ((tran sursol nz1 (0.5*vz1)) et
  81. (tran (sursol plus (0.5*vz1)) nz12 (0.5*vz12)) et
  82. (tran sursol nz1 ((-0.5)*vz1)) et
  83. (tran (sursol plus ((-0.5)*vz1)) nz12 ((-0.5)*vz12)))
  84. coul rouge;
  85. *
  86. SI FLAGVISU;
  87. trac (volsol et volint1 et volint2);
  88. FINSI;
  89. *
  90. surint = (d nz1 p2 p2h)
  91. et (d nz1 p2 p2b);
  92. surint2 = (d nz1 p3 p3h)
  93. et (d nz1 p3 p3b);
  94. surext = (d nz1 p4 p4h) et (d nz12 p4h p4h2)
  95. et (d nz1 p4 p4b)
  96. et (d nz12 p4b p4b2);
  97. surext2 = (d nz1 p5 p5h) et (d nz12 p5h p5h2)
  98. et (d nz1 p5 p5b)
  99. et (d nz12 p5b p5b2);
  100. surver1 = (surbas1 plus (0.5*vz1))
  101. et (surbas1 plus (-0.5*vz1));
  102. surver2 = (surbas2 plus (0.5*(vz1 plus vz12)))
  103. et (surbas2 plus (-0.5*(vz1 plus vz12)));
  104. *
  105. elim 0.0001
  106. (surint et surint2 et surext et surext2
  107. et volint1 et volint2 et surver1 et surver2);
  108. surextb = surext plus (0. 0.);
  109. surint2b = surint2 plus (0. 0.);
  110. elim 0.0001
  111. (surint2b et surextb et volsol);
  112. *
  113. SI FLAGVISU;
  114. trac (surint et surext et surint2 et surext2
  115. et volint1 et volint2 et volsol);
  116. FINSI;
  117. *
  118. *
  119. *
  120. * Modèle de fluide
  121. *
  122. ro1 = 800.;
  123. cf1 = 343.;
  124. *
  125. mod1 = MODE (volint1) liquide lqu4;
  126. mat1 = MATE mod1 rho ro1 rorf ro1
  127. cson cf1 cref cf1 lcar 1.D-3 g 0.;
  128. mod2 = MODE (volint2) liquide lqu4;
  129. mat2 = MATE mod2 rho ro1 rorf ro1
  130. cson cf1 cref cf1 lcar 1.D-3 g 0.;
  131. rig1 = rigi (mod1 et mod2) (mat1 et mat2);
  132. mas1 = mass (mod1 et mod2) (mat1 et mat2);
  133. *
  134. blfl = bloq UR UT UZ
  135. (surint et surext2 );
  136. blflI = bloq IUR IUT IUZ
  137. (surint et surext2 );
  138. blfl2 = bloq UZ
  139. (surver1 et surver2);
  140. blfl2I = bloq IUZ
  141. (surver1 et surver2);
  142. *
  143. *
  144. opti elem rac2;
  145. mrac2 = racc 0.0001 (surext)
  146. (surextb);
  147. mrac1 = racc 0.0001 (surint2)
  148. (surint2b);
  149. *
  150. modrac1 = MODE mrac1 mecanique liquide rac2;
  151. matrac1 = cara modrac1 liqu volint1;
  152. *
  153. modrac2 = MODE mrac2 mecanique liquide rac2;
  154. matrac2 = cara modrac2 liqu volint2;
  155. *
  156. modliq = MODE (volint1 et volint2) liquide lqu4;
  157. matliq = MATE (modliq et modrac1 et modrac2) rho ro1 rorf ro1
  158. cson cf1 cref cf1 lcar 1. g 0.;
  159. *
  160. bl1 = bloq UR UT UZ (surext2 et surint);
  161. bl2ur = bloq UR (surextb et surint2b);
  162. bl2ut = bloq UT (surextb et surint2b);
  163. bl2uz = bloq UZ (surextb et surint2b);
  164. bl1I = bloq IUR IUT IUZ (surext2 et surint);
  165. bl2Iur = bloq IUR (surextb et surint2b);
  166. bl2Iut = bloq IUT (surextb et surint2b);
  167. bl2Iuz = bloq IUZ (surextb et surint2b);
  168. dep1ur = depi bl2ur 1.;
  169. dep1ut = depi bl2ut -1.;
  170. Ktot = rigi (modrac1 et modrac2 et modliq)
  171. (matrac1 et matrac2 et matliq) ;
  172. Mtot = mass (modrac1 et modrac2 et modliq)
  173. (matrac1 et matrac2 et matliq);
  174. *
  175. Mimpe = 'IMPE' Mtot 1. 'MASSE';
  176. Kimpe = 'IMPE' Ktot 'RAIDEUR';
  177. *
  178. SI FLAGLONG;
  179. prfreq = prog 10 pas 10. 1000.;
  180. SINON;
  181. prfreq = prog 100 pas 100. 1000.;
  182. FINSI;
  183. *
  184. prFR = prog;
  185. prFT = prog;
  186. repeter lab1 (dime prfreq);
  187. fr1 = extr prfreq &lab1;
  188. ImpTot = Kimpe et (((2.*pi*fr1)**2)*Mimpe)
  189. et bl1 et bl2ur et bl2ut et bl2uz et
  190. bl1I et bl2Iur et bl2Iut et bl2Iuz et blfl et blflI
  191. et blfl2 et blfl2I;
  192. deptot = reso ImpTot (dep1ur et dep1ut);
  193. reatot = reac deptot (bl2ur et bl2ut);
  194. prFR = prfr et (prog (maxi (exco (resu reatot) FR)));
  195. prFT = prft et (prog (maxi (exco (resu reatot) FT)));
  196. fin lab1;
  197. *
  198. evfr = evol manu prfreq prfr;
  199. evft = evol manu prfreq prft;
  200. *
  201. mm = 57.6d-3;
  202. evfrma = evol rouge manu prfreq ((-1.)*((2.*pi*prfreq)**2)*mm);
  203. *
  204. TEST = ABS ((SOMM ((extr evfrma ordo) - (extr evfr ordo)))/
  205. (SOMM (extr evfrma ordo)));
  206. *
  207. SI FLAGVISU;
  208. trac (volint1 et volint2) (exco deptot 'P');
  209. trac (volint1 et volint2) (vecteur reatot 'FR' 'FZ');
  210. dess (evfr et evfrma);
  211. FINSI;
  212. *
  213. SI (TEST >EG 1.d-2);
  214. ERRE 4;
  215. FINSI;
  216. *
  217. tab1 = table;
  218. tab1.'CHAMPS_POINTS_SYMETRIQUES' = TABLE;
  219. tab1.'CHAMPS_POINTS_SYMETRIQUES'. 1 = exco deptot 'P' 'SCAL';
  220. tab1.'EFFORTS' = TABLE;
  221. tab1.'EFFORTS'. 1 = reatot;
  222. mesh3d1 = four2tri (modrac1 et modrac2 et modliq)
  223. tab1 (prog 0. pas 10. 270.) 1;
  224. SI FLAGVISU;
  225. trac mesh3D1 (tab1.'CHAMPS_POINTS_SYME_3D'. 1 );
  226. trac mesh3D1 (vecteur (tab1.'EFFORTS_3D'. 1 ) FX FY FZ);
  227. FINSI;
  228. *
  229. FIN;
  230. *
  231.  
  232.  
  233.  
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  
  244.  
  245.  
  246.  

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