Télécharger la plaquette Cast3M


fr en

Télécharger cubhr2.eso

Retour à la liste

Numérotation des lignes : 

cubhr2

  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