A volume optimization design method of low-harmonic magnetically controlled reactor

Abstract

The present application relates to the field of reactor volume optimization design, in particular to a kind of low harmonic magnetic control reactor volume optimization design method.Solution includes: selecting low harmonic magnetic control reactor parameter variable, including core length, core height, core cross section side length, single side air gap length, winding total number of turns and single side winding layer number;According to the selected low harmonic magnetic control reactor parameter variable, the variable required in optimization process is calculated;Set constraint condition;Optimize the optimal parameter combination of the core volume minimum by improved genetic algorithm optimization solution.The target of the objective function of the present application is to minimize the outer volume of the core, while imposing a penalty on the inductance deviation, and various constraint conditions are set, including wire constraint, winding number of turns constraint, winding thickness constraint, winding insulation and shape constraint, maximum inductance value constraint and peak magnetic flux density constraint, to avoid performance failure caused by simply pursuing volume minimization.The present application is suitable for reactor volume optimization design.

Description

Technical Field

[0001] The present invention relates to the field of reactor volume optimization design, and in particular to a method for volume optimization design of a low-harmonic magnetically controlled reactor. Background Art

[0002] In the field of power system equipment, tuned-inductance series resonance test systems are widely used for AC withstand voltage testing of power transformers, power cables, and other capacitive equipment. Low-harmonic magnetically controlled reactors, with their ability to smoothly adjust their reactance, effectively achieve series resonance with the test equipment, making them a key component of this test system. Given the frequent movement required in practical applications, such equipment places stringent requirements on compactness. However, as a specialized electromagnetic device, low-harmonic magnetically controlled reactors lack a comprehensive design process. Currently, the determination of core window size and winding turns relies primarily on simulation results, lacking systematic design criteria. Existing design methods fail to prioritize reactor size minimization, resulting in widespread redundancy in core length design. This not only wastes material but also increases equipment cost and mobility. Furthermore, the optimization design of low-harmonic magnetically controlled reactors involves a mixture of continuous and discrete variables. This complex nature further complicates achieving minimum design size, necessitating the development of efficient optimization methods to overcome existing technical bottlenecks. Summary of the Invention

[0003] The purpose of the present invention is to overcome the shortcomings of the prior art and provide a volume optimization design method for a low-harmonic magnetic-controlled reactor, thereby realizing a design with minimized volume of the low-harmonic magnetic-controlled reactor.

[0004] The present invention adopts the following technical solutions to achieve the above-mentioned purpose. The present invention provides a volume optimization design method for a low-harmonic magnetically controlled reactor, comprising:

[0005] S1. Select the parameter variables of the low harmonic magnetic controlled reactor, including the core length, core height, core cross-section side length, single-side air gap length, total number of winding turns, and number of single-side winding layers;

[0006] S2. Calculate the variables required in the optimization process based on the selected low harmonic magnetically controlled reactor parameter variables;

[0007] S3. Set the objective function to minimize the outer volume of the core while imposing a penalty on inductance deviation.

[0008] S4. Set constraints;

[0009] S5. Optimize and solve the optimal parameter combination that minimizes the core volume through an improved genetic algorithm.

[0010] Furthermore, step S2 specifically includes:

[0011] Calculate the maximum inductance;

[0012] ;

[0013] Where, represents the total magnetic circuit length, Indicates the core length, Indicates the core height;

[0014] ;

[0015] Where, represents the cross-sectional area of ​​the core, Indicates the cross-sectional side length of the core;

[0016] ;

[0017] Where, Indicates the magnetic resistance of the core part at maximum inductance; represents the air gap length, and the BH magnetization curve is fitted into a two-fold line model. It represents the slope of the linear section of the BH magnetization curve of the core material, that is, the unsaturated magnetic permeability of the core;

[0018] ;

[0019] Where, represents the magnetic resistance of the air gap part, Indicates the magnetic permeability at the air gap, that is, the magnetic permeability of air;

[0020] The maximum inductance is:

[0021] ;

[0022] Where, Indicates the maximum inductance, N indicates the total number of winding turns;

