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 PDFInfo
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design 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
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:
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;
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 quantityWherein N is the node number of the finite element unit, NjIs a function of the shape in the cell,the variable value of the corresponding node is the size;
fifthly, calculating a coefficient matrix of the finite element according to the following formula:
wherein the computational expression of each term is as follows:
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:and (4) solving by adopting a Newton iteration method, wherein the iteration formula of the k step is as follows:wherein:
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
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;an admittance matrix that is reflected back to the transmission line within the nonlinear element;is the reflected voltage reflected back into the nonlinear cell; j. the design is a squareeIs a Jacobian matrix;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: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;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 variableWherein N is the node number of the finite element unit, NjIs a function of the shape in the cell,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:
written in matrix form KeAe′=beNamely:
wherein the computational expression of each term is as follows:
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:the matrix is solved by adopting a Newton iteration method, and the iteration formula of the k step is as follows:wherein:
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
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;an admittance matrix for the transmission lines of fig. 4;is the reflected voltage reflected back into the cell; j. the design is a squareeIs a Jacobian matrix;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: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;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 variableWherein N is the node number of the finite element unit, NjIs a function of the shape in the cell,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:
written in matrix form KeAe′=beNamely:
wherein the computational expression of each term is as follows:
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:the matrix is solved by adopting a Newton iteration method, and the iteration formula of the k step is as follows:wherein:
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
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
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)
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:
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;
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 quantityWherein N is the node number of the finite element unit, NjIs a function of the shape in the cell,the variable value of the corresponding node is the size;
fifthly, calculating a coefficient matrix of the finite element according to the following formula:
wherein the computational expression of each term is as follows:
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:and (4) solving by adopting a Newton iteration method, wherein the iteration formula of the k step is as follows: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
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;an admittance matrix that is reflected back to the transmission line within the nonlinear element;is the reflected voltage reflected back into the nonlinear cell;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.
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)
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)
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)
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 |
-
2018
- 2018-10-13 CN CN201811192762.0A patent/CN109408927B/en active Active
Patent Citations (4)
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)
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 |