Télécharger fofiss.eso

Retour à la liste

Numérotation des lignes :

fofiss
  1. C FOFISS SOURCE CB215821 23/01/25 21:15:14 11573
  2. C
  3. SUBROUTINE FOFISS
  4. C
  5. C=======================================================================
  6. C
  7. C OPERATEUR FORCES INTERNES
  8. C
  9. C CET OPERATEUR NOUS PERMET DE CALCULER LE CHPOIN FOR1
  10. C DES FORCES NODALES RESULTANT DE L'INTEGRATION DU CHA-
  11. C MELEM SIG1 DE CONTRAINTES ISSU PAR EXEMPLE DE L'OPE-
  12. C RATEUR SIGMA ET DU CHAMELEM GRAD1 DE GRADIENTS ( ET
  13. C GRAF1 DE GRAFLEXIS ) ISSU PAR EXEMPLE DE L'OPERATEUR
  14. C GRAD ( ET GRAF ).
  15. C
  16. C SYNTAXE : FOR1=FOFISS MOD1 SIG1 GRAD1 (GRAF1)(CAR1) ;
  17. C
  18. C ENTREES :
  19. C MOD1=OBJET DE TYPE MMODEL
  20. C SIG1=CHAPM /ELMTS DE CONTRAINTES (TYPE MCHAML)
  21. C GRAD1=CHAPM /ELMTS DE GRADIENTS (TYPE MCHAML)
  22. C GRAF1=CHAPM /ELMTS DE GRAFLEXIS (TYPE MCHAML)
  23. C CAR1=CHAMP /ELMTS DE CARACTERISTIQUES (FACULTATIF)
  24. C POUR LES POUTRES ET LES COQUES (TYPE MCHAML)
  25. C SORTIES:
  26. C FOR1 =CHAMP PAR POINT DONNANT LES FORCES NODALES
  27. C CODE JUILLET 87 SUO X.Z.
  28. C PASSAGE AUX NOUVEAUX MCHAMLS PAR P.DOWLATYARI SEP 91
  29. C=======================================================================
  30. C
  31. IMPLICIT INTEGER(I-N)
  32. IMPLICIT REAL*8(A-H,O-Z)
  33. C
  34.  
  35. -INC SMCOORD
  36. -INC PPARAM
  37. -INC CCOPTIO
  38. C
  39. IPMODL=0
  40. IPCHE1=0
  41. IPCHE2=0
  42. IPCHE3=0
  43. IPCHE4=0
  44. IPCHA3=0
  45. IPCHA4=0
  46. C
  47. C **** LECTURE D'UN OBJET MMODEL
  48. C
  49. CALL LIROBJ('MMODEL',IPMODL,1,IRT1)
  50. CALL ACTOBJ('MMODEL',IPMODL,1)
  51. IF(IERR.NE.0)RETURN
  52. C
  53. C **** LECTURE DU PREMIER CHAMP/ELEMENT
  54. C
  55. CALL LIROBJ('MCHAML',IPIN,1,IRT1)
  56. CALL ACTOBJ('MCHAML',IPIN,1)
  57. IF(IERR.NE.0)RETURN
  58. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  59. IF(IR .NE. 1) CALL ERREUR(KER)
  60. IF(IERR .NE. 0) RETURN
  61. C
  62. C **** LECTURE DU DEUXIEME CHAMP/ELEMENT
  63. C
  64. CALL LIROBJ('MCHAML',IPIN,1,IRT1)
  65. CALL ACTOBJ('MCHAML',IPIN,1)
  66. IF(IERR.NE.0) RETURN
  67. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  68. IF(IR .NE. 1) CALL ERREUR(KER)
  69. IF(IERR .NE. 0) RETURN
  70. C
  71. C **** LECTURE EVENTUELLE DU TROISIEME CHAMP
  72. C
  73. CALL LIROBJ('MCHAML',IPIN,0,IRT1)
  74. IF(IERR.NE.0)RETURN
  75. IPCHA3=0
  76. IF (IRT1 .EQ. 1) THEN
  77. CALL ACTOBJ('MCHAML',IPIN,1)
  78. CALL REDUAF(IPIN,IPMODL,IPCHA3,0,IR,KER)
  79. IF(IR .NE. 1) CALL ERREUR(KER)
  80. IF(IERR .NE. 0) RETURN
  81. ENDIF
  82. C
  83. C **** LECTURE EVENTUELLE DU QUATRIEME CHAMP
  84. C
  85. CALL LIROBJ('MCHAML',IPIN,0,IRT1)
  86. IF(IERR.NE.0)RETURN
  87. IPCHA4=0
  88. IF (IRT1 .EQ. 1) THEN
  89. CALL ACTOBJ('MCHAML',IPIN,1)
  90. CALL REDUAF(IPIN,IPMODL,IPCHA4,0,IR,KER)
  91. IF(IR .NE. 1) CALL ERREUR(KER)
  92. IF(IERR .NE. 0) RETURN
  93. ENDIF
  94. C
  95. IF (IPCHA3.NE.0.OR.IPCHA4.NE.0) THEN
  96. CALL RNGCHA(IPCHA3,IPCHA4,'GRAFLEXI',
  97. 1 'CARACTERISTIQUES',IPCHE3,IPCHE4)
  98. IF(IERR.NE.0) RETURN
  99. ENDIF
  100. C
  101. SEGACT,MCOORD
  102. CALL FOFIS1(IPMODL,IPCHE1,IPCHE2,IPCHE3,IPCHE4,IPCHP1,IRET)
  103. SEGDES,MCOORD
  104. IF(IERR.NE.0) RETURN
  105. C
  106. IF(IRET.EQ.1) THEN
  107. CALL ACTOBJ('CHPOINT ',IPCHP1,1)
  108. CALL ECROBJ('CHPOINT ',IPCHP1)
  109. ENDIF
  110.  
  111. END
  112.  
  113.  
  114.  
  115.  

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