* RLS_NO *
r8.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.            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 *
        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 
     5               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
*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.
*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.
*DIODE  * ID * NODE1 * NODE2 *  V0   *   I0  *
*D2         4      1      4      0.1    0.0001 
*EQ  * ID * NODE1 * NODE2 *
EQ      4      1      4 
* FUNCTION V(I)  * I(A), V(V) *
V0=0.1;
I0=0.1e-3;
IL=I0*(-0.99);
[I<= IL] V1(I)=-V0*( (I-IL)/(IL+I0) +log((IL+I0)/I0) );
[IL <I] V2(I)=-V0* log((I+I0)/I0 );
V(I)=V1(I)+V2(I);;
*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
*   C0    *   C1    *   C2    *   C3     *   C4     *  C5    *   C6    *
*    100.0     0.0       0.0       0.0        0.0       0.0       0.0
*  TEXP   *   TCYCLE  *
*    0.0       0.00
* 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.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             2             0
* GLOBAL_MOTION_ID *
     1
* POSITION0  * POSITION1 * NO_RIGID_MAT_IDS *  NO_DEFORM_MAT_IDS *
   0.0002       0.012            1                     1
* RIGID_MAT_IDS *
     2
* 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