Télécharger la plaquette Cast3M


fr en

Télécharger bgfdst.eso

Retour à la liste

Numérotation des lignes : 

bgfdst

  1. C BGFDST    SOURCE    CHAT      05/01/12    21:39:30     5004
  2.       SUBROUTINE BGFDST(IGAU,XE,NGRAF,QSI,
  3.      1                  ETA,HS4,HS5,HS6,BGF)
  4. *
  5. *            CALCUL LA MATRICE BGF RELATIVE AUX
  6. *            EFFETS DE FLEXION
  7. *
  8. *
  9. * ENTREES :
  10. *
  11. *   IGAU :NUMERO DU POINT DE GAUSS
  12. *    XE  :COORDONNEES DES NOEUDS
  13. *   NGRAF:NOMBRE DE COMPOSANTES DE GRADIENTS
  14. *    QSI :COORDONNEES QSI DES POINTS DE GAUSS
  15. *    ETA :COORDONNEES ETA DES POINTS DE GAUSS
  16. * HS4,HS5:FONCTIONS QUI RELIENT LES ROTATIONS (BETA)SK A {U}
  17. * HS6
  18. *
  19. *  SORTIES :
  20. *
  21. *   BGF   :MATRICE BGF
  22. *
  23. *
  24. * AUTEUR : X.Z. SUO FEVRIER 1996
  25. *
  26. ************************************************************************
  27.       IMPLICIT INTEGER(I-N)
  28.       IMPLICIT REAL*8 (A-H,O-Z)
  29.       DIMENSION XE(3,*),QSI(*),ETA(*)
  30.       DIMENSION HS4(*),HS5(*),HS6(*),BGF(NGRAF,*)
  31.       DIMENSION SX(3),SY(3),RL(3),EXX(3),EYY(3)
  32.       DIMENSION HXABS(9),HXORD(9),HYABS(9)
  33.       DIMENSION HYORD(9),B(4,9)
  34. *
  35. *     MATRICE BGF RELATIVE A L'EFFET DE FLEXION
  36. *
  37.       DO 40 K=4,6
  38.        IF(K.EQ.4)THEN
  39.          IJ=1
  40.          I=2
  41.          J=3
  42.        ELSEIF(K.EQ.5)THEN
  43.          IJ=2
  44.           I=3
  45.           J=1
  46.        ELSE
  47.          IJ=3
  48.           I=1
  49.           J=2
  50.        ENDIF
  51.        SX(IJ)=XE(1,I)-XE(1,J)
  52.        SY(IJ)=XE(2,I)-XE(2,J)
  53.        RL(IJ)=SQRT(SX(IJ)*SX(IJ)+SY(IJ)*SY(IJ))
  54.        EXX(IJ)=-SX(IJ)/RL(IJ)
  55.        EYY(IJ)=-SY(IJ)/RL(IJ)
  56.   40   CONTINUE
  57.        AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
  58. *
  59.        CALL DERIVO(SX,SY,RL,QSI(IGAU),ETA(IGAU),HS4,HS5,
  60.      1             HS6,EXX,EYY,HXABS,HXORD,HYABS,HYORD)
  61.        CALL BGRMAT(SX,SY,HXABS,HYABS,HXORD,HYORD,B)
  62. *
  63.        K=2
  64.        KK=0
  65.        DO 327 NPOI=1,3
  66.          DO 527 IY=1,3
  67.            BGF(1,K+IY)=B(1,IY+KK)
  68.            BGF(2,K+IY)=B(2,IY+KK)
  69.            BGF(4,K+IY)=B(3,IY+KK)
  70.            BGF(5,K+IY)=B(4,IY+KK)
  71.   527    CONTINUE
  72.          KK=KK+3
  73.          K=K+6
  74.   327 CONTINUE
  75. *
  76.       RETURN
  77.       END
  78.  
  79.  
  80.  
  81.  

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