Télécharger rcdst.eso

Retour à la liste

Numérotation des lignes :

  1. C RCDST SOURCE CHAT 05/01/13 02:46:04 5004
  2. SUBROUTINE RCDST(XE,NSTRS,LRE,DDHOMU,HS4,HS5,HS6,
  3. 1 REL,BGENE,IOP)
  4. *
  5. * CALCUL LA MATRICE B ET LA MATRICE DE RIGIDITE
  6. * RELATIVE AUX CISAILLEMENTS TRANSVERSES
  7. *
  8. ************************************************************************
  9. *
  10. * ENTREES :
  11. *
  12. * XE : COORDONNEES DES NOEUDS
  13. * EPAIST :EPAISSEUR
  14. * NSTRS :NOMBRE DE COMPOSANTES DE CONTRAINTES
  15. * LRE :NOMBRE DE DEGRES DE LIBERTE
  16. * DDHOMU :MATRICE DE HOOKE
  17. * OP := 0 LA SORTIE EST LA MATRICE DE RIGIDITE
  18. * = 1 LA SORTIE EST LA MATRICE B
  19. *
  20. * SORTIES :
  21. *
  22. * HS4,HS5,HS6 : LES FONCTIONS QUI RELIE LES TERMES (BETA)SK A {U}
  23. * REL :MATRICE DE RIGIDITE
  24. * BGENE :MATRICE DE B
  25. *
  26. * AUTEUR :P.DOWLATYARI JUILLET 91
  27. ************************************************************************
  28. IMPLICIT INTEGER(I-N)
  29. IMPLICIT REAL*8(A-H,O-Z)
  30. DIMENSION XE(3,*),DDHOMU(NSTRS,*),HS4(*),HS5(*)
  31. DIMENSION HS6(*),REL(LRE,*),BGENE(NSTRS,*)
  32. DIMENSION SX(3),SY(3),RL(3),EXX(3),EYY(3)
  33. DIMENSION DCIS(2,2),EKCIS(9,9),BCISA(2,9),HFX(9),HFY(9)
  34. DIMENSION ASS(6),BSS(6),AXX(3),BXX(3),AYY(3),BYY(3)
  35. *
  36. *
  37. DO 10 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. 10 CONTINUE
  57. AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
  58. *
  59. CALL CISAF(XE,EXX,EYY,ÄDHOM,AKR,ASS,AXX,
  60. 1 AYY,BSS,BXX,BYY)
  61. *
  62. CALL CISAR(ASS,BSS,AXX,AYY,BXX,BYY,RL,EXX,EYY,
  63. 1 HS4,HS5,HS6,HFX,HFY)
  64. *
  65. DO 20 IC=1,9
  66. BCISA(1,IC)=HFX(IC)
  67. BCISA(2,IC)=HFY(IC)
  68. 20 CONTINUE
  69. *
  70. IF(IOP.EQ.0)THEN
  71. DCIS(1,1)=DDHOMU(7,7)
  72. DCIS(1,2)=DDHOMU(7,8)
  73. DCIS(2,1)=DDHOMU(8,7)
  74. DCIS(2,2)=DDHOMU(8,8)
  75. *
  76. CALL ZERO(EKCIS,9,9)
  77. CALL BDBST(BCISA,AIR,DCIS,9,2,EKCIS)
  78. *
  79. * REMPLISSAGE DANS REL
  80. *
  81. DO 30 IL=1,9
  82. IL1=(IL-1)/3
  83. IL2=IL+3*IL1+2
  84. DO 30 IC=1,9
  85. IC1=(IC-1)/3
  86. IC2=IC+IC1*3+2
  87. REL(IL2,IC2)=EKCIS(IL,IC)
  88. 30 CONTINUE
  89. ELSE
  90. *
  91. * REMPLISSAGE DANS BGENE
  92. *
  93. DO 40 IL=1,2
  94. IL1=IL+6
  95. DO 40 IC=1,9
  96. IC1=(IC-1)/3
  97. IC2=IC1*3+IC+2
  98. BGENE(IL1,IC2)=BCISA(IL,IC)
  99. 40 CONTINUE
  100. ENDIF
  101. *
  102. RETURN
  103. END
  104.  
  105.  

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