CN117235832A - Method, device, equipment and medium for selecting aeroelastic coupling simulation interface point - Google Patents

Method, device, equipment and medium for selecting aeroelastic coupling simulation interface point Download PDF

Info

Publication number
CN117235832A
CN117235832A CN202311525648.6A CN202311525648A CN117235832A CN 117235832 A CN117235832 A CN 117235832A CN 202311525648 A CN202311525648 A CN 202311525648A CN 117235832 A CN117235832 A CN 117235832A
Authority
CN
China
Prior art keywords
point
grid
array
points
grid point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202311525648.6A
Other languages
Chinese (zh)
Other versions
CN117235832B (en
Inventor
李豪
江盟
孙岩
王昊
孟德虹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202311525648.6A priority Critical patent/CN117235832B/en
Publication of CN117235832A publication Critical patent/CN117235832A/en
Application granted granted Critical
Publication of CN117235832B publication Critical patent/CN117235832B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The application discloses a method, a device, equipment and a medium for selecting a pneumatic elastic coupling simulation interface point, which relate to the field of pneumatic elastic coupling numerical simulation and comprise the following steps: acquiring grid model file information and creating a dynamic array containing all grid point arrays and basic point grid point arrays; selecting a preset number of basic point grid points from all grid points based on a preset basic point selection method, and filling the basic point grid points into a basic point grid point array to obtain an initial basic point grid point array; determining a target grid point based on a preset base point grid point determining method, filling the target grid point into an initial base point grid point array to obtain a target base point grid point array; and judging whether the point selection ending rule is met based on the target base point grid point array, and ending the point selection operation if the point selection ending rule is met. According to the application, the grid points are interpolated, the grid points meeting the preset conditions are selected and added into the basic point grid point array, and the greedy algorithm is utilized to continuously select points and update the points, so that the simplification of the data points of the coupling interface is realized.

Description

Method, device, equipment and medium for selecting aeroelastic coupling simulation interface point
Technical Field
The invention relates to the field of aeroelastic coupling numerical simulation, in particular to the field of coupling data transmission in the aeroelastic coupling numerical simulation, and particularly relates to a method, a device, equipment and a medium for selecting an interface point of the aeroelastic coupling simulation.
Background
In recent years, with the rapid development of computational methods, computational fluid dynamics CFD (Computational Fluid Dynamics) is widely used in various complex flow numerical simulations, and plays an increasingly significant role in scientific research and product performance analysis. The pneumatic/structural coupling numerical simulation technology based on Computational Fluid Dynamics (CFD) and computational structural mechanics (CSD, computational structural dynamics) is also gradually applied to the design of the aerospace craft, and plays a key role in the performance analysis, optimization and evaluation processes.
In the aerodynamic/structural coupling numerical simulation based on CFD/CSD, CFD calculates the flow-around flow field to obtain aerodynamic load of the surface of the aircraft at first, convert the aerodynamic load data into external load input by the structure through the data conversion of the coupling interface and transmit to CSD, then CSD calculates the response of the structure based on the input external load, then convert the displacement, speed, etc. of the structure into displacement and speed of the flow field boundary and transmit to CFD through the data conversion of the coupling interface, CFD calculates grid and boundary condition according to the boundary displacement, speed of input update flow field, calculate the flow-around flow field in the new state, iterate so as to until coupling is stable or converged. The coupling interface data conversion defines the transmission relation of load, displacement and other data between CFD and CSD, and is an important factor affecting the accuracy and calculation efficiency of pneumatic/structural coupling simulation. In the pneumatic/structural coupling numerical simulation, the coupling interface data conversion is described by a matrix, and the conversion matrix is defined and calculated by the interface grid point coordinates of the CSD according to the principle of conservation of virtual work. In order to improve the efficiency of coupling data transmission, a limited number of structural finite element model interface grid points are generally selected during engineering application so as to reduce the scale of a coupling interface conversion matrix.
At present, the selection of the grid points of the interface of the finite element model of the structure is still realized by using a manual selection mode, so that time and labor are wasted, and when the grid points are large in scale, the internal grid points are easy to select, so that the point selection error is caused.
Disclosure of Invention
In view of the above, the application aims to provide a method, a device, equipment and a medium for selecting a pneumatic elastic coupling simulation interface point, which can realize the simplification of coupling interface data points, reduce the workload of manual processing and remarkably improve the overall efficiency of pneumatic/structural coupling analysis. The specific scheme is as follows:
in a first aspect, the application discloses a method for selecting a pneumatic elastic coupling simulation interface point, which comprises the following steps:
acquiring grid model file information, and creating a dynamic array based on the grid model file information; the dynamic array comprises all grid point arrays and basic point grid point arrays;
selecting a preset number of basic point grid points from all grid points in the grid model file information based on a preset basic point selection method, and filling the basic point grid points into the basic point grid point array to obtain an initial basic point grid point array;
determining a target grid point from all the grid points based on a preset base point grid point determining method, and filling the target grid point into the initial base point grid point array to obtain a target base point grid point array;
And judging whether a point selection ending rule is met based on the target base point grid point array, and ending the point selection operation if the point selection ending rule is met.
Optionally, before the obtaining the mesh model file information, the method further includes:
establishing a data structure of grid points; the data structure includes a grid point number, a grid point true coordinate, a grid point prediction coordinate, an error between the grid point true coordinate and the grid point prediction coordinate, a residual.
Optionally, the grid model file information is obtained, and a dynamic array is created based on the grid model file information; the dynamic array includes all grid point arrays and a base point grid point array, including:
acquiring grid model file information, and acquiring the total number of grid points in the grid model file information;
creating the all grid point array based on the data structure and all grid points in the grid model file information;
creating the basic point grid point array based on the data structure and a preset maximum reduced point;
initializing the whole grid point array and the basic point grid point array.
Optionally, before the determining method of the grid point based on the preset base point determines one target grid point from all the grid points, the method further includes:
Acquiring all grid points and corresponding real coordinates of the grid points, and calculating the grid point prediction coordinates corresponding to each grid point by using a preset difference method;
and calculating a difference value between the grid point predicted coordinates of each grid point and the grid point real coordinates to obtain the error.
Optionally, the determining method based on the preset base point and grid points determines a target grid point from the all grid points, and fills the target grid point into the initial base point and grid point array to obtain a target base point and grid point array, including:
the errors corresponding to all grid points are circularly traversed, and the grid point corresponding to the largest error is determined to be a target grid point;
determining the maximum value of the error as the current residual error;
and filling the target grid point, the corresponding target grid point prediction coordinates and the target grid point real coordinates into the initial basic point grid point array, and updating the residual value of the initial basic point grid point array by utilizing the current residual to obtain the target basic point grid point array.
Optionally, the determining whether the setpoint ending rule is satisfied based on the target base point grid point array includes:
Judging whether the current residual is smaller than a preset minimum residual and/or judging whether the grid point number of the target base point grid point array reaches the preset maximum reduced point number;
if the current residual is smaller than the preset minimum residual and/or if the number of grid points of the target base point grid point array reaches the preset maximum reduced point number, judging that the point selection ending rule is met;
if the current residual is greater than or equal to the preset minimum residual and/or if the number of grid points of the target base point grid point array does not reach the preset maximum reduced point number, judging that the point selection ending rule is not met;
and when the point selection ending rule is not met, re-entering the step of circularly traversing the errors corresponding to all grid points.
Optionally, the selecting a preset number of base point grid points from all grid points in the grid model file information based on a preset base point selection method, and filling the base point grid points into the base point grid point array to obtain an initial base point grid point array, including:
selecting a preset number of basic point grid points from all grid points in the grid model file information through a random function;
Acquiring real coordinates of the base point grid points corresponding to the base point grid points and the base point grid point numbers;
and filling the base point grid points, the base point grid point real coordinates and the base point grid point numbers into the base point grid point array to obtain an initial base point grid point array.
In a second aspect, the application discloses a device for selecting a pneumatic elastic coupling simulation interface point, comprising:
the dynamic array creation module is used for acquiring the grid model file information and creating a dynamic array based on the grid model file information; the dynamic array comprises all grid point arrays and basic point grid point arrays;
the base point and grid point selection module is used for selecting a preset number of base point and grid points from all grid points in the grid model file information based on a preset base point selection method, and filling the base point and grid points into the base point and grid point array to obtain an initial base point and grid point array;
the target grid point determining module is used for determining one target grid point from all the grid points based on a preset base point grid point determining method, and filling the target grid point into the initial base point grid point array to obtain a target base point grid point array;
And the ending point selecting module is used for judging whether the point selecting ending rule is met based on the target base point grid point array, and ending the point selecting operation if the point selecting ending rule is met.
In a third aspect, the present application discloses an electronic device, comprising:
a memory for storing a computer program;
a processor for executing the computer program to implement the steps of the aeroelastic coupling simulation interface point selection method as previously disclosed.
In a fourth aspect, the present application discloses a computer-readable storage medium for storing a computer program; wherein the computer program when executed by a processor implements the aeroelastic coupling simulation interface point selection method as previously disclosed.
The application provides a method for selecting a pneumatic elastic coupling simulation interface point, which comprises the following steps: acquiring grid model file information, and creating a dynamic array based on the grid model file information; the dynamic array comprises all grid point arrays and basic point grid point arrays; selecting a preset number of basic point grid points from all grid points in the grid model file information based on a preset basic point selection method, and filling the basic point grid points into the basic point grid point array to obtain an initial basic point grid point array; determining a target grid point from all the grid points based on a preset base point grid point determining method, and filling the target grid point into the initial base point grid point array to obtain a target base point grid point array; and judging whether a point selection ending rule is met based on the target base point grid point array, and ending the point selection operation if the point selection ending rule is met. Therefore, the method and the device select the grid points meeting the preset conditions by interpolating the grid points, add the grid point array of the base points, continuously select the points by using a greedy algorithm and update the grid point array of the base points, realize the simplification of the data points of the coupling interface, reduce the workload of manual processing and obviously improve the overall efficiency of pneumatic/structural coupling analysis.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings that are required to be used in the embodiments or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are only embodiments of the present application, and that other drawings can be obtained according to the provided drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for selecting a pneumatic elastic coupling simulation interface point according to the present application;
FIG. 2 is a flow chart of a particular method for selecting a pneumatic elastic coupling analog interface point in accordance with the present application;
FIG. 3 is a schematic diagram of a real coordinate model of all points on the surface of a wing disclosed by the application;
FIG. 4 is a schematic diagram of a model of true positions and random two points of all points on the surface of a wing disclosed by the application;
FIG. 5 is a schematic diagram of the actual and predicted positions of local grid points of a wing in accordance with the present disclosure;
FIG. 6 is a schematic diagram of the true positions of all points on the surface of a wing and a selected three-point model according to the present disclosure;
FIG. 7 is a simplified model schematic diagram of a surface point of a wing in accordance with the present disclosure;
FIG. 8 is a schematic view of a model of all-point prediction of the wing surface of the present disclosure;
FIG. 9 is a schematic diagram of a pneumatic elastic coupling simulation interface point selection device according to the present application;
fig. 10 is a block diagram of an electronic device according to the present application.
Detailed Description
The following description of the embodiments of the present application will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present application, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
At present, the selection of the grid points of the interface of the finite element model of the structure is still realized by adopting a manual selection mode, so that time and labor are wasted, and when the grid points are large in scale, the internal grid points are easy to select, so that the point selection error is caused. Therefore, the application provides a pneumatic elastic coupling simulation interface point selection method, which can realize the simplification of coupling interface data points, reduce the workload of manual processing and remarkably improve the overall efficiency of pneumatic/structural coupling analysis.
The embodiment of the application discloses a method for selecting a pneumatic elastic coupling simulation interface point, which is shown in fig. 1 and comprises the following steps:
Step S11: and acquiring grid model file information, and creating a dynamic array based on the grid model file information.
In the embodiment, grid model file information is acquired, and a dynamic array is created based on the grid model file information; the dynamic array includes an array of all grid points and an array of base grid points. Specifically, grid model file information is obtained, and the total number of grid points in the grid model file information is obtained; creating the all grid point array based on the data structure and all grid points in the grid model file information; creating the basic point grid point array based on the data structure and a preset maximum reduced point; initializing the whole grid point array and the basic point grid point array. It can be understood that, before acquiring the mesh model file information, a data structure of the mesh points is established; the data structure includes a grid point number, a grid point true coordinate, a grid point prediction coordinate, an error between the grid point true coordinate and the grid point prediction coordinate, a residual.
It will be appreciated that the mesh point data structure is first established, with the mesh point data structure established as empty. The information of the grid point data structure includes a grid point number INDEX i Wherein i ranges from 1 to nv, nv representing the total number of surface grid points; grid point real coordinates (Xi, yi, zi), where i ranges from 1 to nv; grid point prediction coordinates (XXi, YYi, ZZi), where i ranges from 1 to nv; error between real coordinates and predicted coordinates of grid pointsDELTA i Wherein i ranges from 1 to nv; residual RESID.
Further, two dynamic arrays are created to obtain grid model file information. Firstly, obtaining the total point number of the grid by reading the file information of the grid surface modelnvAnd set the maximum reduced point numbernselectThe method comprises the steps of carrying out a first treatment on the surface of the Then, according to the established data structure, creating two dynamic arrays of grid point types, which are all grid point arrays respectivelyxyz[nv]Simple grid point array (i.e. basic point grid point array)select[nselect]Wherein, the method comprises the steps of, wherein,nvrepresenting the total number of grid points,nselectthe number is the preset maximum reduced point number; then initialize the entire grid point arrayxyz[nv]The coordinate information of grid points in the grid surface model file is processedX i ,Y i ,Z i ) Grid point number informationINDEX i Writing an array of all grid pointsxyz[nv]In, and initialize the reduced grid point arrayselect[nselect]。
Step S12: and selecting a preset number of base point grid points from all grid points in the grid model file information based on a preset base point selection method, and filling the base point grid points into the base point grid point array to obtain an initial base point grid point array.
In this embodiment, after mesh model file information is obtained and a dynamic array is created based on the mesh model file information, a preset number of base point mesh points are selected from all mesh points in the mesh model file information based on a preset base point selection method, and the base point mesh points are filled into the base point mesh point array to obtain an initial base point mesh point array. Specifically, selecting a preset number of base point grid points from all grid points in the grid model file information through a random function; acquiring real coordinates of the base point grid points corresponding to the base point grid points and the base point grid point numbers; and filling the base point grid points, the base point grid point real coordinates and the base point grid point numbers into the base point grid point array to obtain an initial base point grid point array.
It will be appreciated that one or more surface grid points are randomly selected as the basis points for the interpolation method, resulting in basis point grid points. The method comprises the following specific steps: randomly selecting k grid points from all surface grid points by a random function, wherein k represents one or more grid points; then adding the selected grid points into a simplified grid point array select [ nselect ], and writing coordinate information and number information of the corresponding grid points into the array select [ nselect ] to obtain an initial base point grid point array; the selected grid points are used as the base points of the interpolation method.
Step S13: and determining a target grid point from all the grid points based on a preset base point grid point determining method, and filling the target grid point into the initial base point grid point array to obtain a target base point grid point array.
In this embodiment, after the base point grid point is filled into the base point grid point array to obtain an initial base point grid point array, a target grid point is determined from all the grid points based on a preset base point grid point determining method, and the target grid point is filled into the initial base point grid point array to obtain a target base point grid point array. Specifically, the errors corresponding to all grid points are circularly traversed, and the grid point corresponding to the largest error is determined as a target grid point; determining the maximum value of the error as the current residual error; and filling the target grid point, the corresponding target grid point prediction coordinates and the target grid point real coordinates into the initial basic point grid point array, and updating the residual value of the initial basic point grid point array by utilizing the current residual to obtain the target basic point grid point array.
Step S14: and judging whether a point selection ending rule is met based on the target base point grid point array, and ending the point selection operation if the point selection ending rule is met.
In this embodiment, after the target grid point is filled into the initial base point grid point array to obtain the target base point grid point array, whether the point selection ending rule is satisfied is determined based on the target base point grid point array, and if the point selection ending rule is satisfied, the point selection operation is ended. It can be understood that, whether the current residual is smaller than a preset minimum residual and/or whether the number of grid points of the target base point grid point array reaches the preset maximum reduced point number is judged; if the current residual is smaller than the preset minimum residual and/or if the number of grid points of the target base point grid point array reaches the preset maximum reduced point number, judging that the point selection ending rule is met; if the current residual is greater than or equal to the preset minimum residual and/or if the number of grid points of the target base point grid point array does not reach the preset maximum reduced point number, judging that the point selection ending rule is not met; and when the point selection ending rule is not met, re-entering the step of circularly traversing the errors corresponding to all grid points.
It is to be understood that the end-of-setpoint determination operation is performed every time after one target grid point is determined from the all grid points based on a preset base point grid point determination method and the target grid point is filled into the initial base point grid point array. Until residual error RESIDLess than or equal to a preset valueRESID min Or the total number of selected points of the target base point grid point array reaches the preset reduced point numbernselectAnd finishing the point selection and obtaining a final reduced point set.
Further, a model file containing reduced point information, all point prediction information and all point real information is obtained. The method comprises the following specific steps: creating a model file, and writing coordinate information of the reduced points in the final reduced point set into the file; then writing all grid point prediction coordinate information into a file; and finally, writing real coordinate information of all grid points, and closing the file.
The invention provides a data point simplifying method of a structure finite element model coupling interface, namely, a mesh point with the largest coordinate error is selected to be added into a simplified point set by interpolating the coordinates of the mesh point of the surface, and the points are continuously selected by a greedy algorithm to obtain the final simplified point set, so that the data point of the coupling interface is simplified. The method comprises the following specific steps: reading grid point information from a structural finite element grid surface model file, creating all grid point arrays to record the grid point information, and creating a base point array to store final simplified point information; randomly selecting a plurality of grid points from all the grid points as basic points of an interpolation method, obtaining predicted coordinates of all the grid points through the interpolation method based on RBF (Radial Basis Function ), and selecting points with the largest error between the true coordinates and the predicted coordinates of all the grid points to be added into a basic point array; and (3) obtaining updated predicted coordinates of all grid points by interpolation of the updated base point set, selecting points with the largest error between the true coordinates and the predicted coordinates of all grid points again, adding the points into the updated base point array, repeatedly greedy selecting points until the error of the selected points is smaller than or equal to a preset value or the total number of the selected points reaches the preset maximum simplified point, ending the base point selection, and enabling the base point set at the moment to be the final simplified point set, thereby realizing the simplification of the data points of the coupling interface. The invention reduces the workload of manual processing and obviously improves the overall efficiency of pneumatic/structural coupling analysis.
The application provides a method for selecting a pneumatic elastic coupling simulation interface point, which comprises the following steps: acquiring grid model file information, and creating a dynamic array based on the grid model file information; the dynamic array comprises all grid point arrays and basic point grid point arrays; selecting a preset number of basic point grid points from all grid points in the grid model file information based on a preset basic point selection method, and filling the basic point grid points into the basic point grid point array to obtain an initial basic point grid point array; determining a target grid point from all the grid points based on a preset base point grid point determining method, and filling the target grid point into the initial base point grid point array to obtain a target base point grid point array; and judging whether a point selection ending rule is met based on the target base point grid point array, and ending the point selection operation if the point selection ending rule is met. Therefore, the method and the device select the grid points meeting the preset conditions by interpolating the grid points, add the grid point array of the base points, continuously select the points by using a greedy algorithm and update the grid point array of the base points, realize the simplification of the data points of the coupling interface, reduce the workload of manual processing and obviously improve the overall efficiency of pneumatic/structural coupling analysis.
Referring to fig. 2, an embodiment of the present invention discloses a method for selecting a pneumatic elastic coupling analog interface point, and compared with the previous embodiment, the present embodiment further describes and optimizes a technical solution.
Step S21: and acquiring grid model file information, and creating a dynamic array based on the grid model file information.
In this embodiment, in a specific real-time manner, for example, the test case is an airfoil surface model file, and the file includes coordinate information of an airfoil surface point and surface triangle unit information. The total number of the surface points of the model is 9793, and the number of the set reduced points is 500. First, a data structure of grid points is established, wherein the data structure comprises grid point numbersINDEX i WhereiniIn the range of 1 tonvnv=9793The method comprises the steps of carrying out a first treatment on the surface of the Grid point true value coordinates [ ]X i ,Y i ,Z i ) WhereiniIn the range of 1 tonvnv=9793; grid point predictive coordinates [ ]XX i ,YY i ,ZZ i ) WhereiniIn the range of 1 tonvnv=9793; error between real coordinates and predicted coordinates of grid pointsDELTA i WhereiniIn the range of 1 tonvnv=9793; residual isRESID. Then, obtaining the total number of grid points of the wing surface model by reading the file information of the wing surface modelnv=9793, and set the maximum reduced point numbernselect=500, as shown in fig. 3, 500 points need to be selected from the full-point model of the wing surface instead of the full-point model; then, according to the established data structure, creating two dynamic arrays of grid point types, which are all grid point arrays respectively xyz[nv]Simplified grid point arrayselect[nselect]Wherein, the method comprises the steps of, wherein,nv=9793 denotes the total number of grid points,nselect=500 is a preset maximum number of selected points; initializing an array of full grid pointsxyz[nv]Actual coordinate information of grid points in the grid surface model fileX i ,Y i ,Z i ) Grid point number informationINDEX i Writing an array of all grid pointsxyz[nv]In, and initialize the reduced grid point arrayselect[nselect]。
Step S22: and selecting a preset number of base point grid points from all grid points in the grid model file information based on a preset base point selection method, and filling the base point grid points into the base point grid point array to obtain an initial base point grid point array.
In this embodiment, in a specific real-time manner, 2 grid points are randomly selected from all the surface grid points by a random function; the selected grid point is then added to the reduced grid point arrayselect[nselect]Obtaining an initial basic point grid point array, and writing coordinate information and number information corresponding to the two grid points into the arrayselect[1]Andselect[2]in (a) and (b); finally, the two grid points are taken as the base points of the interpolation method. As shown in fig. 4, the larger "black sphere" represents the randomly selected 2 dots, and the other smaller "black dots" represent all other grid points.
Step S23: and acquiring all grid points and the corresponding real coordinates of the grid points, and calculating the corresponding predicted coordinates of the grid points by using a preset difference method.
In this embodiment, in a specific implementation, the grid point array is simplified fromselect[nselect]Two points are randomly selected from the (initial array of base point grid points)select[1]Andselect[2]as a base point of the interpolation method; then obtaining predicted coordinates of all grid points through an RBF method and an interpolation method which is not limited to the RBFXX i ,YY i ,ZZ i )。
Step S24: and calculating a difference value between the grid point predicted coordinates of each grid point and the grid point real coordinates to obtain the error.
In this embodiment, a difference between the grid point predicted coordinates and the grid point true coordinates of each of the grid points is calculated to obtain the error. It can be understood that the predicted coordinates of all grid points are obtainedXX i ,YY i ,ZZ i ) And calculate the error between the true coordinates and the predicted coordinates of the grid pointsDELTA i . By readingkInterpolation base points are obtained by utilizing interpolation method which is not limited to RBF, and all grid point prediction coordinates are obtainedXX i ,YY i ,ZZ i ) The method comprises the steps of carrying out a first treatment on the surface of the Through the true coordinates of each grid pointX i ,Y i ,Z i ) And forecast coordinates [ ]XX i ,YY i ,ZZ i ) Calculating the error of each grid pointDELTA i
In a specific embodiment, the true coordinates of each grid point are calculatedX i ,Y i ,Z i ) And forecast coordinates [ ]XX i ,YY i ,ZZ i ) Calculating the error of each grid pointDELTA i . The true position of the local grid point of the wing is different from the predicted position as shown in FIG. 5, wherein the "black sphere" represents the true position of the grid point The "circle" adjacent to the "black sphere" represents the predicted position of the grid point, which represents the difference in position when the point 500 is selected.
Step S25: and determining a target grid point from all the grid points based on a preset base point grid point determining method, and filling the target grid point into the initial base point grid point array to obtain a target base point grid point array.
In this embodiment, a target grid point is determined from all the grid points based on a preset base point grid point determining method, and the target grid point is filled into the initial base point grid point array to obtain a target base point grid point array. It can be appreciated that the reduced set of points is selected by a greedy algorithm to obtain the array of target base point grid points. The method comprises the following specific steps: cycling through the error of each grid point calculated in step fourDELTA i The grid points at this time are grid points not selected as the base points, and an error is foundDELTA i Maximum pointADots (dot)ATo be selected as the firstk+1 base points, residual errorRESIDEqual point(s)AError of (2), i.eRESID=MAX(DELTA i ) And grid points are combinedACoordinate information and residual of (a)RESIDAssignment toselect[k+1]The method comprises the steps of carrying out a first treatment on the surface of the Then, the selected materials are againk+1 points are used again as updated base point set RBFInterpolation method for obtaining updated predicted coordinates of all grid pointsXX i ,YY i ,ZZ i ) And traversing each grid point again to obtain updated errors of each grid pointDELTA i And find updated errorsDELTA i Maximum pointBAt this time, a dotBI.e. the selected firstk+2 points, residual errorRESIDError equal to point BDELTA i And spotBCoordinate information and residual of (a)RESIDAssignment toselect[k+2]。
In a specific embodiment, the calculated error of each grid point is circularly traversedDELTA i Find errorDELTA i Maximum pointADots (dot)AI.e. the 3 rd point selected, also the 1 st point selected by interpolation method, residualRESIDEqual to the grid pointsAError of (2), i.eRESID=MAX(DELTA i ) = 0.595055, as shown in fig. 6, the true positions of all points on the wing surface are calculated by selecting a three-point model, wherein the larger "black square" represents the 1 st point selected by the interpolation method, the larger "black sphere" represents the 2 points selected randomly, and the other smaller "black points" represent all other points. Handle pointACoordinate information and residual of (a)RESIDAssignment toselect[3]The method comprises the steps of carrying out a first treatment on the surface of the Then using the selected 3 points as updated base points, and using RBF interpolation method again to obtain updated predicted coordinates of all grid pointsXX i ,YY i ,ZZ i ) Traversing each point again to obtain updated error of each pointDELTA i And find updated distance error DELTA i Maximum pointBDots (dot)BI.e. the 4 th selected point, grid pointBDistance errorDELTA i The 2 nd point residual error selected as interpolation methodRESIDI.e.RESID=MAX(DELTA i ) = 0.329984, and dotBCoordinate information and residual of (a)RESIDAssignment toselect[4]The method comprises the steps of carrying out a first treatment on the surface of the Finally, the points are greedy selected continuously through the flow until the residual error of the selected pointsRESIDLess than or equal to a preset valueRESID min Or the total number of selected points reaches the preset reduced point numbernselect=500, in this example, the residual error of the 500 th point is selectedRESID=1.20348×10 -5 Is larger than a preset valueRESID min =10 -5 And finally, after the point selection is finished, obtaining a simplified point set.
Further, a model file is created, coordinate true value information of 500 reduced points is written into the file, all grid point prediction coordinate information is written into the file, all grid point real coordinate information is written into the file, and the file is closed, so that the model file containing the reduced points, all point prediction information and all point true value information is obtained. The wing surface point reduced model shown in fig. 7, 500 "black points" represent selected reduced points, and the wing surface full-point interpolation model shown in fig. 8, all "black points" represent interpolation points of all grid points.
Step S26: and judging whether a point selection ending rule is met based on the target base point grid point array, and ending the point selection operation if the point selection ending rule is met.
For the specific content of the above step S26, reference may be made to the corresponding content disclosed in the foregoing embodiment, and no further description is given here.
Therefore, the embodiment of the application creates the dynamic array by acquiring the grid model file information and based on the grid model file information; selecting a preset number of basic point grid points from all grid points in the grid model file information based on a preset basic point selection method, and filling the basic point grid points into the basic point grid point array to obtain an initial basic point grid point array; acquiring all grid points and corresponding real coordinates of the grid points, and calculating the grid point prediction coordinates corresponding to each grid point by using a preset difference method; calculating a difference between the grid point predicted coordinates and the grid point true coordinates of each grid point to obtain the error; determining a target grid point from all the grid points based on a preset base point grid point determining method, and filling the target grid point into the initial base point grid point array to obtain a target base point grid point array; and judging whether a point selection ending rule is met based on the target base point grid point array, ending the point selection operation if the point selection ending rule is met, simplifying data points of a coupling interface, reducing the workload of manual processing, and obviously improving the overall efficiency of pneumatic/structural coupling analysis.
Referring to fig. 9, the embodiment of the application also correspondingly discloses a device for selecting a pneumatic elastic coupling simulation interface point, which comprises:
the dynamic array creation module 11 is used for acquiring grid model file information and creating a dynamic array based on the grid model file information; the dynamic array comprises all grid point arrays and basic point grid point arrays;
a base point and grid point selection module 12, configured to select a preset number of base point and grid points from all grid points in the grid model file information based on a preset base point selection method, and fill the base point and grid points into the base point and grid point array to obtain an initial base point and grid point array;
a target grid point determining module 13, configured to determine a target grid point from all the grid points based on a preset base point grid point determining method, and fill the target grid point into the initial base point grid point array to obtain a target base point grid point array;
and the ending point selecting module 14 is configured to determine whether a point selecting ending rule is satisfied based on the target base point grid point array, and if the point selecting ending rule is satisfied, end the point selecting operation.
It can be seen that the present application includes: acquiring grid model file information, and creating a dynamic array based on the grid model file information; the dynamic array comprises all grid point arrays and basic point grid point arrays; selecting a preset number of basic point grid points from all grid points in the grid model file information based on a preset basic point selection method, and filling the basic point grid points into the basic point grid point array to obtain an initial basic point grid point array; determining a target grid point from all the grid points based on a preset base point grid point determining method, and filling the target grid point into the initial base point grid point array to obtain a target base point grid point array; and judging whether a point selection ending rule is met based on the target base point grid point array, and ending the point selection operation if the point selection ending rule is met. Therefore, the method and the device select the grid points meeting the preset conditions by interpolating the grid points, add the grid point array of the base points, continuously select the points by using a greedy algorithm and update the grid point array of the base points, realize the simplification of the data points of the coupling interface, reduce the workload of manual processing and obviously improve the overall efficiency of pneumatic/structural coupling analysis.
In some embodiments, the dynamic array creation module 11 specifically includes:
a data structure establishing unit for establishing a data structure of the grid points; the data structure comprises a grid point number, a grid point real coordinate, a grid point prediction coordinate, an error between the grid point real coordinate and the grid point prediction coordinate and a residual error;
the file information acquisition unit is used for acquiring the file information of the grid model;
a grid point total number acquisition unit, configured to acquire a grid point total number in the grid model file information;
a total grid point array creation unit configured to create the total grid point array based on the data structure, the total grid points in the grid model file information;
a base point grid point array creation unit, configured to create the base point grid point array based on the data structure and a preset maximum reduced point number;
and the array initializing unit is used for initializing the whole grid point array and the basic point grid point array.
In some specific embodiments, the base point and grid point selecting module 12 specifically includes:
a base point and grid point selecting unit for selecting a preset number of base point and grid points from all grid points in the grid model file information through a random function;
The real coordinate and number acquisition unit is used for acquiring real coordinates of the base point grid points corresponding to the base point grid points and the base point grid point numbers;
and the initial base point grid point array acquisition unit is used for filling the base point grid points, the base point grid point real coordinates and the base point grid point numbers into the base point grid point array to obtain an initial base point grid point array.
In some specific embodiments, the target grid point determining module 13 specifically includes:
a real coordinate acquiring unit, configured to acquire real coordinates of all grid points and the corresponding grid points;
a predicted coordinate calculation unit, configured to calculate the predicted coordinate of the grid point corresponding to each grid point by using a preset difference method;
an error calculation unit configured to calculate a difference between the grid point predicted coordinates and the grid point true coordinates of each of the grid points to obtain the error;
a first error circulation traversing unit, configured to circulate and traverse the errors corresponding to all grid points;
a target grid point determining unit configured to determine the grid point corresponding to the error that is the largest as a target grid point;
A current residual determining unit, configured to determine a value of the error that is the largest as a current residual;
an initial basic point grid point array filling unit, configured to fill the target grid point, the corresponding target grid point prediction coordinates, and the target grid point real coordinates into the initial basic point grid point array;
and the residual value updating unit is used for updating the residual value of the initial base point grid point array by utilizing the current residual so as to obtain the target base point grid point array.
In some embodiments, the ending setpoint module 14 specifically includes:
a current residual error judging unit, configured to judge whether the current residual error is smaller than a preset minimum residual error;
a grid point number judging unit, configured to judge whether the number of grid points of the target base point grid point array reaches the preset maximum reduced point number;
the first setpoint ending judging unit is used for judging that the setpoint ending rule is met if the current residual is smaller than the preset minimum residual and/or if the number of the grid points of the target base point grid point array reaches the preset maximum reduced point number;
a second setpoint ending determining unit, configured to determine that the setpoint ending rule is not satisfied if the current residual is greater than or equal to the preset minimum residual and/or if the number of the grid points of the target base point grid point array does not reach the preset maximum reduced number of points;
A second error cycle traversing unit, configured to reenter the step of cycle traversing the errors corresponding to all grid points when the setpoint ending rule is not satisfied;
and the ending unit is used for ending the point selecting operation if the point selecting ending rule is met.
Further, the embodiment of the application also provides electronic equipment. Fig. 10 is a block diagram of an electronic device 20, according to an exemplary embodiment, and the contents of the diagram should not be construed as limiting the scope of use of the present application in any way.
Fig. 10 is a schematic structural diagram of an electronic device 20 according to 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. Wherein the memory 22 is configured to store a computer program that is loaded and executed by the processor 21 to implement the relevant steps of the aeroelastic coupling analog interface point selection method disclosed in any of the foregoing embodiments. In addition, the electronic device 20 in the present embodiment may be specifically an electronic computer.
In this embodiment, the power supply 23 is configured to provide an 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 an external device, and the communication protocol to be followed is any communication protocol applicable to the technical solution of the present application, which is not specifically limited herein; the input/output interface 25 is used for acquiring external input data or outputting external output data, and the specific interface type thereof may be selected according to the specific application requirement, which is not limited herein.
The memory 22 may be a carrier for storing resources, such as a read-only memory, a random access memory, a magnetic disk, or an optical disk, and the resources stored thereon may include an operating system 221, a computer program 222, and the like, and the storage may be temporary storage or permanent storage.
The operating system 221 is used for managing and controlling various hardware devices on the electronic device 20 and computer programs 222, which may be Windows Server, netware, unix, linux, etc. The computer program 222 may further comprise a computer program capable of performing other specific tasks in addition to the computer program capable of performing the aeroelastic coupling analog interface point selection method performed by the electronic device 20 as disclosed in any of the previous embodiments.
Further, the embodiment of the application also discloses a computer readable storage medium, wherein the medium stores a computer program, and when the computer program is loaded and executed by a processor, the steps of the aeroelastic coupling simulation interface point selection method disclosed in any embodiment are realized.
In this specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, so that the same or similar parts between the embodiments are referred to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant points refer to the description of the method section.
Finally, it is further noted that relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
The method, the device, the equipment and the medium for selecting the aeroelastic coupling simulation interface point provided by the invention are described in detail, and specific examples are applied to the principle and the implementation mode of the invention, and the description of the examples is only used for helping to understand the method and the core idea of the invention; meanwhile, as those skilled in the art will have variations in the specific embodiments and application scope in accordance with the ideas of the present invention, the present description should not be construed as limiting the present invention in view of the above.

