*23456789012345678901234567890123456789012345678901234567890123456789012
************************************************************************
*                                                                       
*  Input Format of EDDY3D-MJ     EDDY3D-MJ r6.1  Oct 14, 1997
*       Sience Solutions International Laboratory Inc.        
*             2-21-7 Naka-cho, Meguro-ku, Tokyo Japan 153          
*             TEL: +81-3-3711-8908, FAX: +81-3-3711-8910     
*             EMail: support@ssil.com               
*   
*
*     File name : input 
*                 This file must exist in the directry where
*                 calculation is executed.
*
*     File format : Basically, free format devided by blank.
*                   0 or 0.0 cannot be replace by blank.
*                   The raws with * at the 1st column are comments.
*                   Limit of raw length is 100.
*
*     List output file is output.
*     Error message is output to stderr.
*     
*     Finite Element  : Linear Hexahedral
*                       Quadratic Hexahedral
*                       Linear Prism
*                       Linear Tetra
*                       Linear Quadrilateral  (2D case)
*                       Linear Triangle       (2D case)
*                                                                       
************************************************************************
** 0. Release No. ******************************************************
* RLS_NO *
r8.3
************************************************************************
** 1. Job task ********(0;OFF, 1;ON) ***********************************
* PRE_PROCESSING * MAKE_SYSTEM_MATRICES *
       1                  1                    
* SOLVE_EQUATION * POST_PROCESSING *
       1                  1
*-----------------------------------------------------------------------
*
*  PRE_CONDITIONING : Input mesh data and process the data.
*                     Generate edge elements.
*                     Process boundary conditions.
*                     Calculate element volume and element matrices.
*
*                     Input file  : input
*                                   pre_geom ( ATLAS format mesh data )
*                                   Nastran_input
*                                            ( NASTRAN format mesh data)
*                         The mesh data must include nodes and elements.
*
*                     Output file : mesh_data
*                                   element_c_matrices
*                                   element_h_matrices
*                                   post_geom (optional)
*
*  MAKE_SYSTEM_MATRICES : Allocate unknown and fixed variables.
*                         Renumber variables.
*                         Assemble system matrices.               
*                      Output file : potentials
*                                    system_matrices
*                                    source_matrix
*
*  PRE_CONDITIONING  : Incomplete choleskey decomposition in linear
*                      case.
*
*                      Output file : system_r_matrix ( transient case )
*                                    system_t_matrix  
*                
*  SOLVE_EQUATION    : Calculate source term.
*                      Solve eqation.
*
*                      Input file  : old_solutions   (restart)
*                                     copied from file solution of
*                                     previous run.
*                      Output file : solutions
*
*  POST_PROCESSING   : Output results.
*
*                      Output file : current         (optional)
*                                    magnetic        (optional)
*                                    magnetic_iron   (optional)
*                                    force_J_B       (optional)
*                                    force_to_nodes  (optional)
*                                              
*  Notice: The output files of previous processes are input to the
*          successive process.
*          The output files are overwritten by the following run.
*          Rename the files if you want.
*
*-----------------------------------------------------------------------
** 2. Calculation *****(0;OFF, 1;ON) ***********************************
* STATIC * STEP * AC * TRANSIENT * MOTION * NON_LINEAR *
    1       0     0       0          0           1
*-----------------------------------------------------------------------
*
*  STATIC : Static magnetic field calculation.
*  STEP   : Static magnetic field calculation just after step exitation
*           of the external field.
*  AC     : AC ( steady ) magnetic and eddy current field calculation.
*  TRANSIENT : Transient magnetic and eddy current field calculation.
*  MOTION : Transient calculation when the external field moves.
*  NON_LINEAR : Static or transient calculation with non_linear magnetic
*               property.
*
*-----------------------------------------------------------------------
*** 3. Gauge condition *************************************************
* POTENTIAL *
    0
* TREE_GAUGE *
   -1 
*-----------------------------------------------------------------------
*                  Total region        Reduced region
* POTENTIAL  = 0 :    A                  A
*            = 1 :    A                  Omega
*            = 2 :    A-Phi              A
*            = 3 :    A-Phi              Omega
* TREE_GAUGE = -1 : No gauging by tree and co-tree. ** recommended **
*            =  0 : Complete gauging by tree and co-tree. Tree is
*                   constructed randomly.
*            =  1 : The trees in x-directions set to be zero.
*            =  2 : The trees in y-directions set to be zero.
*            =  3 : The trees in z-directions set to be zero.
*
*-----------------------------------------------------------------------
************************************************************************
** 4. Order of shape functions *****************************************
* ORDER_OF_SHEPE_FUNCTIONS *
         1
*-----------------------------------------------------------------------
*
*   =1: Linear hexahedral and prism elements. ** recommended **
*   =2: Quadratic hexahedral elements.
*
*-----------------------------------------------------------------------
** 5. ICCG and non-linear iterations ***********************************
* ICCG_CONV * ACCEL_FACTOR * DIV_FACTOR * DIV_ITERATIONS * SOLVER *
    1.e-6        1.02           2.0           10         0
* NON_LINERAR_CONV * RELAXATION * NO_ITERATIONS *
    1.e-3             1.0            50
*-----------------------------------------------------------------------
*
*  ICCG_CONV         : Convergence error condition of ICCG.
*  ACCEL_FACTOR      : Acceleration facter of Incomplete Cholesky Decomposition.
*                      When gauge is fixed, 1.0 is O.K.
*                      When gauge is not fixed, 1.02~1.2
*                      When =0.0, automatically determined.         
*  DIV_FACTOR        : ICCG is decided to diverge if the error continues 
*  DIV_ITERATINS     :  greater than minimum_error*DIV_FACTOR 
*                       DIV_ITERATION times.
*                      IF convergence is very smooth,
*                       DIV_FACTOR=2., DIV_ITERATIONS=10 are recommended.
*                      IF convergence is spiky or waving,
*                       DIV_FACTOR>10., DIV_ITERATIONS>10.
*  SOLVER           :=0; ICCG Method
*                    =1; Bi_CGSTAB Method 
* 
*  NON_LINEAR_CONV   : Non-linear iterations error condition.
*  RELAXATION        : Relaxation factor of non_linaer iteration.
*                      0.5 - 1 is recomended. If non-linear iteration
*                      diverges, input less value.
*  NO_ITERATIONS     : Max numbers of non-linear iterations.
*  
*-----------------------------------------------------------------------
** 6. GAUSS integration points *****************************************
* POINTS_IN_1D * POINT_IN_TRIANGLE  * MAT_POINTS *
         3               7                 1
