Télécharger la plaquette Cast3M


fr en

Télécharger MATFIB.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