Welding bead iron core loss calculation method and device based on finite element

Abstract

The invention provides a weld bead iron core loss calculation method and device based on a finite element, and the method comprises the steps: carrying out the simulation based on an iron core simulation model, and obtaining a model transient calculation result; obtaining a flux density curve corresponding to a non-welding area of the iron core according to a transient calculation result of the model; calculating iron core magnetic hysteresis loss and iron core eddy current loss of the whole iron core non-welding area according to the flux density curve and the bonding sample test BP curve; according to a model transient calculation result, calculating iron core hysteresis loss of the whole iron core welding area and iron core eddy current loss of the whole iron core welding area; and outputting the iron core magnetic hysteresis loss and the iron core eddy-current loss of the whole iron core non-welding area and the iron core magnetic hysteresis loss and the iron core eddy-current loss of the whole iron core welding area. According to the method and the device, the welding stator core loss can be quickly and accurately simulated and calculated, so that the calculation complexity is effectively reduced, the calculation efficiency is improved, the engineering design efficiency is improved, and the engineering design cost is reduced.

Description

Technical Field

[0001] The present application relates to the technical field of iron loss measurement, and in particular to a method and device for calculating weld core loss based on finite element. Background Art

[0002] In automotive drive motors, the core is a core functional component, and its performance directly affects the thermal efficiency and overall performance of the motor. In particular, the iron loss of the core (i.e., energy loss per unit mass) is crucial to improving motor efficiency. Accurately calculating the loss of the welded core has become a key link in improving motor efficiency. The existing calculation method mainly relies on the establishment of a stacking model and a short-circuit bar model, and evaluates the overall iron loss by calculating the eddy current loss of a single stack and estimating the hysteresis loss using a formula; for the welding area, it is simplified to a conductor with a certain conductivity, and only the eddy current loss generated is calculated. However, in practice, it is found that due to the extremely thin silicon steel stacks (about 0.2mm to 0.5mm), extremely small segmentation sizes are required for finite element analysis, resulting in a large number of computational grids, which not only consumes a lot of computing resources, but also significantly prolongs the calculation time. In some cases, the calculation cannot be completed due to resource limitations, which seriously affects the design efficiency and cost control. Summary of the invention

[0003] The purpose of the embodiments of the present application is to provide a method and device for calculating the core loss of a weld based on finite elements, which can quickly and accurately simulate and calculate the core loss of a welded stator, thereby effectively reducing the calculation complexity and improving the calculation efficiency, thereby improving the engineering design efficiency and helping to reduce the engineering design cost.

[0004] The first aspect of the present application provides a method for calculating weld core loss based on finite element method, comprising:

[0005] Establishing a core simulation model; wherein the core simulation model includes a core non-welding area, a core welding area, and an excitation coil winding;

[0006] Perform simulation based on the core simulation model to obtain transient calculation results of the model;

[0007] Obtain the BP curve of the bonding sample test;

[0008] Obtaining a magnetic flux density curve corresponding to the non-welded area of ​​the core according to the transient calculation result of the model; wherein the magnetic flux density curve is a curve of the magnetic flux density changing with time in each three-dimensional grid unit divided into the non-welded area of ​​the core;

[0009] Calculate the core hysteresis loss of the entire non-welded area of ​​the core according to the magnetic density curve and the BP curve of the bonding sample test; and calculate the core eddy current loss of the entire non-welded area of ​​the core according to the magnetic density curve and the BP curve of the bonding sample test;

[0010] Calculate the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area according to the transient calculation result of the model;

[0011] The core hysteresis loss of the entire core non-welding area, the core eddy current loss of the entire core non-welding area, the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area are output.

[0012] Furthermore, the establishment of the core simulation model includes:

[0013] Establish a sample core model that needs to be welded;

[0014] Acquire first configuration data of the excitation coil winding, second configuration data of the core non-welding area, and third configuration data of the core welding area;

[0015] According to the first configuration data, the second configuration data and the third configuration data, corresponding settings are respectively made for the core non-welding area, the core welding area and the excitation coil winding in the sample core model to obtain a core simulation model.

[0016] Furthermore, the first configuration data at least includes setting data of the current plane inlet and outlet of the excitation coil winding, setting conditions of the current plane inlet and outlet, and setting data of the current source excitation; wherein the condition of the current plane inlet and outlet is set as a symmetrical boundary condition;

[0017] The second configuration data at least includes a BH curve of a bonding sample test, a BP curve of a bonding sample test, a core lamination stacking coefficient, a core lamination conductivity, and core lamination eddy current setting data; wherein the core lamination conductivity is 0, and the core lamination eddy current setting data is that the core lamination does not allow eddy current;

[0018] The third configuration data at least includes a BH curve of a bonding sample test, a BP curve of a bonding sample test, electrical conductivity of a welding entity, and eddy current setting data of a welding entity; wherein the electrical conductivity of the welding entity is electrical conductivity of a silicon steel sheet, and the eddy current setting data of the welding entity is allowable eddy current of the welding entity.

[0019] Further, the calculation of the core hysteresis loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test includes:

[0020] Obtaining the amplitude and frequency of the fundamental wave large loop and each harmonic small loop according to the magnetic density curve;

[0021] According to the BP curve of the bonding sample test, the amplitude and frequency of the fundamental wave large loop and each harmonic small loop, the first core total loss relationship is determined; wherein the first core total loss relationship is a relationship between the core hysteresis loss and the core eddy current loss as a function of the magnetic flux density and the frequency;

[0022] Calculating the hysteresis loss of each three-dimensional grid unit according to the first core total loss relationship;

[0023] The core hysteresis loss of the entire core non-welding area is calculated according to the hysteresis loss of each three-dimensional grid unit.

[0024] Further, the calculation of the core eddy current loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test includes:

[0025] Performing FFT analysis on the magnetic density curve to obtain magnetic density amplitudes corresponding to different frequencies;

[0026] According to the BP curve of the bonding sample test and the magnetic flux amplitude corresponding to the different frequencies, the second core total loss relationship is determined; wherein the second core total loss relationship is a relationship between the core eddy current loss density and the magnetic flux and frequency;

[0027] Calculate the core eddy current loss of each three-dimensional grid unit according to the second core total loss relationship;

[0028] The eddy current loss of the core of the entire non-welded area of ​​the core is calculated according to the eddy current loss of the core of each three-dimensional grid unit.

[0029] Further, the calculating of the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area according to the transient calculation result of the model includes:

[0030] Calculate the core hysteresis loss of the entire core welding area according to the transient calculation result of the model and the BP curve of the bonding sample test;

[0031] Calculating the vector magnetic potential of each three-dimensional grid unit according to the transient calculation results of the model and the vector magnetic potential equation;

[0032] Calculating the electric field intensity of each three-dimensional grid unit according to the vector magnetic potential;

[0033] Calculating the current density of each three-dimensional grid unit according to the electric field strength;

[0034] Calculating the eddy current loss density within a unit volume according to the current density;

[0035] The core eddy current loss of the entire core welding area is calculated according to the eddy current loss density within the unit volume.

[0036] A second aspect of the present application provides a device for calculating the core loss of a weld based on finite elements, the device for calculating the core loss of a weld based on finite elements comprising:

[0037] A model building unit, used to build a core simulation model; wherein the core simulation model includes a core non-welding area, a core welding area and an excitation coil winding;

[0038] A simulation unit, used for performing simulation based on the core simulation model to obtain a transient calculation result of the model;

[0039] A first acquisition unit is used to acquire a BP curve of a bonding sample test;

[0040] A second acquisition unit is used to acquire a magnetic flux density curve corresponding to the non-welding area of ​​the core according to the transient calculation result of the model; wherein the magnetic flux density curve is a curve of the magnetic flux density changing with time in each three-dimensional grid unit divided by the non-welding area of ​​the core;

[0041] A first calculation unit is used to calculate the core hysteresis loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test;

[0042] A second calculation unit is used to calculate the core eddy current loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test;

[0043] A third calculation unit is used to calculate the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area according to the transient calculation result of the model;

[0044] The output unit is used to output the core hysteresis loss of the entire core non-welding area, the core eddy current loss of the entire core non-welding area, the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area.

[0045] Furthermore, the model building unit includes:

[0046] Establishing a subunit for establishing a sample core model to be welded;

[0047] A first acquisition subunit is used to acquire first configuration data of the excitation coil winding, second configuration data of the core non-welding area, and third configuration data of the core welding area;

[0048] The subunit is configured to perform corresponding configurations for the core non-welding area, the core welding area and the excitation coil winding in the sample core model according to the first configuration data, the second configuration data and the third configuration data to obtain a core simulation model.

[0049] Furthermore, the first configuration data at least includes setting data of the current plane inlet and outlet of the excitation coil winding, setting conditions of the current plane inlet and outlet, and setting data of the current source excitation; wherein the condition of the current plane inlet and outlet is set as a symmetrical boundary condition;

[0050] The second configuration data at least includes a BH curve of a bonding sample test, a BP curve of a bonding sample test, a core lamination stacking coefficient, a core lamination conductivity, and core lamination eddy current setting data; wherein the core lamination conductivity is 0, and the core lamination eddy current setting data is that the core lamination does not allow eddy current;

[0051] The third configuration data at least includes a BH curve of a bonding sample test, a BP curve of a bonding sample test, electrical conductivity of a welding entity, and eddy current setting data of a welding entity; wherein the electrical conductivity of the welding entity is electrical conductivity of a silicon steel sheet, and the eddy current setting data of the welding entity is allowable eddy current of the welding entity.

[0052] Furthermore, the first computing unit includes:

[0053] A second acquisition subunit is used to acquire the amplitude and frequency of the fundamental wave large loop and each harmonic small loop according to the magnetic density curve;

[0054] A first determination subunit is used to determine a first core total loss relationship according to the BP curve tested by the bonding sample, the amplitude and frequency of the fundamental wave large loop and each harmonic small loop; wherein the first core total loss relationship is a relationship between the core hysteresis loss and the core eddy current loss as the magnetic flux density and frequency change;

[0055] A first calculation subunit, configured to calculate the hysteresis loss of each three-dimensional grid unit according to the first core total loss relationship;

[0056] The first calculation subunit is further used to calculate the hysteresis loss of the entire core non-welding area according to the hysteresis loss of each three-dimensional grid unit.

[0057] Furthermore, the second computing unit includes:

[0058] An analysis subunit, used for performing FFT analysis on the magnetic density curve to obtain magnetic density amplitudes corresponding to different frequencies;

[0059] The second determination subunit is used to determine a second core total loss relationship according to the BP curve of the bonding sample test and the magnetic flux amplitude corresponding to the different frequencies; wherein the second core total loss relationship is a relationship between the core eddy current loss density and the magnetic flux and frequency;

[0060] A second calculation subunit, configured to calculate the core eddy current loss of each three-dimensional grid unit according to the second core total loss relationship;

[0061] The second calculation subunit is further used to calculate the core eddy current loss of the entire core non-welding area according to the core eddy current loss of each three-dimensional grid unit.

[0062] Furthermore, the third calculation unit is specifically used to calculate the core hysteresis loss of the entire core welding area according to the transient calculation results of the model and the BP curve of the bonding sample test; calculate the vector magnetic potential of each three-dimensional grid unit according to the transient calculation results of the model and the vector magnetic potential equation; calculate the electric field strength of each three-dimensional grid unit according to the vector magnetic potential; calculate the current density of each three-dimensional grid unit according to the electric field strength; calculate the eddy current loss density per unit volume according to the current density; calculate the core eddy current loss of the entire core welding area according to the eddy current loss density per unit volume.

[0063] A third aspect of the present application provides an electronic device, including a memory and a processor, wherein the memory is used to store a computer program, and the processor runs the computer program to enable the electronic device to execute the finite element-based weld core loss calculation method described in any one of the first aspect of the present application.

[0064] A fourth aspect of the present application provides a computer-readable storage medium storing computer program instructions. When the computer program instructions are read and executed by a processor, the finite element-based weld core loss calculation method described in any one of the first aspect of the present application is executed.

[0065] The beneficial effects of the present application are as follows: the method and device can quickly and accurately simulate and calculate the welding stator core loss, thereby effectively reducing the calculation complexity, improving the calculation efficiency, thereby improving the engineering design efficiency, and helping to reduce the engineering design cost. BRIEF DESCRIPTION OF THE DRAWINGS

[0066] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the embodiments of the present application will be briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present application and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0067] Figure 1 A schematic flow chart of a method for calculating weld core loss based on finite element provided in an embodiment of the present application;

[0068] Figure 2A schematic flow chart of another method for calculating weld core loss based on finite element provided in an embodiment of the present application;

[0069] Figure 3 A schematic diagram of a ring-shaped sample core model provided in an embodiment of the present application;

[0070] Figure 4 An example schematic diagram of the current inflow surface and outflow surface of an excitation coil winding provided in an embodiment of the present application;

[0071] Figure 5 An external circuit diagram of an excitation coil winding provided in an embodiment of the present application;

[0072] Figure 6 An example flow chart of a method for calculating welding core loss based on finite element provided in an embodiment of the present application;

[0073] Figure 7 A schematic diagram of the structure of a device for calculating weld core loss based on finite element provided in an embodiment of the present application;

[0074] Figure 8 A schematic structural diagram of another device for calculating weld core loss based on finite element method provided in an embodiment of the present application. DETAILED DESCRIPTION

[0075] The technical solutions in the embodiments of the present application will be described below in conjunction with the drawings in the embodiments of the present application.

[0076] It should be noted that similar reference numerals and letters represent similar items in the following drawings, so once an item is defined in one drawing, it does not need to be further defined and explained in subsequent drawings. At the same time, in the description of this application, the terms "first", "second", etc. are only used to distinguish the description and cannot be understood as indicating or implying relative importance.

[0077] Example 1

[0078] Please see Figure 1 , Figure 1 A schematic flow chart of a method for calculating the core loss of a weld based on finite elements provided in this embodiment. The method for calculating the core loss of a weld based on finite elements includes:

[0079] S101, establishing a core simulation model; wherein the core simulation model includes a core non-welding area, a core welding area, and an excitation coil winding.

[0080] S102, performing simulation based on the core simulation model to obtain transient calculation results of the model.

[0081] S103, obtaining a BP curve of the bonding sample test.

[0082] S104, obtaining a magnetic flux density curve corresponding to the non-welding area of ​​the core according to the transient calculation result of the model; wherein the magnetic flux density curve is a curve of the magnetic flux density changing with time in each three-dimensional grid unit divided into the non-welding area of ​​the core.

[0083] S105, calculating the core hysteresis loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test; and calculating the core eddy current loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test.

[0084] S106, calculating the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area according to the transient calculation result of the model.

[0085] S107, outputting the core hysteresis loss of the entire core non-welding area, the core eddy current loss of the entire core non-welding area, the core hysteresis loss of the entire core welding area, and the core eddy current loss of the entire core welding area.

[0086] In this embodiment, the execution subject of the method may be a computing device such as a computer or a server, and no limitation is made in this embodiment.

[0087] In this embodiment, the execution subject of the method may also be a smart device such as a smart phone, a tablet computer, etc., which is not limited in this embodiment.

[0088] It can be seen that the finite element-based weld core loss calculation method described in this embodiment can quickly and accurately simulate and calculate the welded stator core loss, thereby effectively reducing the calculation complexity and improving the calculation efficiency, thereby improving the engineering design efficiency and helping to reduce the engineering design cost.

[0089] Example 2

[0090] Please see Figure 2 , Figure 2 A schematic flow chart of a method for calculating the core loss of a weld based on finite elements provided in this embodiment. The method for calculating the core loss of a weld based on finite elements includes:

[0091] S201. Establish a sample core model that needs to be welded.

[0092] In this embodiment, the method can give priority to establishing the annular sample core model that needs to be welded, such as Figure 3 As shown. Among them, S1 is the non-welding area of ​​the core, which is made of laminated silicon steel sheets with a lamination coefficient of 0.96; W1 is the welding area of ​​the core, where no solder is used for welding. The silicon steel sheet base material is melted and then cooled again into a whole. Usually, the shape of the welding area is equivalent to a semicircle; C1 is the excitation coil winding.

[0093] S202, obtaining first configuration data of the excitation coil winding, second configuration data of the core non-welding area, and third configuration data of the core welding area.

[0094] In this embodiment, the first configuration data at least includes setting data of the current plane inlet and outlet of the excitation coil winding, setting conditions of the current plane inlet and outlet, and setting data of the current source excitation; wherein the condition of the current plane inlet and outlet is set as a symmetrical boundary condition;

[0095] The second configuration data at least includes a BH curve of a bonding sample test, a BP curve of a bonding sample test, a core lamination stacking coefficient, a core lamination conductivity, and core lamination eddy current setting data; wherein the core lamination conductivity is 0, and the core lamination eddy current setting data is that the core lamination does not allow eddy current;

[0096] The third configuration data at least includes the BH curve of the bonding sample test, the BP curve of the bonding sample test, the electrical conductivity of the welding entity and the eddy current setting data of the welding entity; wherein the electrical conductivity of the welding entity is the electrical conductivity of the silicon steel sheet, and the eddy current setting data of the welding entity is the allowable eddy current of the welding entity.

[0097] In this embodiment, for the related settings of the excitation coil winding C1, the inlet and outlet of the winding current surface can be set first, such as Figure 4 As shown ( Figure 4 Figure 1 shows an example schematic diagram of the current inflow and outflow surfaces of the excitation coil winding. Among them, face1 of the outer ring of C1 is the current inflow inlet, and face2 is the current outflow outlet; face3 of the inner ring of C1 is the current inflow inlet, and face4 is the current outflow outlet. At the same time, it is necessary to set symmetric boundary conditions for face1, face2, face3 and face4 to determine the current flow direction during simulation; by giving the C1 excitation coil winding Coil function (Coil function: the number of turns and resistance of the winding can be set), the inner and outer ring windings of the C1 excitation coil winding are connected into a closed circuit in the external circuit, and a current source excitation I1 is added to simulate the actual excitation coil. Figure 5 , Figure 5 An external circuit diagram of an excitation coil winding is shown.

[0098] In this embodiment, for the iron loss calculation method and related settings of the non-welded area S1 of the core, the BH and BP curves obtained by the "bonded ring sample" test can be first input into the non-welded area S1 of the core. The reason is that the "welded ring sample" is welded on the basis of the "bonded ring sample", which can ensure the rigor and accuracy of the experimental test data; secondly, in order to simulate the effect of core stacking, it is necessary to set the "core stacking coefficient"; at the same time, the iron loss calculation of the non-welded area of ​​the core is calculated according to the calculation idea of ​​the BP loss curve, so the conductivity of the core laminations needs to be set to 0, and the core laminations are set to not allow eddy currents to prevent repeated calculation of the eddy current loss of the core.

[0099] In this embodiment, for the iron loss calculation method and related settings of the core welding area W1, the BH and BP curves obtained by the "bonded ring sample" test can be first input into the core welding area W1. The reason is that the "welded ring sample" is welded on the basis of the "bonded ring sample", which can ensure the rigor and accuracy of the experimental test data; secondly, since the gap between the silicon steel sheets is extremely small, the silicon steel sheet mother sheet can be directly melted and then cooled into a whole, and the weld metal is still the silicon steel sheet mother material, so the electrical conductivity of the weld (core welding area) is set to be the electrical conductivity of the silicon steel sheet; and the weld (core welding area) is composed of the whole The core welding zone is composed of a silicon steel sheet base material, and the entire weld needs to be set to "allow eddy current" to calculate the eddy current loss of the core welding zone through the actual eddy current; since the material of the core welding zone is still the silicon steel sheet base material, this method assumes that the hysteresis loss characteristics of the core welding zone remain unchanged; therefore, in the iron loss calculation method of the core welding zone W1, the core hysteresis loss calculation method is still set to Method_h1, and the core eddy current loss calculation method is Method_e2, which is calculated according to the eddy current loss calculation method that actually considers the skin effect. The purpose of setting "allow eddy current" is to consider the influence of the skin effect.

[0100] S203, performing corresponding settings for the core non-welding area, the core welding area and the excitation coil winding in the sample core model according to the first configuration data, the second configuration data and the third configuration data, to obtain a core simulation model.

[0101] In this embodiment, the core simulation model includes a core non-welding area, a core welding area, and an excitation coil winding.

[0102] In this embodiment, the method divides the iron loss calculation method of the core non-welding area S1 into the core hysteresis loss calculation method Method_h1 and the core eddy current loss calculation method Method_e1; and divides the iron loss calculation method of the core welding area W1 into the core hysteresis loss calculation method Method_h1 (the same as the former) and the core eddy current loss calculation method Method_e2.

[0103] S204, performing simulation based on the core simulation model to obtain transient calculation results of the model.

[0104] In this embodiment, in the calculation process (steps S204 to S210 ) of the core hysteresis loss calculation method Method_h1 , the transient calculation result of the model is obtained first.

[0105] S205, obtaining a BP curve for bonding sample testing.

[0106] In this embodiment, this step is used to obtain the BP curve obtained by inputting the "bonded annular sample" test into the non-welded area S1 of the core. The horizontal axis of the BP curve of iron loss is the magnetic flux density, and the vertical axis is the hysteresis loss and eddy current loss.

[0107] S206. Obtain the magnetic flux density curve corresponding to the non-welding area of ​​the core according to the transient calculation result of the model.

[0108] In this embodiment, the magnetic flux density curve is a curve showing the change of magnetic flux density over time in each three-dimensional grid unit divided by the non-welding area of ​​the core.

[0109] In this embodiment, this step is used to extract the curve B(t) of the magnetic density changing with time in each three-dimensional grid unit Element N (N=1, 2, 3...N) divided in the non-welded area of ​​the core, and the time interval is a complete period T (T=1 / f; f is the fundamental frequency).

[0110] S207, obtaining the amplitude and frequency of the fundamental wave large loop and each harmonic small loop according to the magnetic density curve.

[0111] In this embodiment, the method is used to find the amplitudes B of the fundamental large loop and each harmonic small loop in B(t). k and frequency f k (k=1, 2, 3..., k=1 is a fundamental wave large loop, k=other numbers are harmonic small loops).

[0112] S208. Determine the total loss relationship of the first core according to the amplitude and frequency of the BP curve, fundamental wave large loop and each harmonic small loop tested on the bonding sample.

[0113] In this embodiment, the first total core loss relationship is a relationship between core hysteresis loss and core eddy current loss as a function of magnetic flux density and frequency.

[0114] In this embodiment, the method can fit the relationship between the core hysteresis loss and the core eddy current loss as the magnetic flux density B and the frequency f according to the BP curve of the bonded annular sample:

[0115]

[0116] Among them, Pfe(k) f k Total core loss at frequency; a(B k ) is f k Magnetic field amplitude B corresponding to the frequency k The core hysteresis loss coefficient under the condition of b(B k , f k ) is f k Magnetic field amplitude B corresponding to the frequency k The core eddy current loss coefficient under the condition of a(B k )·f k f k Core hysteresis loss P at frequency h ; f k Core eddy current loss P at frequency e .

[0117] S209, calculating the hysteresis loss of each three-dimensional grid unit according to the first core total loss relationship.

[0118] S210, calculating the core hysteresis loss of the entire core non-welding area according to the hysteresis loss of each three-dimensional grid unit.

[0119] In this embodiment, the method can first calculate the hysteresis loss of each grid unit, and then obtain the core hysteresis loss P of the entire core non-welding area by integration. h :

[0120]

[0121] Among them, N e is the number of grid cells in the three-dimensional grid division; NL is the number of hysteresis loops, k = 1 is the fundamental large loop, k = other is the harmonic small loop; B k is the magnetic flux density amplitude of the kth loop, f k is the kth frequency, k=1 is the fundamental frequency; V e is the volume of the grid cell.

[0122] S211. Perform FFT analysis on the magnetic density curve to obtain magnetic density amplitudes corresponding to different frequencies.

[0123] In this embodiment, the calculation of the core eddy current loss of the core non-welding area S1 (steps S211 to S214) can be performed separately after step S204 and step S206 to calculate the core eddy current loss of the core non-welding area S1.

[0124] In this embodiment, the method can perform FFT analysis on B(t) to obtain different frequencies f k The corresponding magnetic flux density amplitude B k(k=1,2,3…N; k=1 represents fundamental wave, k=other represents harmonic order).

[0125] S212. Determine the total loss relationship of the second core according to the BP curve of the bonding sample test and the magnetic flux density amplitude corresponding to different frequencies.

[0126] In this embodiment, the second total core loss relationship is a relationship between the core eddy current loss density and the magnetic flux density and frequency.

[0127] In this embodiment, the method can fit the BP curve of the bonded annular sample to obtain the relationship between the core hysteresis loss density and the core eddy current loss density as the magnetic flux density B and the frequency f:

[0128]

[0129] Among them, P fe(k) f k The total core loss density at frequency; a(B k ) is f k Magnetic field amplitude B corresponding to the frequency k The core hysteresis loss coefficient under the condition of b(B k , f k ) is f k Magnetic field amplitude B corresponding to the frequency k The core eddy current loss coefficient under the condition of a(B k )·f k f k Core hysteresis loss density at frequency; f k The core eddy current loss density at frequency.

[0130] S213, calculating the core eddy current loss of each three-dimensional grid unit according to the second core total loss relationship.

[0131] S214. Calculate the core eddy current loss of the entire core non-welding area according to the core eddy current loss of each three-dimensional grid unit.

[0132] In this embodiment, the method can first calculate the core eddy current loss of each grid unit, and then obtain the core eddy current loss P of the entire core non-welding area by integration. e :

[0133]

[0134] Among them, N e is the number of grid cells in the three-dimensional grid division; N is the highest harmonic order, k = 1 is the fundamental wave, k = other is the harmonic order; B k is the magnetic flux density amplitude of the kth harmonic, f kis the kth harmonic frequency, k=1 is the fundamental frequency; V e is the volume of the grid cell.

[0135] S215. Calculate the core hysteresis loss of the entire core welding area according to the model transient calculation results and the BP curve of the bonding sample test.

[0136] In this embodiment, this step is used to represent the calculation process of the core hysteresis loss in the core welding area, wherein the process is the same as the calculation process of the core hysteresis loss in the core non-welding area.

[0137] S216. Calculate the vector magnetic potential of each three-dimensional grid unit according to the model transient calculation results and the vector magnetic potential equation.

[0138] In this embodiment, during the calculation process of the calculation method Method_e2 of the core eddy current loss of the core welding area W1, a model transient calculation may be performed first.

[0139] In this embodiment, the transient calculation results are applied to the solution process of the vector magnetic potential equation (ie, the magnetic vector potential equation). Specifically, the transient calculation results include the distribution of the electromagnetic field, the change of the current density, and the like.

[0140] In this embodiment, the method can solve the vector magnetic potential A of each grid unit according to the vector magnetic potential equation; the magnetic vector potential equation is as follows:

[0141]

[0142] Where, μ is the core magnetic permeability; A is the vector magnetic potential of a grid unit; σ is the core electrical conductivity; J ext is the external current density.

[0143] S217. Calculate the electric field strength of each three-dimensional grid unit according to the vector magnetic potential.

[0144] In this embodiment, the method can calculate the electric field strength E of a certain grid unit according to the vector magnetic potential A obtained by solving:

[0145]

[0146] S218. Calculate the current density of each three-dimensional grid unit according to the electric field strength.

[0147] In this embodiment, the method can calculate the current density J of a certain grid unit according to the solved electric field strength E:

[0148] J = σE.

[0149] In this embodiment, considering the influence of skin effect, the current density J is a quantity related to the skin depth δ. In the finite element calculation, "allowed eddy current" is set, that is, the influence of skin effect is considered. The law shows that the closer to the surface of the silicon steel sheet, the greater the current density J; or the current density J shows a trend of gradually decaying from the surface of the silicon steel sheet to the inside. When it decays to 0.367 times the maximum eddy current density on the surface, the depth at this time is the skin depth δ, and its expression is:

[0150]

[0151] Among them, f is the frequency of magnetic field alternation; μ is the magnetic permeability of the iron core; σ is the electrical conductivity of the iron core.

[0152] S219. Calculate the eddy current loss density per unit volume based on the current density.

[0153] In this embodiment, the method can calculate the eddy current loss density of a certain grid unit according to the current density J obtained by solving, that is, the eddy current loss density p per unit volume. ve :

[0154]

[0155] Among them, p ve is the eddy current loss density per unit volume of a certain grid cell.

[0156] S220, calculating the core eddy current loss of the entire core welding area according to the eddy current loss density within the unit volume.

[0157] In this embodiment, the method can be based on the loss density p per unit volume. ve , the eddy current loss P in the core welding area is obtained by integration e :

[0158]

[0159] Among them, P e is the eddy current loss in the core welding area; V e is the volume of the grid unit, and Ne is the number of grid units.

[0160] S221, outputting the core hysteresis loss of the entire core non-welding area, the core eddy current loss of the entire core non-welding area, the core hysteresis loss of the entire core welding area, and the core eddy current loss of the entire core welding area.

[0161] In summary, the core loss value considering the weld was further obtained through finite element rapid calculation, and the effectiveness of this method was verified through experimental tests.

[0162] For the experimental test, the specific parameters of the experimental samples are as follows: the outer diameter of the annular sample silicon steel sheet is 60mm, the inner diameter is 30mm, the number of silicon steel sheets is 24, the thickness of a single sheet is 0.2mm, 2 symmetrical welds, argon arc welding, and the weld depth is about 5.0mm; the primary coil of the test coil has 100 turns and the secondary coil has 20 turns; the test instrument is TD8520, which is specially used to test the broadband AC magnetic measurement system of electrical steel sheets; due to the low magnetic flux density at the edge of the actual stator core weld (generally below 1.0T), when the frequency and magnetic flux density of the experimental test are low, the loss obtained by the test is low, and the loss data deviation is large due to the test interference, so the measured loss at higher frequency and higher magnetic flux density (the loss value is larger, and the error caused by interference is smaller) is taken here for comparative analysis to ensure the accuracy of the experimental test data.

[0163] Therefore. This method takes the simulation and experimental data of 500Hz / 0.9T and 1000Hz / 0.9T for comparative analysis. According to the existing technical method, the finite element segmentation simulation is performed. Since the thickness of the single piece is 0.2mm, the segmentation size of the non-welding area is set to 0.1mm, and the segmentation size of the welding area is set to 0.2mm. The total number of grids is about 24 million, as shown in Table 1. The existing computer memory resources are insufficient and the calculation cannot be run. This is also a major disadvantage of the existing technology; and according to this method, the finite element segmentation simulation is performed, the segmentation size of the non-welding area is set to 0.8mm, and the segmentation size of the welding area is set to 0.2mm. The total number of grids is about 490,000, and the calculation time is about 1.5 hours. The number of grids is reduced by nearly 50 times, which greatly improves the calculation efficiency.

[0164] By performing simulation calculations according to this method, we can finally obtain the results shown in Table 2. By comparing with the experimental test results, it is found that the iron loss value obtained by the calculation model of this method is basically consistent with the experimental test value, with an error fluctuating around 6%, which basically meets the requirements of engineering applications and verifies the effectiveness of this method.

[0165] Table 1 Comparison of finite element mesh numbers and simulation time used in the prior art and this embodiment

[0166]

[0167] Table 2 Comparison of simulated iron loss and experimental iron loss of this method

[0168]

[0169] Please see Figure 6 , Figure 6 An example flow chart of a method for calculating welding core loss based on finite element is shown.

[0170] In this embodiment, the execution subject of the method may be a computing device such as a computer or a server, and no limitation is made in this embodiment.

[0171] In this embodiment, the execution subject of the method may also be a smart device such as a smart phone, a tablet computer, etc., which is not limited in this embodiment.

[0172] It can be seen that the implementation of the finite element-based weld core loss calculation method described in this embodiment can improve the "non-welding area modeling" into a 3D modeling of an entire entity, and use the BP curve calculation method to calculate the iron loss; at the same time, on the basis of the unchanged modeling method of the welding area, a conductor with the same conductivity is set, and the welding area is set as a magnetic conductive overlap, so as to calculate the hysteresis loss of the welding area. Based on this, this method can improve the efficiency of welding iron loss calculation while ensuring the accuracy of iron loss calculation, thereby helping to improve engineering design efficiency and help reduce engineering design costs.

[0173] Example 3

[0174] Please see Figure 7 , Figure 7 The schematic diagram of the structure of a device for calculating the core loss of a weld based on finite element provided in this embodiment. Figure 7 As shown, the finite element-based weld core loss calculation device includes:

[0175] The model building unit 310 is used to build a core simulation model; wherein the core simulation model includes a core non-welding area, a core welding area, and an excitation coil winding;

[0176] A simulation unit 320 is used to perform simulation based on the core simulation model to obtain a transient calculation result of the model;

[0177] A first acquisition unit 330 is used to acquire a BP curve of a bonding sample test;

[0178] The second acquisition unit 340 is used to acquire the magnetic flux density curve corresponding to the non-welding area of ​​the core according to the transient calculation result of the model; wherein the magnetic flux density curve is a curve of the magnetic flux density changing with time in each three-dimensional grid unit divided into the non-welding area of ​​the core;

[0179] The first calculation unit 350 is used to calculate the core hysteresis loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test;

[0180] The second calculation unit 360 is used to calculate the eddy current loss of the core in the non-welded area of ​​the entire core according to the magnetic density curve and the BP curve of the bonding sample test;

[0181] The third calculation unit 370 is used to calculate the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area according to the transient calculation result of the model;

[0182] The output unit 380 is used to output the core hysteresis loss of the entire core non-welding area, the core eddy current loss of the entire core non-welding area, the core hysteresis loss of the entire core welding area, and the core eddy current loss of the entire core welding area.

[0183] In this embodiment, the explanation of the weld core loss calculation device based on finite element method can refer to the description in Embodiment 1 or Embodiment 2, which will not be further elaborated in this embodiment.

[0184] It can be seen that the finite element-based weld core loss calculation device described in this embodiment can accurately calculate the welding core loss, effectively reduce the calculation complexity, improve the calculation efficiency, improve the engineering design efficiency, and reduce the engineering design cost.

[0185] Example 4

[0186] Please see Figure 8 , Figure 8 The schematic diagram of the structure of a device for calculating the core loss of a weld based on finite element provided in this embodiment. Figure 8 As shown, the finite element-based weld core loss calculation device includes:

[0187] The model building unit 310 is used to build a core simulation model; wherein the core simulation model includes a core non-welding area, a core welding area, and an excitation coil winding;

[0188] A simulation unit 320 is used to perform simulation based on the core simulation model to obtain a transient calculation result of the model;

[0189] A first acquisition unit 330 is used to acquire a BP curve of a bonding sample test;

[0190] The second acquisition unit 340 is used to acquire the magnetic flux density curve corresponding to the non-welding area of ​​the core according to the transient calculation result of the model; wherein the magnetic flux density curve is a curve of the magnetic flux density changing with time in each three-dimensional grid unit divided into the non-welding area of ​​the core;

[0191] The first calculation unit 350 is used to calculate the core hysteresis loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test;

[0192] The second calculation unit 360 is used to calculate the eddy current loss of the core in the non-welded area of ​​the entire core according to the magnetic density curve and the BP curve of the bonding sample test;

[0193] The third calculation unit 370 is used to calculate the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area according to the transient calculation result of the model;

[0194] The output unit 380 is used to output the core hysteresis loss of the entire core non-welding area, the core eddy current loss of the entire core non-welding area, the core hysteresis loss of the entire core welding area, and the core eddy current loss of the entire core welding area.

[0195] As an optional implementation, the model building unit 310 includes:

[0196] Establishing subunit 311, used for establishing a sample core model to be welded;

[0197] A first acquisition subunit 312 is used to acquire first configuration data of the excitation coil winding, second configuration data of the core non-welding area, and third configuration data of the core welding area;

[0198] The setting subunit 313 is used to perform corresponding settings for the core non-welding area, the core welding area and the excitation coil winding in the sample core model according to the first configuration data, the second configuration data and the third configuration data to obtain a core simulation model.

[0199] In this embodiment, the first configuration data at least includes setting data of the current plane inlet and outlet of the excitation coil winding, setting conditions of the current plane inlet and outlet, and setting data of the current source excitation; wherein the condition of the current plane inlet and outlet is set as a symmetrical boundary condition;

[0200] The second configuration data at least includes a BH curve of a bonding sample test, a BP curve of a bonding sample test, a core lamination stacking coefficient, a core lamination conductivity, and core lamination eddy current setting data; wherein the core lamination conductivity is 0, and the core lamination eddy current setting data is that the core lamination does not allow eddy current;

[0201] The third configuration data at least includes the BH curve of the bonding sample test, the BP curve of the bonding sample test, the electrical conductivity of the welding entity and the eddy current setting data of the welding entity; wherein the electrical conductivity of the welding entity is the electrical conductivity of the silicon steel sheet, and the eddy current setting data of the welding entity is the allowable eddy current of the welding entity.

[0202] As an optional implementation, the first calculation unit 350 includes:

[0203] The second acquisition subunit 351 is used to acquire the amplitude and frequency of the fundamental wave large loop and each harmonic small loop according to the magnetic density curve;

[0204] The first determination subunit 352 is used to determine the first core total loss relationship according to the amplitude and frequency of the BP curve, fundamental large loop and each harmonic small loop tested by the bonding sample; wherein the first core total loss relationship is the relationship between the core hysteresis loss and the core eddy current loss as a function of the magnetic flux density and the frequency;

[0205] A first calculation subunit 353, configured to calculate the hysteresis loss of each three-dimensional grid unit according to a first core total loss relationship;

[0206] The first calculation subunit 353 is further used to calculate the hysteresis loss of the entire core non-welding area according to the hysteresis loss of each three-dimensional grid unit.

[0207] As an optional implementation, the second calculation unit 360 includes:

[0208] The analysis subunit 361 is used to perform FFT analysis on the magnetic density curve to obtain the magnetic density amplitude corresponding to different frequencies;

[0209] The second determination subunit 362 is used to determine the second core total loss relationship according to the BP curve of the bonding sample test and the magnetic flux amplitude corresponding to different frequencies; wherein the second core total loss relationship is the relationship between the core eddy current loss density and the magnetic flux and frequency;

[0210] A second calculation subunit 363, configured to calculate the core eddy current loss of each three-dimensional grid unit according to a second core total loss relationship;

[0211] The second calculation subunit 363 is further used to calculate the core eddy current loss of the entire core non-welding area according to the core eddy current loss of each three-dimensional grid unit.

[0212] As an optional implementation, the third calculation unit 370 is specifically used to calculate the core hysteresis loss of the entire core welding area according to the model transient calculation results and the bonding sample test BP curve; calculate the vector magnetic potential of each three-dimensional grid unit according to the model transient calculation results and the vector magnetic potential equation; calculate the electric field strength of each three-dimensional grid unit according to the vector magnetic potential; calculate the current density of each three-dimensional grid unit according to the electric field strength; calculate the eddy current loss density per unit volume according to the current density; calculate the core eddy current loss of the entire core welding area according to the eddy current loss density per unit volume.

[0213] In this embodiment, the explanation of the weld core loss calculation device based on finite element method can refer to the description in Embodiment 1 or Embodiment 2, which will not be further elaborated in this embodiment.

[0214] It can be seen that the finite element-based weld core loss calculation device described in this embodiment can accurately calculate the welding core loss, effectively reduce the calculation complexity, improve the calculation efficiency, improve the engineering design efficiency, and reduce the engineering design cost.

[0215] An embodiment of the present application provides an electronic device, including a memory and a processor, wherein the memory is used to store a computer program, and the processor runs the computer program to enable the electronic device to execute the finite element-based weld core loss calculation method in Example 1 or Example 2 of the present application.

[0216] An embodiment of the present application provides a computer-readable storage medium storing computer program instructions. When the computer program instructions are read and executed by a processor, the finite element-based weld core loss calculation method in Embodiment 1 or Embodiment 2 of the present application is executed.

[0217] In several embodiments provided in the present application, it should be understood that the disclosed devices and methods can also be implemented in other ways. The device embodiments described above are merely schematic. For example, the flowcharts and block diagrams in the accompanying drawings show the possible architecture, functions and operations of the devices, methods and computer program products according to multiple embodiments of the present application. In this regard, each box in the flowchart or block diagram can represent a module, a program segment or a part of a code, and the module, a program segment or a part of a code contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the box can also occur in a different order from the order marked in the accompanying drawings. For example, two consecutive boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram and / or flowchart, and the combination of boxes in the block diagram and / or flowchart can be implemented with a dedicated hardware-based system that performs a specified function or action, or can be implemented with a combination of dedicated hardware and computer instructions.

[0218] In addition, the functional modules in the various embodiments of the present application may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.

[0219] If the functions are implemented in the form of software function modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application can be essentially or partly embodied in the form of a software product that contributes to the prior art. The computer software product is stored in a storage medium, including several instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in each embodiment of the present application. The aforementioned storage medium includes: various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0220] The above description is only an embodiment of the present application and is not intended to limit the scope of protection of the present application. For those skilled in the art, the present application may have various changes and variations. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present application should be included in the scope of protection of the present application. It should be noted that similar reference numerals and letters represent similar items in the following drawings, so once an item is defined in one drawing, it does not need to be further defined and explained in the subsequent drawings.

