Télécharger la plaquette Cast3M


fr en

Télécharger a1ti3d.eso

Retour à la liste

Numérotation des lignes : 

a1ti3d

  1. C A1TI3D    SOURCE    FD218221  24/02/07    21:15:02     11834          
  2.       subroutine a1ti3d(NBRINC,NDIMG,ak,bk,am,JEA,ngf,AA,ERR1,affiche)
  3.  
  4. c     calcul du syteme d equations pour le tir visco THM lineaire avec une inclusion      
  5.  
  6.       implicit real*8 (a-h,o-z)
  7.       implicit integer (i-n)
  8.  
  9.       integer NDIMG,ERR1
  10. c     Jacobienne dVe/dVA
  11.       real*8 JEA(NDIMG,NDIMG+12)
  12. c     coeff de fluage      
  13.       integer NBRINC
  14. c     fluage
  15.       real*8 ak(NDIMG),bk(NDIMG),am(NDIMG)
  16.       integer ngf
  17.       real*8 aa(ngf,ngf+1) 
  18.       logical affiche
  19.  
  20.  
  21.       if((ngf.lt.12).or.(NDIMG.ne.12)) then
  22.         print*,'Redimensionner ngf>12 dans incl3d pour a1ti3d'
  23.         print*,'ngf:',ngf
  24.         print*,'ndimg:',ndimg
  25.         err1=1
  26.         return
  27.       end if
  28.  
  29.       if(affiche) then
  30.           print*,'dans a1ti3d, AK BK AM'
  31.           do i=1,12
  32.                 write(*,'(3E10.3)') AK(i),BK(i),AM(i)
  33.           end do
  34.           print*,'Dans a1ti3d, JEA pour une inclusion'
  35.           do i=1,12
  36.             write (*,11) (JEA(i,j),j=1,20) 
  37. 11              format (20E10.3)
  38.           end do
  39.       end if
  40.  
  41.       t2 = 0.2D1 * ak(1) + am(1)
  42.       t7 = 0.2D1 * ak(2) + am(2)
  43.       t11 = 0.2D1 * ak(3) + am(3)
  44.       t15 = 0.2D1 * ak(4) + am(4)
  45.       t19 = 0.2D1 * ak(5) + am(5)
  46.       t23 = 0.2D1 * ak(6) + am(6)
  47.       t27 = 0.2D1 * ak(7) + am(7)
  48.       t31 = 0.2D1 * ak(8) + am(8)
  49.       t35 = 0.2D1 * ak(9) + am(9)
  50.       t39 = 0.2D1 * ak(10) + am(10)
  51.       t43 = 0.2D1 * ak(11) + am(11)
  52.       t47 = 0.2D1 * ak(12) + am(12)
  53.       AA(1,1) = 0.1D1 - t2 * JEA(1,1) / 0.2D1
  54.       AA(1,2) = -t7 * JEA(1,2) / 0.2D1
  55.       AA(1,3) = -t11 * JEA(1,3) / 0.2D1
  56.       AA(1,4) = -t15 * JEA(1,4) / 0.2D1
  57.       AA(1,5) = -t19 * JEA(1,5) / 0.2D1
  58.       AA(1,6) = -t23 * JEA(1,6) / 0.2D1
  59.       AA(1,7) = -t27 * JEA(1,7) / 0.2D1
  60.       AA(1,8) = -t31 * JEA(1,8) / 0.2D1
  61.       AA(1,9) = -t35 * JEA(1,9) / 0.2D1
  62.       AA(1,10) = -t39 * JEA(1,10) / 0.2D1
  63.       AA(1,11) = -t43 * JEA(1,11) / 0.2D1
  64.       AA(1,12) = -t47 * JEA(1,12) / 0.2D1
  65.       AA(2,1) = -t2 * JEA(2,1) / 0.2D1
  66.       AA(2,2) = 0.1D1 - t7 * JEA(2,2) / 0.2D1
  67.       AA(2,3) = -t11 * JEA(2,3) / 0.2D1
  68.       AA(2,4) = -t15 * JEA(2,4) / 0.2D1
  69.       AA(2,5) = -t19 * JEA(2,5) / 0.2D1
  70.       AA(2,6) = -t23 * JEA(2,6) / 0.2D1
  71.       AA(2,7) = -t27 * JEA(2,7) / 0.2D1
  72.       AA(2,8) = -t31 * JEA(2,8) / 0.2D1
  73.       AA(2,9) = -t35 * JEA(2,9) / 0.2D1
  74.       AA(2,10) = -t39 * JEA(2,10) / 0.2D1
  75.       AA(2,11) = -t43 * JEA(2,11) / 0.2D1
  76.       AA(2,12) = -t47 * JEA(2,12) / 0.2D1
  77.       AA(3,1) = -t2 * JEA(3,1) / 0.2D1
  78.       AA(3,2) = -t7 * JEA(3,2) / 0.2D1
  79.       AA(3,3) = 0.1D1 - t11 * JEA(3,3) / 0.2D1
  80.       AA(3,4) = -t15 * JEA(3,4) / 0.2D1
  81.       AA(3,5) = -t19 * JEA(3,5) / 0.2D1
  82.       AA(3,6) = -t23 * JEA(3,6) / 0.2D1
  83.       AA(3,7) = -t27 * JEA(3,7) / 0.2D1
  84.       AA(3,8) = -t31 * JEA(3,8) / 0.2D1
  85.       AA(3,9) = -t35 * JEA(3,9) / 0.2D1
  86.       AA(3,10) = -t39 * JEA(3,10) / 0.2D1
  87.       AA(3,11) = -t43 * JEA(3,11) / 0.2D1
  88.       AA(3,12) = -t47 * JEA(3,12) / 0.2D1
  89.       AA(4,1) = -t2 * JEA(4,1) / 0.2D1
  90.       AA(4,2) = -t7 * JEA(4,2) / 0.2D1
  91.       AA(4,3) = -t11 * JEA(4,3) / 0.2D1
  92.       AA(4,4) = 0.1D1 - t15 * JEA(4,4) / 0.2D1
  93.       AA(4,5) = -t19 * JEA(4,5) / 0.2D1
  94.       AA(4,6) = -t23 * JEA(4,6) / 0.2D1
  95.       AA(4,7) = -t27 * JEA(4,7) / 0.2D1
  96.       AA(4,8) = -t31 * JEA(4,8) / 0.2D1
  97.       AA(4,9) = -t35 * JEA(4,9) / 0.2D1
  98.       AA(4,10) = -t39 * JEA(4,10) / 0.2D1
  99.       AA(4,11) = -t43 * JEA(4,11) / 0.2D1
  100.       AA(4,12) = -t47 * JEA(4,12) / 0.2D1
  101.       AA(5,1) = -t2 * JEA(5,1) / 0.2D1
  102.       AA(5,2) = -t7 * JEA(5,2) / 0.2D1
  103.       AA(5,3) = -t11 * JEA(5,3) / 0.2D1
  104.       AA(5,4) = -t15 * JEA(5,4) / 0.2D1
  105.       AA(5,5) = 0.1D1 - t19 * JEA(5,5) / 0.2D1
  106.       AA(5,6) = -t23 * JEA(5,6) / 0.2D1
  107.       AA(5,7) = -t27 * JEA(5,7) / 0.2D1
  108.       AA(5,8) = -t31 * JEA(5,8) / 0.2D1
  109.       AA(5,9) = -t35 * JEA(5,9) / 0.2D1
  110.       AA(5,10) = -t39 * JEA(5,10) / 0.2D1
  111.       AA(5,11) = -t43 * JEA(5,11) / 0.2D1
  112.       AA(5,12) = -t47 * JEA(5,12) / 0.2D1
  113.       AA(6,1) = -t2 * JEA(6,1) / 0.2D1
  114.       AA(6,2) = -t7 * JEA(6,2) / 0.2D1
  115.       AA(6,3) = -t11 * JEA(6,3) / 0.2D1
  116.       AA(6,4) = -t15 * JEA(6,4) / 0.2D1
  117.       AA(6,5) = -t19 * JEA(6,5) / 0.2D1
  118.       AA(6,6) = 0.1D1 - t23 * JEA(6,6) / 0.2D1
  119.       AA(6,7) = -t27 * JEA(6,7) / 0.2D1
  120.       AA(6,8) = -t31 * JEA(6,8) / 0.2D1
  121.       AA(6,9) = -t35 * JEA(6,9) / 0.2D1
  122.       AA(6,10) = -t39 * JEA(6,10) / 0.2D1
  123.       AA(6,11) = -t43 * JEA(6,11) / 0.2D1
  124.       AA(6,12) = -t47 * JEA(6,12) / 0.2D1
  125.       AA(7,1) = -t2 * JEA(7,1) / 0.2D1
  126.       AA(7,2) = -t7 * JEA(7,2) / 0.2D1
  127.       AA(7,3) = -t11 * JEA(7,3) / 0.2D1
  128.       AA(7,4) = -t15 * JEA(7,4) / 0.2D1
  129.       AA(7,5) = -t19 * JEA(7,5) / 0.2D1
  130.       AA(7,6) = -t23 * JEA(7,6) / 0.2D1
  131.       AA(7,7) = 0.1D1 - t27 * JEA(7,7) / 0.2D1
  132.       AA(7,8) = -t31 * JEA(7,8) / 0.2D1
  133.       AA(7,9) = -t35 * JEA(7,9) / 0.2D1
  134.       AA(7,10) = -t39 * JEA(7,10) / 0.2D1
  135.       AA(7,11) = -t43 * JEA(7,11) / 0.2D1
  136.       AA(7,12) = -t47 * JEA(7,12) / 0.2D1
  137.       AA(8,1) = -t2 * JEA(8,1) / 0.2D1
  138.       AA(8,2) = -t7 * JEA(8,2) / 0.2D1
  139.       AA(8,3) = -t11 * JEA(8,3) / 0.2D1
  140.       AA(8,4) = -t15 * JEA(8,4) / 0.2D1
  141.       AA(8,5) = -t19 * JEA(8,5) / 0.2D1
  142.       AA(8,6) = -t23 * JEA(8,6) / 0.2D1
  143.       AA(8,7) = -t27 * JEA(8,7) / 0.2D1
  144.       AA(8,8) = 0.1D1 - t31 * JEA(8,8) / 0.2D1
  145.       AA(8,9) = -t35 * JEA(8,9) / 0.2D1
  146.       AA(8,10) = -t39 * JEA(8,10) / 0.2D1
  147.       AA(8,11) = -t43 * JEA(8,11) / 0.2D1
  148.       AA(8,12) = -t47 * JEA(8,12) / 0.2D1
  149.       AA(9,1) = -t2 * JEA(9,1) / 0.2D1
  150.       AA(9,2) = -t7 * JEA(9,2) / 0.2D1
  151.       AA(9,3) = -t11 * JEA(9,3) / 0.2D1
  152.       AA(9,4) = -t15 * JEA(9,4) / 0.2D1
  153.       AA(9,5) = -t19 * JEA(9,5) / 0.2D1
  154.       AA(9,6) = -t23 * JEA(9,6) / 0.2D1
  155.       AA(9,7) = -t27 * JEA(9,7) / 0.2D1
  156.       AA(9,8) = -t31 * JEA(9,8) / 0.2D1
  157.       AA(9,9) = 0.1D1 - t35 * JEA(9,9) / 0.2D1
  158.       AA(9,10) = -t39 * JEA(9,10) / 0.2D1
  159.       AA(9,11) = -t43 * JEA(9,11) / 0.2D1
  160.       AA(9,12) = -t47 * JEA(9,12) / 0.2D1
  161.       AA(10,1) = -t2 * JEA(10,1) / 0.2D1
  162.       AA(10,2) = -t7 * JEA(10,2) / 0.2D1
  163.       AA(10,3) = -t11 * JEA(10,3) / 0.2D1
  164.       AA(10,4) = -t15 * JEA(10,4) / 0.2D1
  165.       AA(10,5) = -t19 * JEA(10,5) / 0.2D1
  166.       AA(10,6) = -t23 * JEA(10,6) / 0.2D1
  167.       AA(10,7) = -t27 * JEA(10,7) / 0.2D1
  168.       AA(10,8) = -t31 * JEA(10,8) / 0.2D1
  169.       AA(10,9) = -t35 * JEA(10,9) / 0.2D1
  170.       AA(10,10) = 0.1D1 - t39 * JEA(10,10) / 0.2D1
  171.       AA(10,11) = -t43 * JEA(10,11) / 0.2D1
  172.       AA(10,12) = -t47 * JEA(10,12) / 0.2D1
  173.       AA(11,1) = -t2 * JEA(11,1) / 0.2D1
  174.       AA(11,2) = -t7 * JEA(11,2) / 0.2D1
  175.       AA(11,3) = -t11 * JEA(11,3) / 0.2D1
  176.       AA(11,4) = -t15 * JEA(11,4) / 0.2D1
  177.       AA(11,5) = -t19 * JEA(11,5) / 0.2D1
  178.       AA(11,6) = -t23 * JEA(11,6) / 0.2D1
  179.       AA(11,7) = -t27 * JEA(11,7) / 0.2D1
  180.       AA(11,8) = -t31 * JEA(11,8) / 0.2D1
  181.       AA(11,9) = -t35 * JEA(11,9) / 0.2D1
  182.       AA(11,10) = -t39 * JEA(11,10) / 0.2D1
  183.       AA(11,11) = 0.1D1 - t43 * JEA(11,11) / 0.2D1
  184.       AA(11,12) = -t47 * JEA(11,12) / 0.2D1
  185.       AA(12,1) = -t2 * JEA(12,1) / 0.2D1
  186.       AA(12,2) = -t7 * JEA(12,2) / 0.2D1
  187.       AA(12,3) = -t11 * JEA(12,3) / 0.2D1
  188.       AA(12,4) = -t15 * JEA(12,4) / 0.2D1
  189.       AA(12,5) = -t19 * JEA(12,5) / 0.2D1
  190.       AA(12,6) = -t23 * JEA(12,6) / 0.2D1
  191.       AA(12,7) = -t27 * JEA(12,7) / 0.2D1
  192.       AA(12,8) = -t31 * JEA(12,8) / 0.2D1
  193.       AA(12,9) = -t35 * JEA(12,9) / 0.2D1
  194.       AA(12,10) = -t39 * JEA(12,10) / 0.2D1
  195.       AA(12,11) = -t43 * JEA(12,11) / 0.2D1
  196.       AA(12,12) = 0.1D1 - t47 * JEA(12,12) / 0.2D1
  197.  
  198.  
  199.  
  200.  
  201.       if(affiche) then
  202.         print*, 'Dans a1ti3d : matrice de couplage'
  203.         do i=1,12
  204.             write (*,'(12E10.3)') (AA(i,j),j=1,12)
  205.         end do      
  206.       end if
  207.  
  208.       return
  209.       end
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales