Télécharger mtcmrs.eso

Retour à la liste

Numérotation des lignes :

  1. C MTCMRS SOURCE CHAT 05/01/13 01:53:51 5004
  2. CCC
  3. C **********************************************************************
  4. CCC
  5. SUBROUTINE MTC_MRSMAC (KMAT,NDIMK,SIG,NDIMS,VINTCON,VINTCAP,LCON,
  6. . lcap,xmat,nescri,ues,nnumer,deltax,kerre)
  7. IMPLICIT INTEGER(I-N)
  8. integer ndims,kerre,nescri,nnumer,ues,ndimk,
  9. . i,ifl,ndimx,kmod,ndimv
  10. real*8 sig(ndims),vintcon,vintcap,Kmat(*),xmat(*),deltax,x(7),
  11. . lcon,lcap,lam
  12. real*8 exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,exl,xnn,aaa,
  13. . bbb,pc0,exr,cca,phi,alp,voi1,voi2,voi3,voi4,voi5,cpr,
  14. . cpm,cfr
  15. common /MRScone/ exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
  16. . exl,xnn,aaa,bbb
  17. common /MRScapp/ pc0,exr,cca,phi,alp,voi1,voi2,voi3,voi4,voi5,
  18. . cpr,cpm,cfr
  19. do i=1,7
  20. x(i)=0.D0
  21. enddo
  22. kerre=0
  23. call CaracmateMAC(xmat)
  24. if ((lcon.le.0.D0).and.(lcap.le.0.D0)) then
  25. if (nescri.eq.1) write(ues,*)' Modelo MAC: paso elastico'
  26. ifl=1
  27. call MatHok(Kmat,ndimk,ndims,ifl)
  28. goto 100
  29. endif
  30. if ((lcon.gt.0.D0).and.(lcap.gt.0.D0)) then
  31. write(*,*)' Modelo MAC: no puede corner'
  32. return
  33. endif
  34. do i=1,ndims
  35. x(i)=sig(i)
  36. enddo
  37. if ((lcon.gt.0.D0).and.(lcap.lt.0.D0)) then
  38. kmod=21
  39. lam=lcon
  40. else
  41. kmod=22
  42. lam=lcap
  43. endif
  44. if ((exr.lt.0.D0).and.(exv.lt.0.D0)) then
  45. if (nescri.eq.1) write(ues,*)' Modelo MAC: plast perfecta'
  46. write(*,*)' Formulacion no disponible'
  47. kerre=1
  48. c ndimv=0
  49. c call MTC0(Kmat,ndimk,sig,ndims,ndimv,lam,
  50. c . kmod,nescri,ues,nnumer,deltax,kerre)
  51. else if ((exr.lt.0.D0).and.(exv.gt.0.D0)) then
  52. if (nescri.eq.1) write(ues,*)' Modelo MAC: cone hardening'
  53. write(*,*)' Formulacion no disponible'
  54. kerre=1
  55. else if ((exr.gt.0.D0).and.(exv.lt.0.D0)) then
  56. if (nescri.eq.1) write(ues,*)' Modelo MAC: cap hardening '
  57. write(*,*)' Formulacion no disponible'
  58. kerre=1
  59. else if ((exr.gt.0.D0).and.(exv.gt.0.D0)) then
  60. if (nescri.eq.1) write(ues,*)' Modelo MAC: cone-cap harde'
  61. x(ndims+1)=vintcon
  62. x(ndims+2)=vintcap
  63. ndimx=ndims+2
  64. call MTC21(Kmat,ndimk,x,ndimx,lam,ndims,
  65. . kmod,nescri,ues,nnumer,deltax,kerre)
  66. endif
  67. 100 continue
  68. return
  69. end
  70.  
  71.  
  72.  
  73.  

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