Turn-to-turn short circuit fault deduction method, system, equipment and medium

Abstract

The invention belongs to the technical field of power system faults, and discloses a turn-to-turn short circuit fault deduction method, system, equipment and medium, and the turn-to-turn short circuit fault deduction method comprises the steps: determining the critical strain value of the fracture of an insulating material based on the axial and radial tensile test data of a winding, and building a mechanical characteristic database of the dependent variable and the insulation fracture; establishing a Poisson boundary value equation set according to the magnetic vector potential and the current density, solving to obtain the magnetic induction intensity, establishing a circuit equation in a short circuit state based on the primary side voltage of the winding, and analyzing the linear relation between the deformation state of the winding and the electromagnetic load to obtain the maximum axial and radial deformation quantity of the winding; and based on the mechanical characteristic database, the turn-to-turn insulation state and the turn-to-turn short circuit condition are judged by combining the axial and radial maximum deformation quantity of the winding. According to the method, when the winding deformation is small, the occurrence condition of the turn-to-turn short circuit fault can be accurately deduced.

Description

Technical Field

[0001] The present invention relates to the technical field of power system faults, and in particular to a turn-to-turn short circuit fault deduction method, system, equipment and medium. Background Art

[0002] With the continuous expansion of the scale of power systems and the continuous growth of load demand, the transformer, as a key equipment in the power system, has a direct impact on the safe and stable operation of the entire power grid due to its operational reliability. Especially in the case of power grid failure, the transformer will be subjected to the impact of short-circuit current, which will cause the transformer winding to deform and have a significant cumulative effect; with the increase in the number of short-circuit impacts, the amount of winding deformation will gradually increase, and in severe cases, it may cause the transformer to have an inter-turn short-circuit fault, causing equipment damage and system operation interruption.

[0003] At present, the research on transformer winding deformation mainly focuses on the application of multi-physics field coupling technology, which simulates and analyzes the deformation process of windings under the impact of short-circuit current by establishing a coupling model of electromagnetic field, thermal field and mechanical field. This method can calculate the cumulative deformation of windings under multiple impacts, and judge the possible instability of transformers by the degree of winding deformation, providing theoretical support for transformer design and operation and maintenance.

[0004] However, the existing technology has obvious limitations: on the one hand, when the calculated winding deformation is small, it is not possible to directly determine whether the transformer will have an instability failure; on the other hand, the existing technology mainly focuses on the overall deformation characteristics of the winding and cannot effectively deduce the situation where the transformer has a turn-to-turn short-circuit failure, which makes the prevention and control of transformer turn-to-turn short-circuit failures face great challenges.

[0005] Therefore, how to provide a solution for transformer inter-turn short-circuit fault simulation based on material mechanical properties and multi-physics field simulation technology has become an urgent problem to be solved. Summary of the invention

[0006] The embodiments of the present invention provide a turn-to-turn short-circuit fault deduction method, system, device and medium to solve the problem in the prior art that when the calculated winding deformation is small, it is impossible to directly judge whether the transformer will have an instability fault and it is impossible to effectively deduce the situation where the transformer has a turn-to-turn short-circuit fault.

[0007] In order to have a basic understanding of some aspects of the disclosed embodiments, a brief summary is given below. This summary is not intended to be a general review, nor is it intended to identify key / important components or to delineate the scope of protection of these embodiments. Its only purpose is to present some concepts in a simple form as a preface to the detailed description that follows.

[0008] According to a first aspect of an embodiment of the present invention, a turn-to-turn short circuit fault deduction method is provided.

[0009] In one embodiment, the turn-to-turn short circuit fault deduction method includes:

[0010] Based on the axial tensile test data and radial tensile test data of the winding, the critical strain value of the insulation material rupture is determined, and according to the critical strain value of the insulation material rupture, a database of strain and mechanical properties of insulation rupture is established;

[0011] A group of Poisson boundary value equations is established according to the magnetic vector potential and current density, and the magnetic induction intensity is obtained by solving the group of Poisson boundary value equations; a circuit equation under a short-circuit state is established based on the primary side voltage of the winding, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed in combination with the magnetic induction intensity, and the maximum axial deformation and the maximum amplitude deformation of the winding are obtained;

[0012] Based on the mechanical property database of strain and insulation rupture, combined with the maximum axial deformation and maximum amplitude deformation of the winding, the inter-turn insulation status and inter-turn short circuit situation are judged.

[0013] In one embodiment, a Poisson boundary value equation group is established according to the magnetic vector potential and the current density, and the magnetic induction intensity is obtained by solving the Poisson boundary value equation group; a circuit equation under a short-circuit state is established based on the primary side voltage of the winding, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed in combination with the magnetic induction intensity, and the maximum axial deformation variable and the maximum amplitude deformation variable of the winding are obtained, including:

[0014] Based on the magnetic vector potential and current density, a Poisson boundary value equation group is established, the magnetic vector value of each node in the winding space is obtained by solving the Poisson boundary value equation group, and the magnetic induction intensity is calculated according to the magnetic vector value;

[0015] According to the primary side voltage of the winding, a circuit equation under the short-circuit state is established, and the maximum short-circuit current is calculated by solving the circuit equation under the short-circuit state;

[0016] The electromagnetic load is calculated based on the magnetic induction intensity and the circuit equation under the short-circuit state. The linear relationship between the deformation state of the winding and the electromagnetic load is analyzed by constructing the deformation variable equation to obtain the maximum axial deformation and the maximum radial deformation of the winding.

[0017] In one embodiment, a Poisson boundary value equation group is established based on the magnetic vector potential and the current density, the Poisson boundary value equation group is solved to obtain the magnetic vector value of each node in the winding space, and the magnetic induction intensity is calculated according to the magnetic vector value, including:

[0018] Based on the component relationship between magnetic vector potential and current density, a Poisson boundary value equation system including boundary conditions and magnetic resistivity is established.

[0019] The Poisson boundary value equations are solved by numerical iteration method to obtain the magnetic vector value of each node.

[0020] According to the magnetic vector value of each node, the corresponding magnetic induction intensity is calculated to obtain the leakage magnetic components in all directions in the winding space.

[0021] In one embodiment, the expression of the Poisson boundary value equations including boundary conditions and magnetic resistivity is:

[0022]

