CN115659765A - Cable joint temperature field calculation method and device, electronic equipment and storage medium - Google Patents
Cable joint temperature field calculation method and device, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN115659765A CN115659765A CN202211587703.XA CN202211587703A CN115659765A CN 115659765 A CN115659765 A CN 115659765A CN 202211587703 A CN202211587703 A CN 202211587703A CN 115659765 A CN115659765 A CN 115659765A
- Authority
- CN
- China
- Prior art keywords
- flow rate
- temperature
- air flow
- field
- cable
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
The invention discloses a method and a device for calculating a temperature field of a cable joint, electronic equipment and a storage medium, wherein the method comprises the following steps: acquiring initial physical parameters of a cable to be measured, calculating the air flow rate of the next period according to the initial physical parameters through a POD (POD) order reduction technology, judging whether the air flow rate is within an error range or not based on a flow rate threshold value and a flow rate error, if not, setting the air flow rate of the next period as an initial value, returning to the step of calculating the air flow rate, if so, calculating the field temperature of the joint through a windward finite element method, judging whether the field temperature is within the error range or not based on a temperature threshold value and a temperature error, if so, generating a temperature field calculation result, if not, setting the air flow rate of the next period as the initial value, and returning to the step of calculating the air flow rate. The method is beneficial to solving the technical problem that numerical value oscillation occurs when the conventional temperature field calculation method is used for calculating the heat convection problem, and the calculation efficiency of the cable joint temperature field is improved.
Description
Technical Field
The invention relates to the technical field of cable joint temperature fields, in particular to a method and a device for calculating a cable joint temperature field, electronic equipment and a storage medium.
Background
Along with the continuous promotion of high tension transmission line cabling rate, the use quantity of power cable intermediate head also increases by a wide margin. Compared with a cable body, the cable joint is the weakest link in a circuit due to more complex structure. According to the Chinese yearbook statistics data, if the cable fault caused by external force damage is not considered, about 64 percent of the operating cables are faults caused by the intermediate joints, and the type of the faults is basically insulation faults. In an actual operation environment of the cable, the air flow rate plays an important role in heat dissipation of a cable joint, so that it is necessary to consider the air flow rate in real time and calculate a temperature field.
The existing numerical calculation method for the flow field and the temperature field of the 110kV cable and the cable joint mainly comprises a finite volume method, a least square finite element method, a finite element method and the like. When the finite volume method is used for calculating a temperature field, the requirement on a grid is low, the calculation speed is high, but the continuity of the heat flux density of a fluid-solid coupling interface needs to be considered, and the solid-liquid boundary and the fuel tank wall need to be processed separately, which increases the complexity of calculation. The rigidity matrix formed in the least square finite element method calculation is a symmetrical and positive fixed sparse matrix, the numerical stability is high, and when the heat transfer problem is calculated, the flow-solid interface does not need to be processed independently, and the energy conservation equation is automatically satisfied. The formats of the control equations discrete by the least square finite element method are uniform, and the universality of a computer program is strong in the face of different problems. However, the condition number of the rigidity matrix formed in the calculation is large, and improper treatment can greatly increase the iteration number and even cause large numerical deviation. The rigidity matrix formed in the finite element method in the calculation is a sparse matrix, the boundary condition processing is convenient, the adaptability to irregular areas is good, and the method is suitable for the calculation of the temperature field of the multi-layer medium with boundaries and has strong universality. However, when solving the heat transfer problem, the numerical solution may have non-physical oscillation phenomena.
Therefore, in order to improve the calculation efficiency of the cable joint temperature field and solve the technical problem that numerical value oscillation occurs when the existing temperature field calculation method is used for calculating the heat convection problem, it is necessary to construct a cable joint temperature field calculation method.
Disclosure of Invention
The invention provides a cable joint temperature field calculation method, a cable joint temperature field calculation device, electronic equipment and a storage medium, and solves the technical problem that numerical value oscillation occurs when the conventional temperature field calculation method is used for calculating the heat convection problem.
In a first aspect, the present invention provides a method for calculating a temperature field of a cable joint, including:
s1, acquiring initial physical property parameters of a cable to be detected; the initial physical parameters comprise joint temperature, ambient temperature and air flow rate at the initial moment;
s2, calculating the air flow rate of the cable to be tested in the next period according to the joint temperature, the ambient temperature and the air flow rate at the initial moment by using a POD (power POD) order reduction technology;
s3, judging whether the air flow rate of the next period is within an error range or not based on a preset flow rate threshold and a preset flow rate error; if yes, executing step S4; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the step S2;
s4, calculating the field temperature of the integral joint in the cable to be detected by a windward finite element method;
s5, judging whether the field temperature is within an error range or not based on a preset temperature threshold and a preset temperature error; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the step S2; and if so, generating calculation result data of the joint temperature field of the cable to be tested based on the field temperature and the corresponding air flow rate.
Optionally, the step S3 includes:
step S31, calculating to obtain a first difference value between the air flow rate of the next period and the preset flow rate threshold value;
step S32, judging whether the first difference value is smaller than the preset flow speed error; if yes, determining that the air flow rate of the next period is within an error range, and executing a step S4; if not, determining that the air flow rate of the next period is not in the error range, and executing step S33;
and step S33, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to the step S2.
Optionally, the step S5 includes:
step S51, calculating to obtain a second difference value between the field temperature and the preset temperature threshold value;
step S52, judging whether the second difference value is smaller than the preset temperature error; if so, determining that the field temperature is within an error range, and generating calculation result data of the joint temperature field of the cable to be detected; if not, determining that the field temperature is not within the error range, and executing a step S53;
and step S53, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to step S2.
Optionally, the step S52 includes:
step S521, determining whether the second difference is smaller than the preset temperature error; if yes, determining that the field temperature is within an error range, and executing step S522; if not, determining that the field temperature is not in the error range, and executing a step S53;
and S522, summarizing the field temperature and the corresponding air flow rate to generate calculation result data of the joint temperature field of the cable to be tested.
In a second aspect, the present invention provides a cable joint temperature field calculation apparatus, including:
the acquisition module is used for acquiring initial physical property parameters of the cable to be detected; the initial physical parameters comprise joint temperature, environment temperature and air flow rate at the initial moment;
the first calculation module is used for calculating the air flow rate of the cable to be tested in the next period according to the joint temperature, the ambient temperature and the air flow rate at the initial moment by using a POD (power POD) order reduction technology;
the first judgment module is used for judging whether the air flow rate of the next period is within an error range or not based on a preset flow rate threshold value and a preset flow rate error; if yes, executing a second calculation module; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the first calculation module;
the second calculation module is used for calculating the field temperature of the integral joint in the cable to be measured by a windward finite element method;
the second judgment module is used for judging whether the field temperature is within an error range or not based on a preset temperature threshold and a preset temperature error; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the first calculation module; and if so, generating calculation result data of the joint temperature field of the cable to be detected based on the field temperature and the corresponding air flow rate.
Optionally, the first determining module includes:
the first difference submodule is used for calculating to obtain a first difference between the air flow rate of the next period and the preset flow rate threshold;
the first judgment sub-module is used for judging whether the first difference value is smaller than the preset flow rate error or not; if yes, determining that the air flow rate of the next period is within an error range, and executing a second calculation module; if not, determining that the air flow rate of the next period is not in the error range, and executing a first return submodule;
and the first return submodule is used for setting the air flow rate of the next period as the initial value of the air flow rate and returning to execute the first calculation module.
Optionally, the second determining module includes:
the second difference submodule is used for calculating to obtain a second difference between the field temperature and the preset temperature threshold;
the second judgment sub-module is used for judging whether the second difference value is smaller than the preset temperature error or not; if so, determining that the field temperature is within an error range, and generating calculation result data of the joint temperature field of the cable to be detected; if not, determining that the field temperature is not within the error range, and executing a second return submodule;
and the second return submodule is used for setting the air flow rate of the next period as the initial value of the air flow rate and returning to execute the first calculation module.
Optionally, the second determination submodule includes:
the judging unit is used for judging whether the second difference value is smaller than the preset temperature error or not; if yes, determining that the field temperature is within an error range, and executing a generating unit; if not, determining that the field temperature is not in the error range, and executing a second return submodule;
and the generating unit is used for summarizing the field temperature and the corresponding air flow rate thereof and generating calculation result data of the joint temperature field of the cable to be tested.
In a third aspect, the present application provides an electronic device comprising a processor and a memory, wherein the memory stores computer readable instructions, and when the computer readable instructions are executed by the processor, the steps of the method as provided in the first aspect are executed.
In a fourth aspect, the present application provides a storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the method as provided in the first aspect above.
According to the technical scheme, the invention has the following advantages: the invention provides a cable joint temperature field calculation method, which includes the steps of S1, obtaining initial physical parameters of a cable to be measured, wherein the initial physical parameters include joint temperature, environment temperature and air flow rate at an initial moment, S2, calculating air flow rate of the cable to be measured in a next period according to the joint temperature, the environment temperature and the air flow rate at the initial moment through a POD (POD-level-reduction) technology, S3, judging whether the air flow rate of the next period is in an error range or not based on a preset flow rate threshold value and a preset flow rate error, if so, executing S4, if not, setting the air flow rate of the next period as an initial value of the air flow rate, returning to the step S2, step S4, calculating field temperature of an integral joint in the cable to be measured through a windward finite element method, if not, judging whether the field temperature is in the error range or not based on the preset temperature threshold value and the preset temperature error, setting the initial value of the field temperature of the cable to be measured as the initial value, returning to the step S2, if so, generating a result that the air flow rate of the cable joint temperature field temperature is in the cable to be measured, and calculating the problem of the existing air flow rate calculation method is solved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the description of the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without inventive labor.
FIG. 1 is a flowchart illustrating a first embodiment of a method for calculating a temperature field of a cable joint according to the present invention;
FIG. 2 is a flowchart illustrating a second embodiment of a method for calculating a temperature field of a cable connector according to the present invention;
fig. 3 is a schematic diagram illustrating a calculation result of a temperature field of a cable joint structure according to the present invention;
fig. 4 is a block diagram of a cable connector temperature field calculation apparatus according to an embodiment of the present invention.
Detailed Description
The embodiment of the invention provides a cable joint temperature field calculation method, a cable joint temperature field calculation device, electronic equipment and a storage medium, which are used for solving the technical problem that numerical value oscillation occurs when the conventional temperature field calculation method is used for calculating the heat convection problem.
In order to make the objects, features and advantages of the present invention more obvious and understandable, 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, and it is obvious that the embodiments described below are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, fig. 1 is a flowchart of a first embodiment of a method for calculating a temperature field of a cable connector according to the present invention, including:
step S101, acquiring initial physical property parameters of a cable to be detected; the initial physical parameters comprise joint temperature, environment temperature and air flow rate at the initial moment;
step S102, calculating the air flow rate of the cable to be tested in the next period according to the joint temperature, the ambient temperature and the air flow rate at the initial moment by using a POD (power POD) order reduction technology;
step S103, judging whether the air flow rate of the next period is within an error range or not based on a preset flow rate threshold value and a preset flow rate error; if yes, go to step S104; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the step S102;
in this embodiment of the present invention, in step S1031, a first difference between the air flow rate of the next cycle and the preset flow rate threshold is calculated, in step S1032, it is determined whether the first difference is smaller than the preset flow rate error, if so, the air flow rate of the next cycle is determined to be within the error range, step S104 is executed, if not, the air flow rate of the next cycle is determined not to be within the error range, step S1033 is executed, the air flow rate of the next cycle is set as an initial value of the air flow rate, and step S2 is executed again.
Step S104, calculating the field temperature of the integral joint in the cable to be tested by using a windward finite element method;
step S105, judging whether the field temperature is within an error range or not based on a preset temperature threshold and a preset temperature error; if so, generating calculation result data of the joint temperature field of the cable to be detected based on the field temperature and the corresponding air flow rate; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the step S102;
in this embodiment of the present invention, in step S1051, a second difference between the field temperature and the preset temperature threshold is calculated, in step S1052, it is determined whether the second difference is smaller than the preset temperature error, if yes, the field temperature is determined to be within an error range, and calculation result data of the joint temperature field of the cable to be tested is generated, and if no, the field temperature is determined not to be within the error range, step S1053 is executed, in step S1053, the air flow rate of the next period is set as an initial value of the air flow rate, and the step S102 is executed again.
In the method for calculating the temperature field of the cable joint provided by the embodiment of the invention, through the step S101, the initial physical parameters of the cable to be measured are obtained, the initial physical parameters include the joint temperature, the ambient temperature and the air flow rate at the initial moment, the step S102 is to calculate the air flow rate of the cable to be measured in the next period according to the joint temperature, the ambient temperature and the air flow rate at the initial moment through the POD order reduction technology, the step S103 is to judge whether the air flow rate of the next period is within the error range based on the preset flow rate threshold and the preset flow rate error, if so, the step S104 is to be executed, if not, the air flow rate of the next period is set as the initial value of the air flow rate, the step S102 is to be executed again, the step S104 is to calculate the field temperature of the integral joint in the cable to be measured through the windward finite element method, the step S105 is to judge whether the field temperature is within the error range based on the preset temperature threshold and the preset temperature error, if not, the field temperature and the air flow rate of the cable to be calculated based on the current air flow rate calculation method for solving the problem that the calculation of the cable to be measured has occurred in the heat transfer field.
In a second embodiment, referring to fig. 2, fig. 2 is a flowchart illustrating a method for calculating a temperature field of a cable connector according to the present invention, including:
step S201, acquiring initial physical property parameters of a cable to be detected; the initial physical parameters comprise joint temperature, ambient temperature and air flow rate at the initial moment;
step S202, calculating the air flow rate of the cable to be tested in the next period according to the joint temperature, the environment temperature and the air flow rate at the initial moment by using a POD (power POD) order reduction technology;
in the embodiment of the invention, according to the joint temperature, the ambient temperature and the air flow rate at the initial time of the initial physical parameters, the air flow rate of the cable to be tested in the next period is calculated by the POD reduction technology.
In a specific implementation, the principle of the POD order reduction method based on singular value decomposition is as follows:
a snapshot matrix of
Then there is a diagonal matrixDOrthogonal matrix
So that
. Then it can be deduced
The results are as follows:
so vector
Has a length of
Wherein
Is composed of
Is determined by the characteristic value of (a),
is composed of
The feature vector of (2). Note the book
Is composed of
Unit direction vector of (2) and records
It can be derived that:
then the
Called a matrixBThe singular value of (a). Will be provided with
The eigenvalues of (2) are arranged from large to small, and the singular values are
Written in a matrix form, the calculation formula is specifically as follows:
and matrixBVAlso is provided with
A non-zero vector, memory
，
According to the above calculation
Is calculated by
Give a result of
Thus, a vector is derived
、
Also orthogonal, the calculation formula is specifically:
the above formula is transformed into
This will be a non-square matrixBSingular value decomposition is performed and the matrix
The quadratic root of the eigenvalue is the matrixBSingular values of, matrix
Is the matrixBThe singular vectors of (a). The matrix of singular vectors at this time is the reduced subspace composed of POD orthogonal bases required for the reduction.
Defining the solution of projection of sample data to the reduced subspace as:
wherein the content of the first and second substances,
is a collection of sample data;
is the orthonormal basis for the reduced order subspace.
D eigenvectors are selected as orthogonal bases of the reduced order subspace, and the matrix formed by the d eigenvectors is
，
Then the solution vector corresponding to the discrete nodes of the field domain at a certain time
The construction using the selected orthogonal basis weights is specifically:
wherein, the first and the second end of the pipe are connected with each other,
corresponding to the coordinates of the respective POD orthogonal bases.
Step S203, calculating to obtain a first difference value between the air flow rate of the next period and the preset air flow rate threshold;
step S204, judging whether the first difference value is smaller than the preset flow speed error; if yes, determining that the air flow rate of the next period is within the error range, and executing step S206; if not, determining that the air flow rate of the next period is not in the error range, and executing step S205;
in this embodiment of the present invention, when the first difference is smaller than the preset flow rate error, it is determined that the air flow rate of the next cycle is within the error range, step S206 is performed, and when the first difference is not smaller than the preset flow rate error, it is determined that the air flow rate of the next cycle is not within the error range, step S205 is performed.
Step S205, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to step S202;
step S206, calculating the field temperature of the integral joint in the cable to be tested by a windward finite element method;
in the embodiment of the invention, the field temperature of the integral joint in the cable to be tested is calculated by adopting an upwind finite element method.
In a specific implementation, the windward finite element method principle is as follows:
the control equation of the transient temperature field is specifically as follows:
wherein the content of the first and second substances,Tis the temperature of point (x, y);
is a measure of the thermal conductivity of the material,
heat generation rate per unit volume; in the case of a non-heat-source region,
then 0 is taken. The boundary conditions of the whole field and the coupling conditions between the fluid and the solid are specifically as follows:
wherein, the first and the second end of the pipe are connected with each other,
is the temperature of the zone boundary;
is the air medium temperature;
is the convective heat transfer coefficient;
is the outer boundary of the field;
、
respectively representing the fluid side and solid side boundaries on the fluid-solid interface;
、
representing the thermal conductivity of the solid material and the fluid material, respectively. The 1 st equation is a class 1 boundary condition, the 2 nd equation is a class 2 and 3 boundary condition, and the 3 rd and 4 th equations are a fluid-solid coupling interface condition.
The invention adopts the shape function of 9 node quadrilateral units to calculate the temperature field, and the calculation formula of the temperature field in 1 unit is specifically as follows:
wherein the content of the first and second substances,
is the temperature value corresponding to the ith node in the unit e.
The finite element form of the energy conservation equation is specifically:
wherein the content of the first and second substances,Wrepresenting a weight function;
indicating the heat flow density. In calculating the convection diffusion equation, theThe flow terms can cause the solution process to be unstable. After the flow terms are subjected to numerical discretization, the coupling between adjacent nodes is eliminated, so that poor numerical stability and non-physical oscillation can be caused in the solving process.
Therefore, the invention uses the finite element method of the windward format to add the convection term as the correction quantity into the weight function, and the format has good numerical stability, can effectively eliminate the non-physical numerical oscillation, and can not increase the calculation quantity. The weight function form of the windward finite element is specifically as follows:
wherein the content of the first and second substances,
is a windward factor;
is the characteristic length of the cell. The asymmetry of the weight function increases the contribution to the flow term without contributing to the time and diffusion terms, thereby reducing non-physical oscillations due to greater diffusion. The equation for solving the transient temperature field is specifically as follows:
wherein the content of the first and second substances,
characterizing the boundary satisfying the 2 nd and 3 rd boundary conditions;
is ambient temperature. Next, derivative with time
The treatment is carried out, wherein,Tis a matrix of node temperatures at different times. Dispersing the time items and introducing parameters
Then, the temperature matrix of the k +1 th time step is specifically:
the final formula is as follows:
wherein, the first and the second end of the pipe are connected with each other,
by the above formula, the temperature field distribution at each time can be obtained in combination with the boundary conditions.
Step S207, calculating to obtain a second difference value between the field temperature and the preset temperature threshold value;
step S208, judging whether the second difference value is smaller than the preset temperature error; if so, determining that the field temperature is within an error range, and generating calculation result data of the joint temperature field of the cable to be detected; if not, determining that the field temperature is not within the error range, and executing step S209;
in an alternative embodiment, the step S208 includes:
step S2081, judging whether the second difference value is smaller than the preset temperature error; if yes, determining that the field temperature is within the error range, and executing step S2082; if not, determining that the field temperature is not within the error range, and executing step S209;
and S2082, summarizing the field temperature and the corresponding air flow rate, and generating calculation result data of the joint temperature field of the cable to be tested.
In the embodiment of the invention, when the second difference is smaller than the preset temperature error, the field temperature is determined to be within the error range, the field temperature and the corresponding air flow rate are summarized, and calculation result data of the joint temperature field of the cable to be measured are generated; when the second difference is not less than the preset temperature error, it is determined that the field temperature is not within the error range, and step S209 is performed.
Step S209, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to step S202;
in this embodiment of the present invention, when the second difference is not less than the preset temperature error, it is determined that the field temperature is not within the error range, the air flow rate of the next period is set as the initial value of the air flow rate, and the process returns to step S202.
In specific implementation, please refer to fig. 3, fig. 3 is a schematic diagram of a calculation result of a temperature field of a cable joint structure according to the present invention, in which 301 is a cable core, 302 is an XLPE cable, 303 is a water blocking tape, 304 is an aluminum sheath, 305 is an outer sheath, a ordinate 306 is temperature/deg.c, a abscissa 307 is time/s, a cable model of a reference calculation model is YJLW 03-Z64/110kV 630mm2, an initial air flow rate is set to 2m/s, and parameters of a cable joint are specifically shown in the following table:
according to the parameters, the calculation is performed by combining the cable joint temperature field calculation method provided by the invention, and the obtained temperature field calculation result is specifically shown in fig. 3.
In the method for calculating the temperature field of the cable joint provided by the embodiment of the invention, through the step S1, the initial physical parameters of the cable to be measured are obtained, the initial physical parameters comprise the joint temperature, the environment temperature and the air flow rate at the initial moment, the step S2, the air flow rate of the cable to be measured in the next period is calculated according to the joint temperature, the environment temperature and the air flow rate at the initial moment through the POD order reduction technology, the step S3, whether the air flow rate of the next period is in an error range is judged based on a preset flow rate threshold and a preset flow rate error, if yes, the step S4 is executed, if no, the air flow rate of the next period is set as the initial value of the air flow rate, the step S2 and the step S4 are executed again, the field temperature of the integral joint in the cable to be measured is calculated through the windward finite element method, the step S5, whether the field temperature is in the error range is judged based on the preset temperature threshold and the preset temperature error, if no, the field temperature of the field temperature and the air flow rate of the cable to be measured are set as the initial value, the calculation method for calculating the air flow rate of the cable joint temperature, and the air flow rate calculation method for calculating the current air flow rate of the cable to be measured by the existing technology, the calculation method for calculating the problem of the cable to be measured can be solved.
Referring to fig. 4, fig. 4 is a block diagram of a cable connector temperature field calculating device according to an embodiment of the present invention, including:
an obtaining module 401, configured to obtain an initial physical property parameter of a cable to be tested; the initial physical parameters comprise joint temperature, ambient temperature and air flow rate at the initial moment;
a first calculating module 402, configured to calculate, by using a POD order reduction technology, an air flow rate of the cable to be tested in a next period according to the joint temperature, the ambient temperature, and the air flow rate at the initial time;
a first judging module 403, configured to judge whether the air flow rate of the next cycle is within an error range based on a preset flow rate threshold and a preset flow rate error; if yes, execute the second calculation module 404; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the first calculation module 402;
a second calculating module 404, configured to calculate a field temperature of an integral joint in the cable to be measured by using a windward finite element method;
a second determining module 405, configured to determine whether the field temperature is within an error range based on a preset temperature threshold and a preset temperature error; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the first calculation module; if yes, generating calculation result data 402 of the joint temperature field of the cable to be tested based on the field temperature and the corresponding air flow rate.
In an optional embodiment, the first determining module 403 includes:
the first difference submodule is used for calculating to obtain a first difference between the air flow rate of the next period and the preset flow rate threshold;
the first judgment submodule is used for judging whether the first difference value is smaller than the preset flow rate error or not; if yes, determining that the air flow rate of the next period is within an error range, and executing a second calculation module; if not, determining that the air flow rate of the next period is not in the error range, and executing a first return submodule;
and the first return submodule is used for setting the air flow rate of the next period as the initial value of the air flow rate and returning to execute the first calculation module.
In an optional embodiment, the second determining module 405 includes:
the second difference submodule is used for calculating to obtain a second difference between the field temperature and the preset temperature threshold;
the second judgment submodule is used for judging whether the second difference value is smaller than the preset temperature error or not; if so, determining that the field temperature is within an error range, and generating calculation result data of the joint temperature field of the cable to be detected; if not, determining that the field temperature is not within the error range, and executing a second return submodule;
and the second return submodule is used for setting the air flow rate of the next period as the initial value of the air flow rate and returning to execute the first calculation module.
In an optional embodiment, the second determination sub-module includes:
the judging unit is used for judging whether the second difference value is smaller than the preset temperature error or not; if yes, determining that the field temperature is within an error range, and executing a generating unit; if not, determining that the field temperature is not within the error range, and executing a second return submodule;
and the generation unit is used for summarizing the field temperature and the corresponding air flow rate thereof and generating calculation result data of the joint temperature field of the cable to be detected.
An embodiment of the present invention further provides an electronic device, which includes a memory and a processor, where the memory stores a computer program, and when the computer program is executed by the processor, the processor is enabled to execute the steps of the cable connector temperature field calculation method according to any of the above embodiments.
The embodiment of the present invention further provides a computer storage medium, on which a computer program is stored, where the computer program is executed by the processor to implement the steps of the method for calculating a temperature field of a cable joint according to any one of the above embodiments.
It can be clearly understood by those skilled in the art that, for convenience and simplicity of description, the specific working processes of the above-described systems, apparatuses and units may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the several embodiments provided in the present application, it should be understood that the method, apparatus, electronic device and storage medium disclosed in the present invention may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the units is only one type of logical functional division, and other divisions may be realized in practice, for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of a software functional unit.
The integrated unit, if implemented as a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a readable storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned readable storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-only memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
The above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; although the present invention has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A cable joint temperature field calculation method is characterized by comprising the following steps:
s1, acquiring initial physical property parameters of a cable to be detected; the initial physical parameters comprise joint temperature, ambient temperature and air flow rate at the initial moment;
s2, calculating the air flow rate of the cable to be tested in the next period according to the joint temperature, the ambient temperature and the air flow rate at the initial moment by using a POD (power POD) order reduction technology;
s3, judging whether the air flow rate of the next period is within an error range or not based on a preset flow rate threshold and a preset flow rate error; if yes, executing step S4; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the step S2;
s4, calculating the field temperature of the integral joint in the cable to be detected by a windward finite element method;
s5, judging whether the field temperature is within an error range or not based on a preset temperature threshold and a preset temperature error; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the step S2; and if so, generating calculation result data of the joint temperature field of the cable to be tested based on the field temperature and the corresponding air flow rate.
2. The cable joint temperature field calculation method according to claim 1, wherein the step S3 includes:
step S31, calculating to obtain a first difference value between the air flow rate of the next period and the preset flow rate threshold value;
step S32, judging whether the first difference value is smaller than the preset flow speed error; if yes, determining that the air flow rate of the next period is within an error range, and executing a step S4; if not, determining that the air flow rate of the next period is not in the error range, and executing a step S33;
and step S33, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to the step S2.
3. The cable joint temperature field calculation method according to claim 1, wherein the step S5 includes:
step S51, calculating to obtain a second difference value between the field temperature and the preset temperature threshold value;
step S52, judging whether the second difference value is smaller than the preset temperature error; if so, determining that the field temperature is within an error range, and generating calculation result data of the joint temperature field of the cable to be detected; if not, determining that the field temperature is not within the error range, and executing a step S53;
and S53, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to the step S2.
4. The cable joint temperature field calculation method according to claim 3, wherein the step S52 includes:
step S521, judging whether the second difference value is smaller than the preset temperature error; if yes, determining that the field temperature is within the error range, and executing step S522; if not, determining that the field temperature is not in the error range, and executing a step S53;
and S522, summarizing the field temperature and the corresponding air flow rate to generate calculation result data of the joint temperature field of the cable to be tested.
5. A cable joint temperature field calculation device, comprising:
the acquisition module is used for acquiring initial physical property parameters of the cable to be detected; the initial physical parameters comprise joint temperature, environment temperature and air flow rate at the initial moment;
the first calculation module is used for calculating the air flow rate of the cable to be tested in the next period according to the joint temperature, the ambient temperature and the air flow rate at the initial moment by using a POD (power POD) order reduction technology;
the first judgment module is used for judging whether the air flow rate of the next period is within an error range or not based on a preset flow rate threshold value and a preset flow rate error; if yes, executing a second calculation module; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the first calculation module;
the second calculation module is used for calculating the field temperature of the integral joint in the cable to be measured by a windward finite element method;
the second judgment module is used for judging whether the field temperature is within an error range or not based on a preset temperature threshold and a preset temperature error; if not, setting the air flow rate of the next period as the initial value of the air flow rate, and returning to execute the first calculation module; and if so, generating calculation result data of the joint temperature field of the cable to be tested based on the field temperature and the corresponding air flow rate.
6. The cable joint temperature field calculation device according to claim 5, wherein the first determination module includes:
the first difference submodule is used for calculating to obtain a first difference between the air flow rate of the next period and the preset flow rate threshold;
the first judgment sub-module is used for judging whether the first difference value is smaller than the preset flow rate error or not; if yes, determining that the air flow rate of the next period is within an error range, and executing a second calculation module; if not, determining that the air flow rate of the next period is not in the error range, and executing a first return submodule;
and the first return submodule is used for setting the air flow rate of the next period as the initial value of the air flow rate and returning to execute the first calculation module.
7. The cable joint temperature field calculation device of claim 5, wherein the second determination module comprises:
the second difference submodule is used for calculating to obtain a second difference between the field temperature and the preset temperature threshold;
the second judgment submodule is used for judging whether the second difference value is smaller than the preset temperature error or not; if so, determining that the field temperature is within an error range, and generating calculation result data of the joint temperature field of the cable to be detected; if not, determining that the field temperature is not within the error range, and executing a second return submodule;
and the second return submodule is used for setting the air flow rate of the next period as the initial value of the air flow rate and returning to execute the first calculation module.
8. The cable joint temperature field calculation device of claim 7, wherein the second determination submodule includes:
the judging unit is used for judging whether the second difference value is smaller than the preset temperature error or not; if yes, determining that the field temperature is within an error range, and executing a generating unit; if not, determining that the field temperature is not in the error range, and executing a second return submodule;
and the generation unit is used for summarizing the field temperature and the corresponding air flow rate thereof and generating calculation result data of the joint temperature field of the cable to be detected.
9. An electronic device comprising a processor and a memory storing computer readable instructions that, when executed by the processor, perform the method of any one of claims 1-4.
10. A storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, performs the method according to any of claims 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211587703.XA CN115659765B (en) | 2022-12-12 | 2022-12-12 | Cable joint temperature field calculation method and device, electronic equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211587703.XA CN115659765B (en) | 2022-12-12 | 2022-12-12 | Cable joint temperature field calculation method and device, electronic equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115659765A true CN115659765A (en) | 2023-01-31 |
| CN115659765B CN115659765B (en) | 2023-04-14 |
Family
ID=85017487
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211587703.XA Active CN115659765B (en) | 2022-12-12 | 2022-12-12 | Cable joint temperature field calculation method and device, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115659765B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001351919A (en) * | 2000-06-05 | 2001-12-21 | Nec Corp | Wiring fault analysis method |
| JP2006250689A (en) * | 2005-03-10 | 2006-09-21 | Tokyo Electric Power Co Inc:The | Cable conductor temperature estimation method considering movement of air in cable tunnel, cable conductor temperature estimation system, and cable conductor temperature estimation program |
| CN106066212A (en) * | 2016-05-27 | 2016-11-02 | 三峡大学 | A kind of cable conductor temperature indirect measurement method |
| CN109255192A (en) * | 2018-09-21 | 2019-01-22 | 国网电力科学研究院武汉南瑞有限责任公司 | A kind of emulated computation method of Transformer Winding Temperature Rise characteristic |
| CN114781886A (en) * | 2022-04-24 | 2022-07-22 | 广东电网有限责任公司 | Cable current-carrying capacity evaluation method and device, electronic equipment and storage medium |
-
2022
- 2022-12-12 CN CN202211587703.XA patent/CN115659765B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001351919A (en) * | 2000-06-05 | 2001-12-21 | Nec Corp | Wiring fault analysis method |
| JP2006250689A (en) * | 2005-03-10 | 2006-09-21 | Tokyo Electric Power Co Inc:The | Cable conductor temperature estimation method considering movement of air in cable tunnel, cable conductor temperature estimation system, and cable conductor temperature estimation program |
| CN106066212A (en) * | 2016-05-27 | 2016-11-02 | 三峡大学 | A kind of cable conductor temperature indirect measurement method |
| CN109255192A (en) * | 2018-09-21 | 2019-01-22 | 国网电力科学研究院武汉南瑞有限责任公司 | A kind of emulated computation method of Transformer Winding Temperature Rise characteristic |
| CN114781886A (en) * | 2022-04-24 | 2022-07-22 | 广东电网有限责任公司 | Cable current-carrying capacity evaluation method and device, electronic equipment and storage medium |
Non-Patent Citations (8)
| Title |
|---|
| JIANGJUN RUAN ET AL.: "《Real-time temperature estimation of three-core medium-voltage cable joint on SVR》" * |
| 傅晨钊,李清泉,王世山,李彦明: "电缆绕组变压器绕组温度场的二维数值计算" * |
| 刘刚;靳艳娇;马永强;孙立鹏;池骋;: "基于混合法的油浸式变压器二维瞬态温度场仿真" * |
| 周蠡;彭迎;李靖;陈然;张毅军;但林烽;黄雄峰;: "基于反演的电力电缆缆芯温度间接测量研究" * |
| 张森;章健;: "汽车通风盘式制动器的流固热多物理场耦合建模与分析" * |
| 李萌;牛胜锁;刘玉芹;赵可为;: "基于多场耦合模型的海底电缆载流量和温度场计算研究" * |
| 潘国友,杨乃伏,顾飞: "高炉煤气和炉料温度场迎风有限元法解析" * |
| 谢裕清;李琳;宋雅吾;王帅兵;: "油浸式电力变压器绕组温升的多物理场耦合计算方法" * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115659765B (en) | 2023-04-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bu et al. | Probabilistic analysis of small-signal stability of large-scale power systems as affected by penetration of wind generation | |
| CN106021189B (en) | A kind of more telegram in reply cable stable state temperature liter acquisition methods adapting to a variety of boundary conditions | |
| CN103279639B (en) | Receiving end Network Voltage Stability overall process Situation Assessment based on response and preventing control method | |
| Alves et al. | Continuation fast decoupled power flow with secant predictor | |
| Flatabo et al. | A method for calculation of margins to voltage instability applied on the Norwegian system for maintaining required security level | |
| CN108233400A (en) | A kind of more feed-in interaction factor computational methods of meter and hvdc control mode | |
| CN106021676B (en) | A kind of more telegram in reply cable stable state temperature rise acquisition methods based on transfer matrix | |
| CN109581178A (en) | Determine the method and apparatus of voltage source inverter IGBT module junction temperature | |
| CN107526872B (en) | Method for calculating thermal stress and deformation of 500kV ultrahigh-voltage cable | |
| CN107871037A (en) | A kind of 35kV dry reactors internal temperature evaluation method | |
| CN108956184A (en) | A kind of transformer forced oil air-cooled cooler radiating efficiency quantitative determination method and device | |
| CN111123041B (en) | Cable sheath fault positioning method based on temperature characteristic | |
| CN111241676A (en) | Modeling method, device and equipment based on current limiting type superconducting cable and storage medium | |
| CN115659765B (en) | Cable joint temperature field calculation method and device, electronic equipment and storage medium | |
| CN115713019A (en) | Method and device for calculating steady-state temperature rise of two-dimensional single-partition winding of power transformer | |
| CN117217052A (en) | Quick calculation method and related device for transient electric heating coupling field of converter transformer valve side sleeve | |
| CN109543246A (en) | A kind of method and system of determining transformer DC magnetic bias thermo parameters method | |
| Wang et al. | Research on digital twin modeling method of transformer temperature field based on POD | |
| CN114784804A (en) | Power system analysis method and device, electronic equipment and storage medium | |
| CN110705066B (en) | Dynamic simulation method for gas-electricity coupling park comprehensive energy system based on projection integration | |
| CN114091623A (en) | Transient voltage partitioning method and device based on multi-binary-table-AP clustering | |
| CN110969292A (en) | Topological graph-based power system state measuring and calculating method and device and electronic equipment | |
| CN107918697A (en) | The threephase cable transient state temperature rise acquisition methods of skin temperature comb laying are not depended on | |
| CN114679132B (en) | Photovoltaic inverter electric heating working condition simulation method and device and storage medium | |
| CN105870936A (en) | Probabilistic load flow-based SVC equipment location method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |