The present disclosure relates to an analysis of the surface mounted permanent magnet (SPM) machines by using 2D meshless method, which belongs to the field of electromagnetic computation.
Permanent magnet (PM) machines have attracted more and more attention due to high reliability, high efficiency and high-power density. PM machines have been widely used in high-end fields such as electric vehicles, and aerospace. Meanwhile, an eccentric rotor shape is widely used in a SPM machine due to its excellent performance. Therefore, an appropriate electromagnetic analysis method is very important, as it directly affects design efficiency and operational performance. Mesh based numerical methods, such as finite difference method (FDM), finite element method (FEM) and boundary element method (BEM), have been quite sophisticated and successfully solved many engineering problems. However, these methods are not perfect in all cases, especially when dealing with motion and deformation. In these cases, the mesh may deform severely, which significantly affects the accuracy of the solution. It is necessary to reconstruct the mesh in these situations. However, reconstruction in complex areas is relatively difficult. In order to solve these problems, meshless methods have gradually developed. Unlike mesh-based methods, meshless methods only focus on the information of nodes. There is no connection between nodes, which is conducive to deal with moving components. Therefore, the meshless method has great potential and advantages in solving the transient magnetic field of PM machines.
The purpose of the present disclosure is to provide a method for analyzing the electromagnetic of the SPM machines, which mainly includes discretizing the solution area of the machine by nodes layout, constructing the support region, converting the partial equations into linear equations to solve. To achieve the above objects, the technical solution adopted by the present disclosure is a meshless method for analyzing SPM machines, which includes the following steps:
Furthermore, in step 1, nodes are arranged in each solution region of the machine, at the interface between regions and at the boundary; The sub regions to be solved include stator, slot, airgap and PM; Due to the fact that the rotor core of SPM machine is generally unsaturated, the Neumann boundary can be applied for the outer surface of the rotor. Thus, there is no need to solve the rotor core region, which can improve computational efficiency.
Furthermore, in step 2, it is necessary to select any point within the sub region as the central node and find a certain number of adjacent nodes closest to the central node to form the support region; All nodes that make up the support region must be within the same sub region.
Furthermore, the specific process of step 3 is as follows: the second-order Taylor expansion is performed at the central nodes for the all nodes in the support region. Then the residual function is constructed, and the corresponding algebraic equations are obtained according to the extreme value principle; The derivative values of node can be expressed as a linear combination of function values of nodes in the support region by solving the algebraic equations.
Further, in step 4, the partial differential equation satisfied by each sub regions are converted into algebraic equations. The PM region and the air gap region satisfy the Laplace's equation. The slot region satisfies the Poisson equation. And the stator core satisfies the nonlinear partial differential equation.
Furthermore, in step 5, two support regions are formed for the nodes at the interface. Then the corresponding equations can be obtained based on the continuity conditions. The nodes distributed on the rotor face satisfy the Neumann boundary condition. The nodes distributed on the stator face satisfy the Dirichlet boundary condition.
Furthermore, in step 6, a System of Linear Algebraic Equations can be obtained; The coefficient matrix of the algebraic equations depends on the permeability, node coordinates and weight function; The source vector of the algebraic equation depends on the current density in the winding and the remanence of PM.
Furthermore, in step 7, the calculations are carried out at different rotor positions to obtain the flux leakage of each tooth in an electrical cycle. The flux line and flux density of each node can be further obtained based on the vector potential of each node by solving algebraic equations. According to the electromagnetic calculation constraints of the machine, it is possible to obtain the flux linkage that flows through each stator tooth of the machine at one time during the electrical angle period. The next rotor position is recalculated to obtain the flux of each tooth during the electrical angle period. Based on this, the three-phase flux linkage and electromotive force are obtained. If it is under load, it can be used to calculate the machine output torque.
The SPM machine of the present disclosure is a 12 slot/10 pole three-phase machine, which can be divided into four parts: stator, air gap, rotor, and shaft; The stator includes the stator yoke, stator teeth, stator slot, and armature winding; The armature slot is a flat bottom slot; The fractional slot concentrated winding configuration is adopted; The rotor is cylindrical and PMs are attached to its surface. The material of PM is neodymium iron boron grade N42UH; The shape of the PM is sector shaped and evenly distributed on the circumference of the rotor; The material of stator core and rotor core is silicon steel sheet DW310_35; The air gap is between the stator and rotor, with a thickness of 1.5 mm; The machine shaft is made of nonmagnetic material, which is a solid cylindrical.
The present disclosure has the following beneficial effects.
The following part will provide a clear and complete description of the technical solution in the present disclosure in conjunction with the drawings.
In order to explain the beneficial effects of the present disclosure more simply and clearly, a specific SPM machine will be described in detail below;
As shown in
The specific process is: the Taylor expansion is performed at the central node for all nodes in the support region. Then the expression of the remainder is obtained. Then the residual function can be achieved by multiplying the weight functions give the residual function. Finally, the corresponding algebraic equation is obtained according to the extreme value principle, and the derivative value of each node can be expressed as a linear combination of the node function values in the support region by solving the algebraic equation.
As shown in
where A0 represents the vector potential value of the center node, An represents the vector potential value of the nth node around the center node, Δxn=xn−x0 represents the difference in x coordinates, Δyn=yn−y0 represents the difference in y coordinates, and o(p3) represents a higher-order infinitesimal, So, the residual function can be defined as:
where Ao represents the vector potential value of the center node, An represents the vector potential value of the nth node around the center node, Δxn=xn−x0 represents the difference between the x-axis of the nth node around it and the central node, Δyn=yn−y0 represents the difference between the y-axis of the nth node around it and the central node, p=(Δxn2+Δyn2) represents the distance between two nodes, and w is the weight function. If the surrounding nodes are further away from the center node, the weight function value will be lower. Obviously, the residual function value should be as small as possible, according to the extreme value principle. There is:
where, R represents the residual function defined by (2). In this way, a set of algebraic equations can be obtained:
By solving (4), the derivative values of node (xo,yo) can be expressed as a linear combination of function values of nodes in the support region:
D
5×1=(F5×5−1E5×(N+1))U(N+1)×1=K5×(N+1)U(N+1)×1 (5)
The PM region and air gap region satisfy the Laplace's equation. The slot area satisfies the Poisson equation. And the stator core satisfies the two-dimensional nonlinear partial differential equation.
The vector potential of PM region 1-4 and air gap region 1-5 in
where A represents the vector potential, and x and y represent coordinate variables, respectively;
The discrete format of formula (6) can be written as:
where A, represents the vector potential value of the i-th node in the support region, k3,i+1 represents the element in the third row and i+1 -th column of matrix K in formula (5) of step 3, and k4,i+1 represents the element in the fourth row and i+1 -th column of matrix K in formula (5). The slot regions 1-3 satisfy Poisson's equation:
where, A represents the vector potential, x and y represent coordinate variables, respectively, μ0 is the magnetic permeability in vacuum, and J is the current density in the slot. Similarly, its discrete format can be written as:
where Ai represents the vector potential value of the i-th node in the support region, k3,i+1 represents the element in the third row and i+1-th column of matrix K in formula (5) of step 3, k4,i+1 represents the element in the fourth row and i+1-th column of matrix K in formula (5), μ0 is the magnetic permeability in vacuum, and J is the current density in the slot. The stator core region meets the second-order nonlinear partial differential equation:
where A represents the vector potential function, x and y represent coordinate variables, and v is the magnetic permeability in the iron core. Similarly, its discrete format can be written as:
where Ai represents the vector magnetic potential value of the i-th node in the support region, vi represents the magnetic permeability value of the i-th node in the support region, v0 represents the magnetic permeability value of the center node, k1,i+1 represents the element in the first row and i+1-th column of matrix K in formula (5), k2,i+1 represents the element in the second row and i+1 -th column of matrix K in formula (5), k3,i+1 represents the element in the third row and i+1-th column of matrix K in formula (5), and k4,i+1 represents the element in the fourth row and i+1-th column of matrix K in formula (5).
Step 5: The additional processing is required for the nodes at the interface and the boundary. The nodes at the interface need to meet continuity conditions while the nodes at the boundary need to meet corresponding boundary conditions.
For nodes at the interface between two sub regions, it is necessary to use this node as the central node to construct the support region in each sub region. Then the corresponding equations based on the continuity conditions of magnetic field are obtained. The nodes on the rotor boundary satisfy Neumann boundary condition. The nodes on the outer surface of the stator satisfy the Dirichlet boundary condition.
Due to the discontinuity of the derivative of the vector potential at the interface between two regions, the nodes distributed there need to be processed according to continuity conditions. Firstly, as shown in
The vector potential at the interface of two permanents magnets satisfies the following relationship:
where nx and ny are the x-direction component and y-direction component of the unit tangent vector at the interface, AI represents the vector potential in the first PM, AII represents the vector potential in the second PM, μ is the permeability of the PM, and Hcx and Hcy are the x-direction and y-direction components of the coercive force of the PM, respectively. According to step 4, the derivative of equation (12) is expressed as a linear combination of function values in the support region, so the discrete format can be expressed as:
where AI represents the vector potential in the first PM, AII represents the vector potential in the second PM, μ is the permeability of the PM, Hcx and Hcy are the x-direction and y-direction components of the coercive force of the PM, respectively, kI1,i+1 represents the coefficient in the first region which determined by the element in the first row and i+1-th column of matrix K in formula (5), kI2,i+1 represents the coefficient in the first region which determined by the element in the second row and i+1-th column of matrix K in formula (5), represents the coefficient in the second region which determined by the element in the first row and i+1-th column of matrix K in formula (5), kII2,i+1 represents the coefficient in the second region which determined by the element in the second row and i+1-th column of matrix K in formula (5). At the interface between the PM and the air gap region, the vector magnetic potential satisfies the following relationship:
where APM represents the vector potential in the PM, and Aair represents the vector potential in the air gap, μPM is the permeability of PM, μ0 is the air permeability, while Hcx and Hcy are the x-direction and y-direction components of the coercive force of the PM, respectively. Its discrete format can be written as:
APM,i represents the vector potential value of the i-th node in the PM, Aair,i represents the vector potential value of the i-th node in the air gap, μPM is the permeability of pm, μ0 is the air permeability, Hcx and Hcy are the x-direction and y-direction components of the coercive force of the PM, respectively, kPM,1,i+1 represents the coefficient in the PM region, which is determined by the element in the first row, i+1 column of matrix K determined by formula (5), kPM,2,i+1 represents the coefficient in the PM region which determined by the element in the second row, i+1-th column of matrix K determined by formula (5), kair,1,i+1 represents the coefficient in the air gap region which determined by the element in the first row, i+1 column of matrix K determined by formula (5), kair,1,i+1 represents the coefficient in the air gap region which determined by the element in the second row, i+1-th column of matrix K determined by formula (5); At other interfaces, the vector potential satisfies the following relationship:
where nx and ny are the x-direction component and y-direction component of the unit tangent vector at the interface, AI and AII are functions of vector potential in regions I and II, μ1 and μ2 are the magnetic permeability of the two solving regions. The discrete format of equation (16) can be expressed as:
where nx and ny are the x-direction component and y-direction component of the unit tangent vector at the interface, AI and AII are functions of vector potential in regions I and II, μ1 and μ2 are the magnetic permeability of the two solving regions, kI1,i+1 represents the coefficient in the first region which determined by the element in the first row, i+1-th column of matrix K determined by formula (5), kI2,i+1 represents the coefficient in the first region which determined by the element in the second row, i+1 column of matrix K determined by formula (5), kII1,i+1 represents the coefficient in the second region which determined by the element in the first row, i+1-th column of matrix K determined by formula (5), kII2,i+1 represents the coefficient in the second region which determined by the element in the second row, i+1-th column of matrix K determined by formula (5); For nodes at the boundary, corresponding boundary conditions can be given. Dirichlet boundary condition can be applied for the nodes surface of the stator 1-1:
Astator=0 (18)
where Astator represents the vector potential value of the outer surface of the stator. Neumann boundary condition can be given for the nodes distributed on the surface of rotor 1-2.
where n represents the unit tangent vector of the boundary, μPm represents the permeability of PM, Arotor represents the vector potential of the outer surface of the rotor, and Hc represents the coercive force; Its discrete form can be expressed as:
where nx and ny are the x direction component and y direction component of the unit tangent vector on the outer surface of the rotor respectively, μPM represents the permeability of the PM, Arotor,i represents the vector potential value of the node on the outer surface of the rotor, Hcx and Hcy are the x-direction and y-direction components of the coercive force of PM, k1,i+1 represents the element in the first row and i+1 -th column of matrix K in formula (5), and k2,i+1 represents the element in the second row and i+1 -th column of matrix K in formula (5).
The coefficient matrix G of the algebraic equations depends on the permeability, node coordinates and weight function; The source matrix S depends on the current density in the slot area and the magnetization of the magnet.
According to steps 4 and 5, operating on each node can get an algebraic equation. By combining them together, a matrix equation can be obtained:
GA=S (21)
where G is the coefficient matrix, which depends on the node coordinates, weight function, and magnetic permeability; A is the column vector of vector potential to be solved, and S is the source matrix, which depends on the current density and PM magnetization. Since the permeability in the iron core is not constant, the equations (21) are nonlinear algebraic equations. Equation (21) can be solved using the successive linearization method. Its iterative scheme can be expressed as:
v
i=(1−f1)vk+1+f1vk+1(Bk−1)
A
k=(1−f2)Ak−1+f2(Gk)−1Sk (22)
where Ak is the vector potential vector obtained from the kth calculation, Ak−1 is the vector potential vector obtained from the k−1st calculation, f1 and f2 are relaxation factors used to control the convergence speed, vk is the permeability value used for the kth calculation, vk−1 is the permeability value used for the k−1st calculation, vk−1 (Bk−1) is the permeability value obtained from the B-H curve based on the flux density Bk−1 calculated from the k−1st calculation.
Step 7: Based on the vector potential of each node solved in Step 6, the flux density distribution and the flux lines can be obtained; According to the electromagnetic calculation constraints of the machine, electromagnetic parameters such as the back electromotive force and electromagnetic torque of the machine can be calculated.
The vector potential of each node can be obtained through step 6, and the flux density of each node can be further calculated through vector potential; The flux line can be obtained by drawing the contour line of the vector potential; By using the node vector potential difference, the flux through the stator teeth of the machine at a predetermined time during the electrical angle cycle can be obtained; Based on this, the electromagnetic parameters such as three-phase flux linkage and induced electromotive force of the machine are obtained; If it is under load, it can be used to calculate the output torque of machine. In order to verify the accuracy and reliability of the 2D meshless analysis method for analyzing SPM machine proposed by the disclosure, the results obtained by the FEA commercial software are used for comparative verification.
In summary, a 2D meshless method for analyzing the SPM machine of the present disclosure includes discretizing the solution region of the machine by nodes distribution; For each discrete node in the region, the closest nodes to the central node to are searched form a support region; Based on Taylor expansion and weighted least squares principle, the derivative values of each discrete node's vector potential can be approximated as a linear combination of the vector potential values of each node in the support region; In this way, the partial differential equation satisfied by vector potential can be converted into algebraic equations; For nodes at the interface between two regions, the additional processing based on continuity conditions is required; The nodes at the boundary will meet the corresponding boundary conditions; The algebraic equations can be obtained by performing such operations on nodes; The number of algebraic equations is equal to the number of discrete nodes; The vector potential of each node can be obtained by solving the algebraic equations, and then the distribution of flux lines and flux density can be obtained. According to the electromagnetic calculation constraints of the machine, the parameters such as the back electromotive force and electromagnetic torque of the machine can be obtained. The present disclosure conducts meshless modeling analysis for SPM machines, and the proposed scheme can provide reference research for this type of SPM machines.
Although the present disclosure has been disclosed in preferred embodiments, the embodiments are not intended to limit the present disclosure. Without departing from the spirit and scope of the present disclosure, any equivalent changes or embellishments made shall fall within the scope of protection limited by the claims attached to this application.
Number | Date | Country | Kind |
---|---|---|---|
202111160699.4 | Sep 2021 | CN | national |
This application is the national phase entry of International Application No. PCT/CN2021/125118, filed on Oct. 21, 2021, which is based upon and claims priority to Chinese Patent Application No. 202111160699.4, filed on Sep. 30, 2021, the entire contents of which are incorporated herein by reference.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2021/125118 | 10/21/2021 | WO |