[0023] Wherein, Ω represents the solution domain of the winding space; z represents the axial coordinate in the cylindrical coordinate system; r represents the radial coordinate in the cylindrical coordinate system; θ represents the angular coordinate in the cylindrical coordinate system; t represents time; A θ represents magnetic vector potential; A represents magnetic vector potential; σ′ represents electrical conductivity; J sθ represents the current density J s The θ component of the θ component; Γ1 represents the first type of boundary; Γ2 represents the second type of boundary; H t represents its tangent vector; n represents the normal direction; t0 represents the initial time; ν′ represents the magnetic reluctance; f0(z,r) represents the functional relationship between the magnetic vector potential at time t0 and the coordinates z and r;

[0024] The expression of the magnetic vector value of each node is:

[0025] B=▽×A

[0026] In the formula, A represents the magnetic vector value of each node; ▽ represents the vector differential operator; B represents the isotropic leakage magnetic component in the winding space;

[0027] The expression of the isotropic leakage magnetic component in the winding space is:

[0028]

[0029] Wherein, B represents the isotropic leakage magnetic component of the winding space; x represents the abscissa in the rectangular coordinate system; y represents the ordinate in the rectangular coordinate system.

[0030] In one embodiment, according to the primary side voltage of the winding, a circuit equation under a short-circuit state is established, and by solving the circuit equation under the short-circuit state, the maximum short-circuit current is calculated, which includes:

[0031] According to the primary voltage of the winding, combined with the equivalent resistance and the equivalent leakage inductance, the circuit equation under the short-circuit state is established;

[0032] Solve the circuit equation under short-circuit state to obtain the steady-state component and transient component of the short-circuit current;

[0033] Based on the steady-state component of the short-circuit current, the maximum short-circuit current is calculated.

[0034] In one embodiment, the expression of the circuit equation in the short-circuit state is:

[0035]

[0036] Where, u1 represents the primary side voltage of the winding; U 1N represents the effective value of the primary voltage; w represents the frequency of the sine wave; t represents the time; α represents the initial phase of the voltage at the time of short circuit; R k Represents the equivalent resistance of a short-circuit circuit; i k Indicates short-circuit current; L k Indicates the equivalent leakage inductance of a short circuit;

[0037] The expression for the maximum short-circuit current is:

[0038]

[0039] In the formula, i k Indicates the maximum short-circuit current; I k Indicates the effective value of short-circuit current.

[0040] In one embodiment, based on the magnetic induction intensity and the circuit equation under the short-circuit state, the electromagnetic load is calculated, and by constructing a deformation variable solution equation, the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed, and the maximum axial deformation and the maximum amplitude deformation of the winding are obtained, including:

[0041] According to the magnetic induction intensity and the circuit equation under the short-circuit state, the electromagnetic load acting on the winding under the short-circuit state is calculated;

[0042] Based on the electromagnetic load acting on the winding under the short-circuit state, the deformation equation of the winding is established;

[0043] Based on the deformation equation of the winding, the equilibrium equation of the winding under state change is constructed, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed;

[0044] According to the analysis results of the linear relationship between the deformation state of the winding and the electromagnetic load, the deformation characteristics of the winding under the electromagnetic load are calculated, including the maximum axial deformation and the maximum radial deformation of the winding.

[0045] According to a second aspect of an embodiment of the present invention, a turn-to-turn short circuit fault deduction system is provided.

[0046] In one embodiment, the turn-to-turn short circuit fault deduction system includes:

[0047] A mechanical property database unit is used to determine the critical strain value of insulation material rupture based on the axial tensile test data and radial tensile test data of the winding, and to establish a mechanical property database of strain amount and insulation rupture according to the critical strain value of insulation material rupture;

[0048] The multi-physics field calculation unit is used to establish a Poisson boundary value equation group according to the magnetic vector potential and the current density, and obtain the magnetic induction intensity by solving the Poisson boundary value equation group; establish the circuit equation under the short-circuit state based on the primary side voltage of the winding, and analyze the linear relationship between the deformation state of the winding and the electromagnetic load in combination with the magnetic induction intensity, and obtain the maximum axial deformation and the maximum amplitude deformation of the winding;

[0049] The turn-to-turn short-circuit diagnostic unit is used to judge the turn-to-turn insulation status and turn-to-turn short-circuit situation based on the mechanical property database of strain and insulation rupture, combined with the maximum axial deformation and maximum amplitude deformation of the winding.

[0050] In one embodiment, the multi-physics computing unit includes:

[0051] The magnetic field distribution calculation module is used to establish a Poisson boundary value equation group based on the magnetic vector potential and current density, solve the Poisson boundary value equation group to obtain the magnetic vector value of each node in the winding space, and calculate the magnetic induction intensity based on the magnetic vector value;

[0052] A short-circuit current calculation module is used to establish a circuit equation under a short-circuit state according to the primary side voltage of the winding, and calculate the maximum short-circuit current by solving the circuit equation under the short-circuit state;

[0053] The winding deformation calculation module is used to calculate the electromagnetic load based on the magnetic induction intensity and the circuit equation under the short-circuit state, and to analyze the linear relationship between the deformation state of the winding and the electromagnetic load by constructing the deformation variable equation to obtain the maximum axial deformation and the maximum amplitude deformation of the winding.

[0054] In one embodiment, a Poisson boundary value equation group is established based on the magnetic vector potential and the current density, the Poisson boundary value equation group is solved to obtain the magnetic vector value of each node in the winding space, and the magnetic induction intensity is calculated according to the magnetic vector value, including:

[0055] Based on the component relationship between magnetic vector potential and current density, a Poisson boundary value equation system including boundary conditions and magnetic resistivity is established.

[0056] The Poisson boundary value equations are solved by numerical iteration method to obtain the magnetic vector value of each node.

[0057] According to the magnetic vector value of each node, the corresponding magnetic induction intensity is calculated to obtain the leakage magnetic components in all directions in the winding space.

[0058] In one embodiment, the expression of the Poisson boundary value equations including boundary conditions and magnetic resistivity is:

[0059]

[0060] Wherein, Ω represents the solution domain of the winding space; z represents the axial coordinate in the cylindrical coordinate system; r represents the radial coordinate in the cylindrical coordinate system; θ represents the angular coordinate in the cylindrical coordinate system; t represents time; A θ represents magnetic vector potential; A represents magnetic vector potential; σ′ represents electrical conductivity; J sθ represents the current density J s The θ component of the θ component; Γ1 represents the first type of boundary; Γ2 represents the second type of boundary; H t represents its tangent vector; n represents the normal direction; t0 represents the initial time; ν′ represents the magnetic reluctance; f0(z,r) represents the functional relationship between the magnetic vector potential at time t0 and the coordinates z and r;

[0061] The expression of the magnetic vector value of each node is:

[0062] B=▽×A

[0063] In the formula, A represents the magnetic vector value of each node; ▽ represents the vector differential operator; B represents the isotropic leakage magnetic component in the winding space;

[0064] The expression of the isotropic leakage magnetic component in the winding space is:

[0065]

[0066] Wherein, B represents the isotropic leakage magnetic component of the winding space; x represents the abscissa in the rectangular coordinate system; y represents the ordinate in the rectangular coordinate system.

[0067] In one embodiment, according to the primary side voltage of the winding, a circuit equation under a short-circuit state is established, and by solving the circuit equation under the short-circuit state, the maximum short-circuit current is calculated, which includes:

[0068] According to the primary voltage of the winding, combined with the equivalent resistance and the equivalent leakage inductance, the circuit equation under the short-circuit state is established;

[0069] Solve the circuit equation under short-circuit state to obtain the steady-state component and transient component of the short-circuit current;

[0070] Based on the steady-state component of the short-circuit current, the maximum short-circuit current is calculated.

[0071] In one embodiment, the expression of the circuit equation in the short-circuit state is:

[0072]

[0073] Where, u1 represents the primary side voltage of the winding; U1N represents the effective value of the primary voltage; w represents the frequency of the sine wave; t represents the time; α represents the initial phase of the voltage at the time of short circuit; R k Represents the equivalent resistance of a short-circuit circuit; i k Indicates short-circuit current; L k Indicates the equivalent leakage inductance of a short circuit;

[0074] The expression for the maximum short-circuit current is:

[0075]

[0076] In the formula, i k Indicates the maximum short-circuit current; I k Indicates the effective value of short-circuit current.

[0077] In one embodiment, based on the magnetic induction intensity and the circuit equation under the short-circuit state, the electromagnetic load is calculated, and by constructing a deformation variable solution equation, the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed, and the maximum axial deformation and the maximum amplitude deformation of the winding are obtained, including:

[0078] According to the magnetic induction intensity and the circuit equation under the short-circuit state, the electromagnetic load acting on the winding under the short-circuit state is calculated;

[0079] Based on the electromagnetic load acting on the winding under the short-circuit state, the deformation equation of the winding is established;

[0080] Based on the deformation equation of the winding, the equilibrium equation of the winding under state change is constructed, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed;

[0081] According to the analysis results of the linear relationship between the deformation state of the winding and the electromagnetic load, the deformation characteristics of the winding under the electromagnetic load are calculated, including the maximum axial deformation and the maximum radial deformation of the winding.

[0082] According to a third aspect of an embodiment of the present invention, a computer device is provided.

[0083] In one embodiment, the computer device includes a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps of the above method when executing the computer program.

[0084] According to a fourth aspect of embodiments of the present invention, a computer-readable storage medium is provided.

[0085] In one embodiment, the computer-readable storage medium stores a computer program, which implements the steps of the above method when executed by a processor.

[0086] The technical solution provided by the embodiment of the present invention may have the following beneficial effects:

[0087] The present invention proposes a deduction technology for turn-to-turn short-circuit faults based on material mechanical properties and multi-physical field simulation technology, which is of great significance for guiding the internal structure design of transformers. At the same time, the present invention breaks through the limitation of the prior art that instability faults cannot be judged under conditions of small deformation, and deeply analyzes the local deformation characteristics of the winding through multi-physical field simulation technology, effectively solving the technical problem that the prior art only focuses on the overall deformation but cannot accurately deduce turn-to-turn short-circuit faults, and provides a new technical approach for preventive maintenance and reliability design of transformers.

[0088] It is to be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0089] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the invention and, together with the description, serve to explain the principles of the invention.

[0090] Figure 1 is a flow chart of a method for deducing a turn-to-turn short-circuit fault according to an exemplary embodiment;

[0091] Figure 2 is a specific implementation diagram of a turn-to-turn short-circuit fault deduction method according to an exemplary embodiment;

[0092] Figure 3 is a schematic diagram of an insulating paper wrapped winding tensile test of a turn-to-turn short circuit fault deduction method according to an exemplary embodiment;

[0093] Figure 4 It is a transformer electromagnetic-structural field sequential coupling calculation flow chart of a turn-to-turn short circuit fault deduction method according to an exemplary embodiment;

[0094] Figure 5 is a structural block diagram of a turn-to-turn short-circuit fault deduction system according to an exemplary embodiment;

[0095] Figure 6 The diagram is a schematic structural diagram of a computer device according to an exemplary embodiment. DETAILED DESCRIPTION

[0096] The following description and accompanying drawings fully illustrate the specific embodiments of this article so that those skilled in the art can practice them. Parts and features of some embodiments may be included in or replace parts and features of other embodiments. The scope of the embodiments of this article includes the entire scope of the claims, as well as all available equivalents of the claims. Herein, the terms "first", "second", etc. are only used to distinguish one element from another, without requiring or implying any actual relationship or order between these elements. In fact, the first element can also be called the second element, and vice versa. Moreover, the terms "include", "comprise" or any other variant thereof are intended to cover non-exclusive inclusion, so that the structure, device or equipment including a series of elements includes not only those elements, but also other elements that are not explicitly listed, or also include elements inherent to such structure, device or equipment. In the absence of more restrictions, the elements defined by the sentence "including one..." do not exclude the existence of other identical elements in the structure, device or equipment including the elements. Each embodiment is described in a progressive manner herein, and each embodiment focuses on the differences from other embodiments, and the same and similar parts between the embodiments can be referred to each other.

[0097] The terms "longitudinal", "lateral", "upper", "lower", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside", etc. in this document indicate the orientation or position relationship based on the orientation or position relationship shown in the drawings, and are only for the convenience of describing this document and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operate in a specific orientation, and therefore cannot be understood as a limitation on the present invention. In the description of this document, unless otherwise specified and limited, the terms "installed", "connected", and "connected" should be understood in a broad sense, for example, it can be a mechanical connection or an electrical connection, it can also be the internal communication of two elements, it can be a direct connection, or it can be an indirect connection through an intermediate medium. For ordinary technicians in this field, the specific meanings of the above terms can be understood according to specific circumstances.

[0098] As used herein, the term "plurality" means two or more than two, unless otherwise specified.

[0099] In this document, the character " / " indicates that the preceding and following objects are in an "or" relationship. For example, A / B means: A or B.

[0100] In this article, the term "and / or" is a description of the association relationship between objects, indicating that three relationships may exist. For example, A and / or B means: A or B, or, A and B.

