Télécharger trac3d.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : trac3d.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. ******************************************************************
  6. * *
  7. * Mots-clés : Vibrations, calcul modal, coque, 2D Fourier, *
  8. * Reconstruction 3D *
  9. * *
  10. * TEST TRAC3D *
  11. * *
  12. * *
  13. * P3 *
  14. * | *
  15. * | *
  16. * | *
  17. * | *
  18. * | *
  19. * | *
  20. * | *
  21. * ______________| *
  22. * P1 P2 *
  23. * Soit un cylindre dont on calcule la fréquence propre du *
  24. * deuxieme mode de flexion en mode fourier 1 puis 3. On *
  25. * construit alors les champs de déplacement tridimensionnels *
  26. * ainsi que le maillage 3D. on superpose les deux calculs pour *
  27. * obtenir un champs de deplacement 3D et une déformée. *
  28. * *
  29. ******************************************************************
  30. *
  31. OPTION ECHO 1;
  32. GRAPH='N';
  33. SAUT PAGE ;
  34. OPTI NORM AUTO ;
  35. OPTI MODE FOUR NOHARM ;
  36. OPTI ELEM SEG2;
  37. *
  38. *================================================================*
  39. * GEOMETRY *
  40. *==============================================================*
  41. *
  42. P1 = 0. 0. ;
  43. P2 = 3.5 0.0 ;
  44. P3 = 3.5 7.0 ;
  45. *
  46. LI1 = P1 DROIT 3 P2 ;
  47. STEEL = LI1 ET (P2 DROIT 8 P3) ;
  48. *
  49. *==============================================================*
  50. * MODELE - MATERIAU - RIGIDITE - MASSE *
  51. *==============================================================*
  52. MODCOQ=MODE STEEL MECANIQUE COQ2 ;
  53. MATCOQ=MATE MODCOQ RHO 8.0E3 YOUN 2.E11 NU 0.3 ;
  54. CARCOQ=CARA MODCOQ EPAI 0.035 ;
  55. MATCOQ=MATCOQ ET CARCOQ;
  56. *==============================================================*
  57. * BOUNDARY CONDITION *
  58. *==============================================================*
  59. CL1 = BLOQ DEPL ROTA LI1 ;
  60. *==============================================================*
  61. *============= MODE FOURIER 1 ================================*
  62. *==============================================================*
  63. OPTI MODE FOUR 1 ;
  64. RIG1 = RIGI MODCOQ MATCOQ ;
  65. MAS1 = MASS MODCOQ MATCOQ ;
  66. *==============================================================*
  67. * CALCULATION OF THE FREQUENCIES *
  68. * AND *
  69. * EXTRACTION OF SOME RESULTS *
  70. *==============================================================*
  71. RESUL = VIBR proche (prog 150.) (CL1 ET RIG1 ) MAS1 IMPR ;
  72. FREQ1 = RESUL.MODES. 1 . FREQUENCE ;
  73. DEP1 = RESUL. MODES. 1 . DEFORMEE_MODALE ;
  74. DEP1 = NNOR DEP1 AVEC (MOTS UR) ;
  75. DEP3D1 = CREER_3d STEEL DEP1 360. 20 1 ;
  76. *==============================================================*
  77. *============= MODE FOURIER 3 ================================*
  78. *==============================================================*
  79. OPTI MODE FOUR 3 ;
  80. RIG1 = RIGI MODCOQ MATCOQ ;
  81. MAS1 = MASS MODCOQ MATCOQ ;
  82. *==============================================================*
  83. * CALCULATION OF THE FREQUENCIES *
  84. * AND *
  85. * EXTRACTION OF SOME RESULTS *
  86. *==============================================================*
  87. RESUL = VIBR PROC (PROG 100.) (CL1 ET RIG1 ) MAS1 IMPR ;
  88. FREQ1 = RESUL.MODES. 1 . FREQUENCE ;
  89. DEP1 = RESUL. MODES. 1 . DEFORMEE_MODALE ;
  90. DEP1 = NNOR DEP1 AVEC (MOTS UR) ;
  91. DEP3D3 = CREER_3d STEEL DEP1 360. 20 3 ;
  92.  
  93. ELIM DEP3D1.MAILLAGE DEP3D3.MAILLAGE 1.E-3 ;
  94. DEP3DT = DEP3D1.DEPLACEMENT + DEP3D3.DEPLACEMENT ;
  95. DEF3DT= DEFO DEP3DT DEP3D1.MAILLAGE 1 ;
  96. TITR 'MAILLAGE 3D MODE FOURIER 1 ET 3 MODE AXIAL 2' ;
  97. SI ( EGA GRAPH 'O') ;
  98. TRAC (0. 12. 30.) CACH DEF3DT ;
  99. FINSI;
  100. POI1 = 0. 3.5 7. ;
  101. POI2 = DEP3D1.MAILLAGE POIN PROC POI1 ;
  102. DEPPOI2 = REDU DEP3DT POI2 ;
  103. UXPOI2 = EXCO UX DEPPOI2 ;
  104. POI22 = EXTR UXPOI2 MAIL ;
  105. PP1 = POI22 POIN PROC (0 0 0) ;
  106. UXPOI22 = EXTR UXPOI2 SCAL PP1 ;
  107. ERR1=100*(ABS(0.14 - UXPOI22)/0.14);
  108. SI (ERR1 < 5 );
  109. ERRE 0;
  110. SINON ;
  111. ERRE 5;
  112. FINSI;
  113. FIN ;
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  

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