Télécharger dstgr1.eso

Retour à la liste

Numérotation des lignes :

  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