[0101] It should be understood that, although the various steps in the flow chart are displayed in sequence according to the indication of the arrows, these steps are not necessarily executed in sequence according to the order indicated by the arrows. Unless there is a clear description in this article, the execution of these steps is not strictly limited in order, and these steps can be executed in other orders. Moreover, at least a portion of the steps in the figure may include multiple sub-steps or multiple stages, and these sub-steps or stages are not necessarily executed at the same time, but can be executed at different times, and the execution order of these sub-steps or stages is not necessarily to be carried out in sequence, but can be executed in turn or alternately with other steps or at least a portion of the sub-steps or stages of other steps.

[0102] Each module in the device or system of the present application can be implemented in whole or in part by software, hardware, or a combination thereof. The above modules can be embedded in or independent of a processor in a computer device in the form of hardware, or can be stored in a memory in a computer device in the form of software, so that the processor can call and execute operations corresponding to the above modules.

[0103] In the absence of conflict, the embodiments of the present invention and the features of the embodiments may be combined with each other.

[0104] Figure 1-Figure 2 An embodiment of a turn-to-turn short-circuit fault deduction method of the present invention is shown.

[0105] In this optional embodiment, the turn-to-turn short circuit fault deduction method includes:

[0106] Step S101, based on the axial tensile test data and the radial tensile test data of the winding, determine the critical strain value of the insulation material rupture, and establish a strain amount and a mechanical property database of insulation rupture according to the critical strain value of the insulation material rupture;

[0107] Specifically, the present invention realizes the deduction of turn-to-turn short-circuit faults based on the mechanical properties of materials and multi-physics field simulation technology. Taking the transformer as an example, the specific implementation steps of the deduction are divided into three steps, including ① determination of the mechanical properties of the insulating paper wrapped winding: the test platform supports axial and radial stretching of the insulating paper wrapped winding by a specific distance; ② calculation of the winding deformation under multi-physics field coupling: based on the stress-strain curve of the winding, the elastic model quantity and other mechanical characteristic parameters, the finite element simulation calculation method is used to calculate the transient state of the winding deformation under short-circuit impact; ③ turn-to-turn short-circuit diagnosis: based on the calculation results of the winding deformation using multi-physics fields and the determination of the mechanical properties of the insulating paper wrapped winding, it is determined whether the insulating material on the winding surface is broken.

[0108] Specifically, turn-to-turn short-circuit faults caused by winding deformation are mostly caused by winding insulation rupture due to excessive deformation. Therefore, it is necessary to determine the deformation amount that causes the insulation material to rupture by judging the mechanical properties of the insulating paper wrapped winding. In this embodiment, axial and lateral tensile tests are performed on insulating paper wrapped windings of the same specifications as those of actual transformer manufacturing. The occurrence of obvious cracks in the insulating paper is used as the basis for judgment, and a mechanical property database of strain and insulation rupture is constructed to facilitate subsequent judgment of its mechanical properties.

[0109] Specifically, in the above embodiment, the axial tensile test data and the radial tensile test data of the winding are obtained by performing a pre-performed tensile test on the insulating paper wrapped winding, such as Figure 3 As shown in the figure, a is a clamping plate used to simulate the role of a pad; b is an insulating paper-wrapped winding; F a F is the axial pressure exerted by the outside world on the winding; r It is the radial pressure applied to the winding by the outside world. Through the tensile test, the insulating paper wrapped winding is stretched axially and radially for a specific distance, so as to measure the mechanical properties of the insulating paper wrapped winding.

[0110] Step S103, establishing a Poisson boundary value equation group according to the magnetic vector potential and the current density, and obtaining the magnetic induction intensity by solving the Poisson boundary value equation group; establishing a circuit equation under a short-circuit state based on the primary side voltage of the winding, and analyzing the linear relationship between the deformation state of the winding and the electromagnetic load in combination with the magnetic induction intensity, and obtaining the maximum axial deformation variable and the maximum amplitude deformation variable of the winding;

[0111] Step S105 : Based on the strain and the mechanical property database of insulation rupture, combined with the maximum axial deformation and the maximum radial deformation of the winding, the inter-turn insulation state and the inter-turn short circuit situation are determined.

[0112] Specifically, based on the mechanical property database measured in the experiment and the maximum axial and lateral deformation of the transformer winding calculated based on the multi-physics field simulation method, the diagnosis of the transformer inter-turn insulation status and the deduction of the inter-turn short-circuit situation can be achieved.

[0113] In this optional embodiment, a Poisson boundary value equation group is established according to the magnetic vector potential and the current density, and the magnetic induction intensity is obtained by solving the Poisson boundary value equation group; a circuit equation under a short-circuit state is established based on the primary side voltage of the winding, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed in combination with the magnetic induction intensity, and the maximum axial deformation variable and the maximum amplitude deformation variable of the winding are obtained, including:

[0114] Step S1031, establishing a Poisson boundary value equation group based on the magnetic vector potential and the current density, solving the Poisson boundary value equation group to obtain the magnetic vector value of each node in the winding space, and calculating the magnetic induction intensity according to the magnetic vector value;

[0115] Step S1033, establishing a circuit equation under a short-circuit state according to the primary side voltage of the winding, and calculating the maximum short-circuit current by solving the circuit equation under the short-circuit state;

[0116] Step S1035: Calculate the electromagnetic load based on the magnetic induction intensity and the circuit equation under the short-circuit state, and analyze the linear relationship between the deformation state of the winding and the electromagnetic load by constructing a deformation variable equation to obtain the maximum axial deformation and the maximum radial deformation of the winding.

[0117] Specifically, in the above embodiment, step ② of calculating the winding deformation under multi-physical field coupling includes three parts: (1) transformer magnetic field calculation; (2) transformer short-circuit current calculation; and (3) transformer winding deformation calculation.

[0118] In this optional embodiment, a Poisson boundary value equation group is established based on the magnetic vector potential and the current density, the Poisson boundary value equation group is solved to obtain the magnetic vector value of each node in the winding space, and the magnetic induction intensity is calculated according to the magnetic vector value, including:

[0119] Step S10311, based on the component relationship between the magnetic vector potential and the current density, establish a Poisson boundary value equation group including boundary conditions and magnetic resistivity;

[0120] Specifically, the transformer magnetic field calculation is based on the magnetic vector potential A θ As the intermediate quantity, a Poisson boundary value equation is constructed to obtain the magnetic field distribution of the transformer. The expression of the Poisson boundary value equation group including boundary conditions and magnetic reluctance is:

[0121]

