Télécharger dplas3d.eso

Retour à la liste

Numérotation des lignes :

  1. C DPLAS3D SOURCE AM 12/02/27 23:39:38 7289
  2. SUBROUTINE dplas3D(x1,x2,P2,pi2,ag,fc,fb,
  3. . parahot3,idimpara3,dplas,lcp)
  4.  
  5. c This subroutine calculates the INCREMENT OF PLASTIC STRAIN corresponding to
  6. c the plastic multipliers Dlambdat = x1 and Dlambdac = x2
  7.  
  8. IMPLICIT REAL*8 (A-B,D-H,O-Z)
  9. implicit integer (I-K,M,N)
  10. implicit logical (L)
  11. implicit character*10 (C)
  12.  
  13. dimension sigmai(6)
  14. dimension pi2(6)
  15. dimension dgc(6),dgt(6),dplas(6),P2(6,6)
  16. dimension parahot3(idimpara3)
  17.  
  18. i0 = 0
  19. i1 = 1
  20. i2 = 2
  21. i3 = 3
  22. i4 = 4
  23. i5 = 5
  24. i6 = 6
  25. i7 = 7
  26.  
  27. r2 = 2.
  28.  
  29. Dlambdat = x1
  30. Dlambdac = x2
  31.  
  32. c take the effective stress at the previous converged step
  33. do iloc=i1,i6
  34. sigmai(iloc) = parahot3(idimpara3-48+iloc)
  35. end do
  36.  
  37. call dsiggc(sigmai,dgc,ag,P2,pi2,i1,i6,fc,fb,lcp)
  38. call dsigRank(sigmai,dgt,i3,i6,parahot3,idimpara3,rkappait,lcp)
  39.  
  40. do iloc=i1,i6
  41. dplas(iloc) = Dlambdat*dgt(iloc) + Dlambdac*dgc(iloc)
  42. end do
  43.  
  44. RETURN
  45. END
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  

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