CN112967401A - Three-dimensional space electric field calculation method based on finite element method - Google Patents

Three-dimensional space electric field calculation method based on finite element method Download PDF

Info

Publication number
CN112967401A
CN112967401A CN202110412305.3A CN202110412305A CN112967401A CN 112967401 A CN112967401 A CN 112967401A CN 202110412305 A CN202110412305 A CN 202110412305A CN 112967401 A CN112967401 A CN 112967401A
Authority
CN
China
Prior art keywords
electric field
calculating
equation
finite element
dimensional space
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.)
Pending
Application number
CN202110412305.3A
Other languages
Chinese (zh)
Inventor
李立峰
付炜平
梁利辉
池城
李秦
刘云鹏
何旺龄
韩笑
张胤禄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Super High Voltage Branch Of State Grid Hebei Electric Power Co ltd
State Grid Corp of China SGCC
Original Assignee
North China Electric Power University
Maintenance Branch of State Grid Hebei Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North China Electric Power University, Maintenance Branch of State Grid Hebei Electric Power Co Ltd filed Critical North China Electric Power University
Priority to CN202110412305.3A priority Critical patent/CN112967401A/en
Publication of CN112967401A publication Critical patent/CN112967401A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to the technical field of electric field research, in particular to a method for calculating a three-dimensional space electric field based on a finite element method. The method comprises the steps of carrying out three-dimensional modeling and space simulation on a three-dimensional space electric field, carrying out calculation analysis on an electric field boundary value, carrying out calculation analysis on a minimum electrostatic energy boundary value problem, subdividing a simulation model, carrying out calculation analysis on a unit energy value, carrying out multivariate quadratic function extreme value calculation, calculating an electric field intensity vector and the like. The invention designs a computer engineering software which adopts a finite element method as a principle to carry out simulation modeling of a three-dimensional space electric field, disperses the simulation electric field into a plurality of units to carry out energy calculation, obtains a more accurate approximate result of the electric field intensity by calculating the approximate solution of each unit, realizes the intensity analysis of a virtual electric field, and further not only can provide theoretical basis and support for the line structure design of a power transmission line, but also can provide a technical method for inhibiting the electrostatic induction and transient electric shock influence of the built line.

Description

