Télécharger la plaquette Cast3M


fr en

Télécharger dbco33.eso

Retour à la liste

Numérotation des lignes : 

dbco33

  1. C DBCO33    SOURCE    BP208322  17/03/01    21:16:56     9325           
  2.       SUBROUTINE DBCO33(XE,DD,XDDL,EPS2,SHP,BPSS,XEL,XLOC,XSTRS)
  3. C======================================================================
  4. C
  5. C  AJOUTE DANS LES CONTRAINTES XSTRS LA PARTIE DHO * EPS2
  6. C                  QUE POUR LES ELEMENTS COQ3
  7. C  ENTREE
  8. C    XE(3,3)           =COORDONNEES GLOBALES DU COQ3
  9. C    DD(6,6)            =MATRICE DE HOOKE
  10. C    XDDL(18)          =D.D.L.GLOBAUX DE LA FORMULATION
  11. C    XSTRS(6)           =CONTRAINTES  DHO*(EPS1        )
  12. C  TRAVAIL
  13. C    EPS2(3)            = EPSILON 2
  14. C    SHP(3,9)           = FONCTIONS DE FORME
  15. C    XEL(3,3)           = COORDONNEES LOCALES
  16. C    XLOC(18)         = D.D.L. LOCAUX
  17. C  SORTIE
  18. C    XSTRS(6)           =CONTRAINTES DHO*(EPS1 + EPS2 )
  19. C
  20. C   CODE EBERSOLT AOUT  85
  21. C=======================================================================
  22.       IMPLICIT INTEGER(I-N)
  23.       IMPLICIT REAL*8(A-H,O-Z)
  24.       DIMENSION XE(3,*),DD(6,*),XSTRS(*),XDDL(*)
  25.       DIMENSION EPS2(*),SHP(3,*),BPSS(3,*)
  26.       DIMENSION XLOC(*),XEL(3,*)
  27. C
  28. C
  29.       DATA XZER,UNDEMI,UNTIER/0.D0,.5D0,.333333333333333333D0/
  30. C
  31.       DWDX=XZER
  32.       DWDY=XZER
  33. C
  34.       CALL VPAST(XE,BPSS)
  35.       CALL VCORLC(XE,XEL,BPSS)
  36. *      WRITE(6,*) ' BPSS', ((BPSS( IOU,IYU),IOU=1,3),IYU=1,3)
  37. *      WRITE(6,*) 'XEL' , ((XEL(iou,iyu),iou=1,3), iyu=1,3)
  38.       CALL MATVEC(XDDL,XLOC,BPSS,6)
  39.       CALL COQ3SH(UNTIER,UNTIER,XEL,SHP)
  40. *      write(6,*) ((shp(iou,iyu), iou=1,3),iyu=1,9)
  41. *      write(6,*) (xloc(iou),iou=1,18)
  42.       DO 200 IA=1,3
  43.       IA1 = 2 +(IA-1)*6
  44.       IA2 =    (IA-1)*3
  45.       DO 200 IB=1,3
  46.       DWDX =SHP(2,IA2+IB)*XLOC(IA1+IB) +DWDX
  47.       DWDY =SHP(3,IA2+IB)*XLOC(IA1+IB) +DWDY
  48.   200 CONTINUE
  49. C
  50.       EPS2(1)= DWDX*DWDX*UNDEMI
  51.       EPS2(2)= DWDY*DWDY*UNDEMI
  52.       EPS2(3)= DWDX*DWDY
  53. **      WRITE(6,888) EPS2(1),EPS2(2),EPS2(3),DWDX,DWDY
  54. *  888 FORMAT(4X,' DEBUGG EPS2 ',5E12.5)
  55. C
  56. C     MULTIPLICATION PAR LA MATRICE DE HOOKE
  57. C
  58.       DO 300 IB=1,3
  59.       DO 300 IC=1,3
  60.       XSTRS(IB)=XSTRS(IB)+ DD(IB,IC)*EPS2(IC)
  61.   300 CONTINUE
  62. C
  63.       RETURN
  64.       END
  65.  
  66.  
  67.  
  68.  

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