Télécharger inosc1.eso

Retour à la liste

Numérotation des lignes :

  1. C INOSC1 SOURCE CHAT 05/01/13 00:38:32 5004
  2. CAP +AMAX en argument
  3. SUBROUTINE INOSC1(IPT,IPG,DFREQ,XSI,RMAX,TMAX,AMAX)
  4. CAP
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. C
  8. C =====================================================
  9. C =
  10. C CALCUL DU MAXIMUM DE LA REPONSE A UN SIGNAL D'UN =
  11. C OSCILLATEUR AVEC CONDITIONS INITIALES EN VITESSE =
  12. C ET EN DEPLACEMENT NULLES. =
  13. C =
  14. C CREATION : 03/06/87 =
  15. C PROGRAMMEUR : MALAVAL =
  16. C =
  17. C =====================================================
  18. C
  19.  
  20. -INC PPARAM
  21. -INC CCOPTIO
  22. -INC SMLREEL
  23. -INC CCREEL
  24. LOGICAL PASCON
  25. MLREE1=IPT
  26. SEGACT MLREE1
  27. MLREE2=IPG
  28. SEGACT MLREE2
  29. NN=MLREE1.PROG(/1)
  30. C
  31. C LE SIGNAL EST-IL A PAS CONSTANT ?
  32. C
  33. PASCON=.TRUE.
  34. H=MLREE1.PROG(2)-MLREE1.PROG(1)
  35. DO 10 J=1,NN-2
  36. H1=MLREE1.PROG(J+2)-MLREE1.PROG(J+1)
  37. TEST=H1-H
  38. IF (TEST.GT.(1.E-6*H)) THEN
  39. PASCON=.FALSE.
  40. GOTO 11
  41. ENDIF
  42. 10 CONTINUE
  43. C
  44. C CALCUL DE CONSTANTES UTILES
  45. C
  46. 11 W=2*XPI*DFREQ
  47. W2=W*W
  48. A=W*SQRT(1-XSI*XSI)
  49. B=2*XSI/W
  50. C=(2*XSI*XSI-1)/A
  51. D=XSI*W
  52. E=W2/A
  53. F=D/A
  54. VITEP=0.D0
  55. DEPP=0.D0
  56. RMAX=0.D0
  57. CAP
  58. ACCP=0.D0
  59. AMAX=0.D0
  60. CAP
  61. C TMAX=0.D0
  62. C CREATION DE LA BOUCLE
  63. C
  64. S=MLREE2.PROG(1)
  65. H1=MLREE1.PROG(1)
  66. DO 100 I=1,NN-1
  67. C
  68. C DIFFERENCE ENTRE 2 ACCELERATIONS SUCCESSIVES : G
  69. C
  70. T=MLREE2.PROG(I+1)
  71. G=T-S
  72. IF (PASCON.AND.I.NE.1) GOTO 12
  73. C
  74. C DIFFERENCE ENTRE 2 TEMPS SUCCESSIFS : H
  75. C
  76. H2=MLREE1.PROG(I+1)
  77. H=H2-H1
  78. P=EXP(-D*H)
  79. WH=A*H
  80. Q=COS(WH)
  81. R=SIN(WH)
  82. C
  83. C CALCUL EXPLICITE DU DEPLACEMENT
  84. C
  85. PSA=P/A
  86. UNSW2=1/W2
  87. QPFR=Q+F*R
  88. *
  89. 12 O=G/H
  90. DEP=PSA*((A*Q+D*R)*DEPP+R*VITEP)
  91. |+UNSW2*(P*(O*(B*Q+C*R)-S*QPFR)+T-B*O)
  92. VITEP=PSA*((-W2)*R*DEPP+((-D)*R+A*Q)*
  93. |VITEP)+UNSW2*(P*((-O)*QPFR+S*E*R)+O)
  94.  
  95. DEPP=DEP
  96. CAP
  97. ACCP=-W2*DEP-2*D*VITEP
  98. ABSO1=ABS(ACCP)
  99. IF (ABSO1.GT.AMAX) AMAX=ABSO1
  100. CAP
  101. ABSO=ABS(DEPP)
  102. IF (ABSO.GT.RMAX) THEN
  103. RMAX=ABSO
  104. TMAX=MLREE1.PROG(I)
  105. ENDIF
  106. S=T
  107. H1=H2
  108. C
  109. 100 CONTINUE
  110. SEGDES MLREE1
  111. SEGDES MLREE2
  112. RETURN
  113. END
  114.  
  115.  
  116.  
  117.  
  118.  

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