Télécharger la plaquette Cast3M


fr en

Télécharger compou.eso

Retour à la liste

Numérotation des lignes : 

compou

  1. C COMPOU    SOURCE    PV        17/12/08    21:16:50     9660           
  2. *
  3.       SUBROUTINE COMPOU(IB,MWRKXE,ililuc,iwrk53)
  4. *
  5. ****************************************************
  6. *  ENTREES
  7. ****************************************************
  8. *
  9. *  MFR   : formulation de l'élément
  10. *  MELE  : numéro élément fini
  11. *  NBNN  : nombre de noeuds par élément (WRK4)
  12. *  IB    : numéro de l'élément
  13. *  XEL   : coordonnées des points de l'élément (MWRK2)
  14. *
  15. ****************************************************
  16. *  SORTIES
  17. ****************************************************
  18. *
  19. ****************************************************
  20.       IMPLICIT INTEGER(I-N)
  21.       IMPLICIT REAL*8(A-H,O-Z)
  22. *
  23.  
  24. -INC PPARAM
  25. -INC CCOPTIO
  26. *
  27. -INC SMCHAML
  28. * segment deroulant le mcheml
  29. -INC DECHE
  30. *
  31.       SEGMENT MWRKXE
  32.         REAL*8 XEL(3,NBNN)
  33.       ENDSEGMENT
  34. c
  35.       liluc = ililuc
  36.       wrk53 = iwrk53
  37. C
  38. C  XEL est rempli systematiquement dans le sous-programme appelant
  39. c
  40. c    special poutres et tuyaux  sauf timoschenko
  41. c
  42.       IF ((MFR.EQ.7.OR.MFR.EQ.13).AND.MELE.NE.84) THEN
  43. C
  44.         IF (IDIM.EQ.3) THEN
  45. C
  46.             XL=SQRT ( (XEL(1,2)-XEL(1,1))**2 +
  47.      &                (XEL(2,2)-XEL(2,1))**2 +
  48.      &                (XEL(3,2)-XEL(3,1))**2   )
  49.             pilnec = liluc(11,2)
  50. *******            segact pilnec*nomod
  51.             mran = pilobl(/2)
  52.             deche = pilobl(6,mran)
  53. *********            segact deche*nomod
  54.             melval = ieldec
  55.             DSIGP = VELCHE(1,IB) - VELCHE(2,IB)
  56.             deche = pilobl(2,mran)
  57. ********            segact deche*nomod
  58.             melval = ieldec
  59. *********            segact melval*mod
  60.             VELCHE(1,IB)= DSIGP / XL
  61.             VELCHE(2,IB)=VELCHE(1,IB)
  62.             deche = pilobl(5,mran)
  63. ***********            segact deche*nomod
  64.             melval = ieldec
  65.             DSIGP = VELCHE(2,IB) - VELCHE(1,IB)
  66.             deche = pilobl(3,mran)
  67. ***********            segact deche*nomod
  68.             melval = ieldec
  69. ************            segact melval*mod
  70.             VELCHE(1,IB)= DSIGP / XL
  71.             VELCHE(2,IB)=VELCHE(1,IB)
  72. C
  73.         ELSE
  74. C  CAS 2D
  75.             XL=SQRT ( (XEL(1,2)-XEL(1,1))**2 +
  76.      &                (XEL(2,2)-XEL(2,1))**2 )
  77.             pilnec = liluc(11,2)
  78. ***********            segact pilnec*nomod
  79.             mran = pilobl(/2)
  80.             deche = pilobl(3,mran)
  81. *************            segact deche*nomod
  82.             melval = ieldec
  83.             DSIGP = VELCHE(1,IB) - VELCHE(2,IB)
  84. *
  85.             deche = pilobl(2,mran)
  86. ************            segact deche*nomod
  87.             melval = ieldec
  88. ***************            segact melval*mod
  89.             VELCHE(1,IB)= DSIGP / XL
  90.             VELCHE(2,IB)=VELCHE(1,IB)
  91. C
  92.         ENDIF
  93. C
  94.       ENDIF
  95. c
  96.       RETURN
  97.       END
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  

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