Télécharger la plaquette Cast3M


fr en

Télécharger dkthsh.eso

Retour à la liste

Numérotation des lignes : 

dkthsh

  1. C DKTHSH    SOURCE    CHAT      05/02/25    21:15:10     5035
  2.       SUBROUTINE DKTHSH(tabuv,tabw,tabrot,DJAC,XSTRS,REL,epai)
  3. C---------------------------------------------------------------------
  4. C   CALCULE LA MATRICE  KSIGMA   POUR L ELEMENT DKT
  5. C
  6. C      ENTREES
  7. C         tabuv(6,3)    = FONCTIONS U U,X U,Y V V,X V,Y
  8. C         tabw(6,9)   = FONCTIONS W W,X  W,Y
  9. C         tabrot(4,9) = FONCTIONS BETAx,x BETAx,y BETA y,x BETAy,y
  10. C         DJAC        = JACOBIEN
  11. C         XSTRS(6)    = CONTRAINTES AU POINT DE GAUSS
  12. C         REL (18,18) = MATRICE KSIGMA
  13. C         epai        = Epaisseur
  14. C      SORTIE
  15. C         REL (18,18) = MATRICE KSIGMA
  16. C
  17. C   RE-ECRIT  LE 2/4/90
  18. C---------------------------------------------------------------------
  19.       IMPLICIT INTEGER(I-N)
  20.       IMPLICIT REAL*8(A-H,O-Z)
  21.       DIMENSION REL(18,*),XSTRS(*),tabuv(6,*),tabw(6,9),tabrot(4,9)
  22.       DATA UNTIER/.33333333333333333D0/
  23.       DATA UNDOUZ/.08333333333333333D0/
  24.       xNXX=XSTRS(1)
  25.       xNYY=XSTRS(2)
  26.       xNXY=XSTRS(3)
  27.       xNXX2=XSTRS(1) * UNDOUZ * epai * epai
  28.       xNYY2=XSTRS(2) * UNDOUZ * epai * epai
  29.       xNXY2=XSTRS(3) * UNDOUZ * epai * epai
  30.  
  31.       DO 16  J=1,3
  32.       JJ=(J-1)*6+1
  33.       DO 17  I=1,3
  34.       II=(I-1)*6+1
  35.         S= xNXX*tabuv(2,I)*tabuv(2,J) +
  36.      .     xNYY*tabuv(3,I)*tabuv(3,J) +
  37.      .     xNXY*(tabuv(2,I)*tabuv(3,J)+tabuv(3,I)*tabuv(2,J))
  38.       S=S*DJAC
  39.       REL(II,JJ)    =REL(II,JJ)    +S
  40.       REL(II+1,JJ+1)=REL(II+1,JJ+1)+S
  41.    17 CONTINUE
  42.    16 CONTINUE
  43.  
  44.       DO 20 L=1,3
  45.       DO 21 M=1,3
  46.       DO 22 I=1,3
  47.       DO 23 K=1,3
  48.          S= xNXX*tabw(2,3*L-3+I)*tabw(2,3*M-3+K) +
  49.      .      xNYY*tabw(3,3*L-3+I)*tabw(3,3*M-3+K) +
  50.      .      xNXY*tabw(2,3*L-3+I)*tabw(3,3*M-3+K) +
  51.      .      xNXY*tabw(3,3*L-3+I)*tabw(2,3*M-3+K) +
  52.      .      xNXX2*tabrot(1,3*L-3+I)*tabrot(1,3*M-3+K) +
  53.      .      xNXX2*tabrot(3,3*L-3+I)*tabrot(3,3*M-3+K) +
  54.      .      xNYY2*tabrot(2,3*L-3+I)*tabrot(2,3*M-3+K) +
  55.      .      xNYY2*tabrot(4,3*L-3+I)*tabrot(4,3*M-3+K) +
  56.      .      xNXY2*tabrot(1,3*L-3+I)*tabrot(2,3*M-3+K) +
  57.      .      xNXY2*tabrot(2,3*L-3+I)*tabrot(1,3*M-3+K) +
  58.      .      xNXY2*tabrot(3,3*L-3+I)*tabrot(4,3*M-3+K) +
  59.      .      xNXY2*tabrot(4,3*L-3+I)*tabrot(3,3*M-3+K)
  60.          S=S*DJAC
  61.          REL(6*L-4+I,6*M-4+K) = REL(6*L-4+I,6*M-4+K) + S
  62.    23 CONTINUE
  63.    22 CONTINUE
  64.    21 CONTINUE
  65.    20 CONTINUE
  66.  
  67.       RETURN
  68.       END
  69.  
  70.  

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