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  23/06/21    21:15:03     11690          
  2.       SUBROUTINE DEBIT
  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.       IMPLICIT INTEGER(I-N)
  23.       IMPLICIT REAL*8(A-H,O-Z)
  24.  
  25. -INC PPARAM
  26. -INC CCOPTIO
  27.  
  28. -INC SMCHPOI
  29.  
  30.       CHARACTER*(4) MOSCAL,MOT1
  31.       CHARACTER*(8) TYPO
  32.  
  33.       REAL*8 XXX
  34.  
  35.       DATA MOSCAL/'SCAL'/
  36.  
  37.       MOT1 = '    '
  38.       TYPO = '        '
  39.  
  40.       IPTFP = 0
  41. *
  42. *     ON LIT UN MMODEL
  43. *
  44.       TYPO = 'MMODEL  '
  45.       CALL LIROBJ(TYPO,IPMODL,1,IRETOU)
  46.       IF (IERR.NE.0) RETURN
  47. *
  48. *     ON LIT SOIT UN FLOTTANT ET UN MAILLAGE ,SOIT UN CHPOINT
  49. *
  50.       CALL LIRREE(XXX,0,IRETOU)
  51.       IF (IERR.NE.0) RETURN
  52.       IF (IRETOU.NE.0) THEN
  53.          TYPO = 'MAILLAGE'
  54.          CALL LIROBJ(TYPO,IPTGEO,1,IRETOU)
  55.          IF (IERR.NE.0) RETURN
  56.          CALL ECRREE(XXX)
  57.          CALL ECRCHA(MOSCAL)
  58.          CALL ECRENT(1)
  59.          CALL ECROBJ(TYPO,IPTGEO)
  60.          CALL MANUCH
  61.          IF (IERR.NE.0) RETURN
  62.          TYPO = 'CHPOINT '
  63.          CALL LIROBJ(TYPO,ICHPOP,1,IRETOU)
  64.          IF (IERR.NE.0) RETURN
  65.       ELSE
  66.          TYPO = 'CHPOINT '
  67.          CALL LIROBJ(TYPO,ICHPOP,1,IRETOU)
  68.          IF (IERR.NE.0) RETURN
  69.       ENDIF
  70. *
  71.       CALL LIRCHA(MOT1,0,IRETOU)
  72.       IF (IERR.NE.0) RETURN
  73. *
  74.       CALL FDEBIT(ICHPOP,IPMODL,IPTFP,MOT1,IRET)
  75.       IF (IRET.EQ.0) RETURN
  76. *
  77. *        LE NUMERO DE L HARMONIQUE EST PRIS EGAL A NIFOUR
  78. *        POUR TOUTES LES COMPOSANTES DU CHPOINT
  79. *
  80.       MCHPOI=IPTFP
  81.       SEGACT MCHPOI*MOD
  82. *         ATTRIBUTION D'UNE NATURE DISCRETE AU CHPO SORTANT
  83.       NAT = MAX ( 1 , JATTRI(/1) )
  84.       NSOUPO = IPCHP(/1)
  85.       SEGADJ MCHPOI
  86.       JATTRI(1) = 2
  87.  
  88.       DO IA=1,IPCHP(/1)
  89.          MSOUPO=IPCHP(IA)
  90.          SEGACT MSOUPO*mod
  91.          DO NC=1,NOHARM(/1)
  92.             NOHARM(NC)=NIFOUR
  93.          ENDDO
  94.       ENDDO
  95.  
  96.       TYPO = 'CHPOINT '
  97.       CALL ACTOBJ(TYPO,IPTFP,1)
  98.       CALL ECROBJ(TYPO,IPTFP)
  99.  
  100. C      RETURN
  101.       END
  102.  
  103.  
  104.  

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