Sources Esope | Cast3M

defpou
C DEFPOU    SOURCE    OF166741  25/02/21    21:15:50     12166          
      SUBROUTINE DEFPOU(MFR,MELE,MELEME,IB,WRK4,IVASTF)
*
****************************************************
*  ENTREES
****************************************************
*
*  MFR   : formulation de l'élément
*  MELE  : numéro élément fini
*  MELEME: pointeur sur le maillage
*  NBNN  : nombre de noeuds par élément (WRK4)
*  IB    : numéro de l'élément
*  IVASTF: pointeur sur un segment mptval de contraintes
*
****************************************************
*  SORTIES
****************************************************
*
*  XE    : coordonnées des points de l'élément (WRK4)
*
****************************************************
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
 
-INC SMCHAML
-INC SMELEME
-INC SMCOORD
 
-INC TMPTVAL
 
      SEGMENT WRK4
        REAL*8 XE(3,NBNN)
      ENDSEGMENT
c
c    special poutres et tuyaux  sauf timoschenko
c
      IF((MFR.EQ.7.OR.MFR.EQ.13).AND.MELE.NE.84) THEN
         NBNN=XE(/2)
         CALL DOXE(XCOOR,IDIM,NBNN,NUM,IB,XE)
         XL=SQRT ( (XE(1,2)-XE(1,1))**2 +
     &             (XE(2,2)-XE(2,1))**2 +
     &             (XE(3,2)-XE(3,1))**2   )
         MPTVAL=IVASTF
         MELVAL=IVAL(6)
         SIG16=VELCHE(1,IB)
         SIG26=VELCHE(2,IB)
         MELVAL=IVAL(2)
         VELCHE(1,IB)=(SIG16 - SIG26)/XL
         VELCHE(2,IB)=VELCHE(1,IB)
         MELVAL=IVAL(5)
         SIG15=VELCHE(1,IB)
         SIG25=VELCHE(2,IB)
         MELVAL=IVAL(3)
         VELCHE(1,IB)=(SIG25 - SIG15)/XL
         VELCHE(2,IB)=VELCHE(1,IB)
      ENDIF
 
      RETURN
      END