STRESS FIELD DISTRIBUTION IN MAGNETOELASTIC PRESSURE FORCE SENSOR STRESS FIELD DISTRIBUTION IN MAGNETOELASTIC PRESSURE FORCE SENSOR

Modelling of the stress field distribution in a magnetoelastic sensor core is proposed. The components of stress field are carried out by a numerical solution of boundary–value problem for stress field. The obtained


Introduction
The magnetoelastic pressure force sensor operates on the Villari effect principle. The Villari effect, also called a stress-magnetization, is based on the change of permeability of a ferromagnetic body exposed to a mechanical force, which results in mechanical stress within the body. The permeability increment is proportional to that component of mechanical stress, which is collinear to magnetic field vectors [1].
The modelling of the magnetic field distribution in the core of such a sensor with respect to the variation of the permeability due to mechanical stress requires to determine the stress field distribution.
In the general case, the planar stress field is unambiguously determined by three components, two perpendicular components and a shear component [2]. These components can be appointed by various methods. In this paper a method based on a numerical solution of boundary-value problem is proposed. This task was solved by professional code PDE Toolbox in package MATLAB.

Magnetoelastic pressure force sensor
The arrangement of the magnetoelastic pressure force sensor is depicted in Fig. 1. The sensor core consists of fifty ferromagnetic transformer sheets, thickness of each being 0.5mm. Each sheet contains four holes of the radius of 1 mm. The conductors of primary (10 turns/0.35mm) and secondary windings (8 turns/0.25mm) are led through these holes collaterally. The sensor is assigned for the measuring of nominal pressure force 120 kN that is equivalent to 100 MPa. A simplifying representation of the sensor core is in Fig. 2. If Cartesian coordinate system is applied to such a simplified form of sensor core (see Fig. 2 a nonzero component only in the z-direction. This component produces magnetic field whose vectors have nonzero components only in the directions x and y [3]. Moreover, the dimensions of the sheet are designed in such a way that the sheet can be divided into four squared shapes, each with a circle hole in the middle. A conductor (the cross-section of which is the same as the cross-section of the primary windings conductors) replaces the primary winding. One of such shapes is called the integrative sensor element (Fig. 3left). With respect to the symmetry of the sheet it is sufficient to determine the magnetic field distribution in the integrative sensor element only. The stress field distribution is necessary to determine only in a part of the integrative sensor element that is made of a ferromagnetic material (see Fig. 3right).
For the modelling of the magnetic field distribution in the sensor core with respect to the variation of material permeability due to mechanical stress it is necessary to investigate only the xdirection stress ( x ) and y-direction stress ( y ) in the integrative sensor element (see Fig. 4).

PDE system formulation for plane stress
The determination of stress field in the integrative sensor element can be formulated as a boundary-value problem in terms of the displacement components. The first of all it is essential to formulate the PDE system for balance of force in terms of the displacement components with the boundary conditions.
For plane stress, in the body made of the homogeneous isotropic material, the balance of force equations are [2]: , where x is the x-direction stress, y being the y-direction stress, xy being the shear stress and X, Y being the volume forces.
The stress components are closely related to the strain components and these relations are defined by Hooke law (3) where E is Young's modulus, v being Poisson's ratio, The strain components can be expressed by the displacements, u the displacement in the x-direction and v the displacement in the y-direction (4) where x is the x -direction strain, y being the y -direction strain, ␥ xy being the shear strain.
Combining the equations (1) -(4) and assuming that there are no volume forces, we arrived at the PDE system for balance of force in terms of the displacement components. The system takes form From mathematical point of view, the equations (5) and (6) represent a system of two PDE of elliptic type.
The following boundary conditions for the system PDE can be given: G For each point on the boundary is specified the value of the displacement This boundary condition is called Dirichlet boundary condition.  This boundary condition is called Neumann boundary condition.

Boundary-value problem formulation
The boundary-value problem formulation also requires determining the domain, in which the stress field will be solved. The integrative sensor element has the squared shape with a circle hole in the middle (see Fig. 4). It can be considered to be a planar body with the following properties: dimension 12-by-12 mm, thickness 0.5 mm and the radius of the circle hole 1 mm. The element is made of ferromagnetic transformer sheet, whose parameters are: E ϭ 1.86.105 MN/m 2 , v ϭ 0.3. The element is exposed to a constant continuous pressure stress p ϭ Ϫp x e x (e x being the unit vector in the direction of x-axis).
The domain ⍀ u,v , in which the stress field will be solved, is surrounded by five boundaries Fig. 5). The boundary conditions for the displacement components are following: G The boundary ⌫ (2) u,v is subjected to a pressure p ϭ Ϫp x e x : , .
The investigated boundary-value problem was solved by using a professional code PDE Toolbox. The toolbox is based on the finite element method (FEM). The solution was made by using the adaptive mode option, which enables to refine the mesh in areas where the gradient of the solution is large in order to increase the accuracy of the solution [7]. The obtained results (for generated triangular mesh consisting of 42 112 nodes and 83 456 triangles) at the pressure force p x ϭ 8,23 MPa are depicted in Figs. 6 -11. In these figures are plotted the contour lines, e.g. the lines connected points of equal value of the depicted quantity in the domain ⍀ u,v . In Figs. 6 and 7 are visualised the displacements u(x, y) in the xdirection and v(x, y) in the y -direction. The x -direction strain and the y -direction strain are depicted in Figs. 8 and 9. The xdirection stress and the y -direction stress are shown in Figs. 10 and 11. Fig. 6 The displacement u in the x -direction Fig. 7 The displacement v in the y -direction Fig. 8 The strain x in the x -direction Fig. 9 The strain y in the y -direction Fig. 10 The stress x in the x -direction Fig. 11 The stress y in the y -direction

Conclusion
The paper is aimed at assignment of stress field distribution in the magnetoelastic pressure force sensor. The task was solved as the boundary-value problem in terms of the displacement components by the professional code PDE Toolbox in package MATLAB. The toolbox provides a powerful and flexible environment not only for the solution of elliptic, parabolic and hyperbolic partial differential equations, but it also handles the solution of a system of partial differential equations. The normal components of the stress field and also the components of the strain and displacement were visualised.