Methods, systems, equipment, and media for simulating the hysteresis and loss characteristics of magnetic powder cores.

Abstract

This application belongs to the technical field of hysteresis modeling of ferromagnetic materials, and discloses a method, system, device, and medium for simulating the hysteresis and loss characteristics of magnetic powder cores. The method includes: solving the eddy current field expression inside the metal spherical particles based on Maxwell's equations; multiplying the calculation result of the eddy current field expression by the total number of magnetic powder core particles to obtain the total eddy current loss inside the particles; dividing the surface of the magnetic powder core into equally spaced segments and setting the height of each segment. h Based on the eddy current path conduction coefficient, parameter optimization is performed using the nonlinear least squares method, combined with the skin effect, to determine the expression for the interparticle eddy current field. Based on the total eddy current loss within the particles and the interparticle loss of the metal sphere calculated according to the interparticle eddy current field expression, combined with the Preisach static hysteresis model, a simulation model of the hysteresis and loss characteristics of magnetic powder cores is constructed. This application can simulate the hysteresis and loss characteristics of granular magnetic powder core materials under high-frequency sinusoidal and non-sinusoidal excitation.

Description

Technical Field

[0001] This application relates to the field of hysteresis model technology for ferromagnetic materials, and in particular to a method, system, device and medium for simulating the hysteresis and loss characteristics of magnetic powder cores. Background Technology

[0002] Currently, methods for analyzing the hysteresis and loss characteristics of tape-type amorphous nanocrystals are relatively mature. However, research on how to accurately predict the hysteresis and loss characteristics of granular magnetic powder cores is still limited. The difficulty in accurately predicting the hysteresis and loss of magnetic powder cores lies in the unsolvable eddy current loss. This is because magnetic powder cores operate under high-frequency and complex excitation for extended periods. On the one hand, this results in extremely complex changes in the eddy current field within the metal spheres. On the other hand, during the pressure-bearing fabrication process of soft magnetic composite magnetic powder cores, the external insulating layer of the particles is prone to breakage, leading to direct contact between the metal particles and the formation of inter-particle eddy current conduction. These conditions make it difficult to predict and assess the eddy current field distribution and loss characteristics of granular magnetic powder cores. Summary of the Invention

[0003] To address the aforementioned issues, this application provides a method, system, device, and medium for simulating the hysteresis and loss characteristics of magnetic powder cores, thereby simulating the hysteresis and loss characteristics of granular magnetic powder core materials under high-frequency sinusoidal and non-sinusoidal excitation.

[0004] According to the first technical solution of this application, a method for simulating the hysteresis and loss characteristics of magnetic powder cores is provided, the method comprising:

[0005] The eddy current field expression inside the metal sphere particles is solved based on Maxwell's equations. The calculation result of the eddy current field expression is multiplied by the total number of magnetic powder particles to obtain the total eddy current loss inside the particles.

[0006] The surface of the magnetic powder core is divided into equally spaced segments, and the height of each segment is set. h and eddy path transmission coefficient k Based on the nonlinear least squares method h and k Parameter optimization was performed, and the expression for the interparticle eddy field was determined by combining the skin effect.

[0007] Based on the calculated total eddy current loss inside the particle and the inter-particle loss of the metal sphere calculated from the expression of the inter-particle eddy current field, a first simulation curve and a second simulation curve are constructed. The first simulation curve and the second simulation curve are compared with the measured loss curves of the magnetic powder core under high-frequency sinusoidal and multi-harmonic excitation to determine the simulation accuracy.

[0008] Furthermore, the expression for the vortex field is:

[0009] ;

[0010] In the formula, H cl This represents the vortex field inside the metal sphere. Represents electrical conductivity. B Represents magnetic flux density. t Represents time, a This represents the total number of particles inside the magnetic powder core. R Represents the radius of the metal sphere particles.

[0011] Furthermore, methods for solving the eddy current field expression inside metal spherical particles based on Maxwell's equations include:

[0012] Based on Maxwell's equations, the fundamental governing equations for the vortex field inside the metal sphere are obtained:

[0013] ;

[0014] In the formula, It is a second-order partial differential operator;

[0015] Expanding the second-order partial differential operator in the coordinate system into the sum of the second-order partial derivatives along the three coordinate axes, we obtain the three-dimensional Laplace operator expansion:

[0016] ;

[0017] In the formula, x , y and z These represent the x, y, and z directions of the coordinate system, respectively.

[0018] Based on the fundamental governing equations of the vortex field inside the metal sphere and the three-dimensional Laplace operator expansion, the governing equations of the vortex field inside the metal sphere in spherical coordinates are obtained as follows:

[0019] ;

[0020] In the formula, r Represents the distance from the center of the ball;

[0021] The control equations of the eddy current field inside the metal sphere in the spherical coordinate system are solved to obtain the expression of the eddy current field inside the metal sphere particles.

[0022] Furthermore, the expression for the interparticle eddy field is:

[0023] ;

[0024] In the formula, H cl-outer-i This represents the eddy current field between particles. w Represents angular frequency. mThe permeability of the magnetic powder core material ,d h Represents skin depth, R Represents the radius of the metal sphere particles. t Represents time, R ( h / d h ) represents the skin effect coefficient. Representative and d h Regarding the initial phase angle, Represents the output lag angle. i This represents the phase parameter related to the output hysteresis angle;

[0025] Furthermore, the expression for the interparticle eddy field is determined by the following method:

[0026] Based on setting the height of each magnetic powder core segment h The number of segments in the magnetic powder core is determined, and the surface of the magnetic powder core is divided into equally spaced segments.

[0027] By aligning the external alternating magnetic field strength parallel to the Z-axis, the fundamental governing equations for the interparticle eddy current field are established, as follows:

[0028] ;

[0029] In the formula, y i Represents the distance of each segment, - h / 2≤ y i ≤ h / 2, J i Represents current density. B i Representing the i The magnetic field strength of the segment, H i The amplitude of the external alternating magnetic field strength. w Represents angular frequency. t Represents time, Represents electrical conductivity. H i ( t () represents the intensity of the external alternating magnetic field. e Represents the natural constant;

[0030] Based on the fundamental governing equations of the interparticle eddy current field, the expression for the total magnetic field strength is obtained by solving:

[0031] ;

[0032] In the formula, H (y i , t () represents the total magnetic field strength. v Represents wave vector;

[0033] The expressions for skin depth and wave vector are determined as follows:

[0034] ;

[0035] In the formula, d h Represents skin depth, f Represents frequency;

[0036] Integrating both sides of the expression for the total magnetic field strength yields the average magnetic field strength. H avg :

[0037] ;

[0038] In the formula, R ( h / d h ) represents the skin effect coefficient. represent d h The relevant initial phase angle;

[0039] Based on the average magnetic field strength, the fundamental expression for the eddy current field of the segmented section is determined as follows:

[0040] ;

[0041] In the formula, H cl-outer-i This represents the eddy current field between particles. Representative and d h Regarding the initial phase angle, Represents the output lag angle. i This represents the phase parameter related to the output hysteresis angle;

[0042] R ( h / d h ), and i The expression is:

[0043] ;

[0044] In the formula, ch Represents the hyperbolic cosine function. sh Represents the hyperbolic sine function;

[0045] Introducing eddy current path transmission coefficient k The basic expression for the segmented cross-section eddy field is modified to obtain the expression for the interparticle eddy field.

[0046] Furthermore, based on the aforementioned expression for the interparticle eddy current field, the interparticle loss of the metal sphere is calculated using the following formula:

[0047] ;

[0048] In the formula, W cl-outer-i This represents the cross-sectional eddy current loss of each magnetic powder core segment. W cl-oute This represents the total interparticle eddy current loss of the magnetic ring sample. B Represents magnetic flux density. i Segment index representing the magnetic powder core. n This represents the number of segments in the magnetic powder core.

[0049] Furthermore, the standard for the simulation accuracy is:

[0050] Under high-frequency sinusoidal and multi-harmonic excitation, the error between simulated loss and measured loss is controlled within 5%.

[0051] According to the second technical solution of this application, a simulation system for the hysteresis and loss characteristics of magnetic powder cores is provided, the system comprising:

[0052] The particle internal loss calculation module is configured to solve the eddy current field expression inside the metal sphere particle based on Maxwell's equations, and multiply the calculation result of the eddy current field expression by the total number of magnetic powder core particles to obtain the total eddy current loss inside the particle.

[0053] The interparticle eddy current field calculation module is configured to divide the surface of the magnetic powder core into equally spaced segments, and set the height of each segment of the magnetic powder core. h and eddy path transmission coefficient k Based on the nonlinear least squares method h and k Parameter optimization was performed, and the expression for the interparticle eddy field was determined by combining the skin effect.

[0054] The simulation accuracy verification module is configured to construct a simulation model of the hysteresis and loss characteristics of the magnetic powder core based on the total eddy current loss inside the particle and the inter-particle loss of the metal sphere calculated according to the expression of the inter-particle eddy current field, combined with the Preisach static hysteresis model. Under high-frequency sinusoidal and multi-harmonic excitation, the predicted curve output by the simulation model of the hysteresis and loss characteristics of the magnetic powder core is compared with the measured curve to determine the simulation accuracy.

[0055] According to the third technical solution of this application, an electronic device is provided, the electronic device comprising: a memory for storing a computer program; and a processor for executing the computer program to implement the method described above.

[0056] According to the fourth technical solution of this application, a non-transitory computer-readable storage medium storing instructions is provided, which, when executed by a processor, performs the method described above.

[0057] The simulation methods, systems, equipment, and media for the hysteresis and loss characteristics of magnetic powder cores according to the various schemes in this application have at least the following technical effects:

[0058] This application, considering the particle properties of soft magnetic composite magnetic powder cores, divides the eddy current loss of the magnetic powder core into eddy current loss inside the metal spherical particles and eddy current loss between particles. For the eddy current loss inside the metal spherical particles, Maxwell's equations are established to solve for the eddy current field expression inside the metal spherical particles; for the eddy current loss between metal spherical particles, the surface of the magnetic powder core is divided into equally spaced segments, and the height of each segment is defined. h and eddy current transmission coefficient k Based on the principle of nonlinear least squares, for height h and eddy current transmission coefficient k Parameter optimization was performed, and the skin effect caused by high-frequency action was considered to derive an expression for the eddy current field between metal spherical particles. Finally, combined with the Preisach static hysteresis model, a simulation method for the hysteresis and loss characteristics of magnetic powder cores that takes into account the spatial eddy current distribution characteristics at the particle scale was established. Based on this model, the measured loss curves of magnetic powder cores under high-frequency sinusoidal and multi-harmonic excitation were compared with those of the magnetic powder cores. The loss calculation error was controlled at about 5%, verifying the effectiveness of the proposed dynamic hysteresis model.

[0059] Additional aspects and advantages of this application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of this application. Attached Figure Description

[0060] Figure 1 A flowchart illustrating a method for simulating the hysteresis and loss characteristics of a magnetic powder core, provided in an embodiment of this application;

[0061] Figure 2 A schematic diagram of the internal eddy current field of a metal sphere provided in an embodiment of this application;

[0062] Figure 3 A schematic diagram of an equally spaced segmented vortex field provided in an embodiment of this application;

[0063] Figure 4The following is a comparison chart of simulated and measured values ​​of hysteresis loss provided in the embodiments of this application; wherein, (a) is a prediction chart with a fundamental frequency of 500Hz, a harmonic order of 4, a harmonic content of 30%, and a magnetic flux density of 0.4T; (b) is a prediction chart with a fundamental frequency of 500Hz, a harmonic order of 8, a harmonic content of 30%, and a magnetic flux density of 0.6T.

[0064] Figure 5 This is a structural diagram of a magnetic powder core hysteresis and loss characteristic simulation system provided in an embodiment of this application. Detailed Implementation

[0065] To enable those skilled in the art to better understand the technical solution of this application, the application will be described in detail below with reference to the accompanying drawings and specific embodiments.

[0066] One aspect of this application provides a method for simulating the hysteresis and loss characteristics of magnetic powder cores. This method is applicable to calculating the hysteresis loss of granular magnetic powder cores under high-frequency sinusoidal and non-sinusoidal excitation. The method can be executed by a hysteresis loss measurement and calculation device for ferromagnetic materials. This device can be implemented in hardware and / or software and can be configured in a computer device. Please refer to... Figure 1 This is a flowchart of a method for simulating the hysteresis and loss characteristics of a magnetic powder core, provided in an embodiment of this application. The method includes the following steps S10 to S30.

[0067] S10: Solve the eddy current field expression inside the metal sphere particles based on Maxwell's equations, and multiply the calculation result of the eddy current field expression by the total number of magnetic powder core particles to obtain the total eddy current loss inside the particles.

[0068] In some embodiments, please refer to Figure 2 This is a schematic diagram of the eddy current field inside a metal sphere particle provided in an embodiment of this application. The derivation of the expression for the eddy current field inside the metal sphere particle is as follows:

[0069] Assume the particle size of the metal spheres is 2. R And it is located in a rectangular coordinate system XYZ, with an externally applied alternating magnetic field. H a(t) Parallel to the Z-axis, the vortex field generated inside the metal sphere particle is in a counterclockwise direction. A schematic diagram of the vortex field inside the particle is shown below. Figure 2 As shown.

[0070] Based on Maxwell's equations, the fundamental governing equations for the eddy current field inside a metal sphere can be obtained (the equation relating the eddy current field to the rate of change of magnetic induction, derived from Maxwell's equations), which describes the eddy current field inside the metal sphere. H clThe second-order partial differential equations that satisfy the equations reflect the relationship between the eddy current field and the changes in conductivity and magnetic induction intensity over time. The basic governing equations of the eddy current field inside the metal sphere are as shown in equation (1):

[0071] (1);

[0072] In equation (1), H cl This represents the vortex field inside the metal sphere. Represents electrical conductivity. B Magnetic flux density t Represents time.

[0073] Based on equation (1), the second-order partial differential operator Expanding the Laplace operator in a Cartesian coordinate system into the sum of the second-order partial derivatives along the three coordinate axes yields the three-dimensional Laplace operator expansion, which is used for subsequent transformation equations in coordinate systems adapted to the geometric features of the metal sphere, such as spherical coordinates. This three-dimensional Laplace operator expansion is shown in equation (2).

[0074] (2);

[0075] Combining equations (1) and (2), the governing equations for the vortex field inside the metal sphere in spherical coordinates can be obtained. By transforming the Laplace operator in rectangular coordinates to spherical coordinates, for example, by utilizing the transformation relationship between spherical and rectangular coordinates and combining the spherical symmetry geometry of the metal sphere, the equations can be simplified to obtain the equations only relating to the radial distance (i.e., the distance from the center of the sphere). r The second-order partial differential equation of ) is more in line with the geometric characteristics analysis of the vortex field distribution inside the metal sphere. The governing equation of the vortex field inside the metal sphere in this spherical coordinate system is as shown in equation (3):