[0221] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any technician familiar with the technical field can easily think of changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

[0222] It should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the existence of other identical elements in the process, method, article or device including the elements.

Claims

1. A method for calculating weld core loss based on finite element, characterized in that: include: Establishing a core simulation model; wherein the core simulation model includes a core non-welding area, a core welding area, and an excitation coil winding; Perform simulation based on the core simulation model to obtain transient calculation results of the model; Obtain the BP curve of the bonding sample test; Obtaining a magnetic flux density curve corresponding to the non-welded area of ​​the core according to the transient calculation result of the model; wherein the magnetic flux density curve is a curve of the magnetic flux density changing with time in each three-dimensional grid unit divided into the non-welded area of ​​the core; Calculate the core hysteresis loss of the entire non-welded area of ​​the core according to the magnetic density curve and the BP curve of the bonding sample test; and calculate the core eddy current loss of the entire non-welded area of ​​the core according to the magnetic density curve and the BP curve of the bonding sample test; Calculate the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area according to the transient calculation result of the model; The core hysteresis loss of the entire core non-welding area, the core eddy current loss of the entire core non-welding area, the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area are output.

2. The method for calculating weld core loss based on finite element according to claim 1, characterized in that: The core simulation model is established, comprising: Establish a sample core model that needs to be welded; Acquire first configuration data of the excitation coil winding, second configuration data of the core non-welding area, and third configuration data of the core welding area; According to the first configuration data, the second configuration data and the third configuration data, corresponding settings are respectively made for the core non-welding area, the core welding area and the excitation coil winding in the sample core model to obtain a core simulation model.

