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

Abstract

The invention discloses a multidisciplinary coupling surface extensible communication data coding method, decoding method, device and equipment, and relates to the field of multidisciplinary coupling numerical simulation, and the method comprises the steps: coding a plurality of to-be-communicated coupling data values of any physical quantity type and type of a plurality of coupling interfaces and coupling data basic information into a coupling data array, data communication on a plurality of coupling interfaces can be completed through one-time calling of the message passing interface, the calling frequency of the message passing interface is greatly reduced, coupling data communication can be carried out on physical data in subject solvers of different numerical simulation methods according to the same encoding and decoding format, coupling data types can be randomly expanded, and the coupling data communication efficiency is improved. The development difficulty of multidisciplinary coupling analysis software is reduced, and the efficiency of multidisciplinary coupling analysis 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 scalable communication data encoding method, decoding method, device and equipment. Background Art

[0002] At present, the partition coupling solution strategy for each discipline can use mature single-discipline solvers, thus greatly reducing the difficulty of solver development and being applicable to extremely complex application scenarios, becoming the main way to achieve multidisciplinary coupling simulation. In the partition coupling solution strategy, different discipline solvers run independently, and the iterative update of different physical field data is achieved through the exchange of interface information. The coupling of multiple components through multiple interfaces to exchange coupling information has become a common phenomenon in multidisciplinary coupling analysis scenarios. For example, in the aerodynamic / structural coupling analysis of aircraft main wings, ailerons, tail wings and other components, the fluid solver and the structural solver need to exchange data from multiple wing coupling interfaces.

[0003] In the analysis of multiple coupling interfaces, there are many types of coupling data physical quantities and data types. At the same time, the storage locations of coupling data of different physical quantity types on the discrete grid are inconsistent, resulting in the need to use different types of points as interpolation base points when interpolating coupling data. Therefore, when communicating coupling data on multiple coupling interfaces based on the message passing interface, it is necessary not only to communicate the coupling data itself, but also to communicate basic information such as the type, type and distribution location of the coupling data in the discrete grid to other discipline solvers. However, the message passing interface method can only send and receive a single type of data at a time. When there are many types and 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 scalable communication data encoding method, decoding method, device and equipment, which can reduce the call of the message passing interface and improve the efficiency of multidisciplinary coupling simulation analysis. The specific scheme is as follows:

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

[0006] Determine the physical quantity types and physical quantity quantities of all coupling data involved in coupling communication in the multidisciplinary coupling simulation analysis problem, calculate the number of coupling interfaces involved in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface, and determine the physical field dimension of the coupling data;

[0007] Assigning a corresponding type association symbol to each physical quantity of each physical quantity type, assigning a corresponding storage location symbol to each physical quantity stored in the grid unit surface of the coupling interface, assigning a corresponding data type symbol to each physical quantity of scalar type and / or vector type data type, establishing a coding order for each coupling interface, establishing a coupling data coding order for each physical quantity type on the current coupling interface, and establishing a grid node coding order in the current coupling interface for the location information of the grid node where each physical quantity is located on the current coupling interface;

[0008] Encode the storage location symbol and the data type symbol to obtain a first encoding result, traverse and encode all grid node information of a single coupling interface in the grid node encoding order to obtain a corresponding second encoding result, and construct a third encoding result of coupling data encoding characterizing a single physical quantity type on a single coupling interface according to the order of the first encoding result and the second encoding result;

[0009] Encoding the physical quantity quantity into a coupling data sequence, and filling the type association symbol into the encoded coupling data sequence, and then encoding the third encoding result into the encoded coupling data sequence according to the coupling data encoding order, so as to obtain a fourth encoding result after coupling data encoding representing all physical quantity types on a single coupling interface; wherein the coupling data sequence is constructed based on the values ​​of each physical quantity;

[0010] The basic information of the coupling data including the physical field dimension, the number of coupling interfaces, and the coding order of the coupling interfaces is encoded to obtain a basic information coding result, and the fourth coding result is filled after the basic information coding result to obtain a target coding result.

[0011] Optionally, the calculating the number of coupling interfaces of the coupling interfaces involved in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface includes:

[0012] According to the grid information of the subject solver for solving the multidisciplinary coupling simulation analysis problem, all boundaries of the corresponding current grid are traversed, and the number of coupling interfaces of the coupling interfaces participating in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface are calculated according to the coupling interface boundary conditions set when the current grid is generated.

[0013] Optionally, the step of assigning a corresponding type association symbol to each physical quantity of each physical quantity type includes:

[0014] The target associated symbols corresponding to the physical quantity types are screened from the preset physical quantity type association table to obtain a type associated symbol sequence constructed with the target associated symbols.

[0015] Optionally, the step of establishing a coupling data coding sequence of each physical quantity type on the current coupling interface includes:

[0016] According to the arrangement order of the physical quantity type-association symbols in the preset physical quantity type association table, the target arrangement order of each of the physical quantity types on the current coupling interface is obtained, so as to use the target arrangement order as the coupling data encoding order of each of the physical quantity types on the current coupling interface.

[0017] Optionally, the assigning a corresponding storage location symbol to each of the physical quantities stored in the grid unit surface of the coupling interface, and assigning a corresponding data type symbol to each of the physical quantities of scalar and / or vector data types, comprises:

[0018] If the physical quantity is stored on a grid node in a grid unit surface of the coupling interface, the storage location symbol assigned to the physical quantity is 1, otherwise, the storage location symbol assigned to the physical quantity is 0;

[0019] If the type of the physical quantity is a scalar type, the data type symbol assigned to the physical quantity is 0; if the type of the physical quantity is a vector type, the data type symbol assigned to the physical quantity is 1.

[0020] Optionally, the establishing of the coding sequence of each coupling interface includes:

[0021] A coding sequence of each coupling interface is established according to the structural topological grid type and / or the non-structural topological grid type, and a corresponding coupling interface association symbol is allocated to each coupling interface.

[0022] Optionally, after obtaining the target encoding result, the method further includes:

[0023] According to the communication mapping relationship between the communication domains of each subject solver, the target encoding result is communicated to other subject solvers.

[0024] In a second aspect, the present application discloses a multidisciplinary coupling surface scalable communication data decoding method, comprising:

[0025] Obtaining a target coding result encoded by a subject solver using the multi-disciplinary coupling surface scalable communication data coding method;

[0026] Initially decoding the target encoding result to sequentially obtain the physical field dimension of the coupling data, the number of coupling interfaces, and the encoding order of the coupling interfaces;

[0027] According to the coding sequence of the coupling interfaces, the fourth coding result of each coupling interface is decoded twice to obtain the number of physical quantities, and the step of decoding the fourth coding result of each coupling interface according to the coding sequence of the coupling interfaces is executed twice for the same number of times as the number of physical quantities to obtain a type association symbol sequence, and the third coding result corresponding to each type association symbol in the type association symbol sequence is determined according to a preset physical quantity type association table;

[0028] The third encoding result is decoded to obtain a storage location symbol and a data type symbol as well as information of each grid node on a single coupling interface to obtain decoded coupling data.

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

[0030] An information determination module is used to determine the physical quantity type and physical quantity quantity of all coupling data involved in coupling communication in a multidisciplinary coupling simulation analysis problem, and calculate the number of coupling interfaces of the coupling interfaces involved in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface, and determine the physical field dimension of the coupling data;

[0031] An information allocation module, used to allocate a corresponding type association symbol to each physical quantity of each physical quantity type, allocate a corresponding storage location symbol to each physical quantity stored in a grid unit surface of the coupling interface, allocate a corresponding data type symbol to each physical quantity of a scalar and / or vector data type, establish a coding sequence for each coupling interface, establish a coupling data coding sequence for each physical quantity type on a current coupling interface, and establish a grid node coding sequence in a current coupling interface based on the position information of the grid node where each physical quantity is located on the current coupling interface;

[0032] A first encoding module is used to encode the storage location symbol and the data type symbol to obtain a first encoding result, traverse and encode all grid node information of a single coupling interface in the grid node encoding order to obtain a corresponding second encoding result, and construct a third encoding result of coupling data encoding representing a single physical quantity type on a single coupling interface according to the order of the first encoding result and the second encoding result;

[0033] a second encoding module, for encoding the physical quantity quantity into a coupling data sequence, and filling the type association symbol into the encoded coupling data sequence, and then encoding the third encoding result into the encoded coupling data sequence according to the coupling data encoding order, so as to obtain a fourth encoding result after coupling data encoding representing all types of physical quantities on a single coupling interface; wherein the coupling data sequence is constructed based on the values ​​of each physical quantity;

[0034] The third encoding module is used to encode the basic information of the coupling data including the physical field dimension, the number of coupling interfaces, and the encoding order of the coupling interfaces to obtain a basic information encoding result, and fill the fourth encoding result after the basic information encoding result to obtain a target encoding result.

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

[0036] Memory, used to store computer programs;

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

[0038] It can be seen that the present application discloses a multidisciplinary coupling surface scalable communication data encoding method, including: determining the physical quantity type and physical quantity quantity of all coupling data involved in the coupling communication in the multidisciplinary coupling simulation analysis problem, and calculating the number of coupling interfaces of the coupling interfaces involved in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface, and determining the physical field dimension of the coupling data; assigning a corresponding type association symbol to each physical quantity of each physical quantity type, assigning a corresponding storage location symbol to each physical quantity stored in the grid unit surface of the coupling interface, assigning a corresponding data type symbol to each physical quantity of scalar and / or vector data types, establishing a coding order for each coupling interface, establishing a coupling data coding order for each physical quantity type on the current coupling interface, and establishing a grid node coding order in the current coupling interface for the position information of the grid node where each physical quantity is located on the current coupling interface; encoding the storage location symbol and the data type symbol to obtain a first The encoding result is obtained by traversing and encoding all grid node information of a single coupling interface in accordance with the grid node encoding order to obtain a corresponding second encoding result, and constructing a third encoding result of coupling data encoding representing a single physical quantity type on a single coupling interface according to the order of the first encoding result and the second encoding result; encoding the number of physical quantities into a coupling data sequence, and filling the type association symbol into the encoded coupling data sequence, and then encoding the third encoding result into the encoded coupling data sequence in accordance with the coupling data encoding order to obtain a fourth encoding result after coupling data encoding representing all physical quantity types on a single coupling interface; wherein the coupling data sequence is constructed based on the values ​​of each physical quantity; encoding the coupling data basic information including the physical field dimension, the number of coupling interfaces, and the encoding order of the coupling interfaces to obtain a basic information encoding result, and filling the fourth encoding result into the basic information encoding result to obtain a target encoding result. It can be seen that through the allocation process of type association symbols and storage location symbols, it is possible to process coupled data of various physical quantity types and different data types (scalar type, vector type), meet the needs of complex multidisciplinary coupling analysis scenarios such as fluid-solid coupling and fluid-solid thermal coupling, reflect the processing ability of various data types, and enhance flexibility. By establishing the coding order of the coupling interface, the coding order of the coupling data, and the coding order of the grid nodes, different grid structures (structural topological grids and non-structural topological grids) and the storage locations (nodes or face centers) of physical quantities in the grid unit surface of the coupling interface are further considered. The coding order is established for these differences, so that the method can adapt to different grid space discretization methods and physical quantity storage methods of various discipline solvers, and improve adaptability.Finally, based on the above information and through the preset coding method, the basic information of the physical quantities and coupling information of various physical quantity types of each coupling interface are uniformly encoded, and then the data communication on multiple coupling interfaces can be completed through a single message passing interface call. The number of calls to the message passing interface is greatly reduced, and the problem of increased time cost caused by multiple calls to the interface is reduced, especially when there are a large number of calculation processes, which effectively avoids the problem of communication time accounting for too large a proportion of the overall calculation time, thereby improving the efficiency of multidisciplinary coupling simulation analysis. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] 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.

[0040] Figure 1 A flow chart of a multi-disciplinary coupling surface scalable communication data encoding method disclosed in this application;

[0041] Figure 2 A fluid solver coupling data encoding process disclosed in this application;

[0042] Figure 3 A flow chart of a multi-disciplinary coupling surface scalable communication data decoding method disclosed in this application;

[0043] Figure 4 A fluid solver coupling data decoding process disclosed in this application;

[0044] Figure 5 This is a schematic diagram of the structure of a multi-disciplinary coupling surface expandable communication data encoding device disclosed in this application;

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

[0046] 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.

[0047] With the rapid improvement of computer hardware and software capabilities and the continuous development and improvement of disciplinary numerical simulation 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 to provide high-precision prediction data for product design. However, with the continuous improvement of the sophistication and requirements of modern product design, a large number of 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.

[0048] At present, the partition coupling solution strategy for each discipline can use mature single-discipline solvers, thus greatly reducing the difficulty of solver development and being applicable to extremely complex application scenarios, becoming the main way to achieve multidisciplinary coupling simulation. In the partition coupling solution strategy, different discipline solvers run independently, and the iterative update of different physical field data is achieved through the exchange of interface information. The coupling of multiple components through multiple interfaces to exchange coupling information has become a common phenomenon in multidisciplinary coupling analysis scenarios. For example, in the aerodynamic / structural coupling analysis of aircraft main wings, ailerons, tail wings and other components, the fluid solver and the structural solver need to exchange data from multiple wing coupling interfaces.

[0049] In the analysis of multiple coupling interfaces, there are many types of coupling data physical quantities and data types. At the same time, the storage locations of coupling data of different physical quantity types on the discrete grid are inconsistent, resulting in the need to use different types of points as interpolation base points when interpolating coupling data. Therefore, when communicating coupling data on multiple coupling interfaces based on the message passing interface, it is necessary not only to communicate the coupling data itself, but also to communicate basic information such as the type and type of coupling data and the distribution location in the discrete grid to other subject solvers. However, the message passing interface method can only send and receive a single type of data each time. When there are many types and 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 account for a large proportion of the overall calculation time, thereby reducing the efficiency of multi-disciplinary coupling simulation analysis. At the same time, it is also necessary to communicate the basic information of the coupling data on each coupling interface and the coupling data value to other subject solvers at the same time, so that the subject solver can interpolate the coupling data value according to the basic information of the coupling data, linear weighting and other data processing techniques.

[0050] To this end, the present invention provides a multidisciplinary coupling surface scalable communication data encoding and decoding scheme, which can reduce the call of the message passing interface and improve the efficiency of multidisciplinary coupling simulation analysis.

[0051] Reference Figure 1 As shown, the embodiment of the present invention discloses a multi-disciplinary coupling surface scalable communication data encoding method, including:

[0052] Step S11: Determine the physical quantity type and physical quantity quantity of all coupling data involved in coupling communication in the multidisciplinary coupling simulation analysis problem, calculate the number of coupling interfaces of the coupling interfaces involved in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface, and determine the physical field dimension of the coupling data.

[0053] In view of the start-up / structural coupling simulation analysis problem of the aircraft, an aerodynamic / structural coupling analysis case folder of the aircraft is generated, which contains folders such as the fluid solver, structural solver, coupling data interpolation module and configuration files of the coupling framework. The solver folder stores the mesh and configuration files of the model respectively, and the coupling interpolation folder saves the interpolation algorithm parameter setting file. For the fluid solver setting, first use the mesh generation software to generate the mesh of the aircraft's external flow field, set all wing surfaces as coupling interfaces, 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, etc.), mesh deformation method (RBF-TFI, (Radial Basis Function-Transfinite Interpolation, radial basis function-transfinite interpolation)) and solution method (discrete format, turbulence model, etc.) in the fluid solver configuration file. Then set the structural solver configuration file, use the structural pre-processing software to input the structural physical property parameters, set all wings as coupling interfaces, and use the modal analysis function to calculate and export the first four eigenfrequencies and modal coordinates. In the coupling interpolation algorithm configuration file, set the interpolation method RBF, interpolation radius and other parameters. In the coupling framework configuration file, design the coupling analysis type, time step, maximum number of iterations, solver name involved in the coupling analysis, number of cores and configuration file address and other parameters.

[0054] In this embodiment, the type and quantity of physical quantities are first determined, and the type and quantity of all data physical quantities involved in the coupling communication in the multidisciplinary coupling analysis are determined according to the type of multidisciplinary coupling analysis problem and the coupling iterative mechanism between disciplines. pFor example, in fluid-solid coupling analysis, the fluid needs to transfer the calculated force on the coupling interface to the solid so that the solid can perform static or dynamic analysis. The solid needs to transfer the calculated displacement after deformation to the fluid so that the fluid can update its shape and calculate the force under the new state. Therefore, the types and quantities of physical quantities involved in coupling communication in fluid-solid coupling include: physical field dimensions , the number of coupled surface mesh nodes , grid node coordinates , Node aerodynamics , node displacement .

[0055] Specifically, in the aerodynamic / structural coupling analysis of the aircraft wing, after the fluid solver solves a single coupling time step, the node aerodynamic force on the coupling interface is passed to the structural solver for dynamic analysis; after the structural solver solves a single coupling time step, the displacement on the coupling interface is passed to the fluid solver for fluid mesh deformation. Therefore, the types of physical quantities involved in the coupling communication of the fluid solver include: node coordinates, node aerodynamic force, and setting the number of physical quantities. Equal to 2; the types of physical quantities involved in the coupling communication of the structural solver include: node coordinates, node displacements, and setting the number of physical quantities Equals 2.

[0056] In this embodiment, according to the grid information of the subject solver for solving the multidisciplinary coupling simulation analysis problem, all boundaries of the corresponding current grid are traversed, and the number of coupling interfaces of the coupling interfaces participating in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface are calculated according to the coupling interface boundary conditions set when the current grid is generated. It can be understood that according to the grid information of the subject solver, all boundaries of the solver grid are traversed, and the number of coupling interfaces participating in the multidisciplinary coupling analysis is calculated according to the coupling interface boundary type or boundary name. and the number of mesh nodes on each coupling interface For example, in multidisciplinary analysis problems with boundary surface coupling, such as aerodynamic structure coupling analysis of aircraft, there are usually multiple wing-like coupling interfaces and a large number of grid nodes are distributed on each coupling interface. In the grid file, all boundaries are traversed, and the number of coupling interfaces is calculated according to the boundary type or boundary name set for the coupling interface when the grid is generated, and the number of nodes is calculated on each coupling interface. When encoding and decoding coupling data, the number of data of a single physical quantity needs to be determined by the number of coupling interfaces and the number of nodes on each coupling interface.

[0057] Specifically, in the aerodynamic / structural coupling analysis of the aircraft wing, in the fluid solver, all grid blocks of the flow field grid are traversed, and the wing coupling interface data is calculated according to the coupling interface boundary conditions set when the grid is generated. is 3, which are the main wing, horizontal tail and vertical tail respectively. At the same time, the number of coupling interface nodes is calculated according to the difference between the maximum and minimum values ​​of the boundary node numbers of the grid block where the coupling interface is located. , No. The number of nodes per wing is recorded as The structural solver uses the modal method to solve the structural control equations. The wing structure mesh unit type 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 3, the number of nodes on the coupling interface is also calculated through the grid file .

[0058] In this embodiment, the physical field dimension of the coupling data is determined according to the type of multi-physics field coupling analysis problem and the coupling data transmission method. For example, for most coupled problems, the physics dimension The value is 3, but for some special problems, such as plate coupling problems, since the thickness direction is much smaller than the other two directions, only the coordinates in the length and width directions are used to construct the transfer matrix, and the force and displacement are also only calculated using the thickness direction, that is, the plane normal. At this time, the physical field dimension The value is 2.

[0059] Specifically, in the fluid solver, the wing fluid grid is a three-dimensional structured grid, so the physical field dimension data in the fluid solver is set 3; 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 component 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.

[0060] Step S12: assigning a corresponding type association symbol to each physical quantity of each physical quantity type, assigning a corresponding storage location symbol to each physical quantity stored in the grid unit surface of the coupling interface, assigning a corresponding data type symbol to each physical quantity of scalar and / or vector data type, establishing a coding order for each coupling interface, establishing a coupling data coding order for each physical quantity type on the current coupling interface, and establishing a grid node coding order in the current coupling interface based on the position information of the grid nodes where each physical quantity is located on the current coupling interface.

[0061] In this embodiment, each target association symbol corresponding to each physical quantity type is screened from the preset physical quantity type association table to obtain a type association symbol sequence constructed with the target association symbol. It can be understood that the preset physical quantity type association table is as shown in Table 1.

[0062] Table 1

[0063]

[0064] According to the physical quantity types and physical quantity type association symbols of the fluid solver participating in the multidisciplinary coupling analysis as shown in Table 1, the corresponding physical quantity type symbols are matched according to the node coordinates and the node aerodynamic physical quantities, so the physical quantity symbol sequence of the fluid solver is {b, h}; similarly, according to the physical quantity types of the structural solver participating in the coupling analysis, the physical quantity symbol sequence is matched to {b, i}.

[0065] According to the physical laws and coupling iteration strategies in multidisciplinary coupling analysis, the subject solver needs to communicate multiple types of data with each subject solver. In order to achieve the communication of multiple types of physical quantity data to each solver based on the same coding format, without forcibly stipulating the coding order of physical quantities, and to improve the flexibility of coding, an association symbol is assigned to each type of physical quantity, and the physical quantity association symbol is encoded into the coupling data. When the subject solver decodes the coupling data, it can decode the specified type of coupling data according to the association symbol and the association relationship in Table 1. For example: in fluid-solid coupling analysis, the physical quantity types of the fluid solver coupling communication include: node coordinates, node forces, and node displacements. According to the association relationship in Table 1, the symbols of the physical quantity types are assigned respectively: {b, h, j}, and the physical quantity type association symbols are encoded into the coupling data. When the structural solver decodes the coupling data, it can obtain the coupling data type information of the fluid solver coupling interface according to the symbol sequence and Table 1.

[0066] In this embodiment, if the physical quantity is stored on a grid node in the grid unit surface of the coupling interface, the storage location symbol assigned to the physical quantity is 1, otherwise, the storage location symbol assigned to the physical quantity is 0; if the type of the physical quantity is a scalar type, the data type symbol assigned to the physical quantity is 0. It can be understood that in multidisciplinary coupling analysis, there are multiple grid space discretization methods and physical quantity types in each discipline solver, resulting in differences in the storage methods of physical quantities on the coupling interface of the coupled boundary grid unit; at the same time, when the coupled data is interpolated, the physical quantities at different positions will use interpolation matrices constructed with different types of base points. Therefore, the storage location information of the physical quantity must be included in the encoded data. Although there are great differences in the storage of physical quantity data in the unit, it can be divided into two types: face center storage and node storage. Face center storage refers to the storage of physical quantity data at one or more points inside the unit surface, such as: the centroid of the unit surface, the Lagrange interpolation point of the unit surface; node storage refers to the storage of physical quantity data at all nodes of the unit surface. When encoding coupled data, a symbol representing the storage location is assigned to each type of physical quantity to mark the physical quantity storage location. For example, when storing coupled data in a computer, a bit is used to mark the physical quantity storage location. 1 represents that the physical quantity is stored at the center of the face, and 0 represents that it is stored at the node. In this way, the corresponding storage location symbol is assigned according to the storage location.

[0067] Specifically, in the aerodynamic / structural coupling analysis of the aircraft wing, the physical quantity coupling data of the fluid solver and the structural solver participating in the coupling analysis are all located on the nodes of the unit surface, so the data storage location signs of all physical quantities of the fluid / structural solver are set to 0.

[0068] In this embodiment, if the type of the physical quantity is vector type, the data type symbol assigned to the physical quantity is 1. It can be understood that in multidisciplinary coupling analysis, when a discipline solver communicates coupling data to another discipline solver, the coupling data at a single point can be divided into scalar data and vector data in terms of dimension. For example: during fluid-solid thermal coupling analysis, the physical quantities that the fluid solver needs to transmit to the structural solver include: node coordinates, node forces, node temperatures, and unit surface heat flux density, etc. Node coordinates, node forces, and unit surface heat flux density are all vector data, and node temperature is scalar data. Since there are multiple types of physical quantity data in the discipline solver, the number of decoded data for different types of coupling data is inconsistent when decoding the coupling data. In order to mark the type of coupling data for a single physical quantity, a symbol representing the data type is assigned to each type of physical quantity, and the number of scalar coupling data on each coupling interface in the discipline solver is obtained. The number of coupled data with the vector type For example, when using a computer to store coupling data for fluid-solid-thermal coupling analysis, a bit is used to mark the type of physical quantity data. 1 represents the physical quantity as vector data, and 0 represents scalar data. The number of scalar data in the fluid solver 1, the number of vector data It is 3. In this way, the corresponding data type symbol is assigned according to the type to which it belongs.

[0069] Specifically, the node coordinates, node aerodynamic forces, and node displacement data are all vector data, so the data type symbols of the physical quantity coupling data in the fluid solver and the structural solver are all set to 1. Number of scalar physical quantity types in the fluid / structural solver 0, the number of vector physical quantity types is 2.

[0070] In this embodiment, the coding order of each coupling interface is established according to the structural topological grid type and / or the non-structural topological grid type, and a corresponding coupling interface associated symbol is assigned to each coupling interface. It can be understood that, according to the obtained coupling interface, the coding order of the coupling interface coding is established, and an associated symbol is assigned to each coupling interface. For example: the coupling interface sequence of the structural topological grid is in the order of the number of the grid block where the coupling interface is located from small to large. If there are multiple coupling interfaces on the same grid block, they can be arranged in order from small to large according to the number of the coupling interface in the grid block, and each coupling interface is assigned a number starting from 0 as an associated symbol; the coupling interface sequence of the non-structural topological grid is in the order of the arrangement of the boundary conditions in the subject solver grid, and the name of the coupling interface in the grid file is used as an associated symbol, that is, the coding order of the coupling interface is obtained.

[0071] Specifically, since the fluid solver and the structural solver both contain three coupling interfaces, the coding order for establishing the coupling interfaces is: main wing, horizontal tail, vertical tail, and the associated symbol sequence of the coupling interfaces is established respectively: {0, 1, 2}.

[0072] In this embodiment, according to the arrangement order of the physical quantity type-association symbol in the preset physical quantity type association table, the target arrangement order of each physical quantity type on the current coupling interface is obtained, so as to use the target arrangement order as the coupling data encoding order of each physical quantity type on the current coupling interface. It can be understood that the order of coupling data encoding of different physical types on a single coupling interface is established according to the arrangement order of the physical quantity types from top to bottom in Table 1. When encoding the physical quantity of a single coupling interface, the physical quantity of the coupling interface is encoded in sequence according to the established physical quantity encoding order. For example: when a typical aircraft is conducting fluid-solid thermal coupling analysis, all wings of the aircraft are used as coupling interfaces, including: main wing, horizontal tail, vertical tail. The types of physical quantities involved in the coupling analysis in the fluid solver include: node number, node coordinates, face center coordinates, node aerodynamic force, node displacement, node temperature, face center heat flux, etc., so the physical quantity encoding order established is {a, b, c, h, j, k, l}.

[0073] Specifically, according to the types of physical quantities involved in the coupling analysis in the established fluid solver and structural solver and their order in Table 1, the physical quantity coding order in the fluid solver is established as: node coordinates, node aerodynamic force; the physical quantity coding order in the structural solver is established as: node coordinates, node displacement.

[0074] In this embodiment, the coding order of the grid nodes in the current coupling interface is established as follows: In multidisciplinary coupling analysis, a single coupling interface usually contains coupling data of multiple physical quantity types on the same node set. In order to correspond to the same node order when encoding or decoding different types of coupling data, the order of the nodes or face centers of a single coupling interface is established. For example: in a three-dimensional structural topological grid, the nodes in a single grid block are numbered in the format of a triple (I, J, K). The node number I value in the coupling interface in the I direction is the same, and the node coding order is sorted from small to large according to the values ​​of J and K respectively. The coding order of the coupling interface nodes in the J and K directions is the same as that in the I direction. There are four nodes on the unit face where the face center is located. The face center coding order is based on the node with the smallest value in the J and K directions among the four nodes. The reference nodes of all face centers constitute a reference node set. The order of face centers is based on the order of each node in the node coding order in the reference node set as the face center coding order. In the unstructured topological grid, the order of the coupling interface nodes is coded in ascending order according to the global number of the nodes in the grid file; the order of the coding of the unit face centers of the coupling interface is coded in ascending order according to the global number of the unit face where they are located.

[0075] Specifically, all the physical quantities involved in the coupling analysis of the fluid solver and the structural solver are stored on the nodes of the unit surface, so it is only necessary to establish the coding order of a single coupling interface node. The flow field grid file format is the structured grid plot3D format, and the order of the nodes is in the order of the grid block where the node is located, the number of the grid surface on the grid block, and the number of the node on the grid surface from small to large. First, traverse all the grid blocks and their surfaces, collect the grid blocks and their grid surfaces containing the coupling interface boundary conditions, and then sort the boundary grid blocks in order from small to large by the number of the grid blocks, and the grid surfaces are in order of the order of the grid blocks in which they are located; the order of the nodes in the grid surface is in order from small to large according to the number values ​​of the nodes in the two directions of the grid surface, for example: on the grid surface in the I direction, the nodes are in the order of fixing the value of J first, then in the order of K from small to large, and then increasing the value of J. The order of the structural solver nodes is in the order of the node arrangement order in the finite element model file.

[0076] In this embodiment, in the multidisciplinary coupling iterative analysis, after the discipline solver advances a coupling time step in the physical time domain, it traverses all coupling interfaces of the solver to obtain a certain physical quantity data value. For example: when conducting aerodynamic / structural coupling analysis of a wing based on a first-order loose coupling strategy method, after the fluid solver advances a coupling time step, it traverses all unit surfaces on the wing surface, and obtains the unit surface aerodynamic force after performing area integration on the pressure at the center of the face, and distributes data processing methods such as linear weighting of the unit surface aerodynamic force to each node of the unit surface, and obtains the aerodynamic force data values ​​of all nodes on the wing surface, that is, the coupling data values.

[0077] Specifically, the fluid solver uses the unsteady flow method to solve and advance the flow field by one time step in physical time. The pressure area is integrated on the coupling interface unit surface to obtain the aerodynamic force on the unit surface. The linear weighted method is used to distribute the aerodynamic force of the unit surface to each node of the unit surface. All unit surfaces are traversed, and the aerodynamic forces on the same nodes are accumulated to obtain the aerodynamic data of all nodes. The structural solver uses the aerodynamic force of the wing surface node as an external load and solves the structural control equation based on the modal method. The coordinate data and displacement data of the structural surface nodes in the next time step are obtained as the coupling data value. The coupling data value of the current time step on each coupling interface is obtained to construct a coupling data sequence based on each coupling data value.

[0078] Step S13: Encode the storage location symbol and the data type symbol to obtain a first encoding result, traverse and encode all grid node information of a single coupling interface in the grid node encoding order to obtain a corresponding second encoding result, and construct a third encoding result of the coupling data encoding characterizing a single physical quantity type on a single coupling interface according to the order of the first encoding result and the second encoding result.

[0079] In this embodiment, the encoding of single physical type coupling data refers to encoding the coupling data of a single physical type on a single coupling interface. The encoding format of the single physical type coupling data is: storage location symbol, data type symbol, coupling data value. The encoding order of the coupling data value is encoded into the coupling data in sequence according to the established node / face center order. When the data type is vector data, the component encoding order of the coupling data of each point is divided into the following according to the physical field dimension: When it is 2, the encoding order is x, y; when When it is 3, the encoding order is x, y, z. For example: in a typical aircraft aerodynamic / structural coupling analysis, the node coordinate data of the main wing surface in the three-dimensional unstructured grid fluid solver is encoded and the coupling data is stored in computer memory. The node coordinates are stored on the nodes in the discrete grid, so the data storage location symbol is 0; the node coordinates are vector data, so the data type symbol is 1. Allocate a byte of memory space and store the data storage location symbol and the data type symbol in the first two bits. Establish a node encoding order based on the order from small to large in the coordinate sequence of the grid file for all nodes on the main wing. And fill the node coordinate components into the memory in sequence according to the established node order, and the coordinates of each node are in the order of x, y, z.

[0080] Specifically, in the wing aerodynamic / structural coupling analysis, in order to maintain the consistency of the encoding or decoding of each solver during the coupling analysis process, a multi-step prediction-correction coupling strategy format with second-order coupling accuracy is adopted. The fluid solver or structural solver first encodes the physical quantity data storage location symbol and data type symbol set in step five and step six; then, according to the coupling interface node encoding sequence established in step ten, traverse all the nodes of the coupling interface, and encode the node coordinates after the data storage location symbol and data type symbol according to the established node encoding sequence, and the coordinate components of each node are encoded in the order of x, y, and z; similarly, traverse all the nodes on the coupling interface, encode the node aerodynamic force or node displacement after the coordinate data according to the established node encoding sequence, and the aerodynamic force component of each node is encoded in the order of x, y, and z.

[0081] Step S14: Encode the physical quantity quantity into a coupling data sequence, and fill the type association symbol into the encoded coupling data sequence, and then encode the third encoding result into the encoded coupling data sequence according to the coupling data encoding order to obtain a fourth encoding result after coupling data encoding representing all physical quantity types on a single coupling interface; wherein the coupling data sequence is constructed based on the values ​​of each physical quantity.

[0082] In this embodiment, the coupling data encoding format of a single coupling interface is: physical quantity quantity , physical quantity type symbol sequence, scalar physical quantity type number , Number of vector physical quantity types , all physical quantity data values. First, the physical quantity Encode into the data sequence; then fill the physical quantity type associated symbol sequence on the coupling interface into the data sequence according to the established coupling data encoding sequence, and then encode the established single physical type encoded coupling data into the coupling data according to the physical quantity encoding sequence. For example: in the aircraft aerodynamic / structural coupling analysis, establish the coupling interface encoding sequence: main wing, horizontal tail, vertical tail, and the physical quantity type of each coupling interface is: node coordinates, face center coordinates, node aerodynamic force, node temperature, face center heat flux density, etc. Number of scalar physical quantity types 1, the number of vector physical quantity types is 4, and the encoding content of each coupling interface is: the number of physical quantities is 5, the type symbol sequence is {a, b, c, h, j, k, l}, and the number of scalar physical quantity types is , Number of vector physical quantity types , nodal coordinate coupling data, face-center coordinate coupling data, nodal aerodynamic coupling data, nodal temperature coupling data, face-center heat flux density coupling data.

[0083] Specifically, encoding the data of a single coupling interface is mainly divided into three steps: First, the fluid solver or structural solver converts the physical quantity Encoded into the coupling data sequence; then filled into the coupling data sequence according to the associated symbol sequence of the node coordinates, node aerodynamic force or node displacement established in step 4, and finally encoded into the coupling data sequence according to the physical quantity coding order of the established single physical type encoding, forming a coupling data sequence of a single coupling interface flow field solver or structural solver of the main wing, horizontal tail wing and vertical tail wing.

[0084] Step S15: Encode the basic information of the coupling data including the physical field dimension, the number of coupling interfaces, and the coding order of the coupling interfaces to obtain a basic information coding result, and fill the fourth coding result after the basic information coding result to obtain a target coding result.

[0085] In this embodiment, the solver coupling data encoding format is: physical field dimension , Number of coupling interfaces , coupling interface symbol sequence, all coupling interface data values. First, the acquired physical field dimension data into the coupling data; then, the number of coupling interfaces obtained Encode into the coupling data. Encode the established coupling interface symbol sequence into the coupling data, and encode the completed coupling data of a single coupling interface into the coupling data in sequence according to the coupling interface coding order. At this point, the coupling data encoding of all coupling interfaces of the subject solver is completed. For example: in the aircraft aerodynamic / structural coupling analysis, establish the coupling interface coding order: main wing, horizontal tail, vertical tail, and establish the coupling interface associated symbol sequence: {0, 1, 2}. The coupling data after encoding a single coupling interface is established. According to this step, the solver coding format encodes the content of the fluid solver coupling data as follows: physical field dimension is 3, the number of coupling interfaces is 3, the coupling interface symbol sequence is {0, 1, 2}, main wing coupling data, horizontal tail coupling data, vertical tail coupling data.

[0086] Specifically, the physical dimension data of the calculation , the number of coupling interfaces calculated , the established coupling interface symbol sequence {0, 1, 2} and other basic coupling data information are encoded; then the established coupling data of each wing is encoded, and the coupling data of the main wing, horizontal tail, and vertical tail are filled in after the basic coupling data information. So far, the fluid solver and the structural solver have established the coupling data of the coupling interface of each wing. Figure 2 As shown, the first is the single physical type coupling data encoding format, specifically, the node coordinate data: , , , , Here and Labels for identifying locations and attributes, , , etc. are specific scalar values, i.e., node coordinate data. Node aerodynamic data: , , , , Similar to scalar data, and It's a label. , , The aerodynamic data are the values ​​on different coordinate axes, and then the single interface coupling data encoding format. Then the single interface coupling data encoding format is transitioned from the above-mentioned single physical type coupling data encoding format to the single interface coupling data encoding format, and finally to the fluid solver coupling data encoding format, and the single coupling interface data on each main wing coupling interface, horizontal tail coupling interface and vertical tail coupling interface are encoded as a whole to obtain the corresponding encoding results.

[0087] In this embodiment, after obtaining the target coding result, it also includes: communicating the target coding result to other subject solvers according to the communication mapping relationship between the communication domains of each subject solver. It can be understood that after the coupling data is encoded, the coupling data encoded by the current solver is communicated to other subject solvers according to the communication mapping relationship between the subject solver communication domains established by multi-disciplinary coupling. For example: during fluid-solid coupling analysis, after the fluid solver encodes the node coordinates and the node aerodynamic force, it communicates the coupling data to the structural solver through the solver communication mapping relationship. The structural solver receives the coupling data, calculates the aerodynamic force interpolation matrix based on the fluid node coordinates, completes the node aerodynamic force interpolation, and generates the node aerodynamic force on the surface of the structure.

[0088] It can be seen that the present application discloses a multidisciplinary coupling surface scalable communication data encoding method, including: determining the physical quantity type and physical quantity quantity of all coupling data involved in the coupling communication in the multidisciplinary coupling simulation analysis problem, and calculating the number of coupling interfaces of the coupling interfaces involved in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface, and determining the physical field dimension of the coupling data; assigning a corresponding type association symbol to each physical quantity of each physical quantity type, assigning a corresponding storage location symbol to each physical quantity stored in the grid unit surface of the coupling interface, assigning a corresponding data type symbol to each physical quantity of scalar and / or vector data types, establishing a coding order for each coupling interface, establishing a coupling data coding order for each physical quantity type on the current coupling interface, and establishing a grid node coding order in the current coupling interface for the position information of the grid node where each physical quantity is located on the current coupling interface; encoding the storage location symbol and the data type symbol to obtain a first The encoding result is obtained by traversing and encoding all grid node information of a single coupling interface in accordance with the grid node encoding order to obtain a corresponding second encoding result, and constructing a third encoding result of coupling data encoding representing a single physical quantity type on a single coupling interface according to the order of the first encoding result and the second encoding result; encoding the number of physical quantities into a coupling data sequence, and filling the type association symbol into the encoded coupling data sequence, and then encoding the third encoding result into the encoded coupling data sequence in accordance with the coupling data encoding order to obtain a fourth encoding result after coupling data encoding representing all physical quantity types on a single coupling interface; wherein the coupling data sequence is constructed based on the values ​​of each physical quantity; encoding the coupling data basic information including the physical field dimension, the number of coupling interfaces, and the encoding order of the coupling interfaces to obtain a basic information encoding result, and filling the fourth encoding result into the basic information encoding result to obtain a target encoding result. It can be seen that through the allocation process of type association symbols and storage location symbols, it is possible to process coupled data of various physical quantity types and different data types (scalar type, vector type), meet the needs of complex multidisciplinary coupling analysis scenarios such as fluid-solid coupling and fluid-solid thermal coupling, reflect the processing ability of various data types, and enhance flexibility. By establishing the coding order of the coupling interface, the coding order of the coupling data, and the coding order of the grid nodes, different grid structures (structural topological grids and non-structural topological grids) and the storage locations (nodes or face centers) of physical quantities in the grid unit surface of the coupling interface are further considered. The coding order is established for these differences, so that the method can adapt to different grid space discretization methods and physical quantity storage methods of various discipline solvers, and improve adaptability.Finally, based on the above information and through the preset coding method, the basic information of the physical quantities and coupling information of various physical quantity types of each coupling interface are uniformly encoded, and then the data communication on multiple coupling interfaces can be completed through a single message passing interface call. The number of calls to the message passing interface is greatly reduced, and the problem of increased time cost caused by multiple calls to the interface is reduced, especially when there are a large number of calculation processes, which effectively avoids the problem of communication time accounting for too large a proportion of the overall calculation time, thereby improving the efficiency of multidisciplinary coupling simulation analysis.

[0089] Reference Figure 3 As shown, the embodiment of the present invention discloses a multi-disciplinary coupling surface scalable communication data decoding method. Compared with the previous embodiment, this embodiment further illustrates and optimizes the technical solution. Specifically:

[0090] Step S21: obtaining a target coding result encoded by a subject solver using the multi-disciplinary coupling surface scalable communication data coding method.

[0091] 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 obtains the target encoding result after the structural encoder encodes, and the structural solver obtains the coupling data after the fluid encoder encodes. The subject solvers receive the encoded coupling data, that is, the target encoding result.

[0092] Step S22: Initially decode the target coding result to obtain the physical field dimension of the coupling data, the number of coupling interfaces, and the coding order of the coupling interfaces in sequence.

[0093] In this embodiment, the coupling data is decoded according to the encoding format of the coupling data. First, the first two data are obtained from the target encoding result, which are: physical field dimension data Data on the number of coupling interfaces ; Then get data, i.e. the coupling interface symbol sequence. Then loop Number of times, The next cycle is based on the correspondence between the established coupling interface and its symbol. The coupling interface symbols are restored to the corresponding coupling interfaces, and the first coupling interface is decoded from the encoded coupling data. Coupling data of a coupling interface.

[0094] Specifically, after the structural solver receives the encoded coupling data transmitted by the fluid solver, it starts to decode the coupling data of the fluid solver. The decoding operation is divided into three steps: First, the first and second data are obtained from the coupling data in turn: the physical field dimension data of the flow field Data on the number of coupling interfaces ; Then, obtain the flow field coupling interface symbol sequence {0, 1, 2} from the coupling data; Finally, decode the coupling interfaces one by one according to the obtained coupling interface symbol sequence. Similarly, after the fluid solver receives the encoded coupling data transmitted by the structural solver, it starts to decode the coupling data of the structural solver. Similarly, first obtain the physical field dimension data of the structural solver Data on the number of coupling interfaces ; Then, the structural coupling interface symbol sequence {0, 1, 2} is obtained, and finally the coupling interfaces are decoded one by one according to the obtained coupling interface symbol sequence. The decoding process of the fluid solver is as follows Figure 4 shown.

[0095] Step S23: perform secondary decoding on the fourth coding result of each coupling interface according to the coding order of the coupling interface to obtain the number of physical quantities, and jump to execute the step of secondary decoding of the fourth coding result of each coupling interface according to the coding order of the coupling interface the same number of times as the number of physical quantities to obtain a type-associated symbol sequence, and determine the third coding result corresponding to each type-associated symbol in the type-associated symbol sequence according to a preset physical quantity type association table.

[0096] In this embodiment, when obtaining the When the data of a coupling interface is obtained, the data of the current coupling interface is decoded. First, the first data is obtained from the encoded coupling data, that is, the number of physical quantities of the coupling interface. ; Then obtain the physical quantity type symbol sequence, and convert the symbol sequence into a physical quantity type sequence according to Table 1. The last loop Number of times, The next loop obtains Coupling data for each physical quantity type.

[0097] Specifically, when the fluid solver or structural solver decodes the data of a single coupling interface (main wing / horizontal tail / vertical tail) in sequence, first, the physical quantity of the coupling interface is obtained from the coupling data sequence. ; then loop Each time a physical quantity type symbol is obtained, a physical quantity symbol sequence of the coupling interface is formed, and the symbol sequence is converted into a physical quantity type sequence according to Table 1. times, each time obtaining a coupling data sequence of a physical quantity type on the coupling interface.

[0098] Step S24: Decoding the third encoding result to obtain a storage location symbol and a data type symbol as well as information of each grid node on a single coupling interface to obtain decoded coupling data.

[0099] In this embodiment, when decoding a single physical quantity data of a single coupling interface, the coupling data storage position symbol is first obtained from the encoded coupling data. If the symbol is 0, the total number of coupling interface nodes is obtained from the coupling data. , and set the number of points equal ; If the symbol is 1, the total number of coupling interface face centers is obtained from the coupling data , and set the number of points equal Then get the data type sign of the coupled data. If the sign is 0, the data component of each point If the symbol is 1, the data component of each point for . Calculate the number of data of this physical quantity . Then obtain from the coupling data Data, according to Each data corresponds to a point, the data is mapped to each node or face center, and the physical quantity data is assigned to the physical quantity of the subject solver.

[0100] Specifically, the decoding operation of coupling data of a physical quantity type on a single coupling interface can be divided into three main steps: First, obtain the first two data from the coupling data: the data storage location symbol and the data type symbol, which are 0 and 1 for the wing aerodynamic / structural coupling analysis problem respectively; Since the coupling data is stored on the grid nodes and all physical quantities are vector data, the total number of coupling interface nodes is obtained And the coupled data component at each point for Calculate the total amount of data for this physical quantity . Then obtain from the coupling data data to form a coupled data set. Each data corresponds to a point, the data set is distributed to each node, and the decoded coupled data set is assigned to the physical quantity of the current subject solver.

[0101] In this embodiment, the fluid solver takes the fluid node coordinates and the structural node coordinates as the base points, and establishes the displacement interpolation matrix G based on the radial basis function (RBF). The structural node displacement is taken as a known quantity, 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 takes the fluid node coordinates and the structural node coordinates as the base points, and establishes the aerodynamic interpolation matrix F based on the principle of virtual work. The fluid node aerodynamic force is taken as a known quantity, and the matrix-vector product is performed with the aerodynamic interpolation matrix to obtain the aerodynamic force of the structural grid wing surface node.

[0102] After obtaining the node displacements of the fluid grid wing surface, the grid deformation methods such as RBF-TFI are used to calculate the internal node coordinates of the grid. The steps of using 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 based on 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. In this way, with the help of mesh deformation methods such as RBF-TFI, based on the node displacement information obtained after decoding, the internal node coordinates of the mesh can be accurately calculated, and the geometric information of the mesh (such as node coordinates, centroid coordinates, volume, area, area vector) and Jacobi transformation coefficient matrix can be recalculated. This provides an effective solution for dealing with the mesh deformation problem involved in complex multidisciplinary coupling analysis, which can better adapt to the changes in mesh morphology caused by coupling and ensure the accuracy and reliability of the analysis.

[0103] When the number of coupling iteration steps does not exceed the set maximum number of coupling iteration steps, the physical time accumulates a coupling time step, and the coupling iteration step increases by 1. Jump to execute the fluid and structure discipline solver time advancement until the maximum number of coupling analysis steps is reached. In this way, in each coupling iteration, by accurately decoding to obtain the required data, after completing the corresponding calculation and analysis steps, when the number of coupling iteration steps does not exceed the set maximum number of coupling iteration steps, it can achieve physical time accumulation of a coupling time step, the coupling iteration step increases by 1, and jump to execute the fluid and structure discipline solver time advancement until the maximum number of coupling analysis steps is reached. This good adaptability to the coupling iteration process allows multidisciplinary coupling analysis to be carried out in an orderly manner according to the predetermined iteration rules, ensuring the integrity and effectiveness of the entire analysis process.

[0104] It can be seen that the initial decoding of the coupling data encoding format can accurately obtain key information such as the physical field dimension of the coupling data, the number of coupling interfaces, and the encoding order of the coupling interfaces. The subsequent steps further decode the encoding results of each part on this basis, so as to completely and accurately restore all kinds of information contained in the original encoding, including the number of physical quantities, type-related symbol sequences, storage location symbols, data type symbols, and information of each grid node, etc., ensuring that the encoded data can be accurately restored during the decoding process, providing an accurate data basis for multidisciplinary coupling analysis.

[0105] Reference Figure 5 As shown, the present invention also discloses a multi-disciplinary coupling surface scalable communication data encoding device, including:

[0106] The information determination module 11 is used to determine the physical quantity type and physical quantity quantity of all coupling data involved in the coupling communication in the multidisciplinary coupling simulation analysis problem, and calculate the number of coupling interfaces of the coupling interfaces involved in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface to determine the physical field dimension of the coupling data;

[0107] An information allocation module 12 is used to allocate a corresponding type association symbol to each physical quantity of each physical quantity type, allocate a corresponding storage location symbol to each physical quantity stored in the grid unit surface of the coupling interface, allocate a corresponding data type symbol to each physical quantity of scalar type and / or vector type data type, establish a coding order for each coupling interface, establish a coupling data coding order for each physical quantity type on the current coupling interface, and establish a grid node coding order in the current coupling interface based on the position information of the grid node where each physical quantity is located on the current coupling interface;

[0108] A first encoding module 13 is used to encode the storage location symbol and the data type symbol to obtain a first encoding result, traverse and encode all grid node information of a single coupling interface in the grid node encoding order to obtain a corresponding second encoding result, and construct a third encoding result of coupling data encoding representing a single physical quantity type on a single coupling interface according to the order of the first encoding result and the second encoding result;

[0109] A second encoding module 14 is used to encode the physical quantity quantity into a coupling data sequence, fill the type association symbol into the encoded coupling data sequence, and then encode the third encoding result into the encoded coupling data sequence according to the coupling data encoding order to obtain a fourth encoding result after coupling data encoding representing all physical quantity types on a single coupling interface; wherein the coupling data sequence is constructed based on the values ​​of each physical quantity;

[0110] The third encoding module 15 is used to encode the basic information of the coupling data including the physical field dimension, the number of coupling interfaces, and the encoding order of the coupling interfaces to obtain a basic information encoding result, and fill the fourth encoding result after the basic information encoding result to obtain a target encoding result.

[0111] It can be seen that the present application discloses determining the physical quantity type and the number of physical quantities of all coupling data involved in coupling communication in a multidisciplinary coupling simulation analysis problem, and calculating the number of coupling interfaces of the coupling interfaces involved in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface, and determining the physical field dimension of the coupling data; assigning a corresponding type association symbol to each physical quantity of each physical quantity type, assigning a corresponding storage location symbol to each physical quantity stored in the grid unit surface of the coupling interface, assigning a corresponding data type symbol to each physical quantity of scalar and / or vector data type, establishing a coding order for each coupling interface, establishing a coupling data coding order for each physical quantity type on the current coupling interface, and establishing a grid node coding order in the current coupling interface for the position information of the grid node where each physical quantity is located on the current coupling interface; encoding the storage location symbol and the data type symbol to obtain a first coding result, and according to the grid The node coding sequence traverses and encodes all grid node information of a single coupling interface to obtain a corresponding second coding result, and constructs a third coding result of coupling data coding representing a single physical quantity type on a single coupling interface according to the order of the first coding result and the second coding result; encodes the number of physical quantities into a coupling data sequence, and fills the type association symbol into the encoded coupling data sequence, and then encodes the third coding result into the encoded coupling data sequence according to the coupling data coding sequence to obtain a fourth coding result after coupling data coding representing all physical quantity types on a single coupling interface; wherein the coupling data sequence is constructed based on the values ​​of each physical quantity; encodes the coupling data basic information including the physical field dimension, the number of coupling interfaces, and the coding order of the coupling interfaces to obtain a basic information coding result, and fills the fourth coding result after the basic information coding result to obtain a target coding result. It can be seen that through the allocation process of type association symbols and storage location symbols, it is possible to process coupling data of various physical quantity types and different data types (scalar type, vector type), meet the needs of complex multidisciplinary coupling analysis scenarios such as fluid-solid coupling and fluid-solid thermal coupling, reflect the processing ability of multiple data types, and enhance flexibility. By establishing the coding order of the coupling interface, the coding order of the coupling data, and the coding order of the grid nodes, different grid structures (structural topological grids and non-structural topological grids) and the storage locations (nodes or face centers) of the physical quantities in the grid unit surface of the coupling interface are further considered. The coding order is established for these differences, so that the method can adapt to different grid space discretization methods and physical quantity storage methods of various disciplinary solvers, and improve adaptability. Finally, based on the above information and through the preset coding method, the basic information of the physical quantities and coupling information of various physical quantity types of each coupling interface are uniformly encoded, and then the data communication on multiple coupling interfaces can be completed through a single message passing interface call.The number of calls to the message passing interface is greatly reduced, reducing the problem of increased time cost caused by multiple calls to the interface. Especially when there are a large number of computing processes, it effectively avoids the problem of communication time accounting for too large a proportion of the overall computing time, thereby improving the efficiency of multidisciplinary coupled simulation analysis.

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

[0113] Figure 6 The 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.

[0114] 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.

[0115] 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.

[0116] 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.

[0117] 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.

[0118] 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.

[0119] 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.

[0120] 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.

[0121] 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.

[0122] 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 multidisciplinary coupling surface scalable communication data encoding method, characterized in that: include: Determine the physical quantity types and physical quantity quantities of all coupling data involved in coupling communication in the multidisciplinary coupling simulation analysis problem, calculate the number of coupling interfaces involved in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface, and determine the physical field dimension of the coupling data; Assigning a corresponding type association symbol to each physical quantity of each physical quantity type, assigning a corresponding storage location symbol to each physical quantity stored in the grid unit surface of the coupling interface, assigning a corresponding data type symbol to each physical quantity of scalar type and / or vector type data type, establishing a coding order for each coupling interface, establishing a coupling data coding order for each physical quantity type on the current coupling interface, and establishing a grid node coding order in the current coupling interface for the location information of the grid node where each physical quantity is located on the current coupling interface; Encode the storage location symbol and the data type symbol to obtain a first encoding result, traverse and encode all grid node information of a single coupling interface in the grid node encoding order to obtain a corresponding second encoding result, and construct a third encoding result of coupling data encoding characterizing a single physical quantity type on a single coupling interface according to the order of the first encoding result and the second encoding result; Encoding the physical quantity quantity into a coupling data sequence, and filling the type association symbol into the encoded coupling data sequence, and then encoding the third encoding result into the encoded coupling data sequence according to the coupling data encoding order, so as to obtain a fourth encoding result after coupling data encoding representing all physical quantity types on a single coupling interface; wherein the coupling data sequence is constructed based on the values ​​of each physical quantity; The basic information of the coupling data including the physical field dimension, the number of coupling interfaces, and the coding order of the coupling interfaces is encoded to obtain a basic information coding result, and the fourth coding result is filled after the basic information coding result to obtain a target coding result.

2. The multidisciplinary coupling surface scalable communication data encoding method according to claim 1, characterized in that: The calculation of the number of coupling interfaces involved in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface includes: According to the grid information of the subject solver for solving the multidisciplinary coupling simulation analysis problem, all boundaries of the corresponding current grid are traversed, and the number of coupling interfaces of the coupling interfaces participating in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface are calculated according to the coupling interface boundary conditions set when the current grid is generated.

3. The multidisciplinary coupling surface scalable communication data encoding method according to claim 1, characterized in that: The step of assigning a corresponding type association symbol to each physical quantity of each physical quantity type includes: The target associated symbols corresponding to the physical quantity types are screened from the preset physical quantity type association table to obtain a type associated symbol sequence constructed with the target associated symbols.

4. The multidisciplinary coupling surface scalable communication data encoding method according to claim 3 is characterized in that: The step of establishing the coupling data coding sequence of each physical quantity type on the current coupling interface includes: According to the arrangement order of the physical quantity type-association symbols in the preset physical quantity type association table, the target arrangement order of each of the physical quantity types on the current coupling interface is obtained, so as to use the target arrangement order as the coupling data encoding order of each of the physical quantity types on the current coupling interface.

5. The multidisciplinary coupling surface scalable communication data encoding method according to claim 1, characterized in that: The step of assigning a corresponding storage location symbol to each of the physical quantities stored in the grid unit surface of the coupling interface, and assigning a corresponding data type symbol to each of the physical quantities of scalar and / or vector data types, comprises: If the physical quantity is stored on a grid node in a grid unit surface of the coupling interface, the storage location symbol assigned to the physical quantity is 1, otherwise, the storage location symbol assigned to the physical quantity is 0; If the type of the physical quantity is a scalar type, the data type symbol assigned to the physical quantity is 0; if the type of the physical quantity is a vector type, the data type symbol assigned to the physical quantity is 1.

6. The multidisciplinary coupling surface scalable communication data encoding method according to claim 1, characterized in that: The coding sequence of establishing each of the coupling interfaces includes: A coding sequence of each coupling interface is established according to the structural topological grid type and / or the non-structural topological grid type, and a corresponding coupling interface association symbol is allocated to each coupling interface.

7. The multidisciplinary coupling surface scalable communication data encoding method according to claim 1, characterized in that: After obtaining the target encoding result, the method further includes: According to the communication mapping relationship between the communication domains of each subject solver, the target encoding result is communicated to other subject solvers.

8. A multidisciplinary coupling surface scalable communication data decoding method, characterized in that: include: Obtaining a target encoding result encoded by a subject solver using the multi-disciplinary coupling surface scalable communication data encoding method according to any one of claims 1 to 7; Initially decoding the target encoding result to sequentially obtain the physical field dimension of the coupling data, the number of coupling interfaces, and the encoding order of the coupling interfaces; According to the coding sequence of the coupling interfaces, the fourth coding result of each coupling interface is decoded twice to obtain the number of physical quantities, and the step of decoding the fourth coding result of each coupling interface according to the coding sequence of the coupling interfaces is executed twice for the same number of times as the number of physical quantities to obtain a type association symbol sequence, and the third coding result corresponding to each type association symbol in the type association symbol sequence is determined according to a preset physical quantity type association table; The third encoding result is decoded to obtain a storage location symbol and a data type symbol as well as information of each grid node on a single coupling interface to obtain decoded coupling data.

9. A multidisciplinary coupling surface scalable communication data encoding device, characterized in that: include: An information determination module is used to determine the physical quantity type and physical quantity quantity of all coupling data involved in coupling communication in a multidisciplinary coupling simulation analysis problem, and calculate the number of coupling interfaces of the coupling interfaces involved in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface, and determine the physical field dimension of the coupling data; An information allocation module, used to allocate a corresponding type association symbol to each physical quantity of each physical quantity type, allocate a corresponding storage location symbol to each physical quantity stored in a grid unit surface of the coupling interface, allocate a corresponding data type symbol to each physical quantity of a scalar and / or vector data type, establish a coding sequence for each coupling interface, establish a coupling data coding sequence for each physical quantity type on a current coupling interface, and establish a grid node coding sequence in a current coupling interface based on the position information of the grid node where each physical quantity is located on the current coupling interface; A first encoding module is used to encode the storage location symbol and the data type symbol to obtain a first encoding result, traverse and encode all grid node information of a single coupling interface in the grid node encoding order to obtain a corresponding second encoding result, and construct a third encoding result of coupling data encoding representing a single physical quantity type on a single coupling interface according to the order of the first encoding result and the second encoding result; a second encoding module, for encoding the physical quantity quantity into a coupling data sequence, and filling the type association symbol into the encoded coupling data sequence, and then encoding the third encoding result into the encoded coupling data sequence according to the coupling data encoding order, so as to obtain a fourth encoding result after coupling data encoding representing all types of physical quantities on a single coupling interface; wherein the coupling data sequence is constructed based on the values ​​of each physical quantity; The third encoding module is used to encode the basic information of the coupling data including the physical field dimension, the number of coupling interfaces, and the encoding order of the coupling interfaces to obtain a basic information encoding result, and fill the fourth encoding result after the basic information encoding result to obtain a target encoding result.

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