·Table of Contents ·Computer Processing and Simulation | The Numerical Modelling and residual field computation of flat ferromagnet at its magnetization by inhomogeneous pulse magnetic fieldMatyuk V.F., Churilo V.R., Strelyukhin A.V.Institute of Applied Physics of National Academy of Sciences of Belarus Minsk, Republic Belarus Contact |
By model construction the next assumptions were made: the test piece initially demagnetized is, the material isotropic is, the magnetic ductility absent is. The dependence of magnetic properties on the external field by the known approximations for static magnetic characteristics were determined. By solution the discrete model as in the space (the model of ferromagnet with piece-permanent approximation of magnetization on the elements of division), as in time is used.
The general rules of distribution of magnetic field in the article by the system of Maxwell equations can be described the solving of what in our case to the nonlinear vector equation of parabolic type can be reduced
(1) |
where the magnetic field strength is; m_{d}(H) the incremental permeability is;l the electrical conductivity of the material of disk is.
Because the cylindrical symmetry of the problem its further solution in the cylindrical coordinates system (H_{j} = 0) was made. Then for components (H_{r},H_{z}) the equation (1) becomes
(2) |
As in the process of pulse magnetization the field on the surface of ferromagnet not only by the field of solenoid, but by the field of magnetization of the article itself and by the field of eddy currents determined is it was suggested to use the next initial and boundary conditions allowing to take into account the contribution in the field on the boundary of the three above discussed components.
The initial condition
(3) |
The boundary conditions
(4) |
where G^{*} is the portion of the boundary G of the rectangle that haven't the side r = 0; H^{*}_{r}, H^{*}_{z} is the field on the boundary G^{*} including the field of solenoid ^{c}, field of ferromagnet ^{M} and field of eddy currents ^{d }.
Consider that the field of solenoid ^{c} at every instant of time t_{l} known is, and the components ^{M} and ^{d } from the solution (2) on the preceding time interval t_{l-1} determined are.
Then
(5) |
For determination the magnetic state of separate elements of division in dependence on the increase or the decrease of the field strength H in that element the equations for m _{d}(H) from paper [3] were used.
The components ^{c} in the elements of division by the Biot-Savare-Laplace Law [4] determine.
For computation the value of magnetization in each element of division in every time interval the nonlinear integral equation [4] was solved
(6) |
where κ the magnetic susceptibility is; V_{M} the volume of the magnetized material (element) is; Q the observation point is; N the source point (the point determining the coordinates of element of division of ferromagnet) is; the radius-vector from the source point (N) into the observation point (Q) is.
In connection with the nonlinear dependence of magnetic susceptibility on the magnetic field strength the solution (6) is realized as reduction to the system of linear equations and using the method of iteration.
The computation of contribution of magnetic field strength on the boundary due to the eddy currents on the base of finite differential approximation of equation is performed
(7) |
and the components of magnetic field using the Biot-Savare-Laplace Law were determined.
On the computed in the end of pulse remanent magnetization values using equations
(8) |
where Z_{1}, Z_{2}, R_{1}, R_{2} the geometrical coefficients are; n = 1,2,...N_{1} the number of the elements volume is,
the components of the residual field strength in the observation regions were determined.
Fig 1: The distribution of normal H_{rz}(O) and tangential H_{r r}() components of residual field at pulse magnetization along the disk surface (the number of division was 5´50) in the obser-vation region from the side of magnetization (computation --, experiment -O-, --):
H_{m} = 5·10^{5} A/m; a) - t_{0} = 0,85 ms; t_{3} = 1 ms; b) - t_{0} = 0,85 ms; t_{3} = 10 ms; c) - t_{0} = 0,2 ms; t_{3} = 6 ms; d) - t_{0} = 1,2 ms; t_{0} == 6 ms; H_{m} is the pulse amplitude; t_{0} , t_{3} are the duration of the leading and the trailing edge of magnetizing pulse |
The received results show that the considered model the real distribution of residual field magnetization after applying the nonuniform pulse magnetic field to the flat ferromagnet in the form of a disk represents. For more precise quantitative correspondence the further development of the method of computation, the increase of the number of division of space and time grid, specification of more precise computation of iteration procedures and use of more powerfull PC are required.
Thus the general rules of residual magnetization field distribution computed using the suggested method in type and in value coincide with the experimental data with acceptable for practice accuracy. This confirms the legitimacy of use the developed model for computation the magnetic state of ferromagnetic disk in the nonuniform pulse field of the superimposed solenoid.
© AIPnD , created by NDT.net | |Home| |Top| |