Abstract
The solenoid valve used on the hydraulic roof support in the coal mine is an intrinsically safe solenoid. It requires the solenoid valve to achieve large electromagnetic force at low current. The load can not be driven if the current is too low, while a large temperature rise and energy loss would occur in the form of heat if the current is too high. In order to solve this contradictory problem, the load force of the solenoid valve was analyzed first. Then the ANSYS models of the solenoid were established to carry out the simulation jobs. The effects of the three parameters, the sleeve length, the seat length and the boss height, on the electromagnetic forces and magnetic fields were researched so as to optimize the three parameters. The static and dynamic characteristics of the optimized solenoid valve were also studied. An experimental setup was established to verify the simulation results finally. Results show that the simulation values are in good agreement with the experimental values. The peak value of starting current and holding current is small and the solenoid responses rapidly. This study proves the accuracy and reliability of the simulation models and the methods. It also provides valuable references for the design and optimization of the intrinsically safe solenoid valve.
Introduction
Solenoid valve is an important electro-mechanical conversion component which is widely used for automatic control in the aerospace equipments, marine equipments, mining equipments and other industrial machines. Generally speaking, the most important indicator of an solenoid valve is the electromagnetic force [1–4] which is affected by many factors such as the current, coil number, volume and so on [5,6]. However, there are strict restrictions on the solenoid valve used in the coal mine for the safety reasons. For example, the current is not allowed to exceed 200 mA, and the voltage is 12 V, so it requires large electromagnetic force at low power and is difficult to design the intrinsically safe solenoid valve.
Many scholars investigated on the force calculation of the solenoid valve in the past. Some scholars used the theoretical model to calculate the electromagnetic force and found ways to increase it, but the models were simplified a lot and it leaded to big errors [7,8]. Some others used double coils to increase the electromagnetic force [9], but this type of solenoid valve is much larger and takes up more space, so it is not suitable for the hydraulic support used in the coal mine. Some ones used the materials with good magnetic permeability to increase the electromagnetic forces [10–12], but it is hard for these materials to meet the harsh environment in the coal mine because it is wet and corrosive. Furthermore, with the modern development of computer analysis technology, many scholars use the finite element method to improve the distribution of magnetic fields in order to shorten the design process, reduce the cost and improve the design efficiency [13–16]. However, most people focused too much on the magnetic force, and they ignored the static and dynamic characteristics of the solenoid.
In this paper, the finite element method is used to optimize the key structural parameters, and the static and dynamic characteristics of the optimized solenoid are also studied.
The models
Structure and force analysis
Figure 1 shows the electro-hydraulic valve used on the hydraulic support in the coal mine. The electro-hydraulic valve is consisted of two parts, one is the solenoid valve, and the other one is a hydraulic valve. The output force of the solenoid valve can make the hydraulic valve open and close. The detailed working principle is as the following: when the coil of the solenoid is energized, the armature moves to the right under the action of the electromagnetic force and pulls the lower part of the lever, then the output force of the solenoid valve is amplified through the lever. The amplified force, overcoming the resistant forces such as the spring force of the hydraulic valve, makes the moving components move to the right. Thus, the pressure port P is connected to port A. In this case, the hydraulic valve begins to work. On the contrary, when the coil of the solenoid valve is de-energized, the armature of the solenoid valve and the moving components of the hydraulic valve both move to the initial left side by the respective springs. Then port P is disconnected to port A. The hydraulic valve stops work.

The structure of the electro-hydraulic valve.
The electromagnetic force is so important that it directly depends whether the moving components of the hydraulic valve can move reliably. The opening actions can be divided into two stages as the following.
(1) The initial stage
The moving components begin to move at this stage. The electromagnetic force must satisfy the following inequality (1).
When the moving components move to the right side at the maximum opening state, the electromagnetic force meets the following force inequality (3).
It should be noted that the pressure force is not taken into consideration in the above equations because the pressure force is balanced in the structure design of the hydraulic valve. Above all, a lever of 10 times magnification is used to increase the output force of the solenoid valve. The output force of the solenoid valve should at least be 4 N.
The maximum opening of the hydraulic valve is 0.2 mm, so the output displacement of the solenoid valve should be at least 2 mm. On the basis of this, an intrinsically safe solenoid valve is designed. The two dimensional sectional view of the solenoid valve is shown in Fig. 2. The three parameters, the sleeve length l, the seat length m, the boss height n, affect the magnetic fields and the electromagnetic force, so they would be optimized.