*-----------------------------------------------------------------------
*
*  POINTS_IN_1D : Number of gauss integration in one derection.
*                 allowable values -- 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12
*                 recommented value - 3
*  POINT_IN_TRIANGLE : Number of integration points in triangle.
*                 allowable values -- 1, -3, 3, 4, 6, 7, 12
*                 recommended value -- 7 
*
*-----------------------------------------------------------------------
************************************************************************
** 7. Initial value  ***************************************************
* INITIAL_STEP *
     -1
*-----------------------------------------------------------------------
*
* INITIAL_STEP = -1: Zero initial value.
*              > 0 : Step number in file old_solution, which is copyed
*                    from file solution by previous run. 
*
* * The file solution starts from 1 step in static calculation. No
*    initial value is output.
* * The file solution starts from 0 step in transient calculation. In
*    the case, intial value is output to 0 step.
*   
*-----------------------------------------------------------------------
** 8. Calculation steps ************************************************
* NO_STEPS * INITIAL_TIME(s) *  DELTA_TIME(s) *
      4         0.0              1.0 
*-----------------------------------------------------------------------
* 
* NO_STEPS     : Number of calculation steps in static and transient
*                calculation.
*                The 0 step is not included.
* INITIAL_TIME : Initial time for the 0 step. 
* DELTA_TIME   : Time step.
*
*  The last step is recommended more one step than you wants in
*  transient calculation.
*
*  Calculation start from step 1 at static calculation.  
*  The data has no meaning in AC calculation, but the data is necessary.
*
*-----------------------------------------------------------------------
** 9. Output steps *****************************************************
* In transient and static calculation. ***
* INITIAL_STEP * LAST_STEP * STEP_INTERVAL * TIME_DIFFERENCE *  
       1              4          1                3
* In AC calculation ***
* INITIAL_PHASE * LAST_PHSE * PHASE_INTERVAL * 
*       0.0            90.        90.
*-----------------------------------------------------------------------
* 
*   INITIAL_STEP : Initial step of output in transient calculation.
*                  Must be greater than 0.
*   LAST_STEP    : Last step.
*   STEP_INTERVAL: Interval of time step.
*                  Must be greater than 0.
*   TIME_DIFFERENCE :No_points for time defferentiation.
*                    =2 or 3
*                    When =2, output is at mid_points of time steps. 
*
* IN AC calculation.
*   INITIAL_PHASE  : Initial phase (degree) of output.
*   LAST_PHASE     : Last phase (degree). INITIAL_PHASE<LAST_PHASE.
*   PHASE_INTERVAL : Interval of phase step (degree), must be greater
*                    than 0.
*  
* The steps effects on the following print and file output.
*
*-----------------------------------------------------------------------
** 10. Pre-post-files output *******************************************
* INPUT_MESH_FILE * ( 0;No file * 1;IDEAS *2;ATLAS * 3;NASTRAN * 4;KSWAD)
        2
* POST_DATA_FILE ** ELEM_OUT * NODE_OUT *NUMBER_OUTPUT_MATS *
        5               1          1           2
* MAT_ID_NOS *
       11    12     
* MESH * CURRENT * MAGNETIC * FORCE_J_B* FORCE_NODAL * DISP * ELEM *(0;OFF, 1;ON) ***
   1        0         1            0          1          0     0
*-----------------------------------------------------------------------
*
* INPUT_MESH_FILE =0 : Automatic mesh generation of a block. No input
*                      mesh file is neaded.
*                 =1 : I-DEAS format input mesh file IDEAS_mesh.
*                 =2 : ATLAS format input mesh file pre_geom.
*                 =3 : NASTRAN format input mesh file NASTRAN_mesh.
*                 =4 : KSWAD format input mesh file KSWAD_mesh.
*
* POST_DATA_FILE  =0 : No file output.
*                 =1 : Output IDEAS Univesal file-format files.
*                 =2 : Output ATLAS file-format files.
*                 =4 : Output KSWAD file-format files.
*
* MESH    : Output mesh data into file post_geom.
*
* CURRENT : Output current and Joule loss distributions into file
*            current.
*
* MAGNETIC: Output magnetic field distributions in file magnetic.
*           =2 --> The magnetic field in non_magnetic and magnetic
*                  materials is output separately into files magnetic
*                  and magnetic_iron, respectively.
*           >0 --> Induced + applied field. 
*           <0 --> Only induced magnetic field.
*
* FORCE_J_B : Output force distribution calculated by Lorents force when
*             no magnetic material is includes. Into file force_J_B.
*             Unit N/M**3. Only in ATLAS format.
*
* FORCE_NODAL : Output force distribution calculated by nodal force
*               formulation into file force_to_nodes.
*               Unit N. Only in ATLAS format.
*
*-----------------------------------------------------------------------
** 11. Print output ****************************************************
** 11.1 Print output_options *(0;OFF, 1;ON) ****************************
*   MESH   *   A   *   V   *   B   *   B_INTEG *    J_Q  * 
     0         0       0       0         0           0
*  Q_AVERAGE *  Q_TOTAL  *  MAG_FLUX  * CUR_FLUX *
        0         0            1            0
* FORCE_J_B * FORCE_NODAL *
        0         1
