Identification of the Equivalent Permeability of Step-Lap Joints of Transformer E.Napieralska Juszczak Universite Artois, France
COMPUTATION OF THE EQUIVALENT CHARACTERISTICS OF LAMINATED ANISOTROPIC MAGNETIC CORES
Presentation of the structure Layer 1 Layer Layer 3 Layer 4 Air gap 90 0
Method of establishing equivalent characteristics for various structures Using the homogenization technique makes it possible to replace real 3D structures by simpler homogeneous D. In the homogenization technique, two concepts are of importance: macrostructure and microstructure. In the considered problem, the macrostructure is a complete assembly of layers in the overlap, while the microstructure is a repeatable structure of two or more layers made of sheets with the different rolling direction.
The relation between the flux density vector in the repetitive structures r B = dx 1 n r dy H i= 1 B dx dy h i h i The calculation of components Bix and Biy is based on the total flux penetrating through the limiting surfaces of the volume V.
The calculation is based on the total flux penetrating through the limiting surfaces of the volume V. φ y B ix ' + '' = φ ix φ ix dy h i φ x φ x φ y axis x B iy ' φ ' φ iy = + ' iy dx h i
The total energy of the macrostructure results from the sum of the magnetic energy W µi stored in each volume V i of the layer i r r V B H layi i i W = µi r n ν = = V layi r B B i ( ν Fe Fe Fe B, W ) µ min = nb ν + n ν W ( ) B, α i i i 0 0 B 0 r ( µ, ) µ min = Wµi i Bi The unknowns are flux densities vectors Bi in layers and the angles
The starting values for B1x, B1y, Bx and By follow an initial assumption that magnetic flux goes entirely through the laminations avoiding the gaps. Moreover, for small values of induction the flux goes mainly along the rolling direction, whereas for larger values The minimization task was solved using Hook Jeeves' method.
To calculate the equivalent curve B-H the value of flux density B was changed in the interval 0-Bsat A single point of the curve B-H was determined accordmg to the following algorithm: 1) To calculate B1,B,alfa1,alfa B alfa from the minimization task ) To calculate ν 1 ( B 1, α 1 ), ν ( B, α ) 3) To calculate the homogenized reluctivity ν ( r ) [ ( ) ( ) ] ν1 B1, α1 B1 n1 + ν B, α B n B, W = μ min ( n 4) To calculate field intensity. 1 + n ) B min
D model of the core corner 1 3 1 9 1011 5 6 7 5 6 5 4 3 10 11 7 1 14 13 8 3 1 5 6 7 7 6 5
Family of characteristic B-H for the directions 0 for the 4 layers
. Equivalent static B/H curves for the air-gap (y;1c;1g) at four different anisotropy angles B [T] 1.6 14 1.4 1. 1 0.8 0.6 0.4 0. 0 y;1c;1g 0 y;1c;1g 30 y;1c;1g ; 60 y;1c;1g 90 0 500 1000 1500 H [A/m]
The A two-turn experimental search coil for set up measuringmagneticflux alongy axis y 0 0 x 1 0 0 0 0 0 0 0 Layer 1 A two-turnsearch coilfor measuringmagneticflux alongx axis 0 0 0 0 0 0 0
Measured characteristics for the whole structure (before and after drilling of the holes and installation of the search coils). 1.5 1.1 1 111 1.11 1.1 B (T) 1.09 108 1.08 0.5 1.07 60 65 70 75 80 85 original lamination (no holes) after drilling of the holes with search coils in position 0 0 00 400 600 800 1000 100 H (A/m)
3D simulation
Eddy current distributions in the x y cross-section of the lamination
Distribution of Bx along gpath 1 0.40 0.35 0.30 Bx[T] 0.5 0.0 3D D measured 140 130 110 90 path 1 0.15 y [mm] 70 path 0.10 50 30 0.05 0.00 0 0 40 60 80 100 10 140 160 x[mm] 10 0 1T 1 T 10 30 50 71 91 111 131 151 x [mm]
Distribution of Bx along gpath 0.4 0.35 140 130 110 0.3 90 path 1 0.5 y [mm] 70 50 path Bx[T] 0. 0.15 3D D measured 30 10 0 1 T 1 T 10 30 50 71 91 111 131 151 x [mm] 0.1 005 0.05 0 0 0 40 60 80 100 10 140 160 x[mm]
Distribution of By along path 1 0.4 0.35 By[T] 0.3 0.5 0. 0.15 0.1 3D D measured y [mm] 140 130 110 90 path 1 path 70 50 30 1 T 10 0 1T 10 30 50 71 91 111 131 151 x [mm] 0.05 0 0 0 40 60 80 100 10 140 160-0.0505 x[mm]
Distribution of By along path 0.5 0.45 0.4 140 130 110 By[T] 0.35 0.3 0.5 0. 0.15 3D D measured y [mm] 90 path 1 path 70 50 30 1 T 10 0 1T 10 30 50 71 91 111 131 151 x [mm] 0.1 0.05 0-0.05 005 0 0 40 60 80 100 10 140 160 x[mm]
Fluxes in three consecutive layers Φ 1x z x y Φ Φ Φ Φ 1y y 1 sxy y Φ x 3y Φ 3x Φ 3 sxy Φ z Φ Φ Φ z y Φ x d1 (conductor) d (air gap + insulations) d3 (conductor)
The flux density components in a sample of dimensions δx δy δz may be found as: B ix = 1 ' ix φ + δ y δ φ z ' ' ix ; B iy = 1 ' iy φ + δ x δ φ z '' iy ; B iz = 1 ' iz '' iz φ + φ iz iz y x δ δ ;
Models geometry model with two laminations the 90 step lap joint
Model 1 A series of simulations were executed: For excitations with the flux 4e-6 Wb (excitation ) and 35e-6 Wb (excitation 1) the repartitions of the flux density and of the eddy currents were simulated. The real conductivity was taken into account. For the excitation ti with the flux 4e-6 Wb and 35e-6 Wb the simulations neglecting the eddy currents were executed. For excitations with the flux 4e-6 Wb and 35e-6WB the anisotropic conductivity was applied. Zero conductivity in the normal direction. It permits the analysis of the eddy currents inducted by the normal flux. For excitations with the flux 4e-6 Wb and 35e-6WB the anisotropic i conductivity it was applied. Zero conductivity it in the tangential direction. It permits the analysis of the eddy currents inducted by the longitudinal flux. The same simulations were executed for isotropic materials (ThyssenKrupp Electrical Steel Company NO 0.35 mm).
Variation of B x along x in the area of the lamination (excitation 1and ) 1.8 1.6 B x [T] without eddy currents 1.4 with eddy currents 1. 1 0.8 with eddy currents without eddy currents 0.6 0.4 0. 0-0. x [mm] 80 85 90 95 100 105 110
Variation of B z along x in the area of the lamination (excitation 35e-6 Wb and 4e-6 Wb. 0.16 0.14 B z [T] without eddy currents 0.1 0.1 0.08 with eddy currents 0.06 with eddy currents 0.04 0.0 x [mm] 0 80 85 90 95 100 105 110
Variation of B x along z in the area of the lamination 0.0 0.18 Bz [T] 016 0.16 0.14 0.1 without eddy currents 0.10 with eddy currents 0.08 0.06 0.04 0.0 z [mm] 000 0.00-0.05 0.00 0.05 0.10 0.15 0.0 0.5 0.30 0.35 0.40
The repartition of the flux density (left) and eddy currents (right) in the overlapping area ωt=0 ωt=30
The repartition of the flux density (left) and eddy currents (right) in the overlapping area ωt=60
The repartition of the flux density vectors (left) and eddy currents (right) for the different time instants make it possible to analyze the comportment of the phenomena s in the overlapping area. All the flux passes there from one lamination to another. This results in a significant normal flux in addition to the longitudinal flux. Both cause the eddy currents in the x-y and y-z plans. The inducted eddy currents are sources of own fluxes. The result is the existence of the flux density components in the y direction and eddy currents in x direction. Before the overlapping in the lamination only the longitudinal flux exists. Inside the overlapping area longitudinal flux decreases and the normal increases or eddy currents will decrease in the y-z plane and increase in the x-y plane. This results in an excitation of the significant x component of the eddy currents and of the turning of the flux in y direction (Fig. 9). The simulations with the assumption of zero conductivity values in x direction those with zero conductivity in z direction (fig 9) confirms the presented conclusions. Eddy current is shifted in phase relative to the flux density of a little less than 90.
Conclusions Existence of both normal and longitudinal fluxes should be taken into account in the regions where the flux passes from one lamination to another. The eddy currents induced by the fluxes modify the field distribution in the structure and should be taken into account. The internal air-gaps greater than 0.1 mm have a powerful impact on the field distribution, the insulation between the laminations of 0.01mm have a negligible effect on the passing of the flux.
The normal flux passing of the one sheet to another does not change direction immediately after the entrance to the lamination, the transition is progressively. The main purpose of these 3D simulations was to prepare a mathematical model that makes possible calculation of eddy currents distribution in thin ferromagnetic sheet with given flux. The estimation of eddy currents distribution ib i and their influence on electromagnetic field distribution, for both the longitudinal and normal fluxes, have been done.
We remark that the normal and longitudinal fluxes induce eddy currents. The induced eddy currents modify the field distribution in the structure and should be taken into account. In the prepared model existence of both normal and longitudinal fluxes should be taken into the consideration. The analysis of the step-lap p structure (model ) permits to conclude that the internal air-gaps higher than 0.1mm have an important influence on the field distribution. On the contrary the isolation between the laminations of 0.01mm has a negligible effect and can be omitted in the mathematical model.
Very important fact was concluded, d the direction i of the normal flux from one sheet to another one does not change immediately after the entrance of the lamination, and the transition is progressive. Therefore is not possible to introduce a simple boundary condition stated that flux turns immediately after entering into sheet. To take into account the eddy currents the superposition of D and 1D problem is needed