CN109408927A - A kind of two-dimentional magnetostatic field parallel finite element method accelerated method based on black box transmission line model - Google Patents
A kind of two-dimentional magnetostatic field parallel finite element method accelerated method based on black box transmission line model Download PDFInfo
- Publication number
- CN109408927A CN109408927A CN201811192762.0A CN201811192762A CN109408927A CN 109408927 A CN109408927 A CN 109408927A CN 201811192762 A CN201811192762 A CN 201811192762A CN 109408927 A CN109408927 A CN 109408927A
- Authority
- CN
- China
- Prior art keywords
- linear
- finite element
- transmission line
- circuit
- iteration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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 kind of finite element solving methods of two-dimension non linearity magnetostatic model based on black box transmission line model, and described method includes following steps: one, determining variable to be solved and solve domain;Two, rz rectangular coordinate system is established;Three, governing equation in two-dimensional axial symmetric nonlinear magnetostatic field and boundary condition formula and one differential equation group of composition are listed;Four, subnetting is carried out to solution domain;Five, the coefficient matrix of finite element is calculated;Six, equivalent black box circuit model is established;Seven, one section of transmission line line segment is inserted between black box circuit model and linear circuit;Eight, Transmission line method iteration is carried out;Nine, repeat iteration, until iteration result converges to fixed error, end is solved.The present invention can concurrently carry out solution calculating to required problem, and can be used in complicated finite element meshing model solution, to solve the problems, such as that bring solves time length, low efficiency when Newton iteration method solving finite element nonlinear problem.
Description
Technical field
The invention belongs to Numerical Calculation of Electromagnetic Fields technical field, it is related to a kind of the limited of two-dimension non linearity magnetostatic model
First method for solving, and in particular to a kind of Finite Element Equations by two-dimension non linearity magnetostatic field are converted into black box circuit of equal value and ask
Topic, and realized by transmission line iteration methods and the method solved is accelerated parallel to the finite element of two-dimension non linearity magnetostatic field.
Background technique
FInite Element is most widely used numerical computation method in industrial circle, is used by many commercial simulation softwares.
However, being with traditional Newton iteration method with increasingly complicating for solving model and being increasing for subnetting number of unit
The nonlinear finite element method for solving of core, which is faced with, solves time-consuming serious problem, this emulation for being directly related to product is ground
The speed and efficiency of hair.
The core of the solution of Finite Element Equations is solution system of linear equations, and for nonlinear problem, it is traditional
The each step of Newton iteration method is employed to the global matrix that new iteration result regenerates finite element model, with model point
The time for the consumption such as the continuous increase of net, the dimension of global matrix constantly become larger, and the LU of each step matrix is decomposed can increase accordingly
Greatly, the overall solution time may form geometric increase as the change of subnetting is close.It can be avoided using transmission line iteration methods
Duplicate LU decomposable process solves the time.However, traditional transmission line iteration methods is only used for having for single order triangular element
It limits in meta-model, is not used to complicated subnetting model of element.Therefore, it is necessary to study a kind of new alternative manner, Lai Tigao
The ability of the scope of application of Transmission line method and reply complex model, to solve that Newton iteration method solving finite element is non-linear to ask
Bring solves the problem of long time, low efficiency when topic.
Summary of the invention
The object of the present invention is to provide a kind of the limited of two-dimension non linearity magnetostatic model based on black box transmission line model
First method for solving, this method can concurrently carry out solution calculating to required problem, and can be used in complicated finite element fraction
In pessimistic concurrency control solves, including triangular element and quadrilateral units etc., so that it is non-thread to solve Newton iteration method solving finite element
Property problem when bring the problem of solving time long, low efficiency.
The purpose of the present invention is what is be achieved through the following technical solutions:
A kind of finite element solving method of the two-dimension non linearity magnetostatic model based on black box transmission line model, including such as
Lower step:
One, it determines variable to be solved and solves domain, in which: the variable to be solved is that a two-dimensional axial symmetric is non-
The magnetic potential A of linear magnetostatic field, the domain that solves is the region where two-dimensional axial symmetric nonlinear magnetostatic field;
Two, a rz rectangular coordinate system is established;
Three, governing equation in two-dimensional axial symmetric nonlinear magnetostatic field and boundary condition formula and one differential side of composition are listed
Journey group, in which:
The governing equation are as follows:
In formula, J is current density variable, and μ '=r μ, μ is the magnetic conductivity of triangular unit, and A '=rA, A are magnetic potential, and r is cross
Coordinate, z are ordinate;
The boundary condition formula are as follows:
Γ1: A '=0;
In formula, m is to solve borderline normal vector;
Four, subnetting is carried out to domain is solved using subnetting program, it is discrete for triangular element or quadrangle list that domain will be solved
Member, in each finite element unit, variableIn formula, n is the interstitial content of finite element unit, NjFor list
Shape function in member,For the size of the variate-value of corresponding node;
Five, the coefficient matrix of finite element is calculated according to the following formula:
Wherein, the calculation expression of each single item is as follows:
Six, by the matrix K in step 5eRegard the admittance matrix of circuit, b aseRegard electric current source vector as, establishes equivalent
Black box circuit model;
Seven, one section of transmission line line segment is inserted into the connecting wire any position between black box circuit model and linear circuit;
Eight, Transmission line method iteration is carried out:
Voltage signal is incident in Linear Network, the circuit, i.e. equation group: (Y are solvedlinear+YTL) A '=b+
2ViYTL, in entire iterative process, matrix (Ylinear+YTL) remain unchanged, a matrix is executed in the step 1 of iteration
LU operation splitting, subsequent step without executing again, after the completion of solution, calculates reflected voltage Vr=A '-Vi;
Voltage signal is incident in nonlinear finite element unit, the circuit, i.e. Nonlinear System of Equations are solved:It is solved using Newton iteration method, the iterative formula of kth step is as follows:Wherein:
Solving equations in each non-linear unit are individually put into independent calculating core and solve, and realize parallel
It calculates, after the completion of solution, calculates reflected voltage again
In formula, YlinearAdmittance matrix for the linear circuit being incident in Linear Network;YTLTo be incident on Linear Network
The admittance matrix of interior transmission line circuit;A ' is node voltage to be asked;VrFor the reflected voltage for being reflected into nonlinear circuit;ViFor
It is incident on the incident voltage of linear circuit;KeFor the coefficient matrix of non-linear unit;To be reflected back the biography in non-linear unit
The admittance matrix of defeated line;To be reflected back the reflected voltage in non-linear unit;JeFor Jacobian matrix;Repeatedly for kth time
For when Jacobian matrix;F is the intermediate variable in calculating;Ae′For the magnetic potential value of node each in unit.
Nine, the iteration for repeating step 8, until iteration result converges to fixed error, end is solved.
Compared with the prior art, the present invention has the advantage that
1, discrete modeling is carried out to magnetostatic field using finite element method, by the unit coefficient of each unit in finite element
Matrix is equivalent to the black box circuit model of a sealing, and FEM equations are equivalent to a nonlinear electricity to realize
Road network problem.
2, it is inserted into transmission line line segment between the black-box model of non-linear unit and the network of linear unit, utilizes transmission
The iteration rule of line, realizes transmission line interation, to realize the solution of Magnetostatic Field Problems.
3, in transmission line iterative process, each non-linear black-box model is independent from each other, its solution calculating can
To be calculated using individual computing resource (CPU GPU core), parallel accelerate is realized.
4, it is able to carry out the iterative solution of finite element using Transmission line method and black-box model, is capable of handling Various Complex
Finite element model, such as triangular element and quadrilateral units.
5, during iterative solution, global matrix Y be able to maintain it is constant, in Matrix Solving process, using LU point
Solution, it is only necessary to carry out LU decomposition in the first step of calculating, due to LU decompose it is general occupy Matrix Solving 95% or so when
Between, it will substantially reduce the iteration time of each step in this way.
6, non-linear unit can be isolated using Transmission line method with linear solution region, realizes parallel acceleration effect
Fruit.
Detailed description of the invention
Fig. 1 is black box circuit model;
Fig. 2 is the black box circuit network being inserted into after transmission line;
Fig. 3 is that voltage is incident on the equivalent circuit diagram after Linear Network;
Fig. 4 is that voltage reflection returns the equivalent circuit diagram after non-linear unit;
Fig. 5 is iterative solution flow chart;
Fig. 6 is the structure chart of certain model contactor;
Fig. 7 is the sectional view of certain model contactor;
Fig. 8 is contactor simplified model;
Fig. 9 is contactor computation model;
Figure 10 is subnetting as a result, (a) triangular element, (b) quadrilateral units;
Figure 11 is the Distribution of Magnetic Field situation for solving to obtain using Transmission line method under different conditions;
Figure 12 is that single step calculates time acceleration effect.
Specific embodiment
Further description of the technical solution of the present invention with reference to the accompanying drawing, and however, it is not limited to this, all to this
Inventive technique scheme is modified or replaced equivalently, and without departing from the range of technical solution of the present invention, should all be covered in this hair
In bright protection scope.
Specific embodiment 1: present embodiments provide for a kind of two-dimension non linearity based on black box transmission line model is quiet
The finite element solving method of magnetic field model, described method includes following steps:
One, it determines variable to be solved and solves domain, variable to be solved is a two-dimensional axial symmetric nonlinear magnetostatic field
Magnetic potential A, two-dimensional axial symmetric nonlinear magnetostatic field generates by the electric current in hot-wire coil, and each element around hot-wire coil is equal
For ferromagnetic material, solving domain is the region where two-dimensional axial symmetric nonlinear magnetostatic field.
Two, a rz rectangular coordinate system is established, r is horizontal axis, and z is the longitudinal axis.
Three, governing equation in two-dimensional axial symmetric nonlinear magnetostatic field and boundary condition formula and one differential side of composition are listed
Journey group, governing equation are as follows:Wherein, J is current density variable, and μ '=r μ, μ is
The magnetic conductivity of triangular unit, A '=rA, A are magnetic potential, and r is abscissa, and z is ordinate, boundary condition are as follows: Γ1: A '=0;M is to solve borderline normal vector.
Four, subnetting is carried out to domain using subnetting program, it is discrete for triangular element or quadrangle that domain will be solved
Unit.In each finite element unit, variableWherein, n is the interstitial content of finite element unit, NjFor
Shape function in unit,For the size of the variate-value of corresponding node.
Five, the coefficient matrix of finite element is calculated.According to Galerkin method, integral, the side of obtaining are weighted to governing equation
Journey:
Write as matrix form KeAe′=be, i.e., are as follows:
Wherein, the calculation expression of each single item is as follows:
Six, by the matrix K in step 5eRegard the admittance matrix of circuit, b aseRegard electric current source vector as, establishes equivalent
Black box circuit model, as shown in Figure 1.
It seven, is to be connected between black-box model and linear circuit by conducting wire, between black-box model and linear circuit
These conducting wires any position be inserted into one section of transmission line line segment, as shown in Fig. 2.
Eight, Transmission line method iteration is carried out.
Voltage signal is incident in Linear Network, as shown in figure 3, solving the circuit, i.e. equation group: (Ylinear+YTL)A′
=b+2ViYTL, in entire iterative process, matrix (Ylinear+YTL) remain unchanged, first moment is executed in the step 1 of iteration
The LU operation splitting of battle array, subsequent step without executing again, after the completion of solution, calculates reflected voltage Vr=A '-Vi。
Voltage signal is incident in nonlinear finite element unit, as shown in figure 4, solving the circuit, i.e. nonlinear equation
Group:The matrix is solved using Newton iteration method, and the iterative formula of kth step is as follows:Wherein:
Since each non-linear unit is isolation between each other, the solving equations in each unit are individually put
It is solved into independent calculating core, realizes parallel computation.After the completion of solution, reflected voltage is calculated again
In formula, YlinearFor the admittance matrix of linear circuit in Fig. 3;YTLIt is led for transmission line circuit in Fig. 3 dotted line frame
Receive matrix;A ' is node voltage to be asked;VrFor the reflected voltage for being reflected into nonlinear circuit;ViTo be incident on entering for linear circuit
Radio pressure;KeFor the coefficient matrix of non-linear unit;For the admittance matrix of transmission line in Fig. 4;To be reflected back in unit
Reflected voltage;JeFor Jacobian matrix;Jacobian matrix when iteration secondary for kth.
It is constantly iterated, until iteration result converges to fixed error, end is solved.Solve flow chart such as Fig. 5.
Specific embodiment 2: present embodiments provide for a kind of two-dimension non linearity based on black box transmission line model is quiet
The finite element solving method of magnetic field model, the method specific implementation step are as follows:
One, it determines variable to be solved and solves domain, variable to be solved is a two-dimensional axial symmetric nonlinear magnetostatic field
Magnetic potential A, two-dimensional axial symmetric nonlinear magnetostatic field generates by the electric current in hot-wire coil, and each element around hot-wire coil is equal
For ferromagnetic material, solving domain is the region where two-dimensional axial symmetric nonlinear magnetostatic field.Present embodiment is directed to certain model axis pair
Structure contactor is claimed to be calculated, the outside drawing of the model is as shown in fig. 6, sectional view is as shown in Figure 7.It is simplified,
And model, simplified model is as shown in figure 8, computation model is as shown in Figure 9.
Two, a rz rectangular coordinate system is established, r is horizontal axis, and z is the longitudinal axis.
Three, governing equation in two-dimensional axial symmetric nonlinear magnetostatic field and boundary condition formula and one differential side of composition are listed
Journey group, governing equation are as follows:Wherein, J is current density variable, and μ '=r μ, μ is
The magnetic conductivity of triangular unit, A '=rA, A are magnetic potential, and r is abscissa, and z is ordinate, boundary condition are as follows: Γ1: A '=0;M is to solve borderline normal vector.
Four, subnetting is carried out to domain using subnetting program, it is discrete for triangular element or quadrangle that domain will be solved
Unit, the results are shown in Figure 10 for subnetting.In each finite element unit, variableWherein, n is finite element
The interstitial content of unit, NjFor the shape function in unit,For the size of the variate-value of corresponding node.
For triangular element,Δ is the area of triangular element,
For quadrilateral units,
Five, the coefficient matrix of finite element is calculated.According to Galerkin method, integral, the side of obtaining are weighted to governing equation
Journey:
Write as matrix form KeAe′=be, i.e., are as follows:
Wherein, the calculation expression of each single item is as follows:
Six, by the matrix K in step 5eRegard the admittance matrix of circuit, b aseRegard electric current source vector as, establishes equivalent
Black box circuit model, as shown in Figure 1.
Seven, one section of transmission line line segment is inserted between black-box model and linear circuit, as shown in Figure 2.
Eight, Transmission line method iteration is carried out.
Voltage signal is incident in Linear Network, as shown in figure 3, solving the circuit, i.e. equation group: (Ylinear+YTL)A′
=b+2ViYTL, then, calculate reflected voltage Vr=A '-Vi。
Voltage signal is incident in nonlinear finite element unit, as shown in figure 4, solving the circuit, i.e. nonlinear equation
Group:The matrix is solved using Newton iteration method, and the iterative formula of kth step is as follows:Wherein:
After the completion of solution, reflected voltage is calculated again
It is constantly iterated, until iteration result converges to fixed error, end is solved.Solve flow chart such as Fig. 5.
Solve to obtain the Distribution of Magnetic Field situation under different conditions using Transmission line method, as a result as shown in Figure 11.
The different solving model of table 1
Case 1-5 is triangular element subnetting;Case 6-10 is quadrilateral units subnetting.
The error of 2 N-R iterative method of table and BB-TLM iterative method of the present invention compares
No.1-5 is to appoint that the error of five points is taken to compare in 5 model of table 1case;No.6-10 is 1 case of table, 10 mould
Type, which is appointed, takes the error of five points to compare.
Case10 model in table 1 is accelerated using 1,4,8,16,20 core cpu respectively, compares Newton iteration
The single step of method and Transmission line method calculates the time, and as seen from Figure 12, it is faster than Newton iteration method several that the single step of Transmission line method calculates the time
Times.
Case10 model in table 1 is accelerated using 1,4,8,16,20 core cpu respectively, compares Newton iteration
The total evaluation time of method and Transmission line method, when core cpu increases, the acceleration effect of Transmission line method is more than Newton iteration method
Obviously.
Newton iteration method is respectively adopted for the model in table 1 for table 3 and Transmission line method solves (CPU core is 20)
As shown in Table 3, when subnetting model increases, the calculating time of Transmission line method is faster than Newton iteration method.
Claims (1)
1. a kind of finite element solving method of the two-dimension non linearity magnetostatic model based on black box transmission line model, it is characterised in that
Described method includes following steps:
One, it determines variable to be solved and solves domain, in which: the variable to be solved is a two-dimensional axial symmetric nonlinear static
The magnetic potential A in magnetic field, the domain that solves is the region where two-dimensional axial symmetric nonlinear magnetostatic field;
Two, a rz rectangular coordinate system is established;
Three, governing equation in two-dimensional axial symmetric nonlinear magnetostatic field and boundary condition formula and one differential equation group of composition are listed,
Wherein:
The governing equation are as follows:
In formula, J is current density variable, and μ '=r μ, μ is the magnetic conductivity of triangular unit, and A '=rA, A are magnetic potential, and r is abscissa,
Z is ordinate;
The boundary condition formula are as follows:
Γ1: A '=0;
In formula, m is to solve borderline normal vector;
Four, subnetting is carried out to domain is solved using subnetting program, it is discrete for triangular element or quadrilateral units to solve domain,
In each finite element unit, variableIn formula, n is the interstitial content of finite element unit, NjIn unit
Shape function,For the size of the variate-value of corresponding node;
Five, the coefficient matrix of finite element is calculated according to the following formula:
Wherein, the calculation expression of each single item is as follows:
Six, by the matrix K in step 5eRegard the admittance matrix of circuit, b aseRegard electric current source vector as, establishes equivalent black box
Circuit model;
Seven, one section of transmission line line segment is inserted into the connecting wire any position between black box circuit model and linear circuit;
Eight, Transmission line method iteration is carried out:
Voltage signal is incident in Linear Network, the circuit, i.e. equation group: (Y are solvedlinear+YTL) A '=b+2ViYTL, whole
In a iterative process, matrix (Ylinear+YTL) remain unchanged, the LU operation splitting of a matrix is executed in the step 1 of iteration,
Subsequent step without executing again, after the completion of solution, calculates reflected voltage Vr=A '-Vi;
Voltage signal is incident in nonlinear finite element unit, the circuit, i.e. Nonlinear System of Equations are solved:It is solved using Newton iteration method, the iterative formula of kth step is as follows:Solving equations in each non-linear unit are individually put into independent meter
It calculates and is solved in core, realize parallel computation, after the completion of solution, calculate reflected voltage again
In formula, YlinearAdmittance matrix for the linear circuit being incident in Linear Network;YTLTo be incident on the biography in Linear Network
The admittance matrix of defeated line circuit;A ' is node voltage to be asked;VrFor the reflected voltage for being reflected into nonlinear circuit;ViTo be incident on
The incident voltage of linear circuit;KeFor the coefficient matrix of non-linear unit;For the transmission line that is reflected back in non-linear unit
Admittance matrix;To be reflected back the reflected voltage in non-linear unit;Jacobian matrix when iteration secondary for kth;
Nine, the iteration for repeating step 8, until iteration result converges to fixed error, end is solved.
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 true CN109408927A (en) | 2019-03-01 |
CN109408927B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112016232A (en) * | 2020-08-31 | 2020-12-01 | 中国原子能科学研究院 | Tear finite element process processing method and system |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050228601A1 (en) * | 2003-06-05 | 2005-10-13 | Fujitsu Limited | Device and method for calculating electro-magnetic field intensity by cooperation of circuit analyses and electro-magnetic wave analyses |
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 |
US20180246154A1 (en) * | 2015-09-30 | 2018-08-30 | Osaka University | Method and program for calculating potential, current, and peripheral electromagnetic field in electric circuit |
-
2018
- 2018-10-13 CN CN201811192762.0A patent/CN109408927B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050228601A1 (en) * | 2003-06-05 | 2005-10-13 | Fujitsu Limited | Device and method for calculating electro-magnetic field intensity by cooperation of circuit analyses and electro-magnetic wave analyses |
US20180246154A1 (en) * | 2015-09-30 | 2018-08-30 | Osaka University | Method and program for calculating potential, current, and peripheral electromagnetic field in electric circuit |
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 |
---|
LI BING等: ""Influence of electromagnetic pulse on electromagnetic relay by using equivalent circuit model"", 《26TH INTERNATIONAL CONFERENCE ON ELECTRICAL CONTACTS (ICEC 2012)》 * |
周建美: ""各向异性地层中可控源电磁法一维全参数反演及三维有限体积正演算法研究"", 《中国博士学位论文全文数据库 基础科学辑》 * |
张旭锋: ""传输线理论及电磁兼容计算的半解析方法研究"", 《中国博士学位论文全文数据库 信息科技辑》 * |
张雷: ""基于巨磁阻抗效应的弱磁传感器技术研究"", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112016232A (en) * | 2020-08-31 | 2020-12-01 | 中国原子能科学研究院 | Tear finite element process processing method and system |
Also Published As
Publication number | Publication date |
---|---|
CN109408927B (en) | 2022-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jafari et al. | Using dragonfly algorithm for optimization of orthotropic infinite plates with a quasi-triangular cut-out | |
Song et al. | Stability, steady‐state bifurcations, and Turing patterns in a predator–prey model with herd behavior and prey‐taxis | |
CN103902785B (en) | One kind is based on polynary probabilistic structural finite element model updating method | |
Shivanian | Meshless local Petrov–Galerkin (MLPG) method for three-dimensional nonlinear wave equations via moving least squares approximation | |
Wei et al. | Dynamical analysis of the generalized Sprott C system with only two stable equilibria | |
CN106777827A (en) | A kind of electromechanical electromagnetism hybrid simulation method and system | |
CN110275733A (en) | The GPU parallel acceleration method of phonon Boltzmann equation is solved based on finite volume method | |
CN105005675A (en) | Composite insulator electric field optimization method based on multi-target genetic algorithm | |
Liu et al. | Bifurcation of a predator–prey model with disease in the prey | |
CN108921288A (en) | Neural network activates processing unit and the neural network processor based on the device | |
CN110502785A (en) | A kind of Three-dimensional Time Domain calculates the electromagnetism numerical method of waveguide S parameter | |
CN107038302A (en) | A kind of semi-analytical calculation method in the Eddy Current Nondestructive Testing magnetic field containing columnar defects | |
CN109408927A (en) | A kind of two-dimentional magnetostatic field parallel finite element method accelerated method based on black box transmission line model | |
CN106777729A (en) | A kind of algorithms library simulation and verification platform implementation method based on FPGA | |
CN110362877A (en) | The environment electromagnetics analysis of scattering method of uncertain factor | |
CN106649939B (en) | The method for solving of 2D Axisymmetric Nonlinear magnetostatic model based on transmission line interation | |
Nicolas et al. | Improved adaptive mesh refinement for conformal hexahedral meshes | |
De Su et al. | Numerical solution for a kind of nonlinear telegraph equations using radial basis functions | |
Yin et al. | Multi-step Prediction Algorithm of Traffic Flow Chaotic Time Series Based on Volterra Neural Network. | |
CN103942376A (en) | Collision correction method for real-time hair motion simulation | |
CN106651979A (en) | Projection-based space curved surface quadrangle grid dissection method | |
Moirangthem et al. | GIS integrated automation of a near real-time power-flow service for electrical grids | |
CN106096176A (en) | The modeling method of the antenna model of the pigtail of shielded cable and device | |
Monsegny et al. | Fast marching method in seismic ray tracing on parallel GPU devices | |
CN107609274B (en) | Two-dimensional static magnetic field parallel finite element method based on transmission line and level scheduling method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |