Télécharger advar.eso

Retour à la liste

Numérotation des lignes :

  1. C ADVAR SOURCE PV 07/11/23 21:15:16 5978
  2. C ADVAR SOURCE ROUX G-M 03/02/24 21:15:01 4587
  3. SUBROUTINE ADVAR(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,DSPT,EPSVPT,
  4. & DEPST,VARPT,XMAT,NSTRS,NVARI,IFOUR,INPLAS,MFR)
  5.  
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. DIMENSION SIG(*),EPSV(*),VAR(*),SIG1(*),EPSV1(*),VAR1(*)
  8. DIMENSION DSPT(*),EPSVPT(*),VARPT(*),XMAT(*),DEPST(*)
  9. real*8 lamb,mu
  10. DATA UN/1.D0/
  11. DATA DEUX/2.D0/
  12.  
  13. C*******************************************************
  14. C ENTREE DE ADVAC
  15. C*******************************************************
  16.  
  17.  
  18. C----------------------------------------
  19. C mise à jour des déformations plastiques
  20. C----------------------------------------
  21. YOU=XMAT(1)
  22. XNU=XMAT(2)
  23.  
  24. DO 60 I=1,NSTRS
  25. EPSV1(I)=EPSV(I)+TAU*EPSVPT(I)
  26. 60 CONTINUE
  27.  
  28. C----------------------------
  29. C mise à jour des contraintes
  30. C----------------------------
  31. LAMB = YOU*XNU/((UN+XNU)*(UN-DEUX*XNU))
  32. MU = YOU/(DEUX*(UN+XNU))
  33. TRACE=EPSVPT(1)+EPSVPT(2)+EPSVPT(3)
  34. DO 70 I=1,3
  35. SIG1(I) = SIG(I) + TAU*(DSPT(I)-LAMB*TRACE-DEUX*MU*EPSVPT(I))
  36. 70 CONTINUE
  37. DO 80 I=4,NSTRS
  38. SIG1(I) = SIG(I) + TAU*(DSPT(I)-DEUX*MU*EPSVPT(I))
  39. 80 CONTINUE
  40.  
  41. C-----------------------------------
  42. C mise à jour des variables internes
  43. C-----------------------------------
  44. DO 90 I=1,NVARI
  45. VAR1(I)=VAR(I)+TAU*VARPT(I)
  46. VAR1(I)=VAR(I)+TAU*VARPT(I)
  47. 90 CONTINUE
  48.  
  49. RETURN
  50. END
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  

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