Télécharger devlfb.eso

Retour à la liste

Numérotation des lignes :

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

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