[0023] Calculate the minimum inductance;

[0024] ;

[0025] Where, Indicates the magnetic resistance of the core part at minimum inductance;

[0026] ;

[0027] Where, represents the magnetic resistance of the orthogonal core, and the BH magnetization curve is fitted into a two-fold line model. It represents the slope of the saturation section of the BH magnetization curve of the core material, that is, the saturation permeability of the core;

[0028] The minimum inductance is:

[0029] ;

[0030] The peak flux density in the maximum inductance state and the minimum inductance state is calculated;

[0031] The core flux and the core effective flux density in the maximum inductance state are:

[0032] ;

[0033] ;

[0034] In the formula, the core flux in the maximum inductance state is represented by, the core effective flux density in the maximum inductance state is represented by, I 0 the effective value of the winding current in the maximum inductance state is represented by;

[0035] The peak flux density in the maximum inductance state is:

[0036] ;

[0037] In the formula, the peak flux density in the maximum inductance state is represented by;

[0038] The core flux and the core effective flux density in the minimum inductance state are:

[0039] ;

[0040] ;

[0041] In the formula, the core flux in the minimum inductance state is represented by, the core effective flux density in the minimum inductance state is represented by, I 1 the effective value of the winding current in the minimum inductance state is represented by;

[0042] The peak flux density in the minimum inductance state is:

[0043] ;

[0044] In the formula, the peak flux density in the minimum inductance state is represented by;

[0045] The inductance variation range is calculated;

[0046] ;

[0047] Where, is the inductance variation range, 、 are the maximum inductance and the minimum inductance respectively.

[0048] Furthermore, the objective function is:

[0049] ;

[0050] ;

[0051] ;

[0052] Where, represents the objective function, represents the outer volume of the core, Indicates that a penalty is imposed on the inductance deviation. represents the target inductance.

[0053] Furthermore, step S4 specifically includes:

[0054] Wire Constraints:

[0055] Select the appropriate wire size based on the given current as follows:

[0056] ;

[0057] Where, Indicates the effective value of the reactor's rated voltage; f Indicates the rated frequency of the reactor, usually 50 Hz.

[0058] According to the calculated maximum current , look up the table to obtain the wire cross-sectional area, and select the corresponding wire specification according to the wire cross-sectional area;

[0059] Winding turns constraint:

[0060] Number of turns per layer n and wire diameter The relationship is: ;

[0061] The winding wrapping position cannot exceed 90% of the core height;

[0062] Winding thickness constraints:

[0063] ;

[0064] ;

[0065] Where, Indicates the distance between the outer side of the winding and the surface of the core. represents the margin coefficient, Indicates the winding thickness, Indicates the frame thickness, m indicates the number of winding layers on one side, Indicates the thickness of the insulation between layers;

[0066] Inter-winding insulation and shape constraints: ;

[0067] Maximum inductance constraint: ;

[0068] Peak flux density constraint:

[0069] ;

[0070] ;

[0071] Where, It is the dividing point of the two broken line model of the BH magnetization curve.

[0072] Furthermore, step S5 specifically includes:

[0073] Initialize the population:

[0074] For continuous variables 、 、 and , generates random values ​​within a defined range;

[0075] For discrete variables N and m, generate random integers within the defined integer range;

[0076] choose:

[0077] First, use the elite strategy to select the best-fit individuals from the current population as elite individuals;

[0078] Then, the tournament selection method is used to randomly select multiple individuals from the remaining individuals to compete, and the individual with the best fitness is selected as the parent;

[0079] cross:

[0080] For continuous variables 、 、 and , using simulated binary crossover, during which the genes of two parent individuals are exchanged and mixed to produce new offspring;

[0081] For discrete variables N and m, single-point crossover is used to randomly select the crossover point and exchange the corresponding parts of the parent individuals to produce new offspring;

[0082] Mutations:

[0083] In the new generation, for continuous variables 、 、 and ,Using polynomial mutation, continuous variables undergo small-scale random changes within their defined range with a certain set probability. During the mutation process, the gene value is adjusted according to the distribution index to generate new mutant individuals;

[0084] For discrete variables N and m, uniform mutation is used to achieve mutation by randomly selecting a new integer value within the allowed range to generate new mutant individuals;

[0085] Combine new mutant individuals with elite individuals to obtain a new generation of population;

[0086] termination:

[0087] An early stopping strategy is adopted. When the fitness improvement is less than the improvement threshold after the set number of iterations is reached, the optimization process is terminated in advance to obtain the optimal parameter combination that minimizes the core volume.

[0088] The beneficial effects of the present invention are:

[0089] The present invention sets an objective function, the goal of which is to minimize the outer volume of the iron core, while imposing a penalty on inductance deviation, and sets various constraints, including wire constraints, winding turn constraints, winding thickness constraints, inter-winding insulation and shape constraints, and maximum inductance value constraints, to ensure that the core functions of the inductor meet the standards, avoid performance failure caused by simply pursuing volume minimization, and prevent design parameters from violating physical laws or process limits.

[0090] In order to process different types of decision variables, the present invention adopts a hybrid coding strategy, that is, processing continuous variables and discrete variables at the same time, and using simulated binary crossover for continuous variables, and single-point crossover for discrete variables. The simulated binary crossover can maintain the excellent properties of continuous variables and has adaptive search capabilities. If the current optimal solution is small and the inductance meets the standard, the simulated binary crossover will conduct a fine search near it to avoid blind disturbances.

[0091] The present invention adopts polynomial mutation, controls the disturbance amplitude through the mutation exponent, and fine-tunes the parameters in small steps, which is suitable for refined optimization and can also automatically suppress the parameters from crossing the boundary after mutation. BRIEF DESCRIPTION OF THE DRAWINGS

[0092] Figure 1 This is a flow chart of a volume optimization design method for a low-harmonic magnetically controlled reactor provided by an embodiment of the present invention;

[0093] Figure 2This is a schematic diagram of the AC port-shaped core structure of a low-harmonic magnetically controlled reactor provided by an embodiment of the present invention;

[0094] Figure 3 This is a schematic diagram of the structure of the U-shaped iron core of the DC port of the low harmonic magnetically controlled reactor provided by an embodiment of the present invention;

[0095] Figure 4 1 is a schematic diagram of the winding structure of a low-harmonic magnetically controlled reactor provided by an embodiment of the present invention;

[0096] Figure 5 1 is a schematic diagram of a low harmonic magnetically controlled reactor model before optimization provided by an embodiment of the present invention;

[0097] Figure 6 Schematic diagram of an optimized low-harmonic magnetically controlled reactor model provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0098] To make the objectives, technical solutions and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention.

[0099] The present invention provides a method for optimizing the volume design of a low harmonic magnetic controlled reactor. Figure 1 As shown, specifically including:

[0100] S1. Select the parameter variables of the low harmonic magnetic controlled reactor, including the core length, core height, core cross-section side length, single-side air gap length, total number of winding turns, and number of single-side winding layers;

[0101] S2. Calculate the variables required in the optimization process based on the selected low harmonic magnetically controlled reactor parameter variables;

[0102] Calculate the maximum inductance;

[0103] ;

[0104] Where, represents the total magnetic circuit length, Indicates the core length, Indicates the core height;

[0105] ;

[0106] Where, represents the cross-sectional area of ​​the core, Indicates the cross-sectional side length of the core;

[0107] ;

[0108] Where, Indicates the magnetic resistance of the core part at maximum inductance; represents the air gap length; the BH magnetization curve is fitted into a two-fold line model, It represents the slope of the linear section of the BH magnetization curve of the core material, that is, the unsaturated magnetic permeability of the core;

[0109] ;

[0110] Where, represents the magnetic resistance of the air gap part, Indicates the magnetic permeability at the air gap, that is, the magnetic permeability of air;

[0111] The maximum inductance is:

[0112] ;

[0113] Where, Indicates the maximum inductance, N indicates the total number of winding turns;

[0114] Calculate the minimum inductance;

