Télécharger comrot.eso

Retour à la liste

Numérotation des lignes :

  1. C COMROT SOURCE PV 17/12/08 21:16:51 9660
  2.  
  3. SUBROUTINE COMROT(WRK53,IB,IPTR1,MWRKXE,WRK54)
  4. *
  5. **********************************************************
  6. * ENTREES
  7. **********************************************************
  8. *
  9. * CMATE : nom du matériau élastique
  10. * MFR : formulation de l'élément
  11. * NBBN : nombre de noeuds de l'élément
  12. * IB : numéro de l'élément
  13. * MELE : numéro élément fini
  14. * IPTR1 : pointeur sur un segment MINTE2
  15. *
  16. **********************************************************
  17. * SORTIES
  18. **********************************************************
  19. * TXR : cosinus directeurs des axes locaux de l'élément massif
  20. * (WTRAV) a remplir si IPTR1 n'est pas NUL !
  21. *
  22. **********************************************************
  23. *
  24. IMPLICIT INTEGER(I-N)
  25. IMPLICIT REAL*8(A-H,O-Z)
  26. *
  27. -INC CCOPTIO
  28. *
  29. -INC SMINTE
  30. -INC DECHE
  31. *
  32. SEGMENT MWRKXE
  33. REAL*8 XEL(3,NBNN)
  34. ENDSEGMENT
  35. c
  36. c calcul des axes locaux
  37. c
  38. IF (IPTR1.GT.0) THEN
  39. C* IF((CMATE.EQ.'ORTHOTRO'.OR.CMATE.EQ.'ANISOTRO'.OR.
  40. C* & CMATE.EQ.'UNIDIREC').AND.(MFR.EQ.1
  41. C* & .OR.MFR.EQ.31.OR.MFR.EQ.33)) THEN
  42. MINTE2=IPTR1
  43. NBSH=MINTE2.SHPTOT(/2)
  44. NBNN = XEL(/2)
  45. CALL RLOCAL(XEL,MINTE2.SHPTOT,NBSH,NBNN,TXR)
  46. if (nbsh.eq.-1) then
  47. call erreur(525)
  48. return
  49. endif
  50. ENDIF
  51. *
  52. CPPU IF(LUNI1)THEN
  53. CPPU NBNN = XEL(/2)
  54. CPPU NLG=NUMGEO(MELE)
  55. CPPU CALL RESHPT(1,NBNN,NLG,MELE,0,IPTR1,IRT1)
  56. CPPU MINTE2=IPTR1
  57. CPPU SEGACT MINTE2
  58. CPPU NBSH=MINTE2.SHPTOT(/2)
  59. CPPU CALL RLOCAL (XEL,MINTE2.SHPTOT,NBSH,NBNN,TXR)
  60. CPPU ENDIF
  61. *
  62. RETURN
  63. END
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70.  

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