The two dimensional sectional view of the solenoid valve.
The solenoid valve is geometrically axial symmetric. In order to reduce the amount of calculation and the number of iterations, only 1/2 two dimensional model is selected for calculation. Its simulation model in ANSYS is shown in Fig. 3. Then different parts of the solenoid valve should be defined with different material properties. The coil and rod material should be copper which is antimagnetic and can be directly selected in ANSYS. The armature should use the material with good magnetic permeability, so the electrician pure iron is selected. The seat and the cover are selected as SA1010 steel which is a type of low carbon steel. The air gap as well as the external fields is filled with air, of which the magnetic permeability is defined as 1. Then the model is meshed wish triangle and quadrilateral meshes. Next, the load should be applied on the model. Here, the current density, 3330000 A∕m2, is applied to the coil. Then the calculation begins.

The ANSYS model of the solenoid valve.
The purpose of the experimental tests is to measure the static and dynamic characteristics of the solenoid valve. The static characteristics can be reflected by two items, one is the displacement-force curve which indicates the load capacity of the armature at different positions, and the other one is the current-force curve which reveals the load capacity of the intrinsically safe solenoid valve at different excitation currents. The dynamic characteristics can be reflected by the transient current or the transient displacement of the solenoid valve under step voltage excitations. It represents the rapid response capability of the solenoid valve.
The principle of the experimental system is shown in Fig. 4. The solenoid valve is fixed on the supporting platform, and the moving platform driven by the combination of a servo motor and a lead screw-nut mechanism moves to change the relative position between the solenoid valve and the moving platform. A load transducer which can be used to measure the electromagnetic force of the solenoid valve and an eddy current displacement sensor which can be used to test the air gap or the displacement of the solenoid valve, are mounted on the moving platform. A grating scale is used to feedback the position of the moving platform. A current sensor is used to test the transient current of the coil. All the moving actions and the data acquisition are controlled by the NI host computer and the Labview software. Figure 5 shows the corresponding experimental setup. It is noted that each experimental data in this study is the average value of three tested groups.

Principle of the experimental system.

The experimental setup.
Effects of the key parameters
(1) the influence of sleeve length l
In the simulation models, all the parameters are kept the same except the sleeve length l which varies from 0 to 23 mm. The initial electromagnetic forces at different sleeve length are obtained as Fig. 6 shows. It can be seen that the electromagnetic force rises slowly when l increases from 0 to 20 mm and falls down sharply when l increases from 20 mm to 23 mm. The reason can be found out through the analysis of the distribution of magnetic lines of flux at different sleeve length in Fig. 7. It is clear that when l increases from 0 to 20 mm, more and more magnetic lines pass through the sleeve and the armature, so the electromagnetic force rises slowly. But when l exceeds 20 mm, the distance between the sleeve and the seat is smaller than the distance between the armature and the seat. So many magnetic lines pass through the gap between the sleeve and the seat rather than through the working air gap between the armature and the seat. Therefore, it is a matter of course that the electromagnetic force becomes sensitive to the sleeve length when l > 20 mm and falls down rapidly with the increase of l. So it can be conclude that the electromagnetic force reaches to the maximum value which is 1.7 N when l = 20 mm.

Electromagnetic force at different sleeve length.

The distribution of magnetic lines of flux at different sleeve length.
Suitable distance between the sleeve and the seat not only affect the distribution of the magnetic lines of flux, but also brings the largest electromagnetic force.
(2) the influence of seat length m
There are 11 lengths selected as the simulation models in this part. The seat length m ranges from 12.6 mm to 32.6 mm, one type every 2 mm. The obtained relationship between the electromagnetic force and the seat length in ANSYS is revealed in Fig. 8. It is clear that the electromagnetic force increases with the seat length first and then decreases. When m = 22.6 mm, where the working air gap is in the middle of the coil, the electromagnetic force is the highest which is 4.7 N. Figure 9 shows the three maps of distribution of magnetic lines of flux. It can be seen that when m increases from 12.6 mm to 32.6 mm, the number of leakage magnetic lines of flux from the armature increases while number the leakage magnetic lines of flux from the seat decreases, and the total number of leakage magnetic lines of flux decreases first and then increases, so the total number of magnetic lines of flux passing through the working gap increases first and then decreases. When m = 22.6 mm, the total leakage flux is the least and the magnetic flux passing through the working air gap is the most, so the electromagnetic force is the biggest.

Electromagnetic force at different seat length.

