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    CB215821  26/10/26    21:15:14     12682          
  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.       K=K+6
  55.   30  CONTINUE
  56. *
  57. *     MATRICE B RELATIVE A L'EFFET DE FLEXION
  58. *
  59.       DO 40 K=4,6
  60.        IF(K.EQ.4)THEN
  61.          IJ=1
  62.          I=2
  63.          J=3
  64.        ELSEIF(K.EQ.5)THEN
  65.          IJ=2
  66.           I=3
  67.           J=1
  68.         ELSE
  69.          IJ=3
  70.           I=1
  71.           J=2
  72.         ENDIF
  73.        SX(IJ)=XE(1,I)-XE(1,J)
  74.        SY(IJ)=XE(2,I)-XE(2,J)
  75.        RL(IJ)=SQRT(SX(IJ)*SX(IJ)+SY(IJ)*SY(IJ))
  76.        EXX(IJ)=-SX(IJ)/RL(IJ)
  77.        EYY(IJ)=-SY(IJ)/RL(IJ)
  78.   40   CONTINUE
  79.        AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
  80. *
  81.        CALL DERIVO(SX,SY,RL,QSI(IGAU),ETA(IGAU),HS4,HS5
  82.      1              ,HS6,EXX,EYY,HXABS,HXORD,HYABS,HYORD)
  83.        CALL BMATO(SX,SY,HXABS,HYABS,HXORD,HYORD,B)
  84.       K=2
  85. *
  86.       KK=0
  87.       DO 60 NPOI=1,3
  88.       DO 70 IX=1,3
  89.       DO 80 IY=1,3
  90.       BGENE(3+IX,K+IY)=B(IX,IY+KK)
  91.   80  CONTINUE
  92.   70  CONTINUE
  93.       KK=KK+3
  94.       K=K+6
  95.   60  CONTINUE
  96. *
  97.       RETURN
  98.       END
  99.  
  100.  
  101.  
  102.  

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