Télécharger defpou.eso

Retour à la liste

Numérotation des lignes :

  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)
  45. CALL DOXE(XCOOR,IDIM,NBNN,NUM,IB,XE)
  46. XL=SQRT ( (XE(1,2)-XE(1,1))**2 +
  47. & (XE(2,2)-XE(2,1))**2 +
  48. & (XE(3,2)-XE(3,1))**2 )
  49. MPTVAL=IVASTF
  50. MELVAL=IVAL(6)
  51. SIG16=VELCHE(1,IB)
  52. SIG26=VELCHE(2,IB)
  53. MELVAL=IVAL(2)
  54. VELCHE(1,IB)=(SIG16 - SIG26)/XL
  55. VELCHE(2,IB)=VELCHE(1,IB)
  56. MELVAL=IVAL(5)
  57. SIG15=VELCHE(1,IB)
  58. SIG25=VELCHE(2,IB)
  59. MELVAL=IVAL(3)
  60. VELCHE(1,IB)=(SIG25 - SIG15)/XL
  61. VELCHE(2,IB)=VELCHE(1,IB)
  62. ENDIF
  63. c
  64. RETURN
  65. END
  66.  
  67.  
  68.  
  69.  
  70.  

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