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. OPTI EPSI LINEAIRE;
  32. *
  33. * MAILLAGE
  34. *
  35. P1 = 0. 0. 0. ;P2 = 0. 70. 0. ;P3 = 0. 140. 0. ;
  36. L1 = D 1 P1 P2 ; L2 = D 1 P2 P3 ;LTOT = L1 ET L2 ;
  37. *
  38. * MODELE
  39. *
  40. MOD1 = MODELE LTOT MECANIQUE ELASTIQUE BARR ;
  41. *
  42. * CARACTERISTIQUES MATERIELLES ET GEOMETRIQUES
  43. *
  44. MAT1 = MATER MOD1 YOUN 1.E4 NU 0.3 RHO 0.1 ;
  45. CAR1 = CARAC MOD1 SECT 1. ;
  46. *
  47. * MATRICE DE RIGIDITE
  48. *
  49. RIG1 = RIGIDITE MOD1 (MAT1 ET CAR1) ;
  50. *
  51. * MATRICE DE MASSE
  52. *
  53. MAS1 = MASSE MOD1 (MAT1 ET CAR1) ;
  54. *
  55. * CHARGEMENT
  56. *
  57. LI1 = PROG 0. 0.0125 1.0125 1.025 3. ;
  58. LI2 = PROG 0. 1. 1. 0. 0. ;
  59. F1 = FORCE (0. 100. 0.) P1 ;
  60. EV = EVOL MANU T LI1 F(T) LI2 ;
  61. CHA1 = CHAR 'FORC' F1 EV ;
  62. *
  63. * PREPARATION DE LA TABLE DE DONNEES POUR DYNAMIC
  64. *
  65. LIS1 = PROG 1.9125 2.0625 ;
  66. TAB1 = TABLE ;
  67. TAB1.'DEPL' = MANU CHPO LTOT 3 UZ 0. UY 0. UX 0. ;
  68. TAB1.'VITE' = MANU CHPO LTOT 3 UZ 0. UY 0. UX 0. ;
  69. TAB1.'CHAR' = CHA1 ;
  70. TAB1.'RIGI' = RIG1 ;
  71. TAB1.'MASS' = MAS1 ;
  72. TAB1.'FREQ' = 20. ;
  73. TAB1.'INST' = LIS1 ;
  74. TAB2 = DYNAMIC TAB1 ;
  75. *
  76. * COMPARAISON AVEC LES RESULTATS ANALYTIQUES
  77. *
  78. DEP1 = TAB2 . 1 . DEPL ;
  79. U1 = EXTR DEP1 P1 UY ; U1TH = 10.89972 ;
  80. ERR1 = ((U1 - U1TH) / (ABS U1TH)) * 100. ;
  81. DEP2 = TAB2 . 2 . DEPL ;
  82. U2 = EXTR DEP2 P1 UY ; U2TH = 11.73228 ;
  83. ERR2 = ((U2 - U2TH) / (ABS U2TH)) * 100. ;
  84. MESS ' DEPLACEMENT VERTICAL DE LA BASE T= 1.9125' ;
  85. 'MESS''THEORIE 'U1TH'CASTEM 'U1 '%ERREUR 'ERR1 ;
  86. MESS ' DEPLACEMENT VERTICAL DE LA BASE T= 2.0625' ;
  87. 'MESS''THEORIE 'U2TH'CASTEM 'U2 '%ERREUR 'ERR2 ;
  88. *
  89. SI ((MAXI ABS (PROG ERR1 ERR2 )) < 5);
  90. ERRE 0;
  91. SINON;
  92. ERRE 5;
  93. FINSI;
  94. FIN ;
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  

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