[0115] ;

[0116] Where, Indicates the magnetic resistance of the core part at minimum inductance;

[0117] ;

[0118] Where, represents the magnetic resistance of the orthogonal core, and the BH magnetization curve is fitted into a two-fold line model. It represents the slope of the saturation section of the BH magnetization curve of the core material, that is, the saturation permeability of the core;

[0119] Then the minimum inductance is:

[0120] ;

[0121] Calculate the peak flux density at maximum and minimum inductance;

[0122] The core magnetic flux and the effective magnetic flux density of the core under the maximum inductance state are:

[0123] ;

[0124] ;

[0125] Where, Indicates the core magnetic flux at maximum inductance, Indicates the effective magnetic flux density of the core under the maximum inductance state, I 0 Indicates the effective value of the winding current under the maximum inductance state;

[0126] The peak magnetic flux density in the maximum inductance state is:

[0127] ;

[0128] wherein, denotes the peak magnetic flux density in the maximum inductance state;

[0129] The core flux and the core effective flux density in the minimum inductance state are:

[0130] ;

[0131] ;

[0132] wherein, denotes the core flux in the minimum inductance state, denotes the core effective flux density in the minimum inductance state, I 1 denotes the effective value of the winding current in the minimum inductance state;

[0133] The peak magnetic flux density in the minimum inductance state is:

[0134] ;

[0135] wherein, denotes the peak magnetic flux density in the minimum inductance state;

[0136] The inductance variation range is calculated;

[0137] ;

[0138] wherein, is the inductance variation range, , are the maximum inductance and the minimum inductance, respectively.

[0139] S3, a target function is set, and the optimization target is to minimize the outer volume of the core while imposing a penalty on the inductance deviation;

[0140] ;

[0141] ;

[0142] ;

[0143] wherein, denotes the target function, denotes the outer volume of the core, denotes the penalty imposed on the inductance deviation, Inductance of target.

[0144] S4, set constraints;

[0145] Wire constraints:

[0146] According to the given current, select the corresponding wire specification, as follows:

[0147] ;

[0148] In the formula, Indicates the effective value of the reactor rated voltage; f Indicates the reactor rated frequency, generally 50 Hz.

[0149] According to the calculated maximum current , look up the wire cross-sectional area, and select the corresponding wire specification according to the wire cross-sectional area;

[0150] Winding turn constraint:

[0151] The relationship between the number of turns n of each layer and the wire diameter ; ;

[0152] The position of the winding cover cannot exceed 90% of the core height;

[0153] Winding thickness constraint:

[0154] ;

[0155] ;

[0156] In the formula, Indicates the distance from the outside of the winding to the surface of the core, Indicates the margin coefficient, Indicates the winding thickness, Indicates the frame thickness, m indicates the number of layers of single-sided winding, Indicates the interlayer insulation thickness, and the winding structure of the low harmonic adjustable reactor is shown in Figure 4 .

[0157] Winding insulation and shape constraint: ;

[0158] Maximum inductance constraint: ;

[0159] Peak magnetic flux density constraint:

[0160] ;

[0161] ;

[0162] Where, It is the dividing point of the two broken line model of the BH magnetization curve.

[0163] S5. Optimize and solve the optimal parameter combination that minimizes the core volume through an improved genetic algorithm.

[0164] Initialize the population:

[0165] For continuous variables 、 、 and , generates random values ​​within a defined range;

[0166] For discrete variables N and m, generate random integers within the defined integer range;

[0167] choose:

[0168] First, use the elite strategy to select the best-fit individuals from the current population as elite individuals;

[0169] Then, the tournament selection method is used to randomly select multiple individuals from the remaining individuals to compete, and the individual with the best fitness is selected as the parent;

[0170] The above process helps to maintain the diversity of the population and avoid falling into local optimality. The advantage of tournament selection is that it is simple and can better avoid the degeneration of elite individuals.

[0171] cross:

[0172] For continuous variables 、 、 and , using simulated binary crossover, during which the genes of two parent individuals are exchanged and mixed to produce new offspring;

