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. OPTI EPSI LINEAIRE;
  38. *
  39. *================================================================*
  40. * GEOMETRY *
  41. *==============================================================*
  42. *
  43. P1 = 0. 0. ;
  44. P2 = 3.5 0.0 ;
  45. P3 = 3.5 7.0 ;
  46. *
  47. LI1 = P1 DROIT 3 P2 ;
  48. STEEL = LI1 ET (P2 DROIT 8 P3) ;
  49. *
  50. *==============================================================*
  51. * MODELE - MATERIAU - RIGIDITE - MASSE *
  52. *==============================================================*
  53. MODCOQ=MODE STEEL MECANIQUE COQ2 ;
  54. MATCOQ=MATE MODCOQ RHO 8.0E3 YOUN 2.E11 NU 0.3 ;
  55. CARCOQ=CARA MODCOQ EPAI 0.035 ;
  56. MATCOQ=MATCOQ ET CARCOQ;
  57. *==============================================================*
  58. * BOUNDARY CONDITION *
  59. *==============================================================*
  60. CL1 = BLOQ DEPL ROTA LI1 ;
  61. *==============================================================*
  62. *============= MODE FOURIER 1 ================================*
  63. *==============================================================*
  64. OPTI MODE FOUR 1 ;
  65. RIG1 = RIGI MODCOQ MATCOQ ;
  66. MAS1 = MASS MODCOQ MATCOQ ;
  67. *==============================================================*
  68. * CALCULATION OF THE FREQUENCIES *
  69. * AND *
  70. * EXTRACTION OF SOME RESULTS *
  71. *==============================================================*
  72. RESUL = VIBR proche (prog 150.) (CL1 ET RIG1 ) MAS1
  73. TBAS IMPR ;
  74. FREQ1 = RESUL.MODES. 1 . FREQUENCE ;
  75. DEP1 = RESUL. MODES. 1 . DEFORMEE_MODALE ;
  76. NNOR DEP1 AVEC (MOTS UR) ;
  77. DEP3D1 = CREER_3d STEEL DEP1 360. 20 1 ;
  78. *==============================================================*
  79. *============= MODE FOURIER 3 ================================*
  80. *==============================================================*
  81. OPTI MODE FOUR 3 ;
  82. RIG1 = RIGI MODCOQ MATCOQ ;
  83. MAS1 = MASS MODCOQ MATCOQ ;
  84. *==============================================================*
  85. * CALCULATION OF THE FREQUENCIES *
  86. * AND *
  87. * EXTRACTION OF SOME RESULTS *
  88. *==============================================================*
  89. RESUL = VIBR PROC (PROG 100.) (CL1 ET RIG1 ) MAS1
  90. TBAS IMPR ;
  91. FREQ1 = RESUL.MODES. 1 . FREQUENCE ;
  92. DEP1 = RESUL. MODES. 1 . DEFORMEE_MODALE ;
  93. NNOR DEP1 AVEC (MOTS UR) ;
  94. DEP3D3 = CREER_3d STEEL DEP1 360. 20 3 ;
  95.  
  96. ELIM DEP3D1.MAILLAGE DEP3D3.MAILLAGE 1.E-3 ;
  97. DEP3DT = DEP3D1.DEPLACEMENT + DEP3D3.DEPLACEMENT ;
  98. DEF3DT= DEFO DEP3DT DEP3D1.MAILLAGE 1 ;
  99. TITR 'MAILLAGE 3D MODE FOURIER 1 ET 3 MODE AXIAL 2' ;
  100. SI ( EGA GRAPH 'O') ;
  101. TRAC (0. 12. 30.) CACH DEF3DT ;
  102. FINSI;
  103. POI1 = 0. 3.5 7. ;
  104. POI2 = DEP3D1.MAILLAGE POIN PROC POI1 ;
  105. DEPPOI2 = REDU DEP3DT POI2 ;
  106. UXPOI2 = EXCO UX DEPPOI2 ;
  107. POI22 = EXTR UXPOI2 MAIL ;
  108. PP1 = POI22 POIN PROC (0 0 0) ;
  109. UXPOI22 = EXTR UXPOI2 SCAL PP1 ;
  110. ERR1=100*(ABS(0.14 - UXPOI22)/0.14);
  111. SI (ERR1 < 5 );
  112. ERRE 0;
  113. SINON ;
  114. ERRE 5;
  115. FINSI;
  116. FIN ;
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  

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