_{1}

^{*}

As permanent magnet motors and generators produce torque, vibration occurs through the small air gap due to the alternating magnetic forces created by the rotating permanent magnets and the current switching of the coils. The magnetic force can be calculated from the flux density by finite element methods and the Maxwell stress tensor in cy-lindrical coordinates. In this paper the magnetic flux density, the magnetic force and the torque of a real three dimen-sional brushless DC motor are simulated using Maxwell 3 D V 11.1.

Magnetic force analysis is important in determining not only the torque as the output of motor, but also the source of vibration in machine itself. Vibration is induced in permanent magnet DC motors and generators as shown in

The analysis of magnetic force has been addressed by a number of investigators. Initially, papers mainly focused on the calculation of the torque as the output of motor. Marinescu et al. [

In this paper, the characteristics of torque and magnetic force acting on the rotor are analyzed as the rotation of the rotor. The magnetic force in the motor with three poles and six teeth were calculated from the flux density by finite element methods and the Maxwell stress tensor. The component and the characteristics of the magnetic force and magnetic torque in a small air gap were analyzed with the introduction of cylindrical coordinate.

The magnetic field generated by a brushless DC motor is governed by the set of Maxwell’s equations. Introducing the scalar potential into the Maxwell equations, with some mathematics, gives a single partial differential equation for the scalar potential.

The field intensity due to the current in the winding and the permanent magnet may always be calculated directly by the Biot-Savart law and the magnetic dipole moment per unit volume, respectively. This equation, like the Poisson's equation for electrostatic fields can be solved using finite element method. Maxwell 3D V11.1 a FEM solver for magnetic field was used to calculate the scalar potential and the magnetic field intensity.

A non uniform distributed force per unit area at the interface between two materials is calculated by use of the Maxwell stress tensor. Since the strain imposed on the material due to magnetostriction is small enough to neglect changes in the density, it can be assumed that the change of permeability is negligible. Thus in tensor notation we have:

where is the Maxwell stress tensor. is the magnetic flux density which is obtained by the multiplication of permeability to the magnetic flux density. From the expression given by Woodson and Melcher [

The normal and the tangential traction can be decomposed as follows:

Since, the magnetic traction can be simplified with the introduction of the cylindrical coordinate on the rotor. Along the air gap, the normal and the tangential traction for tooth face are as follows:

But for tooth side which is perpendicular to the air gap the normal and the tangential traction have the following form:

The torque produced for one position can easily be derived from the integration of the shear force along the small air gap with the fact the field distribution inside a closed surface in air remain unchanged if the external sources are removed and replaced by currents and poles on the surface:

The finite element method has a solid theoretical foundation. It is based on mathematical theorems that guarantee an asymptotic increase of the accuracy of the field calculation towards the exact solution as the size of the finite elements used in the solution process decreases. For time domain solutions the spatial discretization of the problem must be refined in a manner coordinated with the time steps of the calculation according to estimated time constants of the solution (such as magnetic diffusion time constant). Maxwell 3D V11.1 solves the electromagnetic field problems by solving Maxwell’s equations in a finite region of space with appropriate boundary conditions and-when necessary-with user-specified initial conditions in order to obtain a solution with guaranteed uniqueness. In order to obtain the set of algebraic equations to be solved, the geometry of the problem is discretized automatically into tetrahedral elements. All the model solids are meshed automatically by the mesher. The assembly of all tetrahedra is referred to as the finite element mesh of the model or simply the mesh. Inside each tetrahedron, the unknowns characteristic for the field being calculated are represented as polynomials of second order. Thus, in regions with rapid spatial field variation, the mesh density needs to be increased for good solution accuracy (see also adaptive mesh refinement).

Solving an electromagnetic field problem is always based on solving Maxwell’s equations. However the process of obtaining the solution is typically based on solving a second order consequence of Maxwell’s equations with the consideration of applicable constitutive equations. At the same time-as a rule-a subset of complete Maxwell’s equations is considered according to characteristic aspects of the application. Thus for reasons of efficiency of the solution, applications are classified as electrostatic, magnetostatic, frequency domain or time domain and as a consequence a specific type of solver is used in each case. This allows the users to obtain the solution with the desired accuracy but always within the limits of the fundamental assumptions made when the application was classified along the lines of the above mentioned criteria. The guide lines that can be used to correctly identify the type of solution to use are mentioned in the following paragraphs. The unknowns for each type of solution can be different, depending on the formulation.

The solution process for 3D transient applications poses several challenges not apparent for other solver types. For instance, the behavior of fields is more complex than it is for static or steady-state applications, and a special finite element mesh structure is required in order for the software to accurately represent the physics. For 3D transient applications, the behavior of the fields is more complex than it is for static or steady-state applications. A diffusion of the magnetic field into the materials occurs in 3D transient situations. The distribution of the magnetic field inside objects typically has a number of spatial harmonics, which usually means the time step used in the analysis should be less 7 (sometimes much less) than the magnetic diffusion time constant. These time constants depend upon the geometry of objects and also upon their respective material properties. Since eddy currents are usually considered in conductive objects, a special finite element mesh structure is required in order to accurately capture the physics. In general, a careful planning of the (manual) meshing process is required in order to achieve an accurate solution with the available hardware resources. The three dimensional view of the Brushless DC motor is shown in

It is assumed that the square voltages with amplitude of 311 V and frequency of 60 Hz are applied to the brushless DC motor as shown in

A method for the analysis of the magnetic force and the torque in a brushless DC motor has been presented from the finite element method and Maxwell stress tensor, so that the characteristics of the magnetic force and the torque can be predictable. Shear force existing in a small air gap produces the torque. The permanent magnet and the teeth geometry produces the reluctant torque, which increases the magnetic field of the teeth corner with the direction to move the rotor to the equilibrium position. The commutating torque is produced by the interaction of the permanent magnet and the current which increases the magnetic field concentration of the teeth corner in the moving direction and decreases it in the opposite direction and the simulation software has analyzed nicely the performance of the brushless DC motor in the transient state.