Télécharger la plaquette Cast3M


fr en

Télécharger dimen3.eso

Retour à la liste

Numérotation des lignes : 

dimen3

  1. C DIMEN3    SOURCE    PV090527  24/10/30    21:15:03     12061          
  2.       SUBROUTINE DIMEN3 (IPRIGI,NBRINC)
  3. ************************************************************************
  4. *
  5. *                             D I M E N 3
  6. *                             -----------
  7. *
  8. * FONCTION:
  9. * ---------
  10. *
  11. *     DONNER LE NOMBRE D'INCONNUES DU PROBLEME PHYSIQUE ASSOCIE A UN
  12. *     OBJET 'RIGIDITE'.
  13. *
  14. * MODE D'APPEL:
  15. * -------------
  16. *
  17. *     CALL DIMEN3 (IPRIGI,NBRINC)
  18. *
  19. * PARAMETRES:  (E)=ENTREE   (S)=SORTIE
  20. * -----------
  21. *
  22. *     IPRIGI  ENTIER    (E)  POINTEUR D'UNE 'RIGIDITE'.
  23. *     NBRINC  ENTIER    (S)  NOMBRE D'INCONNUES DE LA 'RIGIDITE'.
  24. *
  25. * REMARQUE :
  26. * ----------
  27. *
  28. *     IL N'EST PAS NECESSAIRE DE FACTORISER L'OBJET RIGIDITE POUR 
  29. *     CONNAITRE LE NOMBRE D'INCONNUES DU PROBLEME
  30. *     APPEL A ASSEM1/ASNS1 PLUTOT QUE TRIANG/LDMT1
  31. *
  32. * AUTEUR, DATE DE CREATION:
  33. * -------------------------
  34. *
  35. *     PASCAL MANIGOT     26 AOUT 1985
  36. *
  37. * LANGAGE:
  38. * --------
  39. *
  40. *     ESOPE + FORTRAN77
  41. *
  42. ************************************************************************
  43. *
  44.       IMPLICIT INTEGER(I-N)
  45. -INC PPARAM
  46. -INC CCOPTIO
  47. -INC SMMATRI
  48. -INC SMRIGID
  49. -INC CCREEL
  50. C
  51.       SEGMENT,IPOS(NNOE1)
  52. C
  53.       INTEGER INSYM
  54.       INSYM = 0
  55.       MRIGID = IPRIGI
  56.       SEGACT,MRIGID
  57.       NRG = IRIGEL(/1)
  58.       NBR = IRIGEL(/2)
  59. C     bp, 2014-11-27 : petit ajout pour les rigidites vides
  60.       IF(NBR.EQ.0) THEN
  61.         NBRINC=0
  62.         SEGDES,MRIGID
  63.         RETURN
  64.       ENDIF
  65.       IPMATR = ICHOLE
  66. C
  67.       IF(NORINC.GT.0  .AND. NORIND.GT.0) THEN
  68.          INSYM = 1
  69.       ENDIF
  70.       IF (NRG.GE.7) THEN
  71.          DO  9 IN = 1,NBR
  72.             IANTI=IRIGEL(7,IN)
  73.             IF(IANTI.GT.0) THEN
  74.                INSYM = 1
  75.             ENDIF
  76.   9      CONTINUE
  77.       ENDIF
  78. C
  79.       IF (IPMATR.EQ.0) THEN
  80.         CALL ECRCHA('NOOP')
  81.         IF (INSYM.EQ.0) THEN
  82.           CALL ASSEM1(IPRIGI,MMATRX,INUINX,ITOPOX,IMINIX,IPOX,INCTRX,
  83.      &                IITOPX)
  84.         ELSE
  85.           CALL ASNS1 (IPRIGI,MMATRX,INUINX,ITOPOX,IMINIX,IPOX,INCTRX,
  86.      &                INCTRZ,IITOPX,ITOPOD,IITOPD,IPODD)          
  87.         ENDIF
  88.         IPOS=IPOX
  89.         SEGACT,IPOS
  90.         NBNNOE=IPOS(/1)
  91.         NBRINC=IPOS(NBNNOE)
  92.       ELSE
  93.         MMATRI = IPMATR
  94.         SEGACT,MMATRI
  95.         MILIGN = IILIGN
  96.         SEGACT,MILIGN
  97.         NBRINC = IPNO(/1)
  98.         SEGDES,MILIGN
  99.         SEGDES,MMATRI
  100.       ENDIF
  101. C
  102.       END
  103.  
  104.  
  105.  

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