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 36 : CALL CQ8LOC(XX,NBNN,SHPELE,TXR,IRR) 37 : C 38 : DO 70 K=1,NBNN 39 : TH(K)=EPAIST 40 : EXC(K)=EXCEN 41 : 70 CONTINUE 42 : C 43 : C---------------------------------------------------------------------- 44 : C BOUCLE SUR LES POINTS DE GAUSS 45 : C---------------------------------------------------------------------- 46 : C 47 : C 48 : DO 100 IA=1,NBPGAU 49 : E3 = DZEGAU(IA) 50 : C------------------------------- 51 : C CALCUL DE LA MATRICE B 52 : C------------------------------- 53 : DO 3 I=1,5 54 : DO 3 J=1,LRE 55 : B(I,J)=XZER 56 : 3 CONTINUE 56 : CALL BCOQ8E(IA,XX,NBNN,TH,EXC,B,DET,E3,SHPCOQ,TXR,IRR) 57 : C 58 : IF(IRR.NE.1) RETURN 59 : C 60 : C CALCUL DES DEFORMATIONS POUR LE POINT DE GAUSS 61 : C 61 : CALL BST(B,XDDL,LRE,NSTRS,XSTRS) 62 : C 63 : C RANGEMENT DANS DEFO(NSTRS*NBPGAU) 64 : C 65 : DO 80 K=1,5 66 : JE=5*(IA-1)+K 67 : DEFO(JE)=XSTRS(K) 68 : 80 CONTINUE 69 : C 70 : 100 CONTINUE 71 : RETURN 72 : END 73 :
© Cast3M 2003 - All rights reserved.
Disclaimer