Télécharger pospre.eso

Retour à la liste

Numérotation des lignes :

  1. C POSPRE SOURCE CHAT 05/01/13 02:18:16 5004
  2. SUBROUTINE POSPRE(MTRAV,NSPT,NAMRT,IG,ID,TE,NT)
  3. C ===============================================
  4. C=======================================================================
  5. C COMPUTATION OF THE RESPONSE SPECTRA FROM A POWER
  6. C SPECTRA
  7. C
  8. C MTRAV=F(NSPT),S(NSPT) ,ETI(NAMRT), T(NI), RES(NAMRT,NI)
  9. C
  10. C F(1:NSPT)
  11. C S(1:NSPT) =FREQUENCY AND P. SPC. VALUES
  12. C NSPT =NUMBER OF P. SPC. VALUES
  13. C
  14. C ETI(1:NAMRT) =DAMPING VALUES
  15. C NAMRT =NUMBER OF DAMPING VALUES
  16. C
  17. C IG =TYPE OF THE RESPONSE (1 ACCELERATION,2 VELOCITY,
  18. C 3 DISPLACEMENT)
  19. C ID =TYPE OF THE MAXIMUM DISTRIBUTION (1 CRAMER,
  20. C 2 NEWGUMB)
  21. C
  22. C TE =TIME DURATION
  23. C
  24. C T(1:NT) =PERIODS AVAILABLE (INPUT)
  25. C T(1:NT) =PERIODS (OUTPUT)
  26. C NT =NUMBER OF PERIOD AVAILABLE (INPUT)
  27. C NT =NUMBER OF PERIOD (OUTPUT)
  28. C -->IF NT=0 T SHOULD BE COMPUTED AND NT<201
  29. C -->IF NT>0 T IS USED
  30. C
  31. C RES(1:NAMRT,1:NT)
  32. C = RESPONSE SPECTRA FOR VARIOUS DAMPING
  33. C VALUES
  34. C
  35. C FROM A.P. CUSTUMIZED BY P.P. 27/2/90 + 2/4/90 + 15/9/91
  36. C ============================
  37. C=======================================================================
  38. IMPLICIT INTEGER(I-N)
  39. IMPLICIT REAL*8(A-H,O-Z)
  40. C======+++++++++++ bon fontionnement erreur PSRS/PRNS +++++++++++ ======
  41. -INC CCOPTIO
  42. SEGMENT,MTRAV
  43. IMPLIED F(NSPT),S(NSPT) ,ETI(NAMRT), T(NI), RES(NAMRT,NI)
  44. ENDSEGMENT
  45. SEGMENT,MMTRA
  46. IMPLIED XLTIME(NNT,3)
  47. ENDSEGMENT
  48. NNT=1
  49. SEGINI,MMTRA
  50. C======+++++++++++ bon fontionnement erreur PSRS/PRNS +++++++++++ ======
  51. C
  52. DPI=8*ATAN(1.D0)
  53. C
  54. C SELECTION OF PERIODS
  55. C (IF NT >0 THE PERIODS ARE SELECTED BY THE USER)
  56. C ===================
  57. IF(NT.EQ.0)THEN
  58. NT=75
  59. TINF=.04D0
  60. UNPXI=EXP((LOG(TE)-LOG(TINF))/(NT-1))
  61. T(1)=TINF
  62. DO 300 I=2,NT
  63. 300 T(I)=T(I-1)*UNPXI
  64. T(NT)=TE
  65. ENDIF
  66. C COMPUTING OF RESPONSE SPECTRA
  67. C =============================
  68. DO 350 J=1,NAMRT
  69. ETA=ETI(J)
  70. DO 351 I=1,NT
  71. FRQ=1/T(I)
  72. C -MOMENT
  73. C ------
  74. C CALL MOMENT(FRQ,ETA,TE,NSPT,F,S, IG, XL0,XL1,XL2)
  75. CALL MOMENT(FRQ,ETA,TE,NSPT,F,S, IG, XLTIME(1,1),
  76. > XLTIME(1,2),XLTIME(1,3))
  77. C -DISTRIBUTION
  78. C ------------
  79. C CALL DISTRI(XL0,XL1,XL2,TE, ID,DPI, VALMAX)
  80. CALL DISTRI(MMTRA,TE, ID,DPI, VALMAX)
  81. RES(J,I)=VALMAX
  82. 351 CONTINUE
  83. 350 CONTINUE
  84. SEGSUP,MMTRA
  85. RETURN
  86. END
  87.  
  88.  

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