Télécharger la plaquette Cast3M


fr en

Télécharger sigp3d.eso

Retour à la liste

Numérotation des lignes : 

sigp3d

  1. C SIGP3D    SOURCE    FD218221  24/02/07    21:15:26     11834          
  2.        subroutine sigp3d(nstype,var,nvari,vnvari,ndtens,
  3.      # iso,raideur66,sig3,V33)
  4.  
  5. c      calcul des contraintes principale a partir des deformations 
  6. c      elastiques stockees dans les vari
  7.  
  8. c      declarations externes
  9.        implicit integer(i-n)
  10.        implicit real*8(a-h,o-z)
  11.  
  12.        integer nstype,nvari,ndtens
  13.        logical iso
  14.        integer vnvari(nstype,2)
  15.        real*8 var(nvari)
  16.        real*8 raideur66(6,6)
  17.        real*8 sig3(3),V33(3,3)
  18.  
  19.        real*8 E3(3),E33(3,3),E6(6)
  20.        real*8 sig6(6),sig33(3,3)
  21.        integer nprec,i,j,k
  22.  
  23. c       print*,'ds sigp3d',nstype,nvari,vnvari,ndtens,iso
  24.  
  25.  
  26. c      recuperation des deformations elastique
  27.        call tvar3d(E3,V33,nprec,var,nvari,vnvari,nstype,ndtens,1,1)
  28.  
  29.        if(iso) then
  30. c           cas des materiaux elastiques isotropes en base principale
  31.             do i=1,3
  32.                 sig3(i)=0.d0
  33.                 do j=1,3
  34.                     sig3(i)=sig3(i)+raideur66(i,j)*E3(j)
  35.                 end do
  36.             end do
  37.        else
  38. c           cas d une elasticite non isotrope en base fixe
  39. c           passage de E en base fixe
  40.             call x12x6(E3,V33,E6)
  41. c           calcul du pseudo vecteur des contraintes
  42.             do i=1,6
  43.                 sig6(i)=0.d0
  44.                 do j=1,6
  45.                     sig6(i)=sig6(i)+raideur66(i,j)*E6(j)
  46.                 end do
  47.             end do
  48. c           valeur principales des contraintes et de la matrice de passage
  49.             call x6x33(sig6,sig33)
  50.             call b3_v33(sig33,sig3,V33)
  51.        end if
  52.        return
  53.        end
  54.  
  55.  

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