Télécharger bdkt2.eso

Retour à la liste

Numérotation des lignes :

  1. C BDKT2 SOURCE CHAT 05/02/25 21:15:02 5035
  2. SUBROUTINE BDKT2(XEL,XDDL,IGAU,BGENE,XSTRS)
  3. C======================================================================
  4. C
  5. C AJOUTE DANS LES DEFORMATIOND XSTRS LA PARTIE EPS2
  6. C QUE POUR LES ELEMENTS DKT
  7. C ENTREE
  8. C XEL(3,3) =COORDONNEES LOCALES DU DKT
  9. C XDDL(18) =D D L LOCAUX DE LA FORMULATION
  10. C IGAU =NUMERO DU POINT DE GAUSS DE 1 A 3
  11. C XSTRS(6) =DEFORMATIONS (EPS1 )
  12. C BGENE(6,18) =MATRICE BGENE reliant deformations et XDDL
  13. C SORTIE
  14. C XSTRS(6) =DEFORMATIONS (EPS1 + EPS2 )
  15. C
  16. C CODE EBERSOLT AOUT 85
  17. C=======================================================================
  18. IMPLICIT INTEGER(I-N)
  19. IMPLICIT REAL*8(A-H,O-Z)
  20. PARAMETER(XZER=0.D0,UNDEMI=.5D0)
  21. DIMENSION XEL(3,*),XSTRS(*),XDDL(*)
  22. DIMENSION SHP(6,9),ROT(9),BGENE(6,18)
  23. C
  24. DUDX=XZER
  25. DUDY=XZER
  26. DVDX=XZER
  27. DVDY=XZER
  28. DWDX=XZER
  29. DWDY=XZER
  30. C
  31. CALL DKTSHP(IGAU,XEL,SHP,DJAC)
  32. C
  33. C ON EXTRAIT LES DEGRES W ROT X ROT Y
  34. C
  35. IACC= 0
  36. IADD= 2
  37. DO 100 IA=1,3
  38. DO 110 IB=1,3
  39. ROT(IACC+IB)=XDDL(IADD+IB)
  40. 110 CONTINUE
  41. IACC=IACC+3
  42. IADD=IADD+6
  43. 100 CONTINUE
  44. C
  45. DO 200 IA=1,9
  46. DWDX =SHP(2,IA)*ROT(IA) +DWDX
  47. DWDY =SHP(3,IA)*ROT(IA) +DWDY
  48. 200 CONTINUE
  49. C
  50. C--Calcul des derivees de u et v
  51. DO 300 IA=0,2
  52. DUDX=DUDX + BGENE(1,1+6*IA)*XDDL(1+6*IA)
  53. DUDY=DUDY + BGENE(3,1+6*IA)*XDDL(1+6*IA)
  54. DVDX=DVDX + BGENE(3,2+6*IA)*XDDL(2+6*IA)
  55. DVDY=DVDY + BGENE(2,2+6*IA)*XDDL(2+6*IA)
  56. 300 CONTINUE
  57. C
  58. C--Mise a jour des deformations au deuxieme ordre
  59. XSTRS(1)=XSTRS(1)+UNDEMI*(DUDX*DUDX+DVDX*DVDX+DWDX*DWDX)
  60. XSTRS(2)=XSTRS(2)+UNDEMI*(DUDY*DUDY+DVDY*DVDY+DWDY*DWDY)
  61. XSTRS(3)=XSTRS(3)+DUDX*DUDY+DVDX*DVDY+DWDX*DWDY
  62. C
  63. RETURN
  64. END
  65.  
  66.  

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