Télécharger la plaquette Cast3M


fr en

Télécharger devlfb.eso

Retour à la liste

Numérotation des lignes : 

devlfb

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

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