Télécharger la plaquette Cast3M


fr en

Télécharger chre66.eso

Retour à la liste

Numérotation des lignes : 

chre66

  1. C CHRE66    SOURCE    FD218221  24/02/07    21:15:04     11834          
  2.       subroutine chre66(H066,H166,P033,TP033,P133,TP133)
  3.  
  4. c     A Sellier 2023 12 18  
  5.  
  6. c     P033 matrice de passage en base d orthotropie
  7. c     TPo33 son inverse
  8. c     Pb33 matrice de passage enbase quelconque 
  9. c     TPb33 son inverse  
  10.  
  11. c     changement de base des matrices de souplesse et de raideur
  12. c     attention la matrice H ne s adresse pas a des gamma pour permettre
  13. c     le meme passage pour C et H
  14.  
  15.       implicit real*8 (a-h,o-z)
  16.       implicit integer (i-n)
  17.  
  18. c     variables globales
  19.       real*8 H066(6,6),H166(6,6)
  20.       real*8 P033(3,3),TP033(3,3),P133(3,3),TP133(3,3)
  21.  
  22. c     variables locales
  23.       real*8 epsb(6),epsf(6),epso(6)
  24.       real*8 sigb(6),sigf(6),sigo(6)
  25.  
  26. c     methode par deformation virtuelle base d orthotropie
  27.  
  28. c     boucle sur les composantes virtuelle en base b
  29.       do i=1,6      
  30.          do j=1,6
  31.              if(j.eq.i) then
  32.                 epsb(j)=1.d0
  33.              else
  34.                 epsb(j)=0.d0
  35.              end if
  36.          end do
  37. c         print*,'epsb',epsb         
  38. c        passage en base fixe
  39.          call chrep6(epsb,TP133,.false.,epsf)
  40. c         print*,'epsf',epsf
  41. c        passage en base d orthotropie
  42.          call chrep6(epsf,P033,.false.,epso)
  43. c         print*,'epso',epso
  44. c        calcul des contrainte en base orthotropie
  45.          do j=1,6
  46.             sigo(j)=0.d0
  47.             do k=1,6
  48.                 sigo(j)=sigo(j)+H066(j,k)*epso(k)
  49.             end do
  50.          end do
  51. c         print*,'sigo',sigo
  52. c        retour des contraintes en base fixe
  53.          call chrep6(sigo,TP033,.false.,sigf)
  54. c         print*,'sigf',sigf
  55. c        retour des contraintes en base b         
  56.          call chrep6(sigf,P133,.false.,sigb)
  57. c         print*,'sigb',sigb         
  58. c        stockage de la colonne i de la matrice de rigidite en base b
  59. c         print*,'colonne', i
  60.          do j=1,6
  61.             H166(j,i)=sigb(j)
  62. c            print*,'ligne',j,H166(j,i)
  63.          end do
  64.       end do
  65. c     matrice H en base b
  66. c      call afic66(H166)  
  67. c      read*      
  68.  
  69.       return
  70.       end
  71.  
  72.  
  73.  
  74.  

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