Télécharger la plaquette Cast3M


fr en

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. *
  37. opti elem cub8                                                   ;
  38. opti echo 0                                                      ;
  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.  
  240.  
  241.  

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