Télécharger d3vpas.eso

Retour à la liste

Numérotation des lignes :

  1. C D3VPAS SOURCE BP208322 18/01/25 21:15:22 9710
  2. C
  3. C D2VPAS SOURCE BP208322 18/01/05 21:15:15 9672
  4. SUBROUTINE D3VPAS(Q1,Q2,Q3,NA1,NPC1,XK,XASM,XM,PDT,T,NPAS,
  5. & FTOTA,FEXA,IFEXA,NPFEXA,NLIAA,NLSA,IPALA,IPLIA,XPALA,XVALA,
  6. & NLIAB,NLSB,NPLB,IDIMB,IPALB,IPLIB,JPLIB,XPALB,XVALB,FTOTB,
  7. & FTOTBA,XPTB,FEXPSM,
  8. & FINERT,IERRD,FTEST,FTOTA0,FTEST2,FTOTB0,WEXT,WINT,
  9. & XABSCI,XORDON,NIP,FTEXB,FEXB,RIGIDE,KTPHI,XCHPFB,
  10. & XOPM1,NB1,NB1K,NB1C,NB1M,B1,B2,B3,B4)
  11. IMPLICIT INTEGER(I-N)
  12. IMPLICIT REAL*8(A-H,O-Z)
  13. *--------------------------------------------------------------------*
  14. * *
  15. * Opérateur DYNE : algorithme de Newmark *
  16. * (Acceleration moyenne par defaut) *
  17. * ________________________________________________ *
  18. * *
  19. * Calcul d'un pas de temps, appel aux s-p spécifiques. *
  20. * *
  21. * Paramètres: *
  22. * *
  23. * es Q1(,) Vecteur des déplacements généralisés *
  24. * es Q2(,) Vecteur des vitesses généralisées *
  25. * es Q3(,) Vecteur des accélérations généralisées *
  26. * es NA1 Nombre total d'inconnues en base A *
  27. * es NPC1 Nombre de pas de calcul - 1 *
  28. * es XK Vecteur des raideurs généralisées *
  29. * es XASM Vecteur des amortissements généralisés *
  30. * es XM Vecteur des masses généralisées *
  31. * e PDT pas de temps courant *
  32. * e T temps courant *
  33. * e NPAS Numéro du pas de temps *
  34. * es FTOTA Forces extérieures totalisées, sur la base A *
  35. * es FEXA Evolution des forces extérieures en base A *
  36. * e FTEXB Evolution des forces extérieures en base B *
  37. * e FEXB Forces extérieures sur la base B, servant au calcul *
  38. * des moments pour les corps rigides. *
  39. * e RIGIDE Vrai si corps rigide, faux sinon *
  40. * es IFEXA Numéro du mode correspondant au point de chargement *
  41. * es NPFEXA Nombre de points de chargement *
  42. * e NLIAA Nombre de liaisons sur la base A *
  43. * e NLSA Nombre de liaisons A en sortie *
  44. * e IPALA Tableau renseignant sur le type de liaison (base A) *
  45. * e IPLIA Tableau contenant les numéros "DYNE" des points *
  46. * e XPALA Tableau contenant les paramètres des liaisons *
  47. * es XVALA Tableau contenant les variables internes de liaison A *
  48. * XPHILB Vecteur des deformees modales *
  49. * e NLIAB Nombre de liaisons sur la base B *
  50. * e NLSB Nombre de liaisons base B en sortie *
  51. * e NPLB Nombre total de points de liaisons (base B) *
  52. * e IDIMB Nombre de directions *
  53. * e IPALB Tableau renseignant sur le type de liaison *
  54. * e IPLIB Tableau contenant les numeros "DYNE" des points *
  55. * e JPLIB Tableau contenant les numeros "GIBI" des points *
  56. * e XPALB Tableau contenant les parametres des liaisons (base B) *
  57. * es XVALB Tableau contenant les variables internes de liaison B *
  58. * FTOTB Forces exterieures totalisees sur la base B *
  59. * e XABSCI Tableau contenant les abscisses de la loi plastique *
  60. * e XORDON Tableau contenant les ordonnees de la loi plastique *
  61. * e NIP Nbr de points dans l'evolution de la loi plastique *
  62. * FTOTBA Forces totales base B projetees base A *
  63. * XPTB Deplacements des points de liaison *
  64. * IBASB Appartenance des points de liaison a une sous-base *
  65. * IPLSB Position du point de liaison dans la sous-base *
  66. * INMSB Nombre de modes dans la sous-base *
  67. * IORSB Position du 1er mode de la sous-base dans ens. modes *
  68. * IAROTA Indique la position des modes de rotation *
  69. * NSB Nombre de sous-bases *
  70. * NPLSB Nombre de points de liaison par sous-base *
  71. * NA2 Nombre d'inconnues dans la sous-base *
  72. * FEXPSM Forces exterieures base B *
  73. * FINERT Forces d'inertie base B *
  74. * IERRD Indicateur d'erreur *
  75. * - FTEST Tableau local FTEST de la subroutine D2VLFA *
  76. * - FTOTA0 Tableau local FTOTA0 de la subroutine D2VLFA *
  77. * - FTEST2 Tableau local FTEST de la subroutine DEVLB1 *
  78. * - FTOTB0 Tableau local FTOTB0 de la subroutine DEVLB1 *
  79. * e,s WEXT travail des forces exterieures
  80. * e,s WINT travail des forces interieures (rigidite et
  81. * amortissement et forces de liaison )
  82. * *
  83. * inspire de D2VPAS *
  84. * *
  85. *--------------------------------------------------------------------*
  86. *
  87. SEGMENT,MTPHI
  88. INTEGER IBASB(NPLB),IPLSB(NPLB),INMSB(NSB),IORSB(NSB)
  89. INTEGER IAROTA(NSB)
  90. REAL*8 XPHILB(NSB,NPLSB,NA2,IDIMB)
  91. ENDSEGMENT
  92. *
  93. INTEGER IFEXA(*),IPALA(NLIAA,*),IPLIA(NLIAA,*)
  94. INTEGER IPALB(NLIAB,*),IPLIB(NLIAA,*),JPLIB(*)
  95. REAL*8 Q1(NA1,*),Q2(NA1,*),Q3(NA1,*)
  96. REAL*8 XVALA(NLIAA,4,*),XPALA(NLIAA,*),XM(NA1,*),XK(NA1,*)
  97. REAL*8 XPALB(NLIAB,*),XVALB(NLIAB,4,*),FEXPSM(NPLB,NPC1,2,*)
  98. REAL*8 XASM(NA1,*),FTOTA(NA1,*),FEXA(NPFEXA,NPC1,*)
  99. REAL*8 FTOTB(NPLB,*),FTOTBA(*),XPTB(NPLB,4,*),FINERT(NA1,*)
  100. REAL*8 WEXT(NA1,2),WINT(NA1,2)
  101. REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*),FEXB(NPLB,2,*)
  102. REAL*8 FTEXB(NPLB,NPC1,2,*),XCHPFB(2,NLIAB,4,NPLB)
  103. REAL*8 XOPM1(NB1,NB1,*),Q2DEMI(NA1),FAMOR(NA1)
  104. *
  105. LOGICAL RIGIDE
  106. LOGICAL LWRITE
  107. *
  108. MTPHI = KTPHI
  109. NSB = XPHILB(/1)
  110. NPLSB = XPHILB(/2)
  111. NA2 = XPHILB(/3)
  112. IVINIT = 0
  113.  
  114. III = 1
  115. *
  116. * Predicteurs des déplacements et vitesses
  117. DO 1 I = 1,NA1
  118. Q1(I,1) = Q1(I,2) + Q2(I,2)*pdt + Q3(i,2)*B1
  119. Q2(I,1) = Q2(I,2) + Q3(I,2)*B2
  120. 1 CONTINUE
  121.  
  122. *
  123. * Totalisation des forces extérieures pour la base A
  124. * pour la fin du pas precedent
  125. CALL D2VFXA(FEXA,IFEXA,FTOTA,NPFEXA,NA1,NPC1,NPAS,FTEXB,FEXB,
  126. & NPLB,IDIMB,RIGIDE)
  127.  
  128. * Ajout des forces de raideur
  129. * F_1 = FEXT_1 - K Q_1
  130. CALL DEVLK0(Q1,XK,FTOTA,NA1,NB1K,III)
  131. * Ajout des forces d'amortissement
  132. * F_1 = FEXT_1 - K Q_1 - C Q_2
  133. CALL DEVLK0(Q2,XASM,FTOTA,NA1,NB1C,III)
  134.  
  135. * Ajout des forces de liaison
  136. IF (NLIAA.NE.0) THEN
  137. CALL D2VLFA(Q1,Q2,FTOTA,NA1,IPALA,IPLIA,XPALA,XVALA,
  138. & NLIAA,PDT,T,NPAS,III,FINERT,IVINIT,FTEST,FTOTA0)
  139. ENDIF
  140. IF (NLIAB.NE.0) THEN
  141. CALL D2VLFB(Q1,FTOTA,NA1,IPALB,IPLIB,XPALB,XVALB,NLIAB,
  142. & XPHILB,JPLIB,NPLB,IDIMB,FTOTB,FTOTBA,XPTB,PDT,T,
  143. & NPAS,IBASB,IPLSB,INMSB,IORSB,NSB,NPLSB,NA2,III,
  144. & FEXPSM,NPC1,IERRD,FTEST2,FTOTB0,
  145. & XABSCI,XORDON,NIP,FEXB,RIGIDE,IAROTA,XCHPFB)
  146. IF (IERRD.NE.0) RETURN
  147. ENDIF
  148.  
  149.  
  150. * Accelerations
  151. c -Cas C ou M pleine ou K...
  152. IF (NB1.NE.1) THEN
  153. * TODO
  154. c -Cas C et M diagonales
  155. ELSE
  156. DO 51 I=1,NA1
  157. Q3(I,1) = FTOTA(I,III) /
  158. & ( XM(I,1) + B3*XASM(I,1) + B4*XK(I,1) - FINERT(I,III) )
  159. 51 CONTINUE
  160. ENDIF
  161.  
  162. * Deplacements et vitesses
  163. DO 7 I = 1,NA1
  164. Q1(I,1) = Q1(I,1) + B4*Q3(I,1)
  165. Q2(I,1) = Q2(I,1) + B3*Q3(I,1)
  166. 7 CONTINUE
  167.  
  168. * calcul des travaux
  169. CALL D2VENE (NA1,III,NPAS,FEXA,Q1,Q2,FTOTA,WEXT,WINT,
  170. & XASM,NPC1,NB1C)
  171.  
  172.  
  173. END
  174.  
  175.  
  176.  

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