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    FANDEUR   22/01/03    21:15:44     11237          
  2. C RIGMOD    SOURCE    PETI      88/08/12    21:14:26
  3.       SUBROUTINE RIGMOD(IMODE,IRIG,IRET)
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8 (A-H,O-Z)
  6. C
  7. C=======================================================================
  8. C  CE SUBROUTINE CALCULE POUR LES MODES IMODE LES PETITES MATRICES
  9. C  RIGIDITE CONTENANT LA MASSE (IRIG=1) OU LA RAIDEUR (IRIG=2)GENERALISE
  10. C  CES MATRICES SONT ASSOCIEES A L ELEMENT QUI CONTIENT LE POINT QUI
  11. C  EST L INDICE DE CHAQUE MODE . LE RESULTAT EST MIS DANS IRET .
  12. C  IRIG=3 PREPARATION DE LA MATRICE D'AMORTISSEMENTS MODAUX
  13. C
  14. C  OPERATEUR AVEC LA SYNTAXE : K = RIGI M;
  15. C                    K : OBJET DE TYPE RIGIDITE
  16. C                    M : OBJET DE TYPE SOLUTION( SOUS TYPE MODE)
  17. C
  18. C  PROGRAMME PAR FARVACQUE
  19. C  APPELE PAR RGBASE RIGI
  20. C  N'APPELLE RIEN.
  21. C=======================================================================
  22. C
  23. C
  24. -INC SMRIGID
  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.       SEGINI xMATRI
  52.       ITYPEL=1
  53.       DO 3 IR=1,NELRIG
  54.       MMODE=ISOLEN(IR)
  55.       SEGACT MMODE
  56. *      SEGINI XMATRI
  57.       IF(IRIG.NE.1) GO TO 4
  58.       RE(1,1,ir)=FMMODD(2)
  59.       GO TO 10
  60.  4    CONTINUE
  61.       OMEG=2.*XPI*FMMODD(1)
  62.       IF(IRIG.EQ.3) GOTO 20
  63.       OMEG=OMEG*OMEG
  64.       RE(1,1,ir)=FMMODD(2)*OMEG
  65.       GOTO 10
  66.    20 RE(1,1,ir)=FMMODD(2)*OMEG*2.
  67.  10   CONTINUE
  68.       SEGDES MMODE
  69. *      IMATTT(IR)=XMATRI
  70.  3    CONTINUE
  71.       SEGDES MSOLUT,xMATRI,MSOLEN
  72.       NLIGRP=1
  73.       NLIGRD=1
  74.       SEGINI DESCR
  75.       NOELEP(1)=1
  76.       NOELED(1)=1
  77.       LISINC(1)='ALFA'
  78.       LISDUA(1)='FALF'
  79.       SEGDES DESCR
  80.       NRIGE=6
  81.       NRIGEL=1
  82.       SEGINI MRIGID
  83.       IFORIG=IFOUR
  84.       COERIG(1)=1.D0
  85.       IMGEO1=0
  86.       IMGEO2=0
  87.       ICHOLE=0
  88.       IRIGEL(1,1)=MELEME
  89.       IRIGEL(2,1)=0
  90.       IRIGEL(3,1)=DESCR
  91.       IRIGEL(4,1)=xMATRI
  92.       IRIGEL(5,1)=NIFOUR
  93.       IRIGEL(6,1)=0
  94.       IF(IRIG.NE.1) GO TO 11
  95.       MTYMAT='MASSE'
  96.       GOTO 13
  97.   11  IF(IRIG.NE.2) GO TO 12
  98.       MTYMAT='RIGIDITE'
  99.       GOTO 13
  100.    12 MTYMAT='AMORMODA'
  101.    13 CONTINUE
  102.       SEGDES MRIGID
  103.       SEGDES MELEME
  104.       IRET=MRIGID
  105.  5000 CONTINUE
  106.       RETURN
  107.       END
  108.  
  109.  
  110.  
  111.  

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