Claims (10)

1. A method for selecting a aeroelastic coupling analog interface point, comprising:
acquiring grid model file information, and creating a dynamic array based on the grid model file information; the dynamic array comprises all grid point arrays and basic point grid point arrays;
selecting a preset number of basic point grid points from all grid points in the grid model file information based on a preset basic point selection method, and filling the basic point grid points into the basic point grid point array to obtain an initial basic point grid point array;
determining a target grid point from all the grid points based on a preset base point grid point determining method, and filling the target grid point into the initial base point grid point array to obtain a target base point grid point array;
and judging whether a point selection ending rule is met based on the target base point grid point array, and ending the point selection operation if the point selection ending rule is met.
2. The method for selecting a aeroelastic coupling simulation interface point according to claim 1, wherein before the step of obtaining the mesh model file information, the method further comprises:
establishing a data structure of grid points; the data structure includes a grid point number, a grid point true coordinate, a grid point prediction coordinate, an error between the grid point true coordinate and the grid point prediction coordinate, a residual.
3. The method for selecting aeroelastic coupling simulation interface points according to claim 2, wherein the method comprises the steps of obtaining grid model file information and creating a dynamic array based on the grid model file information; the dynamic array includes all grid point arrays and a base point grid point array, including:
acquiring grid model file information, and acquiring the total number of grid points in the grid model file information;
creating the all grid point array based on the data structure and all grid points in the grid model file information;
creating the basic point grid point array based on the data structure and a preset maximum reduced point;
initializing the whole grid point array and the basic point grid point array.
4. The method of claim 3, wherein before determining a target grid point from the all grid points based on the preset base point grid point determination method, further comprising:
acquiring all grid points and corresponding real coordinates of the grid points, and calculating the grid point prediction coordinates corresponding to each grid point by using a preset difference method;
And calculating a difference value between the grid point predicted coordinates of each grid point and the grid point real coordinates to obtain the error.
5. The method of claim 4, wherein determining a target grid point from the total grid points based on a preset base point grid point determination method and filling the target grid point into the initial base point grid point array to obtain a target base point grid point array, comprises:
the errors corresponding to all grid points are circularly traversed, and the grid point corresponding to the largest error is determined to be a target grid point;
determining the maximum value of the error as the current residual error;
and filling the target grid point, the corresponding target grid point prediction coordinates and the target grid point real coordinates into the initial basic point grid point array, and updating the residual value of the initial basic point grid point array by utilizing the current residual to obtain the target basic point grid point array.
6. The method of claim 5, wherein determining whether the setpoint end rule is satisfied based on the target setpoint grid point array comprises:
Judging whether the current residual is smaller than a preset minimum residual and/or judging whether the grid point number of the target base point grid point array reaches the preset maximum reduced point number;
if the current residual is smaller than the preset minimum residual and/or if the number of grid points of the target base point grid point array reaches the preset maximum reduced point number, judging that the point selection ending rule is met;
if the current residual is greater than or equal to the preset minimum residual and/or if the number of grid points of the target base point grid point array does not reach the preset maximum reduced point number, judging that the point selection ending rule is not met;
and when the point selection ending rule is not met, re-entering the step of circularly traversing the errors corresponding to all grid points.
7. The method according to any one of claims 1 to 6, wherein the selecting a preset number of base point grid points from all grid points in the grid model file information based on a preset base point selection method, and filling the base point grid points into the base point grid point array to obtain an initial base point grid point array includes:
Selecting a preset number of basic point grid points from all grid points in the grid model file information through a random function;
acquiring real coordinates of the base point grid points corresponding to the base point grid points and the base point grid point numbers;
and filling the base point grid points, the base point grid point real coordinates and the base point grid point numbers into the base point grid point array to obtain an initial base point grid point array.
8. A aeroelastic coupling analog interface point selection device, comprising:
the dynamic array creation module is used for acquiring the grid model file information and creating a dynamic array based on the grid model file information; the dynamic array comprises all grid point arrays and basic point grid point arrays;
the base point and grid point selection module is used for selecting a preset number of base point and grid points from all grid points in the grid model file information based on a preset base point selection method, and filling the base point and grid points into the base point and grid point array to obtain an initial base point and grid point array;
the target grid point determining module is used for determining one target grid point from all the grid points based on a preset base point grid point determining method, and filling the target grid point into the initial base point grid point array to obtain a target base point grid point array;
And the ending point selecting module is used for judging whether the point selecting ending rule is met based on the target base point grid point array, and ending the point selecting operation if the point selecting ending rule is met.
9. An electronic device, comprising:
a memory for storing a computer program;
a processor for executing the computer program to implement the steps of the aeroelastic coupling simulation interface point selection method as claimed in any one of claims 1 to 7.
10. A computer-readable storage medium storing a computer program; wherein the computer program when executed by a processor implements the aeroelastic coupling simulation interface point selection method according to any of claims 1 to 7.
CN202311525648.6A 2023-11-16 2023-11-16 Method, device, equipment and medium for selecting aeroelastic coupling simulation interface point Active CN117235832B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311525648.6A CN117235832B (en) 2023-11-16 2023-11-16 Method, device, equipment and medium for selecting aeroelastic coupling simulation interface point

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311525648.6A CN117235832B (en) 2023-11-16 2023-11-16 Method, device, equipment and medium for selecting aeroelastic coupling simulation interface point

