Manual pages | 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

 
        Description : 
        -------------
        In 3D magnetostatics, calculation of both the total and
        the reduced potential by the two-potential method.           
        Isoptropic or transversely isotropic materials                          
        Contents :
        ----------

        TABB     table containing the following indices written in
                 words

        'DPHI'  =  GEO1 mesh of the reduced potential region  
  optional :    
        'RIGCON' =   RIG1  stiffness of reduced potential zone built
          as follows :  
  
               | MOD1 = GEO1 'MODE' 'THERMIQUE' 'ISOTROPE' 
               | RIG1 = 'CONDUC' MOD2  'K'  MU0
               
            RIG1 may be assembly of standard and super element matrix 
          if not given it will be built in the procedure ) as a normal
          stiffness)
                  
        'SEPPHI'=  GEO2 separating surface (pertaining to DPHI)  
                   and directed according to the outside normal
        'ORIG'  =  POINT of the separating surface where reduced and
                   total potentials are requested to be equal
                   
        'POTSYM'= TAB1 
                  with TAB1 table containing tables describing
                  SEPPHI and its surroundings for the calculation
                  and the applying of the potential jump

        TAB1.I  = TABLE STN ;
                different forms of boundary conditions may be entered
         by giving support and type of condition or directly by giving
         the stiffness.   
                |STN.'LGEO' = GEO3 line type mesh (trace of SEPPHI
        1)      |             in a plane of symmetry of the potentials 
                |             if the trace exists. The potential jump 
                |             is null if ORIG was chosen on this line)  
                |             
                |STN.'MTYP' = MOT1   'TBLOQ' OR 'RENSE' condition 
                             also applying to the surface GEO4
                         if TBLOQ  jump put to zero on GEO3 nodes 
                         if RENSE  identite relation put on GEO3 nodes 

        (if SEPPHI has several traces in a plane of symmetry of the
        potentials) : 

                |STN.'SGEO  = GEO4 surface type mesh (equipotential
                |             surface connected with the potential
                |             jump)
                |            GEO3 is a section of its contour 
                
         2)     |STN.'IMPOSE' = CHPO chpoint of imposed potential jump
                 
                            
                            
        Repeat the sequence for I = 1 at N if necessary.

        'MUAIR' =  FLOT1   vacuum permeability
        
        'DPSI'  = GEO4  entire mesh of total potential domain 
         
        |'AIRPSI' = GEO5   mesh of air  of total potential domain  
        |                         ( part of GEO4 )                              
  OU    |'RIGCSPSI' =   RIG2 stiffness of  GEO5 built as follows:
  
               | MOD2 = GEO5 'MODE' 'THERMIQUE' 'ISOTROPE' 
               | RIG2 = 'CONDUC' MOD2  'K'  MUAIR
           RIG2 may be assembly of standard and super element matrix 
          ( if not given it will be built in the procedure ) as a
          normal stiffness)
                
   
   
        
        'TABNUSEC' = TABMAT   table containing description 
                ferro-magnetic domain(s) built as follows :

             |MODI = MODE GEOI THERMIQUE ISOTROPE (ORTHOTROPE)                  
             |TABMAT.MODI = STN ;                                    
             |STN.'EV1'  = evolution MU1(H)  sortie  voir H_B        
  Imat fois  |    if transversely isotropic material 
             |STN.'EV1'  =  MU2(H) direction 1                   
             |STN.'EV2'  =  MU3(H) directions  2 and 3 
             |STN.'DIR1'  = P1 (point type )      
             |STN.'DIR2'  = P2 (point type )  P2 orthogonal to P1 
             
       Repeat the sequence for I = 1 at N  material if necessary.
                                

               
        'BLOQUE' =  general boundary condition 
        
        'BIOT' = CHP1  contains the coil field on DPHI defined 
                     at least on SEPPHI  (can be built with operator
                     BIOT . pay attention that you must pass the H 
                     field  and not B field)    

         
         'SOLIN'   in the presence of this optional word, the linear
                   solution will be calculated straightaway, otherwise
                   it will be calculated in MAG_NLIN.
                   
                   
         'ISTEP'   Integer 1  computation done in a single run
                        or 2  computation in 2 steps :
                         a)  (first computation of H field followed 
                by a SAUVER TABB ;
                         b)  make a RESTI and enter in POT_SCAL and
                MAG_NLIN            
                              
        output data :

        The potential is in TABB.'POTENTIEL'.
        This procedure constructs and stores the data required  
        for a non linear running into MAG_NLIN in TABB.