Télécharger la plaquette Cast3M


fr en

Télécharger ftran_test.dgibi

Retour à la liste

Numérotation des lignes : 

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

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