[0122] Wherein, Ω represents the solution domain of the winding space; z represents the axial coordinate in the cylindrical coordinate system; r represents the radial coordinate in the cylindrical coordinate system; θ represents the angular coordinate in the cylindrical coordinate system; t represents time; A θ represents magnetic vector potential; A represents magnetic vector potential; σ′ represents electrical conductivity; J sθ represents the current density J s The θ component of the θ component; Γ1 represents the first type of boundary; Γ2 represents the second type of boundary; H t represents its tangent vector; n represents the normal direction; t0 represents the initial time; ν′ represents the magnetic reluctance, and v represents magnetic resistance; f0(z,r) represents the functional relationship between the magnetic vector potential and the coordinates z and r at time t0.

[0123] Step S10313, using a numerical iteration method to solve the Poisson boundary value equations to obtain the magnetic vector value of each node;

[0124] It should be noted that the Newton-Raphson iteration method is a numerical method for solving nonlinear equations, also known as the Newton iteration method. It approximates the roots of the equations through iteration. This is a prior art and will not be elaborated on here.

[0125] Specifically, in this embodiment, the Newton-Raphson iteration method is used for analytical calculation to obtain the magnetic vector value A of each node. The expression of the magnetic vector value of each node is:

[0126] B=▽×A

[0127] Where A represents the magnetic vector value of each node; ▽ represents the vector differential operator; and B represents the isotropic leakage magnetic component in the winding space.

[0128] Step S10315: Calculate the corresponding magnetic induction intensity according to the magnetic vector value of each node to obtain the leakage magnetic components in all directions in the winding space.

[0129] Specifically, the B value at the corresponding position in the winding space solution area can be obtained from the expression of the magnetic vector value of each node; thus, the isotropic leakage magnetic component of the winding space can be obtained. The expression of the isotropic leakage magnetic component of the winding space is:

[0130]

[0131] Wherein, B represents the isotropic leakage magnetic component of the winding space; x represents the abscissa in the rectangular coordinate system; y represents the ordinate in the rectangular coordinate system.

[0132] In this optional embodiment, according to the primary side voltage of the winding, a circuit equation under a short-circuit state is established, and by solving the circuit equation under a short-circuit state, the maximum short-circuit current is calculated, which includes:

[0133] Step S10331, establishing a circuit equation under a short-circuit state according to the primary side voltage of the winding, combined with an equivalent resistance and an equivalent leakage inductance;

[0134] Specifically, the transformer short-circuit voltage calculation, α is the initial voltage phase at the time of short circuit, then the transformer primary side voltage in the grid is If the equivalent resistance and leakage inductance are R k With L k , then the expression of the circuit equation under the short-circuit state is:

[0135]

[0136] Where, u1 represents the primary side voltage of the winding; U 1N represents the effective value of the primary voltage; w represents the frequency of the sine wave; t represents the time; α represents the initial phase of the voltage at the time of short circuit; R kRepresents the equivalent resistance of a short-circuit circuit; i k Indicates short-circuit current; L k Represents the equivalent leakage inductance of a short-circuit circuit.

[0137] Step S10333, solving the circuit equation under the short-circuit state to obtain the steady-state component and transient component of the short-circuit current;

[0138] Specifically, in i k | t＝0 = 0, the short-circuit current i k The initial value (time t = 0) is zero, so:

[0139]

[0140] In the formula, is the steady-state component of the short-circuit current, is the transient component of the short-circuit current; I k Indicates the effective value of the transient component of the circuit current.

[0141] Step S10335: Calculate the maximum short-circuit current based on the steady-state component of the short-circuit current.

[0142] Specifically, when α t＝0 =0, the maximum short-circuit current can be obtained. The expression of the maximum short-circuit current is:

[0143]

[0144] In the formula, i k Indicates the maximum short-circuit current; I k Indicates the effective value of short-circuit current.

[0145] In this optional embodiment, based on the magnetic induction intensity and the circuit equation under the short-circuit state, the electromagnetic load is calculated, and by constructing the deformation variable solution equation, the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed, and the maximum axial deformation variable and the maximum amplitude deformation variable of the winding are obtained, including:

[0146] Step S10351, calculating the electromagnetic load acting on the winding in the short-circuit state according to the magnetic induction intensity and the circuit equation in the short-circuit state;

[0147] Specifically, the transformer winding deformation is calculated. Based on the transformer magnetic field calculation and transformer short-circuit current calculation, the electromagnetic load on the winding under short-circuit conditions can be calculated:

[0148]

[0149] Where, f represents the electromagnetic load on the winding under short circuit; J represents the current density; B represents the magnetic leakage component in the winding space; i represents the unit vector in the x-axis direction of the rectangular coordinate system in the vector cross multiplication operation; j represents the unit vector in the y-axis direction of the rectangular coordinate system in the vector cross multiplication operation; k represents the unit vector in the z-axis direction of the rectangular coordinate system in the vector cross multiplication operation; J x represents the x-axis component of the current density J in the rectangular coordinate system; J y represents the component of the current density J on the y-axis in the rectangular coordinate system; J z represents the component of the current density J on the z-axis in the rectangular coordinate system; B x represents the component of the leakage magnetic component B on the x-axis in the rectangular coordinate system; B y It represents the component of leakage magnetic component B on the y-axis in the rectangular coordinate system; B z It represents the component of the leakage magnetic component B on the z-axis in the rectangular coordinate system.

[0150] Step S10353, based on the electromagnetic load acting on the winding in the short-circuit state, establish a definite solution equation for the deformation of the winding;

[0151] Specifically, under the action of electromagnetic load, the deformation equation generated by the winding is:

[0152] [K e ]{S0}={P0}

[0153] In the formula, [K e ] is the stiffness matrix of the winding; {S0} is the structural deformation vector of the winding; {P0} is the electromagnetic load vector of the winding.

[0154] Step S10355: Based on the deformation equation of the winding, a balance equation of the winding under state change is constructed, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed;

[0155] Specifically, when the value of {S0} is small, the equilibrium equation of the winding under the state change (P, S, σ) is:

[0156] [K e +K σ (σ)]{ΔS}={ΔP}

[0157] In the formula, [K σ (σ)] is the initial stress matrix of the winding; σ is the change in stress value caused by the deformation variable {S0}; {ΔS} represents the change in the deformation vector of the winding structure; {ΔP} represents the change in the node vector of the electromagnetic load received by the winding structure.

[0158] Step S10357: Calculate the deformation characteristics of the winding under the electromagnetic load, including the maximum axial deformation and the maximum radial deformation of the winding, based on the linear relationship analysis result between the deformation state of the winding and the electromagnetic load.