3. The method for calculating weld core loss based on finite element according to claim 2, characterized in that: The first configuration data at least includes setting data of the current plane inlet and outlet of the excitation coil winding, setting conditions of the current plane inlet and outlet, and setting data of the current source excitation; wherein the condition of the current plane inlet and outlet is set as a symmetrical boundary condition; The second configuration data at least includes a BH curve of a bonding sample test, a BP curve of a bonding sample test, a core lamination stacking coefficient, a core lamination conductivity, and core lamination eddy current setting data; wherein the core lamination conductivity is 0, and the core lamination eddy current setting data is that the core lamination does not allow eddy current; The third configuration data at least includes a BH curve of a bonding sample test, a BP curve of a bonding sample test, electrical conductivity of a welding entity, and eddy current setting data of a welding entity; wherein the electrical conductivity of the welding entity is electrical conductivity of a silicon steel sheet, and the eddy current setting data of the welding entity is allowable eddy current of the welding entity.

4. The method for calculating weld core loss based on finite element according to claim 1, characterized in that: The calculating the core hysteresis loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test includes: Obtaining the amplitude and frequency of the fundamental wave large loop and each harmonic small loop according to the magnetic density curve; According to the BP curve of the bonding sample test, the amplitude and frequency of the fundamental wave large loop and each harmonic small loop, the first core total loss relationship is determined; wherein the first core total loss relationship is a relationship between the core hysteresis loss and the core eddy current loss as a function of the magnetic flux density and the frequency; Calculating the hysteresis loss of each three-dimensional grid unit according to the first core total loss relationship; The core hysteresis loss of the entire core non-welding area is calculated according to the hysteresis loss of each three-dimensional grid unit.

