Electromagnetic field calculation method and device for metal part, equipment, medium and product

Abstract

The invention relates to the technical field of electromagnetic field calculation, and discloses an electromagnetic field calculation method and device for a metal part, equipment, a medium and a product, and the method comprises the steps: obtaining a discretized electric field mathematical model and a discretized magnetic field mathematical model which correspond to the metal part, and each of the electric field mathematical model and the magnetic field mathematical model comprises an update coefficient; an included angle formed by the metal surfaces on the two sides of a target edge is detected, and the target edge is the edge, flush with the divided grid lines, in the at least one edge; determining a correction coefficient according to the included angle; according to the correction coefficient, updating coefficients corresponding to all grid surfaces covered by the target edge are corrected, and a corrected electric field mathematical model and a corrected magnetic field mathematical model are obtained; and performing time domain iteration on the corrected electric field mathematical model and the corrected magnetic field mathematical model to obtain an electromagnetic field distribution result of the metal part. According to the method, the singular field near the edge of the metal part is corrected, so that the electromagnetic field distribution result can be determined more efficiently and accurately.

Description

Technical Field

[0001] The present invention relates to the technical field of electromagnetic field calculation, and particularly relates to a method, device, equipment, medium and product for calculating the electromagnetic field of metal components. Background Art

[0002] With the increasing integration of electronic devices, the problem of electromagnetic interference between metal components has become increasingly prominent. By calculating the electromagnetic field distribution of metal components, it is possible to predict and optimize the layout and shielding measures of metal components in advance, reduce the electromagnetic interference between different metal components, and ensure the normal operation of electronic devices.

[0003] However, traditional electromagnetic field calculation methods (such as the finite difference method, the finite element method or the boundary element method, etc.) have low accuracy in the electromagnetic field results obtained when dealing with special scenarios such as singular fields. Specifically, since singular fields are usually accompanied by very sharp local changes, standard calculation grids are difficult to accurately capture these changes, resulting in low accuracy of the calculated electromagnetic field results, and even numerical divergence may occur. Among them, a singular field refers to the phenomenon that the electromagnetic field is discontinuous near the edges where the metal components mutate.

[0004] Currently, the accuracy of the electromagnetic field results calculated under singular fields is mainly improved by refining the grid. However, extremely fine grid division consumes a large amount of computing resources, greatly reducing the computing efficiency and making it difficult to meet the requirements of large-scale or real-time calculations. Summary of the Invention

[0005] In view of this, the present invention provides a method, device, equipment, medium and product for calculating the electromagnetic field of metal components to improve the problem of low accuracy of the electromagnetic field results obtained by traditional electromagnetic field calculation methods when dealing with special scenarios such as singular fields.

[0006] In a first aspect, the present invention provides a method for calculating the electromagnetic field of metal components. The metal component has at least one edge. The method includes: obtaining a discretized electric field mathematical model and a magnetic field mathematical model corresponding to the metal component, where both the electric field mathematical model and the magnetic field mathematical model include update coefficients; detecting the angle formed by the metal surfaces on both sides of the target edge, where the target edge is the edge that is flush with the divided grid lines among the at least one edge; determining a correction coefficient according to the angle; correcting the update coefficients corresponding to all grid faces covered by the target edge according to the correction coefficient to obtain a corrected electric field mathematical model and a corrected magnetic field mathematical model; performing time-domain iteration on the corrected electric field mathematical model and the corrected magnetic field mathematical model to obtain the electromagnetic field distribution result of the metal component.

[0007] The electromagnetic field calculation method for metal components provided in this embodiment, after obtaining the discretized electric field mathematical model and magnetic field mathematical model corresponding to the metal component, detects the included angle formed by the metal surfaces on both sides of the target edge, determines the correction coefficient according to the included angle, and then corrects the update coefficients corresponding to all grid faces covered by the target edge according to the correction coefficient to obtain the corrected electric field mathematical model and the corrected magnetic field mathematical model. Finally, time-domain iteration is performed on the corrected electric field mathematical model and the corrected magnetic field mathematical model to obtain the electromagnetic field distribution result of the metal component. In this embodiment, the electromagnetic field mathematical model is corrected by the correction coefficient determined by the included angle, which can correct the singular field near the metal edge and improve the stability and accuracy of the electromagnetic field distribution result. Moreover, compared with the traditional correction method of increasing the number of grids, the present application uses included angle correction, which can reduce the calculation amount, achieve a large correction effect with a minimal calculation cost, and significantly improve the calculation efficiency and calculation accuracy.

[0008] In an alternative embodiment, determining the correction coefficient according to the included angle includes: determining the correction coefficient according to the included angle and the following formula:

[0009]

[0010] where, v1 represents the correction coefficient, and α represents the included angle.

[0011] In an alternative embodiment, correcting the update coefficients corresponding to all grid faces covered by the target edge according to the correction coefficient includes: determining the correction amount according to the correction coefficient and the following formula:

[0012]

[0013] where, C s represents the correction amount, and v1 represents the correction coefficient;

[0014] Correct the update coefficients corresponding to all grid faces covered by the target edge according to the correction amount.

[0015] In this embodiment, after determining the correction coefficient, determining the correction amount based on the correction coefficient and the formula, and then correcting the update coefficients corresponding to all grid faces covered by the target edge according to the correction amount can improve the accuracy of the correction, thereby further improving the accuracy of the electromagnetic field distribution result.

