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 - ((BETA1)**(2)))**(2)) + (((2)*(EP_OSC1))*(BETA1)**(2))))));
  45. B1 = (((-1)*(OM_FOR1/OM_OSC1))*(((AMP_FOR1)/(K_OSC1))* (((1 - ((BETA1)**(2)))/ (((1 - ((BETA1)**(2)))**(2)) + (((2)*(EP_OSC1))*(BETA1)**(2)))))));
  46. *
  47. * Definition des plages temporelle et frequentielle
  48. DT1 = (1)/((2)*(FMAX1));
  49. LT1 = PROG 0. PAS DT1 TMAX1;
  50. *
  51. * Determination du deplacement et de la vitesse
  52. OM_OSC2 = ((OM_OSC1)/(PI))*(180);
  53. OM_FOR2 = ((OM_FOR1)/(PI))*(180);
  54. LU1 = ((EXP(((-1)*(EP_OSC1))*(LT1)))*(((A1)*(COS((OM_OSC2)*(LT1)))) + (((B1)*(SIN((OM_OSC2)*(LT1))))))) + (((AMP_FOR1)/(K_OSC1))* ((((1 - ((BETA1)**(2)))*(SIN((OM_FOR2)*(LT1)))) - (((EP_OSC1)* ((BETA1)**(2)))*(COS((OM_FOR2)*(LT1)))))/((((1 - ((BETA1)**(2)))**(2)) + (((2)*(EP_OSC1))*(BETA1)**(2))))));
  55. LV1 = (((EXP(((-1)*(EP_OSC1))*(LT1)))*(OM_OSC1))*(((A1)*((-1)* (SIN((OM_OSC2)*(LT1))))) + (((B1)*(COS((OM_OSC2)*(LT1))))))) + ((((AMP_FOR1)*(OM_FOR1))/(K_OSC1))*((((1 - ((BETA1)**(2)))* (COS((OM_FOR2)*(LT1)))) + (((EP_OSC1)*((BETA1)**(2)))* (SIN((OM_FOR2)*(LT1)))))/((((1 - ((BETA1)**(2)))**(2)) + (((2)*(EP_OSC1))*(BETA1)**(2))))));
  56. *
  57. * Evolution de la force
  58. LF1 = (AMP_FOR1)*(SIN((OM_FOR2)*(LT1)));
  59. EV_FOR1 = EVOL (ROUG) MANU 'Temps [s]' LT1 'Forc [N]' LF1;
  60. *
  61. EV_DEP1 = EVOL (BLEU) MANU 'Temps [s]' LT1 'Depl [m]' LU1;
  62. EV_VIT1 = EVOL (VERT) MANU 'Temps [s]' LT1 'Depl [m]' LV1;
  63. SI ('NEG' GRAP1 'N');
  64. DESS EV_FOR1 TITR 'Evolution de la force';
  65. DESS EV_DEP1 TITR 'Evolution du deplacement';
  66. DESS EV_VIT1 TITR 'Evolution de la vitesse';
  67. FINSI;
  68. *
  69. F_TTR1 = FTRAN EV_DEP1 EV_FOR1 FMAX1 1;
  70. *
  71. SI ('NEG' GRAP1 'N');
  72. DESS F_TTR1 TITR 'Fonction de Transfert';
  73. FINSI;
  74. *
  75. II1 VAL_X DENS1 = MAXI F_TTR1;
  76. *
  77. * Test sur la frequence de l'oscillateur
  78. SI ((ABS(VAL_X - FR_OSC1)) > 0.00001);
  79. ERREUR 5;
  80. FINSI;
  81. *
  82. FIN;
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  

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