1 : C BLMAST SOURCE CHAT 05/01/12 21:40:59 5004 2 : SUBROUTINE BLMAST(IGAU,MFR,NBNO,LRE,IFOU,NDDL,NN,XEL, 3 : 1 SHPTOT,SHP,BLX,BLY,BLZ,BLT,DJAC) 4 : C======================================================================= 5 : C 6 : C CALCULE LES MATRICES BL POUR LES ELEMENTS LIQUIDES 7 : C 8 : C ROUTINE FORTRAN PUR 9 : C======================================================================= 10 : C INPUT 11 : C IGAU=NUMERO DU POINT DE GAUSS 12 : C MFR =NUMERO DE LA FORMULATION 13 : C NBNO=NOMBRE DE NOEUDS 14 : C LRE =NOMBRE DE COLONNES DE LA MATRICE B 15 : C IFOU=IFOUR DE CCOPTIO 16 : C NDDL=NOMBRE DE DEGRES DE LIBERTE 17 : C NN =NUMERO DU MODE DE FOURIER 18 : C XEL =COORDONNEES LOCALES DE L ELEMENT 19 : C SHPTOT(6,NBNO,NBGAU)=FONCTIONS DE FORMES ET DERIVEES 20 : C ZONE DE TRAVAIL 21 : C SHP(6,NBNO)=TABLEAU DE TRAVAIL 22 : C OUTPUT 23 : C DJAC=JACOBIEN 24 : C BLX = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / X 25 : C BLY = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / Y 26 : C BLZ = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / Z 27 : C BLT = MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME / TETA 28 : C======================================================================= 29 : IMPLICIT INTEGER(I-N) 30 : IMPLICIT REAL*8(A-H,O-Z) 31 : DIMENSION XEL(3,1),SHP(6,1),SHPTOT(6,NBNO,1) 32 : DIMENSION BLX(NDDL,1),BLY(NDDL,1),BLZ(NDDL,1),BLT(NDDL,1) 32 : CALL ZERO(BLX,NDDL,LRE) 32 : CALL ZERO(BLY,NDDL,LRE) 32 : CALL ZERO(BLZ,NDDL,LRE) 32 : CALL ZERO(BLT,NDDL,LRE) 33 : IFR=IFOU+4 34 : DIM3=1.D0 35 : RR=0.D0 36 : GOTO (666,10,10,20,30,40) ,IFR 37 : C 38 : C ELEMENTS LIQUIDES BIDIM CONTRAINTES PLANES OU DEFRMTNS PLANES 39 : C 40 : 10 CONTINUE 41 : DO 101 NP=1,NBNO 42 : SHP(2,NP)=SHPTOT(2,NP,IGAU) 43 : SHP(3,NP)=SHPTOT(3,NP,IGAU) 44 : 101 CONTINUE 44 : CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,2,DIM3,RR,DJAC) 45 : K=1 46 : DO 102 NP=1,NBNO 47 : BLX(2,K+1)=SHP(2,NP) 48 : BLY(2,K+1)=SHP(3,NP) 49 : 102 K=K+2 50 : GOTO 666 51 : C 52 : C ELEMENTS LIQUIDES BIDIM AXISYMETRIQUE 53 : C 54 : 20 CONTINUE 55 : DO 201 NP=1,NBNO 56 : SHP(1,NP)=SHPTOT(1,NP,IGAU) 57 : SHP(2,NP)=SHPTOT(2,NP,IGAU) 58 : SHP(3,NP)=SHPTOT(3,NP,IGAU) 59 : 201 CONTINUE 59 : CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,2,DIM3,RR,DJAC) 60 : K=1 61 : DO 202 NP=1,NBNO 62 : BLX(2,K+1)=SHP(2,NP) 63 : BLY(2,K+1)=SHP(3,NP) 64 : 202 K=K+2 65 : GOTO 666 66 : C 67 : C ELEMENTS LIQUIDES BIDIM FOURIER 68 : C 69 : 30 CONTINUE 70 : DO 301 NP=1,NBNO 71 : SHP(1,NP)=SHPTOT(1,NP,IGAU) 72 : SHP(2,NP)=SHPTOT(2,NP,IGAU) 73 : SHP(3,NP)=SHPTOT(3,NP,IGAU) 74 : 301 CONTINUE 74 : CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,2,DIM3,RR,DJAC) 75 : K=1 76 : DO 302 NP=1,NBNO 77 : BLX(2,K+1)= SHP(2,NP) 78 : BLY(2,K+1)= SHP(3,NP) 79 : BLT(2,K+1)= -SHP(1,NP)*NN/RR 80 : 302 K=K+2 81 : GOTO 666 82 : C 83 : C ELEMENTS LIQUIDES TRIDIM 84 : C 85 : 40 CONTINUE 86 : DO 401 NP=1,NBNO 87 : SHP(2,NP)=SHPTOT(2,NP,IGAU) 88 : SHP(3,NP)=SHPTOT(3,NP,IGAU) 89 : SHP(4,NP)=SHPTOT(4,NP,IGAU) 90 : 401 CONTINUE 90 : CALL DEVOLU(XEL,SHP,MFR,NBNO,IFOU,NN,3,DIM3,RR,DJAC) 91 : K=1 92 : DO 402 NP=1,NBNO 93 : BLX(2,K+1)=SHP(2,NP) 94 : BLY(2,K+1)=SHP(3,NP) 95 : BLZ(2,K+1)=SHP(4,NP) 96 : 402 K=K+2 97 : GOTO 666 98 : C 99 : 666 CONTINUE 100 : RETURN 101 : END 102 : 103 : 104 :
© Cast3M 2003 - All rights reserved.
Disclaimer