[0016] In an alternative embodiment, correcting the update coefficients corresponding to all grid faces covered by the target edge according to the correction amount includes: multiplying the correction amount and the update coefficients corresponding to all grid faces covered by the target edge in sequence to correct the update coefficients corresponding to all grid faces covered by the target edge.

[0017] In an alternative embodiment, before detecting the angle formed by the metal surfaces on both sides of the target edge, the method further includes: traversing all the edges of the metal component; sequentially determining whether each edge is flush with the grid lines, and determining the edge flush with the grid lines as the target edge.

[0018] In an alternative embodiment, the metal component is a semiconductor chip, an electronic device in a vehicle, or an electronic device in a terminal device.

[0019] In a second aspect, the present invention provides an electromagnetic field calculation device for a metal component, including: an acquisition module for acquiring a discretized electric field mathematical model and a magnetic field mathematical model corresponding to the metal component, wherein both the electric field mathematical model and the magnetic field mathematical model include update coefficients; a detection module for detecting the angle formed by the metal surfaces on both sides of the target edge, wherein the target edge is the edge flush with the divided grid lines among at least one edge; a determination module for determining a correction coefficient according to the angle; a correction module for correcting the update coefficients corresponding to all the grid faces covered by the target edge according to the correction coefficient to obtain a corrected electric field mathematical model and a corrected magnetic field mathematical model; and an iteration module for performing time-domain iteration on the corrected electric field mathematical model and the corrected magnetic field mathematical model to obtain the electromagnetic field distribution result of the metal component.

[0020] In an alternative embodiment, the determination module includes: a determination unit for determining a correction coefficient according to the angle and the following formula:

[0021]

[0022] wherein, v1 represents the correction coefficient, and α represents the angle.

[0023] In an alternative embodiment, the correction module includes: a first processing unit for determining a correction amount according to the correction coefficient and the following formula:

[0024]

[0025] wherein, C s represents the correction amount, and v1 represents the correction coefficient;

[0026] a correction unit for correcting the update coefficients corresponding to all the grid faces covered by the target edge according to the correction amount.

[0027] In an alternative embodiment, the correction unit includes: a correction subunit for sequentially multiplying the correction amount by the update coefficients corresponding to all the grid faces covered by the target edge to correct the update coefficients corresponding to all the grid faces covered by the target edge.

[0028] In an alternative embodiment, the apparatus further includes: a traversal module configured to traverse all edges of the metal component; and a processing module configured to sequentially determine whether an edge is flush with the grid line and determine the edge flush with the grid line as the target edge.

[0029] In an alternative embodiment, the metal component is a semiconductor chip, an electronic device in a vehicle, or an electronic device in a terminal device.

[0030] In a third aspect, the present invention provides a computer device, including: a memory and a processor, which are communicatively connected to each other. The memory stores computer instructions, and the processor executes the computer instructions to perform the method according to the first aspect or any corresponding embodiment thereof.

[0031] In a fourth aspect, the present invention provides a computer-readable storage medium storing computer instructions for causing a computer to execute the method according to the first aspect or any corresponding embodiment thereof.

[0032] In a fifth aspect, the present invention provides a computer program product including computer instructions for causing a computer to execute the method according to the first aspect or any corresponding embodiment thereof. Description of the Drawings

[0033] In order to more clearly illustrate the specific embodiments of the present invention or the technical solutions in the related art, the following will briefly introduce the drawings required for use in the description of the specific embodiments or the related art. Obviously, the drawings in the following description are some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0034] Figure 1 is a schematic flowchart of a method for calculating the electromagnetic field of a metal component according to an embodiment of the present invention;

[0035] Figure 2 is a schematic cross-sectional view of a metal component according to an embodiment of the present invention;

[0036] Figure 3 is a schematic diagram of the electromagnetic field position during electromagnetic field polarization according to an embodiment of the present invention;

[0037] Figure 4 is a schematic flowchart of another method for calculating the electromagnetic field of a metal component according to an embodiment of the present invention;

[0038] Figure 5 is a schematic diagram of a metal component adsorbed at the grid according to an embodiment of the present invention;

[0039] Figure 6 It is a schematic flowchart of another electromagnetic field calculation method for metal components according to an embodiment of the present invention;

[0040] Figure 7 It is another simulation schematic diagram of a metal component according to an embodiment of the present invention;

[0041] Figure 8 It is a structural block diagram of an electromagnetic field calculation device for a metal component according to an embodiment of the present invention;

[0042] Figure 9 It is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. Specific embodiments

[0043] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0044] The electromagnetic field calculation method for metal components provided by the present invention is used to calculate the electromagnetic field distribution result of a metal component with edges. Among them, the metal component can be a semiconductor chip, an electronic device in a vehicle, or an electronic device in a terminal device. The edge can refer to the edge part of the metal component, or it can be the edge in the metal component that may cause singular behaviors of the electromagnetic field (such as discontinuity or sharp change in field strength, etc.). For example, the right-angle edge of a semiconductor chip, the edge of a display screen, etc.

[0045] Exemplarily, the semiconductor chip can be a Central Processing Unit (CPU), a Graphics Processing Unit (GPU), a Microcontroller Unit (MCU), or a storage chip, etc. The electronic device in the vehicle can be an Engine Control Unit (ECU), a Body Control Module (BCM), or a distance sensor, etc. The terminal device can be a device such as a mobile phone, a tablet computer, or a smart watch, etc. The electronic device in the terminal device can be a display screen, a camera, or a wireless communication module, etc.

