Numérotation des lignes :

eplc3d
C EPLC3D    SOURCE    PV090527  23/02/13    21:15:06     11592                SUBROUTINE EPLC3D(pt0,pc0,M,s1,s2,s3,lcc,Vdep)c     Ecoulement plastique plastique CamCaly (A.Sellier 2021/04/22)      implicit real*8 (a-h,o-z)      implicit integer (i-n)       real*8 pt0,pc0,M,s1,s2,s3,lcc,Vdep(3) C       t1 = M ** 2C       t2 = pc0 ** 2C       t5 = t2 * pc0 * t1 / 0.12D2C       t6 = pt0 * t1C       t7 = 0.9D1 * t6C       t8 = s1 + s2 + s3C       t9 = 0.2D1 * t1 * t8C       t11 = 0.3D1 * s2C       t12 = 0.3D1 * s3C       t16 = 0.4D1 / 0.9D1 * t1 * t8C       t18 = 0.2D1 / 0.3D1 * s2C       t19 = 0.2D1 / 0.3D1 * s3C       t24 = pt0 ** 2C       t25 = 0.2D1 / 0.3D1 * t1 * t8C       t32 = 0.3D1 * s1C       t37 = 0.2D1 / 0.3D1 * s1C       Vdep(1) = (t5 + t2 * (t7 + t9 + 0.6D1 * s1 - t11 - t12) / 0.36D2 +C      # pc0 * (t6 + t16 - 0.4D1 / 0.3D1 * s1 + t18 + t19) * pt0 / 0.4D1 +C      # (t6 + t25 + 0.2D1 * s1 - s2 - s3) * t24 / 0.12D2) * lccC       Vdep(2) = (t5 + t2 * (t7 + t9 - t32 + 0.6D1 * s2 - t12) / 0.36D2 +C      # pc0 * (t6 + t16 + t37 - 0.4D1 / 0.3D1 * s2 + t19) * pt0 / 0.4D1 +C      # t24 * (t6 + t25 - s1 + 0.2D1 * s2 - s3) / 0.12D2) * lccC       Vdep(3) = (t5 + t2 * (t7 + t9 - t32 - t11 + 0.6D1 * s3) / 0.36D2 +C      # pc0 * (t6 + t16 + t37 + t18 - 0.4D1 / 0.3D1 * s3) * pt0 / 0.4D1 +C      # (t6 + t25 - s1 - s2 + 0.2D1 * s3) * t24 / 0.12D2) * lcc       t1 = sqrt(0.3D1)      t3 = sqrt(0.4D1)      t4 = t3 * t1 * lcc      t6 = (pc0 - pt0) ** 2      t7 = t6 ** 2      t10 = s1 + s2 + s3 + 0.3D1 / 0.2D1 * pc0 + 0.3D1 / 0.2D1 * pt0      t11 = t10 ** 2      t12 = M ** 2      t13 = t12 ** 2      t15 = s1 ** 2      t19 = s2 ** 2      t23 = s3 ** 2      t27 = sqrt((t13 * t11 + 0.81D2 / 0.2D1 * t15 + 0.81D2 / 0.2D1 * (-     #s2 - s3) * s1 + 0.81D2 / 0.2D1 * t19 - 0.81D2 / 0.2D1 * s2 * s3 +      #0.81D2 / 0.2D1 * t23) * t7)      t28 = 0.1D1 / t27      t29 = t6 * t28      t30 = t12 * t10      t32 = 0.9D1 / 0.2D1 * s2      t33 = 0.9D1 / 0.2D1 * s3      t38 = 0.9D1 / 0.2D1 * s1      Vdep(1) = (t30 + 0.9D1 * s1 - t32 - t33) * t29 * t4 / 0.6D1      Vdep(2) = t6 * (t30 - t38 + 0.9D1 * s2 - t33) * t28 * t4 / 0.6D1      Vdep(3) = (t30 - t38 - t32 + 0.9D1 * s3) * t29 * t4 / 0.6D1              if(Vdep(1).ne.vdep(1)) then        print*,'Pb dans eplc3d'        print*,'pt0,pc0,M,s1,s2,s3,lcc,Vdep'        print*,pt0,pc0,M,s1,s2,s3,lcc,Vdep      end if        return      end   

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