1 : C BBGFDK SOURCE CHAT 05/01/12 21:36:26 5004 2 : SUBROUTINE BBGFDK(QSI,ETA,GEOM,BB) 3 : C|===============================================================| 4 : C| | 5 : C| EVALUATION DE LA MATRICE BB(4,9) AU POINT QSI ETA | 6 : C| POUR L'ELEMENT DKT | 7 : C| | 8 : C| CODE SUO X.Z. MARS 86 | 9 : C|===============================================================| 10 : IMPLICIT INTEGER(I-N) 11 : IMPLICIT REAL*8 (A-H,O-Z) 12 : DIMENSION BXP(9),BYP(9),BXE(9),BYE(9) 13 : DIMENSION GEOM(*),BB(4,*) 14 : C 15 : DATA UN/1.0D0/,DEUX/2.0D0/,QUATRE/4.0D0/,SIX/6.0D0/ 16 : C 17 : CL4=GEOM(1) 18 : CL5=GEOM(2) 19 : CL6=GEOM(3) 20 : SL4=GEOM(4) 21 : SL5=GEOM(5) 22 : SL6=GEOM(6) 23 : D4= GEOM(7) 24 : D5= GEOM(8) 25 : D6= GEOM(9) 26 : SU4=GEOM(10) 27 : SU5=GEOM(11) 28 : SU6=GEOM(12) 29 : VJ11=GEOM(13) 30 : VJ12=GEOM(14) 31 : VJ21=GEOM(15) 32 : VJ22=GEOM(16) 33 : DETJ=GEOM(17) 34 : PS2=UN-DEUX*QSI 35 : ET2=UN-DEUX*ETA 36 : SU5E=SU5*ET2 37 : C---- DERIVEE DE BETAX(I) /QSI 38 : BXP(1)=CL6*PS2+ETA*(CL5-CL6) 39 : BXP(2)=-ETA*D5 40 : BXP(3)=-QUATRE+SIX*(QSI+ETA)-ETA*SU5 41 : BXP(4)=-CL6*PS2+ETA*(CL4+CL6) 42 : BXP(5)=ETA*D4 43 : BXP(6)=-DEUX+SIX*QSI+ETA*SU4 44 : BXP(7)=-ETA*(CL5+CL4) 45 : BXP(8)=ETA*(D4-D5) 46 : BXP(9)=-ETA*(SU5-SU4) 47 : C---- DERIVEE DE BETAY(I) /QSI 48 : BYP(1)=ETA*SL5 49 : BYP(2)=UN-ETA*SU5 50 : BYP(3)=-BXP(2) 51 : BYP(4)=ETA*SL4 52 : BYP(5)=-UN+ETA*SU4 53 : BYP(6)=-BXP(5) 54 : BYP(7)=-ETA*(SL5+SL4) 55 : BYP(8)=BXP(9) 56 : BYP(9)=-BXP(8) 57 : C---- DERIVEE DE BETAX(I) /ETA 58 : BXE(1)=-CL5*ET2-QSI*(CL6-CL5) 59 : BXE(2)=D5*ET2-QSI*D5 60 : BXE(3)=-QUATRE+SIX*(QSI+ETA)+SU5E-QSI*SU5 61 : BXE(4)=QSI*(CL4+CL6) 62 : BXE(5)=QSI*D4 63 : BXE(6)=QSI*SU4 64 : BXE(7)=CL5*ET2-QSI*(CL4+CL5) 65 : BXE(8)=D5*ET2+QSI*(D4-D5) 66 : BXE(9)=-DEUX+SIX*ETA+SU5E+QSI*(SU4-SU5) 67 : C---- DERIVEE DE BETAY(I) /ETA 68 : BYE(1)=-SL5*ET2+QSI*SL5 69 : BYE(2)=UN+SU5E-QSI*SU5 70 : BYE(3)=-BXE(2) 71 : BYE(4)=QSI*SL4 72 : BYE(5)=BXE(6) 73 : BYE(6)=-BXE(5) 74 : BYE(7)=SL5*ET2-QSI*(SL4+SL5) 75 : BYE(8)=-UN+SU5E+QSI*(SU4-SU5) 76 : BYE(9)=-BXE(8) 77 : C---- CALCUL DE BB(4,9) 78 : DO 10 I=1,9 79 : BB(1,I)=(VJ11*BXP(I)+VJ12*BXE(I)) 80 : BB(4,I)=(VJ21*BYP(I)+VJ22*BYE(I)) 81 : BB(2,I)=(VJ21*BXP(I)+VJ22*BXE(I)) 82 : BB(3,I)=(VJ11*BYP(I)+VJ12*BYE(I)) 83 : 10 CONTINUE 84 : RETURN 85 : END 86 :
© Cast3M 2003 - All rights reserved.
Disclaimer