[0076] (3);

[0077] In equation (3), r Represents the distance from the center of the ball. e This represents the total number of particles inside the magnetic powder core.

[0078] Solving equation (3) yields the expression for the vortex field inside the metal sphere particles:

[0079] (4);

[0080] The expression for the eddy current field inside the metal sphere can be directly used to calculate the magnitude of the eddy current field inside the metal sphere, and the parameters related to the eddy current field and the particle are established. e , R Quantitative relationship between (etc.) and the rate of change of magnetic induction intensity.

[0081] S20: Divide the surface of the magnetic powder core into equally spaced segments and set the height of each segment. h and eddy path transmission coefficient k Based on the nonlinear least squares method h and k Parameter optimization was performed, and the expression for the interparticle eddy field was determined by combining the skin effect.

[0082] In this embodiment, the eddy current loss between the metal spherical particles is addressed by dividing the surface of the magnetic powder core into equally spaced segments and defining the height of each segment. h and eddy current transmission coefficient k Based on the principle of nonlinear least squares, for height h and eddy current transmission coefficient k By optimizing the parameters and considering the skin effect caused by high-frequency action, an expression for the eddy current field applicable to the particles of a metal sphere is derived.

[0083] In some embodiments, please refer to Figure 3 , Figure 3 A schematic diagram of an equally spaced segmented vortex field provided in this embodiment of the application is shown. The derivation process of the expression for the interparticle vortex field is as follows:

[0084] like Figure 3 As shown, the height of each segment is... h It is divided into a total of n Section. Assume the intensity of the alternating magnetic field in the external environment at this time is... Hi(t) Parallel to the Z-axis, the total magnetic field strength is expressed as H ( y i , t Establish (5):

[0085] (5);

[0086] In equation (5), y i Indicates the distance of each segment, - h / 2≤ y i ≤ h / 2. J i Represents current density. B i This represents the magnetic field strength of that segment. H i This represents the amplitude of the external alternating magnetic field strength. w ω is the angular frequency.

[0087] According to equation (5), the expression for the total magnetic field strength can be obtained:

[0088] (6);

[0089] in, H i =hB i / √2μδ h , m The permeability of the magnetic powder core material ,d h and v The physical meanings of are skin depth and wave vector, respectively, as expressed by equation (7):

[0090] (7);

[0091] Integrating both sides of equation (6) yields the average magnetic field strength for that segment. H avg :

[0092] (8);

[0093] Therefore, the vortex field of this section H cl-outer-i The basic expression is:

[0094] (9);

[0095] In the formula, H cl-outer-i This represents the eddy current field between particles. Representative and d h Regarding the initial phase angle, Represents the output lag angle. i This represents the phase parameter related to the output hysteresis angle.

[0096] R ( h / d h ), t and i The expression is:

[0097] (10);

[0098] In the formula, ch Represents the hyperbolic cosine function. sh This represents the hyperbolic sine function.

[0099] Define the eddy path conduction coefficient, optimize equation (9), and finally determine the interparticle eddy field. H cl-outer-i The expression is given by equation (11):

[0100] (11);

[0101] Based on the nonlinear least squares method and the fast gradient descent algorithm, h、k Two parameters are iteratively optimized to determine the parameters in equation (11). h、k The possible values ​​of the two parameters.

[0102] In some embodiments, based on a determined interparticle eddy field H cl-outer-i The expression - Equation (11) is used to calculate the eddy current loss of each cross section using the following formula. W cl-outer-i Total interparticle eddy current loss of magnetic ring sample W cl-outer :

[0103] (12);

[0104] In the formula, B Represents magnetic flux density. i Segment index representing the magnetic powder core. n This represents the number of segments in the magnetic powder core.

[0105] S30: Based on the total eddy current loss inside the particles and the inter-particle loss of the metal sphere calculated according to the expression of the inter-particle eddy current field, combined with the Preisach static hysteresis model, a simulation model of the hysteresis and loss characteristics of the magnetic powder core is constructed. Under high-frequency sinusoidal and multi-harmonic excitation, the predicted curve output by the simulation model of the hysteresis and loss characteristics of the magnetic powder core is compared with the measured curve to determine the simulation accuracy.

[0106] It should be noted that the Preisach static hysteresis model is a well-known model in the field, and its function is to solve for static hysteresis loss. This embodiment introduces abnormal losses based on the Preisach static hysteresis model to construct a simulation model of the hysteresis and loss characteristics of the magnetic powder core. The abnormal losses include the total eddy current loss inside the particle and the inter-particle loss of the metal spheres calculated in steps S10 and S20. The abnormal losses are obtained by superimposing these two losses. The obtained abnormal losses are then superimposed with the static hysteresis loss output by the Preisach static hysteresis model to obtain the simulated value of the total loss.

[0107] In some embodiments, the experimental platform measurement data is input into the simulation model of the hysteresis and loss characteristics of the magnetic powder core to obtain... Figure 4 The comparison graph showing the simulated and measured values ​​of hysteresis loss demonstrates that this application can accurately simulate the dynamic hysteresis characteristics of granular magnetic powder cores under high-frequency sinusoidal and non-sinusoidal excitation. The calculation error remains around 5% when compared with the measured loss values.

[0108] Another aspect of this application provides a simulation system for the hysteresis and loss characteristics of magnetic powder cores, such as... Figure 5 The diagram shown is a structural diagram of a magnetic powder core hysteresis and loss characteristic simulation system provided in an embodiment of this application. The magnetic powder core hysteresis and loss characteristic simulation system includes:

[0109] The particle internal loss calculation module 501 is configured to solve the eddy current field expression inside the metal sphere particle based on Maxwell's equations, and multiply the calculation result of the eddy current field expression by the total number of magnetic powder core particles to obtain the total eddy current loss inside the particle.

[0110] The interparticle eddy current field calculation module 502 is configured to divide the surface of the magnetic powder core into equally spaced segments and set the height of each segment of the magnetic powder core. h and eddy path transmission coefficient k Based on the nonlinear least squares method h and k Parameter optimization was performed, and the expression for the interparticle eddy field was determined by combining the skin effect.

[0111] The simulation accuracy verification module 503 is configured to construct a simulation model of the hysteresis and loss characteristics of the magnetic powder core based on the total eddy current loss inside the particle and the inter-particle loss of the metal sphere calculated according to the inter-particle eddy current field expression, combined with the Preisach static hysteresis model. Under high-frequency sinusoidal and multi-harmonic excitation, the predicted curve output by the simulation model of the hysteresis and loss characteristics of the magnetic powder core is compared with the measured curve to determine the simulation accuracy.

[0112] It should be noted that the magnetic powder core hysteresis and loss characteristic simulation system provided in the above embodiments and the magnetic powder core hysteresis and loss characteristic simulation method provided in the foregoing embodiments belong to the same concept. The specific way in which each module and unit performs operations has been described in detail in the method embodiments, and will not be repeated here.

[0113] Another aspect of this application provides an electronic device, including: a controller; and a memory for storing one or more programs, which, when executed by the controller, perform the methods described in the various embodiments above.

[0114] Another aspect of this application provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the aforementioned method. This computer-readable storage medium may be included in the electronic device described in the above embodiments, or it may exist independently and not incorporated into the electronic device.

[0115] It should be noted that the computer-readable medium shown in the embodiments of this application can be a computer-readable signal medium or a computer-readable storage medium, or any combination of the two. A computer-readable storage medium can be, for example, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of a computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM), flash memory, optical fiber, portable compact disc read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof. In this application, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device. In this application, a computer-readable signal medium can include a data signal propagated in baseband or as part of a carrier wave, carrying a computer-readable computer program. The transmitted data signal can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. The computer-readable signal medium can also be any computer-readable medium other than a computer-readable storage medium, which can send, propagate, or transmit a program for use by or in connection with an instruction execution system, apparatus, or device. The computer program contained on the computer-readable medium can be transmitted using any suitable medium, including but not limited to wireless, wired, etc., or any suitable combination thereof.

