Télécharger dbco33.eso

Retour à la liste

Numérotation des lignes :

  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