* RLS_NO *
r8.1
* PRE_PROCESSING * MAKE_SYSTEM_MATRICES *
       1                 1          
* SOLVE_EQUATION * POST_PROCESSING *
       1                  1
* STATIC * STEP * AC * TRANSIENT * MOTION * NON_LINEAR * TEMP_DEPEND *
    0       0     0       1           3              1
* POTENTIAL * NODAL_FORCE_OPTION *
    0
* TREE_GAUGE *
   -1 
* NODE_ORDER *
         1
* ICCG_CONV * ACCEL_FACTOR * DIV_FACTOR * DIV_ITERATIONS * SOLVER * THETA * THETA_NETWOEK * THETA_MOTION *
    1.e-6        1.            10.          10                0      0.0        0.5            0.5
* NON_LINERAR_CONV * RELAXATION * NO_ITERATIONS *
    1.e-3             1.0            20
* POINTS_IN_1D * POINT_IN_TRIANGLE * MAT_POINTS * 
       3               7                0
* INITIAL_STEP * DATA_TYPE * INITIAL_PHASE *
     -1
* NO_STEPS * INITIAL_TIME(s) *  DELTA_TIME(s) *
    401          0.             0.5e-3
* INITIAL_STEP * LAST_STEP * STEP_INTERVAL * TIME_DIFF *
     1              400          1                3         
* INPUT_MESH_FILE *  UNIT  *
        5             0
* POST_DATA_FILE * ELEM_OUT * NODE_OUT *NUMBER_OUTPUT_MATS * AVERAGE *
        5               1        1           0                  0
* MESH * CURRENT * MAGNETIC * FORCE_J_B* FORCE_NODAL * DISP * ELEM * HEAT *
    1        0         1            0          0         1            0
*   MESH   *   A   *   V   *   B   *   B_INTEG *    J_Q  *
     0         0       0       0         0           0
*  Q_AVERAGE * HEAT  *  MAG_FLUX  * CUR_FLUX *
        0        0            0            0
* FORCE_J_B * FORCE_NODAL *
        0         0
* GEOMETRY * DELTA_Z(m) or DETLA_THETA(deg)* NO_LAYERS *
     2         1.0
*COORDINATE * NO_COORDINATES *
COORDINATE         0
* FAR_BOUNDARY_CONDITION 
        2
* NO_BE_TERMS * BE_CENTER X(m), Y(m), Z(m) *
       3                   0.     0. -0.01                       
* DISTANCE_JUDGE(m) *
     1.e-6
* NO_BN0_PLANES *
    0
* NO_HT0_PLANES *
    0
* NO_A0_LINE *
    0
* CYCLIC     *  SYMMETRICITY *  SLIDE_NDIV *  SLIDE_MAT *
     0               1             0          11
* NO_MAT_IDS *  EXTEND_TOTAL *  NO_SMAT_IDS 
     6               0            0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID * SIGMA * MU  * PACKING *
     1            0           0         0.     1.0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID * SIGMA * MU  * PACKING *
     2            0           1         2.e5   1000.0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID * SIGMA * MU  * PACKING *
     3            0           1         2.e5   1000.0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID * SIGMA * MU  * PACKING *
     4            0           0         0.0    1.0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID * SIGMA * MU  * PACKING *
     5            0           0         0.0    1.0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID * SIGMA * MU  * PACKING *
     6            0           0         0.0    1.0
*SOURCE *
SOURCE  
*ELMCUR  *SERIES_NO *NO_MAT_IDS * OPTION * IN_ROTOR *
ELMCUR       1          1          1
* MAT_ID  * IN_SURFACE  * OUT_SURFACE  *  CURRENT *    SIGMA  *
    4          1              2           2.48647e6     0.0