5. The method for calculating weld core loss based on finite element according to claim 1, characterized in that: The calculating the core eddy current loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test includes: Performing FFT analysis on the magnetic density curve to obtain magnetic density amplitudes corresponding to different frequencies; According to the BP curve of the bonding sample test and the magnetic flux amplitude corresponding to the different frequencies, the second core total loss relationship is determined; wherein the second core total loss relationship is a relationship between the core eddy current loss density and the magnetic flux and frequency; Calculate the core eddy current loss of each three-dimensional grid unit according to the second core total loss relationship; The eddy current loss of the core of the entire non-welded area of ​​the core is calculated according to the eddy current loss of the core of each three-dimensional grid unit.

6. The method for calculating weld core loss based on finite element according to claim 1, characterized in that: The calculating the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area according to the transient calculation result of the model includes: Calculate the core hysteresis loss of the entire core welding area according to the transient calculation result of the model and the BP curve of the bonding sample test; Calculating the vector magnetic potential of each three-dimensional grid unit according to the transient calculation results of the model and the vector magnetic potential equation; Calculating the electric field intensity of each three-dimensional grid unit according to the vector magnetic potential; Calculating the current density of each three-dimensional grid unit according to the electric field strength; Calculating the eddy current loss density within a unit volume according to the current density; The core eddy current loss of the entire core welding area is calculated according to the eddy current loss density within the unit volume.

