Télécharger HCOMFIB.INC

Retour à la liste

Numérotation des lignes :

  1. c debut romain calcul des contraintes endommages dans les fibres
  2. c reconstruction du tenseur des contraintesde fibres avec endo dans la base fixe
  3. if(ifibre.eq.1) then
  4. c deformation elastique
  5. call x6x33(epstf6,epse133)
  6. c base des fissures
  7. do i=1,3
  8. vecw(1)=vecw33(1,i)
  9. vecw(2)=vecw33(2,i)
  10. vecw(3)=vecw33(3,i)
  11. c deformation normale direction des fissures
  12. epse13(i)=sum(matmul(vecw,epse133)*vecw)
  13. c contribution des fibres zone non endommagee
  14. sigfe(i)=(1.d0-dt3(i))*rhof*
  15. # max(min(Ef*epse13(i)*cos(phimoy(i)*pi/180.d0)**2,
  16. # fuf*cos(phimoy(i)*pi/180.d0)**2),-fuf)
  17. end do
  18. c combinaison avec les fibres dans les fissures
  19. do i=1,3
  20. do j=1,3
  21. if(i.eq.j) then
  22. sigfp33(i,j)=sigfr(i)*dt3(i)+sigfe(i)
  23. else
  24. sigfp33(i,j)=0.d0
  25. end if
  26. end do
  27. end do
  28. sigfx33=matmul(matmul(vecw33,sigfp33),transpose(vecw33))
  29. call x33x6(sigfx33,sfib)
  30. c print*,wplx3(3),"dt3(3) fluendo",dt3(3)
  31. end if
  32. c fin romain fibres

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