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    CB215821  25/04/23    21:15:43     12247          
  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 SMCOORD
  26. -INC SMELEME
  27. -INC SMSOLUT
  28. -INC CCREEL
  29. *-
  30.  
  31. -INC PPARAM
  32. -INC CCOPTIO
  33. C
  34.       IRET=0
  35.       IF(IRIG.NE.1.AND.IRIG.NE.2.AND.IRIG.NE.3) GOTO 5000
  36.       MSOLUT=IMODE
  37.       IF(MSOLUT.EQ.0) GO TO 5000
  38.       SEGACT MSOLUT
  39.       NIPO=MSOLIS(/1)
  40.       MSOLEN=MSOLIS(4)
  41.       MELEME=MSOLIS(3)
  42.       SEGACT MSOLEN
  43.       SEGACT MELEME*MOD
  44.       LTAB=ISOLEN(/1)
  45. C
  46. C  **** NBELEM=NBRE DE POINTS,DE RIGIDITE,D ELEMENTS
  47. C  **** FABRICATION DES PETITES MATRICES
  48. C
  49.       NELRIG=LTAB
  50.       NLIGRP=1
  51.       NLIGRD=1
  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.  

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