DATE     94/01/10

     CETTE PROCEDURE A ETE MISE GRACIEUSEMENT
    A DISPOSITION DE LA COMMUNAUTE  CASTEM2000
  PAR  DELERUYELLE F.   (SOCOTEC-INDUSTRIE à L IPSN/DES)

                                                                         
      Procédure @CDG
      --------------
                                                                         
      XG1 XG2 (XG3)  =  @CDG  MODL1 (CAR1) ;
                                                                         

      Objet :                                                            
      ______
                                                                         
      Cette  procédure calcule les coordonnées du centre de gravité de
      la géométrie contenue dans un modèle.
                                                                         
                                                                         
      Commentaires :
      _____________
                                                                         
      MODL1  : objet modèle (type MMODEL).                              
                                                                         
      CAR1   : champ de caractéristiques géométriques (facultatif)
               (type MCHAML,sous-type CARACTERISTIQUES).
               Si on veut prendre en compte une distribution de masse
               volumique non uniforme, il faut entrer dans CAR1 le champ
               par élément de MATERIAU (pour donner RHO) et le champ de
               CARACTERISTIQUES (pour donner EPAI, SECT, ...).
                                                                         
      XG1    : première coordonnée du centre de gravité (type FLOTTANT).
                                                                         
      XG2    : seconde coordonnée du centre de gravité (type FLOTTANT).
                                                                         
      XG3    : éventuellement, troisième coordonnée du centre de gravité
               en 3 D (type FLOTTANT).
                                                                         
                                                                         
      Exemple :
      ________

      mo1 = mode s3d mecanique elastique isotrope coq3 coq4 ;
      ca1 = cara mo1 epai 0.5 ;
      xg1 xg2 xg3 = @cdg mo1 ca1 ;
      list xg1 ; list xg2 ; list xg3 ;

      Remarques :
      __________

      1) Cette procédure n'a été testée que sur des modèles à formulation
         MECANIQUE.

      2) Pour les coques, le domaine d'intégration est la surface de la
         coque et pour les poutres, le domaine d'intégration est la ligne
         moyenne de la poutre. Si l'on veut intégrer sur le volume de ces
         éléments, il faut donner le champs de caractéristiques géométri-
         ques CAR1 (type MCHAML, sous-type CARACTERISTIQUES).

      3) Dans le cas des coques, même en fournissant CAR1, on ne prend pas
         en compte l'excentrement.

      4) On peut traiter un modèle mécanique contenant plusieurs types
         d'éléments.