*-----------------------------------------------------------------------
*
* MESH     : Mesh data.
* A        : Magnetic vector potentail on edges. Integrated value along
*            edges. 
* V        : Electric scalar potential at nodes. Time Integrated values. 
* B        : Magnetic flux density at element centers.
* B_INTEG  : Magnetic flux density at space points by integration of
*            magnetization. Only in static calculation. The data 11.4 is
*            needed.
* J_Q      : Current density and Joule heat loss at element centers.
* Q_AVERAGE: Joule heat loss averaged in each elements.
* Q_TOTAL  : Total Joule heat loss in all calvculated region.
* MAG_FLUX : Magnetic fluxes through closed lines. The data 11.2 is
*            needed.
* CUR_FLUX : Total currents through surfaces. The data 11.3 is needed.
* FORCE_J_B : Force densities at element centers from Lorents force.
* FORCE_NODAL : Forces in material regions from nodal force formulation.
*
*-----------------------------------------------------------------------
** 11.2. Magnetic flux loops *******************************************
* NO_LOOPS *
      2
* NO_NODES *
      15
* NODE_IDS *
  22733 22771 22809 22847 22885 22923 22961
  22960 22959 22958 22957 22956 22955 22954 22953
      33
  22754 22792 22830 22868 22906 22944 22982 23020 23058 23096 23134 23172
  23210 23209 23208 23207 23206 23205 23204 23203 23202
  23164 23126 23088 23050 23012 22974 22936 22898 22860 22822 22784 22746
*-----------------------------------------------------------------------
* The data is necessary when MAG_FLUX Print out is switched on.
*
*  NO_LOOPS : Number of following flux loops.
*  NO_NODES : Number of nodes connecting the loop.
*  NODE_IDS : Node ID nubmers.
*
*-----------------------------------------------------------------------
** 11.3. Current flux surfaces *****************************************
*number of surfaces *
* NO_SURFACES *
*      2
*  MAT_ID_NOS  *
*    11      12
*-----------------------------------------------------------------------
* The data is necessary when CUR_FLUX Print out is switched on.
*
*  NO_SURFACES : Number surfaces on which normal currents are
*               integrated.
*
*  MAT_ID_NOS  : Material ID numbers of the surface elements.
*                The surface elements data must be included in file
*                pre_geom or NASTRAN_input. 
*
** The surface elements must be coincide with the surfaces of volume
*    elements.
*-----------------------------------------------------------------------
** 11.4. B calculation by integrating nagnetization ********************
* XSYMMETRY * YSYMMETRY * ZSYMMETRY * ROTATION *
*   1             1           -1           1
* NO_POINTS *  NO_MAT_IDS * NO_COIL_SERIES *
*   2             1            1
*   XP(m) *    YP(m)  *      ZP(m)  *
*   0.2         0.           0.0
*   0.3         0.           0.0
* MAT_IDS *
*    1
* COIL_SERIES 
*    1
*-----------------------------------------------------------------------
* The data is necessary when B_INTEG Print out is switched on.
*
*  XSYMMETRY =0: No symmetry to x=0 plane of magnetization.
*            =1: Symmetry to x=0 plane of magnetization.
*           =-1: Anti-symmetry to x=0 plane of magnetization.
*  YSYMMETRY :
*  ZSYMMETRY : Just same as XSYMMETRY.
*  ROTATION  : Number of rotational symmetry.
*             =1 means no rotational symmetry.
*  NO_POINTS : Number of calculation points.
*  NO_MAT_IDS: Number of calculated materials.
*  NO_COIL_SEIRES: Number of calculated_external coils.
*              NO_MAT_IDS + NO_COIL_SERIES >=1.
*  XP, YP, ZP (m) : Coordinates of calculation points.
*  MAT_IDS   : Material IDs in 16 when NO_MAT_IDS !=0.
* COIL_SERIES : Coil series Numbers in 17.1 when NO_COIL_SERIES !=0.
*----------------------------------------------------------------------- 
************************************************************************
** 12. Geometry ********************************************************
* GEOMETRY * DELTA_Z(m) *
     0         1.0
* GEOMETRY * DETLA_THETA(deg) *
*    2         1.0
*-----------------------------------------------------------------------
*
* GEOMETRY = 0 : 3-Dimensional calculation.
*          = 1 : 2-Dimensional, Uniform in z direction.
*          = 2 : Axi-symmetric, Uniform in theta direction.
*          = 3 : Make 3D Axi-symmetric model from 2D data.
*          = 4 : 2D-data (x-y base plane) is expanded to z direction.
*                File mesh_gen is needed. 
*
* The 2-dimentional elements are expanded to 3-dimentional elements.
*
* DELTA_Z (m) : Expansion length to z_direction.
* DELTA_THETA (degree) : Expansion angle to theta_direction.
*
***In 2-dimentional, elements must be defined 
*       at x-y (z=0 ) plane in transitional symmetry (GEOMETRY =1)
*    or at x-z (y=0) in axisymmetric case  ( GEOMETRY =2,3)
*  by linear quadriraleral or triangle.
*  Quadratic elements are not allowed.
***In axisymmteric, the elements connecting to center axis must be
*  quadriaterals.
***The surface number 1 shown in Appendix 1. is defined. Volume elements
*  are generated automatically.
*  
*-----------------------------------------------------------------------
** 12.2 Local coordinate ***********************************************
*COORDINATE * NO_COORDINATE *
COORDINATE         0
*   ID      * TYPE          *    X0    *    Y0     *   Z0    *
*                           *    EX_X  *    EX_Y   *   EX_Z  *
*                           *    EZ_X  *    EZ_Y   *   EZ_Z  *  
*   1           1               0.          0.         0.
*                               1.          0.         0.
*                               0.          0.         1.
*-----------------------------------------------------------------------
*
*    ID           : Coordinate ID No. >0
*    TYPE         :=1 : Cartesian coordinate
*                  =2 : Cylindrical coordinate     
*    X0, Y0, Z0   : Origin global coordinate
*    EX_X, _Y, _Z : Unit local x_directional vector
*    EZ_X, _Y, _Z : Unit local z_directional vector
*-----------------------------------------------------------------------
** 13. Boundary conditions *********************************************
* FAR_BOUNDARY_CONDITION * IMPEDANCE_BNC. * 
* ( 0;At=0(Bn=0), 1;At=free(Bt=0), 2;BE boundary )*
        1                     0               