[0173] For discrete variables N and m, single-point crossover is used to randomly select the crossover point and exchange the corresponding parts of the parent individuals to produce new offspring;

[0174] Mutations:

[0175] In the new generation, for continuous variables 、 、 and ,Using polynomial mutation, continuous variables undergo small-scale random changes within their defined range with a certain set probability. During the mutation process, the gene value is adjusted according to the distribution index to generate new mutant individuals;

[0176] For discrete variables N and m, uniform mutation is used to achieve mutation by randomly selecting a new integer value within the allowed range to generate new mutant individuals;

[0177] Combine new mutant individuals with elite individuals to obtain a new generation of population;

[0178] termination:

[0179] An early stopping strategy is adopted. When the fitness improvement is less than the improvement threshold after the set number of iterations, the optimization process is terminated in advance to avoid unnecessary resource computing consumption and obtain the optimal parameter combination that minimizes the core volume.

[0180] The following is an explanation with reference to specific design cases.

[0181] Design of low harmonic magnetic controlled reactor, rated voltage RMS , , , , the maximum inductance required is 0.1H, and the inductance variation range is .

[0182] The AC port-shaped iron core structure of the low harmonic magnetic control reactor is as follows Figure 2 As shown, the DC port shaped core structure is as follows Figure 3 The following takes the AC square core as an example.

[0183] The bounds on the decision variables are as follows:

[0184] AC port core length: ;

[0185] AC port-shaped core height: ;

[0186] Side length of the core cross section square: ;

[0187] Single-side air gap length: ;

[0188] Total number of turns of AC winding: N∈[200,1000];

[0189] Number of layers of single-sided AC winding: m∈[1,5];

[0190] Genetic algorithm parameters: population size 5000, number of iterations 200, elite ratio 0.1, number of randomly selected competing individuals in the tournament 3, crossover probability 0.8, and mutation probability 0.1.

[0191] Early stopping strategy parameters: improvement threshold 0.001, minimum number of iterations 60.

[0192] Table 1 below is a comparison of parameters before and after optimization:

[0193] Table 1 Comparison of parameters before and after optimization

[0194]

[0195] The low harmonic magnetic controlled reactor models before and after optimization are as follows: Figure 5 and Figure 6 shown.

[0196] The foregoing description is merely a preferred embodiment of the present invention. It should be understood that the present invention is not limited to the form disclosed herein and should not be construed as excluding other embodiments. Rather, the present invention can be used in various other combinations, modifications, and environments and can be modified within the scope of the concept described herein through the above teachings or techniques or knowledge in the relevant field. Modifications and variations made by those skilled in the art that do not depart from the spirit and scope of the present invention are intended to be protected by the appended claims.

Claims

1. A volume optimization design method for a low harmonic magnetic controlled reactor, characterized in that: include: S1. Select the parameter variables of the low harmonic magnetic controlled reactor, including the core length, core height, core cross-section side length, single-side air gap length, total number of winding turns, and number of single-side winding layers; S2. Calculate the variables required in the optimization process based on the selected low harmonic magnetically controlled reactor parameter variables; S3. Set the objective function to minimize the outer volume of the core while imposing a penalty on inductance deviation. S4. Set constraints; S5. Optimizing and solving the optimal parameter combination that minimizes the core volume through an improved genetic algorithm; Step S2 specifically includes: Calculate the maximum inductance; ； Where, represents the total magnetic circuit length, Indicates the core length, Indicates the core height; ； Where, represents the cross-sectional area of ​​the core, Indicates the cross-sectional side length of the core; ； Where, Indicates the magnetic resistance of the core part at maximum inductance, represents the air gap length, and the BH magnetization curve is fitted into a two-fold line model. It represents the slope of the linear section of the BH magnetization curve of the core material, that is, the unsaturated magnetic permeability of the core; ； Where, represents the magnetic resistance of the air gap part, Indicates the magnetic permeability at the air gap, that is, the magnetic permeability of air; The maximum inductance is: ； Where, Indicates the maximum inductance, N indicates the total number of winding turns; Calculate the minimum inductance; ； Where, Indicates the magnetic resistance of the core part at minimum inductance; ； Where, represents the magnetic resistance of the orthogonal core, and the BH magnetization curve is fitted into a two-fold line model. It represents the slope of the saturation section of the BH magnetization curve of the core material, that is, the saturation permeability of the core; Then the minimum inductance is: ； Calculate the peak flux density at maximum and minimum inductance; The core magnetic flux and the effective magnetic flux density of the core under the maximum inductance state are: ； ； Where, Indicates the core magnetic flux at maximum inductance, Indicates the effective magnetic flux density of the core under the maximum inductance state, Indicates the effective value of the winding current under the maximum inductance state; Then the peak flux density at maximum inductance is: ； Where, Indicates the peak magnetic flux density at maximum inductance; The core magnetic flux and the effective magnetic flux density of the core under the minimum inductance state are: ； ； Where, Indicates the core magnetic flux in the minimum inductance state, Indicates the effective magnetic flux density of the core under the minimum inductance state, I 1 Indicates the effective value of the winding current in the minimum inductance state; Then the peak flux density in the minimum inductance state is: ； Where, Indicates the peak magnetic flux density at minimum inductance; Calculate the inductance variation range; ； Where, is the inductance variation range, 、 are the maximum and minimum inductance values ​​respectively.

