Automatic conversion method and system for urban building group load data stream in large-scale numerical simulation

Abstract

The application discloses a kind of urban building group load data stream automatic conversion method and system in large-scale numerical simulation, solve the jamming problem of load data stream transmission in large-scale numerical simulation flow field calculation link and response calculation link, realize the seamless continuation of data stream, so that building group damage assessment method based on multi-degree-of-freedom layer model can fully inherit and utilize CFD fine flow field calculation result, effectively ensure the accuracy and scientificity of building group structure response calculation and damage assessment, blast-resistant / wind-resistant reinforcement design and other series of work, specifically including building group surface measuring point and accessory information module, building group surface load data module, building contour inside and outside identification module, building group surface load direction adjustment module, load data stream type conversion module, load data stream vector operation module, data automatic output module 7 functional modules.

Description

TECHNICAL FIELD

[0001] The present application relates to a data stream conversion method and system, in particular to a large-scale numerical simulation of urban building group load data stream automatic conversion method and system. BACKGROUND

[0002] Urban building group is the main damage target of mid-low altitude / ground large equivalent explosion air shock wave. In the study of damage assessment of explosion shock wave on urban building group, the problem can usually be divided into two main links: shock wave flow field transmission environment calculation (referred to as "flow field calculation") and building group structure response calculation (referred to as "response calculation"). The former belongs to the field of fluid mechanics, aiming to reproduce the propagation and evolution process of shock wave in complex urban space by using advanced computational fluid dynamics (CFD), and to obtain the explosion load borne by the building group; the latter belongs to the field of solid mechanics, aiming to solve the dynamic response of building group under specific excitation load, and to obtain the physical parameters for damage assessment. Although the scientific problems, calculation principles and calculation methods involved in the above two links are significantly different, from the perspective of explosion shock wave damage assessment technology on large-scale urban building group, the two links are successive and indispensable. However, the existing research work or technology is each for war. Either single focus on flow field calculation, such as developing and using various CFD programs, modeling and calculating three-dimensional urban space containing a large number of building groups, extracting explosion load of building group surface, forming load data stream represented by pressure, and each load data reflecting the force condition of building at different positions in a local range. Or single focus on dynamic response research of building group, develop multiple degrees of freedom (MDOF) elastic-plastic layer model suitable for different structure types. However, the load data stream required for layer model response analysis is in the form of concentrated force, which is obviously different from the type of load data stream output by CFD.

[0003] As can be seen, the load output of flow field calculation link cannot be directly used as the load input of response calculation link, and there is a bottleneck between the two main links. So far, there is no available technical solution or method to break through the above bottleneck in the existing public literature, and even there is no related attempt. This makes the building group damage assessment based on multiple degrees of freedom layer model unable to directly and fully use the CFD calculation results, which seriously affects the accuracy and scientificity of building group structure response calculation and damage assessment. SUMMARY

[0004] The application aims to provide a method and system for automatically converting load data flow of urban building groups in large-scale numerical simulation, so as to solve the problem of the bottleneck between the flow field calculation link and the response calculation link, realize seamless connection of load data flow of two key links, and make the building group damage assessment based on the multi-degree-of-freedom layer model directly and fully use the CFD calculation results.

[0005] In order to achieve the above-mentioned purpose, the application adopts the following technical scheme:

[0006] The application discloses a method for automatically converting load data flow of urban building groups in large-scale numerical simulation, which is characterized by comprising the following steps:

[0007] Step 1: obtaining the total number of surface measuring points of each building in the building group to be processed, the side length of the rectangular region represented by the surface measuring points, the surface measuring point pressure load, the total number of pressure loads, the contour centroid coordinates, the coordinates of the two end points of the contour line segment and the coordinates of the midpoint thereof, and making the total number of surface measuring points equal to the total number of pressure loads;

[0008] Step 2: traversing all the buildings in the building group to be processed one by one; for each building, obtaining the concentrated force load at each surface measuring point of the building based on the side length of the rectangular region represented by the surface measuring points and the surface measuring point pressure load, superimposing all the concentrated force loads of the building on the same contour line at the same height, and then superimposing the concentrated force loads on the same contour line according to the floor height of the building, so as to obtain the load data flow A2 of the building; and further obtaining the load data flow A2 of all the buildings;

[0009] Step 3: traversing all the buildings in the building group to be processed one by one; for each building, calculating the inside-outside surface identification data of all the polygonal contours; calculating the directional load data component A5 along the X axis and the directional load data component A6 along the Y axis of all the contour line segments of the building by using the load data flow A2 of the building and the inside-outside surface identification data of the polygonal contours; and further obtaining the directional load data component A5 along the X axis and the directional load data component A6 along the Y axis of all the contour line segments of all the buildings;

