C JBGP3D SOURCE FD218221 24/02/07 21:15:16 11834 subroutine jbgp3d(NBRINC,NINC,NDIMG,NDIMPL,FRAC,GCH,Mbg,bg, # Jsgep,AFFICHE,ERR1) c calcul de Jsgep(i,j)=d (-bgpg(i-1)) / d eps_plastique(j) dans l espace c de minimisation de R implicit real*8 (a-h,o-z) implicit integer (i-n) c dimension de la base des variables generalisees integer NDIMG,ERR1 c nbr max inclusion, nbr reel inclusion, nbr de type de criteres de plast integer NBRINC,NINC,NBTYPL c fractions volumiques des phases real*8 FRAC(0:NBRINC) c indicateur de diffusion produits chimiques dans fissure real*8 GCH(0:NBRINC) c module de Biot du gel real*8 Mbg(0:NBRINC,0:1) c coeff de biot des gels real*8 bg(0:NBRINC,0:1) c Jacobienne des pressions / deformation plastiques real*8 Jsgep(NDIMG,NDIMPL) c fraction inclusion real*8 f logical affiche,affiche_local if(NBRINC.ne.1) then print*,'Nbr inc inattendu ds jbgp3d' err1=1 return end if f=frac(1) t1 = Mbg(1,1) * bg(1,1) t2 = t1 * Gch(1) t3 = 0.1D1 - f t5 = t3 * Gch(1) * t1 t6 = Mbg(0,1) * bg(0,1) t8 = 0.1D1 - 0.3D1 / 0.2D1 * f t10 = t8 * Gch(0) * t6 t12 = 0.1D1 / t3 t14 = t12 * f * Gch(0) * t6 t15 = bg(0,1) ** 2 t16 = Gch(0) * t15 t17 = t8 * t16 t19 = t12 * f * t16 Jsgep(1,1) = t2 Jsgep(1,2) = t2 Jsgep(1,3) = t2 Jsgep(1,4) = 0.0D0 Jsgep(1,5) = 0.0D0 Jsgep(1,6) = 0.0D0 Jsgep(1,7) = t5 Jsgep(1,8) = t5 Jsgep(1,9) = t5 Jsgep(1,10) = 0.0D0 Jsgep(1,11) = 0.0D0 Jsgep(1,12) = 0.0D0 Jsgep(1,13) = 0.0D0 Jsgep(1,14) = 0.0D0 Jsgep(1,15) = 0.0D0 Jsgep(1,16) = 0.0D0 Jsgep(1,17) = 0.0D0 Jsgep(1,18) = 0.0D0 Jsgep(1,19) = 0.0D0 Jsgep(1,20) = 0.0D0 Jsgep(1,21) = 0.0D0 Jsgep(1,22) = 0.0D0 Jsgep(1,23) = 0.0D0 Jsgep(1,24) = 0.0D0 Jsgep(1,25) = 0.0D0 Jsgep(1,26) = 0.0D0 Jsgep(1,27) = 0.0D0 Jsgep(1,28) = 0.0D0 Jsgep(1,29) = 0.0D0 Jsgep(1,30) = 0.0D0 Jsgep(1,31) = 0.0D0 Jsgep(1,32) = 0.0D0 Jsgep(1,33) = 0.0D0 Jsgep(1,34) = 0.0D0 Jsgep(1,35) = 0.0D0 Jsgep(1,36) = 0.0D0 Jsgep(2,1) = t2 Jsgep(2,2) = t2 Jsgep(2,3) = t2 Jsgep(2,4) = 0.0D0 Jsgep(2,5) = 0.0D0 Jsgep(2,6) = 0.0D0 Jsgep(2,7) = t5 Jsgep(2,8) = t5 Jsgep(2,9) = t5 Jsgep(2,10) = 0.0D0 Jsgep(2,11) = 0.0D0 Jsgep(2,12) = 0.0D0 Jsgep(2,13) = 0.0D0 Jsgep(2,14) = 0.0D0 Jsgep(2,15) = 0.0D0 Jsgep(2,16) = 0.0D0 Jsgep(2,17) = 0.0D0 Jsgep(2,18) = 0.0D0 Jsgep(2,19) = 0.0D0 Jsgep(2,20) = 0.0D0 Jsgep(2,21) = 0.0D0 Jsgep(2,22) = 0.0D0 Jsgep(2,23) = 0.0D0 Jsgep(2,24) = 0.0D0 Jsgep(2,25) = 0.0D0 Jsgep(2,26) = 0.0D0 Jsgep(2,27) = 0.0D0 Jsgep(2,28) = 0.0D0 Jsgep(2,29) = 0.0D0 Jsgep(2,30) = 0.0D0 Jsgep(2,31) = 0.0D0 Jsgep(2,32) = 0.0D0 Jsgep(2,33) = 0.0D0 Jsgep(2,34) = 0.0D0 Jsgep(2,35) = 0.0D0 Jsgep(2,36) = 0.0D0 Jsgep(3,1) = t2 Jsgep(3,2) = t2 Jsgep(3,3) = t2 Jsgep(3,4) = 0.0D0 Jsgep(3,5) = 0.0D0 Jsgep(3,6) = 0.0D0 Jsgep(3,7) = t5 Jsgep(3,8) = t5 Jsgep(3,9) = t5 Jsgep(3,10) = 0.0D0 Jsgep(3,11) = 0.0D0 Jsgep(3,12) = 0.0D0 Jsgep(3,13) = 0.0D0 Jsgep(3,14) = 0.0D0 Jsgep(3,15) = 0.0D0 Jsgep(3,16) = 0.0D0 Jsgep(3,17) = 0.0D0 Jsgep(3,18) = 0.0D0 Jsgep(3,19) = 0.0D0 Jsgep(3,20) = 0.0D0 Jsgep(3,21) = 0.0D0 Jsgep(3,22) = 0.0D0 Jsgep(3,23) = 0.0D0 Jsgep(3,24) = 0.0D0 Jsgep(3,25) = 0.0D0 Jsgep(3,26) = 0.0D0 Jsgep(3,27) = 0.0D0 Jsgep(3,28) = 0.0D0 Jsgep(3,29) = 0.0D0 Jsgep(3,30) = 0.0D0 Jsgep(3,31) = 0.0D0 Jsgep(3,32) = 0.0D0 Jsgep(3,33) = 0.0D0 Jsgep(3,34) = 0.0D0 Jsgep(3,35) = 0.0D0 Jsgep(3,36) = 0.0D0 Jsgep(4,1) = 0.0D0 Jsgep(4,2) = 0.0D0 Jsgep(4,3) = 0.0D0 Jsgep(4,4) = 0.0D0 Jsgep(4,5) = 0.0D0 Jsgep(4,6) = 0.0D0 Jsgep(4,7) = 0.0D0 Jsgep(4,8) = 0.0D0 Jsgep(4,9) = 0.0D0 Jsgep(4,10) = t10 Jsgep(4,11) = t10 Jsgep(4,12) = t10 Jsgep(4,13) = 0.0D0 Jsgep(4,14) = 0.0D0 Jsgep(4,15) = 0.0D0 Jsgep(4,16) = 0.0D0 Jsgep(4,17) = 0.0D0 Jsgep(4,18) = 0.0D0 Jsgep(4,19) = 0.0D0 Jsgep(4,20) = 0.0D0 Jsgep(4,21) = 0.0D0 Jsgep(4,22) = 0.0D0 Jsgep(4,23) = 0.0D0 Jsgep(4,24) = 0.0D0 Jsgep(4,25) = t14 Jsgep(4,26) = t14 Jsgep(4,27) = t14 Jsgep(4,28) = 0.0D0 Jsgep(4,29) = 0.0D0 Jsgep(4,30) = 0.0D0 Jsgep(4,31) = 0.0D0 Jsgep(4,32) = 0.0D0 Jsgep(4,33) = 0.0D0 Jsgep(4,34) = 0.0D0 Jsgep(4,35) = 0.0D0 Jsgep(4,36) = 0.0D0 Jsgep(5,1) = 0.0D0 Jsgep(5,2) = 0.0D0 Jsgep(5,3) = 0.0D0 Jsgep(5,4) = 0.0D0 Jsgep(5,5) = 0.0D0 Jsgep(5,6) = 0.0D0 Jsgep(5,7) = 0.0D0 Jsgep(5,8) = 0.0D0 Jsgep(5,9) = 0.0D0 Jsgep(5,10) = t10 Jsgep(5,11) = t10 Jsgep(5,12) = t10 Jsgep(5,13) = 0.0D0 Jsgep(5,14) = 0.0D0 Jsgep(5,15) = 0.0D0 Jsgep(5,16) = 0.0D0 Jsgep(5,17) = 0.0D0 Jsgep(5,18) = 0.0D0 Jsgep(5,19) = 0.0D0 Jsgep(5,20) = 0.0D0 Jsgep(5,21) = 0.0D0 Jsgep(5,22) = 0.0D0 Jsgep(5,23) = 0.0D0 Jsgep(5,24) = 0.0D0 Jsgep(5,25) = t14 Jsgep(5,26) = t14 Jsgep(5,27) = t14 Jsgep(5,28) = 0.0D0 Jsgep(5,29) = 0.0D0 Jsgep(5,30) = 0.0D0 Jsgep(5,31) = 0.0D0 Jsgep(5,32) = 0.0D0 Jsgep(5,33) = 0.0D0 Jsgep(5,34) = 0.0D0 Jsgep(5,35) = 0.0D0 Jsgep(5,36) = 0.0D0 Jsgep(6,1) = 0.0D0 Jsgep(6,2) = 0.0D0 Jsgep(6,3) = 0.0D0 Jsgep(6,4) = 0.0D0 Jsgep(6,5) = 0.0D0 Jsgep(6,6) = 0.0D0 Jsgep(6,7) = 0.0D0 Jsgep(6,8) = 0.0D0 Jsgep(6,9) = 0.0D0 Jsgep(6,10) = t10 Jsgep(6,11) = t10 Jsgep(6,12) = t10 Jsgep(6,13) = 0.0D0 Jsgep(6,14) = 0.0D0 Jsgep(6,15) = 0.0D0 Jsgep(6,16) = 0.0D0 Jsgep(6,17) = 0.0D0 Jsgep(6,18) = 0.0D0 Jsgep(6,19) = 0.0D0 Jsgep(6,20) = 0.0D0 Jsgep(6,21) = 0.0D0 Jsgep(6,22) = 0.0D0 Jsgep(6,23) = 0.0D0 Jsgep(6,24) = 0.0D0 Jsgep(6,25) = t14 Jsgep(6,26) = t14 Jsgep(6,27) = t14 Jsgep(6,28) = 0.0D0 Jsgep(6,29) = 0.0D0 Jsgep(6,30) = 0.0D0 Jsgep(6,31) = 0.0D0 Jsgep(6,32) = 0.0D0 Jsgep(6,33) = 0.0D0 Jsgep(6,34) = 0.0D0 Jsgep(6,35) = 0.0D0 Jsgep(6,36) = 0.0D0 Jsgep(7,1) = 0.0D0 Jsgep(7,2) = 0.0D0 Jsgep(7,3) = 0.0D0 Jsgep(7,4) = 0.0D0 Jsgep(7,5) = 0.0D0 Jsgep(7,6) = 0.0D0 Jsgep(7,7) = 0.0D0 Jsgep(7,8) = 0.0D0 Jsgep(7,9) = 0.0D0 Jsgep(7,10) = t10 Jsgep(7,11) = t10 Jsgep(7,12) = t10 Jsgep(7,13) = 0.0D0 Jsgep(7,14) = 0.0D0 Jsgep(7,15) = 0.0D0 Jsgep(7,16) = 0.0D0 Jsgep(7,17) = 0.0D0 Jsgep(7,18) = 0.0D0 Jsgep(7,19) = 0.0D0 Jsgep(7,20) = 0.0D0 Jsgep(7,21) = 0.0D0 Jsgep(7,22) = 0.0D0 Jsgep(7,23) = 0.0D0 Jsgep(7,24) = 0.0D0 Jsgep(7,25) = t14 Jsgep(7,26) = t14 Jsgep(7,27) = t14 Jsgep(7,28) = 0.0D0 Jsgep(7,29) = 0.0D0 Jsgep(7,30) = 0.0D0 Jsgep(7,31) = 0.0D0 Jsgep(7,32) = 0.0D0 Jsgep(7,33) = 0.0D0 Jsgep(7,34) = 0.0D0 Jsgep(7,35) = 0.0D0 Jsgep(7,36) = 0.0D0 Jsgep(8,1) = 0.0D0 Jsgep(8,2) = 0.0D0 Jsgep(8,3) = 0.0D0 Jsgep(8,4) = 0.0D0 Jsgep(8,5) = 0.0D0 Jsgep(8,6) = 0.0D0 Jsgep(8,7) = 0.0D0 Jsgep(8,8) = 0.0D0 Jsgep(8,9) = 0.0D0 Jsgep(8,10) = t17 Jsgep(8,11) = t17 Jsgep(8,12) = t17 Jsgep(8,13) = 0.0D0 Jsgep(8,14) = 0.0D0 Jsgep(8,15) = 0.0D0 Jsgep(8,16) = 0.0D0 Jsgep(8,17) = 0.0D0 Jsgep(8,18) = 0.0D0 Jsgep(8,19) = 0.0D0 Jsgep(8,20) = 0.0D0 Jsgep(8,21) = 0.0D0 Jsgep(8,22) = 0.0D0 Jsgep(8,23) = 0.0D0 Jsgep(8,24) = 0.0D0 Jsgep(8,25) = t19 Jsgep(8,26) = t19 Jsgep(8,27) = t19 Jsgep(8,28) = 0.0D0 Jsgep(8,29) = 0.0D0 Jsgep(8,30) = 0.0D0 Jsgep(8,31) = 0.0D0 Jsgep(8,32) = 0.0D0 Jsgep(8,33) = 0.0D0 Jsgep(8,34) = 0.0D0 Jsgep(8,35) = 0.0D0 Jsgep(8,36) = 0.0D0 Jsgep(9,1) = 0.0D0 Jsgep(9,2) = 0.0D0 Jsgep(9,3) = 0.0D0 Jsgep(9,4) = 0.0D0 Jsgep(9,5) = 0.0D0 Jsgep(9,6) = 0.0D0 Jsgep(9,7) = 0.0D0 Jsgep(9,8) = 0.0D0 Jsgep(9,9) = 0.0D0 Jsgep(9,10) = t10 Jsgep(9,11) = t10 Jsgep(9,12) = t10 Jsgep(9,13) = 0.0D0 Jsgep(9,14) = 0.0D0 Jsgep(9,15) = 0.0D0 Jsgep(9,16) = 0.0D0 Jsgep(9,17) = 0.0D0 Jsgep(9,18) = 0.0D0 Jsgep(9,19) = 0.0D0 Jsgep(9,20) = 0.0D0 Jsgep(9,21) = 0.0D0 Jsgep(9,22) = 0.0D0 Jsgep(9,23) = 0.0D0 Jsgep(9,24) = 0.0D0 Jsgep(9,25) = t14 Jsgep(9,26) = t14 Jsgep(9,27) = t14 Jsgep(9,28) = 0.0D0 Jsgep(9,29) = 0.0D0 Jsgep(9,30) = 0.0D0 Jsgep(9,31) = 0.0D0 Jsgep(9,32) = 0.0D0 Jsgep(9,33) = 0.0D0 Jsgep(9,34) = 0.0D0 Jsgep(9,35) = 0.0D0 Jsgep(9,36) = 0.0D0 Jsgep(10,1) = 0.0D0 Jsgep(10,2) = 0.0D0 Jsgep(10,3) = 0.0D0 Jsgep(10,4) = 0.0D0 Jsgep(10,5) = 0.0D0 Jsgep(10,6) = 0.0D0 Jsgep(10,7) = 0.0D0 Jsgep(10,8) = 0.0D0 Jsgep(10,9) = 0.0D0 Jsgep(10,10) = t10 Jsgep(10,11) = t10 Jsgep(10,12) = t10 Jsgep(10,13) = 0.0D0 Jsgep(10,14) = 0.0D0 Jsgep(10,15) = 0.0D0 Jsgep(10,16) = 0.0D0 Jsgep(10,17) = 0.0D0 Jsgep(10,18) = 0.0D0 Jsgep(10,19) = 0.0D0 Jsgep(10,20) = 0.0D0 Jsgep(10,21) = 0.0D0 Jsgep(10,22) = 0.0D0 Jsgep(10,23) = 0.0D0 Jsgep(10,24) = 0.0D0 Jsgep(10,25) = t14 Jsgep(10,26) = t14 Jsgep(10,27) = t14 Jsgep(10,28) = 0.0D0 Jsgep(10,29) = 0.0D0 Jsgep(10,30) = 0.0D0 Jsgep(10,31) = 0.0D0 Jsgep(10,32) = 0.0D0 Jsgep(10,33) = 0.0D0 Jsgep(10,34) = 0.0D0 Jsgep(10,35) = 0.0D0 Jsgep(10,36) = 0.0D0 Jsgep(11,1) = 0.0D0 Jsgep(11,2) = 0.0D0 Jsgep(11,3) = 0.0D0 Jsgep(11,4) = 0.0D0 Jsgep(11,5) = 0.0D0 Jsgep(11,6) = 0.0D0 Jsgep(11,7) = 0.0D0 Jsgep(11,8) = 0.0D0 Jsgep(11,9) = 0.0D0 Jsgep(11,10) = t10 Jsgep(11,11) = t10 Jsgep(11,12) = t10 Jsgep(11,13) = 0.0D0 Jsgep(11,14) = 0.0D0 Jsgep(11,15) = 0.0D0 Jsgep(11,16) = 0.0D0 Jsgep(11,17) = 0.0D0 Jsgep(11,18) = 0.0D0 Jsgep(11,19) = 0.0D0 Jsgep(11,20) = 0.0D0 Jsgep(11,21) = 0.0D0 Jsgep(11,22) = 0.0D0 Jsgep(11,23) = 0.0D0 Jsgep(11,24) = 0.0D0 Jsgep(11,25) = t14 Jsgep(11,26) = t14 Jsgep(11,27) = t14 Jsgep(11,28) = 0.0D0 Jsgep(11,29) = 0.0D0 Jsgep(11,30) = 0.0D0 Jsgep(11,31) = 0.0D0 Jsgep(11,32) = 0.0D0 Jsgep(11,33) = 0.0D0 Jsgep(11,34) = 0.0D0 Jsgep(11,35) = 0.0D0 Jsgep(11,36) = 0.0D0 Jsgep(12,1) = 0.0D0 Jsgep(12,2) = 0.0D0 Jsgep(12,3) = 0.0D0 Jsgep(12,4) = 0.0D0 Jsgep(12,5) = 0.0D0 Jsgep(12,6) = 0.0D0 Jsgep(12,7) = 0.0D0 Jsgep(12,8) = 0.0D0 Jsgep(12,9) = 0.0D0 Jsgep(12,10) = t10 Jsgep(12,11) = t10 Jsgep(12,12) = t10 Jsgep(12,13) = 0.0D0 Jsgep(12,14) = 0.0D0 Jsgep(12,15) = 0.0D0 Jsgep(12,16) = 0.0D0 Jsgep(12,17) = 0.0D0 Jsgep(12,18) = 0.0D0 Jsgep(12,19) = 0.0D0 Jsgep(12,20) = 0.0D0 Jsgep(12,21) = 0.0D0 Jsgep(12,22) = 0.0D0 Jsgep(12,23) = 0.0D0 Jsgep(12,24) = 0.0D0 Jsgep(12,25) = t14 Jsgep(12,26) = t14 Jsgep(12,27) = t14 Jsgep(12,28) = 0.0D0 Jsgep(12,29) = 0.0D0 Jsgep(12,30) = 0.0D0 Jsgep(12,31) = 0.0D0 Jsgep(12,32) = 0.0D0 Jsgep(12,33) = 0.0D0 Jsgep(12,34) = 0.0D0 Jsgep(12,35) = 0.0D0 Jsgep(12,36) = 0.0D0 return end