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

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