Télécharger la plaquette Cast3M


fr en

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