[0116] Another aspect of this application provides a computer program product or computer program that includes computer instructions stored in a computer-readable storage medium. A processor of a computer device reads the computer instructions from the computer-readable storage medium and executes the computer instructions, causing the computer device to perform the methods provided in the various embodiments described above.

[0117] According to one aspect of the embodiments of this application, a computer system is also provided, including a Central Processing Unit (CPU), which can perform various appropriate actions and processes based on a program stored in read-only memory (ROM) or a program loaded from storage into random access memory (RAM), such as performing the methods described above. Various programs and data required for system operation are also stored in the RAM. The CPU, ROM, and RAM are interconnected via a bus. Input / output (I / O) interfaces are also connected to the bus.

[0118] For example, a computer system includes a Central Processing Unit (CPU), which can perform various appropriate actions and processes based on programs stored in read-only memory (ROM) or loaded from storage into random access memory (RAM), such as executing the methods described in the above embodiments. The RAM also stores various programs and data required for system operation. The CPU, ROM, and RAM are interconnected via a bus. Input / output (I / O) interfaces are also connected to the bus.

[0119] The following components are connected to the I / O interface: input components including keyboards, mice, etc.; output components including cathode ray tubes (CRTs), liquid crystal displays (LCDs), and speakers; storage components including hard drives; and communication components including network interface cards such as LAN (Local Area Network) cards and modems. The communication components perform communication processing via networks such as the Internet. Drives are also connected to the I / O interface as needed. Removable media, such as disks, optical discs, magneto-optical discs, semiconductor memories, etc., are installed on the drive as needed so that computer programs read from them can be installed into the storage components as required.

