Télécharger d2vtri.eso

Retour à la liste

Numérotation des lignes :

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

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