Sources Esope | Cast3M

dedpta
C DEDPTA    SOURCE    CB215821  16/04/21    21:16:21     8920
CCC
C **********************************************************************
CCC
      SUBROUTINE DER_ENERELAS_DPRAL(X,SIG,NMODEL)
      IMPLICIT INTEGER(I-N)
      integer nmodel,i,j
      real*8  x(3),sig(3),e(3,3)
      real*8                  kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3
      logical                 ogden
      common /elasdata_ogden/ kap,xmu,xmu1,alp1,xmu2,alp2,xmu3,alp3,
     .                        ogden
      if (ogden) goto 10
c Modelo cuadratico
      call der2_enerelas_dpral(x,e,3,nmodel)
      do i=1,3
       sig(i)=0.D0
       do j=1,3
        sig(i)=sig(i)+e(i,j)*x(j)
       enddo
      enddo
c hacer desviador para VMT
      if (nmodel.eq.8) then
       call desviador(sig,e,3)
       call equv(e,sig,3)
      endif
      return
c Modelo ogden
10    continue
      do i=1,3
       sig(i)=EXP(x(i))
      enddo
      e(1,1)=xmu*(xmu1*sig(1)**alp1+xmu2*sig(1)**alp2+xmu3*sig(1)**alp3)
      e(2,1)=xmu*(xmu1*sig(2)**alp1+xmu2*sig(2)**alp2+xmu3*sig(2)**alp3)
      e(3,1)=xmu*(xmu1*sig(3)**alp1+xmu2*sig(3)**alp2+xmu3*sig(3)**alp3)
20    continue
      call desviador(e,sig,3)
      return
      end