Multidisciplinary coupling surface scalable communication data encoding method, decoding method, device and equipment

Abstract

The present application discloses a multidisciplinary coupling surface scalable communication data encoding method, decoding method, device and equipment, which relates to the field of multidisciplinary coupling numerical simulation, including: by encoding multiple coupling data values ​​of arbitrary physical quantity types and categories to be communicated on multiple coupling interfaces and basic information of coupling data into a coupling data array, data communication on multiple coupling interfaces can be completed by calling a message passing interface once, which greatly reduces the number of calls to the message passing interface, and physical data in subject solvers of different numerical simulation methods can carry out coupling data communication in the same encoding and decoding format, thereby realizing arbitrary expansion of coupling data types, reducing the difficulty of multidisciplinary coupling analysis software development, and improving the efficiency of multidisciplinary coupling analysis.

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] Currently, the partitioned coupling solution strategy for each discipline can use mature single-discipline solvers, greatly reducing the difficulty of solver development and being applicable to extremely complex application scenarios, becoming the main way to achieve multidisciplinary coupled simulation. In the partitioned coupling solution strategy, different discipline solvers operate independently, and the iterative update of different physical field data is achieved through the exchange of interface information. The coupling of multiple components through the exchange of coupling information across multiple interfaces has become a common phenomenon in multidisciplinary coupling analysis scenarios. For example, in the aerodynamic / structural coupling analysis of aircraft components such as the main wing, ailerons, and tail, the fluid solver and the structural solver need to exchange data from multiple wing coupling interfaces.

[0003] Analyzing multiple coupling interfaces presents challenges, such as the diverse physical quantity types and data types of the coupled data. Furthermore, the storage locations of the coupled 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 the coupled data. Therefore, when communicating coupled data across multiple coupling interfaces using a message-passing interface, not only the coupled data itself must be communicated, but also basic information such as the coupled data type, type, and distribution location on the discrete grid must be communicated to the other discipline solvers. However, the message-passing interface can only send and receive a single type of data at a time. When the number of coupled data types is large, multiple calls to the message-passing interface are required, requiring the solver to perform multiple communication interface preparations, increasing overall communication time. When the number of computational processes is high, the communication time cost can account for a significant proportion of the overall computational time, thereby reducing the efficiency of multidisciplinary coupled simulation analysis. Summary of the Invention

[0004] In view of this, the purpose of the present invention is to provide a multidisciplinary coupling interface scalable communication data encoding method, decoding method, device and equipment, which can reduce the call of 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 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 a grid unit surface of the coupling interface, assigning a corresponding data type symbol to each physical quantity of a 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 node where each physical quantity is located on the current coupling interface;

[0008] Encoding the storage location symbol and the data type symbol to obtain a first encoding result, traversing and encoding all grid node information of a single coupling interface in the grid node encoding order to obtain a corresponding second encoding result, and constructing a third encoding result representing the coupling data encoding of a single physical quantity type on the 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 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 participating in the multidisciplinary coupling analysis and the number of grid nodes on each coupling interface includes:

[0012] According to the grid information of the discipline 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, assigning a corresponding type association symbol to each physical quantity of each physical quantity type includes:

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

[0015] Optionally, the step of establishing a coupling data coding sequence for 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 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.

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

[0018] If the physical quantity is stored on a grid node in a grid cell surface of the coupling interface, the storage location symbol of the physical quantity is assigned to be 1; otherwise, the storage location symbol of the physical quantity is assigned to be 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, establishing the coding sequence of each coupling interface includes:

[0021] A coding sequence of each coupling interface is established according to the structured topology grid type and / or the unstructured topology 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 coding 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] performing secondary decoding on the fourth coding result of each coupling interface according to the coding order of the coupling interfaces to obtain a physical quantity quantity, and skipping to execute the step of secondary decoding on the fourth coding result of each coupling interface according to the coding order of the coupling interfaces a number of times equal to the physical quantity quantity to obtain a type association symbol sequence, and determining a third coding result corresponding to each type association symbol in the type association symbol sequence 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, so as 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 types and physical quantity quantities of all coupling data involved in the 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;

