Télécharger increm.eso

Retour à la liste

Numérotation des lignes :

  1. C INCREM SOURCE STRU 06/12/14 21:15:18 5611
  2. c Sous-Programme du module MISTRAL1
  3. C --------------------------------------------------------------------------
  4. SUBROUTINE INCREM (TT,FI,FIT,
  5. & SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,
  6. & EPSE,EPSEQ,XX,PSI,
  7. & TTP,FIP,
  8. & SIGP,VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR,
  9. & VEPSE,EPSEQP,XXP,PSIP,
  10. & NDPI,NDVP,NXX,
  11. & DT, KINCR,
  12. & TTF,FIF,FITF,
  13. & SIGF,EPSF,EPSTHF,EPSELF,EPSPF,EPSCRF,
  14. & EPSEF,EPSEQF,XXF,PSIF,
  15. & DELTT,DELFI,DELFIT, DELSIG,DELEPS,DELETH,
  16. & DELEEL,DELEP,DELECR,DELEE,DELEEQ,DELXX,DELPSI)
  17. C --------------------------------------------------------------------------
  18. C Initialisation a 0 des variations DELY des variables Y depuis le
  19. C debut de la periode d'integration, si KINCR = 0.
  20. C Calcul des valeurs YF a l'instant t+DT en fonction des valeurs Y
  21. C a l'instant t et de leurs derivees YP a t ou t+DT/2, si KINCR = 1 ou 2.
  22. C Incrementation des variations DELY des variables Y depuis le debut de la
  23. C periode d'integration en fonction des YP, si KINCR = 2.
  24. C --------------------------------------------------------------------------
  25. IMPLICIT INTEGER (I-N)
  26. IMPLICIT REAL*8 (A-H, O-Z)
  27. c std :
  28. PARAMETER ( NDVPMA = 4 , NXXMAX = 3 )
  29. c std.
  30. DIMENSION SIG(1:*),EPS(1:*),EPSTH(1:*),EPSEL(1:*),EPSCR(1:*),
  31. & EPSP(1:6,0:*),EPSE(0:*),EPSEQ(0:*),XX(1:6,1:*)
  32. DIMENSION SIGP(1:*),VEPS(1:*),VEPSTH(1:*),VEPSEL(1:*),VEPSCR(1:*),
  33. & VEPSP(1:6,0:*),VEPSE(0:*),EPSEQP(0:*),
  34. & XXP(1:6,1:*)
  35. c std :
  36. DIMENSION DSIG(1:6),DEPS(1:6),DEPSTH(1:6),DEPSEL(1:6),DEPSCR(1:6),
  37. & DEPSP(1:6,0:NDVPMA),DEPSE(0:NDVPMA),DEPSEQ(0:NDVPMA),
  38. & DX(1:6,1:NXXMAX)
  39. c std.
  40. DIMENSION SIGF(1:*),EPSF(1:*),EPSTHF(1:*),EPSELF(1:*),EPSCRF(1:*),
  41. & EPSPF(1:6,0:*),EPSEF(0:*),EPSEQF(0:*),
  42. & XXF(1:6,1:*)
  43. DIMENSION DELSIG(1:*),
  44. & DELEPS(1:*),DELETH(1:*),DELEEL(1:*),DELECR(1:*),
  45. & DELEP(1:6,0:*),DELEE(0:*),DELEEQ(0:*),
  46. & DELXX(1:6,1:*)
  47. IF (KINCR.EQ.0) THEN
  48. DELTT = 0.
  49. DELFI = 0.
  50. DELFIT = 0.
  51. DO IJ = 1,6
  52. DELSIG(IJ) = 0.
  53. DELEPS(IJ) = 0.
  54. DELETH(IJ) = 0.
  55. DELEEL(IJ) = 0.
  56. DELECR(IJ) = 0.
  57. END DO
  58. DO IDP = 1-NDPI,NDVP
  59. DO IJ = 1,6
  60. DELEP(IJ,IDP) = 0.
  61. END DO
  62. DELEE(IDP) = 0.
  63. DELEEQ(IDP) = 0.
  64. END DO
  65. DO IXX = 1,NXX
  66. DO IJ = 1,6
  67. DELXX(IJ,IXX) = 0.
  68. END DO
  69. END DO
  70. DELPSI = 0.
  71. RETURN
  72. END IF
  73. DTT = TTP*DT
  74. TTF = TT+DTT
  75. FID = FI
  76. DFI = FIP*DT
  77. FIF = FI+DFI
  78. DFIT = (FID+FIF)/2.*DT
  79. FITF = FIT+DFIT
  80. DO IJ = 1,6
  81. DSIG(IJ) = SIGP(IJ)*DT
  82. SIGF(IJ) = SIG(IJ)+DSIG(IJ)
  83. DEPS(IJ) = VEPS(IJ)*DT
  84. EPSF(IJ) = EPS(IJ)+DEPS(IJ)
  85. DEPSTH(IJ) = VEPSTH(IJ)*DT
  86. EPSTHF(IJ) = EPSTH(IJ)+DEPSTH(IJ)
  87. DEPSEL(IJ) = VEPSEL(IJ)*DT
  88. EPSELF(IJ) = EPSEL(IJ)+DEPSEL(IJ)
  89. DEPSCR(IJ) = VEPSCR(IJ)*DT
  90. EPSCRF(IJ) = EPSCR(IJ)+DEPSCR(IJ)
  91. END DO
  92. DO IDP = 1-NDPI,NDVP
  93. DO IJ = 1,6
  94. DEPSP(IJ,IDP) = VEPSP(IJ,IDP)*DT
  95. EPSPF(IJ,IDP) = EPSP(IJ,IDP)+DEPSP(IJ,IDP)
  96. END DO
  97. DEPSE(IDP) = VEPSE(IDP)*DT
  98. EPSEF(IDP) = EPSE(IDP)+DEPSE(IDP)
  99. DEPSEQ(IDP) = EPSEQP(IDP)*DT
  100. EPSEQF(IDP) = EPSEQ(IDP)+DEPSEQ(IDP)
  101. END DO
  102. DO IXX = 1,NXX
  103. DO IJ = 1,6
  104. DX(IJ,IXX) = XXP(IJ,IXX)*DT
  105. XXF(IJ,IXX) = XX(IJ,IXX)+DX(IJ,IXX)
  106. END DO
  107. END DO
  108. DPSI = PSIP*DT
  109. PSIF = PSI+DPSI
  110. IF (KINCR.EQ.2) THEN
  111. DELTT = DELTT+DTT
  112. DELFI = DELFI+DFI
  113. DELFIT = DELFIT+DFIT
  114. DO IJ = 1,6
  115. DELSIG(IJ) = DELSIG(IJ)+DSIG(IJ)
  116. DELEPS(IJ) = DELEPS(IJ)+DEPS(IJ)
  117. DELETH(IJ) = DELETH(IJ)+DEPSTH(IJ)
  118. DELEEL(IJ) = DELEEL(IJ)+DEPSEL(IJ)
  119. DELECR(IJ) = DELECR(IJ)+DEPSCR(IJ)
  120. END DO
  121. DO IDP = 1-NDPI,NDVP
  122. DO IJ = 1,6
  123. DELEP(IJ,IDP) = DELEP(IJ,IDP)+DEPSP(IJ,IDP)
  124. END DO
  125. DELEE(IDP) = DELEE(IDP)+DEPSE(IDP)
  126. DELEEQ(IDP) = DELEEQ(IDP)+DEPSEQ(IDP)
  127. END DO
  128. DO IXX = 1,NXX
  129. DO IJ = 1,6
  130. DELXX(IJ,IXX) = DELXX(IJ,IXX)+DX(IJ,IXX)
  131. END DO
  132. END DO
  133. DELPSI = DELPSI+DPSI
  134. END IF
  135. RETURN
  136. END
  137.  
  138.  

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