[0010] Step 4: traversing all the buildings in the building group to be processed one by one, and obtaining the equivalent load data flow of all the floors or roofs of each building;

[0011] 4.1: selecting a building from the building group to be processed, superimposing all the directional load data components A5 along the X axis on the same floor or roof of the building to obtain the X-axis directional component, and superimposing all the directional load data components A6 along the Y axis on the same floor or roof to obtain the Y-axis directional component;

[0012] 4.2, Synthesize the X-axis direction component force and Y-axis direction component force of the selected building at the same floor or roof to obtain the equivalent load data stream of all floors or roofs of the building;

[0013] 4.3, Traverse all buildings and execute steps 4.1-4.2 for each building until all equivalent load data streams of the building group to be processed are obtained, and the automatic conversion of the load data stream of the building group to be processed is completed.

[0014] Further, the step 3 specifically comprises:

[0015] 3.1, Select a building from the building group to be processed, and for any polygonal contour of the building, calculate the straight line equation A passing through the two end points of each contour line segment, and calculate the straight line equation B passing through the contour centroid coordinates and the midpoint coordinates of each contour line segment;

[0016] Move the straight line equation A with a slope greater than zero in the negative direction of the X-axis by δ, 0mm<δ<10mm, to obtain the straight line equation C, calculate the intersection point of the straight line equation B and the straight line equation C, and use the ray method to determine whether the intersection point is inside the polygonal contour, if the intersection point is inside the polygonal contour, generate the array [-1, 1]; if the intersection point is outside the polygonal contour, generate the array [1, -1];

[0017] Move the straight line equation A with a slope less than zero in the positive direction of the X-axis by δ, 0mm<δ<10mm, to obtain the straight line equation D, calculate the intersection point of the straight line equation B and the straight line equation D, and use the ray method to determine whether the intersection point is inside the polygonal contour, if the intersection point is inside the polygonal contour, generate the array [1, 1]; if the intersection point is outside the polygonal contour, generate the array [-1, -1];

[0018] 3.2, Traverse all contour line segments of the building and execute step 3.1 for each contour line segment to generate all arrays;

[0019] 3.3, Combine all generated arrays and record as the inside-outside surface identification data of all polygonal contours of the building;

[0020] 3.4, Decompose the load data stream A2 corresponding to the building at the floor or roof contour line segment based on the plane coordinate system corresponding to the surface measurement points of the building to obtain the load data component A3 along the X-axis and the load data component A4 along the Y-axis, respectively;

[0021] 3.5, multiply the load data components A3 along the X axis and A4 along the Y axis of all contour line segments of the building with the corresponding polygon contour inside and outside face identification data in turn, to obtain the directional load data components A5 along the X axis and A6 along the Y axis of all contour line segments of the building;

[0022] 3.6, traverse all buildings in the building group to be processed, and execute steps 3.1-3.5 for each building to obtain the directional load data components A5 along the X axis and A6 along the Y axis of all contour line segments of all buildings.

[0023] Further, step 2 specifically includes:

[0024] 2.1, select a building from the building group to be processed, calculate the rectangular area by the side length of the rectangular area represented by the surface measuring point of the building, multiply the pressure load of each surface measuring point with the corresponding rectangular area to obtain the concentrated force load at each surface measuring point;

[0025] 2.2, superimpose all concentrated force loads of the building at the same height and on the same contour line to obtain the load data stream A1;

[0026] 2.3, superimpose the load data stream A1 of the same contour line of the building according to the floor height of the building to obtain the load data stream A2 of the building;

[0027] 2.4, traverse all buildings in the building group to be processed one by one, and execute steps 2.1-2.3 for each building until the load data stream A2 of all buildings is obtained.

[0028] Further, step 1 specifically includes:

[0029] 1.1, prepare the surface measuring point coordinate and auxiliary information file of the building group to be processed, the surface load data file of the building group to be processed; the surface measuring point coordinate and auxiliary information file of the building group to be processed includes the rectangular area side length represented by the surface measuring point of each building in the building group to be processed, the total number of surface measuring points, the contour centroid coordinates, the contour line segment two end point coordinates and the midpoint coordinates; the surface load data file of the building group to be processed includes the pressure load of the surface measuring point of each building in the building group to be processed, the total number of pressure loads;

[0030] 1.2, extract the total number of surface measuring points of each building from the surface measuring point coordinate and auxiliary information file, and extract the total number of pressure loads of each building from the surface load data file;

[0031] 1.3, compare the total number of surface points with the total number of pressure loads; if they are equal, execute step 2; if they are not equal, return to step 1.1 and update the surface point coordinate and accessory information file of the building group to be processed, the surface load data file of the building group to be processed until they are equal.

[0032] The automatic conversion system of urban building group load data stream in large-scale numerical simulation has the characteristics that it comprises an input layer, a conversion layer and an output layer connected in sequence.

[0033] The input layer comprises a building group surface point and accessory information module and a building group surface load data module.

[0034] The conversion layer comprises a building contour inside and outside surface identification module, a building group surface load direction adjustment module, a load data stream type conversion module and a load data stream vector operation module.

[0035] The output layer comprises a data automatic output module.

[0036] The first output end of the building group surface point and accessory information module is connected with the input end of the building contour inside and outside surface identification module, and the second output end is connected with the first input end of the load data stream type conversion module; the building group surface point and accessory information module is used to obtain the surface point coordinate and its accessory information of the building group to be processed.

[0037] The output end of the building group surface load data module is connected with the second input end of the load data stream type conversion module; the building group surface load data module is used to obtain the surface load data stream of the building group to be processed.

[0038] The output end of the building contour inside and outside surface identification module is connected with the first input end of the building group surface load direction adjustment module; the building contour inside and outside surface identification module is used to generate the inside and outside surface identification data of each building polygon contour.

[0039] The output end of the load data stream type conversion module is connected with the second input end of the building group surface load direction adjustment module; the load data stream type conversion module is used to convert the type of load data stream.

[0040] The output end of the building group surface load direction adjustment module is connected with the input end of the load data stream vector operation module; the building group surface load direction adjustment module is used to generate the directional load data component of all contour lines of all buildings.

[0041] The output end of the load data stream vector operation module is connected with the data automatic output module; the load data stream vector operation module is used to generate all equivalent load data streams of the building group to be processed.

[0042] The data automatic output module is used for automatically outputting and saving the obtained equivalent load data stream result.

[0043] The present application has the following beneficial effects:

[0044] 1. The large-scale numerical simulation urban building group load data stream automatic conversion method and system provided by the present application realizes seamless connection from the flow field calculation link to the response calculation link in large-scale numerical simulation, so that the building group damage assessment method based on the multi-degree-of-freedom layer model can fully inherit and utilize the fine flow field calculation results of CFD, and effectively ensure the accuracy and scientificity of building group structure response calculation and damage assessment, blast resistance / wind resistance reinforcement design and other series of work.

[0045] 2. The large-scale numerical simulation urban building group load data stream automatic conversion method and system provided by the present application strictly follows the basic mechanics principle in the conversion process, and the load data stream obtained after conversion fully reflects the load environment characteristics (such as load intensity, load time change history and load point density, etc.) and building structure characteristics (such as shape contour, geometric size, number of floors and floor height, etc.) faced by the building group, and can be used as scientific and reliable equivalent load input required for building structure layer model response calculation.

[0046] 3. The large-scale numerical simulation urban building group load data stream automatic conversion method and system provided by the present application realizes full-process automation and is simple to operate; only the building group surface measuring points and auxiliary information and the CFD building group surface load data need to be provided by the user, and then one-key batch automatic conversion of all building load data can be completed, without the need for the user to specify other parameters or make any other intervention.

[0047] 4. The large-scale numerical simulation urban building group load data stream automatic conversion method and system provided by the present application is highly efficient in data conversion and has no special requirements for computer hardware. According to test, about 200,000 pieces of building group load data can be automatically converted within 15 minutes by using an ordinary office computer, and the data conversion success rate is as high as 100%.

[0048] 5. The large-scale numerical simulation urban building group load data stream automatic conversion method and system provided by the present application is robust and can automatically convert building group load data streams of different scales, geometric shapes and arrangement modes in large-scale numerical simulation; in addition, the present application is suitable for conversion of building group explosion load data and conversion of building group wind load data. BRIEF DESCRIPTION OF DRAWINGS

[0049] Figure 1 is the large-scale numerical simulation urban building group load data stream automatic conversion system architecture diagram of the present application;

[0050] Figure 2 is a flow chart of an embodiment of the automatic conversion method of urban building group load data stream in large-scale numerical simulation of the present application;

[0051] Figure 3 is a building group layout map in a 2.5km 2 core area of a certain city in an embodiment of the present application;

[0052] Figure 4 is a surface pressure measurement point distribution map of a building group in a 2.5km 2 core area of a certain city in an embodiment of the present application;

[0053] Figure 5 is a typical explosion load data of a building group surface in an embodiment of the present application; wherein (a) is a typical measurement point position diagram of four buildings in the building group; (b) is an explosion pressure load curve diagram at the position of measurement point A; (c) is an explosion pressure load curve diagram at the positions of measurement points B-D;

[0054] Figure 6 is an equivalent concentrated force load curve diagram obtained by conversion at different floors of a certain 48-story high-rise building in an embodiment of the present application. DETAILED DESCRIPTION

[0055] The technical solutions of the present application will be described clearly and completely in combination with the drawings and embodiments. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments of the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

[0056] As shown in Figure 1 , the system architecture of the automatic conversion system of urban building group load data stream in large-scale numerical simulation provided by the present embodiment is composed of three parts of input layer, conversion layer and output layer, and contains seven functional modules. The functional modules are connected with each other and connected with each other, and finally achieve the purpose of automatic conversion of large-scale urban building group load data stream. The functions of the modules and the association relationship between the modules are as follows:

[0057] (1) The input layer contains all the input data required by the present application, which is composed of two functional modules of building group surface measurement point and accessory information module ① and building group surface load data module ②.

[0058] The building group surface measuring point and accessory information module ① reads a data file containing the coordinates of the surface measuring points around the building and related accessory information. The data file uses a commonly used data storage format, such as txt / csv / dat / mat, etc. Each piece of data in the data file represents all the information of a measuring point arranged on the surface of the building, including the number of the measuring point, the number of the building where the measuring point is located, the number of the contour line segment of the building where the measuring point is located, the coordinates of the characteristic points of the contour line segment (the coordinates of the two end points and the coordinates of the midpoint), the coordinates of the contour centroid, the three-dimensional coordinates of the measuring point, the side length of the rectangular region represented by the surface measuring point, the height of the building, and the data of the layers.

[0059] The building group surface load data module ② reads a file containing the surface load data stream of the building group (including the surface measuring point pressure load of each building and the total number of pressure loads). The data file uses a commonly used data storage format, such as txt / csv / dat / mat, etc. The file is a result file automatically output by various CFD programs after calculating the flow field transmission environment of a complex urban space building group. In the result file, the load time history data at each measuring point location is stored in order from small to large according to the measuring point number.

[0060] (2) The conversion layer is the core calculation link for automatically converting the load data stream of a large-scale urban building group, which is composed of four functional modules: the building contour inside and outside identification module ③, the building group surface load direction adjustment module ④, the load data stream type conversion module ⑤, and the load data stream vector operation module ⑥.

[0061] The building contour inside and outside identification module ③ is used to identify the inside and outside of all building contours and generate the inside and outside identification data of all building polygonal contours. Based on the building group surface measuring point and accessory information, this module uses the ray method to distinguish the inside and outside of each polygonal contour of each building one by one, which is used for subsequent adjustment of the direction of the building group surface load.

[0062] The building group surface load direction adjustment module ④ is used to adjust the direction of the building group surface load. Based on the building contour inside and outside identification data obtained by the building contour inside and outside identification module ③, this module adjusts the direction of the concentrated force along the coordinate axis on all polygonal contours of the building group, and generates the directional load data component of the building group polygonal contour on the coordinate axis.

[0063] The load data stream type conversion module 5 is mainly used for converting the type of load data stream. The purpose is to convert the uniform pressure load represented by the mass of measuring points around the building into a series of concentrated force loads located at the contour line segments of the building floor or roof. The module includes the steps of converting the uniform pressure load of the measuring points into the concentrated force load at the location, superimposing the concentrated force loads on the same contour line segment at the same height, superimposing the concentrated force along the floor height on the same contour line segment and concentrating at the location of each floor or roof, etc. The conversion object of the load data stream type conversion module 5 is the building group surface load data in the building group surface load data module 2, i.e. the original load data without load direction adjustment.

[0064] The load data stream vector operation module 6 is a load data stream vector operation module. The module takes the data stream after load data stream type conversion and direction adjustment as the input object, and superimposes and synthesizes the concentrated force loads on different contour line segments at the same floor or roof position of the building along the plane coordinate system according to the operation rule of vector force, and finally generates only one equivalent load data in the form of concentrated force at each floor or roof position.

[0065] (3) The output layer only contains the data automatic output module 7, which is mainly used for automatically outputting and saving the results of the load data stream vector operation module 6 to the local hard disk. The results will be stored in the commonly used data format file, which is convenient for subsequent building structure response calculation.

