Télécharger cq8gra.eso

Retour à la liste

Numérotation des lignes :

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

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