Télécharger la plaquette Cast3M


fr en

Télécharger cmp1de.eso

Retour à la liste

Numérotation des lignes : 

cmp1de

  1. C CMP1DE    SOURCE    CHAT      05/01/12    22:11:00     5004
  2. CCC
  3. C **********************************************************************
  4. CCC
  5.       SUBROUTINE CARAC_MATE_POWDER1_DENSI(XMAT,XDENSI)
  6.       IMPLICIT INTEGER(I-N)
  7.       real*8    xmat(*),xdensi,eta
  8.       real*8            you,xnu,sigy0,siginf,kiso,velo,cpar,mpar
  9.       common /elasdata/ you,xnu
  10.       common /miehdata/         sigy0,kiso,siginf,velo,cpar,mpar
  11.       real*8                  kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3
  12.       logical                 ogden
  13.       common /elasdata_ogden/ kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3,
  14.      .                        ogden
  15.       real*8               nn1,nn2,eta0,aa1,aa2
  16.       common /poder1_data/ nn1,nn2,eta0,aa1,aa2
  17.  
  18.       ogden  =.false.
  19.       you    = xmat( 1)
  20.       xnu    = xmat( 2)
  21.       sigy0  = xmat( 5)
  22.       nn1    = xmat( 6)
  23.       nn2    = xmat( 7)
  24.       eta0   = xmat( 8)
  25.       eta    = xdensi*eta0
  26.       aa1    = (1.D0 - eta**2)/(2.D0 + eta**2)
  27.       if (aa1.lt.1.D-30) aa1=1.D-30
  28.       aa1    = aa1**nn1
  29.       if (aa1.lt.1.D-5) aa1=1.D-5
  30.       aa2    = abs(eta - 0.98*eta0)/(1.D0 - 0.98*eta0)
  31.       if (eta.lt.eta0)  aa2=0.02D0*eta0/(1.D0 - 0.98D0*eta0)
  32.       aa2    = aa2**nn2
  33.  
  34.       siginf = 0.D0
  35.       kiso   = 0.D0
  36.       velo   = 0.D0
  37.       cpar   = 0.D0
  38.       mpar   = 1.D0
  39.       return
  40.       end
  41.  
  42.  
  43.  
  44.  

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