Télécharger hookse.eso

Retour à la liste

Numérotation des lignes :

  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. -INC CCOPTIO
  34. C
  35. DIMENSION VALMAT(*)
  36. DIMENSION DDHOOK(LHOOK,*)
  37. DIMENSION CRIGI(12)
  38. C
  39. C INITIALISATION
  40. C
  41. CALL ZERO(DDHOOK,LHOOK,LHOOK)
  42. C
  43. IF (IGAU.LE.NBGMAT.AND.(IB.LE.NELMAT.OR.NBGMAT.GT.1)) THEN
  44. C
  45. C Cas des poutres et des tuyaux
  46. C
  47. IF (MFR.EQ.7.OR.MFR.EQ.13) THEN
  48. C
  49. IF(IFOU.EQ.2) THEN
  50. C
  51. C LA MATRICE EST PLEINE
  52. C
  53. DDHOOK(1,1)=CRIGI(1)
  54. DDHOOK(1,5)=CRIGI(3)
  55. DDHOOK(5,1)=CRIGI(3)
  56. DDHOOK(1,6)=-CRIGI(2)
  57. DDHOOK(6,1)=-CRIGI(2)
  58. C
  59. DDHOOK(5,5)=CRIGI(6)
  60. DDHOOK(5,6)=-CRIGI(5)
  61. DDHOOK(6,5)=-CRIGI(5)
  62. C
  63. DDHOOK(6,6)=CRIGI(4)
  64. C
  65. DDHOOK(2,2)=CRIGI(8)
  66. DDHOOK(2,4)=-CRIGI(10)
  67. DDHOOK(4,2)=-CRIGI(10)
  68. C
  69. DDHOOK(3,3)=CRIGI(7)
  70. DDHOOK(3,4)=CRIGI(9)
  71. DDHOOK(4,3)=CRIGI(9)
  72. C
  73. DDHOOK(4,4)=CRIGI(11)+CRIGI(12)
  74. C
  75. ELSEIF(IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.-3) THEN
  76. C
  77. C LA MATRICE EST PLEINE
  78. C
  79. DDHOOK(1,1)=CRIGI(1)
  80. DDHOOK(1,3)=-CRIGI(2)
  81. DDHOOK(3,1)=-CRIGI(2)
  82. C
  83. DDHOOK(3,3)=CRIGI(3)
  84. C
  85. DDHOOK(2,2)=CRIGI(4)
  86. C
  87. ELSE
  88. IRET=0
  89. ENDIF
  90. C
  91. C Cas des barres
  92. C
  93. ELSEIF (MFR.EQ.27) THEN
  94. C
  95. YOU=VALMAT(1)
  96. DDHOOK(1,1)=YOU*SECT
  97. C
  98. ELSE
  99. C
  100. C Option non definie
  101. C
  102. IRET=0
  103. ENDIF
  104. ENDIF
  105. C
  106. RETURN
  107. END
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  

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