Télécharger coq8ep.eso

Retour à la liste

Numérotation des lignes :

  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