C JEPG3D SOURCE FD218221 24/02/07 21:15:17 11834 subroutine jepg3d(NBRINC,NINC,NDIMG,NDIMA,NDIMPL,FRAC,GCH,Mbg,bg, # JEA,cg,AFFICHE,ERR1) c calcul du jacobien d eps e / d epsp generalise (NDIMG*NDIMPL) c stocke dans cg() (cas avec une seule inclusion) implicit real*8 (a-h,o-z) implicit integer (i-n) c nbre d inc max et numero de la contrainte du critere actif de traction integer NBRINC,ERR1,NDIMG,NINC,NDIMPL c Jacobienne dVe/dVA real*8 JEA(NDIMG,NDIMA) c fraction des phases real*8 FRAC(0:NBRINC) c indicateur de diffusion produits chimiques dans fissure real*8 GCH(0:NBRINC) c coeff de biot et module de biot des gels real*8 bg(0:NBRINC,0:1),Mbg(0:NBRINC,0:1) c jacobienne Je1epg=cg (cf maple) real*8 cg(NDIMG,NDIMPL) c options d affichage logical affiche,affiche_local c fraction inclusion real*8 f if(NBRINC.ne.1) then print*,'Nbr inc inattendu ds jbgp3d' err1=1 return end if affiche_local=affiche err1=0 f=frac(1) c pour la suite cf fichier maple t2 = Gch(1) * JEA(1,21) * bg(1,1) t7 = 0.1D1 - f t8 = t7 * Gch(1) t9 = t8 * JEA(1,21) * bg(1,1) t13 = JEA(1,24) * bg(0,1) t16 = Gch(0) * (0.1D1 - 0.3D1 / 0.2D1 * f) t17 = t16 * t13 t23 = 0.1D1 / t7 * f * Gch(0) t24 = t23 * t13 t29 = Gch(1) * JEA(2,21) * bg(1,1) t34 = t8 * JEA(2,21) * bg(1,1) t38 = JEA(2,24) * bg(0,1) t39 = t16 * t38 t43 = t23 * t38 t48 = Gch(1) * JEA(3,21) * bg(1,1) t53 = t8 * JEA(3,21) * bg(1,1) t57 = JEA(3,24) * bg(0,1) t58 = t16 * t57 t62 = t23 * t57 t67 = Gch(1) * JEA(4,21) * bg(1,1) t72 = t8 * JEA(4,21) * bg(1,1) t76 = JEA(4,24) * bg(0,1) t77 = t16 * t76 t81 = t23 * t76 t86 = Gch(1) * JEA(5,21) * bg(1,1) t91 = t8 * JEA(5,21) * bg(1,1) t95 = JEA(5,24) * bg(0,1) t96 = t16 * t95 t100 = t23 * t95 t105 = Gch(1) * JEA(6,21) * bg(1,1) t110 = t8 * JEA(6,21) * bg(1,1) t114 = JEA(6,24) * bg(0,1) t115 = t16 * t114 t119 = t23 * t114 t124 = Gch(1) * JEA(7,21) * bg(1,1) t129 = t8 * JEA(7,21) * bg(1,1) t133 = JEA(7,24) * bg(0,1) t134 = t16 * t133 t138 = t23 * t133 t143 = Gch(1) * JEA(8,21) * bg(1,1) t148 = t8 * JEA(8,21) * bg(1,1) t152 = JEA(8,24) * bg(0,1) t153 = t16 * t152 t157 = t23 * t152 t162 = Gch(1) * JEA(9,21) * bg(1,1) t167 = t8 * JEA(9,21) * bg(1,1) t171 = JEA(9,24) * bg(0,1) t172 = t16 * t171 t176 = t23 * t171 t181 = Gch(1) * JEA(10,21) * bg(1,1) t186 = t8 * JEA(10,21) * bg(1,1) t190 = JEA(10,24) * bg(0,1) t191 = t16 * t190 t195 = t23 * t190 t200 = Gch(1) * JEA(11,21) * bg(1,1) t205 = t8 * JEA(11,21) * bg(1,1) t209 = JEA(11,24) * bg(0,1) t210 = t16 * t209 t214 = t23 * t209 t219 = Gch(1) * JEA(12,21) * bg(1,1) t224 = t8 * JEA(12,21) * bg(1,1) t228 = JEA(12,24) * bg(0,1) t229 = t16 * t228 t233 = t23 * t228 cg(1,1) = t2 + JEA(1,1) cg(1,2) = t2 + JEA(1,2) cg(1,3) = t2 + JEA(1,3) cg(1,4) = JEA(1,4) cg(1,5) = JEA(1,5) cg(1,6) = JEA(1,6) cg(1,7) = JEA(1,7) + t9 cg(1,8) = JEA(1,8) + t9 cg(1,9) = JEA(1,9) + t9 cg(1,10) = JEA(1,10) + t17 cg(1,11) = JEA(1,11) + t17 cg(1,12) = JEA(1,12) + t17 cg(1,13) = JEA(1,1) cg(1,14) = JEA(1,2) cg(1,15) = JEA(1,3) cg(1,16) = JEA(1,4) cg(1,17) = JEA(1,5) cg(1,18) = JEA(1,6) cg(1,19) = JEA(1,7) cg(1,20) = JEA(1,8) cg(1,21) = JEA(1,9) cg(1,22) = JEA(1,10) cg(1,23) = JEA(1,11) cg(1,24) = JEA(1,12) cg(1,25) = JEA(1,1) + t24 cg(1,26) = JEA(1,2) + t24 cg(1,27) = JEA(1,3) + t24 cg(1,28) = JEA(1,4) cg(1,29) = JEA(1,5) cg(1,30) = JEA(1,6) cg(1,31) = JEA(1,7) cg(1,32) = JEA(1,8) cg(1,33) = JEA(1,9) cg(1,34) = JEA(1,10) cg(1,35) = JEA(1,11) cg(1,36) = JEA(1,12) cg(2,1) = t29 + JEA(2,1) cg(2,2) = t29 + JEA(2,2) cg(2,3) = t29 + JEA(2,3) cg(2,4) = JEA(2,4) cg(2,5) = JEA(2,5) cg(2,6) = JEA(2,6) cg(2,7) = JEA(2,7) + t34 cg(2,8) = JEA(2,8) + t34 cg(2,9) = JEA(2,9) + t34 cg(2,10) = JEA(2,10) + t39 cg(2,11) = JEA(2,11) + t39 cg(2,12) = JEA(2,12) + t39 cg(2,13) = JEA(2,1) cg(2,14) = JEA(2,2) cg(2,15) = JEA(2,3) cg(2,16) = JEA(2,4) cg(2,17) = JEA(2,5) cg(2,18) = JEA(2,6) cg(2,19) = JEA(2,7) cg(2,20) = JEA(2,8) cg(2,21) = JEA(2,9) cg(2,22) = JEA(2,10) cg(2,23) = JEA(2,11) cg(2,24) = JEA(2,12) cg(2,25) = JEA(2,1) + t43 cg(2,26) = JEA(2,2) + t43 cg(2,27) = JEA(2,3) + t43 cg(2,28) = JEA(2,4) cg(2,29) = JEA(2,5) cg(2,30) = JEA(2,6) cg(2,31) = JEA(2,7) cg(2,32) = JEA(2,8) cg(2,33) = JEA(2,9) cg(2,34) = JEA(2,10) cg(2,35) = JEA(2,11) cg(2,36) = JEA(2,12) cg(3,1) = t48 + JEA(3,1) cg(3,2) = t48 + JEA(3,2) cg(3,3) = t48 + JEA(3,3) cg(3,4) = JEA(3,4) cg(3,5) = JEA(3,5) cg(3,6) = JEA(3,6) cg(3,7) = JEA(3,7) + t53 cg(3,8) = JEA(3,8) + t53 cg(3,9) = JEA(3,9) + t53 cg(3,10) = JEA(3,10) + t58 cg(3,11) = JEA(3,11) + t58 cg(3,12) = JEA(3,12) + t58 cg(3,13) = JEA(3,1) cg(3,14) = JEA(3,2) cg(3,15) = JEA(3,3) cg(3,16) = JEA(3,4) cg(3,17) = JEA(3,5) cg(3,18) = JEA(3,6) cg(3,19) = JEA(3,7) cg(3,20) = JEA(3,8) cg(3,21) = JEA(3,9) cg(3,22) = JEA(3,10) cg(3,23) = JEA(3,11) cg(3,24) = JEA(3,12) cg(3,25) = JEA(3,1) + t62 cg(3,26) = JEA(3,2) + t62 cg(3,27) = JEA(3,3) + t62 cg(3,28) = JEA(3,4) cg(3,29) = JEA(3,5) cg(3,30) = JEA(3,6) cg(3,31) = JEA(3,7) cg(3,32) = JEA(3,8) cg(3,33) = JEA(3,9) cg(3,34) = JEA(3,10) cg(3,35) = JEA(3,11) cg(3,36) = JEA(3,12) cg(4,1) = t67 + JEA(4,1) cg(4,2) = t67 + JEA(4,2) cg(4,3) = t67 + JEA(4,3) cg(4,4) = JEA(4,4) cg(4,5) = JEA(4,5) cg(4,6) = JEA(4,6) cg(4,7) = JEA(4,7) + t72 cg(4,8) = JEA(4,8) + t72 cg(4,9) = JEA(4,9) + t72 cg(4,10) = JEA(4,10) + t77 cg(4,11) = JEA(4,11) + t77 cg(4,12) = JEA(4,12) + t77 cg(4,13) = JEA(4,1) cg(4,14) = JEA(4,2) cg(4,15) = JEA(4,3) cg(4,16) = JEA(4,4) cg(4,17) = JEA(4,5) cg(4,18) = JEA(4,6) cg(4,19) = JEA(4,7) cg(4,20) = JEA(4,8) cg(4,21) = JEA(4,9) cg(4,22) = JEA(4,10) cg(4,23) = JEA(4,11) cg(4,24) = JEA(4,12) cg(4,25) = JEA(4,1) + t81 cg(4,26) = JEA(4,2) + t81 cg(4,27) = JEA(4,3) + t81 cg(4,28) = JEA(4,4) cg(4,29) = JEA(4,5) cg(4,30) = JEA(4,6) cg(4,31) = JEA(4,7) cg(4,32) = JEA(4,8) cg(4,33) = JEA(4,9) cg(4,34) = JEA(4,10) cg(4,35) = JEA(4,11) cg(4,36) = JEA(4,12) cg(5,1) = t86 + JEA(5,1) cg(5,2) = t86 + JEA(5,2) cg(5,3) = t86 + JEA(5,3) cg(5,4) = JEA(5,4) cg(5,5) = JEA(5,5) cg(5,6) = JEA(5,6) cg(5,7) = JEA(5,7) + t91 cg(5,8) = JEA(5,8) + t91 cg(5,9) = JEA(5,9) + t91 cg(5,10) = JEA(5,10) + t96 cg(5,11) = JEA(5,11) + t96 cg(5,12) = JEA(5,12) + t96 cg(5,13) = JEA(5,1) cg(5,14) = JEA(5,2) cg(5,15) = JEA(5,3) cg(5,16) = JEA(5,4) cg(5,17) = JEA(5,5) cg(5,18) = JEA(5,6) cg(5,19) = JEA(5,7) cg(5,20) = JEA(5,8) cg(5,21) = JEA(5,9) cg(5,22) = JEA(5,10) cg(5,23) = JEA(5,11) cg(5,24) = JEA(5,12) cg(5,25) = JEA(5,1) + t100 cg(5,26) = JEA(5,2) + t100 cg(5,27) = JEA(5,3) + t100 cg(5,28) = JEA(5,4) cg(5,29) = JEA(5,5) cg(5,30) = JEA(5,6) cg(5,31) = JEA(5,7) cg(5,32) = JEA(5,8) cg(5,33) = JEA(5,9) cg(5,34) = JEA(5,10) cg(5,35) = JEA(5,11) cg(5,36) = JEA(5,12) cg(6,1) = t105 + JEA(6,1) cg(6,2) = t105 + JEA(6,2) cg(6,3) = t105 + JEA(6,3) cg(6,4) = JEA(6,4) cg(6,5) = JEA(6,5) cg(6,6) = JEA(6,6) cg(6,7) = JEA(6,7) + t110 cg(6,8) = JEA(6,8) + t110 cg(6,9) = JEA(6,9) + t110 cg(6,10) = JEA(6,10) + t115 cg(6,11) = JEA(6,11) + t115 cg(6,12) = JEA(6,12) + t115 cg(6,13) = JEA(6,1) cg(6,14) = JEA(6,2) cg(6,15) = JEA(6,3) cg(6,16) = JEA(6,4) cg(6,17) = JEA(6,5) cg(6,18) = JEA(6,6) cg(6,19) = JEA(6,7) cg(6,20) = JEA(6,8) cg(6,21) = JEA(6,9) cg(6,22) = JEA(6,10) cg(6,23) = JEA(6,11) cg(6,24) = JEA(6,12) cg(6,25) = JEA(6,1) + t119 cg(6,26) = JEA(6,2) + t119 cg(6,27) = JEA(6,3) + t119 cg(6,28) = JEA(6,4) cg(6,29) = JEA(6,5) cg(6,30) = JEA(6,6) cg(6,31) = JEA(6,7) cg(6,32) = JEA(6,8) cg(6,33) = JEA(6,9) cg(6,34) = JEA(6,10) cg(6,35) = JEA(6,11) cg(6,36) = JEA(6,12) cg(7,1) = t124 + JEA(7,1) cg(7,2) = t124 + JEA(7,2) cg(7,3) = t124 + JEA(7,3) cg(7,4) = JEA(7,4) cg(7,5) = JEA(7,5) cg(7,6) = JEA(7,6) cg(7,7) = JEA(7,7) + t129 cg(7,8) = JEA(7,8) + t129 cg(7,9) = JEA(7,9) + t129 cg(7,10) = JEA(7,10) + t134 cg(7,11) = JEA(7,11) + t134 cg(7,12) = JEA(7,12) + t134 cg(7,13) = JEA(7,1) cg(7,14) = JEA(7,2) cg(7,15) = JEA(7,3) cg(7,16) = JEA(7,4) cg(7,17) = JEA(7,5) cg(7,18) = JEA(7,6) cg(7,19) = JEA(7,7) cg(7,20) = JEA(7,8) cg(7,21) = JEA(7,9) cg(7,22) = JEA(7,10) cg(7,23) = JEA(7,11) cg(7,24) = JEA(7,12) cg(7,25) = JEA(7,1) + t138 cg(7,26) = JEA(7,2) + t138 cg(7,27) = JEA(7,3) + t138 cg(7,28) = JEA(7,4) cg(7,29) = JEA(7,5) cg(7,30) = JEA(7,6) cg(7,31) = JEA(7,7) cg(7,32) = JEA(7,8) cg(7,33) = JEA(7,9) cg(7,34) = JEA(7,10) cg(7,35) = JEA(7,11) cg(7,36) = JEA(7,12) cg(8,1) = t143 + JEA(8,1) cg(8,2) = t143 + JEA(8,2) cg(8,3) = t143 + JEA(8,3) cg(8,4) = JEA(8,4) cg(8,5) = JEA(8,5) cg(8,6) = JEA(8,6) cg(8,7) = JEA(8,7) + t148 cg(8,8) = JEA(8,8) + t148 cg(8,9) = JEA(8,9) + t148 cg(8,10) = JEA(8,10) + t153 cg(8,11) = JEA(8,11) + t153 cg(8,12) = JEA(8,12) + t153 cg(8,13) = JEA(8,1) cg(8,14) = JEA(8,2) cg(8,15) = JEA(8,3) cg(8,16) = JEA(8,4) cg(8,17) = JEA(8,5) cg(8,18) = JEA(8,6) cg(8,19) = JEA(8,7) cg(8,20) = JEA(8,8) cg(8,21) = JEA(8,9) cg(8,22) = JEA(8,10) cg(8,23) = JEA(8,11) cg(8,24) = JEA(8,12) cg(8,25) = JEA(8,1) + t157 cg(8,26) = JEA(8,2) + t157 cg(8,27) = JEA(8,3) + t157 cg(8,28) = JEA(8,4) cg(8,29) = JEA(8,5) cg(8,30) = JEA(8,6) cg(8,31) = JEA(8,7) cg(8,32) = JEA(8,8) cg(8,33) = JEA(8,9) cg(8,34) = JEA(8,10) cg(8,35) = JEA(8,11) cg(8,36) = JEA(8,12) cg(9,1) = t162 + JEA(9,1) cg(9,2) = t162 + JEA(9,2) cg(9,3) = t162 + JEA(9,3) cg(9,4) = JEA(9,4) cg(9,5) = JEA(9,5) cg(9,6) = JEA(9,6) cg(9,7) = JEA(9,7) + t167 cg(9,8) = JEA(9,8) + t167 cg(9,9) = JEA(9,9) + t167 cg(9,10) = JEA(9,10) + t172 cg(9,11) = JEA(9,11) + t172 cg(9,12) = JEA(9,12) + t172 cg(9,13) = JEA(9,1) cg(9,14) = JEA(9,2) cg(9,15) = JEA(9,3) cg(9,16) = JEA(9,4) cg(9,17) = JEA(9,5) cg(9,18) = JEA(9,6) cg(9,19) = JEA(9,7) cg(9,20) = JEA(9,8) cg(9,21) = JEA(9,9) cg(9,22) = JEA(9,10) cg(9,23) = JEA(9,11) cg(9,24) = JEA(9,12) cg(9,25) = JEA(9,1) + t176 cg(9,26) = JEA(9,2) + t176 cg(9,27) = JEA(9,3) + t176 cg(9,28) = JEA(9,4) cg(9,29) = JEA(9,5) cg(9,30) = JEA(9,6) cg(9,31) = JEA(9,7) cg(9,32) = JEA(9,8) cg(9,33) = JEA(9,9) cg(9,34) = JEA(9,10) cg(9,35) = JEA(9,11) cg(9,36) = JEA(9,12) cg(10,1) = t181 + JEA(10,1) cg(10,2) = t181 + JEA(10,2) cg(10,3) = t181 + JEA(10,3) cg(10,4) = JEA(10,4) cg(10,5) = JEA(10,5) cg(10,6) = JEA(10,6) cg(10,7) = JEA(10,7) + t186 cg(10,8) = JEA(10,8) + t186 cg(10,9) = JEA(10,9) + t186 cg(10,10) = JEA(10,10) + t191 cg(10,11) = JEA(10,11) + t191 cg(10,12) = JEA(10,12) + t191 cg(10,13) = JEA(10,1) cg(10,14) = JEA(10,2) cg(10,15) = JEA(10,3) cg(10,16) = JEA(10,4) cg(10,17) = JEA(10,5) cg(10,18) = JEA(10,6) cg(10,19) = JEA(10,7) cg(10,20) = JEA(10,8) cg(10,21) = JEA(10,9) cg(10,22) = JEA(10,10) cg(10,23) = JEA(10,11) cg(10,24) = JEA(10,12) cg(10,25) = JEA(10,1) + t195 cg(10,26) = JEA(10,2) + t195 cg(10,27) = JEA(10,3) + t195 cg(10,28) = JEA(10,4) cg(10,29) = JEA(10,5) cg(10,30) = JEA(10,6) cg(10,31) = JEA(10,7) cg(10,32) = JEA(10,8) cg(10,33) = JEA(10,9) cg(10,34) = JEA(10,10) cg(10,35) = JEA(10,11) cg(10,36) = JEA(10,12) cg(11,1) = t200 + JEA(11,1) cg(11,2) = t200 + JEA(11,2) cg(11,3) = t200 + JEA(11,3) cg(11,4) = JEA(11,4) cg(11,5) = JEA(11,5) cg(11,6) = JEA(11,6) cg(11,7) = JEA(11,7) + t205 cg(11,8) = JEA(11,8) + t205 cg(11,9) = JEA(11,9) + t205 cg(11,10) = JEA(11,10) + t210 cg(11,11) = JEA(11,11) + t210 cg(11,12) = JEA(11,12) + t210 cg(11,13) = JEA(11,1) cg(11,14) = JEA(11,2) cg(11,15) = JEA(11,3) cg(11,16) = JEA(11,4) cg(11,17) = JEA(11,5) cg(11,18) = JEA(11,6) cg(11,19) = JEA(11,7) cg(11,20) = JEA(11,8) cg(11,21) = JEA(11,9) cg(11,22) = JEA(11,10) cg(11,23) = JEA(11,11) cg(11,24) = JEA(11,12) cg(11,25) = JEA(11,1) + t214 cg(11,26) = JEA(11,2) + t214 cg(11,27) = JEA(11,3) + t214 cg(11,28) = JEA(11,4) cg(11,29) = JEA(11,5) cg(11,30) = JEA(11,6) cg(11,31) = JEA(11,7) cg(11,32) = JEA(11,8) cg(11,33) = JEA(11,9) cg(11,34) = JEA(11,10) cg(11,35) = JEA(11,11) cg(11,36) = JEA(11,12) cg(12,1) = t219 + JEA(12,1) cg(12,2) = t219 + JEA(12,2) cg(12,3) = t219 + JEA(12,3) cg(12,4) = JEA(12,4) cg(12,5) = JEA(12,5) cg(12,6) = JEA(12,6) cg(12,7) = JEA(12,7) + t224 cg(12,8) = JEA(12,8) + t224 cg(12,9) = JEA(12,9) + t224 cg(12,10) = JEA(12,10) + t229 cg(12,11) = JEA(12,11) + t229 cg(12,12) = JEA(12,12) + t229 cg(12,13) = JEA(12,1) cg(12,14) = JEA(12,2) cg(12,15) = JEA(12,3) cg(12,16) = JEA(12,4) cg(12,17) = JEA(12,5) cg(12,18) = JEA(12,6) cg(12,19) = JEA(12,7) cg(12,20) = JEA(12,8) cg(12,21) = JEA(12,9) cg(12,22) = JEA(12,10) cg(12,23) = JEA(12,11) cg(12,24) = JEA(12,12) cg(12,25) = JEA(12,1) + t233 cg(12,26) = JEA(12,2) + t233 cg(12,27) = JEA(12,3) + t233 cg(12,28) = JEA(12,4) cg(12,29) = JEA(12,5) cg(12,30) = JEA(12,6) cg(12,31) = JEA(12,7) cg(12,32) = JEA(12,8) cg(12,33) = JEA(12,9) cg(12,34) = JEA(12,10) cg(12,35) = JEA(12,11) cg(12,36) = JEA(12,12) return end