Télécharger devtri.eso

Retour à la liste

Numérotation des lignes :

  1. C DEVTRI SOURCE BP208322 17/07/18 21:15:09 9498
  2. SUBROUTINE DEVTRI(Q1,Q2,Q3,NA1,IINS2,FTOTA,XRES,ICHRES,NRES,
  3. & NCRES,NPRES,XVALA,INULA,NLIAA,NLSA,XRESLA,
  4. & XVALB,INULB,NLIAB,NLSB,XRESLB,ILIREA,ILIREB,
  5. & NTVAR,WEXT,WINT,XCHPFB,NPLB)
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. *--------------------------------------------------------------------*
  9. * *
  10. * Operateur DYNE : algorithme de Fu - de Vogelaere *
  11. * ________________________________________________ *
  12. * *
  13. * Transfert des quantites initiales dans le tableau de *
  14. * sauvegarde, si l'on souhaite les garder ... *
  15. * *
  16. * Param}tres: *
  17. * *
  18. * e Q1(.,.) Vecteur des deplacements generalises *
  19. * e Q2(.,.) Vecteur des vitesses generalises *
  20. * e Q3(.,.) Vecteur des accelerations generalises *
  21. * e NA1 Nombre total d'inconnues en base A *
  22. * e IINS2 Numero du pas de sortie *
  23. * e FTOTA Forces exterieures totalisees en base A *
  24. * es XRES Valeurs des variables sauvegardees *
  25. * e ICHRES Indique quelles variables seront sauvegardees *
  26. * e NRES Nombre de variables ( principales et auxilliaires ) *
  27. * e NCRES Nombre de valeurs prises par les variables *
  28. * e NPRES Nombre de pas de sortie *
  29. * e XVALA Tableau contenant les variables internes de liaison *
  30. * e INULA Tableau contenant les rep}res des liaisons *
  31. * e NLIAA Nombre de liaisons sur la base A *
  32. * e NLSA Nombre de liaisons en sortie bas A *
  33. * es XRESLA Valeurs de variables de liaisons sauvegardees base B *
  34. * e XVALB Tableau contenant les variables internes de liaison *
  35. * e INULB Tableau contenant les rep}res des liaisons *
  36. * e NLIAB Nombre de liaisons sur la base B *
  37. * e NLSB Nombre de liaisons en sortie bas B *
  38. * es XRESLB Valeurs de variables de liaisons sauvegardees base B *
  39. * e,s WEXT travail des forces exterieures
  40. * e,s WINT travail des forces interieures (rigidite et
  41. * amortissement et forces de liaison )
  42. * *
  43. * Auteur, date de creation: *
  44. * *
  45. * Lionel VIVAN, le 1 septembre 1989. *
  46. * *
  47. *--------------------------------------------------------------------*
  48. REAL*8 Q1(NA1,*),Q2(NA1,*),Q3(NA1,*)
  49. REAL*8 XVALA(NLIAA,4,*),XRESLA(NLSA,NPRES,*)
  50. REAL*8 XVALB(NLIAB,4,*),XRESLB(NLSB,NPRES,*)
  51. REAL*8 FTOTA(NA1,*),XRES(NRES,NCRES,*)
  52. REAL*8 XCHPFB(2,NLIAB,4,*)
  53. INTEGER ICHRES(*),INULA(*),INULB(*),ILIREA(NLSA,*),ILIREB(NLSB,*)
  54. *
  55. * Sauvegarde des valeurs initiales
  56. *
  57. IF (NRES.NE.0) THEN
  58. KRES = 0
  59.  
  60. c +deplacement
  61. IF (ICHRES(1).GE.1) THEN
  62. KRES = KRES + 1
  63. DO 10 I=1,NA1
  64. XRES(KRES,I,IINS2) = Q1(I,3)
  65. 10 CONTINUE
  66. ENDIF
  67. c +vitesse
  68. IF (ICHRES(2).GE.1) THEN
  69. KRES = KRES + 1
  70. DO 12 I=1,NA1
  71. XRES(KRES,I,IINS2) = Q2(I,3)
  72. 12 CONTINUE
  73. ENDIF
  74. c +deplacement demi-pas
  75. IF (ICHRES(3).GE.1) THEN
  76. KRES = KRES + 1
  77. DO 14 I=1,NA1
  78. XRES(KRES,I,IINS2) = Q1(I,4)
  79. 14 CONTINUE
  80. ENDIF
  81. c +vitesse demi-pas
  82. IF (ICHRES(4).GE.1) THEN
  83. KRES = KRES + 1
  84. DO 16 I=1,NA1
  85. XRES(KRES,I,IINS2) = Q2(I,4)
  86. 16 CONTINUE
  87. ENDIF
  88. c +acceleration
  89. IF (ICHRES(5).GE.1) THEN
  90. KRES = KRES + 1
  91. DO 18 I=1,NA1
  92. XRES(KRES,I,IINS2) = Q3(I,3)
  93. 18 CONTINUE
  94. ENDIF
  95. c +acceleration demi-pas
  96. IF (ICHRES(6).GE.1) THEN
  97. KRES = KRES + 1
  98. DO 20 I=1,NA1
  99. XRES(KRES,I,IINS2) = Q3(I,4)
  100. 20 CONTINUE
  101. ENDIF
  102.  
  103. c +travail exterieur
  104. IF (ICHRES(7).GE.1) THEN
  105. KRES = KRES + 1
  106. DO 60 I=1,NA1
  107. XRES(KRES,I,IINS2) = 0.d0
  108. 60 CONTINUE
  109. ENDIF
  110. c +travail interieur
  111. IF (ICHRES(8).GE.1) THEN
  112. KRES = KRES + 1
  113. DO 70 I=1,NA1
  114. XRES(KRES,I,IINS2) = 0.d0
  115. 70 CONTINUE
  116. ENDIF
  117.  
  118. ENDIF
  119.  
  120. c +liaisons_A
  121. IF (ICHRES(9).GE.1) THEN
  122. DO 30 IL = 1,NLSA
  123. IIL = INULA(IL)
  124. II = 0
  125. DO 32 IV = 1,NTVAR
  126. IF (ILIREA(IL,IV).EQ.1) THEN
  127. II = II + 1
  128. XRESLA(IL,IINS2,II) = XVALA(IIL,3,IV)
  129. ENDIF
  130. 32 CONTINUE
  131. 30 CONTINUE
  132. ENDIF
  133. c +liaisons_B
  134. IF (ICHRES(10).GE.1) THEN
  135. DO 40 IL = 1,NLSB
  136. IIL = INULB(IL)
  137. II = 0
  138. DO 42 IV = 1,NTVAR
  139. IF (ILIREB(IL,IV).EQ.1) THEN
  140. II = II + 1
  141. XRESLB(IL,IINS2,II) = XVALB(IIL,3,IV)
  142. ELSE
  143. IF (ILIREB(IL,IV).EQ.2) THEN
  144. DO 43 IP=1,NPLB
  145. DO 44 ID=1,2
  146. II = II + 1
  147. XRESLB(IL,IINS2,II) = XCHPFB(ID,IIL,3,IP)
  148. 44 CONTINUE
  149. 43 CONTINUE
  150. ENDIF
  151. *
  152. ENDIF
  153. 42 CONTINUE
  154. 40 CONTINUE
  155. ENDIF
  156. *
  157. END
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  

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