[0120] Specifically, according to embodiments of this application, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, embodiments of this application include a computer program product comprising a computer program carried on a computer-readable medium, the computer program including a computer program for performing the methods shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via a communication component, and / or installed from a removable medium. When the computer program is executed by a central processing unit (CPU), it performs various functions defined in the system of this application.

[0121] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of this application. Each block in a flowchart or block diagram may represent a module, segment, or portion of code, which contains one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. For example, two consecutively indicated blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in a block diagram or flowchart, and combinations of blocks in a block diagram or flowchart, may be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.

[0122] The module units described in the embodiments of this application can be implemented in software or hardware, and the described units can also be located in a processor. The names of these units do not necessarily limit the specific unit itself.

[0123] The above embodiments are only used to illustrate this application and are not intended to limit this application. Those skilled in the art can make various changes and modifications without departing from the spirit and scope of this application. Therefore, all equivalent technical solutions also fall within the scope of this application, and the patent protection scope of this application should be defined by the claims.

Claims

1. A method for simulating the hysteresis and loss characteristics of magnetic powder cores, characterized in that, The method includes: The eddy current field expression inside the metal sphere particles is solved based on Maxwell's equations. The calculation result of the eddy current field expression is multiplied by the total number of magnetic powder particles to obtain the total eddy current loss inside the particles. The surface of the magnetic powder core is divided into equally spaced segments, and the height of each segment is set. h and eddy path transmission coefficient k Based on the nonlinear least squares method h and k Parameter optimization was performed, and the expression for the interparticle eddy field was determined by combining the skin effect. Based on the total eddy current loss inside the particles and the inter-particle loss of the metal sphere calculated according to the inter-particle eddy current field expression, combined with the Preisach static hysteresis model, a simulation model of the hysteresis and loss characteristics of the magnetic powder core is constructed. Under high-frequency sinusoidal and multi-harmonic excitation, the predicted curve output by the simulation model of the hysteresis and loss characteristics of the magnetic powder core is compared with the measured curve to determine the simulation accuracy. The expression for the interparticle eddy field is determined using the following method: Based on setting the height of each magnetic powder core segment h The number of segments in the magnetic powder core is determined, and the surface of the magnetic powder core is divided into equally spaced segments. By aligning the external alternating magnetic field strength parallel to the Z-axis, the fundamental governing equations for the interparticle eddy current field are established, as follows: ； In the formula, y i Represents the distance of each segment, - h / 2≤ y i ≤ h / 2, J i Represents current density, B i Representing the i The magnetic field strength of the segment, H i The amplitude of the external alternating magnetic field strength. w Represents angular frequency. t Represents time, Represents electrical conductivity. H i ( t () represents the intensity of the external alternating magnetic field. e Represents the natural constant; Based on the fundamental governing equations of the interparticle eddy current field, the expression for the total magnetic field strength is obtained by solving: ； In the formula, H ( y i , t () represents the total magnetic field strength. v Represents wave vector; The expressions for skin depth and wave vector are determined as follows: ； In the formula, δ h Represents skin depth, f Represents frequency; Integrating both sides of the expression for the total magnetic field strength yields the average magnetic field strength. H avg : ； In the formula, R ( h / δ h ) represents the skin effect coefficient. represent δ h The relevant initial phase angle; Based on the average magnetic field strength, the fundamental expression for the eddy current field of the segmented section is determined as follows: ； In the formula, H cl-outer-i This represents the eddy current field between particles. Represents the output lag angle. θ This represents the phase parameter related to the output hysteresis angle; R ( h / δ h ), and θ The expression is: ； In the formula, ch Represents the hyperbolic cosine function. sh Represents the hyperbolic sine function; Introducing eddy current path transmission coefficient k The basic expression for the segmented cross-section eddy field is modified to obtain the expression for the interparticle eddy field.

