Télécharger la plaquette Cast3M


fr en

Télécharger mtcmrs.eso

Retour à la liste

Numérotation des lignes : 

mtcmrs

  1. C MTCMRS    SOURCE    CHAT      05/01/13    01:53:51     5004
  2. CCC
  3. C **********************************************************************
  4. CCC
  5.       SUBROUTINE MTC_MRSMAC (KMAT,NDIMK,SIG,NDIMS,VINTCON,VINTCAP,LCON,
  6.      .                  lcap,xmat,nescri,ues,nnumer,deltax,kerre)
  7.       IMPLICIT INTEGER(I-N)
  8.       integer ndims,kerre,nescri,nnumer,ues,ndimk,
  9.      .          i,ifl,ndimx,kmod,ndimv
  10.       real*8    sig(ndims),vintcon,vintcap,Kmat(*),xmat(*),deltax,x(7),
  11.      .          lcon,lcap,lam
  12.       real*8    exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,exl,xnn,aaa,
  13.      .          bbb,pc0,exr,cca,phi,alp,voi1,voi2,voi3,voi4,voi5,cpr,
  14.      .          cpm,cfr
  15.       common /MRScone/ exm,qaa,pcc,eww,xk1,xk2,etb,exv,xep,paa,cco,
  16.      .                 exl,xnn,aaa,bbb
  17.       common /MRScapp/ pc0,exr,cca,phi,alp,voi1,voi2,voi3,voi4,voi5,
  18.      .                 cpr,cpm,cfr
  19.       do i=1,7
  20.           x(i)=0.D0
  21.         enddo
  22.       kerre=0
  23.       call CaracmateMAC(xmat)
  24.       if ((lcon.le.0.D0).and.(lcap.le.0.D0)) then
  25.         if (nescri.eq.1) write(ues,*)' Modelo MAC: paso elastico'
  26.         ifl=1
  27.         call MatHok(Kmat,ndimk,ndims,ifl)
  28.         goto 100
  29.       endif
  30.       if ((lcon.gt.0.D0).and.(lcap.gt.0.D0)) then
  31.          write(*,*)' Modelo MAC: no puede corner'
  32.          return
  33.       endif
  34.       do i=1,ndims
  35.          x(i)=sig(i)
  36.       enddo
  37.       if ((lcon.gt.0.D0).and.(lcap.lt.0.D0)) then
  38.         kmod=21
  39.         lam=lcon
  40.       else
  41.         kmod=22
  42.         lam=lcap
  43.       endif
  44.       if ((exr.lt.0.D0).and.(exv.lt.0.D0)) then
  45.         if (nescri.eq.1) write(ues,*)' Modelo MAC: plast perfecta'
  46.         write(*,*)' Formulacion no disponible'
  47.         kerre=1
  48. c         ndimv=0
  49. c         call MTC0(Kmat,ndimk,sig,ndims,ndimv,lam,
  50. c     .             kmod,nescri,ues,nnumer,deltax,kerre)
  51.       else if ((exr.lt.0.D0).and.(exv.gt.0.D0)) then
  52.         if (nescri.eq.1) write(ues,*)' Modelo MAC: cone hardening'
  53.         write(*,*)' Formulacion no disponible'
  54.         kerre=1
  55.       else if ((exr.gt.0.D0).and.(exv.lt.0.D0)) then
  56.         if (nescri.eq.1) write(ues,*)' Modelo MAC: cap hardening '
  57.         write(*,*)' Formulacion no disponible'
  58.         kerre=1
  59.       else if ((exr.gt.0.D0).and.(exv.gt.0.D0)) then
  60.         if (nescri.eq.1) write(ues,*)' Modelo MAC: cone-cap harde'
  61.         x(ndims+1)=vintcon
  62.         x(ndims+2)=vintcap
  63.         ndimx=ndims+2
  64.         call MTC21(Kmat,ndimk,x,ndimx,lam,ndims,
  65.      .             kmod,nescri,ues,nnumer,deltax,kerre)
  66.       endif
  67. 100   continue
  68.       return
  69.       end
  70.  
  71.  
  72.  
  73.  

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