[0159] Specifically, in a small deformation range, the deformation state of the winding is linearly related to the electromagnetic load:

[0160]

[0161] Where μ represents the proportionality coefficient and σ0 represents the stress value.

[0162] Specifically, the deformation characteristics of the winding under the electromagnetic load can be calculated, such as Figure 4 As shown in the figure, when calculating the deformation characteristics of transformer windings in multi-physics fields, the electromagnetic-structural direct coupling method is usually used. The electromagnetic load is a non-contact stress that is directly coupled into the structural field as a load on the structural field.

[0163] Figure 5 An embodiment of a turn-to-turn short-circuit fault deduction system of the present invention is shown.

[0164] In this optional embodiment, the turn-to-turn short circuit fault deduction system includes:

[0165] The mechanical property database unit 201 is used to determine the critical strain value of insulation material rupture based on the axial tensile test data and the radial tensile test data of the winding, and to establish a mechanical property database of strain and insulation rupture according to the critical strain value of insulation material rupture;

[0166] The multi-physics field calculation unit 203 is used to establish a Poisson boundary value equation group according to the magnetic vector potential and the current density, and obtain the magnetic induction intensity by solving the Poisson boundary value equation group; establish a circuit equation under a short-circuit state based on the primary side voltage of the winding, and analyze the linear relationship between the deformation state of the winding and the electromagnetic load in combination with the magnetic induction intensity to obtain the maximum axial deformation and the maximum amplitude deformation of the winding;

[0167] The turn-to-turn short circuit diagnosis unit 205 is used to determine the turn-to-turn insulation state and turn-to-turn short circuit situation based on the mechanical property database of strain and insulation rupture, combined with the maximum axial deformation and the maximum amplitude deformation of the winding.

[0168] In this optional embodiment, the multi-physics field calculation unit 203 includes:

[0169] A magnetic field distribution calculation module (not shown in the figure) is used to establish a Poisson boundary value equation group based on the magnetic vector potential and the current density, solve the Poisson boundary value equation group to obtain the magnetic vector value of each node in the winding space, and calculate the magnetic induction intensity based on the magnetic vector value;

[0170] A short-circuit current calculation module (not shown in the figure) is used to establish a circuit equation under a short-circuit state according to the primary side voltage of the winding, and calculate the maximum short-circuit current by solving the circuit equation under the short-circuit state;

[0171] The winding deformation calculation module (not shown in the figure) is used to calculate the electromagnetic load based on the magnetic induction intensity and the circuit equation under the short-circuit state, and to analyze the linear relationship between the deformation state of the winding and the electromagnetic load by constructing a deformation variable equation to obtain the maximum axial deformation and maximum radial deformation of the winding.

[0172] In this optional embodiment, a Poisson boundary value equation group is established based on the magnetic vector potential and the current density, the Poisson boundary value equation group is solved to obtain the magnetic vector value of each node in the winding space, and the magnetic induction intensity is calculated according to the magnetic vector value, including:

[0173] Based on the component relationship between magnetic vector potential and current density, a Poisson boundary value equation system including boundary conditions and magnetic resistivity is established.

[0174] The Poisson boundary value equations are solved by numerical iteration method to obtain the magnetic vector value of each node.

[0175] According to the magnetic vector value of each node, the corresponding magnetic induction intensity is calculated to obtain the leakage magnetic components in all directions in the winding space.

[0176] In this optional embodiment, the expression of the Poisson boundary value equation system including boundary conditions and magnetic resistivity is:

[0177]

[0178] Wherein, Ω represents the solution domain of the winding space; z represents the axial coordinate in the cylindrical coordinate system; r represents the radial coordinate in the cylindrical coordinate system; θ represents the angular coordinate in the cylindrical coordinate system; t represents time; A θ represents magnetic vector potential; A represents magnetic vector potential; σ′ represents electrical conductivity; J sθ represents the current density J s The θ component of the θ component; Γ1 represents the first type of boundary; Γ2 represents the second type of boundary; H t represents its tangent vector; n represents the normal direction; t0 represents the initial time; ν′ represents the magnetic reluctance; f0(z,r) represents the functional relationship between the magnetic vector potential at time t0 and the coordinates z and r;

[0179] The expression of the magnetic vector value of each node is:

[0180] B=▽×A

[0181] In the formula, A represents the magnetic vector value of each node; ▽ represents the vector differential operator; B represents the isotropic leakage magnetic component in the winding space;

[0182] The expression of the isotropic leakage magnetic component in the winding space is:

[0183]

[0184] Wherein, B represents the isotropic leakage magnetic component of the winding space; x represents the abscissa in the rectangular coordinate system; y represents the ordinate in the rectangular coordinate system.

[0185] In this optional embodiment, according to the primary side voltage of the winding, a circuit equation under a short-circuit state is established, and by solving the circuit equation under a short-circuit state, the maximum short-circuit current is calculated, which includes:

[0186] According to the primary voltage of the winding, combined with the equivalent resistance and the equivalent leakage inductance, the circuit equation under the short-circuit state is established;

[0187] Solve the circuit equation under short-circuit state to obtain the steady-state component and transient component of the short-circuit current;

[0188] Based on the steady-state component of the short-circuit current, the maximum short-circuit current is calculated.

[0189] In this optional embodiment, the expression of the circuit equation in the short-circuit state is:

[0190]

[0191] Where, u1 represents the primary side voltage of the winding; U 1N represents the effective value of the primary voltage; w represents the frequency of the sine wave; t represents the time; α represents the initial phase of the voltage at the time of short circuit; R k Represents the equivalent resistance of a short-circuit circuit; i k Indicates short-circuit current; L k Indicates the equivalent leakage inductance of a short circuit;

[0192] The expression for the maximum short-circuit current is:

[0193]

[0194] In the formula, i k Indicates the maximum short-circuit current; I k Indicates the effective value of short-circuit current.

[0195] In this optional embodiment, based on the magnetic induction intensity and the circuit equation under the short-circuit state, the electromagnetic load is calculated, and by constructing the deformation variable solution equation, the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed, and the maximum axial deformation variable and the maximum amplitude deformation variable of the winding are obtained, including:

[0196] According to the magnetic induction intensity and the circuit equation under the short-circuit state, the electromagnetic load acting on the winding under the short-circuit state is calculated;

[0197] Based on the electromagnetic load acting on the winding under the short-circuit state, the deformation equation of the winding is established;

[0198] Based on the deformation equation of the winding, the equilibrium equation of the winding under state change is constructed, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed;

