Télécharger HMATFIB.INC

Retour à la liste

Numérotation des lignes :

  1. c romain
  2. c ************recup parametres pour fibre3d *******
  3. c indicateur presence de fibres ( romain )
  4. ifibre=int(xmat(NMAT1-2))
  5.  
  6.  
  7. if(ifibre.eq.1) then
  8.  
  9. c print*,"debut recuperation fibres romain",NBRFLU3D,
  10. c #xmat(nbelas3d+92),xmat(nbelas3d+93),xmat(nbelas3d+94)
  11. Df=xmat(NMAT1+1)
  12. Lf=xmat(NMAT1+2)
  13. Hf=xmat(NMAT1+3)
  14. Tmax=xmat(NMAT1+4)
  15. Td=xmat(NMAT1+5)
  16. sk=xmat(NMAT1+6)
  17. fabo=xmat(NMAT1+7)
  18. alphaec=xmat(NMAT1+8)
  19. m0f=xmat(NMAT1+9)
  20. L0f=xmat(NMAT1+10)
  21. muf=xmat(NMAT1+11)
  22. Ef=xmat(NMAT1+12)
  23. lech=xmat(NMAT1+13)
  24. mwf=xmat(NMAT1+14)
  25. fuf=xmat(NMAT1+15)
  26. rtec=xmat(NMAT1+16)
  27. rhof=xmat(NMAT1+17)
  28. eof(1)=xmat(NMAT1+18)
  29. eof(2)=xmat(NMAT1+19)
  30. eof(3)=xmat(NMAT1+20)
  31. vf1(1)=xmat(NMAT1+21)
  32. vf1(2)=xmat(NMAT1+22)
  33. vf1(3)=xmat(NMAT1+23)
  34. vf2(1)=xmat(NMAT1+24)
  35. vf2(2)=xmat(NMAT1+25)
  36. vf2(3)=xmat(NMAT1+26)
  37. fyf=xmat(NMAT1+27)
  38.  
  39. c modif de la resistance a la traction
  40. c Rt00=Rtec*(lech/longr(i))**(1/mw)
  41. c print*,NMAT1+26,nbelas3d,NBRFLU3D,nbparr
  42.  
  43. do i=1,10
  44. cokf1(i)=xmat(NMAT1+NB_PARA_FIBRE_U+i)
  45. cokf2(i)=xmat(NMAT1+NB_PARA_FIBRE_U+10+i)
  46. cowp1(i)=xmat(NMAT1+NB_PARA_FIBRE_U+20+i)
  47. cowp2(i)=xmat(NMAT1+NB_PARA_FIBRE_U+30+i)
  48. cofp(i)=xmat(NMAT1+NB_PARA_FIBRE_U+40+i)
  49. cow03(i)=xmat(NMAT1+NB_PARA_FIBRE_U+50+i)
  50. cowf(i)=xmat(NMAT1+NB_PARA_FIBRE_U+60+i)
  51. end do
  52.  
  53. mink1=xmat(NMAT1+NB_PARA_FIBRE_U+71)
  54. maxk1=xmat(NMAT1+NB_PARA_FIBRE_U+72)
  55. mink2=xmat(NMAT1+NB_PARA_FIBRE_U+73)
  56. maxk2=xmat(NMAT1+NB_PARA_FIBRE_U+74)
  57. minwp1=xmat(NMAT1+NB_PARA_FIBRE_U+75)
  58. maxwp1=xmat(NMAT1+NB_PARA_FIBRE_U+76)
  59. minwp2=xmat(NMAT1+NB_PARA_FIBRE_U+77)
  60. maxwp2=xmat(NMAT1+NB_PARA_FIBRE_U+78)
  61. minfp=xmat(NMAT1+NB_PARA_FIBRE_U+79)
  62. maxfp=xmat(NMAT1+NB_PARA_FIBRE_U+80)
  63. minw03=xmat(NMAT1+NB_PARA_FIBRE_U+81)
  64. maxw03=xmat(NMAT1+NB_PARA_FIBRE_U+82)
  65. minwf=xmat(NMAT1+NB_PARA_FIBRE_U+83)
  66. maxwf=xmat(NMAT1+NB_PARA_FIBRE_U+84)
  67. end if
  68.  
  69. c print*,"fin recuperation fibres romain"
  70. c fin recuperation fibres romain

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