Esope source files

   1 : C KGRAV2    SOURCE    CHAT      09/10/09    21:19:23     6519           
   2 :          
   3 :       SUBROUTINE KGRAV2(IPMAIL,XRG,IPOIN1,IPMATR,IDISS,IPMINT)
   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 COQ4         
  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 : *        ipmint:  pointeur sur le segment d'integration             
  25 : *                                                                     *         
  26 : *   SORTIES :                                                         *         
  27 : *   ________
  28 : *        on remplit les segments xmatri et imatri pour la sous zone                                                             
  29 : *                                                                                  
  30 : *                                                                     *         
  31 : *---------------------------------------------------------------------*         
  32 :       IMPLICIT INTEGER(I-N)
  33 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  34 : -INC CCOPTIO                                                                                                                                        
  35 : -INC SMRIGID                                                                                                                                      
  36 : -INC SMELEME                                                                    
  37 : -INC SMCOORD 
  38 :                                                                                                                                                                                                    
  39 :                                                                               
  40 :       REAL*8 XE(3,4)                                            
  41 :                                                                                                   
  42 :                                                                               
  43 :       MELEME=IPMAIL                                                             
  44 :       NBNN=NUM(/1)                                                              
  45 :       NBELEM=NUM(/2)                                                                                                                       
  46 : *                                                                               
  47 :       xMATRI=IPMATR                                                             
  48 :       NLIGRP=24 
  49 :       NLIGRD=24                                                                                     
  50 :                                                                                                                                                                                                                                          
  51 : C                                                                               
  52 : C     BOUCLE DE CALCUL POUR LES DIFFERENTS ELEMENTS                             
  53 : C                                                                               
  54 :       DO 3027 IB=1,NBELEM                                                       
  55 : C                                                                               
  56 : C     ON CHERCHE LES COORDONNEES DE L ELEMENT IB                                
  57 : C                                                                               
  57 :          CALL DOXE(XCOOR,IDIM,NBNN,NUM,IB,XE)                                                                           
  58 :                                                                  
  59 : *      SEGINI XMATRI                                                             
  60 : *      IMATTT(IB)=XMATRI 
  61 : *      IXMATR = XMATRI                                                                                 
  61 :       CALL CQ4GR1(XE,XRG,IPOIN1,IDISS,xmatri,IPMINT,ib)                                                                                                            
  62 : *      SEGDES XMATRI                                                             
  63 :  3027 CONTINUE                                                                  
  64 :       SEGDES xMATRI                                                                                                                                                                          
  65 :       RETURN                                                                    
  66 :       END                                                                       
  67 :  
  68 :