[0066] The connection mode of all the above functional modules is as follows:

[0067] The first output end of the building group surface measuring point and the attached information module 1 is connected with the input end of the building contour inside and outside surface identification module 3, and the second output end is connected with the first input end of the load data stream type conversion module 5;

[0068] The output end of the building group surface load data module 2 is connected with the second input end of the load data stream type conversion module 5; the output end of the building surface inside and outside surface identification module 3 is connected with the first input end of the building group surface load direction adjustment module 4; the output end of the load data stream type conversion module 5 is connected with the second input end of the building group surface load direction adjustment module 4; the output end of the building group surface load direction adjustment module 4 is connected with the input end of the load data stream vector operation module 6; the output end of the load data stream vector operation module 6 is connected with the data automatic output module 7.

[0069] As shown in Figure 2 The working process of the above-mentioned automatic conversion of the load data stream of the urban building group in large-scale numerical simulation is as follows:

[0070] Step 1, prepare the surface point coordinate and the attached information file of the building group to be processed, the surface load data file of the building group to be processed, as the input layer building group surface point and the attached information module ① and the building group surface load data module ② of the read object file; the surface point coordinate and the attached information file of the building group to be processed include the edge length of the rectangular area represented by the surface point of each building in the building group to be processed, the total number of surface points, the contour centroid coordinates, the coordinates of the two end points and the midpoint coordinates of the contour line segment; the surface load data file of the building group to be processed includes the surface point pressure load of each building in the building group to be processed, the total number of pressure loads;

[0071] Step 2, extract the total number of surface points from the surface point coordinate and the attached information file, and extract the total number of pressure loads from the surface load data file;

[0072] Step 3, compare the total number of surface points with the total number of pressure loads;

[0073] If they are not equal, update the surface point coordinate and the attached information file of the building group to be processed, the surface load data file of the building group to be processed, and execute step 2 until they are equal, and execute step 4;

[0074] If they are equal, execute step 4;

[0075] Step 4, select a building from the building group to be processed through the load data stream type conversion module ⑤, and obtain the load data stream A2 of the building;

[0076] 4.1, and obtain the edge length of the rectangular area represented by each surface point of the building from the surface point coordinate and the attached information file, calculate its area, and obtain the pressure load of each surface point of the building from the surface load data file, multiply each surface point pressure load with the corresponding area to obtain the concentrated force load at each surface point;

[0077] 4.2, superimpose all the concentrated force loads of the building at the same height and on the same contour line to obtain the load data stream A1;

[0078] 4.3, superimpose the load data stream A1 of the same contour line of the building according to the floor height of the building to obtain the load data stream A2 of the building;

[0079] Step 5, traverse all the buildings in the building group to be processed, and execute step 4 for each building until the load data stream A2 of all the buildings is obtained, and execute step 6;

[0080] Step 6, select a building from the building group to be processed by the building surface inside and outside identification module ③ and the building group surface load direction adjustment module ④, and obtain the directional load data component A5 along the X axis and the directional load data component A6 along the Y axis of all contour line segments of the building;

[0081] 6.1, extract the contour centroid coordinates, the coordinates of the two end points of any contour line segment, and the coordinates of the midpoint of the building from the surface measurement point coordinates and the auxiliary information file, calculate the straight line equation A passing through the two end point coordinates, and calculate the straight line equation B passing through the centroid coordinates and the midpoint coordinates;

[0082] Move the straight line equation A with a slope greater than zero in the negative direction of the X axis by δ, 0mm<δ<10mm, to obtain the straight line equation C, calculate the intersection point of the straight line equation B and the straight line equation C, and use the ray method to determine whether the intersection point is inside the polygon contour. If the intersection point is inside the polygon contour, generate the array [-1, 1]; if the intersection point is outside the polygon contour, generate the array [1, -1];

[0083] Move the straight line equation A with a slope less than zero in the positive direction of the X axis by δ, 0mm<δ<10mm, to obtain the straight line equation D, calculate the intersection point of the straight line equation B and the straight line equation D, and use the ray method to determine whether the intersection point is inside the polygon contour. If the intersection point is inside the polygon contour, generate the array [1, 1]; if the intersection point is outside the polygon contour, generate the array [-1, -1];

[0084] 6.2, traverse all contour line segments of the building, and execute step 6.1 for each contour line segment to generate all arrays;

[0085] 6.3, combine all generated arrays, and mark as the inside and outside surface identification data of all polygon contours of the building;

