Télécharger eco_j2.eso

Retour à la liste

Numérotation des lignes :

  1. C ECO_J2 SOURCE CHAT 05/01/12 23:25:03 5004
  2. CCC
  3. C **********************************************************************
  4. CCC
  5. SUBROUTINE ECO_J2 (SIGINI,VARINI,DSIGTR,SIGFIN,VARFIN,
  6. . ddefpl,iplast,ndims,xmat,kerre,
  7. . precis,nitmax,nescri,ues,iterlocal)
  8. IMPLICIT INTEGER(I-N)
  9. integer ndims,kerre,nitmax,nescri,ues,iterlocal,ione
  10. real*8 sigini(ndims),varini(2),ddefpl(ndims),dsigtr(ndims),
  11. . sigfin(ndims),varfin(2),xmat(*),precis
  12. integer iplast,ndimx
  13. real*8 sigtri(6),xtri(7),xini(7),lini,void,resu
  14. real*8 sigy0,kiso,siginf,velo
  15. common /j2j2data/ sigy0,kiso,siginf,velo
  16. void =0.D0
  17. ione =1
  18. kerre =0
  19. iplast=1
  20. do i=1,6
  21. sigtri(i)=0.D0
  22. enddo
  23. do i=1,7
  24. xtri(i)=0.D0
  25. xini(i)=0.D0
  26. enddo
  27. call carac_mate_j2(xmat)
  28. do i=1,ndims
  29. sigtri(i)=sigini(i)+dsigtr(i)
  30. enddo
  31. call yielddJ2(sigtri,ndims,varini,ione,resu)
  32. if (resu.le.0.D0) then
  33. if (nescri.eq.1) write(ues,*)' Paso elastico'
  34. iplast=0
  35. do i=1,ndims
  36. sigfin(i)=sigtri(i)
  37. enddo
  38. varfin(1)=varini(1)
  39. varfin(2)=-1.D0
  40. iterlocal=0
  41. return
  42. endif
  43. if (nescri.eq.1) write(ues,*)' Paso plastico'
  44. do i=1,ndims
  45. xtri(i)=sigtri(i)
  46. xini(i)=sigtri(i)
  47. enddo
  48. lini=0.D0
  49. if ((kiso.eq.0.D0).and.(velo.eq.0.D0)) then
  50. if (nescri.eq.1) write(ues,*)' Plasticidad sin endurecimiento'
  51. call Integra0J2(xtri,xini,ndims,lini,ddefpl,ndims,
  52. . precis,nitmax,nescri,ues,kerre,iterlocal)
  53. varfin(1)=varini(1)+lini
  54. varfin(2)=lini
  55. else
  56. if (nescri.eq.1) write(ues,*)' Plasticidad con endurecimiento'
  57. ndimx=ndims+1
  58. xtri(ndimx)=varini(1)
  59. xini(ndimx)=varini(1)
  60. call Integra1J2(xtri,xini,ndimx,ione,lini,ddefpl,ndims,
  61. . precis,nitmax,nescri,ues,kerre,iterlocal)
  62. varfin(1)=xini(ndimx)
  63. varfin(2)=lini
  64. endif
  65. do i=1,ndims
  66. sigfin(i)=xini(i)
  67. enddo
  68. return
  69. end
  70.  
  71.  
  72.  

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