Télécharger la plaquette Cast3M


fr en

Télécharger dyna12.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  dyna12.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *           Test Dyna12.dgibi: Jeux de données        *
  5. *           ---------------------------------         *
  6. *                                                     *
  7. GRAPH='N';                                                                      
  8. SAUT PAGE ;                                                                     
  9. *******************************************************
  10. *   REPONSE TRANSITOIRE D 'UNE FUSEE -METHODE DIRECT 
  11. *******************************************************
  12. *
  13. * DESCRIPTION DU PROBLEME
  14. * -----------------------
  15. *
  16. * UNE FUSEE EST SOUMISE A UN CHARGEMENT AXIAL POUR UNE 
  17. * DUREE LIMITEE. IL S'AGIT DE TROUVER LES HISTORIQUES 
  18. * DE DEPLACEMENTS AUX EXTREMITES ET AU MILIEU DE LA 
  19. * FUSEE CONSIDERANT LE MOUVEMENT DE CORPS RIGIDE ET 
  20. * LA VIBRATION DE LA STRUCTURE .
  21. *
  22. * REFERENCE :
  23. * ----------
  24. *
  25. * J. S. PRZEMIENIECKI,THEORY OF MATRIX STRUCTURAL 
  26. * ANALYSES (NEW YORK: MCGRAW-HILL BOOK CO.,1968)P.340  
  27. *
  28. *******************************************************
  29. OPTION  ECHO 0 ;
  30. OPTIO DIME 3 ELEM SEG2 ;
  31. *
  32. *  MAILLAGE 
  33. *
  34.  P1 = 0. 0. 0. ;P2 = 0. 70. 0. ;P3 = 0. 140. 0. ;
  35.  L1 = D 1 P1 P2 ; L2 = D 1 P2 P3 ;LTOT = L1 ET L2 ;
  36. *
  37. *  MODELE 
  38. *
  39.  MOD1 = MODELE LTOT MECANIQUE ELASTIQUE BARR ;
  40. *
  41. *  CARACTERISTIQUES MATERIELLES ET GEOMETRIQUES 
  42. *
  43.  MAT1 = MATER MOD1 YOUN 1.E4  NU 0.3 RHO 0.1 ;
  44.  CAR1 = CARAC MOD1 SECT 1. ;
  45. *
  46. *  MATRICE DE RIGIDITE 
  47. *
  48.  RIG1 = RIGIDITE MOD1 (MAT1 ET CAR1) ;
  49. *
  50. *  MATRICE DE MASSE 
  51. *
  52.  MAS1 = MASSE MOD1 (MAT1 ET CAR1)  ;
  53. *
  54. * CHARGEMENT  
  55. *
  56. LI1 = PROG 0. 0.0125 1.0125 1.025  3. ;
  57. LI2 = PROG 0. 1. 1. 0.  0.   ;
  58. F1 = FORCE (0. 100. 0.) P1   ;
  59. EV = EVOL MANU T LI1 F(T) LI2 ;
  60. CHA1 = CHAR 'FORC' F1  EV ;
  61. *
  62. *  PREPARATION DE LA TABLE DE DONNEES POUR DYNAMIC 
  63. *
  64. LIS1 = PROG 1.9125 2.0625  ;
  65. TAB1 = TABLE ;
  66. TAB1.'DEPL' = MANU CHPO LTOT 3 UZ 0. UY 0. UX 0.  ;
  67. TAB1.'VITE' = MANU CHPO LTOT 3 UZ 0. UY 0. UX 0. ;
  68. TAB1.'CHAR' = CHA1 ;
  69. TAB1.'RIGI' = RIG1 ;
  70. TAB1.'MASS' = MAS1 ;
  71. TAB1.'FREQ' = 20. ;
  72. TAB1.'INST' = LIS1 ;
  73. TAB2 = DYNAMIC TAB1 ;
  74. *
  75. * COMPARAISON AVEC LES RESULTATS ANALYTIQUES 
  76. *
  77.  DEP1 = TAB2 . 1 . DEPL ;
  78.  U1 = EXTR DEP1 P1 UY ; U1TH = 10.89972 ;
  79.  ERR1 = ((U1 - U1TH) / (ABS U1TH)) * 100. ;
  80.  DEP2 = TAB2 . 2 . DEPL ;
  81.  U2 = EXTR DEP2 P1 UY ; U2TH = 11.73228  ;
  82.  ERR2 = ((U2 - U2TH) / (ABS U2TH)) * 100. ;
  83. MESS ' DEPLACEMENT VERTICAL DE LA BASE    T= 1.9125'  ;
  84. 'MESS''THEORIE  'U1TH'CASTEM  'U1 '%ERREUR  'ERR1  ;
  85. MESS ' DEPLACEMENT VERTICAL DE LA BASE    T= 2.0625'  ;
  86. 'MESS''THEORIE  'U2TH'CASTEM  'U2 '%ERREUR  'ERR2  ;
  87. *
  88. SI ((MAXI ABS (PROG ERR1 ERR2 )) < 5);
  89. ERRE 0;
  90. SINON;
  91. ERRE 5;
  92. FINSI;
  93. FIN ;
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  

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