Télécharger flslic4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : flslic4.dgibi
  2. ******************************************************************
  3. * *
  4. * test de validation de l' élément de raccord fluide - *
  5. * *
  6. * structure LIC4. *
  7. * *
  8. * calcul du mode n=3 m=1 d'un système de coques *
  9. * *
  10. * concentriques. La coque exterieure est supposée *
  11. * *
  12. * encastrée, la coque interieure est libre. *
  13. * *
  14. * valeurs de comparaison: calcul en mode fourier n=3 *
  15. * *
  16. * Les valeurs sont Fréquence = 47.491 *
  17. * Masse généralisée = 2563.1 *
  18. * *
  19. * On ne modélise qu'un quart de modèle compte tenu des *
  20. * *
  21. * symetries. *
  22. * *
  23. * Pour la coque externe, les conditions sont ux=uy=uz=0 *
  24. * rx=ry=rz=0 *
  25. * Pour la coque interne les conditions sont *
  26. * plan xz SYMETRIE *
  27. * plan yz ANTISYMETRIE *
  28. * uz=0 *
  29. * Pour le fluide les conditions sont *
  30. * plan xz dP/dn=0 *
  31. * plan yz P=0 *
  32. * plan x=0 et x=h dP/dn *
  33. * *
  34. ******************************************************************
  35. *
  36. opti elem cub8 ;
  37. opti echo 0 ;
  38. OPTI EPSI LINEAIRE;
  39. *
  40. ****************************************
  41. * definition des donnees parametriques *
  42. ****************************************
  43. *
  44. ri =1.00; comm rayon de la coque interne ;
  45. re =1.30; comm rayon de la coque externe ;
  46. h =1.00; comm hauteur du modele ;
  47. hc =0.04; comm epaisseur des coques ;
  48. lc0=5.3; comm longueur caracteristique du fluide ;
  49. y0 =2.1e11; comm module d YOUNG des coques ;
  50. roc=7900; comm masse volumique de l acier ;
  51. decouz= 5; comm nombre d elements suivant z ;
  52. decout= 10; comm nombre d elements suivant teta ;
  53. decouf= 5; comm nombre d elements sur le rayon ;
  54. vect= 0. 0. 0.0; comm vect translation pour raccord fluide
  55. structure ;
  56. *
  57. *************************
  58. * definition des points *
  59. *************************
  60. *
  61. pci1= ri 0.00 0.00 ;
  62. pci2= ri 0.00 h ;
  63. pce1= re 0.00 0.00 ;
  64. pce2= re 0.00 h ;
  65. ce1=0.00 0.00 0.00 ;
  66. ce2=0.00 0.00 h ;
  67. *
  68. *************************
  69. * definition des coques *
  70. *************************
  71. *
  72. lci= pci1 d decouz pci2 ;
  73. lce= pce1 d decouz pce2 ;
  74. sci= lci rota decout 90 ce1 ce2 ;
  75. sce= lce rota decout 90 ce1 ce2 ;
  76. *
  77. ************************
  78. * definition du fluide *
  79. ************************
  80. *
  81. lcf1 = lci plus vect ;
  82. lcf2 = lce moins vect ;
  83. sflu1= lcf1 rota decout 90 ce1 ce2 ;
  84. sflu2= lcf2 rota decout 90 ce1 ce2 ;
  85. vflu = volu sflu1 decouf sflu2 ;
  86. *
  87. ***************************************
  88. * definition des elements de raccords *
  89. ***************************************
  90. *
  91. opti elem lia4 ;
  92. raci= liaison 0.001 sflu1 sci ;
  93. race= liaison 0.001 sflu2 sce ;
  94. *
  95. *****************************************
  96. * definition des modeles elements finis *
  97. *****************************************
  98. *
  99. moci = modele sci 'MECANIQUE' 'ELASTIQUE' coq4 ;
  100. moce = modele sce 'MECANIQUE' 'ELASTIQUE' coq4 ;
  101. mocri= modele raci 'MECANIQUE' 'LIQUIDE' lic4 ;
  102. mocre= modele race 'MECANIQUE' 'LIQUIDE' lic4 ;
  103. mocf = modele vflu 'LIQUIDE' lcu8 ;
  104. *
  105. ****************************
  106. * definition des materiaux *
  107. ****************************
  108. *
  109. maci = mate moci youn y0 nu 0.3 rho roc epai hc ;
  110. mace = mate moce youn y0 nu 0.3 rho roc epai hc ;
  111. macf = mate mocf rho 1000 rorf 1000 cson 1435 cref 1435
  112. lcar lc0 g 9.81 ;
  113. macri= mate mocri rho 1000 rorf 1000 cson 1435 cref 1435
  114. lcar lc0 g 9.81 ;
  115. macre= mate mocre rho 1000 rorf 1000 cson 1435 cref 1435
  116. lcar lc0 g 9.81 ;
  117. cacri= carac mocri 'LIQUIDE' vflu ;
  118. cacre= carac mocre 'LIQUIDE' vflu ;
  119. macri= macri et cacri ;
  120. macre= macre et cacre ;
  121. macr = macri et macre ;
  122. *
  123. ************************
  124. * calcul des rigidites *
  125. ************************
  126. *
  127. rigci= rigi moci maci ;
  128. rigce= rigi moce mace ;
  129. rigcf= rigi mocf macf ;
  130. rigcr= rigi (mocri et mocre) macr;
  131. rigt = rigci et rigce et rigcf et rigcr ;
  132. *
  133. ********************************
  134. * calcul des matrices de masse *
  135. ********************************
  136. *
  137. masci= masse moci maci ;
  138. masce= masse moce mace ;
  139. mascf= masse mocf macf ;
  140. mascr= masse (mocri et mocre) macr ;
  141. mast = masci et masce et mascf et mascr ;
  142. *
  143. *****************************************
  144. * definition des conditions aux limites *
  145. *****************************************
  146. *
  147. pcr1= pci1 tour 90 ce1 ce2 ;
  148. pcr2= pci2 tour 90 ce1 ce2 ;
  149. pcr3= pce2 tour 90 ce1 ce2 ;
  150. pnx= sci point plan ce1 pci1 pci2 0.001 ;
  151. vpny= vflu point plan pcr1 pcr2 pcr3 0.001 ;
  152. pnz= sci point plan pcr1 ce1 pci1 0.001 ;
  153. *
  154. blo1= bloque depla rota sce ;
  155. blo2= symt 'DEPL' ce1 pci1 pci2 sci 0.001 ;
  156. blo3= anti 'DEPL' pcr1 pcr2 pcr3 sci 0.001 ;
  157. blo4= symt 'ROTA' ce1 pci1 pci2 sci 0.001 ;
  158. blo5= anti 'ROTA' pcr1 pcr2 pcr3 sci 0.001 ;
  159. blo6= bloque uz sci ;
  160. blop = bloque 'P' 'PI' vpny ;
  161. blot= blo1 et blo2 et blo3 et blo4 et blo5 et blo6 et blop ;
  162. *
  163. ********************************
  164. * calcul de la rigidite totale *
  165. ********************************
  166. *
  167. rigt= rigt et blot ;
  168. *
  169. ************************************
  170. * utilisation de vibre option proc *
  171. ************************************
  172. *
  173. tmx1= vibr proch (prog 40.) rigt mast ;
  174. *
  175. ****************************
  176. * extraction des resultats *
  177. ****************************
  178. *
  179. kfois=1 ;
  180. freq = tmx1.modes.kfois.frequence ;
  181. utot = tmx1.modes.kfois.deformee_modale ;
  182. mn = tmx1.modes.kfois.masse_generalisee ;
  183. *
  184. *************************************************
  185. * recalcul par modification de la normalisation *
  186. * maxi du deplacement ux =1 *
  187. *************************************************
  188. *
  189. uxm = maxi (exco ux utot) ;
  190. unox= utot/uxm ;
  191. udx1= manu 'CHPOINT' (sci et sce et raci et race) 6
  192. ux 1.0 uy 0.0 uz 0.0 rx 0.0 ry 0.0 rz 0.0 nature diffus ;
  193. udx2= manu 'CHPOINT' (raci et race et vflu) 2
  194. 'P' 0.0 'PI' 0.0 nature diffus ;
  195. udx = udx1 + udx2 ;
  196. qx = ytmx unox udx mast ;
  197. mnx = ytmx unox unox mast ;
  198. mnx = mnx ;
  199. *
  200. *****************************************
  201. * comparaison aux valeurs de références *
  202. *****************************************
  203. *
  204. iuer=0.05 ;
  205. tvalth= table ;
  206. tvalth.1=47.491 ;
  207. tvalth.2=2563.1 ;
  208. tvalca= table ;
  209. tvalca.1= freq ;
  210. tvalca.2= mnx ;
  211. ierc=table ;
  212. i=1 ;
  213. iercmax= 0.0 ;
  214. repeter bverif 2 ;
  215. ierc.i= ((tvalth.i)-(tvalca.i)) ;
  216. ierc.i= (ierc.i)/(tvalth.i) ;
  217. si ((ierc.i) < 0.0) ;
  218. xvali= ierc.i ;
  219. ierc.i=abs xvali ;
  220. finsi ;
  221. si (iercmax < (ierc.i)) ;
  222. iercmax= (ierc.i) ;
  223. finsi ;
  224. i=i+1 ;
  225. fin bverif ;
  226. si (iercmax < iuer) ;
  227. mess 'bonne fin du calcul' ;
  228. sinon ;
  229. mess 'ecart maximal vaut ' iercmax 'il y a un probleme' ;
  230. finsi ;
  231. fin ;
  232.  
  233.  
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  

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