[0199] According to the analysis results of the linear relationship between the deformation state of the winding and the electromagnetic load, the deformation characteristics of the winding under the electromagnetic load are calculated, including the maximum axial deformation and the maximum radial deformation of the winding.

[0200] In one embodiment, a computer device is provided. The computer device may be a server, and its internal structure diagram may be as follows: Figure 6 As shown. The computer device includes a processor, a memory and a network interface connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store static information and dynamic information data. The network interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, the steps in the above method embodiment are implemented.

[0201] Those skilled in the art will understand that Figure 6 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present invention, and does not constitute a limitation on the computer device to which the solution of the present invention is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0202] In addition, the present invention also provides a computer device, including a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps in the above method embodiment when executing the computer program.

[0203] In addition, the present invention further provides a computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, the steps in the above method embodiment are implemented.

[0204] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided by the present invention can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory or optical memory, etc. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0205] The present invention is not limited to the structures which have been described above and shown in the drawings, and various modifications and changes may be made without departing from the scope thereof. The scope of the present invention is limited only by the appended claims.

Claims

1. A turn-to-turn short-circuit fault deduction method, characterized in that: The turn-to-turn short-circuit fault deduction method includes: Based on the axial tensile test data and radial tensile test data of the winding, the critical strain value of the insulation material rupture is determined, and according to the critical strain value of the insulation material rupture, a database of strain and mechanical properties of insulation rupture is established; A group of Poisson boundary value equations is established according to the magnetic vector potential and current density, and the magnetic induction intensity is obtained by solving the group of Poisson boundary value equations; a circuit equation under a short-circuit state is established based on the primary side voltage of the winding, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed in combination with the magnetic induction intensity, and the maximum axial deformation and the maximum amplitude deformation of the winding are obtained; Based on the mechanical property database of strain and insulation rupture, combined with the maximum axial deformation and maximum amplitude deformation of the winding, the inter-turn insulation status and inter-turn short circuit situation are judged.

2. A turn-to-turn short-circuit fault deduction method according to claim 1, characterized in that: The Poisson boundary value equations are established according to the magnetic vector potential and the current density, and the magnetic induction intensity is obtained by solving the Poisson boundary value equations; the circuit equation under the short-circuit state is established based on the primary side voltage of the winding, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed in combination with the magnetic induction intensity, and the maximum axial deformation variable and the maximum amplitude deformation variable of the winding are obtained, including: Based on the magnetic vector potential and current density, a Poisson boundary value equation group is established, the magnetic vector value of each node in the winding space is obtained by solving the Poisson boundary value equation group, and the magnetic induction intensity is calculated based on the magnetic vector value; According to the primary side voltage of the winding, a circuit equation under the short-circuit state is established, and the maximum short-circuit current is calculated by solving the circuit equation under the short-circuit state; The electromagnetic load is calculated based on the magnetic induction intensity and the circuit equation under the short-circuit state. The linear relationship between the deformation state of the winding and the electromagnetic load is analyzed by constructing the deformation variable equation to obtain the maximum axial deformation and the maximum radial deformation of the winding.

3. A turn-to-turn short-circuit fault deduction method according to claim 2, characterized in that: The method of establishing a Poisson boundary value equation group based on the magnetic vector potential and the current density, solving the Poisson boundary value equation group to obtain the magnetic vector value of each node in the winding space, and calculating the magnetic induction intensity according to the magnetic vector value includes: Based on the component relationship between magnetic vector potential and current density, a Poisson boundary value equation system including boundary conditions and magnetic resistivity is established; The Poisson boundary value equations are solved by numerical iteration method to obtain the magnetic vector value of each node. According to the magnetic vector value of each node, the corresponding magnetic induction intensity is calculated to obtain the leakage magnetic components in all directions in the winding space.

4. A turn-to-turn short-circuit fault deduction method according to claim 3, characterized in that: The expression of the Poisson boundary value equations including boundary conditions and magnetic resistivity is: Wherein, Ω represents the solution domain of the winding space; z represents the axial coordinate in the cylindrical coordinate system; r represents the radial coordinate in the cylindrical coordinate system; θ represents the angular coordinate in the cylindrical coordinate system; t represents time; A θ represents magnetic vector potential; A represents magnetic vector potential; σ′ represents electrical conductivity; J sθ represents the current density J s The θ component of the θ component; Γ1 represents the first type of boundary; Γ2 represents the second type of boundary; H t represents its tangent vector; n represents the normal direction; t0 represents the initial time; ν′ represents the magnetic reluctance; f0(z,r) represents the functional relationship between the magnetic vector potential at time t0 and the coordinates z and r; The expression of the magnetic vector value of each node is: Where A represents the magnetic vector value of each node; represents the vector differential operator; B represents the isotropic leakage magnetic component in the winding space; The expression of the isotropic leakage magnetic component of the winding space is: Wherein, B represents the isotropic leakage magnetic component of the winding space; x represents the abscissa in the rectangular coordinate system; y represents the ordinate in the rectangular coordinate system.

5. The turn-to-turn short-circuit fault deduction method according to claim 2, characterized in that: The establishing of a circuit equation under a short-circuit state according to the primary side voltage of the winding, and calculating the maximum short-circuit current by solving the circuit equation under the short-circuit state comprises: According to the primary voltage of the winding, combined with the equivalent resistance and the equivalent leakage inductance, the circuit equation under the short-circuit state is established; Solve the circuit equation under short-circuit state to obtain the steady-state component and transient component of the short-circuit current; Based on the steady-state component of the short-circuit current, the maximum short-circuit current is calculated.

6. A turn-to-turn short-circuit fault deduction method according to claim 5, characterized in that: The expression of the circuit equation under the short-circuit state is: Where, u1 represents the primary side voltage of the winding; U 1N represents the effective value of the primary voltage; w represents the frequency of the sine wave; t represents the time; α represents the initial phase of the voltage at the time of short circuit; R k Represents the equivalent resistance of a short-circuit circuit; i k Indicates short-circuit current; L k Indicates the equivalent leakage inductance of a short circuit; The expression of the maximum short-circuit current is: In the formula, i k Indicates the maximum short-circuit current; I k Indicates the effective value of short-circuit current.