[0046] The electromagnetic field calculation method for metal components provided by the present invention can more efficiently and accurately determine the electromagnetic field distribution result by correcting the singular field near the edge of the metal component when calculating the electromagnetic field distribution of the metal component.

[0047] An embodiment of the present invention provides an embodiment of a method for calculating the electromagnetic field of a metal component. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although the logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in a different order than here.

[0048] In this embodiment, a method for calculating the electromagnetic field of a metal component is provided, which can be used in an electromagnetic field calculation device for metal components. The electromagnetic field calculation device for metal components can be a mobile phone, a tablet computer, a computer, etc. Figure 1 It is a schematic flowchart of a method for calculating the electromagnetic field of a metal component according to an embodiment of the present invention, as Figure 1 shown, the method includes the following steps:

[0049] Step S101, obtain the discretized electric field mathematical model and magnetic field mathematical model corresponding to the metal component.

[0050] Among them, both the electric field mathematical model and the magnetic field mathematical model include update coefficients.

[0051] Specifically, the basic characteristics of the electromagnetic field are the changes in the electric field and magnetic field caused by the spatial distribution and motion state of the charge system. The essential laws of the electromagnetic field are described by Maxwell's equations. In the calculation and analysis of the electromagnetic field, an electric field mathematical model and a magnetic field mathematical model can be established through Maxwell's equations, and then calculated and analyzed according to the boundary conditions and physical characteristics to obtain the distribution and change laws of the electromagnetic field in space.

[0052] Maxwell's equations are the basic equations of classical electrodynamics, which unify the theories of electricity and magnetism in mathematical form and reveal the mutual connection and dependence between the electric field and the magnetic field.

[0053] Exemplarily, the discretized electric field mathematical model and magnetic field mathematical model corresponding to the metal component can be determined by means such as the Finite Integration Technique (FIT) algorithm, the finite element method, or the finite difference method.

[0054] The FIT algorithm is based on the integral form of Maxwell's equations. It discretizes the calculation domain into an orthogonal grid, defines the voltage and magnetic flux as unknowns on the grid, and discretizes the continuous Maxwell's equations into a set of algebraic equations. When operating specifically, first divide the calculation area into a series of small cubic blocks, sample the electric field strength at the center of the edges of each cubic block, sample the magnetic field strength at the center of the faces of the cubic block, approximate the integral equation as a discrete summation form through spatial and temporal sampling, and rewrite it in matrix form for solution.

[0055] Exemplarily, the electromagnetic calculation device for a metal component can, based on the user's input operation, obtain the discretized electric field mathematical model and magnetic field mathematical model corresponding to the established metal component, or can use electromagnetic simulation software to determine the discretized electric field mathematical model and magnetic field mathematical model corresponding to the metal component.

[0056] Next, taking the metal component as shown Figure 2 as an example, the discretized electric field mathematical model and magnetic field mathematical model corresponding to the metal component determined by the FIT algorithm will be described.

[0057] Specifically, after determining the metal component to be calculated, computer-aided design software can be used to perform geometric modeling on the electromagnetic system of the metal component, clarify the shape, size of the metal component and its position in the entire calculation region, and determine the electromagnetic characteristic parameters (such as magnetic permeability and permittivity, etc.) of the metal component and the surrounding medium.

[0058] Then, the established geometric model of the metal component is divided into discrete grid cells, and based on the integral form of Maxwell's equations, the integral relationships of the electric field E and the magnetic field H are expressed on the discrete grid cells. By integrating the electric field and the magnetic field on the grid edges and grid surfaces, the continuous Maxwell's equations are transformed into algebraic equations, and this algebraic equation can be used as the discretized electric field mathematical model and magnetic field mathematical model corresponding to the metal component.

[0059] Based on the definitions of the electric field and the magnetic field in the FIT algorithm, the positions of the electromagnetic fields during electric field polarization (magnetic field polarization) can be as shown Figure 3 , where the horizontal arrow represents H y (E y ), the vertical arrow represents H x (E x ), the cross at the grid node represents H z (E z ), and i and j are the identifiers of the divided grids. H x represents the component of the magnetic field H in the x-axis direction, H y represents the component of the magnetic field H in the y-axis direction, H z represents the component of the magnetic field H in the z-axis direction, E x represents the component of the electric field E in the x-axis direction, E y represents the component of the electric field E in the y-axis direction, E z represents the component of the electric field E in the z-axis direction.

[0060] As shown Figure 2 , the cross-section of the metal component is triangular, and the metal component remains unchanged along the z-axis. Therefore, the influence of the z-axis on the field (electric field and magnetic field) values near the metal component can be ignored. The region occupied by the metal component is denoted as SC , the boundary of the metal component is denoted as C, and the region outside the boundary is denoted as S. When the medium in S is a homogeneous medium, the field behavior in the region near the metal component is independent of the properties of the medium and only depends on the angle α of the metal component.

[0061] Therefore, it can be assumed that the region S is completely vacuum. For fields that do not depend on the z-axis, the electromagnetic field equations will be simplified into two sets of uncoupled equations, namely the polarization of the electric field E and the polarization of the magnetic field H.

[0062] The expression of the electric field E determined based on Maxwell's equations can be shown as in formulas (1.a), (1.b), and (1.c):

[0063]

[0064] In formula (1.a), μ0 represents the magnetic permeability in vacuum, which describes the relationship between the magnetic field and the magnetic flux density in vacuum, and its value can be 4π×10 -7 N / A 2 ; represents the partial derivative with respect to time t, which is used to describe the rate of change of a physical quantity with time; represents the partial derivative with respect to the y coordinate. In formula (1.b), represents the partial derivative with respect to the x coordinate. In formula (1.c), ε0 represents the vacuum permittivity, which is used to describe the relationship between the electric field and the electric displacement vector in vacuum, and its value can be approximately 8.854×10 -12 C 2 / (N·m 2 ).

[0065] The expression of the magnetic field H determined based on Maxwell's equations can be shown as in formulas (2.a), (2.b), and (2.c):

[0066]

[0067] There are boundary conditions on the surface C of the metal component, and the boundary conditions are shown as in formula (3):

[0068] n×E = 0 (3)

[0069] In formula (3), n represents the unit normal vector of the surface C of the metal component, which is used to determine the direction of the surface, and n×E = 0 means that on the surface of an ideal metal conductor, the tangential component of the electric field is zero.

[0070] Then, introduce the finite-difference time-domain approximate representation of the electromagnetic field values, and this expression is shown as in formula (4):

[0071] F n (i,j) = F(ih,jh,nk) (4)

[0072] In formula (4), F n (i, j) represents the value of a certain component of the electromagnetic field at the grid position (i, j) at the nth time step. F can represent H x , H y , H z , E x , E y or E z ; h = Δx = Δy represents the grid step size on the x-axis and the grid step size in the y-axis direction, which is used to discretize the continuous space; k = Δt represents the time step size, which is used to discretize the continuous time.

[0073] By performing a difference approximation on formulas (1.a), (1.b), and (1.c) based on formula (4), and performing a difference approximation on formulas (2.a), (2.b), and (2.c), a discretized mathematical model of the electric field and a discretized mathematical model of the magnetic field can be obtained.

[0074] The discretized mathematical model of the magnetic field can be as shown in formulas (5.a), (5.b), and (6.c), and the discretized mathematical model of the electric field can be as shown in formulas (6.a), (6.b), and (5.c):

[0075]

[0076] In formula (5.a), n represents the current time step size, represents the interval of half a time step size, n + 1 represents the next time step size, which is an embodiment of the finite-difference time-domain algorithm for alternately updating the electric field and the magnetic field; i corresponds to the grid position on the x-axis, j corresponds to the grid position on the y-axis, and represent the positions between two integer-numbered grids.

[0077]

[0078] Specifically, in the above-mentioned discretized mathematical models of the electric field and the magnetic field, and represent the update coefficients.

[0079] Step S102, detect the included angle formed by the metal surfaces on both sides of the target edge.

[0080] Among them, the target edge is the edge that is flush with the divided grid lines among at least one edge.

[0081] Exemplarily, after determining the target edge among at least one edge of the metal component, the included angle formed by the metal surfaces on both sides of the target edge can be determined based on the constructed geometric model of the metal component.

[0082] Specifically, after constructing the geometric model of the metal component, the parameters of each face and each edge of the metal component are clear. The included angle can be determined by reading the normal vectors of the two faces connected to the target edge and using the vector dot product formula. Alternatively, the coordinates of the two endpoints of the target edge and the relevant coordinate information of the metal surfaces on both sides of the target edge can be obtained from the constructed geometric model of the metal component. Vectors can be calculated through the coordinates, and then the included angle can be solved using vector operations.

[0083] Step S103: Determine the correction coefficient according to the included angle.

[0084] Specifically, the singular field behavior caused by the target edge will cause sharp changes in the electromagnetic field in the local area (such as electromagnetic field discontinuity). If the above-mentioned discretized electromagnetic field mathematical models (electric field mathematical model and magnetic field mathematical model) are forcibly used for solution calculation, it will lead to large errors in the finally determined electromagnetic field distribution result.

[0085] This application has found through research that the field behavior in the area near the metal component is related to the included angle corresponding to the edge. By correcting the discretized electromagnetic field mathematical model with the correction coefficient determined according to the included angle, the calculation result can be optimized.

[0086] Exemplarily, the first correspondence relationship between the included angle and the correction coefficient can be determined by the designer and pre-configured in the memory of the electromagnetic field calculation device of the metal component. After determining the included angle, the included angle can be input into the first correspondence relationship, and then the output of the first correspondence relationship can be determined as the correction coefficient.

[0087] Step S104: According to the correction coefficient, correct the update coefficients corresponding to all grid faces covered by the target edge to obtain the corrected electric field mathematical model and the corrected magnetic field mathematical model.

[0088] Specifically, after determining the correction coefficient, the electromagnetic field mathematical model can be directly corrected according to the product of the correction coefficient and the update coefficients corresponding to all grid faces covered by the target edge, and then the corrected electric field mathematical model and the corrected magnetic field mathematical model can be obtained.

[0089] Step S105: Perform time-domain iteration on the corrected electric field mathematical model and the corrected magnetic field mathematical model to obtain the electromagnetic field distribution result of the metal component.

[0090] Specifically, after determining the corrected electromagnetic field mathematical model, the unknown quantities of the electric field and the magnetic field can be solved by using solution methods such as Gaussian elimination method or conjugate gradient method to obtain the electromagnetic field distribution result of the metal component.

[0091] Exemplarily, after determining the electromagnetic field distribution result, the electromagnetic field distribution result can also be visually processed (such as drawing the distribution diagrams of the electric field and magnetic field). Through the visual result, the electromagnetic field distribution around the metal component can be intuitively observed, the influence of the metal component on the electromagnetic field can be evaluated, and a basis for further design and optimization can be provided.

[0092] For the electromagnetic field calculation method of the metal component provided in this embodiment, after obtaining the discretized electric field mathematical model and magnetic field mathematical model corresponding to the metal component, the included angle formed by the metal surfaces on both sides of the target edge is detected, and the correction coefficient is determined according to the included angle. Then, the update coefficients corresponding to all the grid faces covered by the target edge are corrected according to the correction coefficient to obtain the corrected electric field mathematical model and the corrected magnetic field mathematical model. Finally, time-domain iteration is performed on the corrected electric field mathematical model and the corrected magnetic field mathematical model to obtain the electromagnetic field distribution result of the metal component. In this embodiment, the correction coefficient determined by the included angle is used to correct the electromagnetic field mathematical model, which can correct the singular field near the metal edge and improve the stability and accuracy of the electromagnetic field distribution result. Moreover, compared with the traditional correction method of increasing the number of grids, the present application uses included angle correction, which can reduce the calculation amount, achieve a large correction effect with a minimal calculation cost, and significantly improve the calculation efficiency and calculation accuracy.

[0093] In this embodiment, an electromagnetic field calculation method for a metal component is provided, which can be used in an electromagnetic field calculation device for a metal component. Figure 4 It is a schematic flowchart of another electromagnetic field calculation method for a metal component according to an embodiment of the present invention, as Figure 4 shown. The method includes the following steps:

[0094] Step S401: Obtain the discretized electric field mathematical model and magnetic field mathematical model corresponding to the metal component.

[0095] For details, please refer to Figure 1 step S101 of the embodiment shown, which will not be elaborated here.

[0096] Step S402: Detect the included angle formed by the metal surfaces on both sides of the target edge.

[0097] For details, please refer to Figure 1 step S102 of the embodiment shown, which will not be elaborated here.

[0098] Step S403: Determine the correction coefficient according to the included angle and formula (7):

[0099]

[0100] In formula (7), v1 represents the correction coefficient, and α represents the included angle. <X

[0101] Specifically, after determining the included angle, the included angle can be substituted into formula (7), and the output result of formula (7) is the correction coefficient.

[0102] Step S404: According to the correction coefficient, correct the update coefficients corresponding to all the grid faces covered by the target edge.

[0103] In some examples, before performing the above step S404, the electromagnetic field calculation method of the metal component further includes: determining the corresponding relationship between the correction coefficient and the correction amount.

[0104] Specifically, as Figure 5 shown, the inflection point of the metal component is adsorbed at the grid of the Yee grid. From the above formula (3), it can be obtained that E z = 0 on the boundary of the metal component, and because of the change of E z at the inflection point, it will cause E z to be non-differentiable near the inflection point. Therefore, the magnetic field components H x and H y at a distance of h / 2 from the inflection point no longer apply to formula (5.a) and formula (5.b). If approximate differentiation is forced, huge errors will be introduced.

[0105] To overcome this error, the magnetic field H is represented and solved in polar coordinates near the inflection point. The expression of the magnetic field in polar coordinates can be as shown in formula (8):

[0106]

[0107] In formula (8), φ represents the azimuthal coordinate in the polar coordinate system, which is used to determine the angle of a point in the plane relative to the reference axis, r represents the radial coordinate in the polar coordinate system, which represents the distance from the origin to a certain point in space, H φ represents the component of the magnetic field H in the direction of the direction angle φ, represents the partial derivative with respect to the radial coordinate r.

[0108] The electric field and magnetic field near the metal component can be approximately expressed as formulas (9) and (10):

[0109]

[0110] In formula (9), c1(t), c2(t) and c n (t) represent the first coefficient functions related to time t, which take different values at different times and are used to adjust the amplitudes of the terms in the electric field and magnetic field expressions; v1, v2 and v n represent exponential constants, which determine the power of r and thus affect the variation law of the electric field and magnetic field components with the radial distance r. In formula (10), Y1(t), Y2(t) and Yn (t) represents a second coefficient function related to time t.

[0111] Then, substitute equations (9) and (10) into equation (8), and let r = h / 2 to obtain the expression of the magnetic field, which is shown in equation (11):

[0112]

[0113] After that, use the corresponding φ values to make H φ fall on the x-axis and y-axis respectively, and the corresponding magnetic field components H x and H y can be obtained. By comparing this expression with equations (5.a) and (5.b), it can be determined that there is an additional coefficient in front of the electric field component E z , and the coefficient in front of the electric field component E z can be determined as the corresponding relationship between the correction coefficient and the correction amount.

[0114] Exemplarily, the corresponding relationship between the correction coefficient and the correction amount can be shown as the following equation (12).

[0115] Among them, the Yee grid is a kind of staggered grid, which arranges the electric and magnetic field components alternately in space and time. In three-dimensional space, the electric field components (E x , E y , E z ) and the magnetic field components (H x , H y , H z ) are located at different grid points respectively. For example, the electric field component E x is located at the face center of the cube unit, while the magnetic field component H x is located at the edge center of the cube unit. This staggered arrangement can accurately describe the spatial distribution and interaction of the electromagnetic field.

[0116] Specifically, the above step S404 includes:

[0117] Step S4041, determine the correction amount according to the correction coefficient and equation (12):

[0118]

[0119] In equation (12), C s represents the correction amount.

[0120] Exemplarily, after determining the correction coefficient, the correction coefficient can be substituted into equation (12), and then the output of equation (12) can be determined as the correction amount.

[0121] Step S4042: According to the correction amount, correct the update coefficients corresponding to all the mesh faces covered by the target edge.

[0122] Exemplarily, the correction amount and the update coefficients corresponding to the mesh faces covered by the target edge can be multiplied in sequence to correct the update coefficients corresponding to the mesh faces covered by the target edge.

[0123] Step S405: Perform time-domain iteration on the corrected electric field mathematical model and the corrected magnetic field mathematical model to obtain the electromagnetic field distribution result of the metal component.

[0124] For details, please refer to Figure 1 Step S105 of the illustrated embodiment, which will not be elaborated here.

[0125] For the electromagnetic field calculation method of the metal component provided in this embodiment, after determining the correction coefficient, the correction amount is determined based on the correction coefficient and the formula, and then the update coefficients corresponding to all the mesh faces covered by the target edge are corrected according to the correction amount, which can improve the accuracy of the correction, thereby further improving the accuracy of the electromagnetic field distribution result.

[0126] In this embodiment, an electromagnetic field calculation method for a metal component is provided, which can be used in an electromagnetic field calculation device for a metal component. Figure 6 It is a schematic flowchart of another electromagnetic field calculation method for a metal component according to an embodiment of the present invention. As Figure 6 shown, the method includes the following steps:

[0127] Step S601: Obtain the discretized electric field mathematical model and magnetic field mathematical model corresponding to the metal component.

[0128] For details, please refer to Figure 1 Step S101 of the illustrated embodiment, which will not be elaborated here.

[0129] Step S602: Traverse all the edges of the metal component.

[0130] Step S603: Determine in sequence whether the edge is flush with the grid line, and determine the edge flush with the grid line as the target edge.

[0131] Specifically, after constructing the geometric model of the metal component, the information of all the edges of the metal component can be obtained, and then it is determined whether the current edge is flush with the grid line. If it is flush, the edge is determined as the target edge, and the following step S604 is executed; if it is not flush, it is determined whether the next edge is flush with the grid line until all the edges of the metal component are determined.

[0132] Exemplarily, it is possible to determine whether an edge is flush with a grid line by calculating the relationship between the coordinates of the edge endpoints and the coordinates of the grid lines. For a two-dimensional grid, if the abscissa or ordinate of the edge endpoints is equal to the corresponding coordinate value of the grid line, it indicates that the edge is flush with the grid line in that direction.

[0133] In other embodiments, it is also possible to process an image containing an edge and grid lines, and use an edge detection algorithm (such as the Canny algorithm) to extract the edge information in the image; then, analyze whether the edges of the extracted edge and grid lines coincide or have a specific parallel relationship to determine whether the edge is flush with the grid line.

[0134] Step S604, detect the included angle formed by the metal surfaces on both sides of the target edge.

[0135] For details, please refer to Figure 1 Step S102 of the illustrated embodiment, which will not be elaborated here.

[0136] Step S605, determine the correction coefficient according to the included angle.

[0137] For details, please refer to Figure 4 Step S403 of the illustrated embodiment, which will not be elaborated here.

[0138] Step S606, according to the correction coefficient, correct the update coefficients corresponding to all grid faces covered by the target edge to obtain a corrected electric field mathematical model and a corrected magnetic field mathematical model.

[0139] For details, please refer to Figure 4 Step S404 of the illustrated embodiment, which will not be elaborated here.

[0140] Step S607, perform time-domain iteration on the corrected electric field mathematical model and the corrected magnetic field mathematical model to obtain the electromagnetic field distribution result of the metal component.

[0141] For details, please refer to Figure 1 Step S105 of the illustrated embodiment, which will not be elaborated here.

[0142] In this embodiment, by traversing all the edges of the metal component to determine whether the edges are flush with the grid lines in turn, it is possible to check each edge of the metal component to avoid omission, thereby ensuring the accuracy of the corrected electromagnetic field mathematical model.

[0143] Exemplarily, the electromagnetic field calculation method for metal components provided in this application is applicable to the case where there are many metal singular edges. For example, as Figure 7As shown, a plurality of stacked semiconductor chips 710 have a large number of metal singular edges (the green edge portions of the chips). In this case, if the method provided in the present application is not adopted, there will be a large error between the finally simulated electromagnetic field distribution result and the measured result at high frequencies. However, after adopting the method provided in the present application, the simulated electromagnetic field distribution result is basically consistent with the measured high-frequency result.

[0144] In this embodiment, an electromagnetic field calculation device for metal components is further provided. This device is used to implement the above-mentioned embodiments and preferred implementation manners, and those that have been described will not be repeated. As used hereinafter, the term "module" can be a combination of software and / or hardware that can achieve a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, implementation in hardware, or a combination of software and hardware is also possible and contemplated.

[0145] This embodiment provides an electromagnetic field calculation device for metal components, as Figure 8 shown, including:

[0146] An acquisition module 801, configured to acquire a discretized electric field mathematical model and a magnetic field mathematical model corresponding to the metal component, wherein both the electric field mathematical model and the magnetic field mathematical model include update coefficients;

[0147] A detection module 802, configured to detect the included angle formed by the metal surfaces on both sides of the target edge, where the target edge is the edge that is flush with the divided grid line among at least one edge;

[0148] A determination module 803, configured to determine a correction coefficient according to the included angle;

[0149] A correction module 804, configured to correct the update coefficients corresponding to all grid faces covered by the target edge according to the correction coefficient, to obtain a corrected electric field mathematical model and a corrected magnetic field mathematical model;

