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 30 : CALL COQ8JC(NOBG,NBNN,E,XX,TH,TXR,SHPCOQ,XJ,DET,IRR) 31 : C 32 : IF(IRR.EQ.-1) RETURN 33 : C* 34 : C* DETERMINATION DES COSINUS DIRECTEURS DES AXES LOCAUX EN CE POINT 35 : C* 36 : DO 10 I=1,3 37 : DO 10 J=1,2 38 : K=3*(J-1)+I 39 : T(K)= XJ(J,I) 40 : 10 CONTINUE 41 : C* 42 : C* PRODUITS VECTORIELS ET NORMALISATIONS 43 : C* 43 : CALL CROSS2(T(1),T(4),T(7),IRR) 43 : CALL CROSS2(T(7),T(1),T(4),IRR) 43 : CALL CROSS2(T(4),T(7),T(1),IRR) 44 : C 45 : RETURN 46 : END 47 :
© Cast3M 2003 - All rights reserved.
Disclaimer