Esope source files

   1 : C CMATEM    SOURCE    CHAT      05/01/12    22:09:13     5004           
   2 : c MODELO MRSLade-Macari
   3 : CCC
   4 : C **********************************************************************
   5 : CCC
   6 :       SUBROUTINE CARACMATEMAC(XMAT)
   7 :       IMPLICIT INTEGER(I-N)
   8 :       real*8    xmat(*)
   9 :       real*8    you,xnu,exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
  10 :      .          exl,xnn,aaa,bbb,pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,
  11 :      .          exm2,eww2,cpr,cpm,cfr
  12 :       common /elasdata/ you,xnu
  13 :       common /MRScone/ exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
  14 :      .                 exl,xnn,aaa,bbb
  15 :       common /MRScapp/ pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,exm2,eww2,
  16 :      .                 cpr,cpm,cfr
  17 :       you = xmat(1)
  18 :       xnu = xmat(2)
  19 :       pcc = xmat(5)
  20 :       paa = xmat(6)
  21 :       qaa = xmat(7)
  22 :       exm = xmat(8)
  23 :       eww = xmat(9)
  24 :       xk1 = xmat(10)
  25 :       xk2 = xmat(11)
  26 :       etb = xmat(12)
  27 :       exv = xmat(13)
  28 :       xep = xmat(14)
  29 :       xnn = xmat(15)
  30 :       cco = xmat(16)
  31 :       exl = xmat(17)
  32 :       pc0 = xmat(18)
  33 :       exr = xmat(19)
  34 :       cca = xmat(20)
  35 :       phi = xmat(21)
  36 :       alp = xmat(22)
  37 : c exv < 0 cone fijo
  38 :       if (exv.le.0.D0) then
  39 :        bbb=1.D0/(abs(exv)*(1.D0+xk1))+xk2*xep/(1.D0+xep)/(1.D0+xep-xk2)
  40 :        aaa=(1.D0+xk1)**(-1.D0/abs(exv))*
  41 :      .     (1.D0-xk2/(1.D0+xep))*etb*dexp(bbb)
  42 :       else
  43 :        bbb=1.D0/(exv*(1.D0+xk1))+xk2*xep/(1.D0+xep)/(1.D0+xep-xk2)
  44 :        aaa=(1.D0+xk1)**(-1.D0/exv)*(1.D0-xk2/(1.D0+xep))*etb*dexp(bbb)
  45 :       endif
  46 :       cpr=(1.D0-alp)*(phi*(1.D0-alp)+alp)/(2.D0*phi*(1.D0-alp)+alp)
  47 :       cpm=(alp*alp+phi*(1.D0-alp*alp))/(2.D0*phi*(1.D0-alp)+alp)
  48 :       cfr=(phi*(1.D0-alp)+alp)*dsqrt(alp/(2.D0*phi*(1.D0-alp)+alp))
  49 :       pcc2=pcc
  50 :       paa2=paa
  51 :       qaa2=qaa
  52 :       exm2=exm
  53 :       eww2=eww
  54 :       return
  55 :       end
  56 :  
  57 :  
  58 :