Télécharger cubhr2.eso

Retour à la liste

Numérotation des lignes :

  1. C CUBHR2 SOURCE CHAT 05/01/12 22:32:40 5004
  2. C CUBHR2 SOURCE LOFT 88/06/23 21:07:05
  3. SUBROUTINE CUBHR2(XGENE,DJAC,VKL1,VKL2,LRE,NST,NBNO,IFOU,REL)
  4. C=======================================================================
  5. C
  6. C INTEGRE LE PRODUIT NT KL N
  7. C POUR LES ELEMENTS HOMOGENEISES TRIH
  8. C
  9. C ROUTINE FORTRAN PUR
  10. C
  11. C=======================================================================
  12. C NBDL = NOMBRE DE DDL PAR NOEUD
  13. C INPUT
  14. C XGENE(NBNO,NST)=MATRICE N
  15. C DJAC =JACOBIEN
  16. C VKL1=(COEFPR**2)/(RHO*C**2)*SFLU/SCEL
  17. C VKL2= EI/((EPS**2)*SCEL)
  18. C LRE=TAILLE DE LA MATRICE DE RIGIDITE
  19. C NST=NOMBRE DE COLONNES DE LA MATRICE XGENE
  20. C NBNO=NOMBRE DE NOEUDS DE L ELEMENT CONSIDERE
  21. C IFOU=IFOUR DE CCOPTIO
  22. C OUTPUT
  23. C REL=MATRICE DE RIGIDITE
  24. C=======================================================================
  25. IMPLICIT INTEGER(I-N)
  26. IMPLICIT REAL*8 (A-H,O-Z)
  27. DIMENSION REL(LRE,*),XGENE((NBNO+NBNO),*)
  28. NBDL=LRE/NBNO
  29. C
  30. C
  31. C TERMES EN P*P
  32. C
  33. IX1=0
  34. IY1=0
  35. DO 300 IX=1,LRE ,NBDL
  36. IX1=IX1 + 1
  37. DO 400 IY=1,IX ,NBDL
  38. IY1=IY1 + 1
  39. REL(IX,IY) = REL(IX,IY) + DJAC*VKL1*XGENE(IX1,1)*XGENE(IY1,1)
  40. REL(IY,IX) = REL(IX,IY)
  41. 400 CONTINUE
  42. IY1=0
  43. 300 CONTINUE
  44. C
  45. C
  46. C TERMES EN (UX,RY )*(UX,RY )
  47. C
  48. IX1=0
  49. IY1=0
  50. DO 302 IX=3,LRE ,NBDL
  51. IX1=IX1 + 1
  52. DO 402 IY=3,IX ,NBDL
  53. IY1=IY1 + 1
  54. REL(IX,IY) = REL(IX,IY) + DJAC*VKL2*XGENE(IX1,2)*XGENE(IY1,2)
  55. REL(IX+1,IY) = REL(IX+1,IY) + DJAC*VKL2*XGENE(IX1+8,2)
  56. #*XGENE(IY1,2)
  57. REL(IX,IY+1) = REL(IX,IY+1) + DJAC*VKL2*XGENE(IX1,2)
  58. #*XGENE(IY1+8,2)
  59. REL(IX+1,IY+1) = REL(IX+1,IY+1) + DJAC*VKL2*XGENE(IX1+8,2)*
  60. #XGENE(IY1+8,2)
  61. REL(IY,IX) = REL(IX,IY)
  62. REL(IY,IX+1) = REL(IX+1,IY)
  63. REL(IY+1,IX) = REL(IX,IY+1)
  64. REL(IY+1,IX+1) = REL(IX+1,IY+1)
  65. 402 CONTINUE
  66. IY1=0
  67. 302 CONTINUE
  68. C
  69. C TERMES EN (UY,RX )*(UY,RX )
  70. C
  71. IX1=0
  72. IY1=0
  73. DO 303 IX=5,LRE ,NBDL
  74. IX1=IX1 + 1
  75. DO 403 IY=5,IX ,NBDL
  76. IY1=IY1 + 1
  77. REL(IX,IY) = REL(IX,IY) + DJAC*VKL2*XGENE(IX1,2)*XGENE(IY1,2)
  78. REL(IX+1,IY) = REL(IX+1,IY) + DJAC*VKL2*XGENE(IX1+8,2)
  79. #*XGENE(IY1,2)*(-1.D0)
  80. REL(IX,IY+1) = REL(IX,IY+1) + DJAC*VKL2*XGENE(IX1,2)
  81. #*XGENE(IY1+8,2)*(-1.D0)
  82. REL(IX+1,IY+1) = REL(IX+1,IY+1) + DJAC*VKL2*XGENE(IX1+8,2)*
  83. #XGENE(IY1+8,2)
  84. REL(IY,IX) = REL(IX,IY)
  85. REL(IY,IX+1) = REL(IX+1,IY)
  86. REL(IY+1,IX) = REL(IX,IY+1)
  87. REL(IY+1,IX+1) = REL(IX+1,IY+1)
  88. 403 CONTINUE
  89. IY1=0
  90. 303 CONTINUE
  91. C
  92. C ENDIF
  93. RETURN
  94. END
  95.  
  96.  
  97.  

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