Télécharger cq8bse.eso

Retour à la liste

Numérotation des lignes :

  1. C CQ8BSE SOURCE CHAT 05/01/12 22:27:08 5004
  2. SUBROUTINE CQ8BSE(XX,NBNN,NBPGAU,LRE,EPAIST,EXCEN,DZEGAU,
  3. + POIGAU,SHPCOQ,SHPELE,WORK,XFORC,IRR)
  4. C=======================================================================
  5. C
  6. C (ROUTINE EQUIVALENTE A COQ8BS AVEC EXCENTREMENT.)
  7. C ENTREES
  8. C XX(3,NBNN): TABLEAU DES COORDONNEES DES NOEUDS
  9. C NBNN : NOMBRE DE NOEUDS
  10. C NBPGAU : NOMBRE DEPOINTS DE GAUSS
  11. C LRE : NOMBRE DE DDL
  12. C EPAIST : EPAISSEUR DE LA COQUE
  13. C EXCEN : EXCENTREMENT DE LA COQUE
  14. C DZEGAU(NBPGAU): COORDONNEES REDUITES DES POINTS DE GAUSS DANS
  15. C L EPAISSEUR
  16. C POIGAU(NBPGAU): POIDS DE GAUSS
  17. C SHPCOQ(6,NBNN,NBPGAU) : FONCTIONS DE FORME ET DERIVEES AUX PTS DE GSS
  18. C SHPELE(6,NBNN,NBNN) : FONCTIONS DE FORME ET DERIVEES AUX NOEUDS
  19. C WORK(5*NBPGAU): 5XNBPGAU CONTRAINTES AUX POINTS DE GAUSS
  20. C SORTIES
  21. C XFORC(LRE) : 6 FORCES AUX NBNN NOEUDS
  22. C
  23. C IRR : INDICATEUR DE SUCCES ( 1 ) , D ECHEC (0 OU-1)
  24. C-----------------------------------------------------------------------
  25. IMPLICIT INTEGER(I-N)
  26. IMPLICIT REAL*8(A-H,O-Z)
  27. PARAMETER(XZER=0.D0,UN=1.D0)
  28. PARAMETER(X577=.577350269189626D0,UNSIX=.166666666666666D0)
  29. DIMENSION XX(3,*),DZEGAU(*),SHPCOQ(6,8,*),SHPELE(6,8,*),WORK(*)
  30. DIMENSION POIGAU(*)
  31. DIMENSION TXR(3,3,8),B(5,48),TH(8),EXC(8),T(9),XFORC(48),XSTRS(5)
  32. C
  33. USX577=UNSIX/X577
  34. C
  35. DO 2 J=1,LRE
  36. XFORC(J)=XZER
  37. 2 CONTINUE
  38. C
  39. C --- DETERMINATION DES AXES LOCAUX AU NOEUD I ---
  40. C
  41. CALL CQ8LOC(XX,NBNN,SHPELE,TXR,IRR)
  42. C
  43. DO 70 K=1,NBNN
  44. TH(K)=EPAIST
  45. EXC(K)=EXCEN
  46. 70 CONTINUE
  47. C----------------------------------------------------------------------
  48. C BOUCLE SUR LES POINTS DE GAUSS
  49. C----------------------------------------------------------------------
  50. DO 100 IA=1,NBPGAU
  51. E3 = DZEGAU(IA)
  52. C-------------------------------
  53. C CALCUL DE LA MATRICE B ET RECUPERATION DES CONTRAITES DU PT DE GAUSS
  54. C-------------------------------
  55. IDEB=(IA-1)*5
  56. DO 3 I=1,5
  57. XSTRS(I)=WORK(IDEB+I)
  58. DO 3 J=1,LRE
  59. B(I,J)=XZER
  60. 3 CONTINUE
  61. CALL BCOQ8E(IA,XX,NBNN,TH,EXC,B,DET,E3,SHPCOQ,TXR,IRR)
  62. C
  63. IF(IRR.NE.1) RETURN
  64. C
  65. C CALCUL DE BSIGMA POUR LE POINT DE GAUSS
  66. C
  67. DJAC=DET*POIGAU(IA)
  68. CALL BSIG8(B,NBNN,XSTRS,DJAC,XFORC)
  69. 100 CONTINUE
  70. C
  71. RETURN
  72. END
  73.  
  74.  

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