Télécharger fofiss.eso

Retour à la liste

Numérotation des lignes :

  1. C FOFISS SOURCE CB215821 16/12/05 21:39:29 9237
  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. -INC CCOPTIO
  35. C
  36. IPMODL=0
  37. IPCHE1=0
  38. IPCHE2=0
  39. IPCHE3=0
  40. IPCHE4=0
  41. IPCHA3=0
  42. IPCHA4=0
  43. C
  44. C **** LECTURE D'UN OBJET MMODEL
  45. C
  46. CALL LIROBJ('MMODEL',IPMODL,1,IRT1)
  47. IF(IERR.NE.0)RETURN
  48. C
  49. C **** LECTURE DU PREMIER CHAMP/ELEMENT
  50. C
  51. CALL LIROBJ('MCHAML',IPIN,1,IRT1)
  52. IF(IERR.NE.0)RETURN
  53. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  54. IF(IR .NE. 1) CALL ERREUR(KER)
  55. IF(IERR .NE. 0) RETURN
  56. C
  57. C **** LECTURE DU DEUXIEME CHAMP/ELEMENT
  58. C
  59. CALL LIROBJ('MCHAML',IPIN,1,IRT1)
  60. IF(IERR.NE.0) RETURN
  61. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  62. IF(IR .NE. 1) CALL ERREUR(KER)
  63. IF(IERR .NE. 0) RETURN
  64. C
  65. C **** LECTURE EVENTUELLE DU TROISIEME CHAMP
  66. C
  67. CALL LIROBJ('MCHAML',IPIN,0,IRT1)
  68. IF(IERR.NE.0)RETURN
  69. IPCHA3=0
  70. IF (IRT1 .EQ. 1) THEN
  71. CALL REDUAF(IPIN,IPMODL,IPCHA3,0,IR,KER)
  72. IF(IR .NE. 1) CALL ERREUR(KER)
  73. IF(IERR .NE. 0) RETURN
  74. ENDIF
  75. C
  76. C **** LECTURE EVENTUELLE DU QUATRIEME CHAMP
  77. C
  78. CALL LIROBJ('MCHAML',IPIN,0,IRT1)
  79. IF(IERR.NE.0)RETURN
  80. IPCHA4=0
  81. IF (IRT1 .EQ. 1) THEN
  82. CALL REDUAF(IPIN,IPMODL,IPCHA4,0,IR,KER)
  83. IF(IR .NE. 1) CALL ERREUR(KER)
  84. IF(IERR .NE. 0) RETURN
  85. ENDIF
  86. C
  87. IF (IPCHA3.NE.0.OR.IPCHA4.NE.0) THEN
  88. CALL RNGCHA(IPCHA3,IPCHA4,'GRAFLEXI',
  89. 1 'CARACTERISTIQUES',IPCHE3,IPCHE4)
  90. IF(IERR.NE.0) RETURN
  91. ENDIF
  92. C
  93. CALL FOFIS1(IPMODL,IPCHE1,IPCHE2,IPCHE3,IPCHE4,IPCHP1,IRET)
  94. IF(IERR.NE.0) RETURN
  95. C
  96. IF(IRET.EQ.1) CALL ECROBJ('CHPOINT',IPCHP1)
  97. C
  98. RETURN
  99. END
  100.  
  101.  
  102.  

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