Télécharger la plaquette Cast3M


fr en

Télécharger coq8ep.eso

Retour à la liste

Numérotation des lignes : 

coq8ep

  1. C COQ8EP    SOURCE    CHAT      05/01/12    22:21:49     5004
  2.       SUBROUTINE COQ8EP(XX,NBNN,NBPGAU,LRE,NSTRS,EPAIST,EXCEN,
  3.      1    DZEGAU,SHPCOQ,SHPELE,XDDL,DEFO,IRR)
  4. C=======================================================================
  5. C
  6. C ENTREES
  7. C     XX(3,NBNN):  TABLEAU DES COORDONNEES DES NOEUDS
  8. C     NBNN    : NOMBRE DE NOEUDS
  9. C     NBPGAU  : NOMBRE DE POINTS DE GAUSS
  10. C     LRE     :    NOMBRE DE DDL
  11. C     NSTRS   :    NOMBRE DE DEFORMATIONS
  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  SHPCOQ(6,NBNN,NBPGAU) :FONCTIONS DE FORME ET DERIVEES AUX PTS DE GSS
  17. C  SHPELE(6,NBNN,NBNN) :FONCTIONS DE FORME ET DERIVEES AUX NOEUDS
  18. C     XDDL(LRE):   DEPLACEMENTS AU NOEUDS
  19. C SORTIES
  20. C  DEFO(5*NBPGAU):5 DEFORMATIONS AUX NBPGAU POINTS DE GAUSS
  21. C     IRR     :    INDICATEUR DE SUCCES ( 1 ) , D ECHEC (0 OU-1)
  22. C-----------------------------------------------------------------------
  23.       IMPLICIT INTEGER(I-N)
  24.       IMPLICIT REAL*8(A-H,O-Z)
  25.       PARAMETER(XZER=0.D0)
  26.       DIMENSION DZEGAU(*),SHPCOQ(6,NBNN,*),SHPELE(6,NBNN,*)
  27.       DIMENSION XX(3,*),DEFO(*),XDDL(*)
  28.       DIMENSION TXR(3,3,8),B(5,48),TH(8),EXC(8),T(9),XSTRS(5)
  29. C
  30.       NDEFO=NSTRS*NBPGAU
  31.       DO 2 J=1,NDEFO
  32.       DEFO(J)=XZER
  33.     2 CONTINUE
  34. C
  35. C   ---   DETERMINATION DES AXES LOCAUX AU NOEUD I   ---
  36. C
  37.       CALL CQ8LOC(XX,NBNN,SHPELE,TXR,IRR)
  38. C
  39.       DO 70 K=1,NBNN
  40.       TH(K)=EPAIST
  41.       EXC(K)=EXCEN
  42.    70 CONTINUE
  43. C
  44. C----------------------------------------------------------------------
  45. C     BOUCLE SUR LES POINTS DE GAUSS
  46. C----------------------------------------------------------------------
  47. C
  48. C
  49.       DO 100 IA=1,NBPGAU
  50.       E3 = DZEGAU(IA)
  51. C-------------------------------
  52. C     CALCUL DE LA MATRICE B
  53. C-------------------------------
  54.       DO 3 I=1,5
  55.       DO 3 J=1,LRE
  56.       B(I,J)=XZER
  57.     3 CONTINUE
  58.       CALL BCOQ8E(IA,XX,NBNN,TH,EXC,B,DET,E3,SHPCOQ,TXR,IRR)
  59. C
  60.       IF(IRR.NE.1) RETURN
  61. C
  62. C     CALCUL DES DEFORMATIONS POUR LE POINT DE GAUSS
  63. C
  64.       CALL BST(B,XDDL,LRE,NSTRS,XSTRS)
  65. C
  66. C     RANGEMENT DANS DEFO(NSTRS*NBPGAU)
  67. C
  68.       DO 80 K=1,5
  69.       JE=5*(IA-1)+K
  70.       DEFO(JE)=XSTRS(K)
  71.    80 CONTINUE
  72. C
  73.   100 CONTINUE
  74.       RETURN
  75.       END
  76.  
  77.  

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