Esope source files

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