* NO_BE_TERMS * BE_CENTER X(m), Y(m), Z(m) *
*       3                  0.0   0.0   0.0 
* DISTANCE_JUDGE(m) * Distance for judgement (m)
      1.e-4
** At=0 (Dirichlet, Bn=0 ) planes
* NO_DIRICHLET_PLANE *
    2
* coefficients ( CX*X + CY*Y + CZ*Z = C )
*     CX        *       CY        *      CZ       *      C(m)   *
      1.0               0.0              0.0             0.0
      0.0               1.0              0.0             0.0
** At=free (Neumann, Bt=0 ) planes
* NO_NEUMANN_PLANE Number
    0
* coefficients ( Cx*X + Cy*Y + Cz*Z = C )
*     Cx         *       Cy        *      Cz       *      C(m)   *
*      0.0               0.0              1.0             0.0
** At=0 lines 
* NO_A_0_LINE
    0
*          start point              *     end point               *
*    SX(m) *    SY(m)  *    SZ(m)  *    EX(m)  *    EY(m)  *    EZ(m)  *
*    0.0        0.0        0.0         0.0         0.0         1.0
*-----------------------------------------------------------------------
*
* FAR_BOUNDARY_CONDITION :
*     =0 : At = 0 (Bn = 0 ) at far boundary.
*     =1 : At = free (Bt = 0) at far boundary.
*     =2 : Boundary elements
*When FAR_BOUNDARY_CONDITION==2,
*  NO_BE_TERMS       : Nomber of expanded terms of bounary elements
*  BE_CENTER         : Center of boundary elements    
*
* IMPEDANCE_BNC    =1: Use impedance boundary condition. 
* DISTANCE_JUDGE (m) : The distance less than this value is judges to be
*     0. It is used in judgement if a node is on a plane or not, etc.
*
* NO_DIRICHLET_PLANE : Number of the planes of Dirichlet boundary
*                      condition ( At = 0 (Bn = 0 ) ).
* NO_NEUMANN_PLANE   : Number of the planes of Neumann boundary
*                      condition ( At = free (Bt = 0)).
* NO_A_O_LINE       : Number of the lines on which At = 0.
*
* CX, CY, CZ, C     : The planes are represented by the equation.
*                       CX * x + CY * y + CZ * z = C (Unit m).
* SX, SY, SZ        : Coordinate of one point on a line.
* EX, EY, EZ        : Coordinate of another point on the line.
*  Unit (m)
*
*-----------------------------------------------------------------------
** 14. Cyclic symmetricity *********************************************
* CYCLIC     *  SYMMETRICITY *  DIVISION *  LINE_ID *
     0               1             0          11
*   DX(m)*   DY(m)*   DZ(m)*
*   0.0      0.0      1.0
*   ANGLE(deg) *
*   180.0
*-----------------------------------------------------------------------
*
*  CYCLIC  =0 : No cyclic symmetry.
*          =1 : Transitional symmetry.
*                 A(x+DX,y+DY,z+DZ)= +(-)A(x,y,z)
*          =2 : Rotational Symmetry.
*                 A(r,z,t+ANGLE)= +(-)A(r,z,t)
*          =3 : Claw symmetry.
*                 Ar(r,-z,t+ANGLE)= +(-)Ar(r,z,t)
*                 Az(r,-z,t+ANGLE)= -(+)Az(r,z,t)
*                 At(r,-z,t+ANGLE)= +(-)At(r,z,t)
*
*  SYMMETRY =0 : Symmetric.
*           =1 : Anti_symmetric, sign in ().
*
*In case MOTION (in 2 ) == 2,
*  DIVISION    : Number of division in displacement direction.
*  LINE_ID     : MAT_ID Number of line or point elements 
*                defining sliding face.
*
* When CYCLIC=1
*  DX, DY, DZ : Periodic distance ( m )
*                 A(x+DX,y+DY,z+DZ)= A(x,y,z)
*
* When CYCLIC=2, 3
*  ANGLE      : Periodic angle (degree)
*
*-----------------------------------------------------------------------
************************************************************************
** 15. Block dat for automatic mesh generation *************************
*
*      Following data is necessary when INPUT_MESH_FILE = 0
** 15.1. Finite elements ***********************************************
* ELEMENT_TYPE *
*     0
*-----------------------------------------------------------------------
*
* ELEMENT =0 : Hexahedral.
*         =1 : Prism.
*         =2 : Tetra
*
*-----------------------------------------------------------------------
** 15.2 Coordinates of geometry ****************************************
* X_AIR_LOWER  * X_MAT_LOWER * X_MAT_UPPER * X_AIR_UPPER *
*      0.0           0.0           0.1           0.22
* Y_AIR_LOWER  * Y_MAT_LOWER * Y_MAT_UPPER * Y_AIR_UPPER *
*      0.0           0.0           0.1           0.22
* Z_AIR_LOWER  * Z_MAT_LOWER * Z_MAT_UPPER * YZAIR_UPPER *
*      0.0            0.0          0.005         0.22
*-----------------------------------------------------------------------
*
* X_AIR_LOWER : x coordinate at lower limit of air region.
* X_MAT_LOWER : x coordinate at lower limit of material region.
* X_MAT_UPPER : x coordinate at upper limit of material region.
* X_AIR_UPPER : x coordinate at upper limit of air region.
*    Unit m.
*    X_AIR_LOWER <= X_MAT_LOWER < X_MAT_UPPER < X_AIR_UPPER.
** The same in the y and z directions.
*
** Material number of material region is 1.
** Material number of air region is 2.
*
*-----------------------------------------------------------------------
** 15.3. Numbers of divisions ******************************************
**  Lower air  *  Conductor        * Upper air *
*  XDIV_AIR_LOWER * XDIV_MAT  * XDIV_AIR_UPPER *
*         0            -6               -6
*  YDIV_AIR_LOWER * YDIV_MAT  * YDIV_AIR_UPPER *
*         0            -6               -6
*  ZDIV_AIR_LOWER * ZDIV_MAT  * ZDIV_AIR_UPPER *
*         0             1               -6
*-----------------------------------------------------------------------
*
* XDIV_AIR_LOWER : Number of division from X_AIR_LOWER to_X_MAT_LOWER
* XDIV_MAT_LOWER : Number of division from X_MAT_LOWER to_X_MAT_UPPER
* XDIV_AIR_UPPER : Number of division from X_MAT_UPPER to_X_AIR_UPPER
*
** The same in the y and z dirextions.
** When the division is positive, mesh is equi-distance.
** When the division is negative, the division ration must be input.
*
*-----------------------------------------------------------------------
** 15.4. Divisin ratios ************************************************
* DIVISION_RATIO *
*      0.02      0.02      0.02      0.02      0.01      0.01
*      0.01      0.01      0.02      0.02      0.03      0.03
*      0.02      0.02      0.02      0.02      0.01      0.01
*      0.01      0.01      0.02      0.02      0.03      0.03
*     0.005      0.01      0.02      0.04      0.06      0.08
*----------------------------------------------------------------------
*
* Division ratios to the negative values correspomding to previous data.
*
*-----------------------------------------------------------------------
************************************************************************
************************************************************************
** 16. Material properties and potential regions ***********************
* NO_MAT_IDS *  EXTEND_TOTAL * NO_SMAT_IDS *
     3               1            0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID * SIGMA *  MU     * PACKING *
     1            1          0          0.0      1.0
    11            0          1          0.0      1.0       1.0
    12            0          1          0.0      1.0       1.0
