Télécharger la plaquette Cast3M


fr en

Télécharger dstgr1.eso

Retour à la liste

Numérotation des lignes : 

dstgr1

  1. C DSTGR1    SOURCE    CHAT      05/01/12    22:56:39     5004
  2.       SUBROUTINE DSTGR1(IGAU,NBNO,XE,NGRA,QSI,ETA,SHPTOT,
  3.      1                  SHP,EXCEN,HS4,HS5,HS6,BGR,DJAC)
  4. *
  5. *            CALCUL LA MATRICE BGR 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. *  EXCEN :EXCENTEMENT DE COQUE
  14. *   NBNO :NOMBRE DE NOEUDS
  15. *  NGRA  :NOMBRE DE COMPOSANTES DE GRADIENTS
  16. *  QSI   :COORDONNEES QSI DES POINTS DE GAUSS
  17. *  ETA   :COORDONNEES ETA DES POINTS DE GAUSS
  18. * SHPTOT :FONCTIONS DE FORME ET LEURS DERIVEE ( EFFET DE MEMBRANE)
  19. *   SHP  :TABLEAU DE TRAVAIL
  20. * HS4,HS5:FONCTIONS QUI RELIENT LES ROTATIONS (BETA)SK A {U}
  21. * HS6
  22. *
  23. *  SORTIES :
  24. *
  25. *   BGR   :MATRICE BGR
  26. *   DJAC  :JACOBIEN
  27. *
  28. *
  29. * AUTEUR : X.Z. SUO FEVRIER 1996
  30. *
  31. ************************************************************************
  32.       IMPLICIT INTEGER(I-N)
  33.       IMPLICIT REAL*8 (A-H,O-Z)
  34.       DIMENSION XE(3,*),QSI(*),ETA(*),SHPTOT(6,3,*),SHP(6,*)
  35.       DIMENSION HS4(*),HS5(*),HS6(*),BGR(NGRA,*)
  36.       DIMENSION SX(3),SY(3),RL(3),EXX(3),EYY(3)
  37.       DIMENSION HXABS(9),HXORD(9),HYABS(9)
  38.       DIMENSION HYORD(9),B(4,9)
  39. *
  40. *     MATRICE BGR RELATIVE A L'EFFET DE MEMBRANE
  41. *
  42.       DO 10 NPOI=1,3
  43.         SHP(1,NPOI)=SHPTOT(1,NPOI,IGAU)
  44.         SHP(2,NPOI)=SHPTOT(2,NPOI,IGAU)
  45.         SHP(3,NPOI)=SHPTOT(3,NPOI,IGAU)
  46.   10  CONTINUE
  47.       CALL JACOBI(XE,SHP,2,3,DJAC)
  48.       K=1
  49.       DO 30 NPOI=1,3
  50.         BGR(1,K  )=SHP(2,NPOI)
  51.         BGR(5,K+1)=SHP(3,NPOI)
  52.         BGR(2,K  )=SHP(3,NPOI)
  53.         BGR(4,K+1)=SHP(2,NPOI)
  54.         K=K+6
  55.   30  CONTINUE
  56. *
  57. *     MODIFICATION SI EXCENTREMENT DE COQUE
  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 BGRMAT(SX,SY,HXABS,HYABS,HXORD,HYORD,B)
  84. *
  85.       K=2
  86.       KK=0
  87.       DO 1027 NPOI=1,NBNO
  88.         DO 1127 IX  =1,3
  89.           BGR(1,K+IX)=EXCEN*B(1,IX+KK)
  90.           BGR(2,K+IX)=EXCEN*B(2,IX+KK)
  91.           BGR(4,K+IX)=EXCEN*B(3,IX+KK)
  92.           BGR(5,K+IX)=EXCEN*B(4,IX+KK)
  93.  1127   CONTINUE
  94.         K=K+6
  95.         KK=KK+3
  96.  1027 CONTINUE
  97. *
  98.       RETURN
  99.       END
  100.  
  101.  
  102.  
  103.  

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