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