* COORDINATE_ID  *  FX(FR)   *  FY(FT)  *  FZ  *       
*     0              1.0         0.5       0.0 
* COORDINATE_ID  *  DX(DR)   *  DY(DT)  *  DZ  *       
*     0              0.0         0.0      1.0 
*      2           1             0         0.0              1.0
* SMAT_ID *  SIGMA   *  MU   *
*     1      5.e7      1.0 
*-----------------------------------------------------------------------
*
* NO_MAT_IDS   : Number of following MAT_ID's.
* NO_SMAT_IDS  : Number of following SMAT_ID's.
* EXTEND_TOTAL : When =1, total potential region is extended to reduced 
*                potential region by a layer. This is needed when 
*                calculate the nodal force.
* MAT_ID       : Matrial ID number. The ID numbers included in mesh data
*                must be input.
* SMAT_ID      : Surface matrial ID number for inpedance surface. 
* POTENTIAL   =0: Elements in total potential region.
*             =1: Elements in reduced potential region.
* B_H_CURVE_ID : When not equal to 0, the material is non_linear and the
*                corresponding B-H curve must be input in 20.
* SIGMA        : Electric conductivity of the material (S/m).
* MU           : Relative magnetic permeability of linear material.
* PACKING_FACTOR   : For non_linear layered iron, the ratio of the iron 
*                    in volume. 
*
** Negative value of SIGMA or MU means anisotropy.
*  Following data are needed.        
* COORDINATE_ID : ID No. of local coordinate defining anisotropy. (12.2)
* FX, FY, FZ    : Conductivity or permeability is multiplied by the factors
* (FR, FT, FZ)           in each direction of local coordinate. 
*
** When PACKING_FACTOR is not 1.0.    
* COORDINATE_ID : ID No. of local coordinate defining packing. (12.2)
* DX, DY, DZ    : Packing direction in the coordinate.
* (DR, DT, DZ)    Isotropic in the normal direction.       
*
** The conductor and magnetic material must be in total potential
*  region.
** The coils defined by the data in 17.1 must be in reduced potential
*  region or outside of calculation region.
*
*-----------------------------------------------------------------------
** 17. Source fields ***************************************************
*SOURCE *
SOURCE  
*-----------------------------------------------------------------------
* Source field data start by SOURCE card and stop by END card.
* Series of field sources are defined by
*             COIL(s),
*             ELMCUR,
*             SDEFCOIL,
*             SUFCUR,
*             MAGNET,
*             CIRCUIT
* No rule to input order of the series.
************************************************************************
** 17.1 Source current or source magnetic field given by field elements*
*COIL   * SERIES_NO  *  TIME_ID  *  NO_ELEMENTS * MOTION_ID *
COIL         1            0             8              
*-----------------------------------------------------------------------
* 
*  TITLE       = COIL
*  SERIES_NO   : Source series no. Multiple coil series are O.K.
*  TIME_ID     : Time data number.
*  NO__ELEMENTS : Number following field elements.
*  MOTION_ID   : ID nunmer of motion in in 19.
*                Necessary when MOTION in 2 is equal 1.
*                Blank or 0 means no motion.
** Following field elements are input with title and the field is added.
*  The end of data of one coil series is shown by END.
*  The fields elements are UNIF, LINE, LOOP, GCE, ARC and DIPO. 
*-----------------------------------------------------------------------
** 17.1.0 Coordinate defining field geometry ***************************
*COORDID * COORDINATE_ID * 
*COORDID    1
*-----------------------------------------------------------------------
* COODINATE_ID : Coordinate ID number (12.2).
*                The fields after thid card are defined 
*                   in the local coordinate.   
*                If no COORDIN card, the global coordinate is used.
** 17.1.1 Uniform field ************************************************
*UNIF  *   BX(T)    *    BY(T)   *   BZ(T)    * A_FORM *
*UNIF        0.           0.           1.         0
*-----------------------------------------------------------------------
*
* BX, BY, BZ (T): Uniform magnetic flux density.
*
* A_FORM : Form of magnetic vector potential. Select the form which
*          satisfies boundary conditions.
*      =0: A = (B x r ) /2, where r is postion vector.
*      =1: Ax = By*z, Ay = Bz*x, Az = Bx*y.
*      =2: Ax =-Bz*y, Ay =-Bx*z, Az =-By*x.
*
*-----------------------------------------------------------------------
** 17.1.2 Straight line current with infinite length *******************
*LINE    * CURRENT(A) *
*            X0(m)    *  Y0(m)  *  Z0(m)   *   DX   *   DY   *   DZ    *
*LINE        2.0e7
*            1.0         0.0       0.0        0.0       1.0       0.0
*LINE        2.0e7
*           -1.0         0.0       0.0        0.0      -1.0       0.0
*LINE        2.0e7
*            0.0         1.0       0.0       -1.0       0.0       0.0
*LINE        2.0e7
*            0.0        -1.0       0.0        1.0       1.0       0.0
*-----------------------------------------------------------------------
*
* CURRENT    : 
* X0, Y0, Z0 : Coordinate of a point on the line.
* DX, DY, DZ : Direction vector of the line.
*
*-----------------------------------------------------------------------
** 17.1.3 Loop current with axis on z-axis *****************************
*LOOP  * CURRENT(A) * RADIUS(m)* CENTER_Z(m)* RADIAL_W(m)*  AXIAL_W(m) *
*LOOP    1.591549e6       1.0         0.0       0.1       0.1        
*-----------------------------------------------------------------------
*
* CURRENT :
* RADIUS  : Radius from axis to center circle of the loop.
* RADIAL_W: Radial width.
* AXIAL_W : Axial width.
*
*-----------------------------------------------------------------------
** 17.1.4. General current element *************************************
*GCE    * CURRENT(A)*
*            XS(m)  *   YS(m)  *  ZS(m)  *  XE(m)  *  YE(m)  *  ZE(m)  *
*            W1X(m) *   W1Y(m) *  W1Z(m) *  W2X(m) *  W2Y(m) *  W2Z(m) *
*GCE        2.0e7
*           1.0         -1.0       0.0       1.0       1.0       0.0
*           0.0          0.0       0.05      0.05      0.0       0.0
GCE              1.0
              0.0285   -0.0145     0.075   0.0285      0.0145   0.075
                 0.0       0.0    0.0483   0.009      0.0        0.0
