Télécharger ftran_test.dgibi

Retour à la liste

Numérotation des lignes :

  1. *******************************************************
  2. ************************************************************************
  3. ************************************************************************
  4. * Cas test de la procedure FTRAN
  5. *******************************************************
  6. * Developpe par :
  7. * Alberto FRAU (alberto.frau[at]cea.fr)
  8. * Benjamin RICHARD (benjamin.richard[at]cea.fr)
  9. *
  10. *******************************************************
  11. * Institution :
  12. * CEA/DEN/DANS/DM2S/SEMT/EMSI
  13. *******************************************************
  14. * Commentaires
  15. *
  16. * On calcule la solution d'un oscillateur 1DDL et on deduit
  17. * la fonction de transfert entre le deplacement au sommet et la
  18. * force appliquée au meme point
  19. *
  20. *******************************************************
  21. *
  22. * Options
  23. OPTI DIME 3 MODE TRID ELEM SEG2;
  24. *
  25. * Sortie graphique
  26. GRAP1 = CHAINE 'N';
  27. *
  28. * Caracteristique oscillatuer
  29. FR_OSC1 = 20.;
  30. OM_OSC1 = ((2)*(PI))*(FR_OSC1);
  31. EP_OSC1 = 0.5;
  32. M_OSC1 = 1000.;
  33. K_OSC1 = ((OM_OSC1)**(2))*(M_OSC1);
  34. *
  35. * Caracteristique force
  36. AMP_FOR1 = 100.;
  37. FR_FOR1 = 2.;
  38. OM_FOR1 = ((2)*(PI))*(FR_FOR1);
  39. FMAX1 = 80.;
  40. TMAX1 = (20.)/(FR_FOR1);
  41. *
  42. * Determination des coef pour la solution analytique
  43. BETA1 = (OM_FOR1)/(OM_OSC1);
  44. A1 = ((-1)*(((AMP_FOR1)/(K_OSC1))*((((2)*(EP_OSC1))*(BETA1))/(((1 -
  45. ((BETA1)**(2)))**(2)) +
  46. (((2)*(EP_OSC1))*(BETA1)**(2))))));
  47. B1 = (((-1)*(OM_FOR1/OM_OSC1))*(((AMP_FOR1)/(K_OSC1))*
  48. (((1 - ((BETA1)**(2)))/
  49. (((1 - ((BETA1)**(2)))**(2)) +
  50. (((2)*(EP_OSC1))*(BETA1)**(2)))))));
  51. *
  52. * Definition des plages temporelle et frequentielle
  53. DT1 = (1)/((2)*(FMAX1));
  54. LT1 = PROG 0. PAS DT1 TMAX1;
  55. *
  56. * Determination du deplacement et de la vitesse
  57. OM_OSC2 = ((OM_OSC1)/(PI))*(180);
  58. OM_FOR2 = ((OM_FOR1)/(PI))*(180);
  59. LU1 = ((EXP(((-1)*(EP_OSC1))*(LT1)))*(((A1)*(COS((OM_OSC2)*(LT1))))
  60. + (((B1)*(SIN((OM_OSC2)*(LT1))))))) + (((AMP_FOR1)/(K_OSC1))*
  61. ((((1 - ((BETA1)**(2)))*(SIN((OM_FOR2)*(LT1)))) - (((EP_OSC1)*
  62. ((BETA1)**(2)))*(COS((OM_FOR2)*(LT1)))))/((((1 -
  63. ((BETA1)**(2)))**(2)) + (((2)*(EP_OSC1))*(BETA1)**(2))))));
  64. LV1 = (((EXP(((-1)*(EP_OSC1))*(LT1)))*(OM_OSC1))*(((A1)*((-1)*
  65. (SIN((OM_OSC2)*(LT1))))) + (((B1)*(COS((OM_OSC2)*(LT1))))))) +
  66. ((((AMP_FOR1)*(OM_FOR1))/(K_OSC1))*((((1 - ((BETA1)**(2)))*
  67. (COS((OM_FOR2)*(LT1)))) + (((EP_OSC1)*((BETA1)**(2)))*
  68. (SIN((OM_FOR2)*(LT1)))))/((((1 - ((BETA1)**(2)))**(2)) +
  69. (((2)*(EP_OSC1))*(BETA1)**(2))))));
  70. *
  71. * Evolution de la force
  72. LF1 = (AMP_FOR1)*(SIN((OM_FOR2)*(LT1)));
  73. EV_FOR1 = EVOL (ROUG) MANU 'Temps [s]' LT1
  74. 'Forc [N]' LF1;
  75. *
  76. EV_DEP1 = EVOL (BLEU) MANU 'Temps [s]' LT1
  77. 'Depl [m]' LU1;
  78. EV_VIT1 = EVOL (VERT) MANU 'Temps [s]' LT1
  79. 'Depl [m]' LV1;
  80. SI ('NEG' GRAP1 'N');
  81. DESS EV_FOR1 TITR 'Evolution de la force';
  82. DESS EV_DEP1 TITR 'Evolution du deplacement';
  83. DESS EV_VIT1 TITR 'Evolution de la vitesse';
  84. FINSI;
  85. *
  86. F_TTR1 = FTRAN EV_DEP1 EV_FOR1 FMAX1 1;
  87. *
  88. SI ('NEG' GRAP1 'N');
  89. DESS F_TTR1 TITR 'Fonction de Transfert';
  90. FINSI;
  91. *
  92. II1 VAL_X DENS1 = MAXI F_TTR1;
  93. *
  94. * Test sur la frequence de l'oscillateur
  95. SI ((ABS(VAL_X - FR_OSC1)) > 0.00001);
  96. ERREUR 5;
  97. FINSI;
  98. *
  99. FIN;
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  

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