2. The method according to claim 1, characterized in that, The expression for the vortex field is: ； In the formula, H cl This represents the vortex field inside the metal sphere. Represents electrical conductivity. B Represents magnetic flux density. t Represents time, a This represents the total number of particles inside the magnetic powder core. R Represents the radius of the metal sphere particles.

3. The method according to claim 2, characterized in that, Methods for solving the eddy current field expression inside metal spherical particles based on Maxwell's equations include: Based on Maxwell's equations, the fundamental governing equations for the vortex field inside the metal sphere are obtained: ； In the formula, It is a second-order partial differential operator; Expanding the second-order partial differential operator in the coordinate system into the sum of the second-order partial derivatives along the three coordinate axes, we obtain the three-dimensional Laplace operator expansion: ； In the formula, x , y and z These represent the x, y, and z directions of the coordinate system, respectively. Based on the fundamental governing equations of the vortex field inside the metal sphere and the three-dimensional Laplace operator expansion, the governing equations of the vortex field inside the metal sphere in spherical coordinates are obtained as follows: ； In the formula, r Represents the distance from the center of the ball; The control equations of the eddy current field inside the metal sphere in the spherical coordinate system are solved to obtain the expression of the eddy current field inside the metal sphere particles.

4. The method according to claim 1, characterized in that, Based on the expression for the interparticle eddy current field, the interparticle loss of the metal sphere is calculated using the following formula: ； In the formula, W cl-outer-i This represents the cross-sectional eddy current loss of each magnetic powder core segment. W cl-oute This represents the total interparticle eddy current loss of the magnetic ring sample. B Represents magnetic flux density. i Segment index representing the magnetic powder core. n This represents the number of segments in the magnetic powder core.

5. The method according to claim 1, characterized in that, The standard for simulation accuracy is: Under high-frequency sinusoidal and multi-harmonic excitation, the error between simulated loss and measured loss is controlled within 5%.

6. A simulation system for the hysteresis and loss characteristics of a magnetic powder core, used to implement the method as described in any one of claims 1 to 5, characterized in that, The system includes: The particle internal loss calculation module is configured to solve the eddy current field expression inside the metal sphere particle based on Maxwell's equations, and multiply the calculation result of the eddy current field expression by the total number of magnetic powder core particles to obtain the total eddy current loss inside the particle. The interparticle eddy current field calculation module is configured to divide the surface of the magnetic powder core into equally spaced segments, and set the height of each segment of the magnetic powder core. h and eddy path transmission coefficient k Based on the nonlinear least squares method h and k Parameter optimization was performed, and the expression for the interparticle eddy field was determined by combining the skin effect. The simulation accuracy verification module is configured to construct a simulation model of the hysteresis and loss characteristics of the magnetic powder core based on the total eddy current loss inside the particle and the inter-particle loss of the metal sphere calculated according to the expression of the inter-particle eddy current field, combined with the Preisach static hysteresis model. Under high-frequency sinusoidal and multi-harmonic excitation, the predicted curve output by the simulation model of the hysteresis and loss characteristics of the magnetic powder core is compared with the measured curve to determine the simulation accuracy.

7. An electronic device, characterized in that, The electronic device includes: Memory, used to store computer programs; A processor for executing the computer program to implement the method as described in any one of claims 1 to 5.

8. A non-transitory computer-readable storage medium storing instructions, characterized in that, When the instructions are executed by the processor, the method according to any one of claims 1 to 5 is performed.

Images