1 : C CQ8BSE SOURCE CHAT 05/01/12 22:27:08 5004 2 : SUBROUTINE CQ8BSE(XX,NBNN,NBPGAU,LRE,EPAIST,EXCEN,DZEGAU, 3 : + POIGAU,SHPCOQ,SHPELE,WORK,XFORC,IRR) 4 : C======================================================================= 5 : C 6 : C (ROUTINE EQUIVALENTE A COQ8BS AVEC EXCENTREMENT.) 7 : C ENTREES 8 : C XX(3,NBNN): TABLEAU DES COORDONNEES DES NOEUDS 9 : C NBNN : NOMBRE DE NOEUDS 10 : C NBPGAU : NOMBRE DEPOINTS DE GAUSS 11 : C LRE : NOMBRE DE DDL 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 POIGAU(NBPGAU): POIDS DE GAUSS 17 : C SHPCOQ(6,NBNN,NBPGAU) : FONCTIONS DE FORME ET DERIVEES AUX PTS DE GSS 18 : C SHPELE(6,NBNN,NBNN) : FONCTIONS DE FORME ET DERIVEES AUX NOEUDS 19 : C WORK(5*NBPGAU): 5XNBPGAU CONTRAINTES AUX POINTS DE GAUSS 20 : C SORTIES 21 : C XFORC(LRE) : 6 FORCES AUX NBNN NOEUDS 22 : C 23 : C IRR : INDICATEUR DE SUCCES ( 1 ) , D ECHEC (0 OU-1) 24 : C----------------------------------------------------------------------- 25 : IMPLICIT INTEGER(I-N) 26 : IMPLICIT REAL*8(A-H,O-Z) 27 : PARAMETER(XZER=0.D0,UN=1.D0) 28 : PARAMETER(X577=.577350269189626D0,UNSIX=.166666666666666D0) 29 : DIMENSION XX(3,*),DZEGAU(*),SHPCOQ(6,8,*),SHPELE(6,8,*),WORK(*) 30 : DIMENSION POIGAU(*) 31 : DIMENSION TXR(3,3,8),B(5,48),TH(8),EXC(8),T(9),XFORC(48),XSTRS(5) 32 : C 33 : USX577=UNSIX/X577 34 : C 35 : DO 2 J=1,LRE 36 : XFORC(J)=XZER 37 : 2 CONTINUE 38 : C 39 : C --- DETERMINATION DES AXES LOCAUX AU NOEUD I --- 40 : C 40 : CALL CQ8LOC(XX,NBNN,SHPELE,TXR,IRR) 41 : C 42 : DO 70 K=1,NBNN 43 : TH(K)=EPAIST 44 : EXC(K)=EXCEN 45 : 70 CONTINUE 46 : C---------------------------------------------------------------------- 47 : C BOUCLE SUR LES POINTS DE GAUSS 48 : C---------------------------------------------------------------------- 49 : DO 100 IA=1,NBPGAU 50 : E3 = DZEGAU(IA) 51 : C------------------------------- 52 : C CALCUL DE LA MATRICE B ET RECUPERATION DES CONTRAITES DU PT DE GAUSS 53 : C------------------------------- 54 : IDEB=(IA-1)*5 55 : DO 3 I=1,5 56 : XSTRS(I)=WORK(IDEB+I) 57 : DO 3 J=1,LRE 58 : B(I,J)=XZER 59 : 3 CONTINUE 59 : CALL BCOQ8E(IA,XX,NBNN,TH,EXC,B,DET,E3,SHPCOQ,TXR,IRR) 60 : C 61 : IF(IRR.NE.1) RETURN 62 : C 63 : C CALCUL DE BSIGMA POUR LE POINT DE GAUSS 64 : C 65 : DJAC=DET*POIGAU(IA) 65 : CALL BSIG8(B,NBNN,XSTRS,DJAC,XFORC) 66 : 100 CONTINUE 67 : C 68 : RETURN 69 : END 70 :
© Cast3M 2003 - All rights reserved.
Disclaimer