CN109408927B - Two-dimensional static magnetic field parallel finite element acceleration method based on black box transmission line model - Google Patents

Two-dimensional static magnetic field parallel finite element acceleration method based on black box transmission line model Download PDF

Info

Publication number
CN109408927B
CN109408927B CN201811192762.0A CN201811192762A CN109408927B CN 109408927 B CN109408927 B CN 109408927B CN 201811192762 A CN201811192762 A CN 201811192762A CN 109408927 B CN109408927 B CN 109408927B
Authority
CN
China
Prior art keywords
solving
iteration
finite element
nonlinear
transmission line
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.)
Active
Application number
CN201811192762.0A
Other languages
Chinese (zh)
Other versions
CN109408927A (en
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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201811192762.0A priority Critical patent/CN109408927B/en
Publication of CN109408927A publication Critical patent/CN109408927A/en
Application granted granted Critical
Publication of CN109408927B publication Critical patent/CN109408927B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)

Abstract

The invention discloses a finite element solving method of a two-dimensional nonlinear static magnetic field model based on a black box transmission line model, which comprises the following steps: firstly, determining variables to be solved and a solution domain; secondly, establishing an rz rectangular coordinate system; thirdly, listing a control equation and a boundary condition formula in the two-dimensional axisymmetric nonlinear static magnetic field and forming a differential equation set; fourthly, carrying out network division on the solution domain; fifthly, calculating a coefficient matrix of the finite element; sixthly, establishing an equivalent black box circuit model; seventhly, inserting a section of transmission line segment between the black box circuit model and the linear circuit; eighthly, iteration of a transmission line method is carried out; and ninthly, repeating iteration until the iteration result converges to a fixed error, and ending the solution. The method can solve and calculate the solved problem in parallel, and can be used for solving a complex finite element networking model, thereby solving the problems of long solving time and low efficiency caused by solving the finite element nonlinear problem by a Newton iteration method.

Description

Two-dimensional static magnetic field parallel finite element acceleration method based on black box transmission line model
Technical Field
The invention belongs to the technical field of electromagnetic field numerical calculation, relates to a finite element solving method of a two-dimensional nonlinear static magnetic field model, and particularly relates to a method for converting a finite element problem of a two-dimensional nonlinear static magnetic field into an equivalent black box circuit problem and realizing parallel accelerated solving of the finite element of the two-dimensional nonlinear static magnetic field by means of a transmission line iteration method.
Background
The finite element method is a numerical calculation method which is most widely applied in the industrial field and is adopted by a plurality of commercial simulation software. However, with the increasing complexity of the solution model and the increasing number of the sub-grid units, the nonlinear finite element solution method using the traditional newton iteration method as the core faces the problem of consuming time for solution, which is directly related to the speed and efficiency of the simulation development of the product.
The core of the finite element problem solution lies in solving a linear equation set, for a nonlinear problem, a new iteration result is used for regenerating a global matrix of a finite element model in each step of a traditional Newton iteration method, the dimension of the global matrix is continuously increased along with the continuous increase of the model sub-network, the time consumed by LU decomposition and the like of each step of the matrix is correspondingly increased, and the total solution time may be increased in a geometric form along with the densification of the sub-network. The transmission line iteration method can avoid repeated LU decomposition process and solving time. However, the conventional transmission line iteration method can only be used in finite element models of first-order triangle elements, and cannot be used in complex networking element models. Therefore, a new iteration method needs to be researched to improve the application range of the transmission line method and the capability of coping with a complex model, so that the problems of long solving time and low efficiency caused by solving the finite element nonlinearity problem by the Newton iteration method are solved.
Disclosure of Invention
The invention aims to provide a finite element solving method of a two-dimensional nonlinear static magnetic field model based on a black box transmission line model, which can solve and calculate the solved problem in parallel and can be used for solving a complex finite element networking model, wherein the complex finite element networking model comprises a triangular unit, a quadrilateral unit and the like, so that the problems of long solving time and low efficiency caused by solving the finite element nonlinear problem by a Newton iteration method are solved.
The purpose of the invention is realized by the following technical scheme:
a finite element solving method of a two-dimensional nonlinear static magnetic field model based on a black box transmission line model comprises the following steps:
firstly, determining variables to be solved and a solution domain, wherein: the variable to be solved is the magnetic potential A of a two-dimensional axisymmetric nonlinear static magnetic field, and the solution domain is the region where the two-dimensional axisymmetric nonlinear static magnetic field is located;
secondly, establishing a rz rectangular coordinate system;
thirdly, listing a control equation and a boundary condition formula in the two-dimensional axisymmetric nonlinear static magnetic field and forming a differential equation set, wherein:
the control equation is:
Figure GDA0001853978190000021
wherein J is a current density variable, μ ═ r μ, μ is the permeability of a triangular unit, a ═ rA, a is the magnetic potential, r is the abscissa, and z is the ordinate;
the boundary condition formula is as follows:
Γ1:A′=0;
Figure GDA0001853978190000031
in the formula, m is a normal vector on a solution boundary;
fourthly, adopting a net separating program to separate the solving domain, dispersing the solving domain into triangular units or quadrilateral units, and in each finite element unit, carrying out variable quantity
Figure GDA0001853978190000032
Wherein N is the node number of the finite element unit, NjIs a function of the shape in the cell,
Figure GDA0001853978190000033
the variable value of the corresponding node is the size;
fifthly, calculating a coefficient matrix of the finite element according to the following formula:
Figure GDA0001853978190000034
wherein the computational expression of each term is as follows:
Figure GDA0001853978190000035
Figure GDA0001853978190000036
sixthly, combining the matrix K in the step fiveeAdmittance matrix viewed as a circuit, beRegarding the current source vector as a current source vector, and establishing an equivalent black box circuit model;
seventhly, inserting a section of transmission line segment at any position of a connecting wire between the black box circuit model and the linear circuit;
eighthly, iteration of a transmission line method is carried out:
the voltage signal is incident into a linear network, and the circuit, namely an equation system is solved: (Y)linear+YTL)A′=b+2ViYTLThroughout the iteration, the matrix (Y)linear+YTL) Keeping unchanged, executing LU decomposition operation of the primary matrix in the 1 st step of iteration, and calculating the reflected voltage V after the solution is completed without executing the subsequent steps againr=A′-Vi
And (3) inputting the voltage signal into a nonlinear finite element unit, and solving the circuit, namely a nonlinear equation system:
Figure GDA0001853978190000041
and (4) solving by adopting a Newton iteration method, wherein the iteration formula of the k step is as follows:
Figure GDA0001853978190000042
wherein:
Figure GDA0001853978190000043
Figure GDA0001853978190000044
solving the equation set in each nonlinear unit is independently put into an independent computing core to be solved, parallel computing is realized, and after the solving is finished, the reflected voltage is calculated again
Figure GDA0001853978190000045
In the formula, YlinearAn admittance matrix that is a linear circuit incident within the linear network; y isTLAn admittance matrix that is a transmission line circuit incident within the linear network; a' is the voltage of the node to be solved; vrIs the reflected voltage reflected to the non-linear circuit; viIs the incident voltage incident to the linear circuit; keA coefficient matrix which is a nonlinear unit;
Figure GDA0001853978190000046
an admittance matrix that is reflected back to the transmission line within the nonlinear element;
Figure GDA0001853978190000047
is the reflected voltage reflected back into the nonlinear cell; j. the design is a squareeIs a Jacobian matrix;
Figure GDA0001853978190000048
is the Jacobian matrix at the k iteration; f is an intermediate variable in the calculation; a. thee′The value of the magnetic potential for each node in the cell.
And ninthly, repeating the iteration of the step eight until the iteration result converges to a fixed error, and ending the solution.
Compared with the prior art, the invention has the following advantages:
1. and (3) carrying out discrete modeling on the static magnetic field by adopting a finite element method, and enabling the element coefficient matrix of each element in the finite element to be equivalent to a sealed black box circuit model, thereby realizing the equivalent of a finite element equation set to be a nonlinear circuit network problem.
2. And a transmission line segment is inserted between the black box model of the nonlinear unit and the network of the linear unit, and transmission line iteration is realized by utilizing the iteration rule of the transmission line, so that the solution of the static magnetic field problem is realized.
3. In the transmission line iteration process, each nonlinear black box model is independent, and the solution calculation of the nonlinear black box model can adopt independent calculation resources (CPU or GPU core) for calculation, so that parallel acceleration is realized.
4. The iterative solution of finite elements can be carried out by adopting a transmission line method and a black box model, and various complex finite element models such as triangular units and quadrilateral units can be processed.
5. In the iterative solution process, the global matrix Y can be kept unchanged, in the matrix solution process, an LU decomposition method is adopted, LU decomposition is only needed to be carried out in the first calculation step, and because LU decomposition generally occupies about 95% of the time of matrix solution, the method can greatly reduce the iteration time of each step.
6. The nonlinear unit and the linear solving area can be isolated by using a transmission line method, and the parallel acceleration effect is realized.
Drawings
FIG. 1 is a black box circuit model;
FIG. 2 is a black box circuit network after insertion of transmission lines;
FIG. 3 is an equivalent circuit diagram after voltage is incident on the linear network;
FIG. 4 is an equivalent circuit diagram after voltage is reflected back to the nonlinear element;
FIG. 5 is a flow chart of an iterative solution;
FIG. 6 is a block diagram of a contactor of a certain type;
FIG. 7 is a cross-sectional view of a type of contactor;
FIG. 8 is a simplified model of a contactor;
FIG. 9 is a contactor calculation model;
FIG. 10 shows the result of the screening process, (a) triangular units, (b) quadrilateral units;
FIG. 11 is a diagram illustrating a magnetic field distribution in different states obtained by solving using a transmission line method;
FIG. 12 is a single step calculation time acceleration effect.
Detailed Description
The technical solutions of the present invention are further described below with reference to the drawings, but the present invention is not limited thereto, and modifications or equivalent substitutions may be made to the technical solutions of the present invention without departing from the scope of the technical solutions of the present invention.
The first embodiment is as follows: the embodiment provides a finite element solving method of a two-dimensional nonlinear static magnetic field model based on a black box transmission line model, which comprises the following steps:
the method comprises the steps of firstly, determining a variable to be solved and a solving domain, wherein the variable to be solved is a magnetic potential A of a two-dimensional axisymmetric nonlinear static magnetic field, the two-dimensional axisymmetric nonlinear static magnetic field is generated by current in an electrified coil, all elements around the electrified coil are made of ferromagnetic materials, and the solving domain is a region where the two-dimensional axisymmetric nonlinear static magnetic field is located.
And secondly, establishing a rz rectangular coordinate system, wherein r is a horizontal axis and z is a vertical axis.
Thirdly, listing a control equation and a boundary condition formula in the two-dimensional axisymmetric nonlinear static magnetic field and forming a differential equation set, wherein the control equation is as follows:
Figure GDA0001853978190000071
wherein J is a current density variable, μ ═ r μ, μ is a magnetic permeability of a triangular unit, a ═ rA, a is a magnetic potential, r is an abscissa, z is an ordinate, and boundary conditions are as follows: gamma-shaped1:A′=0;
Figure GDA0001853978190000072
And m is a normal vector on the solution boundary.
And fourthly, carrying out network division on the solving area by adopting a network division program, and dispersing the solving area into triangular units or quadrilateral units. In each finite element, the variable
Figure GDA0001853978190000073
Wherein N is the node number of the finite element unit, NjIs a function of the shape in the cell,
Figure GDA0001853978190000074
is the magnitude of the variable value of the corresponding node.
And fifthly, calculating a coefficient matrix of the finite element. According to the Galerkin method, carrying out weighted integration on a control equation to obtain an equation:
Figure GDA0001853978190000075
written in matrix form KeAe′=beNamely:
Figure GDA0001853978190000076
wherein the computational expression of each term is as follows:
Figure GDA0001853978190000077
Figure GDA0001853978190000078
sixthly, combining the matrix K in the step fiveeAdmittance matrix viewed as a circuit, beRegarded as a current source vector, an equivalent black box circuit model is established, as shown in fig. 1.
And seventhly, the black box model and the linear circuit are connected through a lead, and a transmission line segment is inserted into any position of the lead between the black box model and the linear circuit, as shown in fig. 2.
Eighthly, iteration of a transmission line method is carried out.
The voltage signal is incident on the linear network and, as shown in fig. 3, the circuit, i.e. the system of equations, is solved: (Y)linear+YTL)A′=b+2ViYTLThroughout the iteration, the matrix (Y)linear+YTL) Keeping unchanged, executing LU decomposition operation of the primary matrix in the 1 st step of iteration, and calculating the reflected voltage V after the solution is completed without executing the subsequent steps againr=A′-Vi
The voltage signal is incident on the nonlinear finite element, and as shown in fig. 4, the circuit, i.e. the nonlinear equation system, is solved:
Figure GDA0001853978190000081
the matrix is solved by adopting a Newton iteration method, and the iteration formula of the k step is as follows:
Figure GDA0001853978190000082
wherein:
Figure GDA0001853978190000083
Figure GDA0001853978190000084
because each nonlinear unit is isolated from each other, the solution of the equation set in each unit is independently put into an independent computing core for solution, and parallel computing is realized. After the solution is completed, the reflected voltage is calculated again
Figure GDA0001853978190000085
In the formula, YlinearIs the admittance matrix of the linear circuit in fig. 3; y isTLAn admittance matrix for the transmission line circuit within the dashed box of fig. 3; a' is the voltage of the node to be solved; vrIs the reflected voltage reflected to the non-linear circuit; viIs the incident voltage incident to the linear circuit; keA coefficient matrix which is a nonlinear unit;
Figure GDA0001853978190000091
an admittance matrix for the transmission lines of fig. 4;
Figure GDA0001853978190000092
is the reflected voltage reflected back into the cell; j. the design is a squareeIs a Jacobian matrix;
Figure GDA0001853978190000093
is the jacobian matrix at the kth iteration.
And continuously iterating until the iteration result converges to a fixed error, and ending the solution. The solving flow chart is shown in figure 5.
The second embodiment is as follows: the embodiment provides a finite element solving method of a two-dimensional nonlinear static magnetic field model based on a black box transmission line model, which comprises the following concrete implementation steps:
the method comprises the steps of firstly, determining a variable to be solved and a solving domain, wherein the variable to be solved is a magnetic potential A of a two-dimensional axisymmetric nonlinear static magnetic field, the two-dimensional axisymmetric nonlinear static magnetic field is generated by current in an electrified coil, all elements around the electrified coil are made of ferromagnetic materials, and the solving domain is a region where the two-dimensional axisymmetric nonlinear static magnetic field is located. In the present embodiment, the calculation is performed for a certain type of axially symmetric contactor, and the appearance view and the cross-sectional view of the model are shown in fig. 6 and 7, respectively. This is simplified and modeled, the simplified model being shown in fig. 8 and the computational model being shown in fig. 9.
And secondly, establishing a rz rectangular coordinate system, wherein r is a horizontal axis and z is a vertical axis.
Thirdly, listing a control equation and a boundary condition formula in the two-dimensional axisymmetric nonlinear static magnetic field and forming a differential equation set, wherein the control equation is as follows:
Figure GDA0001853978190000094
wherein J is a current density variable, μ ═ r μ, μ is a magnetic permeability of a triangular unit, a ═ rA, a is a magnetic potential, r is an abscissa, z is an ordinate, and boundary conditions are as follows: gamma-shaped1:A′=0;
Figure GDA0001853978190000095
And m is a normal vector on the solution boundary.
Fourthly, adopting a network division program to carry out network division on the solving areaThe solution domain is discretized into triangular units or quadrilateral units, and the result of the division is shown in fig. 10. In each finite element, the variable
Figure GDA0001853978190000101
Wherein N is the node number of the finite element unit, NjIs a function of the shape in the cell,
Figure GDA0001853978190000102
is the magnitude of the variable value of the corresponding node.
For the case of a triangular cell,
Figure GDA0001853978190000103
delta is the area of the triangular unit cell,
Figure GDA0001853978190000104
for a unit of a quadrilateral shape,
Figure GDA0001853978190000105
and fifthly, calculating a coefficient matrix of the finite element. According to the Galerkin method, carrying out weighted integration on a control equation to obtain an equation:
Figure GDA0001853978190000106
written in matrix form KeAe′=beNamely:
Figure GDA0001853978190000107
wherein the computational expression of each term is as follows:
Figure GDA0001853978190000108
Figure GDA0001853978190000111
sixthly, combining the matrix K in the step fiveeAdmittance matrix viewed as a circuit, beRegarded as a current source vector, an equivalent black box circuit model is established, as shown in fig. 1.
And seventhly, inserting a transmission line segment between the black box model and the linear circuit, as shown in figure 2.
Eighthly, iteration of a transmission line method is carried out.
The voltage signal is incident on the linear network and, as shown in fig. 3, the circuit, i.e. the system of equations, is solved: (Y)linear+YTL)A′=b+2ViYTLThen, the reflected voltage V is calculatedr=A′-Vi
The voltage signal is incident on the nonlinear finite element, and as shown in fig. 4, the circuit, i.e. the nonlinear equation system, is solved:
Figure GDA0001853978190000112
the matrix is solved by adopting a Newton iteration method, and the iteration formula of the k step is as follows:
Figure GDA0001853978190000113
wherein:
Figure GDA0001853978190000114
Figure GDA0001853978190000115
after the solution is completed, the reflected voltage is calculated again
Figure GDA0001853978190000116
And continuously iterating until the iteration result converges to a fixed error, and ending the solution. The solving flow chart is shown in figure 5.
The magnetic field distribution conditions in different states are obtained by solving with a transmission line method, and the result is shown in fig. 11.
TABLE 1 different solving models
Figure GDA0001853978190000121
case 1-5 is a triangular unit sub-net; case 6-10 is a quadrilateral unit sub-net.
TABLE 2 error comparison of N-R iteration and BB-TLM iteration of the present invention
Figure GDA0001853978190000122
Figure GDA0001853978190000123
No.1-5 is the error comparison of any five points in the case 5 model of Table 1; no.6-10 are error comparisons of five points arbitrarily taken for the case10 model of Table 1.
Comparing the single step calculation time of the newton iteration method and the transmission line method, as shown in fig. 12, the single step calculation time of the transmission line method is several times faster than that of the newton iteration method, compared with the single step calculation time of the newton iteration method and the transmission line method, which is accelerated by using 1, 4, 8, 16, 20 CPU cores in the middle case10 model in table 1.
Compared with the total calculation time of the newton iteration method and the transmission line method, the acceleration effect of the transmission line method is more obvious than that of the newton iteration method when the CPU cores are increased, in the case10 model in table 1, the CPU cores are accelerated by using 1, 4, 8, 16 and 20 respectively.
Table 3 solving the model in Table 1 by Newton's iteration method and transmission line method, respectively (CPU core 20)
Figure GDA0001853978190000131
As can be seen from table 3, when the grid division model is increased, the computation time of the transmission line method is faster than that of the newton iteration method.

