Télécharger cq8lc1.eso

Retour à la liste

Numérotation des lignes :

  1. C CQ8LC1 SOURCE CHAT 05/01/12 22:27:19 5004
  2. SUBROUTINE CQ8LC1(NOBG,NBNN,E,XX,TH,SHPCOQ,TXR,T,IRR)
  3. C=======================================================================
  4. C CALCULE REPERE LOCAL EN UN POINT DE GAUSS POUR COQ8 ET COQ6
  5. C ENTREES
  6. C -------
  7. C NOBG : NUMERO DU POINT DE GAUSS
  8. C NBNN : NOMBRE DE NOEUDS
  9. C E : COORDONNEES REDUITES DU PT DE GAUSS DS L EPAI
  10. C XX(3,NBNN) : TABLEAU DES COORDONNEES DES NOEUDS
  11. C TH(NBNN) : LES EPAISSEURS NODALES
  12. C SHPCOQ(6,NBNN,NBPGAU): VALEURS DES FONCTIONS DE FORME ET
  13. C DE LEURS DERIVEES AUX POINTS DE GAUSS
  14. C TXR(3,3,NBNN) :NBNN MATRICES 3 3 CONTENANT REPERES NOEUDS
  15. C SORTIES
  16. C -------
  17. C T(9) : TABLEAU CONTENANT LE REPERE LOCAL
  18. C IRR : INDICATEUR D ERREUR
  19. C-----------------------------------------------------------------------
  20. C EBERSOLT RECUPERATION BILBO JUILLET 86
  21. C=======================================================================
  22. IMPLICIT INTEGER(I-N)
  23. IMPLICIT REAL*8 (A-H,O-Z)
  24. DIMENSION XX(3,*),T(*),SHPCOQ(6,NBNN,*),TXR(3,3,*),TH(*)
  25. DIMENSION XJ(3,3)
  26. C
  27. DATA XZER/0.D0/
  28. C
  29. C DETERMINATION DU JACOBIEN ET DE SON DETERMINANT AU POINT (R,S,T)
  30. C
  31. CALL COQ8JC(NOBG,NBNN,E,XX,TH,TXR,SHPCOQ,XJ,DET,IRR)
  32. C
  33. IF(IRR.EQ.-1) RETURN
  34. C*
  35. C* DETERMINATION DES COSINUS DIRECTEURS DES AXES LOCAUX EN CE POINT
  36. C*
  37. DO 10 I=1,3
  38. DO 10 J=1,2
  39. K=3*(J-1)+I
  40. T(K)= XJ(J,I)
  41. 10 CONTINUE
  42. C*
  43. C* PRODUITS VECTORIELS ET NORMALISATIONS
  44. C*
  45. CALL CROSS2(T(1),T(4),T(7),IRR)
  46. CALL CROSS2(T(7),T(1),T(4),IRR)
  47. CALL CROSS2(T(4),T(7),T(1),IRR)
  48. C
  49. RETURN
  50. END
  51.  
  52.  

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