Multidisciplinary coupling surface grid point communication data coding method, decoding method, device and equipment

Abstract

The invention discloses a multidisciplinary coupling surface grid point communication data encoding method, decoding method, device and equipment, and relates to the field of multidisciplinary coupling numerical simulation, and the method comprises the steps: encoding various types of coupling data on grid nodes on a plurality of coupling interfaces into a one-dimensional real array, all coupling data communication can be completed through one-time message transmission interface calling, the parallel communication calling frequency is reduced, and the multidisciplinary coupling simulation efficiency is improved.

Description

Technical Field

[0001] The present invention relates to the field of multidisciplinary coupling numerical simulation, and in particular to a multidisciplinary coupling surface grid point communication data encoding method, decoding method, device and equipment. Background Art

[0002] At present, multidisciplinary coupling simulation is mainly achieved through two ways: one is to solve multidisciplinary coupling in a unified method framework, that is, different disciplines use the same grid discretization and the same calculation method; the other is to adopt a separate solution strategy, that is, different disciplines solve independently, and transfer physical information and data between different disciplines through interfaces. The first method has the advantages of high accuracy and good stability, but due to the limitations of the method itself, it is difficult to adapt to complex application scenarios, and the design and development are difficult. It is rarely used in the field of industrial design and is mainly used for academic theoretical research. Although the second method is slightly lower than the first method in accuracy and robustness, it can use a relatively mature single-discipline solver, thereby greatly reducing the difficulty of solver development, and can be applied to extremely complex application scenarios. It is a multidisciplinary coupling simulation method widely used in the current industrial field.

[0003] In the separated coupling solution strategy, different subject solvers run independently, and the iterative update of different physical field data is achieved through the exchange of interface information. When the simulation problem is more complex, the subject solution adopts a multi-process parallel operation mode, which requires parallel communication and transmission of coupling interface data between different processes. At present, the message passing interface is still the main way of parallel data communication in multidisciplinary coupling analysis, but the message passing interface method can only send and receive a single type of data at a time. When there are many types of coupling data, the message passing interface needs to be called multiple times, so that the solver needs to carry out multiple communication interface preparations, resulting in an increase in the overall communication time. When the number of calculation processes is large, the communication time cost will occupy a large proportion of the overall calculation time, thereby reducing the efficiency of multidisciplinary coupling simulation analysis. Summary of the invention

[0004] In view of this, the purpose of the present invention is to provide a multidisciplinary coupling surface grid point communication data encoding method, decoding method, device and equipment, which can communicate multiple types of coupling data once, reduce the number of transmission interface calls, and improve the efficiency of multidisciplinary coupling analysis. The specific scheme is as follows:

[0005] In a first aspect, the present application discloses a multidisciplinary coupling surface grid point communication data encoding method, comprising:

[0006] Determine the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem;

[0007] According to the multidisciplinary coupling analysis problem, a corresponding discipline analysis calculation grid is formed to extract a coupling interface grid from the discipline analysis calculation grid, and the number of coupling interfaces and the number of grid nodes of the coupling interfaces are calculated;

[0008] Performing numerical calculation on the physical quantity to be communicated in the grid nodes of the coupling interface to obtain the value of the physical quantity to be communicated of each of the grid nodes;

[0009] Convert the current number of coupling interfaces, the number of grid nodes of a single coupling interface, and the physical field dimension in the form of integer data into the target number of coupling interfaces, the target number of grid nodes, and the target physical field dimension in the form of real data, respectively, and determine the physical quantity encoding order;

[0010] The target physical field dimension and the target coupling interface quantity are encoded to obtain a first real array, and the target grid node quantity is encoded, and the encoded information is filled into the first real array to obtain a second real array;

[0011] The physical quantity to be communicated is cyclically encoded according to the physical quantity encoding order, and the encoding result is filled into the second real type array to obtain a target real type array.

[0012] Optionally, before determining the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem, the method further includes:

[0013] Generate a calculation example folder for a multidisciplinary coupling analysis problem; wherein the calculation example folder is used to store the folders of each discipline solver and the corresponding coupling framework configuration file;

[0014] When the multidisciplinary coupling analysis problem is obtained, the example folder is called to simulate the multidisciplinary coupling analysis problem to obtain a coupling interface involved in the coupling analysis.

[0015] Optionally, the step of determining the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem includes:

[0016] According to the multidisciplinary coupling analysis problem and the mechanism of iterative exchange between different disciplines, the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated participating in the coupling communication in the multidisciplinary coupling analysis problem are determined.

[0017] Optionally, calculate the number of mesh nodes at the coupling interface, including:

[0018] If there are multiple mesh surfaces in the coupling interface, the number of mesh points of each mesh surface is counted respectively, and the number of mesh points of all the mesh surfaces is cumulatively calculated to obtain the number of mesh nodes of the coupling interface.

[0019] Optionally, the performing numerical calculation on the physical quantity to be communicated in the grid nodes of the coupling interface to obtain the value of the physical quantity to be communicated of each grid node includes:

[0020] Calculating the multidisciplinary coupling analysis problem according to the solvers of each discipline to obtain physical field distribution data;

[0021] The physical field distribution data and the mesh data of the coupling interface are processed by a preset data processing technology to calculate the physical quantity values ​​to be communicated on the mesh nodes of the coupling interface.

[0022] Optionally, the target real array is a one-dimensional coupled data array.

[0023] In a second aspect, the present application discloses a method for decoding communication data of a multidisciplinary coupling surface grid point, comprising:

[0024] Obtain a target real array encoded using the multidisciplinary coupling surface grid point communication data encoding method;

[0025] Extracting target physical field dimensions, target coupling interface quantity, and target grid node quantity from the target real array, and converting the target physical field dimensions, target coupling interface quantity, and target grid node quantity into data format to obtain the physical field dimensions, coupling interface quantity, and grid node quantity in the form of integer data;

[0026] Corresponding communication data is obtained from the encoded physical quantity according to the physical quantity encoding order, the physical field dimension, the number of coupling interfaces and the number of grid nodes.

[0027] Optionally, after obtaining corresponding communication data from the encoded physical quantity according to the physical quantity encoding order, the physical field dimension, the number of coupling interfaces and the number of grid nodes, the method further includes:

[0028] Using the current communication data decoded by any subject solver as control point information, and using the control point information and the initial interpolation function to construct a corresponding target interpolation function;

[0029] The target interpolation function is used to perform interpolation calculation on other communication data except the current communication data to obtain interpolated communication data.

[0030] In a third aspect, the present application discloses a multidisciplinary coupling surface grid point communication data encoding device, comprising:

[0031] An information determination module is used to determine the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem;

[0032] A data calculation module, used to form a corresponding subject analysis calculation grid according to the multidisciplinary coupling analysis problem, to extract a coupling interface grid from the subject analysis calculation grid, and to calculate the number of coupling interfaces and the number of grid nodes of the coupling interfaces;

[0033] A numerical calculation module, used for performing numerical calculation on the physical quantity to be communicated in the grid nodes of the coupling interface to obtain the value of the physical quantity to be communicated of each grid node;

[0034] A form conversion module, used to convert the number of coupling interfaces, the number of grid nodes of a single coupling interface, and the physical field dimension in the current integer data form into the target number of coupling interfaces, the target number of grid nodes, and the target physical field dimension in the real data form, and determine the physical quantity coding order;

[0035] A first encoding module is used to encode the target physical field dimension and the target coupling interface quantity to obtain a first real array, and encode the target grid node quantity, and fill the encoded information into the first real array to obtain a second real array;

[0036] The second encoding module is used to perform cyclic encoding on the physical quantity to be communicated according to the physical quantity encoding order, and fill the encoding result into the second real type array to obtain a target real type array.

[0037] In a fourth aspect, the present application discloses an electronic device, including:

[0038] Memory, used to store computer programs;

[0039] The processor is used to execute the computer program to implement the steps of the aforementioned disclosed method.

[0040] It can be seen that the present application discloses a method for encoding communication data of multidisciplinary coupling surface grid points, including: determining the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem; forming a corresponding discipline analysis calculation grid according to the multidisciplinary coupling analysis problem, extracting the coupling interface grid from the discipline analysis calculation grid, and calculating the number of coupling interfaces and the number of grid nodes of the coupling interface; performing numerical calculations on the physical quantity to be communicated in the grid nodes of the coupling interface to obtain the value of the physical quantity to be communicated of each of the grid nodes; converting the number of coupling interfaces, single and double coupling interfaces in the current integer data form into the physical quantity data type, physical field dimension, and physical quantity to be communicated; and finally, calculating the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data in the multidisciplinary coupling analysis problem ... The number of grid nodes of a coupling interface and the physical field dimension are converted into the target coupling interface number, target grid node number and target physical field dimension in the form of real data, and the physical quantity coding order is determined; the target physical field dimension and the target coupling interface number are encoded to obtain a first real array, and the target grid node number is encoded, and the encoded information is filled into the first real array to obtain a second real array; the physical quantity to be communicated is cyclically encoded according to the physical quantity coding order, and the encoding result is filled into the second real array to obtain the target real array. It can be seen that by uniformly converting the key information such as the number of coupling interfaces in the form of integer data, the number of grid nodes of a single coupling interface, and the physical field dimension into the form of real data, the data type is more consistent in the subsequent coding process, avoiding additional processing steps and potential errors that may be caused by inconsistent data types, thereby improving the overall efficiency of data processing. Encoding according to the physical quantity coding order reduces confusion and uncertainty in the coding process, improves the accuracy and speed of coding, and thus improves the data processing efficiency of the entire multidisciplinary coupling analysis. After a series of encoding processes, the target real array is obtained. This method of integrating multiple types of physical quantities to be communicated into a unified format (target real array) is conducive to data communication between different discipline solvers in multidisciplinary coupling analysis. It enables one data transfer to carry all relevant coupling data, avoiding the communication overhead and time delay that may be caused by multiple transfers of different types of data, thereby enhancing the effect of data communication and improving the overall operating efficiency of multidisciplinary coupling analysis. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying creative work.

[0042] Figure 1A flow chart of a multidisciplinary coupling surface grid point communication data encoding method disclosed in this application;

[0043] Figure 2 A fluid solver coding format diagram disclosed in this application;

[0044] Figure 3 A structural solver coding format diagram disclosed in this application;

[0045] Figure 4 A flow chart of a multidisciplinary coupling surface grid point communication data decoding method disclosed in this application;

[0046] Figure 5 A flowchart of wing flutter analysis based on the encoding and decoding method disclosed in this application;

[0047] Figure 6 This is a schematic diagram of the structure of a multidisciplinary coupling surface grid point communication data encoding device disclosed in this application;

[0048] Figure 7 This is a structural diagram of an electronic device disclosed in this application. DETAILED DESCRIPTION

[0049] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0050] With the rapid improvement of computer hardware capabilities and the continuous development of numerical calculation methods, CAE (Computer-Aided Engineering) industrial software based on numerical calculations has gradually become an important driving force for the transformation of design patterns in the equipment manufacturing field represented by aviation and aerospace, and plays an important role in product design, performance evaluation, and use and maintenance. Among them, the development of single-discipline simulation software has matured and can replace physical simulation or physical experiments in many fields, providing high-precision prediction data for product design. However, with the continuous improvement of the sophistication and requirements of modern product design, many cross-disciplinary coupling problems have emerged, which have become a key problem restricting product performance. The improvement of modern product quality and performance urgently needs the support of high-reliability multidisciplinary coupling simulation analysis tools.

[0051] At present, multidisciplinary coupling simulation is mainly achieved through two ways: one is to solve multidisciplinary coupling in a unified method framework, that is, different disciplines use the same grid discretization and the same calculation method; the other is to adopt a separate solution strategy, that is, different disciplines solve independently, and transfer physical information and data between different disciplines through interfaces. The first method has the advantages of high accuracy and good stability, but due to the limitations of the method itself, it is difficult to adapt to complex application scenarios, and the design and development are difficult. It is rarely used in the field of industrial design and is mainly used for academic theoretical research. Although the second method is slightly lower than the first method in accuracy and robustness, it can use a relatively mature single-discipline solver, thereby greatly reducing the difficulty of solver development, and can be applied to extremely complex application scenarios. It is a multidisciplinary coupling simulation method widely used in the current industrial field.

[0052] In the separated coupling solution strategy, different subject solvers run independently, and the iterative update of different physical field data is achieved through the exchange of interface information. When the simulation problem is more complex, the subject solution adopts a multi-process parallel operation mode, which requires parallel communication and transmission of coupling interface data between different processes. At present, the message passing interface is still the main way of parallel data communication in multidisciplinary coupling analysis, but the message passing interface method can only send and receive a single type of data at a time. When there are many types of coupling data, the message passing interface needs to be called multiple times, so that the solver needs to carry out multiple communication interface preparations, resulting in an increase in the overall communication time. When the number of calculation processes is large, the communication time cost will occupy a large proportion of the overall calculation time, thereby reducing the efficiency of multidisciplinary coupling simulation analysis.

[0053] To this end, the present invention provides a multidisciplinary coupling surface grid point communication data encoding and decoding scheme, which can communicate multiple types of coupling data once, reduce the number of transmission interface calls, and improve the efficiency of multidisciplinary coupling analysis.

[0054] Reference Figure 1 As shown, an embodiment of the present invention discloses a multidisciplinary coupling surface grid point communication data encoding method, comprising:

[0055] Step S11: Determine the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem.

[0056] In this embodiment, before determining the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data in the multidisciplinary coupling analysis problem, it also includes: generating a calculation example folder for the multidisciplinary coupling analysis problem; wherein the calculation example folder is used to store the folders of each discipline solver and the corresponding coupling framework configuration file; when the multidisciplinary coupling analysis problem is obtained, the calculation example folder is called to simulate the multidisciplinary coupling analysis problem to obtain the coupling interface involved in the coupling analysis. It can be understood that the multidisciplinary coupling analysis environment is pre-built. The present invention adopts TRIP fluid solver, modal structure solver, and coupling interpolation algorithm. First, a calculation example folder for wing aerodynamic / structural coupling analysis is generated for the wing aerodynamic coupling analysis problem, and the calculation example folder stores folders such as fluid solver, structure solver, coupling data interpolation module, and configuration files of the coupling framework. The solver folders store the grid files and configuration files of the wing respectively, and the coupling interpolation folder stores the interpolation parameter setting file. For the fluid solver settings, first use the mesh generation software to generate the mesh of the wing's external flow field, set the wing surface as the interface participating in the coupling analysis, and output the flow field mesh in Plot3D format. Set the flow field parameters (angle of attack, Mach number, Reynolds number, etc.), boundary conditions (far field, solid wall, coupling boundary, etc.), mesh deformation method and solution method (time / space discrete format, turbulence model, unsteady iteration steps, time step, etc.) in the fluid solver configuration file. For the structural solver settings, use the structural pre-processing software, input the wing model and structural physical parameters, use the modal analysis function, set the wing surface as the coupling interface, and export the first four eigenfrequencies and modal vibration shape coordinates of the wing after the calculation. Set the interpolation method to IPS, interpolation radius and other parameters in the coupling interpolation algorithm configuration file. Design the coupling analysis type, time step, maximum number of iterations, total number of parallel cores, names of all solvers involved in the coupling analysis, number of cores and configuration file address and other parameters in the coupling framework configuration file.

[0057] In this embodiment, the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated that participate in the coupling communication in the multidisciplinary coupling analysis problem are determined according to the multidisciplinary coupling analysis problem and the mechanism of iterative exchange between different disciplines. It can be understood that all data types participating in the coupling communication in the multidisciplinary coupling analysis are determined according to the multidisciplinary coupling analysis problem and the mechanism of iterative exchange between different disciplines. For example, in fluid-solid coupling analysis, the fluid needs to transfer the calculated force on the coupling interface to the solid to facilitate the solid to carry out static or dynamic analysis, and the solid needs to transfer the calculated deformation displacement to the fluid to facilitate the fluid to update its shape and calculate the force in the new state. Therefore, the physical field dimensions participating in the coupling communication in fluid-solid coupling analysis , physical quantity data types include: grid node coordinates , Node aerodynamics , node displacement , and based on the number of coupled interface mesh surfaces , the number of coupling interface mesh nodes The physical quantity to be communicated is determined.

[0058] Specifically, in the wing aerodynamic coupling analysis problem, after the fluid solver completes the unsteady solution within a single coupling time step, the node aerodynamic interpolation on the coupling interface is passed to the structural solver for dynamic analysis; after the structural solver completes the solution of a single coupling time step, the displacement interpolation on the coupling interface is passed to the fluid solver for fluid mesh deformation. In the data transmission process, it is necessary to predetermine the physical quantity data types of the coupling data to be communicated between the fluid solver and the structural solver. Among them, the physical quantity types of the fluid solver participating in the coupling communication include: node coordinates and node aerodynamic forces; the physical quantity types of the structural solver participating in the coupling communication include: node coordinates and node displacements. Furthermore, since the dimensions of different physical fields determine the number of components of physical quantities, for example, in a three-dimensional physical field ( ), physical quantities such as node coordinates, node aerodynamic forces, and node displacements have three components (e.g. direction), that is, then the corresponding grid node coordinates of the physical quantity to be communicated are , node aerodynamic force , node displacement . In the two-dimensional physical field ( ), such as the plate coupling problem, since the thickness direction is much smaller than the other two directions, only the coordinates of the length and width directions are used to construct the transfer matrix, and the physical quantities only use the thickness direction, i.e., the plane normal, for coupling calculations. These physical quantities have only two components, i.e. Determining the physical field dimension ensures that these physical quantities can be correctly operated during data processing. Taking the value of the physical quantity to be communicated on the coupling interface mesh node as an example, the physical field dimension determines how to calculate the communicated physical quantity based on the physical field distribution data calculated by the subject solver. If In the pressure data processing obtained by fluid calculation, it is necessary to consider the distribution of pressure in three-dimensional space. By calculating the normal vector and area of ​​the grid unit, combined with the pressure data, the size and direction of the pressure load in three-dimensional space are obtained, and then it is reasonably distributed to the nodes. Therefore, the physical field dimensions of the computer involved in the coupled communication are further determined. In the fluid solver, the wing fluid grid belongs to the three-dimensional structured grid, so the physical field dimension data in the fluid solver is set. In the structural solver, although the wing structure finite element model is discretized using two-dimensional quadrilateral shell elements, the node coordinates of each node of the shell element and the components of the node aerodynamic force are still 3. When encoding and decoding, the coupling data of each node is operated according to the number of components, so the physical quantity dimension data is still set. is 3. And the node aerodynamic data and node displacement data of the current time step on each coupling interface that participate in the coupling communication are determined as the physical quantities to be communicated. It should be noted that currently only the physical quantities to be communicated are determined, and the specific values ​​of the physical quantities to be communicated cannot be determined.

[0059] Step S12: forming a corresponding subject analysis calculation grid according to the multidisciplinary coupling analysis problem, extracting a coupling interface grid from the subject analysis calculation grid, and calculating the number of coupling interfaces and the number of grid nodes of the coupling interfaces.

[0060] In this embodiment, according to the definition or marking of the coupling interface in each discipline solver in the multidisciplinary coupling analysis, the coupling interface grid is extracted from the discipline analysis calculation grid, and the number of coupling interfaces is calculated. , calculate the number of mesh nodes on the coupling interface . If there are multiple mesh surfaces in the coupling interface, the number of mesh points on each mesh surface is counted respectively, and the number of mesh points on all mesh surfaces is accumulated to obtain the number of mesh nodes on the coupling interface. It can be understood that in the case where there are multiple mesh surfaces in the coupling interface, the number of mesh points on each mesh surface is counted respectively. , and then accumulate to get the number of mesh nodes on the coupling interface .

[0061] Specifically, in the fluid solver, all grid blocks of the wing flow field grid are traversed, and the wing coupling interface data are calculated according to the coupling interface boundary conditions set when the grid is generated. is 1, and the number of grid nodes on the coupling interface is calculated according to the boundary node number of the grid block where the coupling interface is located. The wing structure grid unit type of the structural solver is quadrilateral shell unit, all nodes of the wing finite element model are selected as coupling interface points, and the coupling interface quantity data is set. Equal to 1, calculate the number of nodes on the coupling interface .

[0062] Step S13: performing numerical calculation on the physical quantity to be communicated in the mesh nodes of the coupling interface to obtain the value of the physical quantity to be communicated of each mesh node.

[0063] In this embodiment, the multidisciplinary coupling analysis problem is calculated according to each discipline solver to obtain physical field distribution data; the physical field distribution data and the mesh data of the coupling interface are calculated by a preset data processing technology to calculate the physical quantity value to be communicated on the mesh node of the coupling interface. It can be understood that according to the physical field distribution data calculated by the discipline solver in the multidisciplinary coupling analysis, the physical quantity value communicated on the mesh node of the coupling interface is calculated by data processing techniques such as block integration and variable conversion. For example, in fluid-solid coupling analysis, the physical field data calculated by the fluid includes physical quantity information such as pressure, density, and velocity. Based on the mesh data of the coupling interface, the normal vector and area of ​​the mesh unit are calculated, and the pressure load size and direction on the mesh unit are obtained by combining the pressure data, and then the force on the unit is distributed to the mesh node by linear weighting. It should be noted that in the multidisciplinary coupling analysis, the physical field distribution data calculated by different discipline solvers have various forms and complex structures. The physical field data calculated by the fluid contains a variety of original physical quantity information such as pressure, density, and velocity. These original physical quantity data cannot be directly used for the calculation of the communication physical quantity value on the mesh node of the coupling interface. Because the communication physical quantities (such as node forces, node displacements, etc.) on the coupling interface mesh nodes need to be determined according to the specific coupling problems and physical mechanisms, the format and units of the original data may not meet the requirements. Data interaction between different disciplines requires that the data have consistent physical meanings and appropriate forms on the coupling interface mesh nodes in order to accurately transmit information and perform subsequent coupling calculations. Therefore, for physical quantities such as pressure obtained by fluid calculation, the mesh unit normal vector and area are calculated based on the coupling interface mesh data, and the pressure load size and direction on the mesh unit can be obtained by block integration combined with the pressure data. This is because the pressure distribution on the unit surface is not uniform. Through block integration, the pressure effect on the unit surface can be integrated to obtain the equivalent pressure load on the entire unit. For example, in a two-dimensional fluid-structure coupling interface mesh unit, the pressure may vary on the unit boundary. By dividing the unit into multiple small blocks (blocks), the pressure on each small block is multiplied by the area corresponding to the small block (integral operation), and the total pressure resultant force on the unit can be obtained. Variable conversion is used to convert the units or representations of different physical quantities into a form suitable for coupling calculations. For example, in some cases, the unit of a physical quantity obtained from one discipline may be one of the units in the International System of Units, while another discipline requires another input unit. The units of physical quantities can be converted through variable conversion. At the same time, variable conversion can also convert physical quantities from one representation to another. For example, from a volume-based physical quantity to a mass-based physical quantity, etc., to meet the needs of communicating physical quantity calculations on the coupled interface mesh nodes.For example, in thermal-structural coupling, the heat transfer rate obtained from the thermal discipline may be in joules per second (J / s), while in the structural discipline, it may need to be converted to watts (W) (because 1W=1J / s). It may also need to be converted into physical quantities related to the temperature change of the structure in order to accurately communicate and calculate physical quantities such as temperature-stress on the coupling interface mesh nodes.

[0064] Specifically, the fluid solver uses the unsteady flow method to advance the flow field by one time step in physical time. Traverse all unit surfaces on the coupling interface, perform area numerical integration of the pressure on the unit surface, and obtain the aerodynamic resultant force on the unit surface. The aerodynamic force of the unit surface is distributed to each node of the unit surface using numerical processing methods such as the linear weighted method. The aerodynamic force on the same node is accumulated to obtain the aerodynamic force data of all nodes on the unit surface. The structural solver uses the aerodynamic force of the surface nodes of the wing finite element model obtained after interpolation as the external load, and solves the structural control dynamics equation based on the modal method to advance the structural motion by one time step, and obtain the displacement data of the surface nodes of the wing finite element model.

[0065] Step S14: convert the current integer data form of the coupling interface number, the grid node number of a single coupling interface, and the physical field dimension into the target coupling interface number, target grid node number, and target physical field dimension in real data form, and determine the physical quantity encoding order.

[0066] In this embodiment, in order to realize the one-time transmission of communication data, all communication data on the coupling interface grid nodes are encoded into a one-dimensional real array, which is mainly realized through three processes. First, type conversion is performed on part of the communication data, and the integer communication data in the communication data is converted into real communication data by adding real data zero or forcing type conversion. The text description string data is converted into an integer data sequence through mapping, and then converted into real communication data. The real communication data is not transformed and is directly assigned and used; secondly, the physical quantity values ​​to be communicated on the coupling interface grid nodes are sorted and encoded, that is, according to the physical quantity data type, the arrangement order of the physical quantity values ​​to be communicated in the one-dimensional real array is determined, that is, the physical quantity encoding order, wherein the physical quantity encoding order in the fluid solver is established as: node coordinates, node aerodynamic force; the physical quantity encoding order in the structure solver is established as: node coordinates, node displacement.

[0067] Specifically, the fluid solver and the structural solver encode the coupling data into a one-dimensional real array in their respective solvers. First, the fluid / structural solver converts the number of mesh nodes on a single coupling interface into By adding real data zero or forcing type conversion, the integer data is converted to real data. The node coordinates, node aerodynamic forces, and node displacement data are real data, so no conversion is required. is equal to 1, so there is no need to establish a coupling interface coding order.

[0068] Step S15: Encode the target physical field dimension and the target coupling interface quantity to obtain a first real array, and encode the target grid node quantity, and fill the encoded information into the first real array to obtain a second real array.

[0069] In this embodiment, the overall information data of the coupled interface grid nodes after conversion (for example, the converted physical field dimensions) is first converted into , the number of coupled interface mesh surfaces after conversion ) is encoded into a one-dimensional real array, that is, the first real array, and then the number of target grid nodes on the coupling interface grid surface is encoded one by one, and the encoded information is filled into the first real array to obtain the second real array.

[0070] Specifically, encoding all coupling data into a one-dimensional array includes: the fluid / structure solver first converts the physical field dimension of the coupling interface mesh node into , the number of coupling interfaces in the solver The overall information data is encoded into a one-dimensional array, that is, the first real array; then the wing coupling interface data is encoded, first the number of nodes on the wing surface Continue encoding into the one-dimensional array to obtain the second real array.

[0071] like Figure 2 As shown, a fluid solver coupling data encoding format diagram is disclosed. First, the physical field dimension , Number of coupling interfaces , Number of nodes , followed by the node coordinates ( )to( ), and finally the node aerodynamic force ( )to( ), similarly, the structural solver coupling data encoding format is as follows Figure 3 As shown, first, the physical field dimension , Number of coupling interfaces , Number of nodes , followed by the node coordinates ( )to( ), and finally the node displacement ( )to( ).

[0072] Step S16: cyclically encode the physical quantity to be communicated according to the physical quantity encoding order, and fill the encoding result into the second real type array to obtain a target real type array.

[0073] In this embodiment, the communication data (physical quantity to be communicated) are cyclically encoded one by one according to the physical quantity encoding order, and the coupling interface cyclically encodes the communication data one by one into the second real type array to obtain a target real type array. The target real type array is a one-dimensional coupling data array.

[0074] Specifically, after obtaining the second real array, the fluid solver encodes the node coordinates and the node aerodynamic force of the wing fluid grid surface into the second real array according to the corresponding physical quantity coding order; the structural solver encodes the node coordinates and the node displacement of the wing finite element model into the second real array, until the coupling interfaces are encoded one by one according to the physical quantity coding order, and the target real array is obtained.

[0075] It can be seen that the present application discloses a method for encoding communication data of multidisciplinary coupling surface grid points, including: determining the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem; forming a corresponding discipline analysis calculation grid according to the multidisciplinary coupling analysis problem, extracting the coupling interface grid from the discipline analysis calculation grid, and calculating the number of coupling interfaces and the number of grid nodes of the coupling interface; performing numerical calculations on the physical quantity to be communicated in the grid nodes of the coupling interface to obtain the value of the physical quantity to be communicated of each of the grid nodes; converting the number of coupling interfaces, single and double coupling interfaces in the current integer data form into the physical quantity data type, physical field dimension, and physical quantity to be communicated; and finally, calculating the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data in the multidisciplinary coupling analysis problem ... The number of grid nodes of a coupling interface and the physical field dimension are converted into the target coupling interface number, target grid node number and target physical field dimension in the form of real data, and the physical quantity coding order is determined; the target physical field dimension and the target coupling interface number are encoded to obtain a first real array, and the target grid node number is encoded, and the encoded information is filled into the first real array to obtain a second real array; the physical quantity to be communicated is cyclically encoded according to the physical quantity coding order, and the encoding result is filled into the second real array to obtain the target real array. It can be seen that by uniformly converting the key information such as the number of coupling interfaces in the form of integer data, the number of grid nodes of a single coupling interface, and the physical field dimension into the form of real data, the data type is more consistent in the subsequent coding process, avoiding additional processing steps and potential errors that may be caused by inconsistent data types, thereby improving the overall efficiency of data processing. Encoding according to the physical quantity coding order reduces confusion and uncertainty in the coding process, improves the accuracy and speed of coding, and thus improves the data processing efficiency of the entire multidisciplinary coupling analysis. After a series of encoding processes, the target real array is obtained. This method of integrating multiple types of physical quantities to be communicated into a unified format (target real array) is conducive to data communication between different discipline solvers in multidisciplinary coupling analysis. It enables one data transfer to carry all relevant coupling data, avoiding the communication overhead and time delay that may be caused by multiple transfers of different types of data, thereby enhancing the effect of data communication and improving the overall operating efficiency of multidisciplinary coupling analysis.

[0076] Reference Figure 4 As shown, an embodiment of the present invention discloses a multidisciplinary coupling surface grid point communication data decoding method, comprising:

[0077] Step S21: Obtain a target real array encoded using the multidisciplinary coupling surface grid point communication data encoding method.

[0078] In this embodiment, according to the communication mapping relationship between the communication domains of the subject solvers established by the multidisciplinary coupling analysis, the fluid solver transmits all the encoded coupling data to the structural solver at one time, and the structural solver transmits all the encoded coupling data to the fluid solver at one time. Therefore, each subject solver can obtain the target real data sent by other subject calculators through the communication mapping relationship. For example: during fluid-solid coupling analysis, the fluid solver encodes the node coordinates and node aerodynamic forces, and then transmits the coupling data to the structural solver through the solver communication mapping relationship. The structural solver receives the coupling data and calculates the aerodynamic interpolation matrix.

[0079] Step S22: extracting the target physical field dimension, the target coupling interface number, and the target grid node number from the target real array, and converting the target physical field dimension, the target coupling interface number, and the target grid node number into data format to obtain the physical field dimension, the coupling interface number, and the grid node number in the form of integer data.

[0080] Step S23: Obtain corresponding communication data from the encoded physical quantity according to the physical quantity encoding order, the physical field dimension, the number of coupling interfaces and the number of grid nodes.

[0081] In this embodiment, after receiving the one-dimensional real coupling data, the subject solver performs a decoding operation on the coupling data and decodes the one-dimensional coupling data into a single physical quantity array. The decoding operation can also be implemented through three processes. First, the overall information data of the coupling interface grid node is obtained from the encoded data, for example: the physical field dimension Number of coupling interfaces . Then loop times, each time obtaining a type of communication data on the coupling interface grid point from the encoded data. According to the arrangement order of the physical quantity data established in the encoding process, the communication data obtained each time is restored to the physical quantity corresponding to the same physical meaning. Finally, when all the coupling data are decoded, the communication data is converted to the original type, and the integer data is converted from real data to integer data by discarding the decimal part or forcing type conversion. The text string of the real sequence is first converted to integer data through real data, and then matched to the character sequence according to the integer data according to the same mapping conversion.

[0082] Specifically, the fluid / structure solver obtains the overall information data of the coupling interface grid nodes from the received one-dimensional coupling data array: physical quantity dimension Number of coupling interfaces ; Then, obtain the wing surface physical quantity data from the coupling data, first obtain the number of coupling interface nodes , and convert it from real data to integer data by discarding the decimal part or forcing type conversion; the fluid solver obtains the quantities of node coordinates, node aerodynamic force, and The structural solver obtains the quantities of And assign each type of coupling data obtained to the solver parameters with the same physical meaning.

[0083] In this embodiment, after obtaining the corresponding communication data from the encoded physical quantity according to the physical quantity encoding order, the physical field dimension, the number of coupling interfaces and the number of grid nodes, it also includes: using the current communication data decoded by any subject solver as control point information, and using the control point information and the initial interpolation function to construct a corresponding target interpolation function; using the target interpolation function to perform interpolation calculations on other communication data except the current communication data to obtain the interpolated communication data. It can be understood that the fluid solver uses the fluid node coordinates contained in itself and the structural node coordinates obtained by decoding as control points, and establishes a displacement interpolation matrix based on the Radial Basis Function (RBF) method. The structural node displacement obtained by decoding is taken as a known vector, and the matrix-vector product is performed with the displacement interpolation matrix to obtain the node displacement of the fluid grid wing surface. Similarly, the structural solver uses the decoded fluid node coordinates and the structural node coordinates contained in it as control points, and establishes the aerodynamic interpolation matrix based on the principle of virtual work. The aerodynamic force of the fluid node obtained by decoding is taken as a known vector, and the matrix-vector product is performed with the aerodynamic force interpolation matrix to obtain the aerodynamic force of the structured grid wing surface node.

[0084] After obtaining the node displacements of the fluid grid wing surface, the grid deformation methods such as RBF-TFI (RBF-TransfiniteInterpolation, radial basis function-transfinite interpolation) are used to calculate the internal node coordinates of the grid. The steps of the RBF-TFI method are as follows: First, the vertices of the grid blocks on the coupling interface are used as control points, and the vertices of the internal grid blocks are used as interpolation base points. The RBF method is used to calculate the vertex displacements of the internal grid blocks according to the displacements of the grid blocks on the coupling interface. Secondly, the vertices of all grid blocks are used as control points, and RBF is used to accurately solve the node displacements of the edges of each grid block. Then, based on the node displacements of the grid block edges, the two-dimensional TFI method based on arc length coordinates is used to quickly calculate the node displacements of each grid block surface. Finally, based on the node displacements of the grid surface of the grid block, the three-dimensional TFI method based on arc length coordinates is used to quickly calculate the displacements of the internal points of the grid block. After the displacements of all nodes are calculated, they are added to the current coordinates of the nodes to obtain the new coordinates of the nodes, and the geometric information of the grid (such as: node coordinates, center of mass coordinates, volume, area, area vector) and Jacobi transformation coefficient matrix are recalculated.

[0085] In this way, after decoding and obtaining the communication data, the current communication data decoded by any subject solver can be used as the control point information, and the corresponding target interpolation function can be constructed in combination with the initial interpolation function. For example, the fluid solver uses the structural node coordinates obtained by decoding and its own fluid node coordinates to construct a displacement interpolation matrix, and the structural solver uses the fluid node aerodynamic force obtained by decoding and its own structural node coordinates to construct an aerodynamic force interpolation matrix. This method of constructing an interpolation function based on decoded data enables different subject solvers to effectively use each other's decoded data for further calculation and analysis, thereby improving the efficiency of data utilization. The interpolated communication data can be obtained by interpolating the communication data other than the current communication data using the constructed target interpolation function. For example, the displacement of the node on the surface of the fluid grid wing is calculated by the displacement interpolation matrix, and the aerodynamic force of the node on the surface of the structural grid wing is calculated by the aerodynamic force interpolation matrix. Interpolation calculation can supplement and improve the missing or refined data based on the relationship between the data on the basis of the existing data, thereby improving the accuracy of multidisciplinary coupling analysis in simulating the actual physical process and more accurately reflecting the physical phenomena and system behavior.

[0086] When the number of coupling iterations does not exceed the set maximum number of coupling iterations, the physical time accumulates a coupling time step, and the number of coupling iterations increases by 1. Jump to execute the fluid and structure solver time advancement until the step required by the maximum number of coupling analysis steps is reached. Figure 5The wing flutter analysis process based on the encoding and decoding method is demonstrated, in which a loosely coupled strategy format with first-order coupling accuracy is adopted. Specifically, the initialization phase: After the start, first parse the coupling configuration file and generate a global communication domain. Perform fluid solver initialization and structural solver initialization respectively. Generate the fluid / structural solver communication domain, and extract the coupling interface data in the two solvers respectively. Flutter analysis phase: Import the subject solver module, and the global communication domain extracts the physical coordinates of all processes. The fluid solver and the structural solver perform time advancement operations respectively. At the same time, the fluid solver establishes a communication domain process mapping relationship based on the physical coordinates and extracts the fluid coupling data on the wing surface. Determine whether the maximum step size is met. If not, perform the following fluid coupling data encoding operations: 1. Perform fluid coupling data communication and send the encoded fluid coupling data. 2. The structural solver receives the fluid coupling data and decodes the fluid coupling data. 3. Interpolate the node aerodynamic force. 4. The structural solver extracts the wing surface structural coupling data and encodes the structural coupling data. 5. Perform structural coupling data communication and send out the encoded structural coupling data. 6. The fluid solver receives the structural coupling data and decodes the structural coupling data. 7. Interpolate the node displacement. 8. The fluid solver deforms the fluid mesh and recalculates the geometric data. Iteration phase: Physical time accumulates a time step. The number of coupling alternative steps is increased by 1, and then the step of judging whether the maximum alternative step is met continues to iterate until the maximum alternative step condition is met, and the process ends.

[0087] It can be seen that by extracting and converting key information (such as target physical field dimension, target coupling interface number, target grid node number, etc.) from the target real array in sequence into integer data form, the basic parameter settings of each physical quantity before encoding can be accurately restored, laying the foundation for the subsequent accurate decoding of physical quantity data. For example, clarifying the physical field dimension can determine the number of components of each physical quantity, and then accurately restore the original form of physical quantities such as node coordinates and node forces, ensuring that the data can restore its original physical meaning and format after decoding, and ensuring the accuracy of the data. In addition, in the decoding and restoration process, the encoded physical quantities of each physical quantity data type are decoded and restored according to the physical quantity encoding order established in the encoding process, so that the decoding operation can be carried out strictly in accordance with the rules during encoding, and the one-dimensional coupling data is accurately decoded into a single physical quantity array. This effectively avoids confusion and errors in the data restoration process, ensures that the communication data can be restored in the correct form and corresponding relationship, and improves the accuracy of data processing.

[0088] Reference Figure 6 As shown, the present invention also discloses a multidisciplinary coupling surface grid point communication data encoding device, including:

[0089] An information determination module 11 is used to determine the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem;

[0090] The data calculation module 12 is used to form a corresponding subject analysis calculation grid according to the multi-disciplinary coupling analysis problem, extract a coupling interface grid from the subject analysis calculation grid, and calculate the number of coupling interfaces and the number of grid nodes of the coupling interfaces;

[0091] A numerical calculation module 13, used for performing numerical calculation on the physical quantity to be communicated in the grid nodes of the coupling interface to obtain the value of the physical quantity to be communicated of each grid node;

[0092] A form conversion module 14 is used to convert the number of coupling interfaces, the number of grid nodes of a single coupling interface, and the physical field dimension in the current integer data form into the target number of coupling interfaces, the target number of grid nodes, and the target physical field dimension in the real data form, and determine the physical quantity coding order;

[0093] A first encoding module 15 is used to encode the target physical field dimension and the target coupling interface quantity to obtain a first real array, and encode the target grid node quantity, and fill the encoded information into the first real array to obtain a second real array;

[0094] The second encoding module 16 is used to perform cyclic encoding on the physical quantity to be communicated according to the physical quantity encoding sequence, and fill the encoding result into the second real type array to obtain a target real type array.

[0095] It can be seen that the present application discloses determining the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem; forming a corresponding discipline analysis calculation grid according to the multidisciplinary coupling analysis problem to extract the coupling interface grid from the discipline analysis calculation grid, and calculating the number of coupling interfaces and the number of grid nodes of the coupling interface; performing numerical calculations on the physical quantities to be communicated in the grid nodes of the coupling interface to obtain the physical quantity values ​​to be communicated of each of the grid nodes; converting the current number of coupling interfaces in integer data form, the number of grid nodes of a single coupling interface, and the physical field dimension into the target number of coupling interfaces, the target number of grid nodes, and the target physical field dimension in real data form, respectively, and determining the physical quantity encoding order; encoding the target physical field dimension and the target number of coupling interfaces to obtain a first real array, encoding the target number of grid nodes, and filling the encoded information into the first real array to obtain a second real array; cyclically encoding the physical quantities to be communicated according to the physical quantity encoding order, and filling the encoding results into the second real array to obtain a target real array. It can be seen that by converting the key information such as the number of coupling interfaces in the form of integer data, the number of grid nodes of a single coupling interface, and the physical field dimension into real data, the data type is more consistent in the subsequent encoding process, avoiding additional processing steps and potential errors caused by inconsistent data types, thereby improving the overall efficiency of data processing. Encoding is performed in the order of physical quantity coding to reduce confusion and uncertainty in the coding process, improve the accuracy and speed of coding, and thus improve the data processing efficiency of the entire multidisciplinary coupling analysis. After a series of encoding processes, the target real array is obtained. This method of integrating multiple types of physical quantities to be communicated into a unified format (target real array) is conducive to data communication between different discipline solvers in multidisciplinary coupling analysis. It enables one data transfer to carry all relevant coupling data, avoiding the communication overhead and time delay that may be caused by multiple transfers of different types of data, thereby enhancing the effect of data communication and improving the overall operating efficiency of multidisciplinary coupling analysis.

[0096] Furthermore, the present application also discloses an electronic device. Figure 7 This is a structural diagram of an electronic device 20 according to an exemplary embodiment. The content in the diagram cannot be considered as any limitation on the scope of use of the present application.

[0097] Figure 7The present invention provides a schematic diagram of the structure of an electronic device 20 provided in an embodiment of the present application. The electronic device 20 may specifically include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 is used to store a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the method disclosed in any of the above embodiments. In addition, the electronic device 20 in this embodiment may specifically be an electronic computer.

[0098] In this embodiment, the power supply 23 is used to provide working voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and the external device, and the communication protocol it follows is any communication protocol that can be applied to the technical solution of the present application, and is not specifically limited here; the input and output interface 25 is used to obtain external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs and is not specifically limited here.

[0099] Among them, the processor 21 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 21 can be implemented in at least one hardware form of DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array). The processor 21 may also include a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also known as a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in the standby state. In some embodiments, the processor 21 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 21 may also include an AI (Artificial Intelligence) processor, which is used to process computing operations related to machine learning.

[0100] In addition, the memory 22, as a carrier for storing resources, can be a read-only memory, a random access memory, a disk or an optical disk, etc. The resources stored thereon can include an operating system 221, a computer program 222, etc., and the storage method can be temporary storage or permanent storage.

[0101] Among them, the operating system 221 is used to manage and control the hardware devices and computer programs 222 on the electronic device 20, so as to realize the operation and processing of the massive data 223 in the memory 22 by the processor 21, which can be Windows Server, Netware, Unix, Linux, etc. In addition to including a computer program that can be used to complete the method performed by the electronic device 20 disclosed in any of the aforementioned embodiments, the computer program 222 can further include a computer program that can be used to complete other specific tasks. In addition to data transmitted from an external device received by the electronic device, the data 223 can also include data collected by its own input and output interface 25, etc.

[0102] Furthermore, the present application also discloses a computer-readable storage medium for storing a computer program; wherein the computer program, when executed by a processor, implements the aforementioned disclosed method. The specific steps of the method can refer to the corresponding contents disclosed in the aforementioned embodiments, and will not be repeated here.

[0103] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0104] Professionals may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented with electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel may use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application. The steps of the method or algorithm described in conjunction with the embodiments disclosed herein can be implemented directly with hardware, a software module executed by a processor, or a combination of the two. The software module can be placed in a random access memory RAM (Random Access Memory), memory, read-only memory ROM (Read Only Memory), electrically programmable EPROM (Electrically Programmable Read Only Memory), electrically erasable programmable EEPROM (ElectricErasable Programmable Read Only Memory), register, hard disk, removable disk, CD-ROM (CompactDisc-Read Only Memory), or any other form of storage medium known in the technical field.

[0105] Finally, it should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0106] The scheme provided by the present invention is introduced in detail above. Specific examples are used in this article to illustrate the principle and implementation mode of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation mode and application scope. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A method for encoding communication data of multidisciplinary coupling surface grid points, characterized in that: include: Determine the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem; According to the multidisciplinary coupling analysis problem, a corresponding discipline analysis calculation grid is formed to extract a coupling interface grid from the discipline analysis calculation grid, and the number of coupling interfaces and the number of grid nodes of the coupling interfaces are calculated; Performing numerical calculation on the physical quantity to be communicated in the grid nodes of the coupling interface to obtain the value of the physical quantity to be communicated of each of the grid nodes; Convert the current number of coupling interfaces, the number of grid nodes of a single coupling interface, and the physical field dimension in the form of integer data into the target number of coupling interfaces, the target number of grid nodes, and the target physical field dimension in the form of real data, respectively, and determine the physical quantity encoding order; The target physical field dimension and the target coupling interface quantity are encoded to obtain a first real array, and the target grid node quantity is encoded, and the encoded information is filled into the first real array to obtain a second real array; The physical quantity to be communicated is cyclically encoded according to the physical quantity encoding order, and the encoding result is filled into the second real type array to obtain a target real type array.

2. The multidisciplinary coupling surface grid point communication data encoding method according to claim 1 is characterized in that: Before determining the physical quantity data type, physical field dimension, and physical quantity of the coupling data to be communicated in the multidisciplinary coupling analysis problem, the method further includes: Generate a calculation example folder for a multidisciplinary coupling analysis problem; wherein the calculation example folder is used to store the folders of each discipline solver and the corresponding coupling framework configuration file; When the multidisciplinary coupling analysis problem is obtained, the example folder is called to simulate the multidisciplinary coupling analysis problem to obtain a coupling interface involved in the coupling analysis.

3. The multidisciplinary coupling surface grid point communication data encoding method according to claim 1 is characterized in that: The determination of the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem includes: According to the multidisciplinary coupling analysis problem and the mechanism of iterative exchange between different disciplines, the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated participating in the coupling communication in the multidisciplinary coupling analysis problem are determined.

4. The multidisciplinary coupling surface grid point communication data encoding method according to claim 1 is characterized in that: Calculate the number of mesh nodes for the coupling interface, including: If there are multiple mesh surfaces in the coupling interface, the number of mesh points of each mesh surface is counted respectively, and the number of mesh points of all the mesh surfaces is cumulatively calculated to obtain the number of mesh nodes of the coupling interface.

5. The multidisciplinary coupling surface grid point communication data encoding method according to claim 1 is characterized in that: The performing numerical calculation on the physical quantity to be communicated in the grid nodes of the coupling interface to obtain the value of the physical quantity to be communicated of each grid node includes: Calculating the multidisciplinary coupling analysis problem according to the solvers of each discipline to obtain physical field distribution data; The physical field distribution data and the mesh data of the coupling interface are processed by a preset data processing technology to calculate the physical quantity values ​​to be communicated on the mesh nodes of the coupling interface.

6. The multidisciplinary coupling surface grid point communication data encoding method according to claim 1, characterized in that: The target real type array is a one-dimensional coupled data array.

7. A method for decoding communication data of multidisciplinary coupling surface grid points, characterized in that: include: Obtain a target real array encoded by the multidisciplinary coupling surface grid point communication data encoding method according to any one of claims 1 to 6; Extracting target physical field dimensions, target coupling interface quantity, and target grid node quantity from the target real array, and converting the target physical field dimensions, target coupling interface quantity, and target grid node quantity into data format to obtain the physical field dimensions, coupling interface quantity, and grid node quantity in the form of integer data; Corresponding communication data is obtained from the encoded physical quantity according to the physical quantity encoding order, the physical field dimension, the number of coupling interfaces and the number of grid nodes.

8. The multidisciplinary coupling surface grid point communication data decoding method according to claim 7 is characterized in that: After obtaining corresponding communication data from the encoded physical quantity according to the physical quantity encoding order, the physical field dimension, the number of coupling interfaces and the number of grid nodes, the method further includes: Using the current communication data decoded by any subject solver as control point information, and using the control point information and the initial interpolation function to construct a corresponding target interpolation function; The target interpolation function is used to perform interpolation calculation on other communication data except the current communication data to obtain interpolated communication data.

9. A multidisciplinary coupling surface grid point communication data encoding device, characterized in that: include: An information determination module, used to determine the physical quantity data type, physical field dimension, and physical quantity to be communicated of the coupling data to be communicated in the multidisciplinary coupling analysis problem; A data calculation module, used to form a corresponding subject analysis calculation grid according to the multidisciplinary coupling analysis problem, to extract a coupling interface grid from the subject analysis calculation grid, and to calculate the number of coupling interfaces and the number of grid nodes of the coupling interfaces; A numerical calculation module, used for performing numerical calculation on the physical quantity to be communicated in the grid nodes of the coupling interface to obtain the value of the physical quantity to be communicated of each of the grid nodes; A form conversion module, used to convert the number of coupling interfaces, the number of grid nodes of a single coupling interface, and the physical field dimension in the current integer data form into the target number of coupling interfaces, the target number of grid nodes, and the target physical field dimension in the real data form, and determine the physical quantity coding order; A first encoding module is used to encode the target physical field dimension and the target coupling interface quantity to obtain a first real array, and encode the target grid node quantity, and fill the encoded information into the first real array to obtain a second real array; The second encoding module is used to perform cyclic encoding on the physical quantity to be communicated according to the physical quantity encoding order, and fill the encoding result into the second real type array to obtain a target real type array.

10. An electronic device, characterized in that: include: Memory, used to store computer programs; A processor, configured to execute the computer program to implement the steps of the method according to any one of claims 1 to 8.

Images