Télécharger bsigco.eso

Retour à la liste

Numérotation des lignes :

  1. C BSIGCO SOURCE CHAT 05/01/12 21:42:54 5004
  2. SUBROUTINE BSIGCO(XXE,DAR,FFE,RPL,FT,TWL,RJ,XR,RC,
  3. * FU,P,FF,XE,YE,ZE)
  4. C
  5. C ***** ELEMENT COQ3 ******
  6. C CE SUB A COMMME BUT DE CALCULER LES FORCES ( EXPRIMEES DANS LE
  7. C REOERE GLOBALE) QUI EQUILIBRENT UN CHAMP DE CONTRAINTES POUR
  8. C
  9. C EN ENTREE - EP : EPAISSEUR DE LA COQUE
  10. C XXE(I,J) : IEME COORDONNEE DU JEME NOEUD DE L'ELEMENT
  11. C DAR(6) : CONTRAINTES - MEMBRANE PUIS FLEXION
  12. C
  13. C EN SORTIE - FFE(I,J) : IEME COMPOSANTE DU CHAMP DE FORCE A
  14. C APPLIQUER AU JEME NOEUD DE L'ELEMENT
  15. C
  16. C CE SUB VIENT DE TRICO VIA BILBO.
  17. C
  18. C LA LONGUEUR DU SEGMENT DE TRAVAIL NECESSAIRE EST 151
  19. C ATTENTION RPL ET FT SONT EN EQUIVALLENCE AINSI QUE FF(1),XE(1)
  20. C FF(4),YE(1) FF(7),ZE(1)
  21. C
  22. IMPLICIT INTEGER(I-N)
  23. IMPLICIT REAL*8(A-H,O-Z)
  24. DIMENSION FFE(6,2),XXE(3,3),DAR(6)
  25. DIMENSION TWL(6)
  26. DIMENSION RL(9,9),RP(9,9),RPL(9,9),JLIM(2,9),XL(6,3),RJ(2,2),
  27. 1 XR(3,9),FF(9),RC(3,3)
  28. DIMENSION P(3,3),XE(3),YE(3),ZE(3)
  29. DIMENSION FT(18)
  30. DIMENSION FU(6)
  31. DATA JLIM/1,1,2,3,2,3,4,4,5, 6,5,6,7,7,8,9,8,9/
  32. DATA RL/ 0., 0., 0., 0. ,0. ,0. , 1. ,0. ,0. ,
  33. 1 1., 0., 0., 0. ,0. ,0. , -1. ,0. ,0. ,
  34. 2 0., 0., 0., 1. ,0. ,0. , -1. ,0. ,0. ,
  35. 3 0., -0.5,0.5,0.,0.5,-0.5, 0., 0., 0.,
  36. 4 0., 0., 0., 0., 0., -0.5, 0., 0., 0.5,
  37. 6 0.,-0.5, 0., 0., 0., 0., 0., 0.5, 0.,
  38. 7 -1., 0.5, -0.5, 1., 0.5, -0.5, 0., 0., 0.,
  39. 8 0., 0., 0., -1., 0., 0.5, 1., 0., 0.5,
  40. 9 1., -0.5, 0., 0., 0., 0., -1., -0.5, 0./
  41. DATA RP/ 1., 0., 0., 0., 0., 0., 0., 0., 0.,
  42. 2 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  43. 3 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  44. 4 0., 0., 0., 1., 0., 0., 0., 0., 0.,
  45. 5 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  46. 6 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  47. 7 0., 0., 0., 0., 0., 0., 1., 0., 0.,
  48. 8 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  49. 9 0., 0., 0., 0., 0., 0., 0., 0., 0./
  50. DATA XL/0.,0.,-1., -0.33333333,0.33333333,0.,
  51. 2 0., -1., 0., 0.33333333,0.,-0.33333333,
  52. 3 0.5,-0.5,-0.5,0.,0.33333333,-0.33333333/
  53. DO 185 IP=1,3
  54. XE(IP)=XXE(1,IP)
  55. YE(IP)=XXE(2,IP)
  56. 185 ZE(IP)=XXE(3,IP)
  57. CALL PASSA(XE,YE,ZE,P,X13,X23,Y13,Y23)
  58. 184 CONTINUE
  59. IELE=3
  60. DO 16 IP=1,3
  61. 16 TWL(IP)=DAR(IP+3)
  62. DO 160 IP=1,3
  63. 160 TWL(3+IP)=DAR(IP)
  64. SO= 0.5D0
  65. FU(1)=SO*(TWL(4)*Y23-TWL(6)*X23)
  66. FU(2)=SO*(-TWL(5)*X23+TWL(6)*Y23)
  67. FU(3)=SO*(-TWL(4)*Y13+TWL(6)*X13)
  68. FU(4)=SO*(TWL(5)*X13-TWL(6)*Y13)
  69. FU(5)=SO*TWL(6)*(X23-X13)
  70. FU(6)=SO*(X23-X13)*TWL(5)
  71. RJ(1,1)=X13
  72. RJ(2,1)=X23
  73. RJ(1,2)=Y13
  74. RJ(2,2)=Y23
  75. DT=X13*Y23-Y13*X23
  76. SURF2= ABS(DT)
  77. DT=-1./ABS(DT)
  78. DO 1 K=1,3
  79. IK=3*K-2
  80. DO 1 I=1,2
  81. II=IK+I
  82. DO 1 J=1,2
  83. IJ=IK+J
  84. 1 RP(IJ,II)=RJ(J,I)
  85. DO 2 J=1,9
  86. K1=JLIM(1,J)
  87. K2=JLIM(2,J)
  88. DO 2 I=4,9
  89. S=0.D0
  90. DO 3 K=K1,K2
  91. 3 S=S+RL(K,I)*RP(K,J)
  92. RPL(I,J)=S
  93. 2 CONTINUE
  94. DO 4 J=1,3
  95. DO 4 I=1,9
  96. S=0.D0
  97. DO 5 K=4,9
  98. L=K-3
  99. 5 S=S+XL(L,J)*RPL(K,I)
  100. 4 XR(J,I)=S
  101. DO 8 I=1,9
  102. 8 FF(I)=0.D0
  103. RC(1,1)=Y23*Y23
  104. RC(2,1)=Y13*Y13
  105. RC(3,1)=-2.*Y13*Y23
  106. RC(1,2)=X23*X23
  107. RC(2,2)=X13*X13
  108. RC(3,2)=-2.*X13*X23
  109. RC(1,3)=-2.*X23*Y23
  110. RC(2,3)=-2.*X13*Y13
  111. RC(3,3)=2.*(X13*Y23+X23*Y13)
  112. DO 9 I=1,3
  113. DO 9 J=1,9
  114. S=0.D0
  115. DO 10 K=1,3
  116. 10 S=S+RC(K,I)*XR(K,J)
  117. 9 FF(J)=FF(J)+S*TWL(I)
  118. K=0
  119. DO 11 I=1,IELE
  120. IP=6*(I-1)
  121. FT(IP+1)=FU(2*I-1)
  122. FT(IP+2)=FU(2*I)
  123. FT(IP+3)=DT*FF(3*I-2)
  124. FT(IP+4)=DT*FF(3*I)
  125. FT(IP+5)=-DT*FF(3*I-1)
  126. FT(IP+6)=0.
  127. 11 CONTINUE
  128. DO 12 I=1,IELE
  129. KP=6*(I-1)
  130. MP=KP+3
  131. DO 13 J=1,3
  132. DO 14 JP=1,3
  133. FFE(J,I)=FFE(J,I)+P(JP,J)*FT(JP+KP)
  134. FFE(J+3,I)=FFE(J+3,I)+P(JP,J)*FT(JP+MP)
  135. 14 CONTINUE
  136. 13 CONTINUE
  137. 12 CONTINUE
  138. RETURN
  139. END
  140.  
  141.  

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