Publications (2)

Publication Number Publication Date
CN117235832A true CN117235832A (en) 2023-12-15
CN117235832B CN117235832B (en) 2024-03-12

Family

ID=89086623

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311525648.6A Active CN117235832B (en) 2023-11-16 2023-11-16 Method, device, equipment and medium for selecting aeroelastic coupling simulation interface point

Country Status (1)

Country Link
CN (1) CN117235832B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115099095A (en) * 2022-06-27 2022-09-23 中国空气动力研究与发展中心计算空气动力研究所 Structure flexibility matrix generation method based on structure finite element model and related components
CN116227043A (en) * 2023-05-10 2023-06-06 中国空气动力研究与发展中心计算空气动力研究所 Aircraft numerical simulation method, system, equipment and computer storage medium
CN116227209A (en) * 2023-03-13 2023-06-06 厦门大学 Multi-dimensional linear difference method for point cloud data, terminal equipment and storage medium
CN116415470A (en) * 2023-06-06 2023-07-11 中国空气动力研究与发展中心计算空气动力研究所 Boundary point extraction method, device and medium for structure finite element
CN116756851A (en) * 2023-06-27 2023-09-15 四川数字经济产业发展研究院 Parameterized grid deformation method and system based on NFFD background grid

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115099095A (en) * 2022-06-27 2022-09-23 中国空气动力研究与发展中心计算空气动力研究所 Structure flexibility matrix generation method based on structure finite element model and related components
CN116227209A (en) * 2023-03-13 2023-06-06 厦门大学 Multi-dimensional linear difference method for point cloud data, terminal equipment and storage medium
CN116227043A (en) * 2023-05-10 2023-06-06 中国空气动力研究与发展中心计算空气动力研究所 Aircraft numerical simulation method, system, equipment and computer storage medium
CN116415470A (en) * 2023-06-06 2023-07-11 中国空气动力研究与发展中心计算空气动力研究所 Boundary point extraction method, device and medium for structure finite element
CN116756851A (en) * 2023-06-27 2023-09-15 四川数字经济产业发展研究院 Parameterized grid deformation method and system based on NFFD background grid

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHENGWEI ZHOU 等: "Direction-of-Arrival Estimation for Coprime Array via Virtual Array Interpolation", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》, vol. 66, no. 22, pages 5956, XP055950239, DOI: 10.1109/TSP.2018.2872012 *
孙岩 等: "径向基/超限插值结构网格运动方法中的变形一致性问题", 《计算机辅助设计与图形学学报》, vol. 30, no. 2, pages 354 - 359 *