*NETWORK * REGION_FACTOR *
NETWORK      360.
*C   * ID * NODE1 * NODE2 * CAPACITANCE * 
C      1      1      2        0.0001       
*SETV  *  ID  * INITIAL_VOLTAGE *
SETV      1         100.
*FEM * ID * NODE1 * NODE2 * SERIED_ID *
FEM    2      2      3        1
*R   * ID * NODE1 * NODE2 * RESISTANCE *
R      3      3      1        70.
*TABLE  * NUMBER *
TABLE    1
*ID  * NO_DATA *
   1    34
* CURRENT * VOLTAGE *
-0.0016	150.560517 
-0.0008	70.560517 
-0.0004	30.560517 
-0.0002	10.560517 
-0.0001	0.560517 
-0.00009	0.230259 
-0.00008	0.160944 
-0.00007	0.120397 
-0.00006	0.091629 
-0.00005	0.069315 
-0.00004	0.051083 
-0.00003	0.035667 
-0.00002	0.022314 
-0.00001	0.010536 
0	0.000000 
0.01	-0.461512 
0.02	-0.530330 
0.03	-0.570711 
0.04	-0.599396 
0.05	-0.621661 
0.06	-0.639859 
0.07	-0.655251 
0.08	-0.668586 
0.09	-0.680351 
0.1	-0.690875 
0.2	-0.760140 
0.3	-0.800670 
0.4	-0.829430 
0.5	-0.851739 
0.6	-0.869968 
0.7	-0.885381 
0.8	-0.898732 
0.9	-0.910509 
1	-0.921044 
*TAB * ID * NODE1 * NODE2 *  TABLE_ID *
TAB    4      1           4              1
*SETV  *  ID  * INITIAL_VOLTAGE *
SETV      4         100.
*L   * ID * NODE1 * NODE2 * INDUCTANCE *
L      5      4      2        0.01 
*END OF NETWORK *
END  
*END of SOURCE *
END
* NO_DATA *
      2
* TIME_ID * OPTION *
      1       11
* FUNCTION f(t) * t(sec) *
f(t)=100.;;
* TIME_ID * OPTION *
      2     3
*  INITIAL_POSITION  * INITIAL_VELOCITY  * POSITION_ERROR *
     10.e-3                0.0               0.01e-3
*  LOWER_LIMIT  * UPPER_LIMIT * LOWER_REFLECTION * UPPER_REFRECTION * 
     0.2e-3         12.0e-3            0.7              0.
*    MASS    * CONST_FORCE  * STATIC_FRICTION*
      90.e-3        0.0          0. 
* NO_SPRINGS * NO_DAMPERS * MATH_EQ *
       0          0               1
* FUNCTION F(x, v)  * x(m), v(m/s), F(N) or x(deg), v(deg/sec), F(N/m) *
mass=90.e-3;
g=-9.8;
k=1.e3;
x0=10.882e-3;
F(x,v) = mass*g -k*(x-x0);;
* MOTION_TYPE   *  REGION_FACTOR  *  NO_MAT_IDS * 
     0                    360.         1
* MAT_IDS *
   2
* MOTION *  NO_MOTIONS *
MOTION           1
*  MOTIN_ID  * X_TIME_ID * Y_TIME_ID * Z_TIME_ID * PHI_TIME_ID * COORD_ID *
     1           0           0             2             0
* GLOBAL_MOTION_ID *
     1
* POSITION0  * POSITION1 * NO_RIGID_MAT_IDS *  NO_DEFORM_MAT_IDS * CONTROL_ELM_MAT_ID *
       0.012      0.0002              2                                     1                                   11
* RIGID_MAT_IDS *
     2   6
* DEFORM_MAT_IDS *
     5
* NO_CURVES *
     1
* BH_CURVE_ID * NO_DATA *
     1            11
*   H_AT/M     *      B_T    *
0.0 		0.00 
500.0 		0.90 
750.0 		1.10 
1000.0 	1.21 
1750.0 	1.40 
3000.0 	1.53 
5000.0 	1.65 
7500.0 	1.71 
17500.0 	1.90 
30000.0 	2.01 
101619.0 	2.10 