*2
*2
*V
*Position
*I
*Force
*RLS_NO
r10.1
* PRE_PROCESSING * MAKE_SYSTEM_MATRICES *
       1                  1          
* SOLVE_EQUATION * POST_PROCESSING *
       1                  1
* STATIC * STEP * AC * TRANSIENT * DISPLACEMENT * NON_LINEAR * TEMP_DEPEND *
    0       0     0       1           3              1
* POTENTIAL * NODAL_FORCE_OPTION *
    0
* TREE_GAUGE *
   -1 
* ORDER_OF_SHEPE_FUNCTIONS *
         1
* ICCG_CONV * ACCEL_FACTOR * DIV_FACTOR * DIV_ITERATIONS * SOLVER * THETA * THETA_NETWOEK * THETA_MOTION *
    1.e-6        1.            100.          100                0        1            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                1
* 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 *
        0               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) * Distance for judgement (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 * ANISOTROPY * IRON_LOSS 
     1        0           0                    0             1         1    0      0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID   *  SIGMA  *       MU       * PACKING * ANISOTROPY * IRON_LOSS 
     2        0           1               2e+005          1000         1    0      0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID   *  SIGMA  *       MU       * PACKING * ANISOTROPY * IRON_LOSS 
     3        0           1               2e+005          1000         1    0      0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID   *  SIGMA  *       MU       * PACKING * ANISOTROPY * IRON_LOSS 
     4        0           0                    0             1         1    0      0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID   *  SIGMA  *       MU       * PACKING * ANISOTROPY * IRON_LOSS 
     5        0           0                    0             1         1    0      0
* MAT_ID  *  POTENTIAL * B_H_CURVE_ID   *  SIGMA  *       MU       * PACKING * ANISOTROPY * IRON_LOSS 
     6        0           0                    0             1         1    0      0
*SOURCE *
SOURCE  
*ELMCUR  *SERIES_NO *NO_MAT_IDS *CURRENT_OP * 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.000000
FEM     1     1     2     1
VPS     2     2     1     2
SETI     2  0.0000000
*END OF NETWORK *
END  
*END of SOURCE *
END
* NO_DATA *
      3
* TIME_ID * OPTION *
      1       11
* FUNCTION f(t) * t(sec) *
f(t)=100.;;
* TIME_ID *	OPTION	*
	    2         4	 
* PSIM_ID * PSIM_OUT_ID *
          1                1
* TIME_ID * OPTION *
      3     3
* INITIAL_POSITION * INITIAL_VELOCITY * POSIITON_ERROR * PSIM_ID * PSIM_FORCE * PSIM_POSITION *
 0.0100000                        0.0000000         0.0000100            2                2                  0
*  LOWER_LIMIT  * UPPER_LIMIT * LOWER_REFLECTION * UPPER_REFRECTION 
     0.2e-3         12.0e-3            0.7              0.
*     0.2e-3         12.0e-3            0.              0.
*    MASS    * CONST_FORCE  * STATIC_FRICTION*
*     90.e-3      -0.882           0. 
      90.e-3        0.0          0. 
* NO_SPRINGS * NO_DAMPERS * MATH_EQ *
       0          0               1
* SPRING_TYPE  * LOWER_RANGE * UPPER_RANGE * 
*       0          -999.          999.
* SPRING_CONST *  EQUIL_POSITION * 
*   1.e3           10.882e-3   
* 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);;
* FORCE_TYPE   *  REGEON_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             3             0
* GLOBAL_MOTION_ID *
     1
*  POSITION0   *  POSITION1  * NO_RIGID_MAT_IDS *  NO_DEFORM_MAT_IDS * CONTROL_VOLUME_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 