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