Télécharger dedpta.eso

Retour à la liste

Numérotation des lignes :

  1. C DEDPTA SOURCE CB215821 16/04/21 21:16:21 8920
  2. CCC
  3. C **********************************************************************
  4. CCC
  5. SUBROUTINE DER_ENERELAS_DPRAL(X,SIG,NMODEL)
  6. IMPLICIT INTEGER(I-N)
  7. integer nmodel,i,j
  8. real*8 x(3),sig(3),e(3,3)
  9. real*8 kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3
  10. logical ogden
  11. common /elasdata_ogden/ kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3,
  12. . ogden
  13. if (ogden) goto 10
  14. c Modelo cuadratico
  15. call der2_enerelas_dpral(x,e,3,nmodel)
  16. do i=1,3
  17. sig(i)=0.D0
  18. do j=1,3
  19. sig(i)=sig(i)+e(i,j)*x(j)
  20. enddo
  21. enddo
  22. c hacer desviador para VMT
  23. if (nmodel.eq.8) then
  24. call desviador(sig,e,3)
  25. call equv(e,sig,3)
  26. endif
  27. return
  28. c Modelo ogden
  29. 10 continue
  30. do i=1,3
  31. sig(i)=EXP(x(i))
  32. enddo
  33. e(1,1)=xmu*(xmu1*sig(1)**alp1+xmu2*sig(1)**alp2+xmu3*sig(1)**alp3)
  34. e(2,1)=xmu*(xmu1*sig(2)**alp1+xmu2*sig(2)**alp2+xmu3*sig(2)**alp3)
  35. e(3,1)=xmu*(xmu1*sig(3)**alp1+xmu2*sig(3)**alp2+xmu3*sig(3)**alp3)
  36. 20 continue
  37. call desviador(e,sig,3)
  38. return
  39. end
  40.  
  41.  
  42.  
  43.  
  44.  
  45.  

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