Télécharger la plaquette Cast3M


fr en

Télécharger ecomrs.eso

Retour à la liste

Numérotation des lignes : 

ecomrs

  1. C ECOMRS    SOURCE    CHAT      05/01/12    23:25:07     5004
  2. CCC
  3. C **********************************************************************
  4. CCC
  5.       SUBROUTINE ECO_MRSMAC (SIGINI,VARINI,DSIGTR,SIGFIN,VARFIN,
  6.      .                      ddefpl,iplast,ndims,xmat,kerre,
  7.      .                      precis,nitmax,nescri,ues,nnumer,deltax,
  8.      .                      iterlocal)
  9. CCC
  10. C      kerre   control de errores                      'KERRE
  11. C          = 0  integracion correcta
  12. C          = 1  sobrepasado el num max de iteraciones/no conv. local
  13. C          = 102  apex
  14. CCC
  15.       IMPLICIT INTEGER(I-N)
  16.       integer ndims,kerre,nitmax,nescri,nnumer,ues,kdummy,ndimv
  17.       integer iplast,kerr,i,iplcon,iplcap,iplapex,ndimx,nmodel
  18.       integer iplcon2,iplcap2,iplapex2,iterlocal
  19.       real*8    sigini(ndims),varini(4),ddefpl(ndims),
  20.      .          dsigtr(ndims),sigfin(ndims),varfin(ndims),
  21.      .          xmat(*),precis,deltax,void
  22.       real*8    sigtri(6),xtri(8),xini(8),lini,siginv(3)
  23.       real*8    exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
  24.      .          exl,xnn,aaa,bbb,pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,
  25.      .          exm2,eww2,cpr,cpm,cfr,you,xnu
  26.       common /elasdata/ you,xnu
  27.       common /MRScone/ exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
  28.      .                 exl,xnn,aaa,bbb
  29.       common /MRScapp/ pc0,exr,cca,phi,alp,pcc2,paa2,qaa2,exm2,eww2,
  30.      .                 cpr,cpm,cfr
  31.       do i=1,3
  32.           siginv(i)=0.D0
  33.         enddo
  34.         do i=1,6
  35.           sigtri(i)=0.D0
  36.         enddo
  37.         do i=1,8
  38.           xtri(i)=0.D0
  39.           xini(i)=0.D0
  40.         enddo
  41.       kerre=0
  42.       iplast=1
  43.       call CaracmateMAC(xmat)
  44.       nmodel=20
  45.       do i=1,ndims
  46.          sigtri(i)=sigini(i)+dsigtr(i)
  47.       enddo
  48.       if (nescri.eq.1) then
  49.          call InvariantesPQT(sigtri,ndims,siginv)
  50.          write(ues,*)'  Sigma trial invari'
  51.          write(ues,'(1x,3(E24.16,1x))')(siginv(i),i=1,3)
  52.       endif
  53.       kdummy=0
  54.       call DeterzonaMAC(sigtri,ndims,varini,iplcon,iplcap,iplapex,
  55.      .                  kdummy,ues)
  56. c ELASTICO
  57. 10    continue
  58.       if (iplapex.ne.0) goto 900
  59.       if ((iplcon.eq.0).and.(iplcap.eq.0)) then
  60.         if (nescri.eq.1) write(*,*)' Paso elastico'
  61.         iplast=0
  62.         do i=1,ndims
  63.           sigfin(i) = sigtri(i)
  64.         enddo
  65.         varfin(1)=varini(1)
  66.         varfin(2)=varini(2)
  67.         varfin(3)=-1.D0
  68.         varfin(4)=-1.D0
  69.         iterlocal=0
  70.         return
  71.       endif
  72. c PLASTICO
  73.          do i=1,ndims
  74.           xtri(i)=sigtri(i)
  75.           xini(i)=sigtri(i)
  76.          enddo
  77.          lini  =0.D0
  78.          if (nescri.eq.1) write(ues,*)' Modelo MAC: cone-cap harde'
  79.          ndimx=ndims+2
  80.          xtri(ndims+1)=varini(1)
  81.          xini(ndims+1)=varini(1)
  82.          xtri(ndims+2)=varini(2)
  83.          xini(ndims+2)=varini(2)
  84. c        prediccion para la var interna del cone:
  85. c        se trata de esquivar una zona de no convergencia del n-r local
  86. c         if ((xini(ndims+1).gt.1.D0).and.(xini(ndims+1).lt.5.D0))
  87. c     .     xini(ndims+1)=5.D0
  88.          if ((iplcon.eq.1).and.(iplcap.eq.0)) then
  89.            if (nescri.eq.1) write(*,*)' Paso plastico: CONE'
  90.            kdummy=21
  91.          elseif ((iplcon.eq.0).and.(iplcap.eq.1)) then
  92.            if (nescri.eq.1) write(*,*)' Paso plastico: CAP'
  93.            kdummy=22
  94.          else
  95.            write(*,*)' Paso plastico: CONE-CAP',iplcon,iplcap
  96.            kerre=1
  97.            return
  98.          endif
  99.          call Integra21(xtri,xini,ndimx,lini,ddefpl,ndims,
  100.      .                   kdummy,precis,nitmax,nescri,ues,nnumer,
  101.      .                   deltax,kerre,iterlocal)
  102.          varfin(1)=xini(ndims+1)
  103.          varfin(2)=xini(ndims+2)
  104.          do i=1,ndims
  105.            sigfin(i)=xini(i)
  106.          enddo
  107.          kdummy=0
  108.          call DeterzonaMAC(xtri,ndims,varfin,iplcon2,iplcap2,iplapex2,
  109.      .                    kdummy,ues)
  110.          if ((iplcon.eq.0).and.(iplcon2.eq.1)) then
  111.             if (nescri.eq.1)
  112.      .         write(*,*)' Cambio de zona de CAP a CONE'
  113.             iplcap=0
  114.             iplcon=1
  115.             goto 10
  116.          endif
  117.          if (kerre.eq.1) return
  118.          if ((iplcon.eq.1).and.(iplcap.eq.0)) then
  119.            varfin(3)=lini
  120.            varfin(4)=-1.D0
  121.          elseif ((iplcon.eq.0).and.(iplcap.eq.1)) then
  122.            varfin(3)=-1.D0
  123.            varfin(4)=lini
  124.          endif
  125.          return
  126. c APEX
  127.   900 continue
  128.       if (nescri.eq.1) write(*,*)' Paso plastico: APEX'
  129.       write(*,*)' Paso plastico: APEX'
  130.       do i=1,3
  131.         sigfin(i)=-xmat(5)
  132.       enddo
  133.       do i=4,ndims
  134.         sigfin(i)=0.d0
  135.       enddo
  136.       varfin(1)=varini(1)
  137.       varfin(2)=varini(2)
  138.       varfin(3)=-2.D0
  139.       varfin(4)=-2.D0
  140.       kerre=1
  141.       return
  142.  1000 continue
  143.       return
  144.       return
  145.       end
  146.  
  147.  
  148.  
  149.  

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