[0086] 6.4, decompose the load data stream A2 corresponding to the building at the floor or roof contour line segment based on the plane coordinate system corresponding to the surface measurement point, to obtain the load data component A3 along the X axis and the load data component A4 along the Y axis, respectively;

[0087] 6.5, multiply the load data components along the X axis and the Y axis of all contour line segments of the building with the polygon contour inside and outside surface identification data in turn, to obtain the directional load data component A5 along the X axis and the directional load data component A6 along the Y axis of all contour line segments of the building;

[0088] Step 7, traverse all buildings in the building group to be processed, execute step 6 for each building to obtain the directional load data component A5 along the X axis and the directional load data component A6 along the Y axis of all contour line segments of all buildings;

[0089] Step 8, select a building through the load data flow vector operation module ⑥, and obtain the equivalent load data flow at all floors or roofs of the building;

[0090] 8.1, vector superimpose all X-axis directional load data components A5 at the same floor or roof to obtain the X-axis directional component force; and vector superimpose all Y-axis directional load data components A6 at the same floor or roof to obtain the Y-axis directional component force;

[0091] 8.2, synthesize the X-axis directional component force and the Y-axis directional component force at the same floor or roof of the building to obtain the equivalent load data flow at all floors or roofs of the building;

[0092] Step 9, traverse all buildings, and perform step 8 on each building until all equivalent load data flows of the building group to be processed are obtained, and the automatic conversion of the load data flow of the building group to be processed is completed.

[0093] It should be particularly emphasized that the present embodiment mainly takes the load data flow in the explosion scene as an example, and the method proposed is not limited to the explosion scene, and is also applicable to the conversion of complex urban space wind field CFD calculation results.

[0094] To further verify the effectiveness of the method provided by the present application, a building group in a 2.5kmx2.5km core area of a city is selected as a typical case, and the conversion process and technical effect of the explosion load data flow of the building group in large-scale numerical simulation are detailed.

[0095] (1) Basic situation

[0096] The region contains 1759 buildings in total, with building floor numbers ranging from 1 to 70, a total of 11142 floors, and single building area ranging from 10.4 to 17904.1m 2 , and the overall layout is shown in Figure 3 . When using the CFD program to calculate the explosion flow field transmission environment of the urban space in the region, in order to extract the reflected shock wave overpressure time history curve of the building group surface, a large number of regular and uniform pressure measuring points are arranged around the building surface. In this case, a total of 193992 measuring points are arranged, as shown in Figure 4 .

[0097] After using the CFD program to complete the calculation of the explosion flow field transmission environment, the explosion load borne by the building surface at each measuring point position can be obtained. The explosion loads of 4 randomly selected measuring points from the 193992 measuring points are shown in Figure 5 .

[0098] As shown in Figure 5It can be seen that the explosion load data at the measuring point is the time history data of the reflected shock wave overpressure, and each curve only represents the uniform pressure load borne by the building surface in the specific rectangular area around the measuring point. Obviously, this type of load data stream cannot be directly used as the load input for the building structure layer model response calculation. By using the load data stream conversion method and system provided in this embodiment, according to the building geometric profile, layer height, number of floors, and layout information of the pressure measuring point, the concentrated force load data stream matched with the building multi-degree-of-freedom layer model can be automatically converted and output.

[0099] (2) Input layer data

[0100] The input layer in the system architecture provided in this embodiment includes two functional modules, which correspond to different data files respectively.

[0101] The “building group surface measuring point and accessory information” file: This file contains rich information such as the number of the surface pressure measuring point around the building group, the building number, the building profile line segment number, the building profile line segment feature point coordinates, the building profile centroid coordinates, the measuring point three-dimensional coordinates, the measuring point coverage area size, the building height, and the layer height. Each row of data represents the position coordinates and accessory information of a measuring point, and the units of the measuring point position coordinates and length in the file are “meters”.

[0102] The “building group surface load data” file: This file is a data file output by the open source CFD after calculating the explosion flow field environment of the building group, and contains the explosion load time history curve of all the surface pressure measuring points around the building.

[0103] (3) Output layer data

[0104] The method proposed in this embodiment converts the load data stream of the input layer, and after successful conversion, a new load data stream file “building group equivalent concentrated force load data” is automatically generated, which is saved in the csv format and contains the equivalent concentrated force load data matched with the total number of floors of the urban building group. In this embodiment, it is statistically found that the file contains 11,142 load data, which is consistent with the total number of floors of the building group in the region.

[0105] To further show the conversion effect of the load data stream, a 48-story high-rise building is randomly selected from the 1,759 buildings in the region, and the equivalent concentrated force load curves of the first floor (bottom floor), the 20th floor, the 40th floor, and the 48th floor (top floor) of the building are extracted from the “building group equivalent concentrated force load data” file, as shown in Figure 6 Due to the shielding of other buildings around, the concentrated force explosion load borne by different floors of the same building is quite different, showing obvious shock characteristics, which is consistent with the actual explosion load variation law of the complex urban space building group.

[0106] Through the test, the corresponding equivalent concentrated force load data of all other buildings are obtained, and the success rate of load data conversion is as high as 100%. It is particularly emphasized that: this case is completed on a Lenovo computer with an 8-core Intel Core i7-9700 processor and 16 GB of RAM. In less than 15 minutes, the surface of the building group with different geometric sizes and contours is automatically converted from 193992 pressure data to 11142 equivalent concentrated force load data, which shows high data conversion efficiency and strong adaptability. The converted load data stream can be directly used as the load input for the multi-degree-of-freedom layer model response calculation of the building structure, and can provide reliable data support for subsequent urban building group damage assessment and blast-resistant reinforcement design series research work.

[0107] The above merely illustrates the specific embodiments of the present application, but the protection scope of the present application is not limited thereto, and any change or replacement within the technical scope disclosed by the present application should be covered within the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A method for automatic conversion of urban building cluster load data stream in large-scale numerical simulation, characterized in that, The method comprises the following steps: Step 1, obtaining the total number of surface measuring points of each building in the building group to be processed, the side length of the rectangular area represented by the surface measuring points, the surface measuring point pressure load, the total number of pressure loads, the profile centroid coordinates, the coordinates of the two end points of the profile segment and the coordinates of the midpoint thereof; making the total number of surface measuring points equal to the total number of pressure loads; Step 2, traversing all the buildings in the building group to be processed one by one; for each building, based on the side length of the rectangular area represented by the surface measuring points and the surface measuring point pressure load, obtaining the concentrated force load at each surface measuring point of the building, superimposing all the concentrated force loads of the building on the same contour line at the same height, and then superimposing all the concentrated force loads on the same contour line according to the floor height of the building, to obtain the load data stream A2 of the building; and further obtaining the load data stream A2 of all the buildings; Step 3, traversing all the buildings in the building group to be processed one by one; for each building, calculating the inside-outside surface identification data of all the polygonal contours thereof; through the load data stream A2 of the building and the inside-outside surface identification data of the polygonal contour, calculating the directional load data component A5 along the X axis and the directional load data component A6 along the Y axis of all the profile segments of the building; and further obtaining the directional load data component A5 along the X axis and the directional load data component A6 along the Y axis of all the profile segments of all the buildings; Step 4, traversing all the buildings in the building group to be processed one by one, and obtaining the equivalent load data stream at all the floors or roofs of each building; 4.1, selecting a building from the building group to be processed, vector superimposing all the directional load data components A5 along the X axis on the same floor or roof of the building to obtain the X axis directional component force, and vector superimposing all the directional load data components A6 along the Y axis on the same floor or roof to obtain the Y axis directional component force; 4.2, synthesizing the X axis directional component force and the Y axis directional component force on the same floor or roof of the selected building to obtain the equivalent load data stream at all the floors or roofs of the building; 4.3, traversing all the buildings, and performing steps 4.1-4.2 on each building, until all the equivalent load data streams of the building group to be processed are obtained, and the automatic conversion of the load data stream of the building group to be processed is completed. Step 3 specifically comprises:

2. The method for automatic conversion of urban building cluster load data stream in large-scale numerical simulation according to claim 1, characterized in that, 3.1, selecting a building from the building group to be processed, for any polygonal contour of the building, calculating the straight line equation A passing through the two end points of each profile segment, and calculating the straight line equation B passing through the profile centroid coordinates and the midpoint coordinates of each profile segment; moving the straight line equation A with a slope greater than zero along the negative direction of the X axis by δ, 0mm<δ<10mm, to obtain the straight line equation C, calculating the intersection point of the straight line equation B and the straight line equation C, and judging whether the intersection point is inside the polygonal contour by using the ray method, if the intersection point is inside the polygonal contour, generating the array [-1, 1]; if the intersection point is outside the polygonal contour, generating the array [1, -1]; ​ Move the straight line equation A with a slope less than zero along the positive direction of the X axis by δ, 0mm<δ<10mm, to obtain a straight line equation D, calculate the intersection point of the straight line equation B and the straight line equation D, and determine whether the intersection point is inside the polygon contour by using the ray method, if the intersection point is inside the polygon contour, generate an array [1, 1]; if the intersection point is outside the polygon contour, generate an array [-1, -1]; 3.2, traverse all contour line segments of the building, and perform step 3.1 for each contour line segment to generate all arrays; 3.3, combine all generated arrays, and record as the inside and outside surface identification data of all polygon contours of the building; 3.4, decompose the load data stream A2 of the building at the floor or roof contour line segment based on the plane coordinate system corresponding to the surface measurement point of the building to obtain the load data component A3 along the X axis and the load data component A4 along the Y axis, respectively; 3.5, multiply the load data component A3 along the X axis and the load data component A4 along the Y axis of all contour line segments of the building with the corresponding polygon contour inside and outside surface identification data in turn to obtain the directional load data component A5 along the X axis and the directional load data component A6 along the Y axis of all contour line segments of the building; 3.6, traverse all buildings in the to-be-processed building group, and perform steps 3.1-3.5 for each building to obtain the directional load data component A5 along the X axis and the directional load data component A6 along the Y axis of all contour line segments of all buildings.

