Télécharger dstgr2.eso

Retour à la liste

Numérotation des lignes :

  1. C DSTGR2 SOURCE CHAT 05/01/12 22:56:47 5004
  2. SUBROUTINE DSTGR2(XE,NGRA,NSTRS,DDHOMU,HS4,HS5,HS6,BGR)
  3. *
  4. * CALCUL LA MATRICE BGR RELATIVE
  5. * AUX CISAILLEMENTS TRANSVERSES
  6. *
  7. ************************************************************************
  8. *
  9. * ENTREES :
  10. *
  11. * XE : COORDONNEES DES NOEUDS
  12. * NSTRS :NOMBRE DE COMPOSANTES DE CONTRAINTES
  13. * NGRA :NOMBRE DE COMPOSANTES DE GRADIENT
  14. * DDHOMU :MATRICE DE HOOKE
  15. *
  16. * SORTIES :
  17. *
  18. * HS4,HS5,HS6 : LES FONCTIONS QUI RELIE LES TERMES (BETA)SK A {U}
  19. * BGR : MATRICE DE BGR
  20. *
  21. * AUTEUR : X.Z. SUO FEVRIER 96
  22. ************************************************************************
  23. IMPLICIT INTEGER(I-N)
  24. IMPLICIT REAL*8(A-H,O-Z)
  25. DIMENSION XE(3,*),DDHOMU(NSTRS,*),HS4(*),HS5(*),HS6(*)
  26. DIMENSION BGR(NGRA,*),SX(3),SY(3),RL(3),EXX(3),EYY(3)
  27. DIMENSION BCISA(2,9),HFX(9),HFY(9)
  28. DIMENSION ASS(6),BSS(6),AXX(3),BXX(3),AYY(3),BYY(3)
  29. *
  30. DO 10 K=4,6
  31. IF(K.EQ.4)THEN
  32. IJ=1
  33. I=2
  34. J=3
  35. ELSEIF(K.EQ.5)THEN
  36. IJ=2
  37. I=3
  38. J=1
  39. ELSE
  40. IJ=3
  41. I=1
  42. J=2
  43. ENDIF
  44. SX(IJ)=XE(1,I)-XE(1,J)
  45. SY(IJ)=XE(2,I)-XE(2,J)
  46. RL(IJ)=SQRT(SX(IJ)*SX(IJ)+SY(IJ)*SY(IJ))
  47. EXX(IJ)=-SX(IJ)/RL(IJ)
  48. EYY(IJ)=-SY(IJ)/RL(IJ)
  49. 10 CONTINUE
  50. AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
  51. SOMM=SX(IJ)+SY(IJ)+RL(IJ)+EXX(IJ)+EYY(IJ)
  52. *
  53. CALL CISAF(XE,EXX,EYY,DDHOMU,AIR,ASS,AXX,
  54. 1 AYY,BSS,BXX,BYY)
  55. *
  56. CALL CISAR(ASS,BSS,AXX,AYY,BXX,BYY,RL,EXX,EYY,
  57. 1 HS4,HS5,HS6,HFX,HFY)
  58. *
  59. DO 20 IC=1,9
  60. BCISA(1,IC)=HFX(IC)
  61. BCISA(2,IC)=HFY(IC)
  62. 20 CONTINUE
  63. *
  64. * REMPLISSAGE DANS BGR
  65. *
  66. DO 40 IL=1,2
  67. IL1=IL+6
  68. IL2=IL*3
  69. DO 40 IA=1,3
  70. IC1=3*(2*IA-1)
  71. IC2=3*IA-2
  72. BGR(IL1,IC1)=BCISA(IL,IC2)
  73. BGR(IL2,IC1+1)=BCISA(IL,IC2+1)
  74. BGR(IL2,IC1+2)=BCISA(IL,IC2+2)
  75. 40 CONTINUE
  76. *
  77. RETURN
  78. END
  79.  
  80.  
  81.  

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