Télécharger la plaquette Cast3M


fr en

Télécharger defpou.eso

Retour à la liste

Numérotation des lignes : 

defpou

  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.  
  26. -INC PPARAM
  27. -INC CCOPTIO
  28. -INC SMCHAML
  29. -INC SMELEME
  30. -INC SMCOORD
  31. *
  32.       SEGMENT MPTVAL
  33.         INTEGER IPOS(NS)    ,NSOF(NS)
  34.         INTEGER      IVAL(NCOSOU)
  35.         CHARACTER*16 TYVAL(NCOSOU)
  36.       ENDSEGMENT
  37. *
  38.       SEGMENT WRK4
  39.         REAL*8 XE(3,NBNN)
  40.       ENDSEGMENT
  41. *
  42. c
  43. c    special poutres et tuyaux  sauf timoschenko
  44. c
  45.          IF((MFR.EQ.7.OR.MFR.EQ.13).AND.MELE.NE.84) THEN
  46.             NBNN=XE(/2)
  47.             CALL DOXE(XCOOR,IDIM,NBNN,NUM,IB,XE)
  48.             XL=SQRT ( (XE(1,2)-XE(1,1))**2 +
  49.      &                (XE(2,2)-XE(2,1))**2 +
  50.      &                (XE(3,2)-XE(3,1))**2   )
  51.             MPTVAL=IVASTF
  52.             MELVAL=IVAL(6)
  53.             SIG16=VELCHE(1,IB)
  54.             SIG26=VELCHE(2,IB)
  55.             MELVAL=IVAL(2)
  56.             VELCHE(1,IB)=(SIG16 - SIG26)/XL
  57.             VELCHE(2,IB)=VELCHE(1,IB)
  58.             MELVAL=IVAL(5)
  59.             SIG15=VELCHE(1,IB)
  60.             SIG25=VELCHE(2,IB)
  61.             MELVAL=IVAL(3)
  62.             VELCHE(1,IB)=(SIG25 - SIG15)/XL
  63.             VELCHE(2,IB)=VELCHE(1,IB)
  64.          ENDIF
  65. c
  66.          RETURN
  67.          END
  68.  
  69.  
  70.  
  71.  
  72.  

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