Télécharger la plaquette Cast3M


fr en

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 ;
  36. OPTI MODE FOUR NOHARM ;
  37. OPTI ELEM SEG2;
  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 IMPR   ;
  73. FREQ1 = RESUL.MODES. 1 . FREQUENCE ;
  74. DEP1 = RESUL. MODES. 1 . DEFORMEE_MODALE ;
  75. DEP1 = NNOR DEP1 AVEC (MOTS UR) ;
  76. DEP3D1 = CREER_3d STEEL DEP1 360. 20 1 ;
  77. *==============================================================*
  78. *============= MODE FOURIER 3  ================================*
  79. *==============================================================*
  80. OPTI MODE FOUR 3 ;
  81. RIG1 = RIGI  MODCOQ  MATCOQ  ;
  82. MAS1 = MASS  MODCOQ  MATCOQ  ;
  83. *==============================================================*
  84. *        CALCULATION OF THE FREQUENCIES                        *
  85. *                 AND                                          *
  86. *        EXTRACTION OF SOME RESULTS                            * 
  87. *==============================================================*
  88. RESUL = VIBR PROC (PROG 100.)  (CL1 ET RIG1 )  MAS1 IMPR   ;
  89. FREQ1 = RESUL.MODES. 1 . FREQUENCE ;
  90. DEP1 = RESUL. MODES. 1 . DEFORMEE_MODALE ;
  91. DEP1 = NNOR DEP1 AVEC (MOTS UR) ;
  92. DEP3D3 = CREER_3d STEEL DEP1 360. 20 3 ;
  93.  
  94. ELIM DEP3D1.MAILLAGE DEP3D3.MAILLAGE 1.E-3 ;
  95. DEP3DT = DEP3D1.DEPLACEMENT + DEP3D3.DEPLACEMENT ;
  96. DEF3DT= DEFO DEP3DT DEP3D1.MAILLAGE 1 ;
  97. TITR 'MAILLAGE 3D MODE FOURIER 1 ET 3 MODE AXIAL 2' ;
  98. SI ( EGA GRAPH  'O') ; 
  99. TRAC (0. 12. 30.) CACH  DEF3DT ;
  100. FINSI;
  101. POI1 = 0. 3.5 7. ;
  102. POI2 = DEP3D1.MAILLAGE POIN PROC POI1 ;
  103. DEPPOI2 = REDU DEP3DT POI2 ;
  104. UXPOI2 = EXCO UX DEPPOI2 ;
  105. POI22 = EXTR UXPOI2 MAIL ;
  106. PP1 = POI22 POIN PROC (0 0 0) ;
  107. UXPOI22 = EXTR UXPOI2 SCAL PP1 ;
  108. ERR1=100*(ABS(0.14 - UXPOI22)/0.14);
  109. SI (ERR1 < 5 );
  110. ERRE 0;
  111. SINON ;
  112. ERRE 5;
  113. FINSI;
  114. FIN ;
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  

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