Sources Esope | Cast3M

cons3d
C CONS3D    SOURCE    FD218221  24/02/07    21:15:07     11834          
       subroutine cons3d(xmat,nmat,vnmat,nstype,var0,varf,nvari,
     # vnvari,ndtens,nbelas,xxk,V33,CC3)
 
c      calcul des coeff de consolidation pour mc3d
 
c      declarations externes
       implicit integer(i-n)
       implicit real*8(a-h,o-z)
 
       integer nmat,nstype,ndtens,nvari,nbelas
       integer vnmat(nstype)
       integer vnvari(nstype,2)
       real*8 xmat(nmat)
       real*8 var0(nvari),varf(nvari)
       real*8 xxk,CC3(3),V33(3,3)
 
       real*8 E3(3),E33(3,3),E6(6),EV6(6),M3(3),M33(3,3),M6(6),MV6(6)
       integer nprec,i
       real*8 xx1,xx2
       real*8 epse_min
       parameter (epse_min=1.0d-6)
 
c      recuperation des deformations elastiques actualisees par l hydratation
       call tvar3d(E3,E33,nprec,varf,nvari,vnvari,nstype,ndtens,1,1)
c      reconstruction du tenseur des deformations de Maxwell en base fixe
       call x12x6(E3,E33,E6)
c      passage de la deformation elastique dans la base souhaitee
       call chrep6(E6,V33,.false.,EV6)
c       print*,'ds cons3d ev6',ev6
 
c      recuperation des deformations de Maxwell actualisees par l hydratation
       call tvar3d(M3,M33,nprec,varf,nvari,vnvari,nstype,ndtens,2,1)
c      reconstruction du tenseur des deformations de Maxwell en base fixe
       call x12x6(M3,M33,M6)
c      passage de la deformation de Maxwell dans la base souhaitee
       call chrep6(M6,V33,.false.,MV6)
c       print*,'ds cons3d mv6',MV6
 
       do i=1,3
        xx1=max(abs(EV6(i)),epse_min)         
        xx2=xxk*xx1        
        CC3(i)=(exp(min(abs(MV6(i))/xx2,25.d0)))/xxk
c        print*,'ds cons3 cc',i,CC3(i)
       end do
       return
       end