Also Published As

Publication number Publication date
CN117235832B (en) 2024-03-12

Similar Documents

Publication Publication Date Title
CN113516248B (en) Quantum gate testing method and device and electronic equipment
JP2022003501A (en) Quantum circuit simulation method, device, facility, storage medium, and program
US8924331B2 (en) System and method for solving quadratic programming problems with bound constraints utilizing a semi-explicit quadratic programming solver
CN116227043B (en) Aircraft numerical simulation method, system, equipment and computer storage medium
CN111185909A (en) Robot operation condition acquisition method and device, robot and storage medium
CN115293623A (en) Training method and device for production scheduling model, electronic equipment and medium
CN114564787B (en) Bayesian optimization method, device and storage medium for target related airfoil design
CN115099095A (en) Structure flexibility matrix generation method based on structure finite element model and related components
CN117235832B (en) Method, device, equipment and medium for selecting aeroelastic coupling simulation interface point
CN107153755B (en) Solving method for shale gas well numerical simulation
CN112100877B (en) Structural rigidity efficient topology optimization method and system
CN112270083B (en) Multi-resolution modeling and simulation method and system
CN110887490A (en) Key frame selection method, medium, terminal and device for laser positioning navigation
CN114818128B (en) Modeling method and optimizing method for ship body local curved surface optimizing neural network
CN112101538B (en) Graphic neural network hardware computing system and method based on memory computing
JP2002297678A (en) Coordinate determination processing method for voxel model, coordinate determination processing program, and recording medium for coordinate determination processing program
CN117078825B (en) Rendering modification method, system, equipment and medium on point cloud data line
CN115017842B (en) Pneumatic data interpolation method and device, electronic equipment and storage medium
CN117875220B (en) Data monitoring method, equipment and medium based on flow field numerical simulation
CN115995279B (en) Material mechanical property evaluation method, device, equipment and readable storage medium
CN111460732B (en) Construction method of nonlinear model of planar motor
JP7331877B2 (en) System computing device and system computing program
CN115995278B (en) Method, device, equipment and readable storage medium for evaluating thermodynamic characteristics of material
CN118095401B (en) Method for accelerating post-state off-track strategy reinforcement learning training of warehouse storage
CN114970052B (en) Simulation method based on node integration, computer equipment and readable storage medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant