Télécharger coq3st.eso

Retour à la liste

Numérotation des lignes :

  1. C COQ3ST SOURCE CHAT 05/01/12 22:21:38 5004
  2. SUBROUTINE COQ3ST(XXE,UUE,CONT,DDHOOK)
  3. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  4. C
  5. C ELEMENT COQ3 CALCUL DES EFFORTS AU CENTRE DE GRAVITE.
  6. C
  7. CCCCCCCCCCC ENTREE
  8. C XXE(3,3) COORDONEES DES 3 NOEUDS DE L ELEMENT
  9. C RANGES X(1,2)=COORDONNEE X DU NOEUD 2
  10. C X(2,2)=COORDONNEE Y DU NOEUD 2
  11. C X(3,2)=COORDONNEE Z DU NOEUD 2
  12. C UUE(18) VECTEUR DEPLACEMENT DE L ELEMENT
  13. C
  14. CCCCCCCCCCC SORTIE
  15. C CONT EFFORTS AUX CENTRE DE GRAVITE
  16. C RANGES 3 EFFORTS PUIS
  17. C 3 MOMENTS DE FLEXION
  18. C S1 EST PORTE PAR LE COTE 1 DE L ELEMENT
  19. C S2 EST PORTE PAR LA DIRECTION PERPENDICULAIRE
  20. C A 1 ET DANS LE PLAN DE L ELEMENT
  21. C
  22. CCCCCCCCCCCC TRAVAIL
  23. C
  24. CCCCCCCCCCCCCCCCAUTRES TABLEUX IL FAUT EN TOUT 160
  25. C
  26. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  27. IMPLICIT INTEGER(I-N)
  28. IMPLICIT REAL*8(A-H,O-Z)
  29. DIMENSION XXE(3,*),UUE(*),CONT(*),DDHOOK(6,*)
  30. DIMENSION XGRAV(3),TW(6),CM(3),TM(3),CONT0(6)
  31. DIMENSION UE(18),XX(3),YY(3),P(3,3),XE(3),YE(3),ZE(3),
  32. 1 DT(6),DF(9)
  33. 1,SK(3,6),T(3,3),Q(3,9),SD(3),CI(3),CE(3)
  34. DATA UNTIER,DETIER/.333333333333333D0,.6666666666666667D0/
  35. DATA QUTIER,CISOIX/1.33333333333333D0,.0833333333333333D0/
  36. DATA UNSIX/0.1666666666666667D0/
  37. C
  38. C
  39. C
  40. PS1=0.5
  41. RAI=1.
  42. DO 2 IP=1,18
  43. 2 UE(IP)=UUE(IP)
  44. DO 4 IP=1,3
  45. XE(IP)=XXE(1,IP)
  46. YE(IP)=XXE(2,IP)
  47. ZE(IP)=XXE(3,IP)
  48. 4 CONTINUE
  49. CALL PASSA(XE,YE,ZE, P,X13,X23,Y13,Y23)
  50. XGRAV(1)=(XE(1)+XE(2)+XE(3))*UNTIER
  51. XGRAV(2)=(YE(1)+YE(2)+YE(3))*UNTIER
  52. XGRAV(3)=(ZE(1)+ZE(2)+ZE(3))*UNTIER
  53. DO 10 IP=1,6
  54. IP1=3*(IP-1)
  55. DO 5 IQ=1,3
  56. 5 XX(IQ)=UE(IQ+IP1)
  57. CALL MULMAT(YY,P,XX,3,1,3)
  58. DO 7 IQ=1,3
  59. 7 UE(IQ+IP1)=YY(IQ)
  60. 10 CONTINUE
  61. DO 8 IP=1,3
  62. IP1=IP-1
  63. IP2=2*IP1
  64. IP3=3*IP1
  65. IP6=6*IP1
  66. DO 9 IQ=1,2
  67. 9 DT(IQ+IP2)=UE(IQ+IP6)
  68. DO 11 IQ=3,5
  69. IQ2=IQ-2
  70. 11 DF(IQ2+IP3)=UE(IQ+IP6)
  71. 8 CONTINUE
  72. S=1./(X13*Y23-X23*Y13)
  73. A1=Y13
  74. A2=-X13
  75. A3=Y23
  76. A4=-X23
  77. SK(1,1)=Y23
  78. SK(1,2)=0.
  79. SK(1,3)=-Y13
  80. SK(1,4)=0.
  81. SK(1,5)=Y13-Y23
  82. SK(1,6)=0.
  83. SK(2,1)=0.
  84. SK(2,2)=-X23
  85. SK(2,3)=0.
  86. SK(2,4)=X13
  87. SK(2,5)=0.
  88. SK(2,6)=X23-X13
  89. SK(3,1)=-X23
  90. SK(3,3)=X13
  91. SK(3,5)=X23-X13
  92. SK(3,2)=Y23
  93. SK(3,4)=-Y13
  94. SK(3,6)=Y13-Y23
  95. T(1,1)=Y23*Y23
  96. T(1,2)=Y13*Y13
  97. T(1,3)=-2.*(Y13*Y23)
  98. T(2,1)=X23*X23
  99. T(2,2)=X13*X13
  100. T(2,3)=-2.*(X13*X23)
  101. T(3,1)=-2.*X23*Y23
  102. T(3,2)=-2.*X13*Y13
  103. T(3,3)=2.*(X13*Y23+X23*Y13)
  104. Q(1,1)=DETIER
  105. Q(1,2)=UNTIER*(2.*A1+A3)
  106. Q(1,3)=UNTIER*(2.*A2+A4)
  107. Q(1,4)=-QUTIER
  108. Q(1,5)=UNTIER*(-A1+2.*A3)
  109. Q(1,6)=UNTIER*(-A2+2.*A4)
  110. Q(1,7)=DETIER
  111. Q(1,8)=-A1+UNTIER*A3
  112. Q(1,9)=-A2+UNTIER*A4
  113. Q(2,1)=-QUTIER
  114. Q(2,2)=UNTIER*(2.*A1-A3)
  115. Q(2,3)=UNTIER*(2.*A2-A4)
  116. Q(2,4)=DETIER
  117. Q(2,5)=UNTIER*(A1+2.*A3)
  118. Q(2,6)=UNTIER*(A2+2.*A4)
  119. Q(2,7)=DETIER
  120. Q(2,8)=-A3+UNTIER*A1
  121. Q(2,9)=-A4+UNTIER*A2
  122. Q(3,1)=-DETIER
  123. Q(3,2)=UNTIER*A1+0.5*A3
  124. Q(3,3)=UNTIER*A2+0.5*A4
  125. Q(3,4)=-DETIER
  126. Q(3,5)=0.5*A1+UNTIER*A3
  127. Q(3,6)=0.5*A2+UNTIER*A4
  128. Q(3,7)=QUTIER
  129. Q(3,8)=-UNSIX*(A1+A3)
  130. Q(3,9)=-UNSIX*(A2+A4)
  131. CALL MULMAT(CM,SK,DT,3,1,6)
  132. CALL MULMAT(SD,Q,DF,3,1,9)
  133. CALL MULMAT(TM,T,SD,3,1,3)
  134. CPLA=S*S
  135. CMEM=RAI*S
  136. DO 306 IP=1,3
  137. CONT0(IP+3)=-TM(IP)*CPLA
  138. CONT0(IP )= CM(IP)*CMEM
  139. 306 CONTINUE
  140. CALL MULMAT(CONT,DDHOOK,CONT0,6,1,6)
  141. RETURN
  142. END
  143.  
  144.  

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