Télécharger hookse.eso

Retour à la liste

Numérotation des lignes :

hookse
  1. C HOOKSE SOURCE BP208322 15/06/22 21:19:17 8543
  2. SUBROUTINE HOOKSE(VALMAT,IB,IGAU,MFR,CRIGI,IFOU,
  3. + NBGMAT,NELMAT,SECT,LHOOK,DDHOOK,IRET)
  4. C
  5. C----------------------------------------------------------------------
  6. C
  7. C Calcul de la matrice de HOOKE dans le cas d'un
  8. C matériau homogeneise
  9. C
  10. C Entrees:
  11. C --------
  12. C VALMAT tableau de materiau
  13. C IB numero de l'element
  14. C IGAU numero du point de Gauss
  15. C MFR numero de formulation
  16. C CRIGI ELEMENT DE REDUCTION DE LA RIGIDITE
  17. C IFOU numero d'harmonique de Fourier
  18. C NBGMAT, NELMAT tailles des tableaux
  19. C SECT SECTION DE L'ELEMENT IB (<> 0 SI MFR.EQ.27)
  20. C LHOOK taille de la matrice de HOOKE
  21. C
  22. C Sorties:
  23. C --------
  24. C DDHOOK matrice de HOOKE
  25. C IRET = 1 si option existante, 0 sinon
  26. C
  27. C---------------------------------------------------------------------
  28. C
  29. IMPLICIT INTEGER(I-N)
  30. IMPLICIT REAL*8(A-H,O-Z)
  31. C
  32. -INC CCHAMP
  33.  
  34. -INC PPARAM
  35. -INC CCOPTIO
  36. C
  37. DIMENSION VALMAT(*)
  38. DIMENSION DDHOOK(LHOOK,*)
  39. DIMENSION CRIGI(12)
  40. C
  41. C INITIALISATION
  42. C
  43. CALL ZERO(DDHOOK,LHOOK,LHOOK)
  44. C
  45. IF (IGAU.LE.NBGMAT.AND.(IB.LE.NELMAT.OR.NBGMAT.GT.1)) THEN
  46. C
  47. C Cas des poutres et des tuyaux
  48. C
  49. IF (MFR.EQ.7.OR.MFR.EQ.13) THEN
  50. C
  51. IF(IFOU.EQ.2) THEN
  52. C
  53. C LA MATRICE EST PLEINE
  54. C
  55. DDHOOK(1,1)=CRIGI(1)
  56. DDHOOK(1,5)=CRIGI(3)
  57. DDHOOK(5,1)=CRIGI(3)
  58. DDHOOK(1,6)=-CRIGI(2)
  59. DDHOOK(6,1)=-CRIGI(2)
  60. C
  61. DDHOOK(5,5)=CRIGI(6)
  62. DDHOOK(5,6)=-CRIGI(5)
  63. DDHOOK(6,5)=-CRIGI(5)
  64. C
  65. DDHOOK(6,6)=CRIGI(4)
  66. C
  67. DDHOOK(2,2)=CRIGI(8)
  68. DDHOOK(2,4)=-CRIGI(10)
  69. DDHOOK(4,2)=-CRIGI(10)
  70. C
  71. DDHOOK(3,3)=CRIGI(7)
  72. DDHOOK(3,4)=CRIGI(9)
  73. DDHOOK(4,3)=CRIGI(9)
  74. C
  75. DDHOOK(4,4)=CRIGI(11)+CRIGI(12)
  76. C
  77. ELSEIF(IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.-3) THEN
  78. C
  79. C LA MATRICE EST PLEINE
  80. C
  81. DDHOOK(1,1)=CRIGI(1)
  82. DDHOOK(1,3)=-CRIGI(2)
  83. DDHOOK(3,1)=-CRIGI(2)
  84. C
  85. DDHOOK(3,3)=CRIGI(3)
  86. C
  87. DDHOOK(2,2)=CRIGI(4)
  88. C
  89. ELSE
  90. IRET=0
  91. ENDIF
  92. C
  93. C Cas des barres
  94. C
  95. ELSEIF (MFR.EQ.27) THEN
  96. C
  97. YOU=VALMAT(1)
  98. DDHOOK(1,1)=YOU*SECT
  99. C
  100. ELSE
  101. C
  102. C Option non definie
  103. C
  104. IRET=0
  105. ENDIF
  106. ENDIF
  107. C
  108. RETURN
  109. END
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  

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