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