3. The method for automatic conversion of urban building cluster load data stream in large-scale numerical simulation according to claim 1, characterized in that, Step 2 specifically includes: 2.1, select a building from the to-be-processed building group, calculate the area of the rectangular region represented by the surface measurement points of the building through the side length of the rectangular region, multiply the pressure load of each surface measurement point with the area of the rectangular region to obtain the concentrated force load at each surface measurement point; 2.2, superimpose all concentrated force loads of the building at the same height and on the same contour line to obtain the load data stream A1; 2.3, superimpose the load data stream A1 of the building at the same contour line according to the floor height of the building to obtain the load data stream A2 of the building; 2.4, traverse all buildings in the to-be-processed building group one by one, and perform steps 2.1-2.3 for each building until the load data stream A2 of all buildings is obtained.

4. The method for automatic conversion of urban building cluster load data stream in large-scale numerical simulation according to claim 1 or 2 or 3, characterized in that, Step 1 specifically includes: 1.1, prepare the surface measurement point coordinates and auxiliary information file of the to-be-processed building group, and the surface load data file of the to-be-processed building group; the surface measurement point coordinates and auxiliary information file of the to-be-processed building group includes the side length of the rectangular region represented by the surface measurement points of each building in the to-be-processed building group, the total number of surface measurement points, the contour centroid coordinates, the contour line segment two end point coordinates and the midpoint coordinates; the surface load data file of the to-be-processed building group includes the pressure load of the surface measurement points of each building in the to-be-processed building group, the total number of pressure loads; 1.2, extract the total number of surface measurement points of each building from the surface measurement point coordinates and auxiliary information file, and extract the total number of pressure loads of each building from the surface load data file; 1.3, compare the total number of surface points with the total number of pressure loads; if they are equal, execute step 2; if they are not equal, return to step 1.1 and update the surface point coordinate and accessory information file of the building group to be processed, the surface load data file of the building group to be processed until they are equal.

5. A system for automatic conversion of urban building cluster load data stream in large scale numerical simulation, characterized in that: The input layer, the conversion layer and the output layer are connected in sequence; The input layer includes a building group surface point and accessory information module and a building group surface load data module; The conversion layer includes a building contour inside and outside surface identification module, a building group surface load direction adjustment module, a load data stream type conversion module and a load data stream vector operation module; The output layer includes a data automatic output module; The first output end of the building group surface point and accessory information module is connected with the input end of the building contour inside and outside surface identification module, and the second output end is connected with the first input end of the load data stream type conversion module; the building group surface point and accessory information module is used to obtain the surface point coordinate and its accessory information of the building group to be processed; The output end of the building group surface load data module is connected with the second input end of the load data stream type conversion module; the building group surface load data module is used to obtain the surface load data stream of the building group to be processed; The output end of the building contour inside and outside surface identification module is connected with the first input end of the building group surface load direction adjustment module; the building contour inside and outside surface identification module is used to generate the inside and outside surface identification data of each building polygon contour; The output end of the load data stream type conversion module is connected with the second input end of the building group surface load direction adjustment module; the load data stream type conversion module is used for the conversion of the load data stream type; The output end of the building group surface load direction adjustment module is connected with the input end of the load data stream vector operation module; the building group surface load direction adjustment module is used to generate the directional load data component of all contour lines of all buildings; The output end of the load data stream vector operation module is connected with the data automatic output module; the load data stream vector operation module is used to generate all equivalent load data streams of the building group to be processed; The data automatic output module is used to automatically output and save the obtained equivalent load data stream results.

Images