2. The volume optimization design method of a low harmonic magnetic controlled reactor according to claim 1, characterized in that: The objective function is: ； ； ； Where, represents the objective function, represents the outer volume of the core, Indicates that a penalty is imposed on the inductance deviation. represents the target inductance.

3. The volume optimization design method of a low harmonic magnetic controlled reactor according to claim 2, characterized in that: Step S4 specifically includes: Wire Constraints: Select the appropriate wire size based on the given current as follows: ； Where, Indicates the effective value of the reactor's rated voltage; f Indicates the rated frequency of the reactor; According to the calculated maximum current , look up the table to obtain the wire cross-sectional area, and select the corresponding wire specification according to the wire cross-sectional area; Winding turns constraint: Number of turns per layer n and wire diameter The relationship is: ; The winding wrapping position cannot exceed 90% of the core height; Winding thickness constraints: ； ； Where, Indicates the distance between the outer side of the winding and the surface of the core. represents the margin coefficient, Indicates the winding thickness, Indicates the frame thickness, m indicates the number of winding layers on one side, Indicates the thickness of the insulation between layers; Inter-winding insulation and shape constraints: ; Maximum inductance constraint: ; Peak flux density constraint: ； ； Where, It is the dividing point of the two broken line model of the BH magnetization curve.

4. The volume optimization design method of a low harmonic magnetic controlled reactor according to claim 3, characterized in that: Step S5 specifically includes: Initialize the population: For continuous variables 、 、 and , generates random values ​​within a defined range; For discrete variables N and m, generate random integers within the defined integer range; choose: First, use the elite strategy to select the best-fit individuals from the current population as elite individuals; Then, the tournament selection method is used to randomly select multiple individuals from the remaining individuals to compete, and the individual with the best fitness is selected as the parent; cross: For continuous variables 、 、 and , using simulated binary crossover, during which the genes of two parent individuals are exchanged and mixed to produce new offspring; For discrete variables N and m, single-point crossover is used to randomly select the crossover point and exchange the corresponding parts of the parent individuals to produce new offspring; Mutations: In the new generation, for continuous variables 、 、 and ,Using polynomial mutation, continuous variables undergo small-scale random changes within their defined range with a certain set probability. During the mutation process, the gene value is adjusted according to the distribution index to generate new mutant individuals; For discrete variables N and m, uniform mutation is used to achieve mutation by randomly selecting a new integer value within the allowed range to generate new mutant individuals; Combine new mutant individuals with elite individuals to obtain a new generation of population; termination: An early stopping strategy is adopted. When the fitness improvement is less than the improvement threshold after the set number of iterations is reached, the optimization process is terminated in advance to obtain the optimal parameter combination that minimizes the core volume.

Images