Télécharger d2vlf2.eso

Retour à la liste

Numérotation des lignes :

  1. C D2VLF2 SOURCE BP208322 15/07/22 21:15:10 8586
  2. C DEVLF2 SOURCE LAVARENN 96/10/30 21:18:02 2349
  3. SUBROUTINE D2VLF2(Q1,Q2,FTOTA,NA1,IPALB,IPLIB,XPALB,XVALB,
  4. & NLIAB,XPHILB,JPLIB,NPLB,IDIMB,FTOTB,FTOTBA,
  5. & XPTB,XDT,NPAS,IBASB,IPLSB,INMSB,IORSB,NSB,
  6. & NPLSB,NA2,IND,FEXPSM,NPC1,IERRD,FTEST2,
  7. & FTOTB0,XABSCI,XORDON,NIP,IAROTA,RIGIDE,
  8. & FEXB,XCHPFB)
  9. IMPLICIT INTEGER(I-N)
  10. IMPLICIT REAL*8(A-H,O-Z)
  11. *--------------------------------------------------------------------*
  12. * *
  13. * Opérateur DYNE : algorithme de Fu - de Vogelaere *
  14. * ________________________________________________ *
  15. * *
  16. * Initialisation pour calcul des forces de choc base B. *
  17. * *
  18. * Paramètres: *
  19. * *
  20. * e Q1(.,.) Vecteur des déplacements généralisés.
  21. * e Q2 vecteur des vitesses generalises *
  22. * es FTOTA Forces extérieures totalisées sur la base A. *
  23. * e NA1 Nombre total d'inconnues en base A. *
  24. * e IPALB Renseigne sur le type de la liaison. *
  25. * e IPLIB Tableau contenant les numéros "DYNE" de la liaison. *
  26. * e IAROTA Indique la position des modes de rotation *
  27. * e XPALB Tableau contenant les paramètres de la liaison. *
  28. * es XVALB Tableau contenant les variables internes des liaisons *
  29. * e NLIAB Nombre de liaisons sur la base B. *
  30. * e XPHILB Tableau des vecteurs propres aux points de liaisons. *
  31. * e JPLIB Tableau contenant les numéros "GIBI" des liaisons. *
  32. * e NPLB Nombre total de points intervenant dans les liaisons. *
  33. * e IDIMB Nombre de directions. *
  34. * e IND Indice du pas. *
  35. * e XABSCI Tableau contenant les abscisses de la loi plastique *
  36. * pour la liaison point-point- ... -plastique *
  37. * e XORDON Tableau contenant les ordonnees de la loi plastique *
  38. * pour la liaison point-point- ... -plastique *
  39. * e RIGIDE Vrai si l'on a un corps rigide *
  40. * *
  41. * - FTEST2 Tableau local FTEST de la subroutine D2VLB1 *
  42. * - FTOTB0 Tableau local FTOTB0 de la subroutine D2VLB1 *
  43. * *
  44. * *
  45. * Auteur, date de création: *
  46. * *
  47. * Bertrand BEAUFILS : le 27 juillet 1990. *
  48. * Ianis Politopoulos 22 /02 /93 DEVRCO ---> DEVCOI *
  49. *--------------------------------------------------------------------*
  50. *
  51. INTEGER IPALB(NLIAB,*),IPLIB(NLIAB,*),JPLIB(*)
  52. INTEGER IBASB(*),IPLSB(*),INMSB(*),IORSB(*),IAROTA(*)
  53. REAL*8 XPALB(NLIAB,*),Q1(NA1,*),Q2(NA1,*),FTOTA(NA1,*)
  54. REAL*8 XVALB(NLIAB,4,*),XPHILB(NSB,NPLSB,NA2,*),XPTB(NPLB,4,*)
  55. REAL*8 XDT(*),FTOTB(NPLB,*),FTOTBA(*),FEXPSM(NPLB,NPC1,2,*)
  56. REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*),FEXB(NPLB,2,*)
  57. REAL*8 XCHPFB(2,NLIAB,4,NPLB)
  58. LOGICAL RIGIDE
  59. *
  60. DO 20 IP = 1,NPLB
  61. DO 30 ID = 1,IDIMB
  62. FTOTB(IP,ID) = 0.D0
  63. 30 CONTINUE
  64. * end do
  65. 20 CONTINUE
  66. * end do
  67. IF (IDIMB.EQ.6) THEN
  68. IDIM = 3
  69. ELSE
  70. IDIM = 2
  71. ENDIF
  72. *
  73. * Recombinaison des déplacements aux points de choc
  74. *
  75. CALL D2VCOI(Q1,Q2,XDT,NA1,XPTB,NPLB,XPHILB,NSB,NPLSB,NA2,IDIMB,
  76. & IBASB,IPLSB,INMSB,IORSB,IND,IAROTA)
  77. *
  78. * Initialisation des tableaux contenants les paramètres de liaison
  79. *
  80. CALL DEVLB2(IPLIB,IPALB,XPALB,XPTB,NLIAB,IND,IDIMB,NPLB,
  81. & XABSCI,XORDON,NIP)
  82. *
  83. * Calcul des forces de choc sur base B
  84. *
  85. CALL D2VLB1(FTOTB,XPTB,IPALB,IPLIB,JPLIB,XPALB,XVALB,NLIAB,
  86. & NPLB,IDIMB,XDT,NPAS,IND,FEXPSM,NPC1,IERRD,
  87. & FTEST2,FTOTB0,XABSCI,XORDON,NIP,XCHPFB)
  88. IF (IERRD.NE.0) RETURN
  89. *
  90. * Calcul des moments
  91. *
  92. IF (RIGIDE) THEN
  93. CALL DEVMOM(FTOTB,Q1,FEXB,XPHILB,IAROTA,IBASB,IPLSB,INMSB,
  94. & NA2,NA1,NSB,NPLSB,NPLB,1,IDIM)
  95. ENDIF
  96. *
  97. * Projection des forces base B sur base A
  98. *
  99. CALL DEVPRO(XPHILB,FTOTB,FTOTBA,IBASB,INMSB,IPLSB,IORSB,NSB,
  100. & NPLSB,NA2,IDIMB,NPLB,NA1)
  101. *
  102. * Ajout des forces projetées aux forces extérieures sur base A
  103. *
  104. DO 10 I = 1,NA1
  105. FTOTA(I,IND) = FTOTA(I,IND) + FTOTBA(I)
  106. 10 CONTINUE
  107. * end do
  108. *
  109. END
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  

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