Notices | Cast3M

$$$$ POT_SCAL NOTICE  CHAT      11/09/12    21:17:41     7124           
                                             DATE     11/09/12

  Procedure  POT_SCAL                  Voir aussi :
        -------------------  
         POT_SCAL  TABB  SOLIN

 
        Objet : 
        -------
        En magnetostatique 3D calcul par la methode a deux potentiels
        le potentiel total et le potentiel reduit.
        Materiaux isotropes ou isotropes transverse .
   
        Commentaires :
        --------------

        TABB     table contenant les indices suivants en toute lettre

        'DPHI'  =  GEO1 maillage  de la zone de potentiel reduit
        
  eventuellement:       
        'RIGCON' =   RIG1 rigidite construite  comme suit:
  
               | MOD1 = GEO1 'MODE' 'THERMIQUE' 'ISOTROPE' 
               | RIG1 = 'CONDUC' MOD1  'K'  MUAIR
                  RIG1 peut  l'assemblage  de rigidites elementaires 
                  et de  rigiditees de super elements  
          (  si absente de la table cette conductivite sera construite 
           - sans super elements) 
          
        'SEPPHI'=  GEO2 surface de separation (appartenant a  DPHI) 
        'ORIG'  =  POINT de la surface de separartion ou on impose
                   l'egalite des potentiels 
                   
                   
        'POTSYM'= TAB1 
                  avec TAB1 table contenant des tables de description 
                  sur SEPPHI et son entourage eventuel pour calcul 
                  et application du saut de potentiel

        TAB1.I  = TABLE STN   (I = 1,2.........N)
                  une table STN  permet de decrire differentes  formes 
         de conditions limites pouvant se presenter  ou de passer
         directement des objets de type rigidite a prendre en compte
         dans la resolution ( conditions de periodicite par exemple)
         3 formes disponibles :

          1)    |STN.'LGEO' = GEO3 maillage type ligne( trace de SEPPHI
                |            dans un plan de symetrie des potentiels 
                |             si elle existe .En general le saut de 
                |             potentiel est nul si on a choisi
                |             ORIG sur cette ligne )
                |STN.'MTYP' = MOT1   'TBLOQ' OU 'RENSE'  condition 
                    s'appliquant aussi a la surface GEO4
                    si TBLOQ  blocage a zero des noeuds de GEO3
                    si RENSE  application relation d ensemble sur GEO3
        eventuellement ( si SEPPHI a plusieurs traces dans un plan de
        symetrie des potentiels) : 

                |STN.'SGEO'  = GEO4 maillage de type surface (surface
                |           equipotentielle liee au saut de potentiel)
                |            GEO3  appartient a son son contour 
                            
                            
           2)   |STN.'IMPOSE' = CHPO chpoint du saut impose localement
                                                           
                            
                           
        'MUAIR' =  FLOT1   permeablite du vide 



        'DPSI'   = GEO4   maillage de la zone de potentiel total
                  (totalite du domaine  air et ferro-magnetique) 
         
        |'AIRPSI' = GEO5   maillage de la zone d' air de DPSI  
        |                      (appartient a DPSI)                              
  OU    |'RIGCSPSI' =   RIG2 rigidite construite  comme suit:
  
               | MOD2 = GEO5 'MODE' 'THERMIQUE' 'ISOTROPE' 
               | RIG2 = 'CONDUC' MOD2  'K'  MUAIR
                  RIG2 peut  l'assemblage  de rigidites elementaires 
                  et de  rigiditees de super elements              
                
   
   
        
        'TABNUSEC' = TABMAT   table stockage de l'ensemble 
              des materiaux  ferro-magnetiques  

             | MODI = MODE GEOI THERMIQUE ISOTROPE (ORTHOTROPE)                 
             |  TABMAT.MODI = STN ;                                    
             |  STN.'EV1'  = evolution MU1(H)  sortie  voir H_B        
  Imat fois  |    si materiau  isotrope transverse (isotrope plan 23)           
             |  STN.'EV2'  = evolution MU2(H)  sortie  voir H_B        
             |  STN.'DIR1' = P1  direction 1 (type point)                       
             |  STN.'DIR2' = P2  direction 2 (type point)
                                                                      
            
                     
        'BLOQUE' = conditions limites du probleme global        
        'BIOT' = CHP1  champ D INDUCTION defini au minimum  sur 
          la surface de separation SEPPHI ( contruit avec l operateur
          BIOT par exemple  attention c'est H et non B )
          
        'ISTEP'  entier 1   le calcul sera fait en une seule fois
                   ou     2   le calcul est fait en deux temps :
                 
              A) cacul du champ d induction   la TABB doit etre
               sauvegardee ( voir SAUVER)
              B) le complement du calcul est lance apres restitution 
              de la sauvegarde en etrant dans POT_SCAL puis MAG_NLIN 
                                                              
   ---------------------------------------------------------------              
   SOLIN    si present  on calcule  la solution  en lineaire 
            Pour un calcul non lineaire apres passage dans POT_SCAL 
            il faut enchainer par la procedure MAGNLIN 
                      
            si absent le premier pas sera  fait dans MAG_NLIN 
            voir ci apres          
   
        Cette procedure construit et archive dans TABB les ingredients
        necessaire a un passage en non lineaire dans MAG_NLIN.                  


        SORTIES :

        le potentiel est dans TABB.'POTENTIEL' 
           POTENTIEL REDUIT  ET POTENTIEL TOTAL  SUR LES ZONES 
           CORRESPONDANTES.
           
           SI ON S INTERESSE A L INDUCTION  VOIR PROCEDURE MKGDT