1 : C BTUYO SOURCE CHAT 05/01/12 21:44:30 5004 2 : SUBROUTINE BTUYO(IGAU,MINTE,WRK1,WRK2,WRK3,XL,DJAC,KERRE) 3 : *----------------------------------------------------------------------- 4 : * 5 : * CALCUL DE B POUR L'ELEMENT TUYO 6 : * 7 : *----------------------------------------------------------------------- 8 : IMPLICIT INTEGER(I-N) 9 : IMPLICIT REAL*8 (A-H,O-Z) 10 : -INC CCOPTIO 11 : -INC SMINTE 12 : * 13 : SEGMENT WRK1 14 : REAL*8 DDHOOK(NSTRS,NSTRS) ,DDHOMU(NSTRS,NSTRS) 15 : REAL*8 REL(LRE,LRE) , XE(3,NBBB) 16 : ENDSEGMENT 17 : * 18 : SEGMENT WRK2 19 : REAL*8 SHPWRK(6,NBNO) ,BGENE(NSTRS,LRE) 20 : ENDSEGMENT 21 : * 22 : SEGMENT WRK3 23 : REAL*8 WORK(LW) 24 : ENDSEGMENT 25 : * 26 : * INITIALISATIONS 27 : * 28 : KERRE=0 29 : NSTRS=DDHOOK(/1) 30 : LRE =REL(/1) 30 : CALL ZERO(BGENE,NSTRS,LRE) 31 : SL=1.D0/XL 32 : SL2 = 2.D0*SL 33 : SL4 = SL2*SL2 34 : * 35 : T=ETAGAU(IGAU) 36 : ST=SIN(T) 37 : CT=COS(T) 38 : S2T=SIN(2.D0*T) 39 : C2T=COS(2.D0*T) 40 : S4T=SIN(4.D0*T) 41 : C4T=COS(4.D0*T) 42 : * 43 : EP=WORK(1) 44 : R =WORK(2)-EP*0.5D0 45 : R2=R*R 46 : GR=WORK(3) 47 : UNSUR=0.D0 48 : IF(GR.NE.0.D0) UNSUR=1.D0/GR 49 : DJAC=XL*0.5D0*R 50 : * 51 : * CERTAINS TERMES SONT MIS EN SOMMEIL 52 : * 53 : * LES MEMBRANES 54 : * 55 : * EPS TET-TET 56 : * 57 : BGENE(1, 7)=SHPTOT(1,1,IGAU)/R 58 : BGENE(1,22)=SHPTOT(1,2,IGAU)/R 59 : * 60 : * EPS PHI-PHI 61 : * 62 : BGENE(2, 1)=SHPTOT(2,1,IGAU)*SL2 63 : BGENE(2,16)=SHPTOT(2,2,IGAU)*SL2 64 : * 65 : BGENE(2, 3)=-SHPTOT(5,3,IGAU)*SL4*R*CT 66 : BGENE(2, 5)= SHPTOT(5,4,IGAU)*SL4*R*CT*XL 67 : BGENE(2,18)=-SHPTOT(5,5,IGAU)*SL4*R*CT 68 : BGENE(2,20)= SHPTOT(5,6,IGAU)*SL4*R*CT*XL 69 : * 70 : BGENE(2, 2)=-SHPTOT(5,3,IGAU)*SL4*R*ST 71 : BGENE(2, 6)=-SHPTOT(5,4,IGAU)*SL4*R*ST*XL 72 : BGENE(2,17)=-SHPTOT(5,5,IGAU)*SL4*R*ST 73 : BGENE(2,21)=-SHPTOT(5,6,IGAU)*SL4*R*ST*XL 74 : * 75 : BGENE(2, 8)=SHPTOT(2,1,IGAU)*SL2*C2T 76 : BGENE(2,23)=SHPTOT(2,2,IGAU)*SL2*C2T 77 : BGENE(2, 9)=SHPTOT(2,1,IGAU)*SL2*C4T 78 : BGENE(2,24)=SHPTOT(2,2,IGAU)*SL2*C4T 79 : BGENE(2,10)=SHPTOT(2,1,IGAU)*SL2*S2T 80 : BGENE(2,25)=SHPTOT(2,2,IGAU)*SL2*S2T 81 : BGENE(2,11)=SHPTOT(2,1,IGAU)*SL2*S4T 82 : BGENE(2,26)=SHPTOT(2,2,IGAU)*SL2*S4T 83 : * 84 : * BGENE(2, 3)=BGENE(2, 3)+SHPTOT(1,3,IGAU)*UNSUR 85 : * BGENE(2, 5)=BGENE(2, 5)-SHPTOT(1,4,IGAU)*UNSUR*XL 86 : * BGENE(2,18)=BGENE(2,18)+SHPTOT(1,5,IGAU)*UNSUR 87 : * BGENE(2,20)=BGENE(2,20)-SHPTOT(1,6,IGAU)*UNSUR*XL 88 : * 89 : BGENE(2,7) =SHPTOT(1,1,IGAU)*CT*UNSUR 90 : BGENE(2,22)=SHPTOT(1,2,IGAU)*CT*UNSUR 91 : * 92 : BGENE(2,12)=SHPTOT(1,1,IGAU)*UNSUR*(CT*C2T+ST*S2T/2.D0) 93 : BGENE(2,27)=SHPTOT(1,2,IGAU)*UNSUR*(CT*C2T+ST*S2T/2.D0) 94 : BGENE(2,13)=SHPTOT(1,1,IGAU)*UNSUR*(CT*C4T+ST*S4T/4.D0) 95 : BGENE(2,28)=SHPTOT(1,2,IGAU)*UNSUR*(CT*C4T+ST*S4T/4.D0) 96 : BGENE(2,14)=SHPTOT(1,1,IGAU)*UNSUR*(CT*S2T-ST*C2T/2.D0) 97 : BGENE(2,29)=SHPTOT(1,2,IGAU)*UNSUR*(CT*S2T-ST*C2T/2.D0) 98 : BGENE(2,15)=SHPTOT(1,1,IGAU)*UNSUR*(CT*S4T-ST*C4T/4.D0) 99 : BGENE(2,30)=SHPTOT(1,2,IGAU)*UNSUR*(CT*S4T-ST*C4T/4.D0) 100 : * 101 : * BGENE(2, 4)=SHPTOT(1,1,IGAU)*R*ST*UNSUR 102 : * BGENE(2,19)=SHPTOT(1,2,IGAU)*R*ST*UNSUR 103 : * 104 : * GAM TET-PHI 105 : * 106 : BGENE(3, 8)=SHPTOT(1,1,IGAU)*(ST*C2T*UNSUR-2.D0*S2T/R) 107 : BGENE(3,23)=SHPTOT(1,2,IGAU)*(ST*C2T*UNSUR-2.D0*S2T/R) 108 : BGENE(3, 9)=SHPTOT(1,1,IGAU)*(ST*C4T*UNSUR-4.D0*S4T/R) 109 : BGENE(3,24)=SHPTOT(1,2,IGAU)*(ST*C4T*UNSUR-4.D0*S4T/R) 110 : BGENE(3,10)=SHPTOT(1,1,IGAU)*(ST*S2T*UNSUR+2.D0*C2T/R) 111 : BGENE(3,25)=SHPTOT(1,2,IGAU)*(ST*S2T*UNSUR+2.D0*C2T/R) 112 : BGENE(3,11)=SHPTOT(1,1,IGAU)*(ST*S4T*UNSUR+4.D0*C4T/R) 113 : BGENE(3,26)=SHPTOT(1,2,IGAU)*(ST*S4T*UNSUR+4.D0*C4T/R) 114 : * 115 : BGENE(3, 4)=-SHPTOT(2,1,IGAU)*SL2*R 116 : BGENE(3,19)=-SHPTOT(2,2,IGAU)*SL2*R 117 : * 118 : * BGENE(3, 1)=SHPTOT(1,1,IGAU)*ST*UNSUR 119 : * BGENE(3,16)=SHPTOT(1,2,IGAU)*ST*UNSUR 120 : * 121 : * BGENE(3, 3)=-SHPTOT(2,3,IGAU)*SL2*R*UNSUR*ST*CT 122 : * BGENE(3, 5)= SHPTOT(2,4,IGAU)*SL2*R*UNSUR*ST*CT*XL 123 : * BGENE(3,18)=-SHPTOT(2,5,IGAU)*SL2*R*UNSUR*ST*CT 124 : * BGENE(3,20)= SHPTOT(2,6,IGAU)*SL2*R*UNSUR*ST*CT*XL 125 : * 126 : * BGENE(3, 2)=-SHPTOT(2,3,IGAU)*SL2*R*UNSUR*ST*ST 127 : * BGENE(3, 6)=-SHPTOT(2,4,IGAU)*SL2*R*UNSUR*ST*ST*XL 128 : * BGENE(3,17)=-SHPTOT(2,5,IGAU)*SL2*R*UNSUR*ST*ST 129 : * BGENE(3,21)=-SHPTOT(2,6,IGAU)*SL2*R*UNSUR*ST*ST*XL 130 : * 131 : BGENE(3,12)=-SHPTOT(2,1,IGAU)*SL2*S2T/2.D0 132 : BGENE(3,27)=-SHPTOT(2,2,IGAU)*SL2*S2T/2.D0 133 : BGENE(3,13)=-SHPTOT(2,1,IGAU)*SL2*S4T/4.D0 134 : BGENE(3,28)=-SHPTOT(2,2,IGAU)*SL2*S4T/4.D0 135 : BGENE(3,14)= SHPTOT(2,1,IGAU)*SL2*C2T/2.D0 136 : BGENE(3,29)= SHPTOT(2,2,IGAU)*SL2*C2T/2.D0 137 : BGENE(3,15)= SHPTOT(2,1,IGAU)*SL2*C4T/4.D0 138 : BGENE(3,30)= SHPTOT(2,2,IGAU)*SL2*C4T/4.D0 139 : * 140 : * LES FLEXIONS 141 : * 142 : * QUI TET-TET 143 : * 144 : BGENE(4,12)=SHPTOT(1,1,IGAU)*3.D0*C2T/R2 145 : BGENE(4,27)=SHPTOT(1,2,IGAU)*3.D0*C2T/R2 146 : BGENE(4,13)=SHPTOT(1,1,IGAU)*15.D0*C4T/R2 147 : BGENE(4,28)=SHPTOT(1,2,IGAU)*15.D0*C4T/R2 148 : BGENE(4,14)=SHPTOT(1,1,IGAU)*3.D0*S2T/R2 149 : BGENE(4,29)=SHPTOT(1,2,IGAU)*3.D0*S2T/R2 150 : BGENE(4,15)=SHPTOT(1,1,IGAU)*15.D0*S4T/R2 151 : BGENE(4,30)=SHPTOT(1,2,IGAU)*15.D0*S4T/R2 152 : * 153 : * 0. POUR LES 2 AUTRES 154 : * 155 : RETURN 156 : END 157 :
© Cast3M 2003 - All rights reserved.
Disclaimer