[0031] an information allocation module, configured 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 cell 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 order for each coupling interface, establish a coupling data coding order for each physical quantity type on a 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;

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

[0033] a second encoding module, configured 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, 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;

[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, comprising:

[0036] Memory, used to store computer programs;

[0037] A 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 number of all coupling data participating in the coupling communication in the multidisciplinary coupling simulation analysis problem, and calculating 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, 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 encoding result, 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, 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. As can be seen, the process of assigning type association symbols and storage location symbols enables the processing of coupled data of various physical quantity types and different data types (scalar and vector), meeting the needs of complex multidisciplinary coupled analysis scenarios such as fluid-structure coupling and fluid-structure thermal coupling. This demonstrates the ability to handle multiple data types and enhances flexibility. By establishing a coding order for the coupling interface, the coupling data, and the grid node coding order, the method further considers different grid structures (structural topological grids and unstructural topological grids) and the storage location of physical quantities in the coupling interface grid element faces (nodes or face centers). By establishing a coding order based on these differences, the method can adapt to the different grid space discretization methods and physical quantity storage methods of various disciplines solvers, thus improving adaptability.Finally, based on this information and using a preset encoding method, the basic information of the physical quantities and coupling information of various types of physical quantities on each coupling interface is uniformly encoded. This allows data communication across multiple coupling interfaces to be completed with a single message-passing interface call. This significantly reduces the number of message-passing interface calls and the time cost associated with multiple interface calls. This effectively avoids the problem of communication time dominating the overall computational time, especially when there are many computational processes, thereby improving the efficiency of multidisciplinary coupled 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 following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are merely embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying any creative work.

[0040] Figure 1 This is a flow chart of a multidisciplinary coupling surface scalable communication data encoding method disclosed in this application;

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

[0042] Figure 3 This is a flow chart of a multidisciplinary coupling surface scalable communication data decoding method disclosed in this application;

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

[0044] Figure 5 This is a schematic structural diagram of a multidisciplinary coupling surface scalable 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 embodiments described 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 efforts are within the scope of protection of the present invention.

[0047] With the rapid advancement of computer hardware and software capabilities and the continuous development and improvement of interdisciplinary numerical simulation methods, numerically-based CAE (Computer-Aided Engineering) industrial software has gradually become a key driver of design paradigm transformation in equipment manufacturing, particularly in the aviation and aerospace sectors, playing a crucial role in product design, performance evaluation, and maintenance. Single-discipline simulation software has matured and, in many areas, can replace physical simulation or testing, providing highly refined predictive data for product design. However, with the increasing sophistication and requirements of modern product design, numerous interdisciplinary coupling issues have emerged, becoming a key challenge hindering product performance. Improving the quality and performance of modern products urgently requires the support of highly reliable multidisciplinary coupled simulation and analysis tools.

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

[0049] Analyzing multiple coupling interfaces presents challenges such as the diverse physical quantity types and data types of coupled data. Furthermore, the storage locations of coupled 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 the coupled data. Therefore, when communicating coupled data across multiple coupling interfaces using message-passing interfaces, not only the coupled data itself must be communicated, but also basic information such as the coupled data type, category, and distribution location on the discrete grid must be communicated to the other disciplinary solvers. However, message-passing interfaces can only send and receive a single type of data at a time. When the number of coupled data types is large, multiple message-passing interfaces must be called, requiring the solver to perform multiple communication interface preparations, increasing overall communication time. When a large number of computational processes are involved, communication time can account for a significant portion of the overall computational time, reducing the efficiency of multidisciplinary coupled simulation analysis. It is also necessary to simultaneously communicate the basic information and coupled data values ​​of each coupling interface to the other disciplinary solvers, allowing them to perform data processing techniques such as interpolation and linear weighting based on the coupled data basic information.

[0050] To this end, the present invention provides a multidisciplinary coupling interface scalable communication data encoding and decoding scheme, which can reduce the call of 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 multidisciplinary 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 the 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.

[0053] For the coupled launch / structure simulation analysis of an aircraft, a case folder for the aircraft's aerodynamic / structural coupling analysis is generated. This folder contains folders for the fluid solver, structural solver, coupled data interpolation module, and configuration files for the coupling framework. The solver folder contains the model mesh and configuration files, respectively, while the coupled interpolation folder contains the interpolation algorithm parameter settings. For the fluid solver setup, first use mesh generation software to generate the vehicle's external flow field mesh, set all wing surfaces as coupled interfaces, and output the flow field mesh in Plot3D format. In the fluid solver configuration file, 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)), and solution method (discrete format, turbulence model, etc.). Then, set the structural solver configuration file, input structural physical properties using structural preprocessing software, set all wings as coupled interfaces, and use modal analysis 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, set 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 number of physical quantities are first determined. According to the type of multidisciplinary coupling analysis problem and the coupling iterative mechanism between disciplines, the types and number of physical quantities of all data involved in the coupling communication in the multidisciplinary coupling analysis are determined. pFor example, in fluid-solid coupling analysis, the fluid needs to transfer the calculated force on the coupling interface to the solid to facilitate the static or dynamic analysis of the solid, and the solid needs to transfer the calculated displacement after deformation to the fluid to facilitate the fluid to update its shape and calculate the force under the new state. Therefore, the types and numbers 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 an aircraft wing, after the fluid solver solves a single coupling time step, the node aerodynamic force on the coupling interface is transferred to the structural solver for dynamic analysis; after the structural solver solves a single coupling time step, the displacement on the coupling interface is transferred to the fluid solver for fluid mesh deformation. Therefore, the physical quantity types 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 physical quantity types of the structural solver participating in the coupling communication include: node coordinates, node displacements, and setting the number of physical quantities Equal to 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 involving boundary surface coupling, such as aircraft aerodynamic-structural coupling analysis, there are typically multiple wing-like coupling interfaces, each with a large number of mesh nodes. In the mesh file, all boundaries are traversed, and the number of coupling interfaces is calculated based on the boundary type or boundary name set for the coupling interface during mesh generation. The number of nodes on each coupling interface is also calculated. When encoding and decoding coupled data, the number of coupling interfaces and the number of nodes on each coupling interface are used to determine the data quantity for a single physical quantity.

[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, representing 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 boundary node numbers of the grid block where the coupling interface is located. , No. The number of nodes of a wing is recorded as The structural solver uses the modal method to solve the structural control equations. The wing structure grid 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 and the components of the node aerodynamic force of each node of the shell element 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 node 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. Therefore, 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 discipline solver needs to communicate multiple types of data with each discipline solver separately. In order to achieve the communication of multiple physical quantity types of data to each solver based on the same coding format, without forcibly specifying the coding order of the 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 discipline 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-structure 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 physical quantity type symbols are assigned: {b, h, j}, and the physical quantity type association symbol is 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 within the mesh cell face 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 scalar, the data type symbol assigned to the physical quantity is 0. It is understood that in multidisciplinary coupling analysis, various disciplinary solvers use a variety of grid space discretization methods and physical quantity types, resulting in differences in the storage of physical quantities within the coupling interface of the coupled boundary mesh cells. Furthermore, when interpolating coupled data, physical quantities at different locations will use interpolation matrices constructed using different types of base points. Therefore, the encoded data must include information about the storage location of the physical quantity. Although the storage of physical quantity data within cells varies significantly, 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 within the cell face, such as the cell face's centroid or the cell face's Lagrange interpolation point; node storage refers to the storage of physical quantity data at all nodes within the cell face. When encoding coupled data, each type of physical quantity is assigned a symbol representing the storage location 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 indicates that the physical quantity is stored at the center of the face, and 0 indicates that it is stored at a 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. Therefore, the data storage location symbols 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: in 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 temperatures are 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 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 is 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, based on the obtained coupling interface, a 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 the order of the number of the coupling interface in the grid block from small to large, 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 the associated symbol, that is, the coding order of the coupling interface is obtained.

[0071] Specifically, since both the fluid solver and the structural solver 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, a target order of physical quantity types on the current coupling interface is obtained based on the order of physical quantity type-association symbols in the preset physical quantity type association table, and the target order is used as the coupling data encoding order for each physical quantity type on the current coupling interface. It is understood that the order of coupling data encoding for different physical quantity types on a single coupling interface is established based on the top-to-bottom order of the physical quantity types in Table 1. When encoding the physical quantities of a single coupling interface, the physical quantities of that coupling interface are encoded sequentially according to the established physical quantity encoding order. For example, when performing fluid-structure thermal coupling analysis on a typical aircraft, all wings of the aircraft are considered coupling interfaces, including the main wing, horizontal tail, and vertical tail. The physical quantity types involved in the coupling analysis in the fluid solver include node numbers, node coordinates, face-centered coordinates, node aerodynamic forces, node displacements, node temperatures, and face-centered heat flux. Therefore, the established physical quantity encoding order 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 forces; the physical quantity coding order in the structural solver is established as: node coordinates, node displacements.

[0074] In this embodiment, the order of mesh node coding in the current coupling interface is established as follows: In multidisciplinary coupling analysis, a single coupling interface typically contains coupling data of multiple physical quantity types on the same node set. To ensure that different types of coupling data correspond to the same node order when encoding or decoding, the order of nodes or face centers in a single coupling interface is established. For example, in a three-dimensional structural topology grid, nodes within a single grid block are numbered using the triplet format (I, J, K). Node numbers within the coupling interface in the I direction have the same I value, and the node coding order is sorted from smallest to largest based on the J and K values, respectively. The node coding order for the coupling interface in the J and K directions is the same as that for the I direction. The element face where the face center resides has four nodes. The face center coding order is based on the node with the smallest J and K values ​​among the four nodes. The reference nodes of all face centers form a reference node set. The face center coding order is based on the order of the nodes in the reference node set in the node coding order. In an 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 element face centers of the coupling interface is coded in ascending order according to the global number of the element face in which they are located.

[0075] Specifically, all physical quantities involved in the coupling analysis between the fluid solver and the structural solver are stored on the nodes of the unit surface, so only the coding order of a single coupling interface node needs to be established. The flow field grid file format is the structured grid plot3D format, and the order of the nodes is based on 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, all grid blocks and their surfaces are traversed to collect the grid blocks and grid surfaces containing the coupling interface boundary conditions. Then, the boundary grid blocks are sorted in order from small to large by the grid block number, and the grid surfaces are sorted in order by the grid block number. The order of the nodes in the grid surface is based on the number value of the node in both directions of the grid surface from small to large. For example, on the grid surface in the I direction, the nodes are first fixed by the value of J, then by the value of K from small to large, and then the value of J is increased. The order of the structural solver nodes is based on the order of the nodes in the finite element model file.

[0076] In this embodiment, in a multidisciplinary coupled 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 specific physical quantity data value. For example, when conducting an aerodynamic / structural coupled analysis of a wing based on a first-order loosely coupled strategy, after the fluid solver advances a coupling time step, it traverses all unit surfaces on the wing surface, performs an area integration of the pressure at the center of the surface, and obtains the unit surface aerodynamic force. Data processing methods such as linear weighting of the unit surface aerodynamic force are assigned to each node of the unit surface to obtain the aerodynamic force data values ​​of all nodes on the wing surface, i.e., the coupled data values.

[0077] Specifically, the fluid solver uses the unsteady flow method to solve, advancing the flow field by one time step in physical time. The pressure area integration is performed on the coupling interface unit surface to obtain the aerodynamic force of 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 node 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 representing the coupling data encoding of 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 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 the order of the established nodes / face centers. 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 at 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 one byte of memory space and store the data storage location symbol and data type symbol in the first two bits. Establish a node encoding order based on the ascending order of the number of all nodes on the main wing in the coordinate sequence of the grid file. And fill the node coordinate components into the memory in sequence according to the established node order, with the coordinates of each node 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 steps five and six; then, according to the coupling interface node coding sequence established in step ten, all nodes of the coupling interface are traversed, and the node coordinates are respectively encoded after the data storage location symbol and data type symbol according to the established node coding sequence. The coordinate components of each node are encoded in the order of x, y, and z. Similarly, all nodes on the coupling interface are traversed, and the node aerodynamic force or node displacement is encoded after the coordinate data according to the established node coding sequence. 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 that represents 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 , the number of vector physical quantity types , all physical quantity data values. First, the physical quantity Encode into the data sequence; then, according to the established coupling data encoding sequence, fill the physical quantity type associated symbol sequence on the coupling interface into the data sequence in turn, 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 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}, the number of scalar physical quantity types , the number of vector physical quantity types , nodal coordinate coupling data, face-centered coordinate coupling data, nodal aerodynamic coupling data, nodal temperature coupling data, face-centered 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 the structural solver converts the physical quantity Encoded into the coupling data sequence; then, the associated symbol sequence of the node coordinates, node aerodynamic force or node displacement established in step 4 is filled into the coupling data sequence, and finally, the coupling data completed by the established single physical type encoding is encoded into the coupling data sequence according to the physical quantity encoding order, 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 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 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 physical field dimension data obtained is encoded into the coupling data; then, the number of coupling interfaces obtained is Encode into the coupling data. Encode the established coupling interface symbol sequence into the coupling data, and encode the completed coupling data of the 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 content of the solver coding format for encoding the fluid solver coupling data is: 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, for the calculated physical dimension data , the number of coupling interfaces calculated , encode the basic information of coupling data such as the established coupling interface symbol sequence {0, 1, 2}; then encode the established coupling data of each wing, and fill the coupling data of the main wing, horizontal tail, and vertical tail in turn after the basic information of coupling data. At this point, the fluid solver and the structural solver have established the coupled data of each wing after the coupling interface encoding. Figure 2 As shown, the first is the encoding format of a single physical type coupling data, specifically, the node coordinate data: 、 、 、 、 Here and Labels used to identify 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 values ​​on different coordinate axes are then encoded in a single interface coupling data format. The single interface coupling data encoding format then transitions from the aforementioned single physical type coupling data encoding format to a single interface coupling data encoding format, and finally to a fluid solver coupling data encoding format. 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 structure surface.

[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 number of all coupling data participating in the coupling communication in the multidisciplinary coupling simulation analysis problem, and calculating 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, 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 encoding result, 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, 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. As can be seen, the process of assigning type association symbols and storage location symbols enables the processing of coupled data of various physical quantity types and different data types (scalar and vector), meeting the needs of complex multidisciplinary coupled analysis scenarios such as fluid-structure coupling and fluid-structure thermal coupling. This demonstrates the ability to handle multiple data types and enhances flexibility. By establishing a coding order for the coupling interface, the coupling data, and the grid node coding order, the method further considers different grid structures (structural topological grids and unstructural topological grids) and the storage location of physical quantities in the coupling interface grid element faces (nodes or face centers). By establishing a coding order based on these differences, the method can adapt to the different grid space discretization methods and physical quantity storage methods of various disciplines solvers, thus improving adaptability.Finally, based on this information and using a preset encoding method, the basic information of the physical quantities and coupling information of various types of physical quantities on each coupling interface is uniformly encoded. This allows data communication across multiple coupling interfaces to be completed with a single message-passing interface call. This significantly reduces the number of message-passing interface calls and the time cost associated with multiple interface calls. This effectively avoids the problem of communication time dominating the overall computational time, especially when there are many computational processes, thereby improving the efficiency of multidisciplinary coupled simulation analysis.

[0089] Reference Figure 3 As shown, the embodiment of the present invention discloses a multidisciplinary coupling surface scalable communication data decoding method. Compared with the previous embodiment, this embodiment further explains 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, based on the communication mapping relationship between the communication domains of the subject solvers established through multidisciplinary coupling analysis, the fluid solver obtains the target encoding result after encoding by the structure encoder, and the structure solver obtains the coupling data after encoding by the fluid encoder. The subject solvers receive the encoded coupling data, which is also the target encoding result.

[0092] Step S22: performing initial decoding on the target coding result to sequentially obtain the physical field dimension of the coupling data, the number of coupling interfaces, and the coding order of the coupling interfaces.

[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, namely the coupling interface symbol sequence. Then loop times, The cycle is based on the corresponding relationship 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 coded coupling data transmitted by the fluid solver, it begins 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, the flow field coupling interface symbol sequence {0, 1, 2} is obtained from the coupling data; finally, the coupling interfaces are decoded one by one according to the obtained coupling interface symbol sequence. Similarly, after the fluid solver receives the coded coupling data transmitted by the structural solver, it starts to decode the coupling data of the structural solver. Similarly, the physical field dimension data of the structural solver is obtained first. 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 on 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 association symbol sequence, and determine the third coding result corresponding to each type association symbol in the type association 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 physical quantity of the coupling interface. ; Then obtain the physical quantity type symbol sequence and convert the symbol sequence into the physical quantity type sequence according to Table 1. The last loop times, The next loop gets the Coupling data of different physical quantity types.

[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, and each time a coupling data sequence of a physical quantity type on the coupling interface is obtained.

[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 location 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 points equal ; If the symbol is 1, the total number of coupling interface face centers is obtained from the coupling data , and set the points equal Then get the data type symbol of the coupled data. If the symbol is 0, the data component of each point If the symbol is 1, the data component of each point for Calculate the number of physical quantity data . Then obtain from the coupled data data, according to Each data corresponds to a point, the data is mapped to each node or surface center, and the physical quantity data is assigned to the physical quantity of the subject solver.

[0100] Specifically, the decoding operation of a physical quantity type coupling data 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 of each point for Calculate the total amount of data for this physical quantity . Then obtain from the coupled 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 coupling data set is assigned to the physical quantity of the current discipline solver.

[0101] In this embodiment, the fluid solver uses the fluid and structural node coordinates as the basis points and establishes a displacement interpolation matrix G based on radial basis functions (RBFs). The structural node displacements are taken as known quantities and matrix-vector products are performed with the displacement interpolation matrix to obtain the displacements of the fluid mesh wing surface nodes. Similarly, the structural solver uses the fluid and structural node coordinates as the basis points and establishes an aerodynamic force interpolation matrix F based on the principle of virtual work. The fluid node aerodynamic forces are taken as known quantities and matrix-vector products are performed with the aerodynamic force interpolation matrix to obtain the aerodynamic forces of the structural mesh wing surface nodes.

[0102] After obtaining the nodal displacements of the fluid mesh wing surface, mesh deformation methods such as RBF-TFI are used to calculate the internal nodal coordinates of the mesh. The RBF-TFI method follows these steps: First, using the vertices of the mesh blocks on the coupling interface as control points and the vertices of the internal mesh blocks as interpolation base points, the RBF method is used to calculate the vertex displacements of the internal mesh blocks based on the mesh block displacements at the coupling interface. Second, using the vertices of all mesh blocks as control points, the RBF method is used to accurately calculate the nodal displacements at the edges of each mesh block. Then, based on the nodal displacements of the mesh block edges, a two-dimensional TFI method based on arc-length coordinates is used to quickly calculate the nodal displacements of each mesh block surface. Finally, based on the nodal displacements of the mesh block surfaces, a three-dimensional TFI method based on arc-length coordinates is used to quickly calculate the displacements of the internal points of the mesh block. Once the displacements of all nodes are calculated, they are added to the current node coordinates to obtain the new node coordinates. The mesh's geometric information (such as node coordinates, center of mass coordinates, volume, area, and area vector) and the Jacobi transform coefficient matrix are then recalculated. In this way, mesh deformation methods such as RBF-TFI can accurately calculate the internal node coordinates of the mesh based on information such as node displacements obtained after decoding, and recalculate the mesh's geometric information (such as node coordinates, center of mass coordinates, volume, area, area vector) and Jacobi transform coefficient matrix. This provides an effective solution for dealing with mesh deformation issues involved in complex multidisciplinary coupled analyses, better adapting to mesh morphological changes caused by coupling, and ensuring 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 one coupling time step, and the coupling iteration step number automatically increases by 1. The process then jumps to executing the fluid and structure solver time advance until the maximum number of coupling analysis steps is reached. In this way, in each coupling iteration, after accurately decoding the required data and 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, the physical time can be accumulated by one coupling time step, the coupling iteration step number automatically increases by 1, and the process jumps to executing the fluid and structure solver time advance 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] This shows that initial decoding of the coupled data's encoding format accurately and sequentially captures key information such as the coupling data's physical field dimensions, the number of coupling interfaces, and the encoding order of the coupling interfaces. Subsequent steps further target decoding of each portion of the encoding results, thereby completely and accurately restoring all types of information contained in the original encoding, including the number of physical quantities, the sequence of type-associated symbols, storage location symbols, data type symbols, and information about each grid node. This ensures that the encoded data can be accurately restored during the decoding process, providing an accurate data foundation for multidisciplinary coupling analysis.

[0105] Reference Figure 5 As shown, the present invention also correspondingly discloses a multidisciplinary coupling surface scalable communication data encoding device, comprising:

[0106] An information determination module 11 is used to determine the physical quantity types and physical quantity quantities of all coupling data involved in the 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;

[0107] an information allocation module 12, configured 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 cell 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 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 configured 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 representing the coupling data encoding of a single physical quantity type on the single coupling interface according to the order of the first encoding result and the second encoding result;

[0109] a second encoding module 14, configured 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, 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;

[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 physical quantity number 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 then encoding the grid node coding order according to the grid node coding result. All grid node information of a single coupling interface is traversed and encoded in node encoding order to obtain a corresponding second encoding result, and a third encoding result of coupling data encoding representing a single physical quantity type on a single coupling interface is constructed according to the order of the first encoding result and the second encoding result; the number of physical quantities is encoded into a coupling data sequence, and the type association symbol is filled into the encoded coupling data sequence, and then the third encoding result is encoded 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; basic coupling data information including the physical field dimension, the number of coupling interfaces, and the encoding order of the coupling interfaces is encoded to obtain a basic information encoding result, and the fourth encoding result is filled after the basic information encoding result to obtain a target encoding result. As can be seen, the process of assigning type association symbols and storage location symbols enables the processing of coupled data of various physical quantity types and different data types (scalar and vector), meeting the needs of complex multidisciplinary coupled analysis scenarios such as fluid-structure coupling and fluid-structure thermal coupling. This demonstrates the ability to handle multiple data types and enhances flexibility. By establishing a coding order for coupling interfaces, coupling data, and mesh node coding, the method further considers different mesh structures (structural and unstructural topological meshes) and the storage locations of physical quantities within the mesh element faces of the coupling interfaces (nodes or face centers). By addressing these differences, the coding order is established, making the method adaptable to the different mesh space discretization methods and physical quantity storage methods of various disciplinary solvers, thus improving adaptability. Finally, based on this information and using a preset coding method, the basic physical quantity and coupling information of various physical quantity types across each coupling interface are uniformly encoded, allowing data communication across multiple coupling interfaces to be completed through a single message passing interface call.The number of calls to the message passing interface is greatly reduced, reducing the time cost increase 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 embodiment of 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 should not be considered as any limitation to the scope of application of the present application.

[0113] Figure 6 This is 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 of the method disclosed in any of the aforementioned embodiments. Furthermore, 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 operating 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. The communication protocol it follows is any communication protocol that can be applied to the technical solution of this 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. 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 resource storage, can be a read-only memory, random access memory, disk or CD, 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] The operating system 221 is used to manage and control the hardware devices and computer programs 222 on the electronic device 20, enabling the processor 21 to calculate and process the massive amount of data 223 in the memory 22. It can be Windows Server, NetWare, Unix, Linux, etc. In addition to including computer programs capable of performing the methods performed by the electronic device 20 as disclosed in any of the aforementioned embodiments, the computer programs 222 may further include computer programs capable of performing other specific tasks. Data 223 may include data received by the electronic device from external devices, as well as data collected by its own input / output interface 25.

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

[0119] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from the other embodiments. Reference can be made to the descriptions of the identical or similar parts between the various embodiments. For the devices disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the descriptions are relatively simple, and the relevant parts can be referred to the descriptions of the methods.

[0120] Professionals may further appreciate that the units and algorithmic steps of each example described in conjunction with the embodiments disclosed herein can be implemented in 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. Professionals and technicians can 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 using 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 (Electric Erasable Programmable Read Only Memory), registers, hard disk, removable disk, CD-ROM (Compact Disc-Read Only Memory), or any other form of storage medium known in the technical field.

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

[0122] The above is a detailed introduction to the solution provided by the present invention. Specific examples are used herein to illustrate the principles and implementation methods 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 ideas. At the same time, for those skilled in the art, according to the ideas of the present invention, there may be changes in the specific implementation methods and application scopes. In summary, the contents of this specification should not be understood as limiting the present invention.

Claims

1. A method for scalable communication data encoding of coupling surfaces in wing aerodynamic / structural coupling simulation analysis, characterized in that: include: Determine the physical quantity type and physical quantity quantity of all coupling data involved in coupling communication in the aircraft wing aerodynamic / structural coupling simulation analysis problem, calculate the number of coupling interfaces involved in the aircraft wing aerodynamic / structural coupling simulation analysis and the number of grid nodes on each coupling interface, and determine the physical field dimension of the coupling data; wherein the coupling interface is the wing surface of the aircraft model; 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 method for scalable communication data encoding of coupling surfaces in wing aerodynamic / structural coupling simulation analysis according to claim 1, characterized in that: The calculation of the number of coupling interfaces involved in the wing aerodynamic / structural coupling simulation analysis of the aircraft and the number of grid nodes on each coupling interface includes: According to the grid information of the subject solver in solving the wing aerodynamic / structural coupling simulation analysis problem of the aircraft, all boundaries of the corresponding current grid are traversed, and the number of coupling interfaces participating in the wing aerodynamic / structural coupling simulation analysis of the aircraft 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 method for scalable communication data encoding of coupling surfaces in wing aerodynamic / structural coupling simulation analysis 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 method for encoding coupling surface scalable communication data in wing aerodynamic / structural coupling simulation analysis according to claim 3, 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 method for scalable communication data encoding of coupling surfaces in wing aerodynamic / structural coupling simulation analysis 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 method for scalable communication data encoding of coupling surfaces in wing aerodynamic / structural coupling simulation analysis 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 method for encoding coupling surface scalable communication data in wing aerodynamic / structural coupling simulation analysis 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 method for decoding coupling surface scalable communication data in wing aerodynamic / structural coupling simulation analysis, characterized in that: include: Obtaining a target encoding result encoded by a subject solver using a coupling surface scalable communication data encoding method in a wing aerodynamic / structural coupling simulation analysis as claimed in any one of claims 1 to 7; Initially decoding the target coding result to sequentially obtain the physical field dimension of the coupling data, the number of coupling interfaces, and the coding order of the coupling interfaces; the coupling interface is the wing surface of the aircraft model; 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 coupling surface expandable communication data encoding device in wing aerodynamic / structural coupling simulation analysis, 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 the aircraft wing aerodynamic / structural coupling simulation analysis problem, and calculate the number of coupling interfaces involved in the aircraft wing aerodynamic / structural coupling simulation analysis and the number of grid nodes on each coupling interface to determine the physical field dimension of the coupling data; wherein the coupling interface is the wing surface of the aircraft model; 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