Sources Esope | Cast3M

chre36
C CHRE36    SOURCE    PV090527  23/01/27    21:15:12     11574          
      SUBROUTINE  CHRE36(sig3,sig6,v33)
c     passage base prin->base fixe
c     (A.Sellier 2021/04/24)
      implicit real*8 (a-h,o-z)
      implicit integer (i-n)      
 
c     passage base prin->base fixe
      real*8 sig3(3),sig6(6),v33(3,3)
 
      integer i,j,k,l
      real*8 v3(3),b33(3,3,3)
      real*8 sig33(3,3)
c     tenseur d orientation      
      do i=1,3
        do j=1,3
            v3(j)=v33(j,i)
        end do
        do l=1,3
            do k=1,3
                b33(i,l,k)=v3(l)*v3(k)
            end do
        end do
      end do
c     tenseur physique      
      do i=1,3
        do j=1,3
           sig33(i,j)=0.d0
        end do
      end do
      do i=1,3
        do l=1,3
            do k=l,3
                 sig33(l,k)=sig33(l,k)+sig3(i)*b33(i,l,k)
            end do
        end do
      end do
c     stockage pseudo vecteur      
      do i=1,3
        sig6(i)=sig33(i,i)
      end do
      sig6(4)=sig33(1,2)
      sig6(5)=sig33(1,3)
      sig6(6)=sig33(2,3)
 
      return
      end