C A1TI3D SOURCE FD218221 24/02/07 21:15:02 11834 subroutine a1ti3d(NBRINC,NDIMG,ak,bk,am,JEA,ngf,AA,ERR1,affiche) c calcul du syteme d equations pour le tir visco THM lineaire avec une inclusion implicit real*8 (a-h,o-z) implicit integer (i-n) integer NDIMG,ERR1 c Jacobienne dVe/dVA real*8 JEA(NDIMG,NDIMG+12) c coeff de fluage integer NBRINC c fluage real*8 ak(NDIMG),bk(NDIMG),am(NDIMG) integer ngf real*8 aa(ngf,ngf+1) logical affiche if((ngf.lt.12).or.(NDIMG.ne.12)) then print*,'Redimensionner ngf>12 dans incl3d pour a1ti3d' print*,'ngf:',ngf print*,'ndimg:',ndimg err1=1 return end if if(affiche) then print*,'dans a1ti3d, AK BK AM' do i=1,12 write(*,'(3E10.3)') AK(i),BK(i),AM(i) end do print*,'Dans a1ti3d, JEA pour une inclusion' do i=1,12 write (*,11) (JEA(i,j),j=1,20) 11 format (20E10.3) end do end if t2 = 0.2D1 * ak(1) + am(1) t7 = 0.2D1 * ak(2) + am(2) t11 = 0.2D1 * ak(3) + am(3) t15 = 0.2D1 * ak(4) + am(4) t19 = 0.2D1 * ak(5) + am(5) t23 = 0.2D1 * ak(6) + am(6) t27 = 0.2D1 * ak(7) + am(7) t31 = 0.2D1 * ak(8) + am(8) t35 = 0.2D1 * ak(9) + am(9) t39 = 0.2D1 * ak(10) + am(10) t43 = 0.2D1 * ak(11) + am(11) t47 = 0.2D1 * ak(12) + am(12) AA(1,1) = 0.1D1 - t2 * JEA(1,1) / 0.2D1 AA(1,2) = -t7 * JEA(1,2) / 0.2D1 AA(1,3) = -t11 * JEA(1,3) / 0.2D1 AA(1,4) = -t15 * JEA(1,4) / 0.2D1 AA(1,5) = -t19 * JEA(1,5) / 0.2D1 AA(1,6) = -t23 * JEA(1,6) / 0.2D1 AA(1,7) = -t27 * JEA(1,7) / 0.2D1 AA(1,8) = -t31 * JEA(1,8) / 0.2D1 AA(1,9) = -t35 * JEA(1,9) / 0.2D1 AA(1,10) = -t39 * JEA(1,10) / 0.2D1 AA(1,11) = -t43 * JEA(1,11) / 0.2D1 AA(1,12) = -t47 * JEA(1,12) / 0.2D1 AA(2,1) = -t2 * JEA(2,1) / 0.2D1 AA(2,2) = 0.1D1 - t7 * JEA(2,2) / 0.2D1 AA(2,3) = -t11 * JEA(2,3) / 0.2D1 AA(2,4) = -t15 * JEA(2,4) / 0.2D1 AA(2,5) = -t19 * JEA(2,5) / 0.2D1 AA(2,6) = -t23 * JEA(2,6) / 0.2D1 AA(2,7) = -t27 * JEA(2,7) / 0.2D1 AA(2,8) = -t31 * JEA(2,8) / 0.2D1 AA(2,9) = -t35 * JEA(2,9) / 0.2D1 AA(2,10) = -t39 * JEA(2,10) / 0.2D1 AA(2,11) = -t43 * JEA(2,11) / 0.2D1 AA(2,12) = -t47 * JEA(2,12) / 0.2D1 AA(3,1) = -t2 * JEA(3,1) / 0.2D1 AA(3,2) = -t7 * JEA(3,2) / 0.2D1 AA(3,3) = 0.1D1 - t11 * JEA(3,3) / 0.2D1 AA(3,4) = -t15 * JEA(3,4) / 0.2D1 AA(3,5) = -t19 * JEA(3,5) / 0.2D1 AA(3,6) = -t23 * JEA(3,6) / 0.2D1 AA(3,7) = -t27 * JEA(3,7) / 0.2D1 AA(3,8) = -t31 * JEA(3,8) / 0.2D1 AA(3,9) = -t35 * JEA(3,9) / 0.2D1 AA(3,10) = -t39 * JEA(3,10) / 0.2D1 AA(3,11) = -t43 * JEA(3,11) / 0.2D1 AA(3,12) = -t47 * JEA(3,12) / 0.2D1 AA(4,1) = -t2 * JEA(4,1) / 0.2D1 AA(4,2) = -t7 * JEA(4,2) / 0.2D1 AA(4,3) = -t11 * JEA(4,3) / 0.2D1 AA(4,4) = 0.1D1 - t15 * JEA(4,4) / 0.2D1 AA(4,5) = -t19 * JEA(4,5) / 0.2D1 AA(4,6) = -t23 * JEA(4,6) / 0.2D1 AA(4,7) = -t27 * JEA(4,7) / 0.2D1 AA(4,8) = -t31 * JEA(4,8) / 0.2D1 AA(4,9) = -t35 * JEA(4,9) / 0.2D1 AA(4,10) = -t39 * JEA(4,10) / 0.2D1 AA(4,11) = -t43 * JEA(4,11) / 0.2D1 AA(4,12) = -t47 * JEA(4,12) / 0.2D1 AA(5,1) = -t2 * JEA(5,1) / 0.2D1 AA(5,2) = -t7 * JEA(5,2) / 0.2D1 AA(5,3) = -t11 * JEA(5,3) / 0.2D1 AA(5,4) = -t15 * JEA(5,4) / 0.2D1 AA(5,5) = 0.1D1 - t19 * JEA(5,5) / 0.2D1 AA(5,6) = -t23 * JEA(5,6) / 0.2D1 AA(5,7) = -t27 * JEA(5,7) / 0.2D1 AA(5,8) = -t31 * JEA(5,8) / 0.2D1 AA(5,9) = -t35 * JEA(5,9) / 0.2D1 AA(5,10) = -t39 * JEA(5,10) / 0.2D1 AA(5,11) = -t43 * JEA(5,11) / 0.2D1 AA(5,12) = -t47 * JEA(5,12) / 0.2D1 AA(6,1) = -t2 * JEA(6,1) / 0.2D1 AA(6,2) = -t7 * JEA(6,2) / 0.2D1 AA(6,3) = -t11 * JEA(6,3) / 0.2D1 AA(6,4) = -t15 * JEA(6,4) / 0.2D1 AA(6,5) = -t19 * JEA(6,5) / 0.2D1 AA(6,6) = 0.1D1 - t23 * JEA(6,6) / 0.2D1 AA(6,7) = -t27 * JEA(6,7) / 0.2D1 AA(6,8) = -t31 * JEA(6,8) / 0.2D1 AA(6,9) = -t35 * JEA(6,9) / 0.2D1 AA(6,10) = -t39 * JEA(6,10) / 0.2D1 AA(6,11) = -t43 * JEA(6,11) / 0.2D1 AA(6,12) = -t47 * JEA(6,12) / 0.2D1 AA(7,1) = -t2 * JEA(7,1) / 0.2D1 AA(7,2) = -t7 * JEA(7,2) / 0.2D1 AA(7,3) = -t11 * JEA(7,3) / 0.2D1 AA(7,4) = -t15 * JEA(7,4) / 0.2D1 AA(7,5) = -t19 * JEA(7,5) / 0.2D1 AA(7,6) = -t23 * JEA(7,6) / 0.2D1 AA(7,7) = 0.1D1 - t27 * JEA(7,7) / 0.2D1 AA(7,8) = -t31 * JEA(7,8) / 0.2D1 AA(7,9) = -t35 * JEA(7,9) / 0.2D1 AA(7,10) = -t39 * JEA(7,10) / 0.2D1 AA(7,11) = -t43 * JEA(7,11) / 0.2D1 AA(7,12) = -t47 * JEA(7,12) / 0.2D1 AA(8,1) = -t2 * JEA(8,1) / 0.2D1 AA(8,2) = -t7 * JEA(8,2) / 0.2D1 AA(8,3) = -t11 * JEA(8,3) / 0.2D1 AA(8,4) = -t15 * JEA(8,4) / 0.2D1 AA(8,5) = -t19 * JEA(8,5) / 0.2D1 AA(8,6) = -t23 * JEA(8,6) / 0.2D1 AA(8,7) = -t27 * JEA(8,7) / 0.2D1 AA(8,8) = 0.1D1 - t31 * JEA(8,8) / 0.2D1 AA(8,9) = -t35 * JEA(8,9) / 0.2D1 AA(8,10) = -t39 * JEA(8,10) / 0.2D1 AA(8,11) = -t43 * JEA(8,11) / 0.2D1 AA(8,12) = -t47 * JEA(8,12) / 0.2D1 AA(9,1) = -t2 * JEA(9,1) / 0.2D1 AA(9,2) = -t7 * JEA(9,2) / 0.2D1 AA(9,3) = -t11 * JEA(9,3) / 0.2D1 AA(9,4) = -t15 * JEA(9,4) / 0.2D1 AA(9,5) = -t19 * JEA(9,5) / 0.2D1 AA(9,6) = -t23 * JEA(9,6) / 0.2D1 AA(9,7) = -t27 * JEA(9,7) / 0.2D1 AA(9,8) = -t31 * JEA(9,8) / 0.2D1 AA(9,9) = 0.1D1 - t35 * JEA(9,9) / 0.2D1 AA(9,10) = -t39 * JEA(9,10) / 0.2D1 AA(9,11) = -t43 * JEA(9,11) / 0.2D1 AA(9,12) = -t47 * JEA(9,12) / 0.2D1 AA(10,1) = -t2 * JEA(10,1) / 0.2D1 AA(10,2) = -t7 * JEA(10,2) / 0.2D1 AA(10,3) = -t11 * JEA(10,3) / 0.2D1 AA(10,4) = -t15 * JEA(10,4) / 0.2D1 AA(10,5) = -t19 * JEA(10,5) / 0.2D1 AA(10,6) = -t23 * JEA(10,6) / 0.2D1 AA(10,7) = -t27 * JEA(10,7) / 0.2D1 AA(10,8) = -t31 * JEA(10,8) / 0.2D1 AA(10,9) = -t35 * JEA(10,9) / 0.2D1 AA(10,10) = 0.1D1 - t39 * JEA(10,10) / 0.2D1 AA(10,11) = -t43 * JEA(10,11) / 0.2D1 AA(10,12) = -t47 * JEA(10,12) / 0.2D1 AA(11,1) = -t2 * JEA(11,1) / 0.2D1 AA(11,2) = -t7 * JEA(11,2) / 0.2D1 AA(11,3) = -t11 * JEA(11,3) / 0.2D1 AA(11,4) = -t15 * JEA(11,4) / 0.2D1 AA(11,5) = -t19 * JEA(11,5) / 0.2D1 AA(11,6) = -t23 * JEA(11,6) / 0.2D1 AA(11,7) = -t27 * JEA(11,7) / 0.2D1 AA(11,8) = -t31 * JEA(11,8) / 0.2D1 AA(11,9) = -t35 * JEA(11,9) / 0.2D1 AA(11,10) = -t39 * JEA(11,10) / 0.2D1 AA(11,11) = 0.1D1 - t43 * JEA(11,11) / 0.2D1 AA(11,12) = -t47 * JEA(11,12) / 0.2D1 AA(12,1) = -t2 * JEA(12,1) / 0.2D1 AA(12,2) = -t7 * JEA(12,2) / 0.2D1 AA(12,3) = -t11 * JEA(12,3) / 0.2D1 AA(12,4) = -t15 * JEA(12,4) / 0.2D1 AA(12,5) = -t19 * JEA(12,5) / 0.2D1 AA(12,6) = -t23 * JEA(12,6) / 0.2D1 AA(12,7) = -t27 * JEA(12,7) / 0.2D1 AA(12,8) = -t31 * JEA(12,8) / 0.2D1 AA(12,9) = -t35 * JEA(12,9) / 0.2D1 AA(12,10) = -t39 * JEA(12,10) / 0.2D1 AA(12,11) = -t43 * JEA(12,11) / 0.2D1 AA(12,12) = 0.1D1 - t47 * JEA(12,12) / 0.2D1 if(affiche) then print*, 'Dans a1ti3d : matrice de couplage' do i=1,12 write (*,'(12E10.3)') (AA(i,j),j=1,12) end do end if return end