Télécharger la plaquette Cast3M


fr en

Télécharger debit.eso

Retour à la liste

Numérotation des lignes : 

debit

  1. C DEBIT     SOURCE    OF166741  24/10/21    21:15:09     12042          
  2.  
  3. *_______________________________________________________________________
  4. *
  5. *                     OPERATEUR DEBIT  ( Inspire de PRESSI )
  6. *
  7. *       FQ = DEBIT MODEL I D MAIL I   (MOT1)
  8. *                        I CHPDEB I
  9. *
  10. *            FQ        CHPOINT CONTENANT LES FORCES NODALES EQUIVALENTES
  11. *            MODEL     OBJET MMODEL ,MASSIF , SUR LEQUEL S'APPLIQUE
  12. *                      LA CONDITION DE DEBIT
  13. *            D         VALEUR ALGEBRIQUE DU DEBIT
  14. *            MAIL      POUR LES MASSIFS ,OBJET MAILLAGE REPRESENTANT
  15. *                      LA FACE SUR LAQUELLE S'APPLIQUE LA CONDITION
  16. *                      DE DEBIT
  17. *            CHPOT     CHPOINT CONTENANT LES VALEURS ALGEBRIQUES
  18. *                      DES DEBITS AUX NOEUDS
  19. *            MOT1      NOM DONNE AU RESULAT ( PAR DEFAUT , CELUI DU MODELE)
  20. *______________________________________________________________________
  21.  
  22.       SUBROUTINE DEBIT
  23.  
  24.       IMPLICIT INTEGER(I-N)
  25.       IMPLICIT REAL*8(A-H,O-Z)
  26.  
  27. -INC PPARAM
  28. -INC CCOPTIO
  29.  
  30. -INC SMCHPOI
  31.  
  32.       CHARACTER*(4) MOSCAL,MOT1
  33.       CHARACTER*(8) TYPO
  34.  
  35.       REAL*8 XXX
  36.  
  37.       DATA MOSCAL/'SCAL'/
  38.  
  39.       MOT1 = '    '
  40.       TYPO = '        '
  41.  
  42.       IPTFP = 0
  43. *
  44. *     ON LIT UN MMODEL
  45. *
  46.       TYPO = 'MMODEL  '
  47.       CALL LIROBJ(TYPO,IPMODL,1,IRETOU)
  48.       IF (IERR.NE.0) RETURN
  49.       CALL ACTOBJ(TYPO,IPMODL,1)
  50.       IF (IERR.NE.0) RETURN
  51. *
  52. *     ON LIT SOIT UN FLOTTANT ET UN MAILLAGE ,SOIT UN CHPOINT
  53. *
  54.       CALL LIRREE(XXX,0,IRETOU)
  55.       IF (IERR.NE.0) RETURN
  56.       IF (IRETOU.NE.0) THEN
  57.          TYPO = 'MAILLAGE'
  58.          CALL LIROBJ(TYPO,IPTGEO,1,IRETOU)
  59.          IF (IERR.NE.0) RETURN
  60.          CALL ECRREE(XXX)
  61.          CALL ECRCHA(MOSCAL)
  62.          CALL ECRENT(1)
  63.          CALL ECROBJ(TYPO,IPTGEO)
  64.          CALL MANUCH
  65.          IF (IERR.NE.0) RETURN
  66.          TYPO = 'CHPOINT '
  67.          CALL LIROBJ(TYPO,ICHPOP,1,IRETOU)
  68.          IF (IERR.NE.0) RETURN
  69.       ELSE
  70.          TYPO = 'CHPOINT '
  71.          CALL LIROBJ(TYPO,ICHPOP,1,IRETOU)
  72.          IF (IERR.NE.0) RETURN
  73.       ENDIF
  74.       CALL ACTOBJ(TYPO,ICHPOP,1)
  75.       IF (IERR.NE.0) RETURN
  76. *
  77.       CALL LIRCHA(MOT1,0,IRETOU)
  78.       IF (IERR.NE.0) RETURN
  79. *
  80.       CALL FDEBIT(ICHPOP,IPMODL,IPTFP,MOT1,IRET)
  81.       IF (IRET.EQ.0) RETURN
  82. *
  83. *        LE NUMERO DE L HARMONIQUE EST PRIS EGAL A NIFOUR
  84. *        POUR TOUTES LES COMPOSANTES DU CHPOINT
  85. *
  86.       MCHPOI=IPTFP
  87.       SEGACT MCHPOI*MOD
  88. *         ATTRIBUTION D'UNE NATURE DISCRETE AU CHPO SORTANT
  89.       NAT = MAX ( 1 , JATTRI(/1) )
  90.       NSOUPO = IPCHP(/1)
  91.       SEGADJ MCHPOI
  92.       JATTRI(1) = 2
  93.  
  94.       DO IA=1,IPCHP(/1)
  95.          MSOUPO=IPCHP(IA)
  96.          SEGACT MSOUPO*mod
  97.          DO NC=1,NOHARM(/1)
  98.             NOHARM(NC)=NIFOUR
  99.          ENDDO
  100.       ENDDO
  101.  
  102.       TYPO = 'CHPOINT '
  103.       CALL ACTOBJ(TYPO,IPTFP,1)
  104.       CALL ECROBJ(TYPO,IPTFP)
  105.  
  106. C      RETURN
  107.       END
  108.  
  109.  
  110.  

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