7. A device for calculating weld core loss based on finite element, characterized in that: The finite element-based weld core loss calculation device comprises: A model building unit, used to build a core simulation model; wherein the core simulation model includes a core non-welding area, a core welding area and an excitation coil winding; A simulation unit, used for performing simulation based on the core simulation model to obtain a transient calculation result of the model; A first acquisition unit is used to acquire a BP curve of a bonding sample test; A second acquisition unit is used to acquire a magnetic flux density curve corresponding to the non-welding area of ​​the core according to the transient calculation result of the model; wherein the magnetic flux density curve is a curve of the magnetic flux density changing with time in each three-dimensional grid unit divided by the non-welding area of ​​the core; A first calculation unit is used to calculate the core hysteresis loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test; A second calculation unit is used to calculate the core eddy current loss of the entire core non-welding area according to the magnetic density curve and the BP curve of the bonding sample test; A third calculation unit is used to calculate the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area according to the transient calculation result of the model; The output unit is used to output the core hysteresis loss of the entire core non-welding area, the core eddy current loss of the entire core non-welding area, the core hysteresis loss of the entire core welding area and the core eddy current loss of the entire core welding area.

8. The device for calculating weld core loss based on finite element according to claim 7, characterized in that: The model building unit comprises: Establishing a subunit for establishing a sample core model to be welded; A first acquisition subunit is used to acquire first configuration data of the excitation coil winding, second configuration data of the core non-welding area, and third configuration data of the core welding area; The subunit is configured to perform corresponding configurations for the core non-welding area, the core welding area and the excitation coil winding in the sample core model according to the first configuration data, the second configuration data and the third configuration data to obtain a core simulation model.

9. An electronic device, characterized in that: The electronic device comprises a memory and a processor, the memory is used to store a computer program, and the processor runs the computer program to enable the electronic device to execute the weld core loss calculation method based on finite element according to any one of claims 1 to 6.

10. A readable storage medium, characterized in that: The readable storage medium stores computer program instructions, and when the computer program instructions are read and executed by a processor, the method for calculating weld core loss based on finite element according to any one of claims 1 to 6 is executed.

Images