GCE              1.0
               0.0145    0.0285    0.075    -0.0145    0.0285     0.075
                 0.0       0.0    0.0483      0.0    0.009       0.0   
GCE              1.0
             -0.0285     0.0145     0.075   -0.0285    -0.0145     0.075
                 0.0       0.0    0.0483   -0.009       0.0       0.0
GCE              1.0
              -0.0145   -0.0285     0.075     0.0145  -0.0285     0.075
                 0.0       0.0    0.0483      0.0   -0.009       0.0  
*-----------------------------------------------------------------------
*
* CURRENT    :
* XS, YS, ZS : Coordinate of the start point of the center line.
* XE, YE, ZE : Coordinate of the end point of the center line.
* W1X, W1Y, W1Z : Half width vector 1.
* W2X, W2Y, W2Z : Half width vector 2.
*
** The direction of W1 x W2 (vector product) is the same direction of
*   center line.
** The magnitudes of vectots W1 and W2 mean the half length of
*   rectangular cross section.
*
*-----------------------------------------------------------------------
** 17.1.5 Arc element **************************************************
*ARC    *  current(A)  *
*          X(m)  *  Y(m)  *   Z(m)  * RADIUS(m)* AXIAL_W(m)*RADIAL_W(m)*
*          ALPHA(deg)* BETA(deg) * PHI1(deg)   * PHI2(deg) *
*ARC       1.591549e6
*          0.0      0.0       0.0      1.0         0.1        0.1
*          0.0      0.0     -180.      180.
ARC              1.00
               0.0145     0.0145    0.075    0.014      0.0966    0.018
                 0.0       0.0       0.        90.
ARC              1.00
              -0.0145     0.0145    0.075    0.014      0.0966    0.018
                 0.0       0.0      90.       180.
ARC              1.00
              -0.0145    -0.0145    0.075    0.014      0.0966    0.018
                 0.0       0.0    -180.      -90.
ARC              1.00
               0.0145    -0.0145    0.075    0.014      0.0966    0.018
                 0.0       0.0     -90.        0.