7. A turn-to-turn short-circuit fault deduction method according to claim 2, characterized in that: The electromagnetic load is calculated based on the magnetic induction intensity and the circuit equation under the short-circuit state, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed by constructing the deformation variable equation to obtain the maximum axial deformation variable and the maximum radial deformation variable of the winding. According to the magnetic induction intensity and the circuit equation under the short-circuit state, the electromagnetic load acting on the winding under the short-circuit state is calculated; Based on the electromagnetic load acting on the winding under the short-circuit state, the deformation equation of the winding is established; Based on the deformation equation of the winding, the equilibrium equation of the winding under state change is constructed, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed; According to the analysis results of the linear relationship between the deformation state of the winding and the electromagnetic load, the deformation characteristics of the winding under the electromagnetic load are calculated, including the maximum axial deformation and the maximum radial deformation of the winding.

8. A turn-to-turn short-circuit fault deduction system, characterized in that: The turn-to-turn short-circuit fault deduction system comprises: A mechanical property database unit is used to determine the critical strain value of insulation material rupture based on the axial tensile test data and radial tensile test data of the winding, and to establish a mechanical property database of strain amount and insulation rupture according to the critical strain value of insulation material rupture; The multi-physics field calculation unit is used to establish a Poisson boundary value equation group according to the magnetic vector potential and the current density, and obtain the magnetic induction intensity by solving the Poisson boundary value equation group; establish the circuit equation under the short-circuit state based on the primary side voltage of the winding, and analyze the linear relationship between the deformation state of the winding and the electromagnetic load in combination with the magnetic induction intensity, and obtain the maximum axial deformation and the maximum amplitude deformation of the winding; The turn-to-turn short-circuit diagnostic unit is used to judge the turn-to-turn insulation status and turn-to-turn short-circuit situation based on the mechanical property database of strain and insulation rupture, combined with the maximum axial deformation and maximum amplitude deformation of the winding.

9. The turn-to-turn short-circuit fault deduction system according to claim 8, characterized in that: The multi-physics field calculation unit comprises: The magnetic field distribution calculation module is used to establish a Poisson boundary value equation group based on the magnetic vector potential and current density, solve the Poisson boundary value equation group to obtain the magnetic vector value of each node in the winding space, and calculate the magnetic induction intensity based on the magnetic vector value; A short-circuit current calculation module is used to establish a circuit equation under a short-circuit state according to the primary side voltage of the winding, and calculate the maximum short-circuit current by solving the circuit equation under the short-circuit state; The winding deformation calculation module is used to calculate the electromagnetic load based on the magnetic induction intensity and the circuit equation under the short-circuit state, and to analyze the linear relationship between the deformation state of the winding and the electromagnetic load by constructing the deformation variable equation to obtain the maximum axial deformation and the maximum amplitude deformation of the winding.

10. The turn-to-turn short-circuit fault deduction system according to claim 9, characterized in that: The method of establishing a Poisson boundary value equation group based on the magnetic vector potential and the current density, solving the Poisson boundary value equation group to obtain the magnetic vector value of each node in the winding space, and calculating the magnetic induction intensity according to the magnetic vector value includes: Based on the component relationship between magnetic vector potential and current density, a Poisson boundary value equation system including boundary conditions and magnetic resistivity is established; The Poisson boundary value equations are solved by numerical iteration method to obtain the magnetic vector value of each node. According to the magnetic vector value of each node, the corresponding magnetic induction intensity is calculated to obtain the leakage magnetic components in all directions in the winding space.

11. The turn-to-turn short-circuit fault deduction system according to claim 10, characterized in that: The expression of the Poisson boundary value equations including boundary conditions and magnetic resistivity is: Wherein, Ω represents the solution domain of the winding space; z represents the axial coordinate in the cylindrical coordinate system; r represents the radial coordinate in the cylindrical coordinate system; θ represents the angular coordinate in the cylindrical coordinate system; t represents time; A θ represents magnetic vector potential; A represents magnetic vector potential; σ′ represents electrical conductivity; J sθ represents the current density J s The θ component of the θ component; Γ1 represents the first type of boundary; Γ2 represents the second type of boundary; H t represents its tangent vector; n represents the normal direction; t0 represents the initial time; ν′ represents the magnetic reluctance; f0(z,r) represents the functional relationship between the magnetic vector potential at time t0 and the coordinates z and r; The expression of the magnetic vector value of each node is: Where A represents the magnetic vector value of each node; represents the vector differential operator; B represents the isotropic leakage magnetic component in the winding space; The expression of the isotropic leakage magnetic component of the winding space is: Wherein, B represents the isotropic leakage magnetic component of the winding space; x represents the abscissa in the rectangular coordinate system; y represents the ordinate in the rectangular coordinate system.

12. The turn-to-turn short-circuit fault deduction system according to claim 9, characterized in that: The establishing of a circuit equation under a short-circuit state according to the primary side voltage of the winding, and calculating the maximum short-circuit current by solving the circuit equation under the short-circuit state comprises: According to the primary voltage of the winding, combined with the equivalent resistance and the equivalent leakage inductance, the circuit equation under the short-circuit state is established; Solve the circuit equation under short-circuit state to obtain the steady-state component and transient component of the short-circuit current; Based on the steady-state component of the short-circuit current, the maximum short-circuit current is calculated.

13. The turn-to-turn short-circuit fault deduction system according to claim 12, characterized in that: The expression of the circuit equation under the short-circuit state is: Where, u1 represents the primary side voltage of the winding; U 1N represents the effective value of the primary voltage; w represents the frequency of the sine wave; t represents the time; α represents the initial phase of the voltage at the time of short circuit; R k Represents the equivalent resistance of a short-circuit circuit; i k Indicates short-circuit current; L k Indicates the equivalent leakage inductance of a short circuit; The expression of the maximum short-circuit current is: In the formula, i k Indicates the maximum short-circuit current; I k Indicates the effective value of short-circuit current.

14. The turn-to-turn short-circuit fault deduction system according to claim 9, characterized in that: The electromagnetic load is calculated based on the magnetic induction intensity and the circuit equation under the short-circuit state, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed by constructing the deformation variable equation to obtain the maximum axial deformation variable and the maximum radial deformation variable of the winding. According to the magnetic induction intensity and the circuit equation under the short-circuit state, the electromagnetic load acting on the winding under the short-circuit state is calculated; Based on the electromagnetic load acting on the winding under the short-circuit state, the deformation equation of the winding is established; Based on the deformation equation of the winding, the equilibrium equation of the winding under state change is constructed, and the linear relationship between the deformation state of the winding and the electromagnetic load is analyzed; According to the analysis results of the linear relationship between the deformation state of the winding and the electromagnetic load, the deformation characteristics of the winding under the electromagnetic load are calculated, including the maximum axial deformation and the maximum radial deformation of the winding.

15. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 7 are implemented.

16. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images