Télécharger la plaquette Cast3M


fr en

Télécharger pospre.eso

Retour à la liste

Numérotation des lignes : 

pospre

  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.  
  42. -INC PPARAM
  43. -INC CCOPTIO
  44.         SEGMENT,MTRAV
  45.           IMPLIED F(NSPT),S(NSPT) ,ETI(NAMRT), T(NI), RES(NAMRT,NI)
  46.         ENDSEGMENT
  47.         SEGMENT,MMTRA
  48.           IMPLIED XLTIME(NNT,3)
  49.         ENDSEGMENT
  50.         NNT=1
  51.         SEGINI,MMTRA
  52. C======+++++++++++ bon fontionnement erreur PSRS/PRNS +++++++++++ ======
  53. C
  54.         DPI=8*ATAN(1.D0)
  55. C
  56. C       SELECTION OF PERIODS
  57. C       (IF NT >0 THE PERIODS ARE SELECTED BY THE USER)
  58. C       ===================
  59.             IF(NT.EQ.0)THEN
  60.         NT=75
  61.         TINF=.04D0
  62.         UNPXI=EXP((LOG(TE)-LOG(TINF))/(NT-1))
  63.         T(1)=TINF
  64.         DO 300 I=2,NT
  65. 300       T(I)=T(I-1)*UNPXI
  66.         T(NT)=TE
  67.             ENDIF
  68. C       COMPUTING OF RESPONSE SPECTRA
  69. C       =============================
  70.         DO 350 J=1,NAMRT
  71.         ETA=ETI(J)
  72.         DO 351 I=1,NT
  73.         FRQ=1/T(I)
  74. C       -MOMENT
  75. C        ------
  76. C       CALL MOMENT(FRQ,ETA,TE,NSPT,F,S, IG, XL0,XL1,XL2)
  77.         CALL MOMENT(FRQ,ETA,TE,NSPT,F,S, IG, XLTIME(1,1),
  78.      >                           XLTIME(1,2),XLTIME(1,3))
  79. C       -DISTRIBUTION
  80. C        ------------
  81. C       CALL DISTRI(XL0,XL1,XL2,TE, ID,DPI, VALMAX)
  82.         CALL DISTRI(MMTRA,TE, ID,DPI, VALMAX)
  83.         RES(J,I)=VALMAX
  84. 351     CONTINUE
  85. 350     CONTINUE
  86.         SEGSUP,MMTRA
  87.         RETURN
  88.         END
  89.  
  90.  

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