*-----------------------------------------------------------------------
*
* CURRENT   :
* X, Y, Z   : Coordinate of the center of arc.
* RADIUS    : Radias from axis to center circle of the arc.
* AXIAL_W   : Axial width.
* RADIAL_W  : Radial width.
* ALPHA     : Eular angle, rotation around z-axis.  First.
* BETA      : Eular angles, rotation around x-axis. Second.
* PHI1      : Start angle from x-axis.
* PHI2      : END angle from x-axis.
*             PHI2 > PHI1.
*   Where axis is fixed on the arc.
*-----------------------------------------------------------------------
** 17.1.6 Dipole field. ************************************************
*DIPO   *  DBZDX(T/m) *
*DIPO       1.0
*-----------------------------------------------------------------------
*
* DBZDX : Gradienet of Bz in direction x.
*         Bx = DBZDX * z, By =0, Bz = DBZDX *x.
*         Ax = 0, Ay = DBZDX* ( x**2-z**2), Az =0.
*-----------------------------------------------------------------------
* END for a COIL series ************************************************
END
************************************************************************
** 17.2 Source current in elements *************************************
*ELMCUR  *SERIES_NO *NO_MAT_IDS *CURRENT_OP  *
*ELMCUR      1          1           1       
* MAT_ID  * IN_SURFACE  * OUT_SURFACE  *  CURRENT * CONDUCTIVITY *
*    1          1              2           1.e6       5.0e7 
*-----------------------------------------------------------------------
*
*  TITLE       = ELMCUR
*  SERIES_NO   : Source series no.
*  NO_MAT_IDS  : Number of following material IDs in magnet region.
*  CURRENT_OP=0: Currents through element surface are input in A.
*            =1: Current densities are input in A/M**2.
*                If =0., resistance ic calculated from CONDUCTIVITY.         
*  MAT_ID      : Material ID number of the current region.
*                Different MAT_IDs are necessary to diffrent following
*                flow conditions to elements.
*  IN_SURFACE  : Element surface number that current flows through in.
*  OUT_SURFACE : Element surface number that current flows through out.
*                Surface numbers are shown Appendix 1.
*  CURRENT     : Total current (A) throgh an element surface when
*                OPTION =0.
*                Current density (A/m**2) when OPTION=1.
* CONDUCTIVITY : Effective conductivity of internal current element.
*                It has meaning when RESISTIVITY is not zero.
*
*** The elements defined in ELMCUR must be in total potential region.
*-----------------------------------------------------------------------
************************************************************************
** 17.3 Suface defined coil        *************************************
*SDEFCOIL  *SERIES_NO *NO_MAT_IDS *
*SDEFCOIL      1          1        
* MAT_ID * S_MAT_ID1 * S_MAT_ID2 * S_MAT_ID3 * S_MAT_ID4 *CURRENT *COND * 
*   1         11         12          13          14         1.e6    5.e7
*------------------------------------------------------------------------
*  TITLE    = SDEFCOIL
*  NO_MAT_IDS  : Number of following material IDs in coil region. 
*  MAT_ID      : Material ID number of the current region.
*  S_MAT_ID1   : Surface Material ID's surrounding the coil region.
*  S_MAT_ID2  
*  S_MAT_ID3
*  S_MAT_ID4
*  CURRENT     : Total current (A) throgh an element surface when
*                OPTION =0.
*  COND : Effective conductivity of internal current element.
*                It has meaning when RESISTIVITY is not zero.
*
*** The elements defined in SDEFCOIL must be in total potential region.
*
*************************************************************************
** 17.4 Source current through conductor sourface **********************
*SUFCUR  *  SERIES_NO * TIME_ID * SURFACE_MAT_ID *  TURN *
*SUFCUR      1         1          1             1.0
*-----------------------------------------------------------------------
*
*  TITLE       = SUFCUR
*  SERIES_NO   : Source series no.
*  TIME_ID     : Time data number.
*  SURFACE_MAT_ID : SURFACE_MAT_ID of surface elements.
*                   Surface element (quadrilateral or triangle) must be 
*                   input in file pre_geom or Nastran_input.
*  TURN        : Current(A) to unit series current. 
** The surface must be on boundary surface of conductor region.
*
*-----------------------------------------------------------------------
************************************************************************
** 17.5 Source magnetization *******************************************
*MAFGNET *  SERIES_NO * TIME_ID *  NO_MAT_IDS *
*MAGNET    1          1        1
* MAT_ID   *  COORD_ID   *   MX(MR)    *    MY(MT)   *     MZ    *
*    1          1                0.         0.            1.0
*-----------------------------------------------------------------------
*  TITLE         = MAGNET
*  SERIES_NO     : Source series number.
*  TIME_ID       : Time data number.
*  NO_MAT_IDS    : Number of following material IDs in magnet region.
*  MAT_ID        : Material ID of magnetized region.
*                  Different MAT_IDs are necessary to diffrent following
*                  magnetization to elements.
*  COORD_ID      : Local coordinate ID No. (12.2)
*                :=0 : Global coordinate.
*  MX, MY, MZ    : Magnetization in local coordinate (Wb/m**2). Cartesian
*  MR, MT, MZ    :                                              Cylindrical
*-----------------------------------------------------------------------
************************************************************************
** 17.6 Circuit connection and power supply ****************************
*CIRCUIT *
CIRCUIT
*  NO_SERIES  *  NO_POWER_SUPPLIES  *  REGION_FACTOR *
        1                1                   4.
*       3                2                   1
* SERIES_IDS *
        1
*       1   2    3
* SERIES_EXTERNAL_INDUCTANCE * 
        0.0
*       0.0
*       0.0 0.0
*       0.0 0.0 0.0  
*  SIRIES_EXTERNAL_RESISTANCE  *
        0.0
*       0.0
*       0.0 0.0
*       0.0 0.0 0.0 
*  CONNECTION_MATRIX *
        1.0 
*       1.0  0.0
*       0.0  1.0
*      -1.0  1.0 
*POWER_SUPPLIES *
* PS_ID    *    TYPE  * TIME_ID *   INITIAL_CURRENT *
     1            0         1                0.0   
*    1            1         1                0.0 
*    2            1         2                0.0          
*-----------------------------------------------------------------------
* 
*  TITLE       = CIRCUIT
*  NO_SERIRS   : Number of series in FEM region.
*  NO_POWER_SUPPLIES  : Number of connected power supplies.
*  REGION_FACTOR : Ratio of real region to model region.
*  SERIES_IDS : Series numbers in 17.1, 17.2, 17.3.
*  SERIES_EXTERNAL_INDUCTANCE (H)
*              : Inductance matrix(lower) connected to the series.
*  SERIES_EXTERNAL_RESISTANCE (OHM) 
*             : RESISTANCE matrix(lower) connected to the series.
*  CONNECTION_MATRIX : Series currents 
*                       =CONNECTION_MATRIX * Power supply currents.
* POWER_SUPPLIES
*   ID   : Power supply ID number.
*   TYPE :=0 Current given.
*         =1 Voltage given.
*   TIME_ID : : Time data number of current or voltage.
*   INITIAL_CURRENT: Initial current in transient calculation.
* END 0f Circuit data **************************************************  
END
* END of field sources *************************************************
END
************************************************************************
** 18. Time data *******************************************************
* NO_DATA
      1