Three-dimensional space electric field calculation method based on finite element method
Technical Field
The invention relates to the technical field of electric field research, in particular to a method for calculating a three-dimensional space electric field based on a finite element method.
Background
The electric shock phenomenon suffered by a human body in the area near the circuit can be divided into a steady-state electric shock and a transient-state electric shock according to different acting time, wherein the transient-state electric shock phenomenon is more common. Common transient electric shocks are easy to occur in metal structures and metal pipelines in the area near the line, and scenes such as clothes hangers, umbrella opening by residents, automobile parking and the like. Because no corresponding environmental protection standard evaluation exists, the treatment of the inductive electricity dispute event of the power transmission line close to the civil house is difficult and troublesome, and is difficult to properly solve. Therefore, the system is required to develop researches on evaluation and suppression measures of transient electric shocks below the overhead transmission line under different voltage classes and different scenes, reveal the electrostatic induction phenomenon and the action process of the electrostatic induction phenomenon with the human body, analyze the feeling characteristics of the human body suffering from the transient electric shocks under different electrostatic induction voltages, and provide measures for suppressing and alleviating the influence of the induction voltages and the transient electric shocks under different scenes on the production and life of residents. However, since the concept of the electric field is virtual, and there is a great risk in performing experimental studies in an actual electric field, the electric field strength can be studied by establishing a three-dimensional electric field. However, the three-dimensional space electric field simulated by adopting computer engineering software has the problem that the electric field intensity is difficult to accurately calculate and control. Therefore, if a more accurate approximate value of the electric field intensity can be calculated, a better research effect can be achieved. The results obtained from the finite element analysis calculation are only approximate, but if the number of the divided elements is very large and reasonable, the obtained results are consistent with the actual situation. At present, there is no perfect method for calculating and analyzing the three-dimensional space electric field by a finite ternary method.
Disclosure of Invention
The invention aims to provide a method for calculating a three-dimensional space electric field based on a finite element method, so as to solve the problems in the background technology.
In order to solve the above technical problem, an object of the present invention is to provide a method for calculating a three-dimensional space electric field based on a finite element method, including the steps of:
s1, selecting COMSOL software based on a finite element method as a principle, and performing three-dimensional modeling and space simulation on the three-dimensional space electric field;
s2, calculating and analyzing the electric field boundary value of the simulated three-dimensional electric field space;
s3, calculating and analyzing the minimum electrostatic energy boundary value problem;
s4, subdividing the simulation model, and calculating and analyzing the unit energy value;
s5, discretizing the functional problem, then performing multivariate quadratic function extremum calculation, and calculating the approximate value required by each node;
and S6, calculating an electric field intensity vector by combining the potential distributions.
As a further improvement of the present technical solution, in S2, the calculation expression of the electric field boundary value is:
Figure BDA0003024345370000021
Figure BDA0003024345370000022
in the electrostatic field, the electric field intensity vector and the electric displacement vector meet the loop theorem and the Gaussian theorem, rho is the charge density of the free body, and equation (1) and equation (2) are differential equations used for electrostatic field analysis in the Maxwell equation set;
will be provided with
Figure BDA0003024345370000023
Substituting formula (2) and considering D ═ E, one can obtain:
Figure BDA0003024345370000024
in the formula (3), ε is a dielectric constant, and an identity is used
Figure BDA0003024345370000025
Equation (3) can be expressed as:
Figure BDA0003024345370000026
the formula (4) is a function of potential
Figure BDA0003024345370000027
And bulk charge density p, dielectric constant epsilon is constant in a region of uniform medium,
Figure BDA0003024345370000028
equation (4) can be simplified to:
Figure BDA0003024345370000029
the expression of equation (5) in the rectangular coordinate system is:
Figure BDA00030243453700000210
in the region where ρ is 0, there are:
Figure BDA00030243453700000211
equation (6) is called the poisson equation of the electrostatic field, and equation (7) is called the laplace equation of the electrostatic field, where the laplace equation is a special case of the poisson equation.
As a further improvement of the present technical solution, in S3, the minimum electrostatic energy margin has three types of problems, and the expression thereof is:
the electrostatic field minimum action principle is adopted, a fixed belt system is positioned in a medium, the surface charge distribution of the fixed belt system enables the synthesized electric field to have minimum electrostatic energy, and the expression of the first-class boundary value problem energy integral is as follows:
Figure BDA0003024345370000031
the electrostatic field problem is therefore equivalent to the functional extremum problem:
Figure BDA0003024345370000032
similarly, the equivalence functional extreme value problems of the second and third types of edge value problems are as follows:
Figure BDA0003024345370000033
therein, functional
Figure BDA0003024345370000034
Quadratic dependence function
Figure BDA0003024345370000035
And partial derivatives thereof, hence the term
Figure BDA0003024345370000036
Is a function of
Figure BDA0003024345370000037
The corresponding variation problem is called a quadratic functional extreme value problem.
As a further improvement of the present technical solution, in S4, the method for performing calculation analysis on the unit energy value includes the following steps:
s4.1, performing secondary functional calculation of all unit energies;
s4.2, sequentially carrying out energy integral calculation of the single unit;
s4.3, calculating unit energy integrals which are independent of x, y and z respectively;
and S4.4, simultaneously combining the equations of all the sub-regions.
As a further improvement of the present technical solution, in S4.1, performing a quadratic functional calculation on all unit energies may be expressed as a sum of all unit energy functions:
Figure BDA0003024345370000041
as a further improvement of the present technical solution, in S4.2, an expression for performing single-unit energy integral calculation is as follows:
Figure BDA0003024345370000042
and the following steps:
Figure BDA0003024345370000043
thus:
Figure BDA0003024345370000044
namely:
Figure BDA0003024345370000045
as a further improvement of the present technical solution, in S4.3, expressions for calculating the unit energy integrals independent of x, independent of y, and independent of z are respectively:
if it is
Figure BDA0003024345370000046
Independent of x, then:
Figure BDA0003024345370000047
wherein:
Figure BDA0003024345370000048
Figure BDA0003024345370000051
wherein S isεThe area enclosed by the subdivision area;
in the same way, if
Figure BDA0003024345370000052
Independently of y, then:
Figure BDA0003024345370000053
Figure BDA0003024345370000054
wherein S isεThe area enclosed by the subdivision area;
in the same way, if
Figure BDA0003024345370000055
Independent of z, then:
Figure BDA0003024345370000056
Figure BDA0003024345370000057
wherein S isεThe area enclosed by the subdivision area.
As a further improvement of the present technical solution, in S4.4, an expression obtained by combining equations in a simultaneous manner is as follows:
from equation (12) -equation (22), it can be derived:
Figure BDA0003024345370000058
by combining equations (23) for all sub-regions concurrently, equation (11) can be written as:
Figure BDA0003024345370000061
wherein
Figure BDA0003024345370000062
And K is the matrix of variables and coefficients after all the equation sets are connected.
As a further improvement of the present technical solution, in S5, the expression of discretizing the functional problem and then performing multivariate quadratic function extremum calculation is as follows:
Figure BDA0003024345370000063
from the theory of extreme values of functions, there are:
Figure BDA0003024345370000064
the formula is as follows:
Figure BDA0003024345370000065
wherein, the formula (27) is a finite element equation, and the equation (27) is solved in sequence to obtain the node of each node
Figure BDA0003024345370000066
The values may be used to obtain the desired approximate solution.
As a further improvement of the present technical solution, in S6, the expression for calculating the electric field strength vector is as follows:
Figure BDA0003024345370000067
the electric field intensity vector E of each potential can be obtained by combining the potential distribution.
Another object of the present invention is to provide an apparatus for calculating a three-dimensional space electric field based on a finite element method, including a processor, a memory, and a computer program stored in the memory and executed on the processor, wherein the processor is configured to implement any of the steps of the method for calculating a three-dimensional space electric field based on a finite element method when the computer program is executed.
It is a further object of the present invention to provide a computer-readable storage medium storing a computer program, which when executed by a processor, implements the steps of any of the above-described methods for calculating a three-dimensional space electric field based on a finite element method.
Compared with the prior art, the invention has the beneficial effects that: in the method for calculating the three-dimensional space electric field based on the finite element method, computer engineering software which takes the finite element method as a principle is adopted to carry out simulation modeling on the three-dimensional space electric field, the finite element method is adopted to calculate the boundary value of the simulation electric field and the minimum electrostatic energy of the electric field, the simulation electric field is dispersed into a plurality of units to carry out energy calculation, and an accurate electric field strength approximate result is obtained by calculating the approximate solution of each unit to realize the intensity analysis of the virtual electric field, so that not only can theoretical basis and support be provided for the circuit structure design of the power transmission line, but also a technical method can be provided for the electrostatic induction and transient electric shock influence suppression of the built circuit, and the method has great significance for building an environment-friendly power grid.
Drawings
FIG. 1 is an overall process flow diagram of the present invention;
FIG. 2 is a flow chart of a partial method of the present invention;
FIG. 3 is a block diagram of an exemplary computer program product of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments 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.
Method embodiment
As shown in fig. 1 to fig. 3, the present embodiment aims to provide a method for calculating a three-dimensional space electric field based on a finite element method, which includes the following steps:
s1, selecting COMSOL software based on a finite element method as a principle, and performing three-dimensional modeling and space simulation on the three-dimensional space electric field;
s2, calculating and analyzing the electric field boundary value of the simulated three-dimensional electric field space;
s3, calculating and analyzing the minimum electrostatic energy boundary value problem;
s4, subdividing the simulation model, and calculating and analyzing the unit energy value;
s5, discretizing the functional problem, then performing multivariate quadratic function extremum calculation, and calculating the approximate value required by each node;
and S6, calculating an electric field intensity vector by combining the potential distributions.
Specifically, the finite element method is a numerical calculation method for solving mathematical equations, is a numerical analysis technology organically combining elastic theory, computational mathematics and computer software, and is a powerful numerical calculation tool for solving the practical engineering problem.
The finite element method is expressed in a matrix form, and is convenient for compiling a computer program. At present, computer engineering software which can be conveniently and directly applied at home and abroad comprises MSC NASTRAN, ADINA, LS-DYNA, ANSYS, ANSOFT, ABAQUS, 2D-sigma, COMSOL and the like, wherein the COMSOL and ANSYS are both suitable for electric field analysis.
In this embodiment, in S2, the calculation expression of the electric field boundary value is:
Figure BDA0003024345370000081
Figure BDA0003024345370000082
in the electrostatic field, the electric field intensity vector and the electric displacement vector meet the loop theorem and the Gaussian theorem, rho is the charge density of the free body, and equation (1) and equation (2) are differential equations used for electrostatic field analysis in the Maxwell equation set;
will be provided with
Figure BDA0003024345370000083
Substituting formula (2) and considering D ═ E, one can obtain:
Figure BDA0003024345370000084
in the formula (3), ε is a dielectric constant, and an identity is used
Figure BDA0003024345370000085
Equation (3) can be expressed as:
Figure BDA0003024345370000086
the formula (4) is a function of potential
Figure BDA0003024345370000087
And bulk charge density p, dielectric constant epsilon is constant in a region of uniform medium,
Figure BDA0003024345370000088
equation (4) can be simplified to:
Figure BDA0003024345370000089
the expression of equation (5) in the rectangular coordinate system is:
Figure BDA00030243453700000810
in the region where ρ is 0, there are:
Figure BDA00030243453700000811
equation (6) is called the poisson equation of the electrostatic field, and equation (7) is called the laplace equation of the electrostatic field, where the laplace equation is a special case of the poisson equation.
In this embodiment, in S3, the minimum electrostatic energy margin has three types of problems, and the expression thereof is:
the electrostatic field minimum action principle is adopted, a fixed belt system is positioned in a medium, the surface charge distribution of the fixed belt system enables the synthesized electric field to have minimum electrostatic energy, and the expression of the first-class boundary value problem energy integral is as follows:
Figure BDA0003024345370000091
the electrostatic field problem is therefore equivalent to the functional extremum problem:
Figure BDA0003024345370000092
similarly, the equivalence functional extreme value problems of the second and third types of edge value problems are as follows:
Figure BDA0003024345370000093
therein, functional
Figure BDA0003024345370000094
Quadratic dependence function
Figure BDA0003024345370000095
And partial derivatives thereof, hence the term
Figure BDA0003024345370000096
Is a function of
Figure BDA0003024345370000097
The corresponding variation problem is called a quadratic functional extreme value problem.
Specifically, the second and third types of boundary condition problems are included in the variation among the requirements for the functional to reach the extremum. And the boundary conditions on the interfaces of different media are also included in the requirement that the functional reaches the extreme value and are automatically met. These conditions are therefore referred to as natural boundary conditions. For the first class of boundary conditions, which need to be given as definite conditions in the variational problem, such boundary conditions are called imposed boundary conditions.
In this embodiment, in S4, the method for performing calculation analysis on the unit energy value includes the following steps:
s4.1, performing secondary functional calculation of all unit energies;
s4.2, sequentially carrying out energy integral calculation of the single unit;
s4.3, calculating unit energy integrals which are independent of x, y and z respectively;
and S4.4, simultaneously combining the equations of all the sub-regions.
In this embodiment, in S4.1, performing the quadratic functional calculation on all unit energies may be expressed as a sum of all unit energy functions:
Figure BDA0003024345370000101
in this embodiment, in S4.2, the expression for performing the single-unit energy integral calculation is as follows:
Figure BDA0003024345370000102
and the following steps:
Figure BDA0003024345370000103
thus:
Figure BDA0003024345370000104
namely:
Figure BDA0003024345370000105
in this embodiment, in S4.3, the expressions for the integral calculation of the unit energy independent of x, independent of y, and independent of z are:
if it is
Figure BDA0003024345370000106
Independent of x, then:
Figure BDA0003024345370000107
wherein:
Figure BDA0003024345370000108
Figure BDA0003024345370000111
wherein S isεThe area enclosed by the subdivision area;
in the same way, if
Figure BDA0003024345370000112
Independently of y, then:
Figure BDA0003024345370000113
Figure BDA0003024345370000114
wherein S isεThe area enclosed by the subdivision area;
in the same way, if
Figure BDA0003024345370000115
Independent of z, then:
Figure BDA0003024345370000116
Figure BDA0003024345370000117
wherein S isεThe area enclosed by the subdivision area.
In this embodiment, in S4.4, the expression for simultaneous combination of equations is:
from equation (12) -equation (22), it can be derived:
Figure BDA0003024345370000118
by combining equations (23) for all sub-regions concurrently, equation (11) can be written as:
Figure BDA0003024345370000121
wherein
Figure BDA0003024345370000122
And K is the matrix of variables and coefficients after all the equation sets are connected.
In this embodiment, in S5, the expression of discretizing the functional problem and then performing multivariate quadratic function extremum calculation is:
Figure BDA0003024345370000123
from the theory of extreme values of functions, there are:
Figure BDA0003024345370000124
the formula is as follows:
Figure BDA0003024345370000125
wherein, the formula (27) is a finite element equation, and the equation (27) is solved in sequence to obtain the node of each node
Figure BDA0003024345370000126
The values may be used to obtain the desired approximate solution.
In this embodiment, in S6, the expression for calculating the electric field strength vector is:
Figure BDA0003024345370000127
the electric field intensity vector E of each potential can be obtained by combining the potential distribution.
Computer program product embodiment
Referring to fig. 3, a schematic structural diagram of a device for calculating a three-dimensional space electric field based on the finite element method according to the present embodiment is shown, where the device includes a processor, a memory, and a computer program stored in the memory and running on the processor.
The processor comprises one or more than one processing core, the processor is connected with the processor through a bus, the memory is used for storing program instructions, and the method for calculating the three-dimensional space electric field based on the finite element method is realized when the processor executes the program instructions in the memory.
Alternatively, the memory may be implemented by any type or combination of volatile or non-volatile memory devices, such as Static Random Access Memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic or optical disks.
In addition, the present invention also provides a computer-readable storage medium, in which a computer program is stored, and the computer program, when executed by a processor, implements the steps of the method for calculating a three-dimensional space electric field based on the finite element method.
Optionally, the present invention further provides a computer program product containing instructions, which when run on a computer, causes the computer to perform the steps of the method for calculating a three-dimensional space electric field based on the finite element method according to the above aspects.
It will be understood by those skilled in the art that all or part of the steps for implementing the above embodiments may be implemented by hardware, or may be implemented by hardware related to instructions of a program, which may be stored in a computer-readable storage medium, such as a read-only memory, a magnetic or optical disk, and the like.
The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, and the preferred embodiments of the present invention are described in the above embodiments and the description, and are not intended to limit the present invention. The scope of the invention is defined by the appended claims and equivalents thereof.

