Télécharger la plaquette Cast3M


fr en

Télécharger eplr3d.eso

Retour à la liste

Numérotation des lignes : 

eplr3d

  1. C EPLR3D    SOURCE    PV090527  23/01/27    21:15:32     11574          
  2.       SUBROUTINE EPLR3D(G,K,s1,s2,s3,f1,f2,f3,Vdep,numcas)
  3. c     Ecoulement plastique Trois Rankine (A.Sellier 2021/04/22)
  4.       implicit real*8 (a-h,o-z)
  5.       implicit integer (i-n)
  6.  
  7.       real*8 G,K,s1,s2,s3,f1,f2,f3,Vdep(3)
  8.       integer numcas
  9.  
  10.       if (numcas.eq.1) then
  11. c       cas 1-3
  12.       Vdep(1) = 0.3D1 * f1 / (0.4D1 * G + 0.3D1 * K)
  13.       Vdep(2) = 0.d0
  14.       Vdep(3) = 0.d0 
  15.  
  16.       else if (numcas.eq.2) then
  17. c       cas 1-2
  18.       Vdep(1) = (0.4D1 * G * f1 + 0.2D1 * G * f2 + 0.3D1 * K * f1 - 0.3D
  19.      #1 * K * f2) / G / (G + 0.3D1 * K) / 0.4D1
  20.       Vdep(2) = (0.2D1 * G * f1 + 0.4D1 * G * f2 - 0.3D1 * K * f1 + 0.3D
  21.      #1 * K * f2) / G / (G + 0.3D1 * K) / 0.4D1
  22.       Vdep(3) = 0.d0
  23.  
  24.       else if(numcas.eq.3) then
  25. c       cas 1-1     
  26.       Vdep(1) = (0.2D1 * G * f1 + 0.2D1 * G * f2 + 0.2D1 * G * f3 + 0.6D
  27.      #1 * K * f1 - 0.3D1 * K * f2 - 0.3D1 * K * f3) / G / K / 0.18D2
  28.       Vdep(2) = (0.2D1 * G * f1 + 0.2D1 * G * f2 + 0.2D1 * G * f3 - 0.3D
  29.      #1 * K * f1 + 0.6D1 * K * f2 - 0.3D1 * K * f3) / G / K / 0.18D2
  30.       Vdep(3) = (0.2D1 * G * f1 + 0.2D1 * G * f2 + 0.2D1 * G * f3 - 0.3D
  31.      #1 * K * f1 - 0.3D1 * K * f2 + 0.6D1 * K * f3) / G / K / 0.18D2
  32.  
  33.       else
  34.         print*,'incoherence dans eplr3d cas:',numcas
  35.         ierr1=1
  36.       end if
  37.  
  38.  
  39.       return
  40.       end
  41.  
  42.  
  43.  

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