Télécharger bgfdst.eso

Retour à la liste

Numérotation des lignes :

  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