Télécharger la plaquette Cast3M


fr en

Télécharger rigmod.eso

Retour à la liste

Numérotation des lignes : 

rigmod

  1. C RIGMOD    SOURCE    PV090527  26/04/30    21:16:21     12529          
  2.       SUBROUTINE RIGMOD(IMODE,IRIG,IRET)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8 (A-H,O-Z)
  5. C
  6. C=======================================================================
  7. C  CE SUBROUTINE CALCULE POUR LES MODES IMODE LES PETITES MATRICES
  8. C  RIGIDITE CONTENANT LA MASSE (IRIG=1) OU LA RAIDEUR (IRIG=2)GENERALISE
  9. C  CES MATRICES SONT ASSOCIEES A L ELEMENT QUI CONTIENT LE POINT QUI
  10. C  EST L INDICE DE CHAQUE MODE . LE RESULTAT EST MIS DANS IRET .
  11. C  IRIG=3 PREPARATION DE LA MATRICE D'AMORTISSEMENTS MODAUX
  12. C
  13. C  OPERATEUR AVEC LA SYNTAXE : K = RIGI M;
  14. C                    K : OBJET DE TYPE RIGIDITE
  15. C                    M : OBJET DE TYPE SOLUTION( SOUS TYPE MODE)
  16. C
  17. C  PROGRAMME PAR FARVACQUE
  18. C  APPELE PAR RGBASE RIGI
  19. C  N'APPELLE RIEN.
  20. C=======================================================================
  21. C
  22. C
  23. -INC SMRIGID
  24. -INC SMCOORD
  25. -INC SMELEME
  26. -INC SMSOLUT
  27. -INC CCREEL
  28. *-
  29.  
  30. -INC PPARAM
  31. -INC CCOPTIO
  32. C
  33.       IRET=0
  34.       IF(IRIG.NE.1.AND.IRIG.NE.2.AND.IRIG.NE.3) GOTO 5000
  35.       MSOLUT=IMODE
  36.       IF(MSOLUT.EQ.0) GO TO 5000
  37.       SEGACT MSOLUT
  38.       NIPO=MSOLIS(/1)
  39.       MSOLEN=MSOLIS(4)
  40.       MELEME=MSOLIS(3)
  41.       SEGACT MSOLEN
  42.       SEGACT MELEME*MOD
  43.       LTAB=ISOLEN(/1)
  44. C
  45. C  **** NBELEM=NBRE DE POINTS,DE RIGIDITE,D ELEMENTS
  46. C  **** FABRICATION DES PETITES MATRICES
  47. C
  48.       NELRIG=LTAB
  49.       NLIGRP=1
  50.       NLIGRD=1
  51.       rigrel=0
  52.       SEGINI xMATRI
  53.       ITYPEL=1
  54.       DO 3 IR=1,NELRIG
  55.       MMODE=ISOLEN(IR)
  56.       SEGACT MMODE
  57. *      SEGINI XMATRI
  58.       IF(IRIG.NE.1) GO TO 4
  59.       RE(1,1,ir)=FMMODD(2)
  60.       GO TO 10
  61.  4    CONTINUE
  62.       OMEG=2.*XPI*FMMODD(1)
  63.       IF(IRIG.EQ.3) GOTO 20
  64.       OMEG=OMEG*OMEG
  65.       RE(1,1,ir)=FMMODD(2)*OMEG
  66.       GOTO 10
  67.    20 RE(1,1,ir)=FMMODD(2)*OMEG*2.
  68.  10   CONTINUE
  69.       SEGDES MMODE
  70. *      IMATTT(IR)=XMATRI
  71.  3    CONTINUE
  72.       SEGDES MSOLUT,xMATRI,MSOLEN
  73.       NLIGRP=1
  74.       NLIGRD=1
  75.       SEGINI DESCR
  76.       NOELEP(1)=1
  77.       NOELED(1)=1
  78.       LISINC(1)='ALFA'
  79.       LISDUA(1)='FALF'
  80.       SEGDES DESCR
  81.       NRIGE=6
  82.       NRIGEL=1
  83.       SEGINI MRIGID
  84.       IFORIG=IFOUR
  85.       COERIG(1)=1.D0
  86.       IMGEO1=0
  87.       IMGEO2=0
  88.       ICHOLE=0
  89.       IRIGEL(1,1)=MELEME
  90.       IRIGEL(2,1)=0
  91.       IRIGEL(3,1)=DESCR
  92.       IRIGEL(4,1)=xMATRI
  93.       IRIGEL(5,1)=NIFOUR
  94.       IRIGEL(6,1)=0
  95.       IF(IRIG.NE.1) GO TO 11
  96.       MTYMAT='MASSE'
  97.       GOTO 13
  98.   11  IF(IRIG.NE.2) GO TO 12
  99.       MTYMAT='RIGIDITE'
  100.       GOTO 13
  101.    12 MTYMAT='AMORMODA'
  102.    13 CONTINUE
  103.       SEGDES MRIGID
  104.       SEGDES MELEME
  105.       IRET=MRIGID
  106.  5000 CONTINUE
  107.       RETURN
  108.       END
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  

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