Claims (10)

1. A three-dimensional space electric field calculation method based on a finite element method is characterized in that: the method comprises the following steps:
s1, selecting COMSOL software based on a finite element method as a principle, and performing three-dimensional modeling and space simulation on the three-dimensional space electric field;
s2, calculating and analyzing the electric field boundary value of the simulated three-dimensional electric field space;
s3, calculating and analyzing the minimum electrostatic energy boundary value problem;
s4, subdividing the simulation model, and calculating and analyzing the unit energy value;
s5, discretizing the functional problem, then performing multivariate quadratic function extremum calculation, and calculating the approximate value required by each node;
and S6, calculating an electric field intensity vector by combining the potential distributions.
2. The finite element method-based method for calculating the electric field in the three-dimensional space according to claim 1, wherein: in S2, the calculation expression of the electric field boundary value is:
Figure FDA0003024345360000011
Figure FDA0003024345360000012
in the electrostatic field, the electric field intensity vector and the electric displacement vector meet the loop theorem and the Gaussian theorem, rho is the charge density of the free body, and equation (1) and equation (2) are differential equations used for electrostatic field analysis in the Maxwell equation set;
will be provided with
Figure FDA0003024345360000013
Substituting formula (2) and considering D ═ E, one can obtain:
Figure FDA0003024345360000014
in the formula (3), ε is a dielectric constant, and an identity is used
Figure FDA0003024345360000015
Equation (3) can be expressed as:
Figure FDA0003024345360000016
the formula (4) is a function of potential
Figure FDA0003024345360000017
And bulk charge density p, dielectric constant epsilon is constant in a region of uniform medium,
Figure FDA0003024345360000018
equation (4) can be simplified to:
Figure FDA0003024345360000019
the expression of equation (5) in the rectangular coordinate system is:
Figure FDA00030243453600000110
in the region where ρ is 0, there are:
Figure FDA0003024345360000021
equation (6) is called the poisson equation of the electrostatic field, and equation (7) is called the laplace equation of the electrostatic field, where the laplace equation is a special case of the poisson equation.
3. The finite element method-based method for calculating the electric field in the three-dimensional space according to claim 2, wherein: in S3, the minimum electrostatic energy margin has three types of problems, and the expression is:
the electrostatic field minimum action principle is adopted, a fixed belt system is positioned in a medium, the surface charge distribution of the fixed belt system enables the synthesized electric field to have minimum electrostatic energy, and the expression of the first-class boundary value problem energy integral is as follows:
Figure FDA0003024345360000022
the electrostatic field problem is therefore equivalent to the functional extremum problem:
Figure FDA0003024345360000023
similarly, the equivalence functional extreme value problems of the second and third types of edge value problems are as follows:
Figure FDA0003024345360000024
therein, functional
Figure FDA0003024345360000028
Quadratic dependence function
Figure FDA0003024345360000025
And partial derivatives thereof, hence the term
Figure FDA0003024345360000026
Is a function of
Figure FDA0003024345360000027
The corresponding variation problem is called a quadratic functional extreme value problem.
4. The finite element method-based method for calculating the electric field in the three-dimensional space according to claim 3, wherein: in S4, the method for performing calculation analysis on the unit energy value includes the following steps:
s4.1, performing secondary functional calculation of all unit energies;
s4.2, sequentially carrying out energy integral calculation of the single unit;
s4.3, calculating unit energy integrals which are independent of x, y and z respectively;
and S4.4, simultaneously combining the equations of all the sub-regions.
5. The finite element method-based method for calculating the electric field in the three-dimensional space according to claim 4, wherein: in S4.1, performing a quadratic functional calculation on all unit energies may be expressed as a sum of all unit energy functions:
Figure FDA0003024345360000031
6. the finite element method-based method for calculating the electric field in the three-dimensional space according to claim 5, wherein: in S4.2, the expression for performing the single-unit energy integral calculation is as follows:
Figure FDA0003024345360000032
and the following steps:
Figure FDA0003024345360000033
thus:
Figure FDA0003024345360000034
namely:
Figure FDA0003024345360000035
7. the finite element method-based method for calculating the electric field in the three-dimensional space according to claim 6, wherein: in S4.3, the expressions for the integral calculation of the unit energy independent of x, independent of y, and independent of z are:
if it is
Figure FDA0003024345360000036
Independent of x, then:
Figure FDA0003024345360000037
wherein:
Figure FDA0003024345360000038
Figure FDA0003024345360000041
wherein S isεThe area enclosed by the subdivision area;
in the same way, if
Figure FDA0003024345360000042
Independently of y, then:
Figure FDA0003024345360000043
Figure FDA0003024345360000044
wherein S isεThe area enclosed by the subdivision area;
in the same way, if
Figure FDA0003024345360000045
Independent of z, then:
Figure FDA0003024345360000046
Figure FDA0003024345360000047
wherein S isεThe area enclosed by the subdivision area.
8. The finite element method-based method for calculating the electric field in the three-dimensional space according to claim 7, wherein: in S4.4, the expression for simultaneous combination of equations is:
from equation (12) -equation (22), it can be derived:
Figure FDA0003024345360000051
by combining equations (23) for all sub-regions concurrently, equation (11) can be written as:
Figure FDA0003024345360000052
wherein
Figure FDA0003024345360000053
And K is the matrix of variables and coefficients after all the equation sets are connected.
9. The finite element method-based method for calculating the electric field in the three-dimensional space according to claim 8, wherein: in S5, the expression of discretizing the functional problem and then performing multivariate quadratic function extremum calculation is:
Figure FDA0003024345360000054
from the theory of extreme values of functions, there are:
Figure FDA0003024345360000055
the formula is as follows:
Figure FDA0003024345360000056
wherein, the formula (27) is a finite element equation, and the equation (27) is solved in sequence to obtain the node of each node
Figure FDA0003024345360000057
The values may be used to obtain the desired approximate solution.
10. The finite element method-based method for calculating the electric field in the three-dimensional space according to claim 9, wherein: in S6, the expression for calculating the electric field strength vector is:
Figure FDA0003024345360000058
the electric field intensity vector E of each potential can be obtained by combining the potential distribution.
CN202110412305.3A 2021-04-16 2021-04-16 Three-dimensional space electric field calculation method based on finite element method Pending CN112967401A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110412305.3A CN112967401A (en) 2021-04-16 2021-04-16 Three-dimensional space electric field calculation method based on finite element method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110412305.3A CN112967401A (en) 2021-04-16 2021-04-16 Three-dimensional space electric field calculation method based on finite element method

Publications (1)

Publication Number Publication Date
CN112967401A true CN112967401A (en) 2021-06-15

Family

ID=76280625

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110412305.3A Pending CN112967401A (en) 2021-04-16 2021-04-16 Three-dimensional space electric field calculation method based on finite element method

Country Status (1)

Country Link
CN (1) CN112967401A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114048589A (en) * 2021-10-19 2022-02-15 中国南方电网有限责任公司超高压输电公司检修试验中心 Electric field determining method and device, electronic equipment and storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412989A (en) * 2013-08-01 2013-11-27 电子科技大学 Parameterized reduced model based three-dimensional electromagnetic field simulation method of periodic structure
CN104834765A (en) * 2014-12-31 2015-08-12 国网电力科学研究院武汉南瑞有限责任公司 Power cable accessory state evaluation method based on finite element electric field analysis

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412989A (en) * 2013-08-01 2013-11-27 电子科技大学 Parameterized reduced model based three-dimensional electromagnetic field simulation method of periodic structure
CN104834765A (en) * 2014-12-31 2015-08-12 国网电力科学研究院武汉南瑞有限责任公司 Power cable accessory state evaluation method based on finite element electric field analysis

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
汪茹 等: "一种基于有限元法的超高压输电线路三维电场计算方法", 电子测量与仪器学报, no. 04, 30 April 2019 (2019-04-30), pages 103 - 110 *
潘静 等: "装设辅助设备输电线路软梯法进等电位电场仿真计算与分析", 电气技术, no. 10, 31 October 2017 (2017-10-31), pages 21 - 26 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114048589A (en) * 2021-10-19 2022-02-15 中国南方电网有限责任公司超高压输电公司检修试验中心 Electric field determining method and device, electronic equipment and storage medium

Similar Documents

Publication Publication Date Title
Rao et al. Estimating the saddle-node bifurcation point of static power systems using the holomorphic embedding method
Clemens et al. Discrete electromagnetism with the finite integration technique
Aristidou et al. Power system dynamic simulations using a parallel two-level schur-complement decomposition
Schöps Multiscale modeling and multirate time-integration of field/circuit coupled problems
De Dios et al. Discontinuous Galerkin methods for the multi-dimensional Vlasov–Poisson problem
Qi et al. Nonlinear model reduction in power systems by balancing of empirical controllability and observability covariances
CN112967401A (en) Three-dimensional space electric field calculation method based on finite element method
Yang et al. A decoupled time-domain simulation method via invariant subspace partition for power system analysis
Villa et al. An implicit three-dimensional fractional step method for the simulation of the corona phenomenon
Qiao et al. An iterative flux tracing method without Deutsch assumption for ion-flow field of AC/DC hybrid transmission lines
CN112924763B (en) Method for evaluating transient electric shock below high-voltage alternating-current transmission line
Liu et al. Faster than real-time simulation: Methods, tools, and applications
Guizán et al. Numerical analysis and safety design of grounding systems in underground compact substations
Amoruso et al. Thundercloud pre-stroke electrostatic modeling
Averkin et al. A parallel electrostatic Particle-in-Cell method on unstructured tetrahedral grids for large-scale bounded collisionless plasma simulations
Gao et al. Element differential method for solving linear and nonlinear electromagnetic problems
Kochetov et al. Stability of full-wave PEEC models: Reason for instabilities and approach for correction
Poursalehi et al. An adaptive mesh refinement approach for average current nodal expansion method in 2-D rectangular geometry
Zhao et al. GPU based parallel matrix exponential algorithm for large scale power system electromagnetic transient simulation
Aristidou et al. Algorithmic and computational advances for fast power system dynamic simulations
Dong et al. Parallel numerical computing of finite element model of conductors and floating potentials
Zhang et al. Parallel Simulation of Power Systems with High Penetration of Distributed Generation Using GPUs and OpenCL
Keßler et al. Vlasov--Poisson System Tackled by Particle Simulation Utilizing Boundary Element Methods
Marcsa et al. Finite element tearing and interconnecting method and its algorithms for parallel solution of magnetic field problems
Aiello et al. Efficient analysis of grounding systems by means of the hybrid FEM–DBCI 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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 050000 66 Zhongsheng Road, Xinhua District, Shijiazhuang, Hebei Province

Applicant after: Super high voltage branch of State Grid Hebei Electric Power Co.,Ltd.

Applicant after: NORTH CHINA ELECTRIC POWER University (BAODING)

Address before: 050022 No.66, Zhongsheng Road, Xinhua District, Shijiazhuang City, Hebei Province

Applicant before: MAINTENANCE BRANCH COMPANY OF STATE GRID HEBEI ELECTRIC POWER Co.

Applicant before: NORTH CHINA ELECTRIC POWER University (BAODING)

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20220922

Address after: 050000 66 Zhongsheng Road, Xinhua District, Shijiazhuang, Hebei Province

Applicant after: Super high voltage branch of State Grid Hebei Electric Power Co.,Ltd.

Applicant after: STATE GRID CORPORATION OF CHINA

Address before: 050000 66 Zhongsheng Road, Xinhua District, Shijiazhuang, Hebei Province

Applicant before: Super high voltage branch of State Grid Hebei Electric Power Co.,Ltd.

Applicant before: NORTH CHINA ELECTRIC POWER University (BAODING)