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