Esope source files

   1 : C DEDPTA    SOURCE    CHAT      05/01/12    22:39:36     5004           
   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 
  14 :       call der2_enerelas_dpral(x,e,3,nmodel)
  15 :       do i=1,3
  16 :        sig(i)=0.D0
  17 :        do j=1,3
  18 :         sig(i)=sig(i)+e(i,j)*x(j)
  19 :        enddo
  20 :       enddo
  21 : c hacer desviador para VMT
  22 :       if (nmodel.eq.8) then   
  22 :        call desviador(sig,e,3)
  22 :        call equv(e,sig,3)
  23 :       endif
  24 :       return
  25 : c Modelo ogden
  26 : 10    continue
  27 :       do i=1,3
  28 :        sig(i)=dexp(x(i))
  29 :       enddo
  30 :       e(1,1)=xmu*(xmu1*sig(1)**alp1+xmu2*sig(1)**alp2+xmu3*sig(1)**alp3)
  31 :       e(2,1)=xmu*(xmu1*sig(2)**alp1+xmu2*sig(2)**alp2+xmu3*sig(2)**alp3)
  32 :       e(3,1)=xmu*(xmu1*sig(3)**alp1+xmu2*sig(3)**alp2+xmu3*sig(3)**alp3)
  33 : 20    continue
  33 :       call desviador(e,sig,3)
  34 :       return
  35 :       end
  36 :  
  37 :  
  38 :