Télécharger la plaquette Cast3M


fr en

Télécharger devpro.eso

Retour à la liste

Numérotation des lignes : 

devpro

  1. C DEVPRO    SOURCE    CHAT      05/01/12    22:47:04     5004
  2.       SUBROUTINE DEVPRO(XPHILB,FTOTB,FTOTBA,IBASB,INMSB,IPLSB,IORSB,
  3.      &                  NSB,NPLSB,NA2,IDIMB,NPLB,NA1)
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6. *--------------------------------------------------------------------*
  7. *                                                                    *
  8. *     Op{rateur DYNE : algorithme de Fu - de Vogelaere               *
  9. *     ________________________________________________               *
  10. *                                                                    *
  11. *     Projection des forces base r{elle sur base A.                  *
  12. *                                                                    *
  13. *     Param}tres:                                                    *
  14. *                                                                    *
  15. * e   XPHILB  Tableau des vecteurs propres aux points de liaisons.   *
  16. * e   FTOTB   Tableau des forces sur base B.                         *
  17. *  s  FTOTBA  Tableau des forces base B projet{es sur base A.        *
  18. * e   IBASB   Indique dans quelle sous base appartient le point de   *
  19. *             liaison.                                               *
  20. * e   INMSB   Nombre de modes par sous base.                         *
  21. * e   IORSB   Donne l'indice du premier mode de la sous base dans    *
  22. *             l'ensemble des modes.                                  *
  23. * e   NSB     Nombre total de sous base.                             *
  24. * e   NPLSB   Nombre total de points intervenant dans les liaisons   *
  25. *             d'une sous base.                                       *
  26. * e   NPLB    Nombre total de points intervenant dans les liaisons.  *
  27. * e   IDIMB   Nombre de ddl retenus.                                 *
  28. * e   NA1     Nombre total d'inconnues en base A.                    *
  29. *                                                                    *
  30. *                                                                    *
  31. *     Auteur, date de cr{ation:                                      *
  32. *                                                                    *
  33. *     Lionel VIVAN, le 22 Septembre 1989.                            *
  34. *                                                                    *
  35. *--------------------------------------------------------------------*
  36. *
  37.       INTEGER IBASB(*),INMSB(*),IPLSB(*),IORSB(*)
  38.       REAL*8  XPHILB(NSB,NPLSB,NA2,*),FTOTB(NPLB,*),FTOTBA(*)
  39. *
  40.       DO 10 IN = 1,NA1
  41.          FTOTBA(IN) = 0.D0
  42.  10      CONTINUE
  43. *     end do
  44.       DO 20 IP = 1,NPLB
  45.          ISB  = IBASB(IP)
  46.          IPLB = IPLSB(IP)
  47.          NA3  = INMSB(ISB)
  48.          INA2 = IORSB(ISB) - 1
  49.          DO 30 IN = 1,NA3
  50.             XRET = 0.D0
  51.             DO 40 ID = 1,IDIMB
  52.                XRET = XRET + XPHILB(ISB,IPLB,IN,ID) * FTOTB(IP,ID)
  53.  40            CONTINUE
  54. *           end do
  55.             IN2 = INA2 + IN
  56.             FTOTBA(IN2) = FTOTBA(IN2) + XRET
  57.  30         CONTINUE
  58. *        end do
  59.  20      CONTINUE
  60. *     end do
  61. *
  62.       END
  63.  
  64.  

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