Claims (1)

1. A finite element solving method of a two-dimensional nonlinear static magnetic field model based on a black box transmission line model is characterized by comprising the following steps:
firstly, determining variables to be solved and a solution domain, wherein: the variable to be solved is the magnetic potential A of a two-dimensional axisymmetric nonlinear static magnetic field, and the solution domain is the region where the two-dimensional axisymmetric nonlinear static magnetic field is located;
secondly, establishing a rz rectangular coordinate system;
thirdly, listing a control equation and a boundary condition formula in the two-dimensional axisymmetric nonlinear static magnetic field and forming a differential equation set, wherein:
the control equation is:
Figure FDA0003317730330000011
wherein J is a current density variable, μ ═ r μ, μ is the permeability of a triangular unit, a ═ rA, a is the magnetic potential, r is the abscissa, and z is the ordinate;
the boundary condition formula is as follows:
Γ1:A′=0;
Γ2:
Figure FDA0003317730330000012
in the formula, m is a normal vector on a solution boundary;
fourthly, adopting a net separating program to separate the solving domain, dispersing the solving domain into triangular units or quadrilateral units, and in each finite element unit, carrying out variable quantity
Figure FDA0003317730330000013
Wherein N is the node number of the finite element unit, NjIs a function of the shape in the cell,
Figure FDA0003317730330000014
the variable value of the corresponding node is the size;
fifthly, calculating a coefficient matrix of the finite element according to the following formula:
Figure FDA0003317730330000021
wherein the computational expression of each term is as follows:
Figure FDA0003317730330000022
Figure FDA0003317730330000023
sixthly, combining the matrix K in the step fiveeAdmittance matrix viewed as a circuit, beRegarding the current source vector as a current source vector, and establishing an equivalent black box circuit model;
seventhly, inserting a section of transmission line segment at any position of a connecting wire between the black box circuit model and the linear circuit;
eighthly, iteration of a transmission line method is carried out:
the voltage signal is incident into a linear network, and the circuit, namely an equation system is solved: (Y)linear+YTL)A′=b+2ViYTLThroughout the iteration, the matrix (Y)linear+YTL) Keeping unchanged, executing LU decomposition operation of the primary matrix in the 1 st step of iteration, and calculating the reflected voltage V after the solution is completed without executing the subsequent steps againr=A′-Vi
And (3) inputting the voltage signal into a nonlinear finite element unit, and solving the circuit, namely a nonlinear equation system:
Figure FDA0003317730330000024
and (4) solving by adopting a Newton iteration method, wherein the iteration formula of the k step is as follows:
Figure FDA0003317730330000025
solving the equation set in each nonlinear unit is independently put into an independent computing core to be solved, parallel computing is realized, and after the solving is finished, the incident voltage incident to the linear circuit is recalculated
Figure FDA0003317730330000026
In the formula, YlinearAn admittance matrix that is a linear circuit incident within the linear network; y isTLAn admittance matrix that is a transmission line circuit incident within the linear network; a' is the voltage of the node to be solved; vrIs the reflected voltage reflected to the non-linear circuit; viIs the incident voltage incident to the linear circuit; keA coefficient matrix which is a nonlinear unit;
Figure FDA0003317730330000031
an admittance matrix that is reflected back to the transmission line within the nonlinear element;
Figure FDA0003317730330000032
is the reflected voltage reflected back into the nonlinear cell;
Figure FDA0003317730330000033
is the Jacobian matrix at the k iteration;
and ninthly, repeating the iteration of the step eight until the iteration result converges to a fixed error, and ending the solution.
CN201811192762.0A 2018-10-13 2018-10-13 Two-dimensional static magnetic field parallel finite element acceleration method based on black box transmission line model Active CN109408927B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811192762.0A CN109408927B (en) 2018-10-13 2018-10-13 Two-dimensional static magnetic field parallel finite element acceleration method based on black box transmission line model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811192762.0A CN109408927B (en) 2018-10-13 2018-10-13 Two-dimensional static magnetic field parallel finite element acceleration method based on black box transmission line model

