Télécharger la plaquette Cast3M


fr en

Télécharger c2d93.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  c2d93.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ***********************************************************************
  5. *
  6. *       2D   AXISYMMETRIC  MAGNETIC FIELD COMPUTATION        
  7. *
  8. *       Formulation  : VECTOR POTENTIAL 
  9. *
  10. *                      NON LINEAR MATERIAL
  11. *
  12. ***********************************************************************                        
  13.  OPTION DIME 2  ELEM  TRI6  COUL VERT echo 1 ;
  14.  GRAPH ='N' ;
  15. * ------------------------Mesh --------------------------------------
  16.  R1 =20. ; R2= 25. ; R3= 27. ; R4= 29. ; R5=39.;R6 = 130.;              
  17.  Z1= 2.5;  Z2=  5. ; Z3= 20. ; Z4= 130.;Z31=30.;Z32=50.;                 
  18.  DI1= 1.; DI2=5.;DI3=60. ;DI4= 10.;NET2= 6;                             
  19.  
  20.   DENS DI1 ;                                                            
  21.  OZ3= 0.  Z3; OZ2= 0.  Z2;OZ1= 0.  Z1;OZ4=0. Z4;                        
  22.  OO= 0.  0.;R0Z0= DI1  0 ;                                              
  23.  R1Z0=  R1  0. ;R2Z0= R2  0. ;R3Z0 = R3  0. ;R4Z0 = R4  0.;             
  24.  R0Z2= DI1 Z2 ; R4Z2 = R4 Z2 ;                                          
  25.   DENS DI2 ;                                                            
  26.  R5Z0=  R5  0. ;R5Z3= R5 Z3 ;R0Z3 =DI2 Z3 ;OZ3 = 0. Z3 ;                
  27.   DENS DI3 ;                                                            
  28.  R6Z0=R6  0.; R6Z4=R6 Z4 ;                                              
  29.   DENS 10 ;                                                             
  30.  R2Z3= R2  Z3 ; R1Z3= R1  Z3; R4Z1 = R4   Z1;                           
  31.  R3Z1= R3  Z1 ;                                                         
  32.  DENSITE DI4 ;                                                          
  33.  R0Z4 = DI4 Z4 ; OZ4 = 0. Z4 ;                                                                                      
  34.  NTRA= 1;                                                                                                                         
  35.  FLAN1 =( D OO R0Z0 D R1Z0 ) TRAN DINI DI1 DFIN DI1 ( OZ1 MOINS OO )    
  36.              TRAN DINI DI1 DFIN DI1 ( OZ2 MOINS OZ1) COUL VERT ;
  37. *-------------------------  COIL SURFACE -------------------------                   
  38.  BOBI  =( D R1Z0 R2Z0 ) TRAN DINI DI1 DFIN DI1 ( OZ1 MOINS OO )        
  39.                                                       COUL BLEU ;
  40. *-----------------------------------------------------------------                                                            
  41.  FLAN2 =(INVE( BOBI COTE 3)) TRAN DINI DI1 DFIN DI1                     
  42.                             ( OZ2 MOINS OZ1) COUL VERT ;                
  43.  FLAN3 =( D R2Z0 R3Z0 ) TRAN DINI DI1 DFIN DI1 ( OZ1 MOINS OO )         
  44.              TRAN DINI DI1 DFIN DI1 ( OZ2 MOINS OZ1) COUL VERT ;        
  45.  FLAN4 =( D R3Z0 R4Z0 ) TRAN DINI DI1 DFIN DI1 ( OZ1 MOINS OO )         
  46.              TRAN DINI DI1 DFIN DI1 ( OZ2 MOINS OZ1) COUL VERT ;        
  47.  
  48. *---------------------IRON ------------------------------------------                 
  49.  FER = (D OZ2 R0Z2 D R4Z2 D 6 R4Z0 D R5Z0 D R5Z3 D R0Z3 D OZ3 D OZ2 )
  50.  SURF PLANE COUL ROUG ;                                                
  51. *---------------------EXTERNAL AIR BOX ------------------------------                
  52.  FLAN5 = ( D   R5Z0 R5Z3  D R6Z4 D R6Z0 D R5Z0 ) SURF PLANE ;
  53. *                                                                                                                                               
  54.  FLAN6 =( D R0Z3  R5Z3 D R6Z4  D R0Z4 D R0Z3 ) SURF PLANE               
  55.                         COUL VERT ;                                     
  56.  TUB6 = (D OZ3 R0Z3 D R0Z4 D OZ4 D OZ3) SURF PLANE COUL vert ;
  57.  AIREXT = (FLAN6 ET TUB6 ET FLAN5 ) coul blan ;
  58.  AIRIN = ((FLAN1 ET FLAN2 ET FLAN3  ET FLAN4 ) coul vert) et BOBI  ;
  59.  TOUT = AIRIN ET AIREXT ET FER ; ELIM .2 TOUT ;
  60. *     
  61. *---------------  BOUNDARIES -----------------------------------------                   
  62. *                                                  
  63.  ENPP = CONTOUR TOUT  ;                                                 
  64.  CEXTR=ENPP POINTS DROITE   (130.  0.) (130. 100.)  .1 ;                
  65.  CEXTH=ENPP POINTS DROITE   ( 0. 130.) (10.  130.)  .1 ;                
  66.  AXE=ENPP POINTS DROITE   ( 0.  0. ) ( 0. 100. )  .1  ;
  67.  AX1 = (AIRIN ET FER) POINT DROITE ( 0.  0. ) ( 0. 100. )  .1  ;
  68.  AX2 = DIFF ( CHAN POI1 AXE) (CHAN POI1 AX1) ;
  69. *                   
  70. *----------------------------------------------------------------------
  71. *        FORMER DESCRIPTION IN  MILLIMETRES
  72. *---------------  SHIFT FOR METERS ------------------------------------
  73. *
  74.  deplacer tout homo  .001  (0. 0.) ;
  75. *
  76.  SI (NEG GRAPH  N ) ;
  77.  TITRE  ' MESH 2D ' ;
  78.  TRAC  tout  ;
  79.  FINSI ;
  80. *----------------------------------------------------------------------    
  81.  MU0= 4. * PI * 1.E-7 ;
  82. *  
  83. *------------------CURENT DESCRIPTION --------------------------------- 
  84.  TABCOUR= TABLE ;                                                       
  85.  DESCOUR TABCOUR 1 BOBI 'AMP' 800.E6 ;                                   
  86.  TABB = TABLE ;
  87. *----------------- AXISYMMETRIC  PROBLEM ------------------------------                                                          
  88.  TABB.'AXI'= VRAI  ;
  89.   NF = FER NBEL ;
  90.    FER1 = FER ELEM ( LECT 1 PAS 1 30 )  ;
  91.    FER2 = FER ELEM ( LECT 31 PAS 1 NF ) ;
  92. *------------------MATERIALS  -----------------------------------------   
  93.  KEVOL = H_B MU0    ;
  94.  TABMAT = TABLE ;                                                        
  95.     OBFER1 = FER1 MODE THERMIQUE ISOTROPE ;
  96.     OBFER2 = FER2 MODE THERMIQUE ISOTROPE ;  
  97.    STN = TABLE ; STN.EV1 = KEVOL   ;
  98.    TABMAT.OBFER1 = STN ;
  99.    STN = TABLE ; STN.EV1 = KEVOL   ;
  100.    TABMAT.OBFER2 = STN ;   
  101.  TABB.'TABNUSEC' = TABMAT ;
  102.  TABB.'MUAIR' =MU0 ;
  103. *------- LINEAR MATERIAL CAN BE TRAITED AS A SUPER ELEMENT OR NOT 
  104. *              
  105.    isuper= 1 ;
  106.  'SI' ('EGA' isuper 1 ) ;
  107.  TABB.'AIRSUP' = AIREXT ;
  108.  TABB.'MAITRES' = ( FER CONTOUR) 'ELEM' 'COMP' R5Z0 OZ3 ;
  109.  TABB.'ENCS' = CEXTR ET CEXTH ET AX2 ;
  110.  TABB.'BLOQUE' = BLOQUER 'T' AX1 ;
  111.  TABB.'AIR' = AIRIN ;
  112.  'SINON' ;
  113.  TABB.'AIR' = (AIRIN ET AIREXT) ;
  114.  TABB.'BLOQUE' = BLOQUER 'T' (CEXTR ET CEXTH  ET AXE ) ;
  115.  'FINSI' ;
  116.  TABB.'COUR'= TABCOUR ;
  117. *                                                 
  118. * ------POTENTIAL  COMPUTATION OR FIRST STEP IF NON LINEAR PROBLEM ---
  119. *------                     
  120.  POT_VECT  TABB 'SOLIN' ;
  121. *-----
  122. *------------------------------------------------------------------
  123.  SOL1 = (TABB.'POTENTIEL' ) ENLEVER LX ;
  124.  RAY = FLAN1 COTE 1 ;
  125.  RAY2 = D 10 (0. .0001) ( .020  .0001) ;
  126. * -------------SOME  POST TRAITMENT   B COMPUTATION ----------------
  127.  BB = INDUCTIO  AIRIN  SOL1  VRAI ;
  128.  PREF = AIRIN POINT PROCHE ( 0. 0. ) ;
  129.  BY0 = EXTR BB 'BY'  PREF  ;
  130.  BY10 = EXTR BB 'BY' ( AIRIN POINT PROCHE ( .010 0.)) ;
  131.  TITRE ' COMPOSANTE BY AVANT ET APRES LISSAGE ' ;
  132.  EVB1  = EVOL ROUG  CHPO  BB    'BY'  RAY  ;
  133. *----------    polynomial  smothing  if  wanted  --------------------
  134.  CHLIS = PROI  POLY  TOUT RAY2   SOL1 1 AXIS ;
  135.  BBY = (EXCO CHLIS 'BY' ) NOMC  'BY' ;
  136.  BY01 = EXTR BBY 'BY' (RAY2 POINT INITIAL) ;
  137.  BAT1  = 800.e6  * 5.e-3 * 5.e-3  * MU0 /10.e-3 ;
  138.  
  139.  option echo 0 ;
  140.  MESS '**************************************************************';
  141.  MESS '*  CIRCULAR COIL  internal radius 20 mm '                      ;
  142.  MESS '*        total    cross section  5*5 mm '                      ;
  143.  MESS '*         SYMMETRY BY HORIZONTAL PLANE  '                      ; 
  144.  MESS '*  AMPERE  mufer >> muo  all AMPERE*TURNS in the GAP of 10 mm '; 
  145.  MESS '*                                                             ';
  146.  MESS '*       J * EP * HAUT = B/ MU0 * e                            ';
  147.  MESS '*   WAITED :   By = '  bat1 '  TESLAS                           '; 
  148.  MESS ' *************************************************************';
  149.  MESS '  NUMERICAL DERIVATION    BY AU CENTRE ' BY0   ;
  150.  MESS '  SMOTHED SOLUTION        BY AU CENTRE ' BY01   ;
  151.  MESS ' *****************************************';
  152. *******************************************************************     
  153.    RAP = ABS ((BY01 - bat1 ) / bat1  );
  154.        SI ( RAP  > .01 ) ; ERREUR 5 ; FINSI ;                            
  155. *******************************************************************     
  156.  SI (NEG  GRAPH  N ) ;
  157.  TITRE ' POTENTIAL   BEFORE AND AFTER   SMOTHING ' ;
  158.  EVPO1 = EVOL ROUG  CHPO  SOL1   'T'   RAY  ;
  159.  EVPO2 = EVOL VERT  CHPO  CHLIS    'A' RAY2  ; 
  160.  dess (evpo1 et evpo2) xbor 0. (R1 * 0.001 * .75); 
  161.  TITRE ' BY COMPONENT  BEFORE AND AFTER SMOTHING  ' ;
  162.  EVB1  = EVOL ROUG  CHPO  BB    'BY'   RAY  ;
  163.  EVB2  = EVOL VERT  CHPO  CHLIS  'BY'  RAY2  ;
  164.  
  165.  dess (evb1 et evb2) xbor 0. (R1 * 0.001 * .75) ;
  166.   FINSI ;
  167. *
  168. **-----------------     NON LINEAR  COMPUTATION  -------------------
  169. *       filling  
  170.  TABB.SOUSTYPE='THERMIQUE' ;
  171.  TABB.CRITERE =1.E-4 ;
  172.  TABB.NITER   =1;
  173.  TABB.'OME' = .1 ;
  174. *-----------  test shortened to 2 iterations ---------------------- 
  175.  TABB.ITERMAX=2;
  176.  TABB.NIVEAU =1;
  177. *
  178. *-------------------------------
  179.   MAG_NLIN  TABB  ;
  180. *-------------------------------
  181. *
  182.  SOL2 = enlever (TABB.'POTENTIEL' ) LX ;
  183.  BB = INDUCTIO  AIRIN SOL2  VRAI ;
  184.  BY02 = EXTR BB 'BY'  PREF ;
  185.  CHLIS = PROI  POLY  TOUT RAY2   SOL2 1 AXIS ;
  186.  BBY = (EXCO CHLIS 'BY' ) NOMC  'BY' ;
  187.  BY03 = EXTR BBY 'BY' (RAY2 POINT INITIAL) ;
  188.  BAT = 2.4106 ;
  189.  MESS ' *****************************************';
  190.  MESS '  EXPECTED SOLUTIONS   AT CENTER   BY  ' BAT      ;
  191.  MESS '  NUMERICAL DERIVATION             BY  ' BY02     ;
  192.  MESS '  SMOTHED SOLUTION                 BY  ' BY03     ;
  193.  MESS ' *****************************************';
  194.     itest = 1 ;
  195. *--------------------GOOD WORKING MESSAGE -----------------------                  
  196.     SI (EGA ITEST 1 ) ;
  197. *------- test   shortened ---------------------------------------
  198.    RAP = ABS ((BY03 - BAT ) / BAT);
  199.        SI ( RAP  > .01 ) ; ERREUR 5 ; FINSI ;
  200. *----------------------------------------------------------------
  201.    SINON ;
  202. *------  COMPUTATION TILL CONVERGENCE OR UP TO 100 MORE ITERATION
  203.  TABB.ITERMAX=100;
  204.   MAG_NLIN  TABB  ;
  205. *----------------------------------------------------------------
  206.  SOL2 = enlever (TABB.'POTENTIEL' ) LX ;
  207.  BB = INDUCTIO  AIRIN SOL2  VRAI ;
  208.  BY02 = EXTR BB 'BY'  PREF ;
  209.  CHLIS = PROI  POLY  TOUT RAY2   SOL2 1 AXIS ;
  210.  BBY = (EXCO CHLIS 'BY' ) NOMC  'BY' ;
  211.  BY03 = EXTR BBY 'BY' (RAY2 POINT INITIAL) ;
  212.  BAT = 2.3468 ; 
  213.  MESS ' *****************************************';
  214.  MESS '  EXPECTED SOLUTION ON AXIS      BY   ' BAT      ;
  215.  MESS '  NUMERICAL DERIVATION  BY AU CENTRE ' BY02     ;
  216.  MESS '  SMOTHED SOLUTION      BY AU CENTRE ' BY03     ;
  217.  MESS ' *****************************************';
  218. * 
  219. *******************************************************************
  220.    RAP = ABS ((BY03 - BAT ) / BAT );
  221.        SI ( RAP  > .01 ) ; ERREUR 5 ; FINSI ;
  222. *******************************************************************
  223.  
  224.   SI (NEG GRAPH N )  ;
  225.  TITRE ' POTENTIAL   BEFORE AND AFTER   SMOTHING ' ;
  226.  EVPO1 = EVOL ROUG  CHPO  SOL2   'T'   RAY  ;
  227.  EVPO2 = EVOL VERT  CHPO  CHLIS    'A' RAY2  ;
  228.  dess (evpo1 et evpo2) xbor 0. (R1 * 0.001 * .75);
  229.  TITRE ' BY COMPONENT  BEFORE AND AFTER SMOTHING ' ;
  230.  EVB1  = EVOL ROUG  CHPO  BB    'BY'  RAY  ;
  231.  EVB2  = EVOL VERT  CHPO  CHLIS  'BY'  RAY2  ;
  232.  dess (evb1 et evb2) xbor 0. (R1 * 0.001 * .75);
  233.  FINSI ;
  234.    FINSI ;
  235.  FIN ;
  236.               ;
  237.  
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  
  244.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales