39 results about "Computational electromagnetics" patented technology

The invention belongs to the technical field of computational electromagnetics, and discloses a numerical calculation method for simulating electromagnetic wave propagation characteristics in time-varying plasma, which is a three-dimensional Newton-HIE-FDTD method for efficiently simulating propagation of electromagnetic waves in time-varying plasma based on a hybrid explicit-implicit finite difference time domain method and Newton's second law of motion. The method inherits the advantages of the Newton-ADE-FDTD method, the time-varying electron density is accurately represented through the time derivative, and the abundant dynamic characteristics of the time-varying plasma can be accurately captured. On the basis, the method also fuses the advantages of an HIE-FDTD method, so that the time step length of the algorithm is only restricted by the sizes of grids in two directions, and compared with other simulation methods based on time-varying plasma, the method provided by the invention is more excellent in calculation efficiency and is more suitable for large-scale popularization and application. The method is especially suitable for analysis and simulation of electromagnetic wave propagation characteristics in time-varying plasma.

The application provides a method for calculating interface polarization charge and current distribution based on reflection coefficient, and relates to the technical field of computational electromagnetics. The method comprises the following steps: establishing a reflection coefficient-dielectric parameter mapping, directly correlating the reflection characteristics of electromagnetic waves with the intrinsic parameters of the medium through a reflection coefficient formula, and laying a physical foundation for polarization analysis; interface electromagnetic field reconstruction, using boundary conditions and constitutive relations to determine the electric field distribution on both sides of the interface and the corresponding polarization intensity. The distribution rule of the electromagnetic field in the medium is determined, the macroscopic electric field is converted into the spatial description of the polarization intensity, and an intermediate variable is provided for charge and current calculation; polarization dynamic response quantification, through the spatial and temporal derivatives of the polarization intensity, the interface polarization charge density and the polarization current density are derived. The theoretical model is converted into actual physical quantities, and the dynamic response mechanism of charges and currents when electromagnetic waves interact with the medium is revealed.

The invention discloses a triangular mesh-based high-order electromagnetic field simulation method, and belongs to the field of computational electromagnetism. In order to solve the problems of significant numerical dispersion and phase error accumulation when an existing triangular mesh time domain finite difference method is used for processing a complex geometric structure, the method comprises the following steps: firstly, introducing a first additional loop on a basic triangular mesh, and combining with a first optimization coefficient to construct a high-order model; and a second additional loop is further introduced, and an optimization model is constructed in combination with a second optimization coefficient, so that the phase error is effectively inhibited under the coarse mesh, and high-precision electromagnetic simulation of a large-scale polygon structure is realized.

The invention relates to the technical field of computational electromagnetism, in particular to a metasurface scattering behavior analysis method and device, equipment and a storage medium. The method comprises the following steps: tracking an incident wave emitted into the metasurface based on geometrical optics; for the tracked incident wave, determining a first relationship between the equivalent surface current and the scattered field based on physical optics; determining an equivalent surface current based on the generalized boundary condition and a surface polarizability model; and determining a scattering field by combining the equivalent surface current and the first relationship by adopting Gordon integral. According to the metasurface scattering behavior analysis method provided by the embodiment of the invention, the surface polarizability model is combined with physical optics and geometrical optics, so that the scattering behavior related to the metasurface target can be effectively analyzed, and meanwhile, the calculation amount is reduced.

The embodiment of the invention provides an electromagnetic scattering calculation method based on a heterogeneous GPU cluster. The method is applied to the field of computational electromagnetism, and comprises the following steps: constructing a hierarchical bounding volume geometric acceleration structure of a target object and adding an edge index, and meanwhile, completely copying and distributing data of the acceleration structure to a video memory of each computational node in a heterogeneous GPU cluster; generating a random ray path starting from an emission source, modeling an electromagnetic scattering process as a ray path set in a probability space, and generating a ray direction through an importance sampling strategy; determining a static scheduling strategy for the ray batch based on the calculation cost of pre-sampling and geometric enhancement, and carrying out non-uniform distribution on task tiles according to the heterogeneous GPU calculation power; and monitoring the execution state of the heterogeneous GPU cluster in real time, and dynamically adjusting the distribution strategy of the residual ray tasks based on execution feedback. According to the method, the calculation overhead can be remarkably reduced while the calculation precision is ensured, and the efficiency and expandability of complex target electromagnetic scattering calculation are improved.

This invention belongs to the field of computational electromagnetics numerical simulation and provides a simulation method for electromagnetic scattering of thin-coated targets based on Schur complement preconditions. First, for the IBC scattering target to be solved, equivalent surface current and equivalent surface magnetic current are introduced as independent unknowns. The electromagnetic scattering solution equation is established using SDIE and then discretized. Next, the impedance matrix is ​​decomposed into near-field and far-field matrices using an octree structure. For the near-field matrix, key operator sub-blocks are explicitly assembled and extracted. Then, approximate operators in block diagonal form and Schur complement operators are constructed and approximately inverted. Simultaneously, the diagonal product of intermediate blocks is pre-calculated and stored based on the approximate operators. Finally, the constructed preconditions are applied to the GMRES iterative solver. When the iteration error reaches the preset convergence criterion, the current and magnetic current distributions on the target surface are output. Further calculations of radiation characteristic parameters such as the target's radar cross section are then performed, achieving low-cost, high-efficiency, and easily parallelizable electromagnetic simulation solutions.

The invention belongs to the field of computational electromagnetism numerical simulation, and particularly relates to a non-uniform conformal FDTD method based on a sub-grid technology. According to the method, non-uniform conformal geometric correction is introduced into the sub-grids, so that high-precision fitting of a curved surface boundary is realized, step approximation errors are reduced, and the calculation precision of field quantities near the boundary is improved; according to the method, bidirectional coupling transfer which simultaneously meets line integral conservation and flux conservation is established on a main-sub grid interface, and a conservation relation required by a non-uniform conformal FDTD method is prevented from being damaged, so that interface numerical value reflection is effectively inhibited, energy non-physical errors are reduced, and long-term stability is improved; and fine grids are only adopted in the local fine structure neighborhood and are matched with time sub-cycle propulsion, so that global grid encryption and limitation of global time step length to a minimum unit are avoided, and the storage and calculation time overhead is remarkably reduced while the precision is ensured. According to the method, geometric modeling precision, interface stability and calculation efficiency are considered, and the method has high engineering application value.

The invention discloses a chaff cloud cluster electromagnetic scattering calculation method based on a binary-voxel method, and belongs to the technical field of computational electromagnetism and electronic countermeasures. According to the method, an incidence voxel coordinate system and a scattering voxel coordinate system are respectively established through two times of dynamic coordinate rotation transformation, voxel division, chaff concentration statistics and attenuation factor calculation are independently executed on a bidirectional path, and accurate modeling of incidence and scattering path attenuation is realized in combination with a linear moment method basic scattering field solution. According to the method, the scattered field calculation precision of the non-uniform chaff cloud cluster in the complex electromagnetic environment is remarkably improved, meanwhile, the calculation complexity is reduced through voxel-level discretization, and the method can be applied to electronic countermeasure efficiency evaluation.

The invention discloses a design method and system of a high-temperature-resistant bearing wave-absorbing / wave-transmitting structure, and belongs to the technical field of calculation of electromagnetic structural strength crossing. Comprising the following steps: carrying out mechanical, thermal and electromagnetic performance characterization on a bearing core layer, and constructing an equivalent transmission line model of the bearing core layer; constructing an equivalent circuit model based on the equivalent transmission line model; importing the structure parameters into the equivalent circuit model, calculating the electrical performance, and researching the influence of the structure parameters on the electromagnetic performance to obtain key design parameters; establishing an optimization model by taking the key design parameters as optimization objects, the electromagnetic performance indexes as optimization targets and the force and thermal performance indexes of the structure as constraint conditions; key design parameters are initialized and imported into the optimization model, deformation and temperature performance indexes of the structure are calculated, and electromagnetic performance indexes under all temperature detection points are obtained; judging whether the thermal performance meets a constraint condition or not and whether the electromagnetic performance index meets an optimization target or not; and if not, modifying the key design parameters until the conditions are met.

The invention discloses a coated target electromagnetic scattering modeling method based on multi-core heterogeneous parallel computing, and belongs to the field of computational electromagnetism and high-performance computing. The method comprises the steps that an MPI + OpenMP + DSP three-level heterogeneous parallel framework is constructed, and optimal configuration of computing resources is achieved through one-to-one binding of an MPI process and a DSP acceleration card and one-to-one mapping of an OpenMP thread and a CPU core; dispersing a to-be-solved coating target by adopting a triangular grid, and defining an RWG primary function to expand surface current; substituting the surface current approximate expansion into an electric field integral equation based on impedance boundary conditions; carrying out impedance matrix calculation by adopting a process-based parallel filling strategy, and carrying out parallel solution on a matrix equation based on a ScaLAPACK library; and finally, calculating a far-region scattering field and a radar scattering sectional area according to the solved current coefficient. According to the method, through an innovative heterogeneous parallel architecture and a calculation task scheduling strategy, the problems of insufficient memory and low calculation efficiency during processing of an electrically large coated target in a traditional method are effectively solved, and an efficient and reliable technical means is provided for electromagnetic stealth design of a complex target.

The invention belongs to the crossing field of computational electromagnetism and artificial intelligence, and particularly provides a three-dimensional metal target electromagnetic simulation method based on a neural network and an equivalence principle. The method is used for solving the technical problems that in three-dimensional metal target electromagnetic scattering characteristic analysis, an equivalent principle algorithm is high in calculation complexity and large in memory consumption in operator matrix construction and scattering characteristic solving stages, and a traditional numerical integration method is low in efficiency. According to the method, the deep neural network is combined with multi-frequency triangular features and a space coupling perception activation mechanism, a prediction model of an equivalent principle operator matrix is constructed, high prediction precision is maintained, computing resources and time required for operator matrix construction and scattering characteristic solving are remarkably reduced, and the method is suitable for large-scale popularization and application. And an efficient, stable and extensible rapid solution scheme is provided for the electromagnetic scattering problem of the complex three-dimensional metal target.

The application discloses a three-dimensional space separation heterogeneous region scattering field simulation method combined with a neural network, and belongs to the technical field of computational electromagnetics and artificial intelligence. The method divides a solution domain into a background domain and an observation domain, and the background domain and the observation domain are separated from each other in space. Cross-domain numerical Green functions and in-domain numerical Green functions are extracted through virtual probe scanning. Physical characteristic vectors corresponding to field source relations are extracted, and training sample sets are constructed in combination with the numerical Green functions. Neural network models are trained by using the training sample sets to respectively fit the cross-domain and in-domain numerical Green functions. The trained neural network models are used to reconstruct numerical Green functions corresponding to to-be-solved foreground scatterers in real time, integral equation solutions of equivalent currents of the foreground scatterers are established, and three-dimensional heterogeneous region scattering fields are calculated in the observation domain. The method can realize fast and accurate calculation of numerical Green functions in a three-dimensional heterogeneous region scattering scene.

The invention relates to the cross technical field of computational electromagnetism and precision permanent magnet device design, in particular to a method for improving the nonuniformity of a magnetic disk space magnetic field, and aims to solve the problems that a method for improving the uniformity of a magnetic field in the prior art is low in optimization efficiency and difficult to break through the uniformity bottleneck caused by angular harmonic waves and a structural coupling effect. According to the scheme, the method comprises the steps that an accurate finite element model containing seven circles of fan-shaped magnets, a parallel magnetic disk, an eddy current disk and an iron support is established in Opera3D, and meshing and solving are carried out; in the first stage, the structural height of the iron stand is optimized through parameterized scanning, and the optimal height parameter is determined; and in the second stage, on the basis of the optimal iron stand height, the outermost two circles of the magnetic disk are set as NdFeB. The non-uniformity deviation (delta B) of the magnetic field of the target spherical area can be reduced by 13%, and the method is suitable for high-precision scenes such as nuclear magnetic resonance and synchrotron radiation light sources and has the advantages of being efficient in design and high in repeatability.

The application belongs to the technical field of computational electromagnetics, and discloses a numerical calculation method for simulating electromagnetic wave propagation characteristics in time-varying plasma, which is a three-dimensional Newton-HIE-FDTD method for efficiently simulating electromagnetic wave propagation in time-varying plasma based on a mixed explicit-implicit finite-difference time-domain method and Newton's second law of motion. The method inherits the advantages of the Newton-ADE-FDTD method, accurately represents the time-varying electron density through the time derivative, and can accurately capture the rich dynamic characteristics of the time-varying plasma. On this basis, the method also integrates the advantages of the HIE-FDTD method, so that the time step of the algorithm is only restricted by the grid size in two directions. Compared with other simulation methods based on time-varying plasma, the method is more efficient in calculation, and is particularly suitable for the analysis and simulation of electromagnetic wave propagation characteristics in time-varying plasma.

The invention discloses an electromagnetic field prediction method based on a deep learning network model, belongs to the technical field of calculation of absorption boundaries in electromagnetics, constructs and trains an Autoformer model, applies a sequence decomposition mechanism and a sparse autocorrelation mechanism to a WCS-FDTD framework based on an electromagnetic model so as to update electromagnetic field data more quickly and accurately, and improves the accuracy of the electromagnetic field prediction. And inputting the obtained historical internal field domain and the electromagnetic field data at the corresponding absorption boundary into a trained neural network model for prediction, and calculating the electromagnetic field data of the internal field domain at the next moment based on the predicted electromagnetic field data at the absorption boundary at the next moment and the electromagnetic field data of the internal field domain at the current moment. And combining the obtained predicted value of the electromagnetic field data at the absorption boundary corresponding to the next moment with the electromagnetic field data of the internal field domain at the next moment, generating input data at the next moment, carrying out iterative prediction, and finally obtaining the electromagnetic field data of the computational domain at each moment.

The invention discloses an unconditionally stable UPML electromagnetic simulation method based on Dosin FDTD, and belongs to the field of computational electromagnetics and numerical analysis, and the method comprises the steps: introducing an anisotropic damping term under a Dosin finite difference time domain (MS-FDTD) framework, embedding the UPML absorption characteristic into a Dosin Maxwell system, and carrying out the UPML electromagnetic simulation. And the consistency of the main domain and the absorption layer in geometric structure and energy transfer is realized. By adopting a Preissmann implicit discrete format and a local one-dimensional operator splitting strategy, a numerical value updating form irrelevant to a time step and a space step is obtained, so that unconditional stable solution is realized. According to the method, the electromagnetic field quantity is decoupled through the intermediate variable of the electric displacement vector and the magnetic induction intensity, and the solving efficiency and the numerical precision are improved. The provided Dosin UPML method can keep energy conservation and structural stability in complex media and long-time simulation, and the electromagnetic wave absorption precision and simulation reliability are remarkably improved.

This invention belongs to the field of numerical simulation in computational electromagnetics, specifically a non-uniform conformal FDTD method based on submesh technology. This invention introduces non-uniform conformal geometric corrections within the submesh to achieve high-precision fitting of curved surface boundaries, reducing step approximation errors and improving the accuracy of field calculations near the boundaries. At the master-submesh interface, a bidirectional coupling transfer is established that simultaneously satisfies line integral conservation and flux conservation, avoiding disruption of the conservation relationships required by the non-uniform conformal FDTD method. This effectively suppresses numerical reflections at the interface, reduces energy non-physical errors, and improves long-term stability. Furthermore, by using fine meshes only in the neighborhood of local microstructures and coordinating with time sub-cycle progression, global mesh refinement and the global time step being limited by the smallest unit are avoided, significantly reducing storage and computation time overhead while maintaining accuracy. This invention balances geometric modeling accuracy, interface stability, and computational efficiency, and has high engineering application value.

The present application belongs to the technical field of block iterative matrix solution in computational electromagnetics, in particular to a block iterative matrix solution method based on matrix dimension reduction and preprocessing technology. The present application uses two technologies of matrix dimension reduction and matrix preprocessing: the right end item matrix is processed by using the matrix dimension reduction, through low rank approximation, which can effectively ensure the accuracy of the solution, significantly reduce the calculation cost, and facilitate the processing of multiple right end linear systems; and the BCOCG algorithm is applied to the preprocessing technology to further improve the matrix behavior and improve the convergence and stability of the algorithm. For the system equation generated by multiple physical excitation sources or multiple scanning angles in electromagnetic problems, the memory cost is shortened from the perspective of short recursive Krylov subspace iteration method, and the block iteration BCOCG algorithm can calculate multiple column right end item vectors at the same time, avoiding the matrix ill-conditioning problem caused by the correlation between columns and columns.

The application discloses a grid generation method for electromagnetic environment simulation on a water surface platform and belongs to the technical field of computational electromagnetics, which comprises the following steps: determining a minimum rectangular area capable of completely covering the shape of a ship platform, and performing segmented function fitting on the outer contour of the ship platform; based on the minimum rectangular area and the segmented outer contour, the number and position coordinates of grid points along the transverse and longitudinal directions of the ship platform are obtained according to a grid division rule, and each grid is indexed by a coordinate number; in the calculation, the grids are called in turn according to the coordinate number, when a certain target grid needs to be used, the position information and area information of the target grid are calculated according to the coordinate number corresponding to the target grid and the grid division rule; and after the current target grid is used up, the grid information of the next coordinate number is covered. The application solves the problem of slow calculation progress or difficulty in calculation caused by the grid generation, storage and calling links in the electromagnetic calculation of a super-electric large-size target.

The invention relates to an unmanned aerial vehicle electromagnetic scattering characteristic solving method and device based on a graph neural network, and belongs to the technical field of computational electromagnetism and radar target recognition, and the method specifically comprises the steps: carrying out the initialization of a system and the discretization of an unmanned aerial vehicle grid; constructing graph structure data based on an RWG primary function, and establishing a mapping relation between an unstructured grid of a physical space and graph data which can be processed by a graph neural network; a GraphSolver model of a graph neural network solver is constructed, and an end-to-end deep neural network model is formed; training and optimizing the model; according to the method, the problems of too long calculation time and huge memory consumption when an existing full-wave numerical calculation method is used for analyzing complex targets such as an unmanned aerial vehicle are solved; and the existing deep learning method is difficult to accurately process the non-uniform triangular mesh of the unmanned aerial vehicle and is easy to lose key geometric details such as wing edges.

The present application belongs to the field of electromagnetic numerical simulation, and particularly relates to a self-adaptive size field construction method for high-curvature model grid generation. The present application generates uniform sampling points and constructs an initial background grid by performing geometric processing on a target geometric model; and calculates discrete central axis points on the basis of the initial background grid, obtains adjacent self-adaptive sizes, calculates model edge curvatures, obtains curve curvature self-adaptive sizes, and fuses the two to construct a self-adaptive size field. On this basis, internal sampling points of a high-curvature region are generated, a high-curvature self-adaptive size is constructed, and the original self-adaptive size field is fused; and finally, a high-quality unstructured self-adaptive grid is generated. The present application improves grid generation efficiency and simulation precision, and provides a basis for high-precision algorithm engineering applications of computational fluid dynamics and computational electromagnetics in the fields of automobile manufacturing, civil engineering, environmental engineering, shipbuilding industry and aviation industry.

This invention belongs to the field of computational electromagnetics, specifically a method for transforming a one-dimensional axial magnetic field into a three-dimensional axisymmetric magnetic field. This invention uses two variable-step-size, four-point numerical differential calculations followed by averaging to approximate the first derivative at the point to be calculated. This avoids the problem that calculating four points once might not improve accuracy, and also avoids the overly complex formula of directly using five points. Furthermore, by using the approximate value of the first derivative at known coordinates combined with Spline curve fitting to estimate the approximate value of the first derivative at the point z to be calculated, the problem of error superposition is effectively avoided, resulting in a higher accuracy for the radial component B of the calculated three-dimensional axisymmetric magnetic field. r The results are closer to the theoretical values. This invention effectively solves the problems of low accuracy and error accumulation in the transformation from one-dimensional axial magnetic field to three-dimensional magnetic field, and provides strong support for magnetic field loading in the design of devices such as vacuum electronic devices, ion sources, and particle accelerators.

Various examples related to microwave dielectric analyzers and their use are provided. In one example, a microwave dielectric analyzer includes a measurement apparatus having a conductive electrode that can couple to a microwave analyzer and processing circuitry that can determine a dielectric characteristic of the dielectric specimen using a reflection coefficient measured by the microwave analyzer. The dielectric characteristic can be determined using a computational electromagnetic model of the measurement apparatus. The reflection coefficient can be measured by the microwave analyzer with the dielectric specimen in contact with the conductive electrode and / or sandwiched between conductive electrodes. The conductive electrodes can be axially aligned, and the second electrode may not be coupled to the microwave analyzer.

The invention relates to a physical optical scattered field calculation method based on pixelated projection, and belongs to the field of computational electromagnetics, and the calculation method comprises the steps: obtaining a three-dimensional CAD model of a target, then carrying out the corresponding setting of parameters, and carrying out the curved surface projection of the target, and generating a projection region; pixel discretization and curvature constraint are carried out on the projection area, a radiation function of pixel points is calculated, and a radiation image is constructed. And calculating an oscillation integral of each pixel point through a closed expression, calculating through a GPU parallel architecture to obtain a total scattering field, and outputting the total scattering field. Based on the basic theory of physical optics (PO), the method can ensure high calculation precision on the premise of meeting PO application conditions. Geometric details and surface characteristics of the target can be described in a pixelated projection mode, so that a complex shape, an edge effect and scattering behaviors of a non-uniform material surface are simulated more accurately, and the applicability of the method to an actual engineering target and the reliability of a calculation result are improved.

The present application belongs to the field of computational electromagnetics, and provides an electromagnetic simulation isosurface extraction method based on space rearrangement and deep learning prediction, which is used to solve the problems of waste of computing resources, memory bottleneck of unstructured grid data and insufficient real-time performance in the prior art. First, the unstructured grid data output by electromagnetic simulation is rearranged in memory by using Morden code, the field points that are adjacent in geometry are forced to be continuous in physical memory, random memory access is converted into streaming memory access, and the bandwidth utilization is greatly improved. Then, a lightweight fully connected neural network is constructed, the field scalar values of the grid cell vertices and the target isosurface are input, and the probability of the grid cell crossing the isosurface is output. The trained network model is used for fast batch reasoning on the rearranged grid data, most of the background space is removed at a very low computing cost, and an active cell candidate table is generated. Finally, accurate geometric interpolation is performed on the active cells to generate high-quality isosurface grids.

The application discloses a target echo classification method and system based on multi-modal data fusion, relates to the technical field of computational electromagnetics, and solves the technical problems of classification accuracy reduction of a model in a noise environment, limited generalization ability, and large performance fluctuation when the distribution of modal data changes; the weight of each feature is analyzed when modeling the time sequence dependence of BiLSTM, global information is reserved through cross-modal fusion after the features are removed, cross fusion focuses on key discriminant features after the features are removed, and both can take into account information integrity and calculation efficiency. The contribution of different modal features to the classification result is dynamically adjusted through the hidden state of GRU. If some modal data is missing, the original fusion path can still use other modal information, and the screening fusion path can depend on the key features of the remaining modal. The balance between redundant feature reservation and key feature screening can reduce the sensitivity of the model to specific noise patterns.

The application discloses a PIN diode modeling and simulation method based on LP-CDI-FDTD, belongs to the field of computational electromagnetics, microwave technology and circuit simulation, and comprises the following steps: firstly, an equivalent lumped parameter model of a semiconductor device is established, and a voltage-current relationship thereof is obtained; then, the relationship is converted into an equivalent current density, and is embedded into a Maxwell curl equation according to the relationship between an electric field and a voltage, so as to construct an electromagnetic-circuit coupling equation set; further, a leapfrog compliant divergence unconditional stable time domain finite difference method is adopted to discretize the coupling equation set, and a discrete updating equation is obtained; finally, iterative calculation is carried out according to the equation, and electromagnetic fields and device states are synchronously updated. The method realizes the inherent strong coupling between full-wave electromagnetic fields and semiconductor physical processes, and breaks through the time step restriction caused by a CFL condition in traditional methods by using the unconditional stability of the algorithm, so that the calculation efficiency can be significantly improved when simulating a device containing a fine structure.

The application discloses a method and device for calculating nominal electric field of direct current transmission based on a firefly algorithm, and relates to the technical field of high-voltage engineering and computational electromagnetics. The method comprises the following steps: obtaining transmission line parameters, setting the number of simulation electric charges, the number and coordinates of matching points and verification points for each conductor, and calling a preset simulation electric charge method calculation framework to solve simulation electric charge quantity vectors and target functions; encoding all simulation electric charge coordinates into optimization variables, randomly generating an initial population and calculating brightness through the firefly algorithm, iteratively optimizing simulation electric charge positions, and finally obtaining a terminal generation optimal simulation electric charge position vector, which is substituted into the calculation framework to calculate a terminal generation optimal simulation electric charge quantity vector, and then the ground nominal electric field is solved. The optimal position of the simulation electric charge is searched through the firefly algorithm, the calculation precision of the ground electric field of the direct current transmission line is improved, and reliable technical support is provided for the ground nominal electric field evaluation under the static electric field condition of the direct current transmission line.

This invention discloses a method and system for integrated load-bearing and wave-transmission optimization design of array antenna windows, relating to the field of computational electromagnetics structural strength cross-section and optimization technology. The invention includes: S1: Constructing a load-bearing and wave-transmission optimization model for the array antenna window, and obtaining initial sample points using the Latin hypercube sampling method; S2: Obtaining the actual response values ​​of the initial sample points using an integrated load-bearing and wave-transmission simulation analysis method for the array antenna window, and constructing a load-bearing and wave-transmission surrogate optimization model for the array antenna window using a Kriging surrogate model. This invention can solve the problems of high time cost and the inability to simultaneously optimize load-bearing / wave-transmission performance during the optimization design of array antenna windows, obtaining an array antenna window structure with optimal wave-transmission performance under load constraints with less time cost, thus reducing the optimization design cycle and cost of array antenna windows.