Télécharger la plaquette Cast3M


fr en

Télécharger dcmp3d.eso

Retour à la liste

Numérotation des lignes : 

dcmp3d

  1. C DCMP3D    SOURCE    PV090527  23/01/27    21:15:21     11574          
  2.         subroutine dcmp3d(iphase,itps,NBRINC,x6,x3,xc6,
  3.      #   v33,v33t,affiche,nomt)
  4.  
  5. c       recuperation des projections normale d un tenseur sur les axes
  6. c       de la base v33 ->x3
  7. c       construction du tenseur complementaire dans la base fixe->xc6
  8.  
  9. c       tables de dimension fixe pour resolution des systemes lineaires 
  10.         implicit real*8 (a-h,o-z)
  11.         implicit integer (i-n)
  12.  
  13.         logical affiche
  14.         character*5 nomt
  15.  
  16.         integer NBRINC,iphase,itps
  17.  
  18.         real*8 x6(-1:NBRINC,6,0:1)
  19.         real*8 x3(-1:NBRINC,3),xc6(-1:NBRINC,6)
  20.         real*8 v33(3,3),v33t(3,3)
  21.  
  22.         real*8 x06(6),x03(3),x16(6)
  23.         integer  j
  24.  
  25.         if((iphase.lt.-1).or.(iphase.gt.NBRINC)) then
  26.           print*,'Pb dans incl3d'
  27.           print*,'dimension NBRINC dans dcmp3d'
  28.           stop
  29.         end if
  30.  
  31.         do j=1,6
  32.             x06(j)=x6(iphase,j,itps)
  33.         end do
  34.         call prjc3d(x06,v33,v33t,x03,x16,.false.)
  35.         do j=1,6
  36.              if(j.le.3) then
  37. c               recuperation des projections            
  38.                 x3(iphase,j)=x03(j)
  39.              end if 
  40. c           recuperation de la partie complémentaire            
  41.             xc6(iphase,j)=x16(j)
  42.         end do
  43.  
  44.  
  45.         if(affiche) then
  46.               print*,'Dans Decomp3D'
  47.               write(*,'(A26,1X,A5,1X,A12,I2)')
  48.      #        'Decomposition du tenseur :',nomt,'de la phase:',iphase 
  49.               do icomp=1,6
  50.                  write(*,'(A5,A1,I2,A2,E10.3,1X,A5,A6,I2,A2,E10.3 )')
  51.      #           nomt,'(',icomp,')=',x6(iphase,icomp,0),
  52.      #           nomt,'_Comp(',icomp,')=',xc6(iphase,icomp)
  53.               end do     
  54.               do icomp=1,3
  55.                     write(*,'(A5,A6,I2,A2,E10.3)') 
  56.      #              nomt,'_PRIN(',icomp,')=',x3(iphase,icomp)
  57.               end do
  58. c              print*,'Valider pour continuer'
  59. c              read*
  60.          end if 
  61.  
  62.         return
  63.         end
  64.  
  65.  
  66.  

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