Télécharger la plaquette Cast3M


fr en

Télécharger bmfdst.eso

Retour à la liste

Numérotation des lignes : 

bmfdst

  1. C BMFDST    SOURCE    CHAT      05/01/12    21:41:26     5004
  2.       SUBROUTINE BMFDST(IGAU,XE,NSTRS1,QSI,ETA,SHPTOT,SHP,HS4,
  3.      1                  HS5,HS6,BGENE,DJAC)
  4. *
  5. *            CALCUL LA MATRICE B  RELATIVE AUX EFFETS DE
  6. *            MEMBRANE ET DE FLEXION
  7. *
  8. *
  9. * ENTREES :
  10. *
  11. *   IGAU :NUMERO DU POINT DE GAUSS
  12. *    XE  :COORDONNEES DES NOEUDS
  13. *  NSTRS1:NOMBRE DE COMPOSANTES DE CONTRAINTES
  14. *  QSI   :COORDONNEES QSI DES POINTS DE GAUSS
  15. *  ETA    COORDONNEES ETA DES POINTS DE GAUSS
  16. * SHPTOT :FONCTIONS DE FORME ET LEURS DERIVEE ( EFFET DE MEMBRANE)
  17. *   SHP  :TABLEAU DE TRAVAIL
  18. * HS4,HS5:FONCTIONS QUI RELIENT LES ROTATIONS (BETA)SK A {U}
  19. * HS6
  20. *
  21. *  SORTIES :
  22. *
  23. *   BGENE :MATRICE B
  24. *   DJAC  :JACOBIEN
  25. *
  26. *
  27. * AUTEUR : P.DOWLATYARI  JUILLET 91
  28. *
  29. ************************************************************************
  30.       IMPLICIT INTEGER(I-N)
  31.       IMPLICIT REAL*8 (A-H,O-Z)
  32.       DIMENSION XE(3,*),QSI(*),ETA(*),SHPTOT(6,3,*),SHP(6,*)
  33.       DIMENSION HS4(*),HS5(*),HS6(*),BGENE(NSTRS1,*)
  34.       DIMENSION SX(3),SY(3),RL(3),EXX(3),EYY(3)
  35.       DIMENSION HXABS(9),HXORD(9),HYABS(9)
  36.       DIMENSION HYORD(9),B(3,9)
  37. *
  38. *      MATRICE B RELATIVE A L'EFFET DE MEMBRANE
  39. *
  40.       DO 10 NPOI=1,3
  41.       SHP(1,NPOI)=SHPTOT(1,NPOI,IGAU)
  42.       SHP(2,NPOI)=SHPTOT(2,NPOI,IGAU)
  43.       SHP(3,NPOI)=SHPTOT(3,NPOI,IGAU)
  44.   10  CONTINUE
  45.       CALL JACOBI(XE,SHP,2,3,DJAC)
  46.       K=1
  47.       DO 30 NPOI=1,3
  48.       BGENE(1,K  )=SHP(2,NPOI)
  49.       BGENE(1,K+1)=0.D0
  50.       BGENE(2,K  )=0.D0
  51.       BGENE(2,K+1)=SHP(3,NPOI)
  52.       BGENE(3,K  )=SHP(3,NPOI)
  53.       BGENE(3,K+1)=SHP(2,NPOI)
  54.   30  K=K+6
  55. *
  56. *     MATRICE B RELATIVE A L'EFFET DE FLEXION
  57. *
  58.       DO 40 K=4,6
  59.        IF(K.EQ.4)THEN
  60.          IJ=1
  61.          I=2
  62.          J=3
  63.        ELSEIF(K.EQ.5)THEN
  64.          IJ=2
  65.           I=3
  66.           J=1
  67.         ELSE
  68.          IJ=3
  69.           I=1
  70.           J=2
  71.         ENDIF
  72.        SX(IJ)=XE(1,I)-XE(1,J)
  73.        SY(IJ)=XE(2,I)-XE(2,J)
  74.        RL(IJ)=SQRT(SX(IJ)*SX(IJ)+SY(IJ)*SY(IJ))
  75.        EXX(IJ)=-SX(IJ)/RL(IJ)
  76.        EYY(IJ)=-SY(IJ)/RL(IJ)
  77.   40   CONTINUE
  78.        AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
  79. *
  80.        CALL DERIVO(SX,SY,RL,QSI(IGAU),ETA(IGAU),HS4,HS5
  81.      1              ,HS6,EXX,EYY,HXABS,HXORD,HYABS,HYORD)
  82.        CALL BMATO(SX,SY,HXABS,HYABS,HXORD,HYORD,B)
  83.       K=2
  84. *
  85.       KK=0
  86.       DO 60 NPOI=1,3
  87.       DO 70 IX=1,3
  88.       DO 80 IY=1,3
  89.       BGENE(3+IX,K+IY)=B(IX,IY+KK)
  90.   80  CONTINUE
  91.   70  CONTINUE
  92.       KK=KK+3
  93.   60  K=K+6
  94. *
  95.       RETURN
  96.       END
  97.  
  98.  
  99.  

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