The distribution of magnetic lines of flux at different seat length.
(3) the influence of boss height n
14 sets of model are chosen to carry out the simulation when the boss height varies from 0.5 mm to 3 mm. The initial electromagnetic forces at different boss height are obtained as Fig. 10 shows. It can be seen that the force increases uniformly with the increase of the boss height. The forces reaches the maximum value when n = 2.4 mm. Then the force drops rapidly when the boss height keeps on increasing. The truth can be revealed through Fig. 11 which the distribution of magnetic lines of flux at different boss height. As the boss height increases, more and more magnetic fluxes pass through the boss and the main working air gap while the number of magnetic fluxes passing through the bottom of the boss becomes less, so the output force increases reasonably. When the boss height goes on increasing, the right side corner of the boss is so close to the seat that many magnetic flux lines pass from the sharp corner to the seat, and the number of the magnetic flux lines passing through the working air gap decreases, so it is a matter of course that the output force drops rapidly. It can be concluded that the electromagnetic force can increases greatly if the boss height increases properly until the right side of the boss is close the seat.

Electromagnetic force at different boss height.

The distribution of magnetic lines of flux at different boss height.
Figure 12 is relationship between the working air gap and the electromagnetic force under different currents. It can be seen from the figure that the simulation results agree well with the experimental results, which proves that the simulation method and the results are reliable. On the whole, the higher the current value, the larger the electromagnetic force under the same air gap. However, it does not mean that the higher the current, the better, because bigger current causes serious temperature rise and too much energy to be wasted in the form of heat during the holding stage of the solenoid valve when the hydraulic valve is working, so the lower current should be selected as long as it meets the initial force requirements. It can be seen that as the size of the air gap increases, the electromagnetic force becomes lower. Since the air gap size designed in this study is 2 mm which is the initial state, more attention should be paid to the situation at 2 mm. It can be seen from the enlarged view that electromagnetic force is greater than 4 N, which meets the force requirement at the initial stage, when the current is 100 mA and 110 mA. When the current is 100 mA, the electromagnetic force is 4.2 N which is a little bigger than the minimum requirement and is safe in practical applications, so 100 mA can be sected as the supplying current. The result is benefit for the circuit design of the coil.

Electromagnetic force at different currents.
The transient responses of current and the displacement of the solenoid valve can be used to evaluate the dynamic characteristics. Figure 13–Fig. 15 are the step responses of current, opening displacement and the closing displacement when the coil is supplied with DC 12 V excitation voltage. It is clear that the experimental values are close to the simulation values, indicating that the simulation models and the methods are accurate and reliable. The maximum current is below 200 mA which is the maximum allowed value for the Intrinsically safe solenoid valve and the response of displacement is rapid. However, in general, the tested current is greater than the simulation current and the experimental displacement lags behind the simulation result, it may be caused by the friction force which is hard to be considered in the ANSYS model. The current in Fig. 13 forms a peak during the rising process due to the excitation impedance of the coil. At last, the current reaches a maximum value and remains. The solenoid valve is in the holding stage during this period. The current is completely dissipated in the form of heat.

Transient current characteristics.

Transient displacement characteristics of the opening process.

Transient displacement characteristics of the closing process.
The ANSYS software is used in this paper to study the effects of the key parameters, such as the sleeve length, the seat length and the boss height, on the electromagnetic force, and these three parameters are optimized. On the basis of this, the static and dynamic characteristics of the optimized solenoid valve are studied through the combination of simulation and experiment. The innovations and contributions of this paper are:
(1) The electromagnetic force does not just simply increase or decrease with the sleeve length, the seat length or the boss height, but it has an optimal value due to the distribution of magnetic fields. Results show that the electromagnetic force reaches the highest value when the three parameter values are 20 mm, 22.5 mm, 2.4 mm, respectively.
(2) The electromagnetic force increases with the current of the coil. When the current is 100 mA, the electromagnetic force can meet the requirement to open the hydraulic valve in the initial state, without too much temperature rise and energy loss at the holding stage.
(3) In the test of dynamic characteristics, the transient test current is slightly higher than the simulation value, and the test displacement lags behind the simulated value. This may caused by the friction force during the tests. However, the simulation results are very close to the test values overall, indicating the accuracy of the simulation model and the methods.
The research results can provide significant references for the design and optimization of the intrinsically safe solenoid.
Footnotes
Acknowledgements
The authors thank the support of Shanxi Science and Technology Department, under project numbers 201701D221139 and 201705D121016. We also thank the support of the National Natural Science Foundation of China, under project number U1610251.
