Télécharger la plaquette Cast3M


fr en

Télécharger cq8bse.eso

Retour à la liste

Numérotation des lignes : 

cq8bse

  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