Accueil | Cast3M

   1 : $$$$ POT_SCAL NOTICE  CHAT      11/09/12    21:17:41     7124           
   2 :                                              DATE     11/09/12
   3 : 
   4 :    Procedure  POT_SCAL                  Voir aussi :
   5 :         -------------------  
   6 :          POT_SCAL  TABB  SOLIN
   7 : 
   8 :  
   9 :         Objet : 
  10 :         -------
  11 :         En magnetostatique 3D calcul par la methode a deux potentiels
  12 :         le potentiel total et le potentiel reduit.
  13 :         Materiaux isotropes ou isotropes transverse .
  14 :    
  15 :         Commentaires :
  16 :         --------------
  17 : 
  18 :         TABB     table contenant les indices suivants en toute lettre
  19 : 
  20 :         'DPHI'  =  GEO1 maillage  de la zone de potentiel reduit
  21 :         
  22 :   eventuellement:       
  23 :         'RIGCON' =   RIG1 rigidite construite  comme suit:
  24 :   
  25 :                | MOD1 = GEO1 'MODE' 'THERMIQUE' 'ISOTROPE' 
  26 :                | RIG1 = 'CONDUC' MOD1  'K'  MUAIR
  27 :                   RIG1 peut  l'assemblage  de rigidites elementaires 
  28 :                   et de  rigiditees de super elements  
  29 :           (  si absente de la table cette conductivite sera construite 
  30 :            - sans super elements) 
  31 :           
  32 :         'SEPPHI'=  GEO2 surface de separation (appartenant a  DPHI) 
  33 :         'ORIG'  =  POINT de la surface de separartion ou on impose
  34 :                    l'egalite des potentiels 
  35 :                    
  36 :                    
  37 :         'POTSYM'= TAB1 
  38 :                   avec TAB1 table contenant des tables de description 
  39 :                   sur SEPPHI et son entourage eventuel pour calcul 
  40 :                   et application du saut de potentiel
  41 : 
  42 :         TAB1.I  = TABLE STN   (I = 1,2.........N)
  43 :                   une table STN  permet de decrire differentes  formes 
  44 :          de conditions limites pouvant se presenter  ou de passer
  45 :          directement des objets de type rigidite a prendre en compte
  46 :          dans la resolution ( conditions de periodicite par exemple)
  47 :          3 formes disponibles :
  48 : 
  49 :           1)    |STN.'LGEO' = GEO3 maillage type ligne( trace de SEPPHI
  50 :                 |            dans un plan de symetrie des potentiels 
  51 :                 |             si elle existe .En general le saut de 
  52 :                 |             potentiel est nul si on a choisi
  53 :                 |             ORIG sur cette ligne )
  54 :                 |STN.'MTYP' = MOT1   'TBLOQ' OU 'RENSE'  condition 
  55 :                     s'appliquant aussi a la surface GEO4
  56 :                     si TBLOQ  blocage a zero des noeuds de GEO3
  57 :                     si RENSE  application relation d ensemble sur GEO3
  58 :         eventuellement ( si SEPPHI a plusieurs traces dans un plan de
  59 :         symetrie des potentiels) : 
  60 : 
  61 :                 |STN.'SGEO'  = GEO4 maillage de type surface (surface
  62 :                 |           equipotentielle liee au saut de potentiel)
  63 :                 |            GEO3  appartient a son son contour 
  64 :                             
  65 :                             
  66 :            2)   |STN.'IMPOSE' = CHPO chpoint du saut impose localement
  67 :                                                            
  68 :                             
  69 :                            
  70 :         'MUAIR' =  FLOT1   permeablite du vide 
  71 : 
  72 : 
  73 : 
  74 :         'DPSI'   = GEO4   maillage de la zone de potentiel total
  75 :                   (totalite du domaine  air et ferro-magnetique) 
  76 :          
  77 :         |'AIRPSI' = GEO5   maillage de la zone d' air de DPSI  
  78 :         |                      (appartient a DPSI)                              
  79 :   OU    |'RIGCSPSI' =   RIG2 rigidite construite  comme suit:
  80 :   
  81 :                | MOD2 = GEO5 'MODE' 'THERMIQUE' 'ISOTROPE' 
  82 :                | RIG2 = 'CONDUC' MOD2  'K'  MUAIR
  83 :                   RIG2 peut  l'assemblage  de rigidites elementaires 
  84 :                   et de  rigiditees de super elements              
  85 :                 
  86 :    
  87 :    
  88 :         
  89 :         'TABNUSEC' = TABMAT   table stockage de l'ensemble 
  90 :               des materiaux  ferro-magnetiques  
  91 : 
  92 :              | MODI = MODE GEOI THERMIQUE ISOTROPE (ORTHOTROPE)                 
  93 :              |  TABMAT.MODI = STN ;                                    
  94 :              |  STN.'EV1'  = evolution MU1(H)  sortie  voir H_B        
  95 :   Imat fois  |    si materiau  isotrope transverse (isotrope plan 23)           
  96 :              |  STN.'EV2'  = evolution MU2(H)  sortie  voir H_B        
  97 :              |  STN.'DIR1' = P1  direction 1 (type point)                       
  98 :              |  STN.'DIR2' = P2  direction 2 (type point)
  99 :                                                                       
 100 :             
 101 :                      
 102 :         'BLOQUE' = conditions limites du probleme global        
 103 :         'BIOT' = CHP1  champ D INDUCTION defini au minimum  sur 
 104 :           la surface de separation SEPPHI ( contruit avec l operateur
 105 :           BIOT par exemple  attention c'est H et non B )
 106 :           
 107 :         'ISTEP'  entier 1   le calcul sera fait en une seule fois
 108 :                    ou     2   le calcul est fait en deux temps :
 109 :                  
 110 :               A) cacul du champ d induction   la TABB doit etre
 111 :                sauvegardee ( voir SAUVER)
 112 :               B) le complement du calcul est lance apres restitution 
 113 :               de la sauvegarde en etrant dans POT_SCAL puis MAG_NLIN 
 114 :                                                               
 115 :    ---------------------------------------------------------------              
 116 :    SOLIN    si present  on calcule  la solution  en lineaire 
 117 :             Pour un calcul non lineaire apres passage dans POT_SCAL 
 118 :             il faut enchainer par la procedure MAGNLIN 
 119 :                       
 120 :             si absent le premier pas sera  fait dans MAG_NLIN 
 121 :             voir ci apres          
 122 :    
 123 :         Cette procedure construit et archive dans TABB les ingredients
 124 :         necessaire a un passage en non lineaire dans MAG_NLIN.                  
 125 : 
 126 : 
 127 :         SORTIES :
 128 : 
 129 :         le potentiel est dans TABB.'POTENTIEL' 
 130 :            POTENTIEL REDUIT  ET POTENTIEL TOTAL  SUR LES ZONES 
 131 :            CORRESPONDANTES.
 132 :            
 133 :            SI ON S INTERESSE A L INDUCTION  VOIR PROCEDURE MKGDT
 134 : 
 135 :         
 136 :         
 137 :