Télécharger la plaquette Cast3M


fr en

Télécharger prtt3d.eso

Retour à la liste

Numérotation des lignes : 

prtt3d

  1. C PRTT3D    SOURCE    PV090527  23/01/27    21:15:59     11574          
  2. c*************************************************************************
  3.       subroutine prtt3d(sige6,sige3,vsige33,vsige33t,
  4.      #siget6,sigec6,sigec3,siget3)
  5.  
  6. c     declaration des varibles externes
  7.       real*8 sige6(6),sige3(3),vsige33(3,3),vsige33t(3,3)
  8.       real*8 siget6(6),sigec6(6),sigec3(3),siget3(3)
  9. c     declaration des varibles locales
  10.       real*8 sige33(3,3)
  11.       real*8 x33(3,3),siget33(3,3),sigec33(3,3),sige33p(3,3),sige6p(6)
  12.       integer i      
  13.  
  14. c     rangement des contraintes effectives en tableau 3*3
  15.       call x6x33(sige6,sige33)
  16. c     diagonalisation contraintes effectives actuelles et valeurs propres par la methode de jacobi
  17.       call b3_v33(sige33,sige3,vsige33)
  18. c     creation de la matrice de passage inverse    
  19.       call traps1(vsige33t,vsige33,3)
  20. c     decomposition des contraintes principales en partie positive et négative dans 
  21. c     la base principale (avec prise en compte des erreurs numeriques de diagonalisation)
  22. c     on suppose sige33p pas tout a fait diagonale par defaut et on utilise
  23. c     que les contraintes normales positves pour faire la partition
  24.       call chre3(sige33p,sige33,vsige33)
  25.       call x33x6(sige33p,sige6p)
  26.       do i=1,3
  27.        siget3(i)=0.5d0*(sige33p(i,i)+abs(sige33p(i,i)))
  28.        sigec3(i)=0.5d0*(sige33p(i,i)-abs(sige33p(i,i)))
  29.        siget6(i)=siget3(i)
  30.        sigec6(i)=sigec3(i)
  31.       end do
  32.       do i=4,6
  33.        siget6(i)=0.d0
  34.        sigec6(i)=sige6p(i)-siget6(i)
  35.       end do
  36. c     stockage des parties positives et negatives en base fixe
  37. c     cas des contraintes de traction
  38.       call x6x33(siget6,x33)
  39.       call chre3(siget33,x33,vsige33t)
  40.       call x33x6(siget33,siget6)
  41. c     cas des contraintes de compression
  42.       call x6x33(sigec6,x33)
  43.       call chre3(sigec33,x33,vsige33t)
  44.       call x33x6(sigec33,sigec6)
  45.       return
  46.       end
  47.  
  48. **************************************************************************************
  49.  
  50.  
  51.  

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