Publications (2)

Publication Number Publication Date
CN109408927A CN109408927A (en) 2019-03-01
CN109408927B true CN109408927B (en) 2022-01-18

Family

ID=65467151

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811192762.0A Active CN109408927B (en) 2018-10-13 2018-10-13 Two-dimensional static magnetic field parallel finite element acceleration method based on black box transmission line model

Country Status (1)

Country Link
CN (1) CN109408927B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112016232B (en) * 2020-08-31 2024-06-14 中国原子能科学研究院 Tear finite element process processing method and system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105631105A (en) * 2015-12-24 2016-06-01 北京航空航天大学 Cable equivalent circuit modeling method containing applied field coupling
CN106649939A (en) * 2016-09-28 2017-05-10 哈尔滨工业大学 Transmission line iteration-based solving method for 2D axial symmetric nonlinear magnetostatic field model
CN106991263A (en) * 2017-05-25 2017-07-28 清华大学 A kind of modeling method of mesolow distribution transformer winding thunder and lightning wideband model
CN107609274A (en) * 2017-09-14 2018-01-19 哈尔滨工业大学 Two-dimentional magnetostatic field parallel finite element method method based on transmission line Yu rank scheduling method

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004109560A1 (en) * 2003-06-05 2004-12-16 Fujitsu Limited Device and method for calculating electromagnetic field intensity by using circuit analysis and electromagnetic wave analysis
KR102066894B1 (en) * 2015-09-30 2020-01-16 오사카 유니버시티 Calculation method and program of electric potential, electric current, and peripheral electromagnetic field in electric circuit

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105631105A (en) * 2015-12-24 2016-06-01 北京航空航天大学 Cable equivalent circuit modeling method containing applied field coupling
CN106649939A (en) * 2016-09-28 2017-05-10 哈尔滨工业大学 Transmission line iteration-based solving method for 2D axial symmetric nonlinear magnetostatic field model
CN106991263A (en) * 2017-05-25 2017-07-28 清华大学 A kind of modeling method of mesolow distribution transformer winding thunder and lightning wideband model
CN107609274A (en) * 2017-09-14 2018-01-19 哈尔滨工业大学 Two-dimentional magnetostatic field parallel finite element method method based on transmission line Yu rank scheduling method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"Influence of electromagnetic pulse on electromagnetic relay by using equivalent circuit model";Li Bing等;《26th International Conference on Electrical Contacts (ICEC 2012)》;20120913;全文 *
"传输线理论及电磁兼容计算的半解析方法研究";张旭锋;《中国博士学位论文全文数据库 信息科技辑》;20120715;全文 *
"各向异性地层中可控源电磁法一维全参数反演及三维有限体积正演算法研究";周建美;《中国博士学位论文全文数据库 基础科学辑》;20140915;全文 *
"基于巨磁阻抗效应的弱磁传感器技术研究";张雷;《中国博士学位论文全文数据库 信息科技辑》;20160115;全文 *

