Télécharger la plaquette Cast3M


fr en

Télécharger SMMATRIK.INC

Retour à la liste

Numérotation des lignes : 

  1. C*******************************************************************            
  2. C                                                                               
  3. C    REPERAGE ET STOKAGE DES MATRICES ELEMENTAIRES  puis assemblees
  4. C    (description par sous-zone associees a un operateur)
  5. C                                                                               
  6. C   COEMTK(I)     : Flottant , coefficient multiplicateur de la matrice
  7. C                                                                               
  8. C   IRIGEL(1,I)   : POINTEUR SUR L'OBJET GEOMETRIE (Inconnues primales)
  9. C   IRIGEL(2,I)   : POINTEUR SUR L'OBJET GEOMETRIE (Inconnues duales)
  10. C   IRIGEL(3,I)   : Non utilise (POINTEUR SUR LE SEGMENT DESCRIPTIF D'UNE
  11. C                                MATRICE ELEMENTAIRE.(SEGMENT DESCR)
  12. C   IRIGEL(4,I)   : POINTEUR SUR LE SEGMENT CONTENANT LES POINTEURS
  13. C                        DES MATRICES DE MRIGIDITE DE CHAQUE ELEMENTS.
  14. C                        (SEGMENT IMATRI)
  15. C   IRIGEL(5,I)   : Non utilise
  16. C   IRIGEL(6,I)   : Non utilise
  17. C   IRIGEL(7,1)   : 0   LA MATRICE EST SYMETRIQUE
  18. C                 : 1   LA MATRICE EST ANTISYMETRIQUE
  19. C                 : 2   LA MATRICE EST NON SYMETRIQUE
  20. C                 : 3   LA MATRICE EST RECTANGULAIRE avec SPGP # SPGD
  21. C                 : 4   LA MATRICE EST type 3 et CCt (on ne stoke que C)
  22. C                 : 5   LA MATRICE EST diagonale
  23. C                 : 6   LA MATRICE EST deja assemblee en morse
  24.  
  25.       SEGMENT MATRIK
  26.       REAL*8  COEMTK(NMATRI)
  27.       INTEGER IRIGEL(NRIGE,NMATRI)
  28.  
  29. C KSIM  =0 matrice symetrique  =2 matrice non symetrique
  30. C KMINC , KMINCP , KMINCD : pointeur sur MINC repartition des inconnues
  31. C totales primales et duales                            PROFKS PKINC
  32. C KIZM  : pointeur sur les connectivites globales
  33. C KISPGT KISPGP KISPGD ; SPG assemble pour inc totales,prim et dua
  34. C KNTTT  KNTTP  KNTTD  ; nb d'inconnues total
  35.  
  36.       INTEGER KSYM,KMINC,KMINCP,KMINCD,KIZM
  37.       INTEGER KISPGT,KISPGP,KISPGD
  38.       INTEGER KNTTT,KNTTP,KNTTD
  39.  
  40. C KIDMAT: pointeur sur stokage bloc IDMAT (Cholevski)    TRIAKS
  41. C KS2B  : pointeur sur second membre(IZA)cree ds PROFKS calcule ds KASMBR
  42. C KMORS : pointeur sur profil Morse(PMORS)               ASSMT (KALMOR)
  43. C KISA  : pointeur sur stokage Morse(IZA)                ASSMT (KALMOR)
  44. C KMRST : pointeur sur profil Morse(PMORS)de AAt         PROFKS(KALMOR)
  45. C KIST  : pointeur sur stokage Morse(IZA) de AAt         PROFKS(KALMOR)
  46. C KCLIM : pointeur sur stokage C lim (CHPT)
  47. C KTRING: 0 pas triangulée 1 triangulée
  48.  
  49. C  nkid=9 :  IDMATP,IDMATD,KS2B,KMORS,KISA,KMRST,KIST,KCLIM,KTRING
  50.       INTEGER KIDMAT(NKID)
  51. C  nkmt=7 :  KMMT,MATRIU,MATRIP,IZDU,IZDP,IZFU,IZFP
  52. C                               (IZA)(IZA)(IZA)(IZA)
  53.       INTEGER KKMMT(NKMT)
  54.       ENDSEGMENT
  55.  
  56. C LIZAFM(NBSOUS,.) description par sous-objet geometrique -> IZAFM
  57. C KSPGP , KSPGD : SPG pour les inconnues primales et duales
  58.       SEGMENT IMATRI
  59.       CHARACTER*8 LISPRI(NBME),LISDUA(NBME)
  60.       INTEGER LIZAFM(NBSOUS,NBME)
  61.       INTEGER KSPGP,KSPGD
  62.       ENDSEGMENT
  63.  
  64. C Stokage matrices elementaires non assemblees (valeurs)
  65.       SEGMENT IZAFM                                                             
  66.       REAL*8  AM(NBEL,NP,MP)
  67.       ENDSEGMENT                                                                
  68.       POINTEUR IPM1.IZAFM,IPM2.IZAFM,IPM3.IZAFM,IPM4.IZAFM
  69.       POINTEUR IPM5.IZAFM,IPM6.IZAFM,IPM7.IZAFM,IPM8.IZAFM
  70.       POINTEUR IPM9.IZAFM
  71.  
  72. C Reperage des inconnues
  73. *  NPT nb de noeud  NBI nb de composantes total  NTT nb total de DDL
  74. *  MPOS(NPT,NBI+1) = 0 si l'inconnue j n est pas defini au noeud i
  75. *            sinon = k rang de cette inconnue pour le noeud i
  76. *  MPOS(i,NBI+1) nb d'inconnues au noeud i
  77. *  NPOS(NPT)     Position de la 1ere inconnue du noeud i
  78. *  NPOS et MPOS  sont donnes ds la numerotation optimisee
  79.  
  80.       SEGMENT MINC
  81.       CHARACTER*8 LISINC(NBI)
  82.       INTEGER     NPOS(NPT+1)
  83.       INTEGER     MPOS(NPT,NBI+1)
  84.       ENDSEGMENT
  85.       POINTEUR MINCP.MINC,MINCD.MINC
  86.  
  87.       SEGMENT PMORS
  88.       INTEGER IA (NTT+1)
  89.       INTEGER JA (NJA)
  90.       ENDSEGMENT
  91.       POINTEUR PMS1.PMORS,PMS2.PMORS
  92.  
  93. C Segment de stokage
  94.       SEGMENT IZA
  95.       REAL*8  A(NBVA)
  96.       ENDSEGMENT
  97.       POINTEUR IZA1.IZA,IZA2.IZA,IZAU.IZA,IZAL.IZA,ISA.IZA
  98.  
  99.       SEGMENT IDMAT
  100. *  KZA(NTT)    Longueur de chaque ligne de la matrice (diag comprise)
  101. *  NUIA(NTT,2) 1/ numero du bloc dans lequel se trouve la ligne i de la
  102. *                 matrice
  103. *              2/ position du debut de la ligne dans le segment IZA - 1
  104.       INTEGER KZA(NTT),NUIA(NTT,2)
  105.       INTEGER NUAN(NPT),NUNA(NPT)
  106.       INTEGER IDIAG
  107.       INTEGER IDESCL(NBLK)
  108.       INTEGER IDESCU(NBLK)
  109.       INTEGER NLDBLK(NBLK+1)
  110.       ENDSEGMENT
  111.  
  112. C*******************************************************************            
  113.  

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