Esope source files

   1 : C KGRAV3    SOURCE    CHAT      09/10/09    21:19:24     6519           
   2 :          
   3 :       SUBROUTINE KGRAV3(IPMAIL,XRG,IPOIN1,IPMATR,IDISS)
   4 : *---------------------------------------------------------------------*         
   5 : *      APPELEE PAR KGRAV1         
   6 : *                                                                     *         
   7 : * calcule la matrice de rigidité liée à la variation de pression à cause
   8 : * d'un mouvement suivant le vecteur IPOIN1.     
   9 : *
  10 : *  éléments COQ3,DKT,DST         
  11 : *                                                                     *         
  12 : *---------------------------------------------------------------------*         
  13 : *                                                                     *         
  14 : *   ENTREES :                                                         *         
  15 : *   ________                                                          *         
  16 : *                                                                     *         
  17 : *        IPMAIL   Pointeur sur un segment  MELEME                     *                             
  18 : *        XRG      coefficient multiplicateur de la rigidité           *
  19 : *        IPOIN1   vecteur(type point), definissant le sens de gravité *
  20 : *        IPMATR   pointeur sur la matrice de rigidité de pesanteur
  21 : *                 de la sous-zone
  22 : *        IDISS : 0 ---> MATRICE SYMETRIQUE
  23 : *                1 ---> MATRICE DISSYMETRIQUE              
  24 : *                                                                     *         
  25 : *   SORTIES :                                                         *         
  26 : *   ________
  27 : *        on remplit les segments xmatri et imatri pour la sous zone                                                             
  28 : *                                                                                  
  29 : *                                                                     *         
  30 : *---------------------------------------------------------------------*         
  31 :       IMPLICIT INTEGER(I-N)
  32 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  33 : -INC CCOPTIO                                                                                                                                        
  34 : -INC SMRIGID                                                                                                                                      
  35 : -INC SMELEME                                                                    
  36 : -INC SMCOORD                                                                                                                                                                                                    
  37 :                                                                               
  38 :       REAL*8 XE(3,3)                                            
  39 :                                                                                                   
  40 :                                                                               
  41 :       MELEME=IPMAIL                                                             
  42 :       NBNN=NUM(/1)                                                              
  43 :       NBELEM=NUM(/2)                                                                                                                       
  44 : *                                                                               
  45 :       xMATRI=IPMATR                                                             
  46 :       NLIGRP=18 
  47 :       NLIGRD=18                                                                                     
  48 :                                                                                                                                                                                                                                          
  49 : C                                                                               
  50 : C     BOUCLE DE CALCUL POUR LES DIFFERENTS ELEMENTS                             
  51 : C                                                                               
  52 :       DO 3027 IB=1,NBELEM                                                       
  53 : C                                                                               
  54 : C     ON CHERCHE LES COORDONNEES DE L ELEMENT IB                                
  55 : C                                                                               
  55 :          CALL DOXE(XCOOR,IDIM,NBNN,NUM,IB,XE)                                                                           
  56 :                                                                  
  57 : *      SEGINI XMATRI                                                             
  58 : *      IMATTT(IB)=XMATRI 
  59 :       IXMATR = XMATRI                                                                                 
  59 :       CALL DKTGR1(XE,XRG,IPOIN1,IDISS,IXMATR,ib)                                                                                                            
  60 : *      SEGDES XMATRI                                                             
  61 :  3027 CONTINUE                                                                  
  62 :       SEGDES xMATRI                                                                                                                                                                          
  63 :       RETURN                                                                    
  64 :       END                                                                       
  65 :  
  66 :