* TIME_ID OPTION
*      1     0
*   C0    *   C1    *   C2    *   C3     *   C4     *  C5    *   C6    *
*    0.0        1.0       0.0       0.0        0.0       0.0       0.0
*  TEXP   *   TCYCLE  *
*    0.0       0.00
* TIME_ID OPTION
     1      1
* NO_TIME_POINTS
     5
*   Time (sec)   *   Value  *
     0.0             0.0
     1.0            1000.
     2.0            3000.
     3.0            4500.
     4.0            5000.
* TIME_ID OPTION
*      1     2
*  AMPLITUDE   *  TCYCLE(s)   *   PHASE(deg)  *
*     1.0           0.001            0.0 
*-----------------------------------------------------------------------
*  
*  NO_DATA         : Number of following time data for field sources.
*
*  TIME_ID         : ID number of the time data.
*  OPTION       = 0: Time variation is given by a equation.
*               = 1: Time variation is given by linear lines.
*               = 2: Time variation is given by cosine wave.
* In case OPTION =0,
*
*  CO, C1(/sec), C2, C3, C4, C5, C6,
*  TXEP(sec), TCYCLE(sec)
*    Time variation is given the equation.
*       F(t)= C0 + C1*t + C2*exp(t/TEXP) 
*           + C3*sin(2pi*t/TCYCLE) + C4*cos(2pi*t/TCYCLE)
*           + C5*exp(t/TEXP)*sin(2pi*t/TCYCLE) 
*           + C6*exp(t/TEXP)*cos(2pi*t/TCYCLE)
*
*      In AC calculation, only C4 and TCYCLE have meaning.
*
* In case OPTION = 1,
*
*   NO_TIME_POINT : Number of following points connecting linear lines.
*
*   TIME (sec) : Times at the points.
*   VALUE      : Magnitude of the applied field at above times.
** At time outside defined time range, the value is linearly 
*    extrapolated.
*
* in case OPTION = 2,
*
*   AMPLITUDE   ;    F(t) = AMPLITUDE*cos(2pi*t/TCYCLE+phase)
*   TCYCLE (sec)
*   PHASE(deg)
*
*     TCYCLE must be equal in all time data in AC calculation. 
*
** The fields is normalized by the data in 17.
*
*-----------------------------------------------------------------------
************************************************************************
** 19. Motion of applied fields ****************************************
* TITLE * NO_MOTION_DATA
*MOTION        1
*  MOTIN_ID  * X_TIME_ID * Y_TIME_ID * Z_TIME_ID * PHI_TIME_ID *
*     1        0         0           0            2
* X_CENTER * Y_CENTER * Z_CENTER *  NX  *  NY  *  NZ  *
*     0.0       0.0        0.0       0.0    0.0    1.0 
* GLOBAL_MOTION_ID *
*      0
*-----------------------------------------------------------------------
* The data is necessary when MOTION( in 2 ) =1 or 2.
* 
*  NO_MOTION_DATA : Numbers of following motion data.
*  MOTIN_ID          : ID number of motion. >0.
*  X_TIME_ID         : Time variatin ID number in 18 for the motion
*                       in x drection of global coordinate(unit m ). 
*  Y_TIME_ID         :  in y drection of global coordinate(unit m ).
*  Z_TIME_ID         :  in z drection of global coordinate(unit m ).
*  PHI_TIME_ID       : Time variatin ID number in 18 for the rotation
*                      (unit degree).
*  X_CENTER, YCENTER, ZCENTER : Global coordinate of rotation center.
*  NX, NY, NZ        : Direction of the rotation axis.
*                      Necessary when PHI_TIME_ID is not 0.
*  GLOBAL_MOTION_ID  : Motion ID number  for the global motion of 
*                      the meshed model (MOTION=1) 
*                      or rotor model (MOTION=2).             
*-----------------------------------------------------------------------
************************************************************************
** 20. B_H curves ******************************************************
* NO_CURVES *
     1
* BH_CURVE_ID * NO_DATA *
     1             39       
*   H_AT/M     *      B_T    *
0.0     0.0
27.     0.01
58.     0.025
100.    0.05
153.    0.10
185.    0.15
205.    0.20
233.    0.30
255.    0.40
285.    0.50
320.    0.60
355.    0.70
405.    0.80
470.    0.90
555.    1.00
673.    1.10
836.    1.20
1065.   1.30
1220.   1.35
1420.   1.40
1720.   1.45
2130.   1.50
2670.   1.55
3480.   1.60
4500.   1.65
5950.   1.70
7650.   1.75
10100.  1.80
13000.  1.85
15900.  1.90
21100.  1.95
26300.  2.00
32900.  2.05
42700.  2.10
61700.  2.15
84300.  2.20
110000. 2.25
135000. 2.30
6262000.  10.0
*-----------------------------------------------------------------------
* The data is necessary when calculation is non_linear.
*
*  NO_CURVES : Number of the following B-H curves.
*
*  BH_CURVE_ID : ID number of B-H curve.
*  NO_DATA     : Number of data points.
*  H_AT/M      : Magnetic intensity H in AT/m
*  B_T         : Magnetic Flux density B in T.
*
** Initial data must be 0.0 0.0.
*-----------------------------------------------------------------------
************************************************************************
** End of data of file input *******************************************
************************************************************************
** Output of time dependent quantities *********************************
* Output file : history
*HISRTORY
* TIME POSITIONX POSITIONY POSITIONZ
* Joule Heating
* QT Qn  
* Forces
* FORCXT FORCYT FOTRCZT
* FORCXn FORCYn FORCZn
* n: Material ID Number. 
*END
