Télécharger la plaquette Cast3M


fr en

Télécharger dbdkt2.eso

Retour à la liste

Numérotation des lignes : 

dbdkt2

  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