Télécharger dbdkt2.eso

Retour à la liste

Numérotation des lignes :

  1. C DBDKT2 SOURCE CHAT 05/02/28 21:15:01 5037
  2. SUBROUTINE DBDKT2(XEL,DD,XDDL,IGAU,XSTRS,SHP1,SHPTOT,BGENE,
  3. * NBNO,LRE,NST)
  4. C======================================================================
  5. C
  6. C AJOUTE DANS LES CONTRAINTES XSTRS LA PARTIE DHO * EPS2
  7. C QUE POUR LES ELEMENTS DKT
  8. C ENTREE
  9. C XEL(3,3) =COORDONNEES LOCALES DU DKT
  10. C DD(6,6) =MATRICE DE HOOKE
  11. C XDDL(18) =D D L LOCAUX DE LA FORMULATION
  12. C IGAU =NUMERO DU POINT DE GAUSS DE 1 A 3
  13. C XSTRS(6) =CONTRAINTES DHO*(EPS1 )
  14. C SORTIE
  15. C XSTRS(6) =CONTRAINTES DHO*(EPS1 + EPS2 )
  16. C
  17. C CODE EBERSOLT AOUT 85
  18. C=======================================================================
  19. IMPLICIT INTEGER(I-N)
  20. IMPLICIT REAL*8(A-H,O-Z)
  21. DIMENSION XEL(3,*),DD(6,*),XSTRS(*),XDDL(*)
  22. DIMENSION EPS2(6),SHP(6,9),ROT(9)
  23. DIMENSION SHP1(6,NBNO),SHPTOT(6,NBNO,*)
  24. DIMENSION BGENE(NST,*)
  25. C
  26. C
  27. DATA XZER,UNDEMI/0.D0,.5D0/
  28. C
  29. DUDX=XZER
  30. DUDY=XZER
  31. DVDX=XZER
  32. DVDY=XZER
  33. DWDX=XZER
  34. DWDY=XZER
  35. eps2(4)=xzer
  36. eps2(5)=xzer
  37. eps2(6)=xzer
  38. C
  39. CALL DKTSHP(IGAU,XEL,SHP,DJAC)
  40. C
  41. C ON EXTRAIT LES DEGRES W ROT X ROT Y
  42. C
  43. IACC= 0
  44. IADD= 2
  45. DO 100 IA=1,3
  46. DO 110 IB=1,3
  47. ROT(IACC+IB)=XDDL(IADD+IB)
  48. 110 CONTINUE
  49. IACC=IACC+3
  50. IADD=IADD+6
  51. 100 CONTINUE
  52. C
  53. DO 200 IA=1,9
  54. DWDX =SHP(2,IA)*ROT(IA) +DWDX
  55. DWDY =SHP(3,IA)*ROT(IA) +DWDY
  56. 200 CONTINUE
  57. C
  58. C--Calcul des derivees de u et v
  59. DO 300 IA=0,2
  60. DUDX=DUDX + BGENE(1,1+6*IA)*XDDL(1+6*IA)
  61. DUDY=DUDY + BGENE(3,1+6*IA)*XDDL(1+6*IA)
  62. DVDX=DVDX + BGENE(3,2+6*IA)*XDDL(2+6*IA)
  63. DVDY=DVDY + BGENE(2,2+6*IA)*XDDL(2+6*IA)
  64. 300 CONTINUE
  65. C
  66. C
  67. EPS2(1)=0.5D0*(DUDX*DUDX+DVDX*DVDX+DWDX*DWDX)
  68. EPS2(2)=0.5D0*(DUDY*DUDY+DVDY*DVDY+DWDY*DWDY)
  69. EPS2(3)=DUDX*DUDY+DVDX*DVDY+DWDX*DWDY
  70. C
  71. C MULTIPLICATION PAR LA MATRICE DE HOOKE
  72. C
  73. DO 700 IB=1,6
  74. DO 700 IC=1,6
  75. XSTRS(IB)=XSTRS(IB)+ DD(IB,IC)*EPS2(IC)
  76. 700 CONTINUE
  77. C
  78. RETURN
  79. END
  80.  
  81.  
  82.  

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