Télécharger bdst2.eso

Retour à la liste

Numérotation des lignes :

  1. C BDST2 SOURCE CHAT 05/01/12 21:38:07 5004
  2. SUBROUTINE BDST2(XEL,XDDL,IGAU,BGENE,CMATE,COSA,SINA,XSTRS)
  3. C======================================================================
  4. C
  5. C AJOUTE DANS LES DEFORMATIONS XSTRS LA PARTIE EPS2
  6. C QUE POUR LES ELEMENTS DST
  7. C ENTREE
  8. C XEL(3,3) =COORDONNEES LOCALES DU DST
  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(8) =DEFORMATIONS EPS1
  12. C BGENE(8,18) =RELIE LES DEFORMATIONS AUX DDL
  13. C CMATE ='ISOTROPE' OU 'ORTHOTRO'
  14. C COSA =COSINUS ET SINUS DU L'ANGLE ENTRE l'AXE
  15. C SINA 1 DE L'ORTHOTROPIE ET L'AXE X
  16. CC SORTIE
  17. C XSTRS(6) =(EPS1 + EPS2 )
  18. C
  19. C CODE EBERSOLT AOUT 85 + P DOWLATYARI OCT. 91 POUR ELEMENT DST
  20. C=======================================================================
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8(A-H,O-Z)
  23. DIMENSION XEL(3,*),XSTRS(*),XDDL(*),BGENE(8,*)
  24. DIMENSION EPS2(3),SHP(6,9),ROT(9)
  25. CHARACTER*8 CMATE
  26. C
  27. C
  28. DATA XZER,UNDEMI,TRDEMI/0.D0,.5D0,1.5D0/
  29. C
  30. BETAX=XZER
  31. BETAY=XZER
  32. GAMAXZ=XZER
  33. GAMAYZ=XZER
  34. C
  35. CALL DKTSHP(IGAU,XEL,SHP,DJAC)
  36. C
  37. C ON EXTRAIT LES DEGRES W ROT X ROT Y
  38. C
  39. IACC= 0
  40. IADD= 2
  41. DO 100 IA=1,3
  42. DO 110 IB=1,3
  43. ROT(IACC+IB)=XDDL(IADD+IB)
  44. 110 CONTINUE
  45. IF(CMATE.NE.'ISOTROPE')THEN
  46. RRX=ROT(IACC+2)*COSA+ROT(IACC+3)*SINA
  47. RRY=ROT(IACC+2)*(-SINA)+ROT(IACC+3)*COSA
  48. ROT(IACC+2)=RRX
  49. ROT(IACC+3)=RRY
  50. ENDIF
  51. IACC=IACC+3
  52. IADD=IADD+6
  53. 100 CONTINUE
  54. C
  55. C LES ROT SONT DES W W,X W,Y AUX 3 NOEUDS
  56. C LES XDDL SONT LES W RX RY AUX 3 NOEUDS
  57. C
  58. CCC = ROT(2)
  59. ROT(2)=-ROT(3)
  60. ROT(3)= CCC
  61. C
  62. CCC = ROT(5)
  63. ROT(5)=-ROT(6)
  64. ROT(6)= CCC
  65. C
  66. CCC = ROT(8)
  67. ROT(8)=-ROT(9)
  68. ROT(9)= CCC
  69. C
  70. DO 200 IA=1,9
  71. BETAX =SHP(2,IA)*ROT(IA) +BETAX
  72. BETAY =SHP(3,IA)*ROT(IA) +BETAY
  73. 200 CONTINUE
  74. C
  75. DO 210 IA=1,18
  76. GAMAXZ=BGENE(7,IA)*XDDL(IA)+GAMAXZ
  77. GAMAYZ=BGENE(8,IA)*XDDL(IA)+GAMAYZ
  78. 210 CONTINUE
  79. C
  80. C ON AJOUTE L'EFFET DE CISAILLEMENT
  81. C
  82. IF(IGAU.EQ.1)THEN
  83. DX=XEL(1,2)-XEL(1,1)
  84. DY=XEL(2,2)-XEL(2,1)
  85. ELSEIF(IGAU.EQ.2)THEN
  86. DX=XEL(1,3)-XEL(1,1)
  87. DY=XEL(2,3)-XEL(2,1)
  88. ELSE
  89. DX=XEL(1,3)-XEL(1,2)
  90. DY=XEL(2,3)-XEL(2,2)
  91. ENDIF
  92. XLONG= DX*DX + DY*DY
  93. XLON2= SQRT(XLONG)
  94. CC=DX/XLON2
  95. SS=DY/XLON2
  96. BETAX=BETAX-TRDEMI*CC*CC*GAMAXZ-TRDEMI*SS*CC*GAMAYZ
  97. BETAY=BETAY-TRDEMI*SS*CC*GAMAXZ-TRDEMI*SS*SS*GAMAYZ
  98. C
  99. C BETA EST EN REALITE -BETA DONC:W,X=BETAX+GAMAXZ
  100. C
  101. C
  102. DWDX = BETAX + GAMAXZ
  103. DWDY = BETAY + GAMAYZ
  104. C
  105. XSTRS(1)= DWDX*DWDX*UNDEMI +XSTRS(1)
  106. XSTRS(2)= DWDY*DWDY*UNDEMI +XSTRS(2)
  107. XSTRS(3)= DWDX*DWDY +XSTRS(3)
  108. C
  109. RETURN
  110. END
  111.  
  112.  

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