Télécharger fofiss.eso

Retour à la liste

Numérotation des lignes :

  1. C FOFISS SOURCE CB215821 19/08/01 21:15:59 10279
  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. CALL ACTOBJ('MMODEL',IPMODL,1)
  48. IF(IERR.NE.0)RETURN
  49. C
  50. C **** LECTURE DU PREMIER CHAMP/ELEMENT
  51. C
  52. CALL LIROBJ('MCHAML',IPIN,1,IRT1)
  53. CALL ACTOBJ('MCHAML',IPIN,1)
  54. IF(IERR.NE.0)RETURN
  55. CALL REDUAF(IPIN,IPMODL,IPCHE1,0,IR,KER)
  56. IF(IR .NE. 1) CALL ERREUR(KER)
  57. IF(IERR .NE. 0) RETURN
  58. C
  59. C **** LECTURE DU DEUXIEME CHAMP/ELEMENT
  60. C
  61. CALL LIROBJ('MCHAML',IPIN,1,IRT1)
  62. CALL ACTOBJ('MCHAML',IPIN,1)
  63. IF(IERR.NE.0) RETURN
  64. CALL REDUAF(IPIN,IPMODL,IPCHE2,0,IR,KER)
  65. IF(IR .NE. 1) CALL ERREUR(KER)
  66. IF(IERR .NE. 0) RETURN
  67. C
  68. C **** LECTURE EVENTUELLE DU TROISIEME CHAMP
  69. C
  70. CALL LIROBJ('MCHAML',IPIN,0,IRT1)
  71. IF(IERR.NE.0)RETURN
  72. IPCHA3=0
  73. IF (IRT1 .EQ. 1) THEN
  74. CALL ACTOBJ('MCHAML',IPIN,1)
  75. CALL REDUAF(IPIN,IPMODL,IPCHA3,0,IR,KER)
  76. IF(IR .NE. 1) CALL ERREUR(KER)
  77. IF(IERR .NE. 0) RETURN
  78. ENDIF
  79. C
  80. C **** LECTURE EVENTUELLE DU QUATRIEME CHAMP
  81. C
  82. CALL LIROBJ('MCHAML',IPIN,0,IRT1)
  83. IF(IERR.NE.0)RETURN
  84. IPCHA4=0
  85. IF (IRT1 .EQ. 1) THEN
  86. CALL ACTOBJ('MCHAML',IPIN,1)
  87. CALL REDUAF(IPIN,IPMODL,IPCHA4,0,IR,KER)
  88. IF(IR .NE. 1) CALL ERREUR(KER)
  89. IF(IERR .NE. 0) RETURN
  90. ENDIF
  91. C
  92. IF (IPCHA3.NE.0.OR.IPCHA4.NE.0) THEN
  93. CALL RNGCHA(IPCHA3,IPCHA4,'GRAFLEXI',
  94. 1 'CARACTERISTIQUES',IPCHE3,IPCHE4)
  95. IF(IERR.NE.0) RETURN
  96. ENDIF
  97. C
  98. CALL FOFIS1(IPMODL,IPCHE1,IPCHE2,IPCHE3,IPCHE4,IPCHP1,IRET)
  99. IF(IERR.NE.0) RETURN
  100. C
  101. IF(IRET.EQ.1) THEN
  102. CALL ACTOBJ('CHPOINT ',IPCHP1,1)
  103. CALL ECROBJ('CHPOINT ',IPCHP1)
  104. ENDIF
  105.  
  106. END
  107.  
  108.  
  109.  

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