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