Esope source files

   1 : C DEFPOU    SOURCE    PV        09/03/12    21:19:10     6325           
   2 :       SUBROUTINE DEFPOU(MFR,MELE,MELEME,IB,WRK4,IVASTF)   
   3 : * 
   4 : ****************************************************
   5 : *  ENTREES  
   6 : ****************************************************  
   7 : *  
   8 : *  MFR   : formulation de l'élément  
   9 : *  MELE  : numéro élément fini  
  10 : *  MELEME: pointeur sur le maillage 
  11 : *  NBNN  : nombre de noeuds par élément (WRK4)  
  12 : *  IB    : numéro de l'élément  
  13 : *  IVASTF: pointeur sur un segment mptval de contraintes 
  14 : * 
  15 : **************************************************** 
  16 : *  SORTIES  
  17 : ****************************************************  
  18 : *  
  19 : *  XE    : coordonnées des points de l'élément (WRK4)  
  20 : * 
  21 : ****************************************************   
  22 :       IMPLICIT INTEGER(I-N)
  23 :       IMPLICIT REAL*8(A-H,O-Z)
  24 : *
  25 : -INC CCOPTIO
  26 : -INC SMCHAML
  27 : -INC SMELEME
  28 : -INC SMCOORD
  29 : *  
  30 :       SEGMENT MPTVAL
  31 :         INTEGER IPOS(NS)    ,NSOF(NS)
  32 :         INTEGER      IVAL(NCOSOU)
  33 :         CHARACTER*16 TYVAL(NCOSOU)
  34 :       ENDSEGMENT
  35 : *  
  36 :       SEGMENT WRK4
  37 :         REAL*8 XE(3,NBNN)
  38 :       ENDSEGMENT
  39 : *  
  40 : c
  41 : c    special poutres et tuyaux  sauf timoschenko
  42 : c
  43 :          IF((MFR.EQ.7.OR.MFR.EQ.13).AND.MELE.NE.84) THEN
  44 :             NBNN=XE(/2) 
  44 :             CALL DOXE(XCOOR,IDIM,NBNN,NUM,IB,XE)
  45 :             XL=SQRT ( (XE(1,2)-XE(1,1))**2 +
  46 :      &                (XE(2,2)-XE(2,1))**2 +
  47 :      &                (XE(3,2)-XE(3,1))**2   )
  48 :             MPTVAL=IVASTF
  49 :             MELVAL=IVAL(6)
  50 :             SIG16=VELCHE(1,IB)
  51 :             SIG26=VELCHE(2,IB)
  52 :             MELVAL=IVAL(2)
  53 :             VELCHE(1,IB)=(SIG16 - SIG26)/XL
  54 :             VELCHE(2,IB)=VELCHE(1,IB)
  55 :             MELVAL=IVAL(5)
  56 :             SIG15=VELCHE(1,IB)
  57 :             SIG25=VELCHE(2,IB)
  58 :             MELVAL=IVAL(3)
  59 :             VELCHE(1,IB)=(SIG25 - SIG15)/XL
  60 :             VELCHE(2,IB)=VELCHE(1,IB)
  61 :          ENDIF
  62 : c 
  63 :          RETURN 
  64 :          END   
  65 :  
  66 :  
  67 :  
  68 :