[0150] An iteration module 805, configured to perform time-domain iteration on the corrected electric field mathematical model and the corrected magnetic field mathematical model, to obtain the electromagnetic field distribution result of the metal component.

[0151] In some alternative implementation manners, the determination module 803 includes:

[0152] A determination unit, configured to determine a correction coefficient according to the included angle and the following formula:

[0153]

[0154] wherein, v1 represents the correction coefficient, and α represents the included angle.

[0155] In some alternative implementation manners, the correction module 804 includes:

[0156] A first processing unit for determining a correction amount according to a correction coefficient and the following formula:

[0157]

[0158] where C s represents the correction amount, and v1 represents the correction coefficient;

[0159] A correction unit for correcting the update coefficients corresponding to all grid faces covered by the target edge according to the correction amount.

[0160] In some alternative embodiments, the correction unit includes:

[0161] A correction subunit for sequentially multiplying the correction amount by the update coefficients corresponding to all grid faces covered by the target edge to correct the update coefficients corresponding to all grid faces covered by the target edge.

[0162] In some alternative embodiments, the apparatus further includes:

[0163] A traversal module for traversing all edges of the metal component;

[0164] A processing module for sequentially determining whether an edge is flush with a grid line and determining the edge flush with the grid line as the target edge.

[0165] In some alternative embodiments, the metal component is a semiconductor chip, an electronic device in a vehicle, or an electronic device in a terminal device.

[0166] The further function descriptions of the above-mentioned various modules and units are the same as those in the corresponding foregoing embodiments, and will not be elaborated herein.

[0167] The electromagnetic field calculation apparatus for the metal component in this embodiment is presented in the form of functional units. Here, the unit refers to an application specific integrated circuit (ASIC), a processor and a memory that execute one or more software or fixed programs, and / or other devices that can provide the above functions.

[0168] This embodiment of the present invention further provides a computer device having the above-mentioned Figure 8 shown electromagnetic field calculation apparatus for the metal component.

[0169] Please refer to Figure 9 , Figure 9 which is a schematic structural diagram of a computer device provided by an alternative embodiment of the present invention. As shown in Figure 9As shown, the computer device includes: one or more processors 910, a memory 920, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. Each component communicates with each other using different buses and can be installed on a common motherboard or installed in other ways as needed. The processor can process instructions executed within the computer device, including instructions stored in the memory or on the memory to display graphical information of the GUI on an external input / output device (such as a display device coupled to the interface). In some alternative embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories. Similarly, multiple computer devices can be connected, and each device provides some necessary operations (for example, as a server array, a set of blade servers, or a multi-processor system). Figure 9 In [the figure], a processor 910 is taken as an example.

[0170] The processor 910 can be a central processing unit, a network processor, or a combination thereof. Among them, the processor 910 can further include a hardware chip. The above hardware chip can be an application-specific integrated circuit, a programmable logic device, or a combination thereof. The above programmable logic device can be a complex programmable logic device, a field programmable gate array, a generic array logic, or any combination thereof.

[0171] Among them, the memory 920 stores instructions executable by at least one processor 910, so that at least one processor 910 executes the method shown in the above embodiments.

[0172] The memory 920 can include a program storage area and a data storage area. Among them, the program storage area can store an operating system and application programs required for at least one function; the data storage area can store data created according to the use of the computer device, etc. In addition, the memory 920 can include high-speed random access memory and can also include non-transitory memory, such as at least one disk storage device, a flash memory device, or other non-transitory solid-state storage devices. In some alternative embodiments, the memory 920 can optionally include a memory remotely set relative to the processor 910, and these remote memories can be connected to the computer device through a network. Examples of the above network include but are not limited to the Internet, an enterprise intranet, a local area network, a mobile communication network, and combinations thereof.

[0173] The memory 920 can include volatile memory, such as random access memory; the memory can also include non-volatile memory, such as flash memory, a hard disk, or a solid-state drive; the memory 920 can also include a combination of the above types of memory.

[0174] The computer device further includes an input device 930 and an output device 940. The processor 910, the memory 920, the input device 930, and the output device 940 may be connected via a bus or other means. Figure 9 Taking connection via a bus as an example.

[0175] The input device 930 can receive input digital or character information, and generate key signal inputs related to the user settings and function controls of the computer device, such as a touch screen, a keypad, a mouse, a trackpad, a touchpad, a pointing stick, one or more mouse buttons, a trackball, a joystick, etc. The output device 940 may include a display device, an auxiliary lighting device (e.g., an LED), and a haptic feedback device (e.g., a vibration motor), etc. The above display device includes, but is not limited to, a liquid crystal display, a light-emitting diode, a display, and a plasma display. In some alternative embodiments, the display device may be a touch screen.

[0176] The embodiment of the present invention also provides a computer-readable storage medium. The method according to the embodiment of the present invention can be implemented in hardware, firmware, or be implemented as computer code that can be recorded on a storage medium, or be implemented as computer code that is originally stored in a remote storage medium or a non-transitory machine-readable storage medium and downloaded through a network and will be stored in a local storage medium, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium can be a magnetic disk, an optical disk, a read-only memory, a random access memory, a flash memory, a hard disk, or a solid-state drive, etc.; further, the storage medium can also include a combination of the above types of memories. It can be understood that a computer, a processor, a microprocessor controller, or programmable hardware includes a storage component that can store or receive software or computer code. When the software or computer code is accessed and executed by the computer, the processor, or the hardware, the method shown in the above embodiments is implemented.

[0177] A part of the present invention can be applied as a computer program product, such as computer program instructions. When executed by a computer, through the operation of the computer, the methods and / or technical solutions according to the present invention can be invoked or provided. Those skilled in the art should be able to understand that the forms in which computer program instructions exist in a computer-readable medium include, but are not limited to, source files, executable files, installation package files, etc. Correspondingly, the ways in which computer program instructions are executed by a computer include, but are not limited to: the computer directly executes the instruction, or the computer compiles the instruction and then executes the corresponding compiled program, or the computer reads and executes the instruction, or the computer reads and installs the instruction and then executes the corresponding installed program. Herein, the computer-readable medium can be any available computer-readable storage medium or communication medium accessible by the computer.

[0178] In the description of the present invention, it should be understood that the orientation or positional relationship indicated by the terms "center", "longitudinal", "transverse", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", "axial", "radial", "circumferential", etc. is based on the orientation or positional relationship shown in the drawings. It is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation. Therefore, it should not be construed as a limitation to the present invention.

[0179] In the description of this specification, the description with reference to the terms "this embodiment", "one embodiment", "some embodiments", "examples", "specific examples", or "some examples", etc. means that the specific features, structures, materials, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described can be combined in a suitable manner in any one or more embodiments or examples. In addition, without contradiction, those skilled in the art can combine and combine the different embodiments or examples described in this specification and the features of different embodiments or examples.

[0180] In addition, the terms "first" and "second" are only used for descriptive purposes and cannot be understood as indicating or implying relative importance or implicitly specifying the quantity of the indicated technical features. Thus, the features defined with "first" and "second" may explicitly or implicitly include at least one of such features. In the description of the present invention, the meaning of "a plurality" is at least two, such as two, three, etc., unless otherwise specifically and clearly defined.

[0181] Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art can make various modifications and variations without departing from the spirit and scope of the present invention, and such modifications and variations fall within the scope defined by the present invention.

Claims

1. A method for calculating the electromagnetic field of a metal component, characterized in that, The metal component has at least one edge, and the method includes: Obtaining a discretized mathematical model of the electric field and a mathematical model of the magnetic field corresponding to the metal component, wherein both the mathematical model of the electric field and the mathematical model of the magnetic field include update coefficients; Detecting an angle formed by the metal surfaces on both sides of the target edge, wherein the target edge is the edge that is flush with the divided grid lines among the at least one edge; Determining a correction coefficient according to the angle; According to the correction coefficient, correcting the update coefficients corresponding to all the grid faces covered by the target edge to obtain a corrected mathematical model of the electric field and a corrected mathematical model of the magnetic field; Performing time-domain iteration on the corrected mathematical model of the electric field and the corrected mathematical model of the magnetic field to obtain the electromagnetic field distribution result of the metal component.

2. The method according to claim 1, wherein The determining the correction coefficient according to the angle includes: Determining the correction coefficient according to the angle and the following formula: wherein, v1 represents the correction coefficient, and α represents the angle.

3. The method according to claim 1, wherein The correcting the update coefficients corresponding to all the grid faces covered by the target edge according to the correction coefficient includes: Determining a correction amount according to the correction coefficient and the following formula: Among them, C s represents the correction amount, and v1 represents the correction coefficient; According to the correction amount, correcting the update coefficients corresponding to all the grid faces covered by the target edge.

4. The method according to claim 3, wherein The correcting the update coefficients corresponding to all the grid faces covered by the target edge according to the correction amount includes: Multiplying the correction amount and the update coefficients corresponding to all the grid faces covered by the target edge in sequence to correct the update coefficients corresponding to all the grid faces covered by the target edge.

5. The method according to any one of claims 1 to 4, characterized in that, Before the detecting the angle formed by the metal surfaces on both sides of the target edge, the method further includes: Traversing all the edges of the metal component; Determining in sequence whether the edge is flush with the grid lines, and determining the edge that is flush with the grid lines as the target edge.

6. The method according to any one of claims 1 to 4, characterized in that, The metal component is a semiconductor chip, an electronic device in a vehicle, or an electronic device in a terminal device.

7. An electromagnetic field calculation device for a metal component, characterized in that, The metal component has at least one edge, and the device includes: An obtaining module, configured to obtain a discretized mathematical model of the electric field and a mathematical model of the magnetic field corresponding to the metal component, wherein both the mathematical model of the electric field and the mathematical model of the magnetic field include update coefficients; A detecting module, configured to detect an angle formed by the metal surfaces on both sides of the target edge, wherein the target edge is the edge that is flush with the divided grid lines among the at least one edge; A determining module, configured to determine a correction coefficient according to the angle; A correcting module, configured to correct the update coefficients corresponding to all the grid faces covered by the target edge according to the correction coefficient to obtain a corrected mathematical model of the electric field and a corrected mathematical model of the magnetic field; An iteration module, configured to perform time-domain iteration on the corrected mathematical model of the electric field and the corrected mathematical model of the magnetic field to obtain the electromagnetic field distribution result of the metal component.

8. A computer device, characterized in that, Includes: A memory and a processor, which are communicatively connected to each other. The memory stores computer instructions, and the processor executes the computer instructions to execute the electromagnetic field calculation method of the metal component according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions for causing a computer to execute the electromagnetic field calculation method of the metal component according to any one of claims 1 to 6.

10. A computer program product, characterized in that, It includes computer instructions for causing a computer to execute the electromagnetic field calculation method of the metal component according to any one of claims 1 to 6.

Images