Télécharger la plaquette Cast3M


fr en

Télécharger devlf2.eso

Retour à la liste

Numérotation des lignes : 

devlf2

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

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