Télécharger la plaquette Cast3M


fr en

Télécharger resdpr.eso

Retour à la liste

Numérotation des lignes : 

resdpr

  1. C RESDPR    SOURCE    PV        22/04/19    16:18:11     11344          
  2.  
  3.       SUBROUTINE RESIDUO_DPRAL(X,R,N,NDIMX,NMODEL,XTRI,SIG,VECM)
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6.       dimension void(1)
  7.       integer n,i,j,ndims,ndimv,nmodel,ndimx
  8.       real*8  x(n),r(n),void,sig(4),vecm(4),xtri(ndimx)
  9.       integer             augla
  10.       real*8                    c
  11.       common /auglagrang1/ augla
  12.       common /auglagrang2/ c
  13.       call zzero(sig,4)
  14.       call zzero(vecm,4)
  15.       void=0.D0
  16.       ndimv=1
  17.       ndims=3
  18.       call der_enerelas_dpral(x,sig,nmodel)
  19.       if (ndimx.eq.3) then
  20.        call yieldd(sig,ndims,void,ndimv,r(n),nmodel)
  21.        call vflsig(sig,ndims,void,ndimv,vecm,nmodel)
  22.       else if (ndimx.eq.4) then
  23.        call yieldd(sig,ndims,x(ndimx),ndimv,r(n),nmodel)
  24.        call vflsig(sig,ndims,x(ndimx),ndimv,vecm,nmodel)
  25.        call vflvar(sig,ndims,x(ndimx),ndimv,vecm(ndimx),nmodel)
  26.       endif
  27.       do i=1,ndimx
  28.          r(i)=x(i)-xtri(i)+ABS(x(n)+c*r(n))*vecm(i)
  29.       enddo
  30. C
  31. c VISCOSIDAD FICTICIA AÑADIDA
  32. c      r(n)=r(n)-1.D-8
  33. c
  34.       return
  35.       end
  36.  
  37.  
  38.  
  39.  
  40.  
  41.  
  42.  

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