Also Published As

Publication number Publication date
CN109408927A (en) 2019-03-01

Similar Documents

Publication Publication Date Title
CN110110413B (en) Structural topology optimization method based on material field reduction progression expansion
JP7133894B2 (en) Data-driven interactive 3D experience
Gyulassy et al. Efficient computation of Morse-Smale complexes for three-dimensional scalar functions
RU2519331C2 (en) Method of calculating physical value, numerical analysis method, programme for calculating physical value, numerical analysis programme, device for calculating physical value and numerical analysis device
CN107526860B (en) VLSI standard unit layout method based on electric field energy modeling technology
CN106682286B (en) Based on etc. geometry analysis methods Functional Gradient Materials modeling method
CN104361185A (en) Wire distribution space automatic generating method for cable virtual design
CN109408927B (en) Two-dimensional static magnetic field parallel finite element acceleration method based on black box transmission line model
WO2024007652A1 (en) Accelerated solving method for large sparse matrix, system, and storage medium
CN108763777B (en) Method for establishing VLSI global layout model based on Poisson equation explicit solution
CN111339688B (en) Method for solving rocket simulation model time domain equation based on big data parallel algorithm
MOON et al. Conservative treatment of boundary interfaces for overlaid grids andmulti-level grid adaptations
Zhou et al. Reduced-order ℒ2–ℒ∞ filtering for singular systems: a linear matrix inequality approach
Burrows et al. Learning overlap optimization for domain decomposition methods
Zhang et al. Adaptive finite element analysis of elliptic problems based on bubble-type local mesh generation
Mavriplis Viscous flow analysis using a parallel unstructured multigrid solver
US8676551B1 (en) Multi-solver simulation of dynamic systems in a modeling environment
CN104809297A (en) Electromagnetic force density transferring method used among special-shaped grids in magnetic field-structure field coupling calculation
Sonar Classical finite volume methods
CN113378279A (en) Building modeling method, building modeling system, electronic equipment and storage medium
CN107609274B (en) Two-dimensional static magnetic field parallel finite element method based on transmission line and level scheduling method
CN109584369B (en) Actual stratum full hexahedron grid generation method and device
CN106598911B (en) A kind of column write method of circuit network state equation
CN106202636B (en) A kind of method and system of micro-capacitance sensor electro-magnetic transient real-time simulation
CN107045488A (en) A kind of surface boundary member computational methods based on Bezier curved surface quadrangle parametric equations

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