Curtain wall deformation condition safety detection method and device based on laser point cloud

Abstract

The invention provides a curtain wall deformation condition safety detection method and device based on laser point clouds, and relates to the technical field of constructional engineering, and the method comprises the steps: 1, carrying out the coordinate system unification processing of a three-dimensional point cloud data set covering the whole surface of a curtain wall, generating a curtain wall surface digital twin model, and carrying out the coordinate system unification processing of a three-dimensional point cloud data set covering the whole surface of the curtain wall in the digital twin model; defining a first detection target point and a second detection target point according to geometrical characteristics and material joint characteristics of the curtain wall sheet; 2, according to the digital twinborn model, surface normal vector distribution characteristics are extracted, actually-measured wind pressure load distribution data and temperature field parameters are fused, a dynamic flexural deformation coupling calculation model is constructed, and global deformation parameters and local deformation parameters are output. According to the invention, the digital twinborn model is constructed through point cloud data processing, deformation and displacement parameters are extracted, the safety coefficient is evaluated, and high-precision, omnibearing and dynamic detection and safety risk early warning of the flatness of the building outer wall are realized.

Description

Technical Field

[0001] The present invention relates to the technical field of construction engineering, and particularly to a safety detection method and device for curtain wall deformation based on laser point cloud. Background Art

[0002] Traditional methods for detecting the flatness of building exterior walls, such as using a straightedge for detection and total station measurement, have problems such as low detection efficiency, strong dependence on manual labor, and insufficient data accuracy, making it difficult to meet the detection requirements of modern large and complex curtain wall structures. Existing detection methods often only focus on the geometric shape changes on the surface of the curtain wall, simply obtaining the flatness value through point cloud data, but ignoring the dynamic effects of environmental factors such as wind pressure and temperature on the curtain wall structure under long-term action.

[0003] For example, during strong winds or drastic temperature changes, the actual stress state and deformation trend of the curtain wall cannot be truly reflected by a single geometric parameter, resulting in deviations in the detection results and being unable to accurately evaluate the safety status of the curtain wall.

[0004] In addition, during the process of detecting data processing and analysis, it is impossible to effectively integrate multi-source data information, making it difficult to achieve a comprehensive and dynamic assessment of the safety of the curtain wall structure. Summary of the Invention

[0005] The technical problem to be solved by the present invention is to provide a safety detection method and device for curtain wall deformation based on laser point cloud, to achieve the problem of comprehensive dynamic assessment, and to improve the accuracy and reliability of detecting the flatness of building exterior walls.

[0006] To solve the above technical problems, the technical solution of the present invention is as follows: In the first aspect, a safety detection method for curtain wall deformation based on laser point cloud, the method includes: Step 1, perform coordinate system unification processing on the three-dimensional point cloud data set covering the entire surface of the curtain wall to generate a digital twin model of the curtain wall surface, and in the digital twin model, define the first detection target point and the second detection target point according to the geometric characteristics and material joint characteristics of the curtain wall panels; Step 2, according to the digital twin model, extract the surface normal vector distribution characteristics, and fuse the measured wind pressure load distribution data and temperature field parameters to construct a dynamic flexural deformation coupling calculation model, and output the global deformation parameters and local deformation parameters; Step 3, according to the local deformation parameters, through phase correlation analysis and wavelet transform, extract the vibration displacement components of the first detection target point and the second detection target point to construct a curtain wall vibration modal feature matrix, and at the same time, analyze the displacement-time series data to identify the over-limit displacement events where the displacement amplitude is greater than or equal to the preset threshold; Step 4: Determine the initial coefficient of the curtain wall structure safety degree according to the global deformation parameters and the occurrence frequency of the over-limit displacement events, and calculate the target correction value based on the local deformation parameters and the dynamic displacement components; Step 5: Superimpose the target correction value on the initial coefficient of the safety degree to generate the corrected coefficient of the curtain wall structure safety degree, realize the dynamic adjustment of the initial coefficient of the safety degree, and output a three-dimensional detection report, including the risk level, the distribution of the corrected safety degree coefficient, and the warning marks associated with the target correction value.

[0007] In a second aspect, a safety detection device for the deformation condition of a curtain wall based on laser point cloud includes: A data processing module, configured to perform coordinate system normalization processing on the three-dimensional point cloud data set covering the entire surface of the curtain wall, generate a digital twin model of the curtain wall surface, and define the first detection target point and the second detection target point in the digital twin model according to the geometric characteristics and material joint characteristics of the curtain wall panels; A model construction module, configured to extract the surface normal vector distribution characteristics according to the digital twin model, and fuse the measured wind pressure load distribution data and the temperature field parameters to construct a dynamic flexural deformation coupling calculation model, and output the global deformation parameters and the local deformation parameters; A feature extraction module, configured to extract the vibration displacement components of the first detection target point and the second detection target point through phase correlation analysis and wavelet transform according to the local deformation parameters to construct a vibration mode feature matrix of the curtain wall. At the same time, analyze the displacement-time series data to identify the over-limit displacement events with the displacement amplitude greater than or equal to the preset threshold; A safety degree evaluation module, configured to determine the initial coefficient of the curtain wall structure safety degree according to the global deformation parameters and the occurrence frequency of the over-limit displacement events, and calculate the target correction value based on the local deformation parameters and the dynamic displacement components; A report generation module, configured to superimpose the target correction value on the initial coefficient of the safety degree to generate the corrected coefficient of the curtain wall structure safety degree, realize the dynamic adjustment of the initial coefficient of the safety degree, and output a three-dimensional detection report, including the risk level, the distribution of the corrected safety degree coefficient, and the warning marks associated with the target correction value.

[0008] In a third aspect, a computing device includes: One or more processors; A storage device, configured to store one or more programs, which when executed by the one or more processors, cause the one or more processors to implement the method described above.

[0009] In a fourth aspect, a computer-readable storage medium stores a program, which when executed by a processor, implements the method described above.

[0010] The above solutions of the present invention at least include the following beneficial effects: By performing coordinate system normalization on the three-dimensional point cloud dataset to generate a digital twin model, the consistency and accuracy of the data are ensured, and the surface morphology of the curtain wall can be truly restored. Defining the detection target points based on the geometric characteristics of the curtain wall panels and the material joint characteristics avoids misjudgment caused by data deviation or improper selection of detection points. When constructing the dynamic flexural deformation coupling calculation model, the distribution characteristics of surface normal vectors, the measured wind pressure load distribution data, and the temperature field parameters are integrated, changing the traditional mode of analyzing single or a small number of factors. This comprehensive multi-factor analysis can comprehensively consider various influences on the curtain wall during actual use, and the output global and local deformation parameters are more in line with the actual situation, thus more accurately reflecting the deformation status of the curtain wall.

[0011] Using phase correlation analysis and wavelet transform to extract the vibration displacement components of the detection target points can sensitively capture the fine vibration characteristics of the curtain wall. The constructed vibration modal characteristic matrix provides detailed data for analyzing the dynamic characteristics of the curtain wall. The analysis of displacement-time series data can timely identify over-limit displacement events, providing a dynamic basis for evaluating the safety of the curtain wall and helping to detect potential safety hazards in advance. Determining the initial safety coefficient based on the global deformation parameters and the occurrence frequency of over-limit displacement events, and calculating the target correction value by combining the local deformation parameters and the dynamic displacement components, realizes the scientific evaluation of the structural safety of the curtain wall. Superimposing the target correction value and the initial coefficient to generate the corrected safety coefficient can be dynamically adjusted according to the actual state of the curtain wall, making the evaluation result more in line with the safety performance of the curtain wall at different stages.

[0012] The output three-dimensional detection report contains rich information such as risk levels, the distribution of corrected safety coefficients, and warning marks, presented in an intuitive three-dimensional visualization form, facilitating managers and technicians to quickly understand the safety status of each part of the curtain wall. The warning marks are associated with the target correction value, which can timely prompt potential risk areas, enabling relevant personnel to quickly take measures for handling, reducing the probability of safety accidents, and ensuring the safety of personnel's lives and property and the normal use of the building. Description of the Drawings

[0013] Figure 1 is a schematic flowchart of a method for safely detecting the deformation status of a curtain wall based on laser point cloud provided by an embodiment of the present invention.

[0014] Figure 2 is a schematic diagram of a device for safely detecting the deformation status of a curtain wall based on laser point cloud provided by an embodiment of the present invention. Detailed Embodiments

[0015] Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided so that the present disclosure can be more thoroughly understood and the scope of the present disclosure can be fully conveyed to those skilled in the art.

[0016] As Figure 1 shown, an embodiment of the present invention provides a safety detection method for the deformation condition of a curtain wall based on laser point cloud. The method includes the following steps: Step 1: Perform coordinate system normalization processing on the three-dimensional point cloud data set covering the entire surface of the curtain wall to generate a digital twin model of the curtain wall surface. In the digital twin model, define the first detection target point and the second detection target point according to the geometric characteristics and material joint characteristics of the curtain wall panels. Step 2: According to the digital twin model, extract the surface normal vector distribution characteristics, and fuse the measured wind pressure load distribution data and temperature field parameters to construct a dynamic flexural deformation coupling calculation model, and output the global deformation parameters and local deformation parameters. Step 3: According to the local deformation parameters, extract the vibration displacement components of the first detection target point and the second detection target point through phase correlation analysis and wavelet transform to construct a vibration mode feature matrix of the curtain wall. At the same time, analyze the displacement-time series data to identify the over-limit displacement events where the displacement amplitude is greater than or equal to the preset threshold. Step 4: Determine the initial coefficient of the curtain wall structure safety degree according to the global deformation parameters and the occurrence frequency of the over-limit displacement events, and calculate the target correction value based on the local deformation parameters and the dynamic displacement components. Step 5: Superimpose the target correction value on the initial coefficient of the safety degree to generate a corrected curtain wall structure safety degree coefficient, realize the dynamic adjustment of the initial coefficient of the safety degree, and output a three-dimensional detection report, including the risk level, the distribution of the corrected safety degree coefficient, and the warning mark associated with the target correction value.

[0017] In the embodiment of the present invention, by performing coordinate system normalization processing on the three-dimensional point cloud data set to generate a digital twin model of the curtain wall surface, and defining the first and second detection target points according to the geometric characteristics and material joint characteristics of the curtain wall panels, the key detection areas are accurately anchored, avoiding the blindness of traditional detection. Compared with the traditional method, it can capture the subtle deformation of the curtain wall more carefully, effectively overcome the problem of large manual measurement errors, and improve the accuracy and reliability of the detection data.

[0018] Integrate the measured wind pressure load distribution data and temperature field parameters to construct a dynamic flexural deformation coupling calculation model, breaking through the limitation of existing technologies that only focus on single geometric deformation parameters. This model can simulate the dynamic effects of environmental factors such as wind pressure and temperature on the curtain wall, output global deformation parameters and local deformation parameters, truly restore the stress and deformation states of the curtain wall under complex environments, make the detection results more consistent with the actual working conditions, effectively avoid detection deviations caused by ignoring environmental factors, and comprehensively reflect the true safety state of the curtain wall.

[0019] Apply phase correlation analysis and wavelet transform technology to extract the vibration displacement components of the detection target points to construct a curtain wall vibration modal feature matrix, and at the same time identify out-of-limit displacement events in combination with displacement-time series data. This process not only achieves accurate capture of the dynamic response of the curtain wall but also can timely detect potential safety hazards. Compared with the situation where traditional methods cannot detect dynamic characteristics and existing technical analysis means are insufficient, it enhances the monitoring ability of the curtain wall operation state and provides a strong guarantee for preventing accidents.

[0020] Determine the initial safety coefficient based on the global deformation parameters and the occurrence frequency of out-of-limit displacement events, and calculate the target correction value in combination with the local deformation parameters and dynamic displacement components to dynamically adjust the safety coefficient. This systematic evaluation method effectively integrates multi-source data information. Compared with the defects of traditional detection lacking an evaluation system and existing technologies being difficult to comprehensively and dynamically evaluate, it can quantitatively evaluate the structural safety of the curtain wall from a global perspective, generate a three-dimensional detection report including risk levels, safety coefficient distributions, and warning marks, and provide an intuitive, comprehensive, and scientific decision-making basis for curtain wall maintenance management, improving the efficiency and accuracy of curtain wall maintenance management.

[0021] In a preferred embodiment of the present invention, in step 1 above, perform coordinate system normalization processing on the three-dimensional point cloud data set covering the entire surface of the curtain wall to generate a digital twin model of the curtain wall surface, and in the digital twin model, define the first detection target point and the second detection target point according to the geometric characteristics and material joint characteristics of the curtain wall panels; the first detection target point is located in the four corner regions of the curtain wall panel and is used to monitor the corner stress concentration effect; the second detection target point is located at the junction of the center and the edge of the curtain wall panel and is used to monitor the temperature difference deformation gradient change, which may include: In the embodiments of the present invention, a three-dimensional point cloud dataset covering the entire surface of the curtain wall is obtained, and a laser scanning device (such as a three-dimensional laser scanner) is used to complete this task. The laser scanning device is placed at a suitable position to scan the curtain wall in all directions. During the scanning, the device emits a laser beam and measures the time it takes to reflect back, thereby obtaining the three-dimensional coordinates of numerous points on the surface of the curtain wall, and finally forming a three-dimensional point cloud dataset. This dataset covers rich geometric information on the surface of the curtain wall, but due to factors such as the position and attitude of the scanning device, the point cloud data of each part is in different coordinate systems. Since the obtained three-dimensional point cloud dataset comes from different scanning positions or devices, there are differences in the coordinate systems of the point cloud data of each part. In order to be able to uniformly process and analyze these data, coordinate system unification processing is required. The specific steps are as follows: Extract representative feature points or feature surfaces in the three-dimensional point cloud dataset, such as the corner points and edge lines of the curtain wall. These features will be used as the reference basis for coordinate transformation. Through the iterative closest point algorithm, a suitable rotation and translation transformation matrix is found to transform the point cloud data in different coordinate systems into the same coordinate system. This can ensure that all point cloud data is processed in the same reference system. After completing the coordinate system unification processing, a digital twin model of the curtain wall surface can be generated based on the three-dimensional point cloud dataset in the unified coordinate system. The digital twin model is a digital mapping of the real curtain wall and can accurately reflect the geometric shape and structural information of the curtain wall.

[0022] Specifically, remove the noise points and outliers in the point cloud data to improve the quality and accuracy of the data. Use the point cloud data to reconstruct the surface model of the curtain wall. Surface reconstruction methods include the moving least squares method. Through these methods, the discrete point cloud data can be converted into a continuous surface model to form a digital twin model of the curtain wall surface. The first detection target points are located in the four corner regions of the curtain wall panels and are used to monitor the corner stress concentration effect. The specific process of determining these target points is as follows: Divide the point cloud data into regular three-dimensional voxel grids. By calculating the average value of the points in each voxel, data downsampling is achieved, reducing the data volume while retaining the features. For each point in the point cloud data, construct a covariance matrix in its neighborhood. Perform eigenvalue decomposition on the covariance matrix, and the eigenvector corresponding to the minimum eigenvalue is the normal vector of the marked point. Use the average curvature algorithm to calculate the curvature value of each point in the point cloud , , where and are two eigenvalues of the covariance matrix. Mark the region where the curvature value > the set threshold and the curvature change rate of adjacent points is greater than or equal to as the potential corner region. The curvature change rate of adjacent points The calculation formula is , where is the curvature value of a certain point within the neighborhood of point , is the distance between two points. By applying the edge detection algorithm, the edge intensity and direction of the point cloud are calculated to identify the edge contour of the slab. The edge intensity , where and are the partial derivatives of the curvature value in the and directions respectively. The edge intersection area is also included in the potential corner area.

[0023] Cluster analysis is performed on all potential corner areas. Using clustering algorithms such as DBSCAN, overlapping areas are merged, and finally the four corner areas of the curtain wall slab are determined. Within each corner area, a spherical range with a radius of centered on the area vertex is delimited. Within this range, 3 - 5 points that are closest to the vertex and have the largest curvature change rate are selected as the first detection target points. Let the distance from point to vertex be , , where,( ) and( ) are the coordinates of point and vertex respectively.

[0024] The second detection target point is located at the junction of the center and the edge of the curtain wall slab for monitoring the change in the temperature difference deformation gradient. The specific process of determining these target points is as follows: According to the geometric shape and size of the curtain wall slab, the geometric center of the slab is calculated. If the point cloud data of the slab is , then the coordinate calculation formula of the geometric center is , where is the total number of points in the point cloud data of the curtain wall slab, represents the coordinates of the th point in the point cloud data of the curtain wall slab. By constructing a Voronoi diagram to divide the point cloud data of the slab, combined with the distance from the point cloud to the geometric center and the shortest distance from the point cloud to the edge, points with a distance from the slab center within 0.6 - 0.8 times the slab radius and a distance to the edge less than a certain threshold The points that meet the following conditions are determined as the boundary region between the center and the edge, i.e., satisfying: and . The point cloud in the boundary region is meshed, and the region is divided into small grids. Within each small grid, the temperature sensitivity coefficient of each point is calculated . According to the type of material used for the curtain wall, the thermal expansion coefficient of building curtain wall materials such as aluminum alloy is approximately between (23.6× -24.0× ) / °C. The thermal expansion coefficient of glass varies depending on the type of glass. The thermal expansion coefficient of ordinary flat glass is approximately (8× -10× ) / °C, while the thermal expansion coefficient of low-expansion glass-ceramics can be as low as close to 0 / °C. By obtaining the specific parameters of the curtain wall material, the accurate value of the material thermal expansion coefficient is determined, and weighted calculation is performed in combination with the direction of the point cloud normal vector , i.e., where is the unit vector in the temperature change direction.

[0025] represents the direction of the point cloud normal vector. Specifically, is the direction vector pointed by the normal vector of the corresponding point in the point cloud data. The normal vector is a vector with magnitude and direction. Here the emphasis is on the property of direction. Through the cross product operation with the unit vector in the temperature change direction and in combination with the material thermal expansion coefficient , the temperature sensitivity coefficient is calculated. Such a calculation method can take into account the relationship between the direction of the point cloud normal vector and the temperature change direction, and thus reflect the difference in the sensitivity of points at different positions to the deformation of the plate under temperature changes.

[0026] After calculating the temperature sensitivity coefficient of each point in each small grid, the point with the largest temperature sensitivity coefficient is selected as the second detection target point to ensure that the selected point can sensitively capture the subtle changes in the deformation gradient of the plate under different temperature difference conditions, thereby effectively monitoring the influence of temperature difference on the flatness of the curtain wall plate.

[0027] In a preferred embodiment of the present invention, in step 2 above, according to the digital twin model, the surface normal vector distribution characteristics are extracted, and the measured wind pressure load distribution data and temperature field parameters are fused to construct a dynamic flexural deformation coupling calculation model, and the global deformation parameters and local deformation parameters are output, which may include: Step 200: According to the surface point cloud data of the digital twin model, perform grid-based block division on the curtain wall surface, and for the point cloud data within each divided block area, extract each target point within the current block and the neighboring point set within a preset radius range. Step 201: Fit the spatial plane of the neighboring point set by the least squares method, calculate the surface normal vector of each target point, and aggregate the normal vectors of each block to generate the surface normal vector distribution characteristics of the curtain wall. Step 202: Generate the measured wind pressure load distribution data by mapping the wind pressure sensor data around the building curtain wall to the corresponding spatial coordinates of the digital twin model; according to the data collected by the temperature sensor and the building structure heat transfer parameters, generate the temperature field parameters through finite element heat transfer analysis, including the temperature spatial distribution matrix and the time-varying gradient sequence. Step 203: Integrate the surface normal vector distribution characteristics of the curtain wall, the measured wind pressure load distribution data, and the temperature field parameters to construct a dynamic flexural deformation coupling calculation model. Step 204: Input the wind pressure load distribution data, the temperature field parameters, and the material elastic modulus into the dynamic flexural deformation coupling calculation model, and set the fixed constraint conditions for the connection nodes between the curtain wall and the main structure. Step 205: Perform dynamic flexural deformation analysis according to the fixed constraint conditions, the wind pressure load distribution data, the temperature field parameters, and the material elastic modulus, and extract the global deformation parameters and the local deformation parameters; the global deformation parameters include the overall waviness and the regional curvature anomaly value; the local deformation parameters include the corner stress gradient of the first detection target point and the peak value of the dynamic flexural displacement at the edge-center junction of the second detection target point.

[0028] In the embodiment of the present invention, based on the surface point cloud data of the digital twin model, first determine the overall range and size of the curtain wall surface. According to the preset grid resolution (for example, divide the curtain wall surface into square grids with a side length of 0.5 meters), use a spatial grid algorithm (such as octree division or uniform grid division) to divide the curtain wall surface into multiple regular divided block areas. Each divided block area forms an independent calculation unit. For the point cloud data within each divided block area, select target points according to certain rules (it can be all points within the divided block area, or select a representative point at a certain interval). For each target point, centered on this point, according to the preset radius (such as 0.3 meters), use a spatial search algorithm (such as KD-Tree search) to find all neighboring points within the radius range in the divided block area, so as to construct the neighboring point set of each target point. This neighboring point set contains the point cloud information of the local area around the target point.

[0029] Step 201: For each target point and its corresponding neighboring point set, the three-dimensional coordinates of the neighboring points ( 、 、 1) As input data, the least squares method is used for spatial plane fitting. The goal of the least squares method is to find a plane equation + + + = 0, such that the sum of the squares of the distances from the neighboring points to the plane is minimized. By solving the corresponding linear equations (matrix equations constructed based on the coordinates of the neighboring points), the coefficients 、 、 、 of the plane equation can be obtained. According to the coefficients of the plane equation, the normal vector of the plane is calculated. The normal vector of the plane can be expressed as ([[]] ), and it is normalized (i.e., divided by the modulus length ) of the normal vector to obtain the unit normal vector, which is the surface normal vector of the target point. The above operations are repeated for all target points within the block region to obtain the surface normal vectors of each target point.

[0030] Aggregate the surface normal vectors of all target points within the block region. Statistical methods, such as calculating the average of the normal vectors, can be used to obtain the representative normal vector of the block region. After performing the same processing on all block regions, the representative normal vectors of each block are integrated to generate the normal vector distribution characteristics of the entire curtain wall surface, which reflect the variation of the normal directions at different positions on the curtain wall surface.

[0031] Step 202, according to the structural characteristics and design requirements of the building curtain wall, reasonably arrange wind pressure sensors at key positions such as the windward side, corners, top, and bottom of the curtain wall. For example, for a high-rise curtain wall building with a rectangular plane, sensors can be evenly arranged at the four corners and the middle area of each floor to ensure coverage of different heights and positions of the curtain wall. These sensors collect the magnitude and direction data of the wind pressure in real time at a fixed sampling frequency (such as once per second) and transmit the data to the data processing center through wireless communication. The data collected by each sensor includes a timestamp, the wind pressure value (unit: Pascal, Pa), the wind pressure direction (expressed in degrees, such as 0° representing the due east direction), and the unique identification ID of the sensor itself.

[0032] Obtain the actual spatial coordinates of each wind pressure sensor, which can be determined by high-precision GPS positioning or measurement based on the building design drawings, using a three-dimensional Cartesian coordinate system ( 、 、 is represented by the axis). At the same time, in the digital twin model, there is also a corresponding coordinate system. To match the sensor coordinates with the digital twin model coordinates, it is necessary to determine the coordinate transformation relationship between the two. If the coordinate systems of the two are inconsistent, operations such as translation, rotation, and scaling may be required. For example, through measurement, it is found that the origin of the coordinate system of the digital twin model is offset by 5 meters in the axis direction from the origin of the actual building coordinate system, axis direction by 3 meters, and there is a rotation angle of 15°. Then, the sensor coordinates need to be adjusted according to the corresponding mathematical transformation formula to make them consistent with the digital twin model coordinates.

[0033] For each point or block area in the digital twin model (assuming the block area is a square with a side length of 0.5 meters), determine the set of wind pressure sensors within a certain range (such as a radius of 5 meters) around it. For each sensor in the set, calculate the Euclidean distance from the sensor to the target point . Then, according to the inverse distance weighted interpolation principle, calculate the wind pressure value at the target point . Suppose there are sensors, and the wind pressure value of each sensor is , then The calculation formula is: = , where is the weight exponent, taking 2, is the index, refers to the Euclidean distance from the th wind pressure sensor to the target point (i.e., the point in the digital twin model where the wind pressure value needs to be calculated). In this formula, the closer the sensor is to the target point, the greater the weight of its wind pressure value in calculating the wind pressure at the target point. For the wind pressure direction, a similar weighted average method is also used for calculation to obtain the wind pressure direction at the target point. After calculating the wind pressure values for all points or block areas, verify the generated measured wind pressure load distribution data. Some known positions that did not participate in the interpolation calculation (such as positions where additional high-precision wind pressure measurement equipment is installed) can be selected, and the calculated wind pressure values are compared with the actual measured values to calculate the error. If the error exceeds the preset threshold (such as 5%), analyze the reasons, which may be unreasonable sensor layout, inaccurate coordinate transformation, or improper parameter settings of the interpolation algorithm, etc. According to the analysis results, adjust the sensor layout, recheck the coordinate transformation relationship, or optimize the interpolation algorithm parameters (such as adjusting the weight exponent ), and then perform the interpolation calculation again until the accuracy requirements are met. Finally, generate accurate and reliable measured wind pressure load distribution data for the entire curtain wall surface, which can describe in detail the size and direction distribution of wind pressure at different positions on the curtain wall surface.

[0034] Temperature data is collected using temperature sensors distributed on the curtain wall surface and inside the building, and at the same time, the heat transfer parameters of the building structure (such as the thermal conductivity, specific heat capacity, heat convection coefficient of materials, etc.) are obtained. Taking these data as inputs, a heat transfer analysis model is established using finite element heat transfer analysis software (such as ANSYS, ABAQUS, etc.). In the model, the mesh is divided according to the geometric shape and material properties of the building structure, and the boundary conditions (such as ambient temperature, solar radiation, etc.) are set. By solving the heat transfer equation (such as Fourier's heat conduction equation), finite element calculations are performed to obtain the temperature field parameters. The temperature field parameters include the temperature spatial distribution matrix (describing the temperature values at different positions on the curtain wall surface at a certain moment) and the time-varying gradient sequence (describing the change trend and gradient of the curtain wall surface temperature over time).

[0035] Step 203: Integrate the normal vector distribution characteristics of the curtain wall surface generated in Step 201, the measured wind pressure load distribution data generated in Step 202, and the temperature field parameters. During the integration process, first, these data are subjected to format conversion and normalization processing to make them have a unified data format and dimension. Then, based on mechanical principles and numerical calculation methods (such as the finite element method), a mathematical model is established to describe the interaction relationship between wind pressure, temperature, and the curtain wall structure. In the model, the normal vector distribution characteristics are used to determine the force direction on the curtain wall surface, the wind pressure load distribution data is used as an external load input, and the temperature field parameters are used to consider the thermal stress caused by temperature changes. By integrating these factors, a dynamic flexural deformation coupling calculation model is constructed, which can simulate the dynamic flexural deformation process of the curtain wall under the combined action of wind pressure and temperature.

[0036] Step 204: Input the measured wind pressure load distribution data, the calculated temperature field parameters, and the elastic modulus of the curtain wall material into the dynamic flexural deformation coupling calculation model. The wind pressure load distribution data is used to apply the external wind load, the temperature field parameters are used to consider the influence of temperature changes on the curtain wall structure, and the material elastic modulus reflects the mechanical properties of the curtain wall material. At the same time, according to the actual connection method between the curtain wall and the main structure, fixed constraint conditions for the connection nodes are set in the model. For example, if the curtain wall is connected to the main structure by bolts, fixed displacement constraints (i.e., restricting the displacement of the node in three directions) are applied at the corresponding connection node positions in the model to accurately simulate the stress state of the curtain wall under actual working conditions.

[0037] Step 205: Define the fixed constraint conditions, which are determined based on the actual connection method between the curtain wall and the building main body. For example, at the bolt-connected parts, the displacement limitations of this position in space need to be set in the analysis to simulate its fixed characteristics; for the hinged connection nodes, specific directions are set to be rotatable, and the displacements in other directions are restricted.

[0038] The wind pressure load distribution data is collected in real time by wind pressure sensors installed on the curtain wall surface, recording the wind pressure magnitude and direction information at different times and positions. The temperature field parameters are obtained through laser point cloud scanning equipment, which can accurately present the spatial distribution of the curtain wall surface temperature and its variation over time. The material elastic modulus is determined based on the material properties of the materials used in the curtain wall (such as glass and metal frames) with reference to material standards, which reflects the material's ability to resist elastic deformation.

[0039] Obtaining the overall waviness: Determine the grid specifications according to the actual size of the curtain wall and the requirements of the detection accuracy. For example, for a curtain wall that is 50 meters long and 30 meters high, if it is selected to be divided into a 10×10 grid, then the size of each grid in the length direction is 5 meters, and in the height direction is 3 meters. During the division process, take the lower left corner vertex of the curtain wall as the coordinate origin (0 0 0), with the horizontal rightward direction as the axis positive direction, the vertical upward direction as the axis positive direction, and the direction perpendicular to the curtain wall surface and outward as the axis positive direction to establish a three-dimensional coordinate system. By equally spacing the coordinate intervals of the axis and the axis, determine the boundary positions of each grid to ensure that the grid evenly covers the entire curtain wall surface. Use a three-dimensional laser scanner to scan the curtain wall in all directions. The scanner emits laser beams, and after the laser encounters the curtain wall surface, it is reflected back to the instrument. By measuring the round-trip time of the laser, calculate the distance from each point on the curtain wall surface to the scanner. Combining the spatial position and attitude information of the scanner, determine the three-dimensional coordinates ( ) of each measurement point, thus forming point cloud data containing a large number of points. During the scanning process, to ensure data integrity, it is necessary to control the scanning angle and movement path of the scanner to ensure that all parts of the curtain wall, including corners and concave-convex areas, can be scanned. After the scanning is completed, classify the point cloud data according to the pre-divided grid areas, extract the point cloud data within each grid, and form independent point cloud data sets.

[0040] For the point cloud data set of each grid area, use the least squares method for plane fitting. The core idea of the least squares method is to find a plane such that the sum of the squares of the perpendicular distances from all the point clouds within the grid to this plane is the smallest. During the specific operation, first assign initial values to the coefficients in the plane equation + + + = 0. Substitute , , Set it as a smaller non-zero random number, such as taking values between -1 and 1, which can avoid special situations in the initial plane (such as being parallel to the coordinate axes), and at the same time provide a reasonable starting point for subsequent optimization; Initially set to 0 to simplify the initial calculation. By continuously adjusting the coefficients , , , values, calculate the perpendicular distance from each point cloud to the plane, and sum the squares of the perpendicular distances of all point clouds. During the process of adjusting the coefficients, use mathematical optimization algorithms (such as the gradient descent method), and according to the change trend of the objective function (the sum of the squares of the perpendicular distances), gradually change the coefficient values according to certain rules (such as the gradient direction and step size). For example, when using the gradient descent method, the amount of each coefficient update is related to the partial derivative of the objective function with respect to the coefficient and the set step size. The step size is a key parameter. If the value is too large, it may cause excessive adjustment of the coefficients and miss the final solution; if the value is too small, the calculation convergence speed will be too slow. In actual operation, the step size will first be set with an empirical value (such as 0.01) and dynamically adjusted according to the change of the objective function during the calculation process.

[0041] Continuously repeat the process of adjusting the coefficients and calculating the objective function value to find the , , , combination that makes this sum of squares reach the minimum value. When the stopping condition is met (such as reaching the preset maximum number of iterations), the plane determined by the coefficient combination at this time is the plane representing the benchmark of the grid area.

[0042] After determining the fitting plane, calculate the perpendicular distance from each point cloud in the grid to this plane. The specific method is to substitute the three-dimensional coordinates of the point cloud ( ) into the distance formula from a point to a plane to obtain the deviation value of each point cloud relative to the fitting plane. These deviation values are positive and negative. A positive value indicates that the point cloud is above the plane, and a negative value indicates that the point cloud is below the plane. The absolute value of the deviation reflects the degree of deviation of the point cloud from the plane. After calculating the deviation values of all point clouds, analyze these deviation values. First, calculate the average value of the deviation values to understand the overall offset trend of the point clouds in the grid area relative to the fitting plane; then calculate the standard deviation of the deviation values to measure the degree of dispersion of the point cloud data relative to the average value. The larger the standard deviation, the more dispersed the point clouds are in the grid area and the more uneven the surface is.

[0043] Determination of regional curvature anomaly values: After collecting the point cloud data of the curtain wall surface using a three-dimensional laser scanner, due to factors such as measurement errors and environmental interference, there may be noise points and outliers in the data. First, identify and remove noise points through statistical analysis methods. For example, calculate the distance between each point and its neighboring points. If the average distance of a certain point to its neighboring points is significantly greater than that of other points, it is determined as a noise point. For outliers, a density-based clustering algorithm can be used to regard the points in the low-density area as outliers and remove them. Before performing NURBS surface fitting, it is necessary to determine the number and distribution of the control points of the surface. According to the complexity and accuracy requirements of the local area of the curtain wall, reasonably arrange the control point grid within the bounding box of the point cloud data. For areas with relatively complex shapes, such as the corners or decorative protrusions of the curtain wall, increase the density of the control points; for relatively flat areas, appropriately reduce the number of control points. Determine the weight factors of the NURBS surface based on the point cloud data and the set control points. The weight factors affect the degree of fitting between the surface and the point cloud data. Empirical pre-judgment assignment is used to assign different initial values to the weight factors of the control points in different regions according to the point cloud density and the general shape of the surface. For example, at the corners where the point cloud is dense, set the initial value of the weight factor to 1.2 - 1.5, so that the surface is more inclined to approach the point cloud data in this area during the initial construction; in relatively flat areas, set the initial value of the weight factor to 0.8 - 1 to avoid shape distortion caused by excessive fitting of the surface.

[0044] Taking the least squares method as an example for optimization, the sum of the squares of the distances from the point cloud to the surface is used as the objective function. During the optimization process, adjust the weight factors according to the change trend of the objective function. Specifically, when adjusting, change the weight factors at a certain step size. The initial value of the step size is set to 0.1, which can not only ensure obvious changes in the weight factors but also avoid missing the optimal solution due to too large an adjustment amplitude. During the calculation process, monitor the change of the objective function value in real time. If the objective function value drops significantly, it indicates that the current step size and adjustment direction are appropriate, and continue to adjust at this step size; if the objective function value drops slowly, then reduce the step size to 0.05 or even smaller to ensure the fineness of the adjustment; if the objective function value rises, change the adjustment direction and appropriately reduce the step size to readjust. Set double termination conditions. One is to set the maximum number of iterations, such as 200300 times. When the number of iterations reaches this upper limit, stop the calculation regardless of whether the objective function value converges; the other is to set the change threshold of the objective function value, such as 0.001. When the difference between the objective function values calculated in two adjacent iterations < this threshold, it is considered that the surface has sufficiently approximated the point cloud data. The finally obtained combination of weight factors is the value that can make the NURBS surface fit the point cloud data as much as possible and is used to construct an accurate local surface model of the curtain wall.

[0045] During the fitting process, it is also necessary to determine the degree of the B-spline basis function. The higher the degree of the basis function, the smoother the surface, but the higher the computational complexity. Select an appropriate degree according to actual needs. For relatively smooth surfaces such as curtain wall surfaces, choose cubic B-spline basis functions, which can not only ensure the smoothness of the surface but also effectively control the amount of calculation. By adjusting the control points, weight factors, and the degree of the basis function multiple times and performing fitting repeatedly, a NURBS surface model that can accurately reflect the local surface shape of the curtain wall is finally obtained. For the generated NURBS surface model, a geometric analysis method is used to calculate the curvature values of each point on the surface. First, determine the parametric coordinates of the points on the surface where the curvature is to be calculated. By evenly dividing the grid on the surface, a series of points are selected as the calculation objects. At each selected point, calculate the first-order derivative and the second-order derivative of the surface. The first-order derivative reflects the tangential direction of the surface at that point, and the second-order derivative is related to the degree of curvature of the surface. By performing derivative operations on the parametric equation of the NURBS surface, these derivative information are obtained.

[0046] Based on the first-order derivative and the second-order derivative, further calculate the normal vector of the surface at that point. The normal vector is perpendicular to the surface. Let the first-order partial derivatives of the surface be and , and the second-order partial derivatives be 、 and , then the unit normal vector of the surface is . Using the normal vector and the derivative information, calculate the mean curvature of this point, where , , , is the dot product of the first-order partial derivative of the surface in the parametric direction with itself, is the dot product of the first-order partial derivative of the surface in the parametric direction and the first-order partial derivative in the parametric direction, is the dot product of the first-order partial derivative of the surface in the parametric direction with itself. Through this formula, the curvature values of each point on the NURBS surface can be obtained, thereby evaluating the degree of curvature at different positions of the surface.

[0047] According to the design drawings of the curtain wall and the mechanical properties of the materials, combined with engineering experience and relevant standards and specifications, a reasonable curvature range is preset as the judgment criterion. Taking a 6-mm thick tempered glass flat curtain wall as an example, based on its elastic modulus of 72 GPa and the allowable The span deformation amount, verified by finite element simulation and engineering practice, sets the normal curvature range to 0.001 - 0.003 . When comparing the curvature values of each point on the calculated surface, a threshold needs to be set to accurately identify outliers. Considering the error of the measuring device ±0.3 mm (such as the accuracy of the laser scanner), the fluctuation range of the elastic modulus of the glass material ±5%, and the process error during construction and installation, the threshold is set to 20% of the upper limit of the standard range. Specifically, when the curvature value of a certain point exceeds 0.0036 (i.e., 0.003 ×1.2), it is determined as an outlier of the regional curvature. This threshold setting mechanism can effectively avoid misjudgment caused by measurement errors or minor fluctuations in material properties. Once an outlier appears, it indicates that there is bending deformation beyond the normal range in the corresponding local area of the curtain wall, which may be due to fatigue damage caused by long-term wind loads on the structure or deviation from the design accuracy requirements during installation. Immediate targeted detailed inspections and mechanical analyses are required.

[0048] Calculation of the angular stress gradient of the first detection target point: According to the stress distribution results obtained from the finite element calculation model, locate the first detection target point and its adjacent points in the corner area of the curtain wall panel. The selection of adjacent points can be based on actual needs. Taking the target point as the center, set a spherical or cubic region, and extract the stress values of all nodes within this region. The stress values include normal stress and shear stress. The finite difference method is used to calculate the partial derivatives of the stress in , , directions. Taking the direction as an example, for the target point and its adjacent points in the and directions, the partial derivative of the stress in the direction is approximately , where is the distance between the two points in the direction. Similarly, the partial derivatives in the and directions are and . Through the partial derivatives calculated in the three directions, according to vector operations, the angular stress gradient vector is obtained. The magnitude of the stress gradient reflects the degree of stress concentration in the corner area; the direction of the stress gradient indicates the trend of stress change.

[0049] Calculation of the peak value of the dynamic flexural displacement at the edge - center junction of the second detection target point: Locate the second detection target point at the junction of the center and the edge of the curtain wall panel from the database of finite element calculation results. In each calculation time step, extract the displacement component data of this target point in 、 、 three directions, and simultaneously obtain the corresponding timestamps. For example, if the finite element calculation uses 0.1 second as a time step, at the 10th time step, it is obtained that the displacement of the target point in direction is 2 mm, direction is -1.5 mm, direction is 0.8 mm, and the timestamp is 1 second. Organize this kind of data into a multi-dimensional array form in chronological order for storage. Use the linear search algorithm to traverse the stored displacement data array. Set the initial maximum displacement value as the combined displacement value of the first time step. The combined displacement calculation formula is where 、 、 are the three-direction displacement components of the first time step. Starting from the second time step, calculate the combined displacement for each time step and compare it with the current maximum displacement value. After traversing the entire array, the finally obtained maximum combined displacement value is the peak value of the dynamic deflection displacement, and the corresponding time point can reflect the moment when the maximum deformation occurs, so as to evaluate the maximum deformation situation generated by the temperature difference and wind force at the junction.

[0050] Suppose there is a high-rise office building with a height of 100 meters and 30 floors, and its exterior facade uses a glass curtain wall structure. When detecting the flatness of the curtain wall: Divide the point cloud data on the curtain wall surface in the digital twin model into grids of 1 meter × 1 meter, and a total of 5000 sub-block areas are obtained. For the point cloud data in each sub-block area, select the point at the center of the sub-block as the target point, and use a preset radius of 0.2 meter to find the adjacent point set with the target point as the center. For example, in a certain sub-block area, 15 adjacent points are found around the target point. For this target point and its adjacent point set, use the least squares method to fit the spatial plane to obtain the plane equation 2x + 3y - 4z + 5 = 0, and calculate and normalize to obtain the surface normal vector of the target point as (0.37, 0.55, -0.74). After performing the same processing on the target points of all sub-block areas, aggregate to obtain the distribution characteristics of the curtain wall surface normal vector.

[0051] A total of 20 wind pressure sensors were installed on the four facades of the building curtain wall. The data collected by the sensors were mapped onto the digital twin model through inverse distance weighted interpolation to generate the measured wind pressure load distribution data. At the same time, data were collected using 10 temperature sensors installed on the curtain wall surface and inside the building. Combining the heat transfer parameters of the building structure, finite element heat transfer analysis was performed through ANSYS to obtain the temperature field parameters. For example, the temperature spatial distribution matrix at a certain moment showed that the temperature at the top of the curtain wall was 28°C and the temperature at the bottom was 25°C, and a time-varying gradient sequence of temperature changing with time was generated. The normal vector distribution characteristics, wind pressure load distribution data, and temperature field parameters were fused, and a dynamic flexural deformation coupling calculation model was constructed based on the finite element method. The wind pressure load distribution data, temperature field parameters, and elastic modulus of the glass material (assumed to be 70 GPa) were input into the model, and the fixed constraint conditions of the connection nodes between the curtain wall and the main structure were set. Through model calculation, global deformation parameters were obtained, such as the overall waviness was 3 mm and there were 3 regional curvature abnormal values; local deformation parameters, such as the corner stress gradient of the first detection target point was 15 MPa / m, and the peak value of the dynamic flexural displacement at the edge-center junction of the second detection target point was 8 mm.

[0052] Through grid-based block processing and extraction of neighboring point sets, the point cloud data on the curtain wall surface was refined. By fitting a plane using the least squares method to calculate the surface normal vector, the geometric feature information of the curtain wall surface could be accurately obtained, improving the accuracy and reliability of data processing compared with traditional methods. The fusion of the measured wind pressure load distribution data and temperature field parameters comprehensively considered the influence of environmental factors on the curtain wall. The real-time collection and accurate mapping of the wind pressure sensor data and temperature sensor data, combined with the generation of temperature field parameters through finite element heat transfer analysis, enabled the model to truly simulate the stress and deformation of the curtain wall under actual working conditions, avoiding the limitations of single-parameter analysis. The constructed dynamic flexural deformation coupling calculation model could comprehensively consider the interaction relationship between wind pressure, temperature, and the curtain wall structure. By setting reasonable constraint conditions and input parameters, accurate simulation of the dynamic flexural deformation of the curtain wall was achieved. The extracted global deformation parameters and local deformation parameters comprehensively described the deformation of the curtain wall from both the overall and local levels. The overall waviness and regional curvature abnormal values reflected the overall flatness and local abnormal deformation of the curtain wall; the corner stress gradient and the peak value of the dynamic flexural displacement at the edge-center junction were analyzed in detail for key parts, helping to detect potential safety hazards in a timely manner and providing specific quantitative indicators for the maintenance and management of the curtain wall.

[0053] In a preferred embodiment of the present invention, in step 3 above, according to the local deformation parameters, through phase correlation analysis and wavelet transform, the vibration displacement components of the first detection target point and the second detection target point are extracted to construct a curtain wall vibration mode feature matrix. At the same time, the displacement-time series data is analyzed to identify over-limit displacement events where the displacement amplitude is greater than or equal to a preset threshold, which may include: Step 300, according to the local deformation parameters, determine the dynamically deformation-sensitive regions corresponding to the first detection target point and the second detection target point in the time-series point cloud data; Step 301, for the point cloud sequence of the corner region of the first detection target point in the sensitive region, through phase correlation analysis, extract the displacement change amount between adjacent time frames, and use wavelet transform to separate the low-frequency thermal expansion displacement component and the high-frequency wind-induced vibration displacement component; the high-frequency wind-induced vibration displacement component is used as the corner vibration displacement component; for the point cloud sequence of the junction region of the second detection target point, through wavelet transform, extract the dynamic displacement component within a preset frequency band; Step 302, spatially and temporally align the corner vibration displacement component of the first detection target point with the dynamic displacement component of the junction of the second detection target point to construct a curtain wall vibration mode feature matrix; Step 303, analyze the displacement-time series data in the vibration mode feature matrix, and preset a displacement amplitude threshold to identify single over-limit events where the displacement amplitude is greater than or equal to the preset displacement amplitude threshold.

[0054] In an embodiment of the present invention, centered on the first detection target point, according to the distribution range of the corner stress gradient, a spatial range including the stress concentration region is delimited. For example, if the corner stress gradient changes significantly within a range of 0.5 meters from the target point, a spherical region with the target point as the center and 0.5 meters as the radius is delimited as the dynamically deformation-sensitive region of the first detection target point. For the second detection target point, based on the influence range of the peak value of its dynamic flexural displacement, a rectangular or irregular-shaped region is delimited at the junction of the center and the edge. Assuming that the peak displacement at the junction where the second detection target point is located mainly affects the surrounding area of 1 meter × 0.8 meters, a rectangular region with corresponding dimensions is delimited. In the time-series point cloud data, the point cloud data sequence falling within the sensitive region is screened by matching the coordinate information of the point cloud with the delimited sensitive region. Since the point cloud data has time-stamp information, these screened point cloud sequences can reflect the deformation of the sensitive region over time.

[0055] Step 301, for the point cloud sequence of the corner region, regard the point cloud data of adjacent time frames as an image (which can be converted into a two-dimensional image by projecting the point cloud onto a specific plane). Phase correlation analysis is based on Fourier transform. First, perform two-dimensional Fourier transform on the images corresponding to the point cloud of two adjacent frames to obtain the amplitude and phase information in the frequency domain. According to the phase correlation theorem, calculate the phase correlation function of the two frames of images Among them, and are the Fourier transform results of two frames of images respectively, is 's complex conjugate, is the frequency coordinate of the corresponding image in the horizontal direction, is the frequency coordinate of the corresponding image in the vertical direction. By finding the peak position of the phase correlation function, the translation amount between two frames of point clouds is determined, that is, the displacement change amount between adjacent time frames. Taking the extracted displacement change amount sequence as the input, a suitable wavelet basis function (such as Daubechies wavelet) is selected for discrete wavelet transform. The wavelet transform decomposes the displacement sequence into different frequency sub-bands. The displacement changes caused by thermal expansion usually have a lower frequency, while the displacement changes caused by wind vibration have a higher frequency. A frequency threshold is set based on historical data statistical analysis. For example, through the spectral analysis of past curtain wall displacement data, it is found that the thermal expansion displacement frequency is mainly between 0 - 0.1 Hz, and the wind vibration displacement frequency is above 0.1 Hz, then the frequency threshold can be set to 0.1 Hz. In the discrete wavelet transform results, the detail components of different layers correspond to different frequency ranges. The higher the layer number, the lower the corresponding frequency. According to the frequency threshold and the frequency characteristics of the wavelet transform, it is determined which detail components and approximation components belong to the thermal expansion displacement component and which belong to the corner vibration displacement component. According to the set frequency threshold, the approximation component and detail components belonging to the thermal expansion displacement component are reconstructed to obtain the thermal expansion displacement component sequence; the detail components belonging to the corner vibration displacement component are reconstructed to obtain the corner vibration displacement component sequence.

[0056] Directly perform wavelet transform on the point cloud sequence in the junction area. According to the frequency range affected by factors such as wind force under the actual working conditions of the curtain wall, a preset frequency band is set (for example, 3 - 10 Hz, which is the common frequency range for general building curtain walls affected by wind vibration). After decomposing the point cloud sequence into different frequency sub-bands through wavelet transform, the dynamic displacement component within the preset frequency band is extracted, and this component contains the displacement change information generated by the dynamic loads such as wind force acting on the junction.

[0057] Step 302. Since there may be slight differences in the displacement data acquisition times of the first and second detection target points and their spatial positions are different, spatio-temporal alignment is required. In terms of time, using the time series of one of the target points as a reference, the displacement data of the other target point is time-calibrated by methods such as linear interpolation to ensure that the displacement data of the two target points corresponds to the same moment on the time axis. In terms of space, according to the three-dimensional coordinates of the two target points in the digital twin model, the displacement data is transformed into a unified coordinate system. The corner vibration displacement components of the first detection target point and the dynamic displacement components at the junction of the second detection target point after spatio-temporal alignment are arranged in chronological order to form a two-dimensional matrix. Each row of the matrix corresponds to a time step, and each column corresponds to the displacement components of the two target points in different directions ( 、 、 ). For example, the first column of the matrix is the vibration displacement component of the first detection target point in the direction, and the second column is the dynamic displacement component of the second detection target point in the direction, and so on. Finally, the vibration mode feature matrix of the curtain wall is constructed.

[0058] Step 303. According to the curtain wall design specifications and safety standards, combined with historical detection data and actual usage conditions, a displacement amplitude threshold is preset. For example, for a certain type of curtain wall, according to its material strength and structural design, the displacement amplitude threshold is set to 5 mm. Analyze the displacement data of each time step in the vibration mode feature matrix, and calculate the combined displacement amplitude of the displacement components in each direction for each time step: , where 、 、 are the displacement components in three directions respectively. If the combined displacement amplitude of a certain time step is greater than or equal to the preset threshold, then this time step is identified as an over-limit displacement event, and information such as the time of the event occurrence and the displacement amplitude is recorded. By analyzing the entire displacement-time series data, the occurrence times, durations and other characteristics of the over-limit displacement events can be statistically analyzed, providing a basis for evaluating the safety of the curtain wall.

[0059] Suppose there is a high-rise commercial building with a unitized glass curtain wall structure on its exterior facade. In a certain detection: Through preliminary calculations, it is found that the corner stress gradient of the first detection target point changes significantly within a region centered at the target point with a radius of 0.4 meters. Therefore, this spherical region is designated as the dynamic deformation sensitive region of the first detection target point; the second detection target point is at the junction of the center and the edge, and the peak influence range of its dynamic flexural displacement is a rectangular region of 1.2 meters × 0.6 meters. This rectangular region is designated as the dynamic deformation sensitive region of the second detection target point. In the time-series point cloud data, the point cloud sequences falling into these two regions are screened to obtain the data for subsequent analysis. For the first detection target point, the point cloud sequences of the corner region in adjacent time frames are projected onto the XY plane and transformed into images for phase correlation analysis, and the displacement change amount between adjacent two frames of point clouds is calculated. Then, discrete wavelet transform is performed using Daubechies wavelets, with the frequency threshold set at 1 Hz. The displacement components below 1 Hz are regarded as thermal expansion displacement components, and those above 1 Hz are extracted as corner vibration displacement components. For the second detection target point, wavelet transform is performed on the point cloud sequence of the junction region, and the dynamic displacement components within the frequency band of 3 - 8 Hz are extracted. Based on the time series of the first detection target point, time interpolation calibration is performed on the displacement data of the second detection target point, and the displacement data of the two target points are converted to the same coordinate system. Then, the corner vibration displacement components of the first detection target point and the dynamic displacement components at the junction of the second detection target point are arranged in chronological order to construct a two-dimensional matrix, completing the construction of the vibration mode characteristic matrix of the curtain wall.

[0060] According to the design requirements of this curtain wall, the preset displacement amplitude threshold is 4 mm. The displacement data in the vibration mode characteristic matrix are analyzed. At a certain time step, the comprehensive displacement amplitude of the first detection target point is calculated to be 5.2 mm, exceeding the threshold. This time step is identified as an over-limit displacement event, and relevant information is recorded.

[0061] Through phase correlation analysis and wavelet transform, the vibration displacement components of the curtain wall under the action of different environmental factors can be accurately extracted from complex point cloud data. Phase correlation analysis can accurately calculate the displacement changes between adjacent time frames, and wavelet transform effectively separates low-frequency thermal expansion and high-frequency wind-induced vibration displacements, avoiding the problem that it is difficult to distinguish the influences of different factors by traditional methods, providing a powerful means for in-depth analysis of the dynamic response of the curtain wall. After spatially and temporally aligning the displacement components of the first and second detection target points, the vibration mode characteristic matrix is constructed, integrating the dynamic information of key positions of the curtain wall. By presetting the threshold to identify over-limit displacement events, abnormal displacement situations occurring during the use of the curtain wall can be detected in a timely manner. Compared with traditional manual inspections, this data-driven automatic identification method is more sensitive and accurate, can early warn of potential safety risks, provide a scientific basis for the maintenance and repair of the curtain wall, reduce the probability of safety accidents, and ensure the safe operation of the building.

[0062] In a preferred embodiment of the present invention, in step 4, according to the global deformation parameter and the occurrence frequency of the over-limit displacement event, the initial coefficient of the curtain wall structure safety degree is determined, and based on the local deformation parameter and the dynamic displacement component, the target correction value is calculated, which may include: Step 400: Obtain the historical cumulative term of wind pressure according to the wind pressure from the starting moment to the current moment and the cumulative attenuation of wind pressure over time; analyze the temperature field data of the curtain wall surface obtained by laser point cloud scanning, and combine the spatial coordinate information of the curtain wall to obtain the temperature gradient reflecting the local thermal stress condition of the curtain wall; determine the fatigue cumulative damage according to the number of cycles and the corresponding fatigue life of the curtain wall under different stress levels; Step 401: Integrate the historical cumulative term of wind pressure, the gradient of the temperature field at the spatial position, and the fatigue cumulative damage to determine the static-related term; Step 402: Determine the maximum displacement of the curtain wall junction under dynamic load, that is, the peak value of the dynamic deflection displacement, and the maximum vibration displacement of the curtain wall corner in high-frequency wind vibration after band-pass filtering, that is, the peak value of the filtered corner vibration displacement, to obtain a comprehensive parameter reflecting the degree of mutual influence between the deformation and vibration of the curtain wall structure under dynamic load; Step 403: Determine the relationship between the static risk and the deformation-vibration sensitivity of the curtain wall structure according to the static-related term and the comprehensive parameter reflecting the degree of mutual influence between the deformation and vibration of the curtain wall structure under dynamic load; obtain the relationship characterizing the coordination between the vibration displacements of the corner and the junction, that is, the dynamic-related term, according to the peak value of the filtered corner vibration displacement and the peak value of the dynamic displacement of the curtain wall junction after the same band-pass filtering; Step 404: Determine the contributions of the static factors and the dynamic factors to the target correction value according to the static-related term and the dynamic-related term, so as to obtain the final target correction value.

[0063] In the embodiment of the present invention, high-precision wind pressure sensors are installed at key positions such as the windward side, corners, top and bottom of the curtain wall, and wind pressure data is collected in real time at a fixed frequency (such as once per second). The wind pressure data collected by the sensors contains magnitude and direction information, and is transmitted to the data center through a wireless communication method to form a sequence of wind pressure changing with time , where is a time variable. The attenuation coefficient is determined according to the climate conditions of the region where the building is located and the historical wind pressure data, and its value range is between 0.02 and 0.1, which is used to describe the attenuation characteristic of wind pressure over time. The wind pressure farther from the current moment has less influence on the structure.

[0064] The surface of the curtain wall is scanned periodically by a laser point cloud scanning device, and the scanning frequency is set according to the change of the ambient temperature. For example, the scanning is densified when the temperature difference between day and night is large or during seasonal alternation. The scanning device combines a high-precision positioning system to obtain the temperature value of each point and its three-dimensional spatial coordinates The temperature gradient is obtained by performing finite difference calculations on the temperature values of adjacent points , which reflects the local thermal stress condition of the curtain wall. Refer to the curtain wall design documents and material property reports to obtain the theoretical fatigue life of the curtain wall structure under different stress levels Meanwhile, install stress sensors at the key stress-bearing parts of the curtain wall to monitor stress changes in real time. When the stress exceeds the set threshold, record one cycle, and accumulate to obtain the actual number of cycles under each stress level .

[0065] Install high-precision displacement sensors at the intersections and corners of the curtain wall to monitor displacement changes under dynamic loads (such as wind force). By setting a band-pass filter, screen out the displacement data in the high-frequency wind vibration frequency band (3 - 10 Hz), and then determine the maximum displacement at the curtain wall intersection under dynamic loads (i.e., the peak value of dynamic flexural displacement), the maximum vibration displacement of the curtain wall corner after frequency band filtering (i.e., the peak value of the filtered corner vibration displacement), and the peak value of the dynamic displacement at the curtain wall intersection after the same frequency band filtering .

[0066] Comprehensively determine the initial safety factor: Monitor curtain walls of different types (such as unit curtain walls, framed curtain walls, point-supported curtain walls, etc.) and different scales (small building curtain walls, large commercial complex curtain walls, high-rise building curtain walls, etc.). These curtain walls cover different geographical locations, service life, and environmental conditions to obtain rich and diverse data samples

[0067] For each data sample, accurately record the following key information. Use high-precision laser scanning equipment to obtain the point cloud data of the curtain wall surface, and calculate the standard deviation of the overall waviness through processing, accurate to 0.01 mm. Through advanced three-dimensional modeling technology combined with the point cloud data, calculate and count the number of regional curvature outliers, and record the degree of deviation of each outlier from the normal range, accurate to 0.001 . Use high-precision displacement sensors installed at the key parts of the curtain wall to count the occurrence frequency of out-of-limit displacement events. The statistical period can be set to one month. For each data sample, combined with the design data and historical maintenance records of the curtain wall, preliminarily determine the corresponding initial safety factor, with the value range set between 0 and 1. The larger the value, the higher the safety level

[0068] After a large number of data samples are collected, the data is first cleaned and preprocessed. Check the integrity of the data, fill in the missing values, verify, correct or eliminate the abnormal data points. Standardize the data such as the standard deviation of the overall waviness, the number and deviation degree of regional curvature outliers, and the occurrence frequency of over-limit displacement events to make them comparable. Use statistical methods, such as correlation analysis, to study the linear or non-linear relationships between various factors and the initial safety coefficient, and preliminarily judge which factors have a greater impact on the initial safety coefficient. Taking multiple linear regression as an example, use the standard deviation of the overall waviness, the number and deviation degree of regional curvature outliers, and the occurrence frequency of over-limit displacement events as independent variables, and the initial safety coefficient as the dependent variable to construct a regression model. By continuously adjusting the model parameters and optimizing the model fitting effect, make the model accurately reflect the internal relationship between various factors and the initial safety coefficient, so as to construct an accurate correspondence table.

[0069] Determination of the initial safety coefficient: Use the measuring equipment and standard measuring methods to accurately obtain the standard deviation of the overall waviness, the number and deviation degree of regional curvature outliers of the curtain wall again, and accurately count the occurrence frequency of over-limit displacement events through displacement sensors. Substitute the data obtained from actual measurement and statistics into the correspondence table established through the above analysis. Use the weighted average method for comprehensive calculation to determine the coefficients and weights corresponding to each factor.

[0070] For example, the coefficient corresponding to the standard deviation of the overall waviness is determined by dividing the interval according to its numerical value. When the standard deviation is less than 0.5 mm, the corresponding coefficient is 0.8 - 1.0; when it is between 0.5 - 1.5 mm, the corresponding coefficient is 0.6 - 0.8; when it is greater than 1.5 mm, the corresponding coefficient is 0.4 - 0.6. The correlation coefficient of regional curvature outliers is determined according to the number and deviation degree of outliers. When the number is small and the deviation degree is low, the coefficient is 0.8 - 1.0; when the number and deviation degree are moderate, the coefficient is 0.6 - 0.8; when the number is large and the deviation degree is high, the coefficient is 0.4 - 0.6. The coefficient corresponding to the occurrence frequency of over-limit displacement events is determined according to the occurrence frequency. When the occurrence frequency within one month is less than 5 times, the coefficient is 0.8 - 1.0; when it is between 5 - 15 times, the coefficient is 0.6 - 0.8; when it is greater than 15 times, the coefficient is 0.4 - 0.6.

[0071] In terms of weights, the weight range of the standard deviation of the overall waviness is 0.3 - 0.5. For curtain walls with high flatness requirements, such as glass curtain walls, this weight can be biased towards 0.5; the weight of the relevant data of regional curvature outliers is 0.2 - 0.4. If the curtain wall structure is complex and local stress problems are likely to occur, the weight can be biased towards 0.4; the weight of the occurrence frequency of over-limit displacement events is 0.2 - 0.4. In the curtain wall scenarios with large wind loads or long service life, the weight can be biased towards 0.4.

[0072] Initial safety factor = coefficient corresponding to the standard deviation of overall waviness × its weight + coefficient related to regional curvature outliers × its weight + coefficient corresponding to the occurrence frequency of over-limit displacement events × its weight. Considering the influence of various factors comprehensively, the accurate initial safety factor of the curtain wall structure is finally determined. 。

[0073] Setting of weight coefficients: Used to adjust the weight of the historical cumulative term of wind pressure in the static-related terms, with a value range of 0.3 - 0.5. When the building is located in a windy area, a relatively larger value can be appropriately taken; Used to adjust the weight of the temperature gradient term, with a value range of 0.2 - 0.4. For curtain wall materials with high thermal sensitivity, a relatively larger value is taken; Reflects the importance degree of the fatigue cumulative damage term, with a value range of 0.2 - 0.3; Used to quantify the relative importance of static-related terms (including historical cumulative wind pressure, temperature gradient, fatigue cumulative damage) when calculating the target correction value, with a value range of 0.6 - 0.8, emphasizing the dominant influence of static factors on the safety of the curtain wall structure; Adjusts the weight of dynamic-related terms (vibration displacement coordination), with a value range of 0.2 - 0.4.

[0074] Positive number Take , used to avoid the situation where the denominator of the formula is zero.

[0075] Using numerical integration methods (such as trapezoidal integration method), divide the time interval , into multiple small intervals, perform weighted calculation on the wind pressure values within each small interval , and then multiply by the weight coefficient to obtain the value of the historical cumulative term of wind pressure . Multiply the obtained temperature gradient by the weight coefficient to obtain the value of the temperature gradient term .

[0076] Add the ratios of the actual number of cycles at different stress levels to the corresponding fatigue life , and then multiply by the weight coefficient to obtain the value of the fatigue cumulative damage term .

[0077] Add the wind pressure historical cumulative term, temperature gradient term, and fatigue cumulative damage term calculated above to obtain the numerator of the static-related term. , combined with the obtained dynamic displacement parameters and calculate the denominator; finally, divide the numerator by the numerator of the static-related term to obtain the value of the static-related term .

[0078] Divide the peak value of the filtered corner vibration displacement by the peak value of the filtered dynamic displacement at the intersection of the curtain wall in the same frequency band , and then multiply by the weight coefficient to obtain the value of the dynamic-related term, which characterizes the coordination of the vibration displacements at the corners and intersections of the curtain wall.

[0079] Multiply the calculated static-related term by the weight coefficient , and add it to the dynamic-related term to obtain the final target correction value , that is . This value is used to correct the initial coefficient of the curtain wall structure safety factor , comprehensively reflecting the influence degrees of static factors and dynamic factors on the curtain wall structure safety. The larger the value, the higher the risk faced by the curtain wall structure.

[0080] The formula combines static factors such as wind pressure historical accumulation, temperature gradient, and fatigue cumulative damage, as well as dynamic factors such as the deformation and vibration relationship of different parts of the curtain wall under dynamic loads, changing the limitations of traditional single-factor or few-factor evaluations, and being able to comprehensively reflect the actual situation of the curtain wall structure under complex environments and loads, improving the accuracy and reliability of the evaluation. It includes considerations of static long-term acting factors (such as long-term wind pressure accumulation and fatigue damage), and incorporates the structural responses under dynamic loads (such as peak vibration displacements and coordination), achieving an evaluation method that combines static and dynamic aspects. It can more realistically simulate the stress and deformation conditions of the curtain wall during actual use and timely discover potential safety hazards.

[0081] By setting different weight coefficients ( , , , , ), according to the specific conditions such as the material properties, structural forms, and use environments of the curtain wall, the contribution degrees of various factors to the target correction value can be flexibly adjusted. This makes the evaluation model have stronger adaptability and is applicable to the safety evaluations of different types of curtain wall structures.

[0082] The relationship between the static risk and the deformation-vibration sensitivity of the curtain wall structure, as well as the relationship such as the vibration displacement synergy at the corners and joints, is clearly determined, and the target correction value is obtained through quantitative calculation. It provides specific quantitative indicators for the safety assessment of the curtain wall structure, facilitating engineers and managers to intuitively understand the safety status of the curtain wall structure and providing a scientific basis for decision-making. Through the calculation of the fatigue cumulative damage term, the fatigue state of the curtain wall structure during long-term use can be tracked, the risk of structural damage caused by fatigue can be pre-warned in advance, which helps to formulate reasonable maintenance and inspection plans, extend the service life of the curtain wall, and reduce the probability of safety accidents.

[0083] In a preferred embodiment of the present invention, in step 5 above, the target correction value is superimposed on the initial safety factor to generate a corrected safety factor of the curtain wall structure, realizing the dynamic adjustment of the initial safety factor, and outputting a three-dimensional detection report, including the risk level, the distribution of the corrected safety factor, and the warning marks associated with the target correction value, which may include: Step 500, integrating the dynamic adjustment information reflected by the target correction value, including local deformation and dynamic displacement components, into the initial safety factor to generate a corrected safety factor of the curtain wall structure; Step 501, according to the corrected safety factor of the curtain wall structure, and referring to the preset risk level classification standard, determine the risk level of the curtain wall; Step 502, map the corrected safety factor onto the digital twin model of the curtain wall surface, process the distribution of the corrected safety factor in each part of the curtain wall, and obtain the distribution status of the corrected safety factor in each part of the curtain wall; Step 503, according to the magnitude and change trend of the target correction value, as well as the association with the distribution status of the safety factor in each part of the curtain wall, add corresponding warning marks in the digital twin model; Step 504, integrate the risk level of the curtain wall, the distribution status of the safety factor in each part of the curtain wall, the warning marks, the basic information of the curtain wall, and the description of the detection process, including the detection time, detection equipment, and detection method, to generate a three-dimensional detection report.

[0084] In the embodiment of the present invention, first, the target correction value is analyzed to clarify the information of local deformation and dynamic displacement components it contains. For example, the target correction value may contain dynamic displacement information such as the peak value of high-frequency vibration displacement at the corners of the curtain wall and the peak value of dynamic flexural displacement at the joints, as well as information related to local deformation such as temperature gradient and fatigue cumulative damage. These information are classified and quantified according to certain rules. Determine the fusion method of the initial safety factor and the target correction value. The fusion strategy is weighted summation, that is, the corrected safety factor of the curtain wall structure + × , where is the initial safety factor, is the fusion weight, and its value range is between 0.5 and 1, which can be adjusted according to actual engineering experience and the degree of emphasis on various factors. This formula means that on the basis of the initial safety factor, the coefficient is dynamically adjusted according to the magnitude of the target correction value and the fusion weight. If the target correction value reflects a higher risk, after calculation by this formula, the corrected safety factor will increase accordingly, and vice versa. For each detection point (or grid unit) of the curtain wall structure, the initial safety factor and the corresponding target correction value are obtained respectively, and the corrected safety factor of the curtain wall structure at each point (or grid unit) is calculated according to the above fusion formula.

[0085] Step 501, before the detection, according to the design specifications, material properties, service life of the curtain wall and relevant industry standards, formulate a preset risk level classification standard. For example, the risk level is divided into four levels: low risk, medium risk, high risk and extremely high risk, and a specific safety factor range is set for each level. Suppose the safety factor range corresponding to low risk is 0.8 - 1.0, medium risk is 0.6 - 0.79, high risk is 0.4 - 0.59, and extremely high risk is 0 - 0.39.

[0086] Compare the corrected safety factor of the curtain wall structure at each detection point (or grid unit) calculated in Step 500 with the preset risk level classification standard. Determine the risk level to which each point (or grid unit) belongs. Then, by counting the number of points (or grid units) in different risk levels and the proportion they account for in the total area of the curtain wall, comprehensively evaluate the risk level of the entire curtain wall. For example, if the area proportion of high-risk and extremely high-risk areas exceeds 30%, it is determined that the risk level of the entire curtain wall is high risk.

[0087] Step 502, perform grid processing on the digital twin model of the curtain wall, and each grid corresponds to an actual area on the curtain wall surface. Map the corrected safety factor of each detection point (or grid unit) calculated in Step 500 to the corresponding grid in the digital twin model. Ensure that the safety factor corresponds accurately to the actual curtain wall part.

[0088] Render the digital twin model using professional 3D visualization software (such as Unity, Unreal Engine, etc.). Assign different colors or transparencies to each grid according to the magnitude of the corrected safety factor. For example, areas with a high safety factor are represented in green, and as the factor decreases, the color gradually transitions to yellow, orange, and finally red for high-risk areas. At the same time, different transparencies can be set to make the distribution of the safety factor more intuitive. In addition, contour lines or isopleths can be added to further display the changing trend of the safety factor on the curtain wall surface. Smooth the mapped safety factor data to avoid discontinuous or abnormal displays caused by detection errors or data fluctuations. Filtering algorithms (such as Gaussian filtering) can be used to process the data to make the distribution of the safety factor smoother and more natural. At the same time, compress and optimize the data to improve the efficiency and performance of the visual display.

[0089] Step 503: Conduct an in-depth analysis of the magnitude and changing trend of the target correction value. Calculate the change rate of the target correction value over a certain time period. For example, calculate the ratio of the difference between the target correction values of two adjacent detections to the previous target correction value to obtain the change rate of the target correction value. At the same time, combine the distribution of the safety factor in various parts of the curtain wall to analyze the correlation between the target correction value and the safety factor. For example, observe whether there is an obvious increasing trend in the target correction value in areas with a low safety factor. Based on the analysis results, formulate warning rules. For example, if the target correction value in a certain area is greater than the preset high-risk threshold and the safety factor is in the medium to low-risk range, or the change rate of the target correction value exceeds a certain percentage (such as 20%), then it is determined that a warning mark needs to be added to this area. The warning rules can be flexibly adjusted according to the actual situation to ensure that potential risks can be discovered promptly and accurately. In the digital twin model, add warning marks to the corresponding areas according to the warning rules. The warning marks can take different forms, such as flashing icons, prominent text prompts, etc. For example, add a red exclamation mark icon in the high-risk area and display a brief risk prompt message beside it, such as "There is a high safety risk here, please check it in time." At the same time, set different priorities for the warning marks so that users can quickly focus on the most urgent risk areas.

[0090] Step 504: Collect the basic information of the curtain wall, including the type of the curtain wall (such as unitized curtain wall, framed curtain wall), area, service life, design unit, etc.; the information describing the detection process, such as detection time, detection equipment (such as laser point cloud scanner, model of wind pressure sensor, etc.), detection methods (such as finite element analysis method, wavelet transform analysis method, etc.); the risk level of the curtain wall determined in Step 501; the distribution of the corrected safety factor obtained in Step 502; the warning marks added in Step 503 and other data. Classify and organize these data to ensure the accuracy and integrity of the data. Design the template of the 3D detection report and determine the structure and content layout of the report. The report can include parts such as cover, table of contents, overview, detailed detection results, conclusions and suggestions.

[0091] Fill the sorted data into the corresponding positions according to the requirements of the report template. For the risk level, display the area proportion of different risk level regions in the form of a chart; for the distribution of the safety factor, insert a 3D visualization image and provide an interactive operation function to facilitate users to view from different angles; for the warning marks, list the detailed warning area locations, risk descriptions and recommended measures and other information. Finally, review and proofread the report. After ensuring that the content is accurate, generate the final 3D detection report. The report can be output in PDF format or 3D interactive file format (such as HTML5 format) for easy viewing and sharing by users.

[0092] Suppose there is a high-rise commercial building with a unitized glass curtain wall structure on its exterior facade. In a certain regular inspection: The initial safety factor is obtained through preliminary calculation = 0.7, the target correction value = 0.15, and the set fusion weight = 0.8. According to the formula + × , the corrected safety factor S of the curtain wall structure is calculated as S = 0.7 + 0.8×0.15 = 0.82. Comparing with the preset risk level classification standard, the low-risk range is 0.8 - 1.0. Since 0.8 ≤ 0.82 ≤ 1.0, it is determined that the risk level of this curtain wall is low risk. Divide the digital twin model evenly into 1000 grid units, and record and mark the corrected safety factor corresponding to each grid one by one. According to the preset rules, mark the grids with the safety factor in the range of 0.8 - 0.9 as "light green", and the grids in the range of 0.9 - 1.0 as "dark green". Process the grid coefficients through a data smoothing algorithm to eliminate the abnormal values caused by detection errors or data fluctuations, and finally form a detailed list of the safety factor distribution of each part of the curtain wall. It can be clearly seen from the list that the safety factors of most regions are in the green marked interval, indicating that the overall safety condition is good.

[0093] In the comparative analysis of the detection data, it was found that in the grid area numbered A-123 in the southeast corner of the curtain wall, the target correction value increased from 0.1 in the previous detection to 0.2 in this detection. After calculation, the change rate was = 100%, significantly exceeding the preset threshold of 20%. At the same time, the safety factor of this area is 0.65, which is in the medium-risk range. According to the established warning rules, a special identification symbol "△" is added to the corresponding grid position in the digital twin model and set to a red flashing state. At the same time, a prompt message "The target correction value in area A-123 in the southeast corner has changed too much, there is a potential risk" is added next to the grid entry in the detection data record form.

[0094] Collect the basic information of the curtain wall, including the type of the curtain wall is unit glass curtain wall, the total area is 5000 square meters, and the service life is 5 years, etc.; detailed record the information of the detection process, such as the detection time is October 1, 2024, and the detection equipment used includes a certain type of laser point cloud scanner, wind pressure sensor, etc., and the detection methods adopted cover technical means such as finite element analysis and wavelet transform. Arrange the risk level assessment results, the distribution list of the safety factors of each grid, the warning mark positions and risk descriptions, etc. according to the standard report template. In the finally generated three-dimensional detection report, the risk level distribution ratio is presented in the form of a data table, the safety factor of each area is described in the form of text combined with a data list, and the grid numbers, specific positions and risk contents of all warning marks are listed in detail.

[0095] By superimposing the target correction value on the initial safety factor, the impacts of various factors such as local deformation and dynamic displacement on the curtain wall safety are comprehensively considered, enabling a more accurate assessment of the actual safety condition of the curtain wall, avoiding the limitations of single-factor assessment, and improving the accuracy and reliability of the assessment results. The dynamic adjustment of the initial safety factor is realized, which can timely reflect the changes in the safety performance of the curtain wall caused by factors such as environmental changes and load effects during its use. This dynamic adjustment mechanism makes the safety assessment more in line with the actual situation and helps to timely detect potential safety hazards. The generated three-dimensional inspection report visually displays the distribution of the corrected safety factor through the digital twin model and adds warning marks. Managers and technicians can quickly and comprehensively understand the safety conditions of various parts of the curtain wall through the three-dimensional visualization image, locate high-risk areas, and improve the efficiency of problem detection and decision-making. Warning marks are added according to the correlation between the target correction value and the safety factor distribution, enabling early warning of potential risks. Users can take corresponding maintenance and repair measures in a timely manner based on the warning information, reducing the probability of safety accidents, ensuring the safety of personnel's lives and property, and also helping to extend the service life of the curtain wall and reduce maintenance costs. The three-dimensional inspection report integrates various aspects of information such as the basic information of the curtain wall, the description of the inspection process, the risk level, the safety factor distribution, and warning marks, providing comprehensive and detailed data for the maintenance, management, and supervision of the curtain wall. It is convenient for different personnel to understand the safety situation of the curtain wall from different perspectives, promoting information sharing and communication.

[0096] As Figure 2 shown, an embodiment of the present invention further provides a safety inspection device for the deformation condition of a curtain wall based on laser point cloud, including: A data processing module, configured to perform coordinate system normalization processing on the three-dimensional point cloud data set covering the entire surface of the curtain wall, generate a digital twin model of the curtain wall surface, and define a first detection target point and a second detection target point in the digital twin model according to the geometric characteristics and material joint characteristics of the curtain wall panels; A model construction module, configured to extract the surface normal vector distribution characteristics according to the digital twin model, fuse the measured wind pressure load distribution data and temperature field parameters, construct a dynamic flexural deformation coupling calculation model, and output global deformation parameters and local deformation parameters; A feature extraction module, configured to extract the vibration displacement components of the first detection target point and the second detection target point through phase correlation analysis and wavelet transform according to the local deformation parameters to construct a curtain wall vibration mode feature matrix. At the same time, analyze the displacement-time series data to identify out-of-limit displacement events with a displacement amplitude greater than or equal to a preset threshold; A safety degree evaluation module, configured to determine the initial safety factor of the curtain wall structure according to the global deformation parameters and the occurrence frequency of out-of-limit displacement events, and calculate the target correction value based on the local deformation parameters and dynamic displacement components; A report generation module, which is used to superimpose the target correction value and the initial safety factor to generate a corrected curtain wall structure safety factor, realize the dynamic adjustment of the initial safety factor, and output a three-dimensional detection report, including the risk level, the distribution of the corrected safety factor, and a warning mark associated with the target correction value.

[0097] It should be noted that this device corresponds to the above method, and all implementation manners in the above method embodiments are applicable to this embodiment and can achieve the same technical effects.

[0098] An embodiment of the present invention also provides a computing device, including: a processor and a memory storing a computer program. When the computer program is run by the processor, it executes the method as described above. All implementation manners in the above method embodiments are applicable to this embodiment and can achieve the same technical effects.

[0099] An embodiment of the present invention also provides a computer-readable storage medium storing instructions. When the instructions are run on a computer, the computer is made to execute the method as described above. All implementation manners in the above method embodiments are applicable to this embodiment and can achieve the same technical effects.

[0100] The above is the preferred embodiment of the present invention. It should be pointed out that for those of ordinary skill in the art, without departing from the principle of the present invention, several improvements and refinements can be made, and these improvements and refinements should also be regarded as the protection scope of the present invention.

Claims

1. A curtain wall deformation safety detection method based on laser point cloud, characterized in that: The method comprises: Step 1: normalize the coordinate system of the three-dimensional point cloud data set covering the entire surface of the curtain wall to generate a digital twin model of the curtain wall surface, and define the first detection target point and the second detection target point in the digital twin model according to the geometric features and material bonding characteristics of the curtain wall panels; Step 2: According to the digital twin model, the surface normal vector distribution characteristics are extracted, and the measured wind pressure load distribution data and temperature field parameters are integrated to construct a dynamic flexural deformation coupling calculation model, and output global deformation parameters and local deformation parameters; Step 3: According to the local deformation parameters, the vibration displacement components of the first detection target point and the second detection target point are extracted through phase correlation analysis and wavelet transform to construct the curtain wall vibration modal characteristic matrix. At the same time, the displacement-time series data is analyzed to identify the over-limit displacement events whose displacement amplitude is greater than or equal to the preset threshold; Step 4: Determine the initial safety coefficient of the curtain wall structure based on the global deformation parameters and the frequency of occurrence of over-limit displacement events, and calculate the target correction value based on the local deformation parameters and dynamic displacement components; Step 5, superimpose the target correction value and the initial safety coefficient to generate a corrected curtain wall structure safety coefficient, realize dynamic adjustment of the initial safety coefficient, and output a three-dimensional detection report, including risk level, corrected safety coefficient distribution and warning marks associated with the target correction value.

2. The curtain wall deformation safety detection method based on laser point cloud according to claim 1 is characterized in that: According to the digital twin model, the surface normal vector distribution characteristics are extracted, and the measured wind pressure load distribution data and temperature field parameters are integrated to build a dynamic flexural deformation coupling calculation model, and output global deformation parameters and local deformation parameters, including: According to the surface point cloud data of the digital twin model, the curtain wall surface is gridded and divided into blocks, and for the point cloud data in each block area, each target point in the current block and the neighboring point set within the preset radius are extracted; The spatial plane of the neighboring point set is fitted by the least square method, the surface normal vector of each target point is calculated, and the normal vector of each block is aggregated to generate the normal vector distribution characteristics of the curtain wall surface; By mapping the wind pressure sensor data around the building curtain wall to the corresponding spatial coordinates of the digital twin model, the measured wind pressure load distribution data is generated; based on the data collected by the temperature sensor and the heat transfer parameters of the building structure, the temperature field parameters, including the temperature spatial distribution matrix and the time-varying gradient sequence, are generated through finite element heat transfer analysis; The normal vector distribution characteristics of the curtain wall surface, the measured wind pressure load distribution data and the temperature field parameters are integrated to construct a dynamic flexural deformation coupling calculation model; Input wind pressure load distribution data, temperature field parameters and material elastic modulus into the dynamic flexural deformation coupling calculation model, and set fixed constraint conditions for the connection nodes between the curtain wall and the main structure; Dynamic flexural deformation analysis is performed based on fixed constraint conditions, wind pressure load distribution data, temperature field parameters and material elastic modulus to extract global and local deformation parameters.

3. The curtain wall deformation safety detection method based on laser point cloud according to claim 2 is characterized in that: The global deformation parameters include the overall waviness and the regional curvature anomaly; the local deformation parameters include the corner stress gradient of the first detection target point and the dynamic flexural displacement peak value at the edge-center junction of the second detection target point.

4. The curtain wall deformation safety detection method based on laser point cloud according to claim 3 is characterized in that: According to the local deformation parameters, the vibration displacement components of the first detection target point and the second detection target point are extracted through phase correlation analysis and wavelet transform to construct the curtain wall vibration modal characteristic matrix. At the same time, the displacement-time series data is analyzed to identify the over-limit displacement events with a displacement amplitude greater than or equal to the preset threshold, including: Determine, according to the local deformation parameters, the dynamic deformation sensitive areas corresponding to the first detection target point and the second detection target point in the time series point cloud data; For the point cloud sequence of the corner area of ​​the first detection target point in the sensitive area, the displacement changes of adjacent time frames are extracted through phase correlation analysis, and the low-frequency thermal expansion displacement component and the high-frequency wind vibration displacement component are separated by wavelet transform; the high-frequency wind vibration displacement component is used as the corner vibration displacement component; for the point cloud sequence of the junction area of ​​the second detection target point, the dynamic displacement component within the preset frequency band is extracted by wavelet transform; The vibration displacement component of the corner of the first detection target point is aligned with the dynamic displacement component of the junction of the second detection target point in time and space to construct the curtain wall vibration modal characteristic matrix; The displacement-time series data in the vibration modal characteristic matrix is ​​analyzed, and a displacement amplitude threshold is preset to identify single over-limit events with a displacement amplitude greater than or equal to the preset displacement amplitude threshold.

5. The curtain wall deformation safety detection method based on laser point cloud according to claim 4 is characterized in that: According to the global deformation parameters and the frequency of over-limit displacement events, the initial safety coefficient of the curtain wall structure is determined, and based on the local deformation parameters and dynamic displacement components, the target correction value is calculated, including: The historical accumulation item of wind pressure is obtained based on the wind pressure from the starting time to the current time and the accumulation of wind pressure attenuation over time; the temperature field data of the curtain wall surface obtained by laser point cloud scanning is analyzed, and the temperature gradient reflecting the local thermal stress condition of the curtain wall is obtained in combination with the spatial coordinate information of the curtain wall; the fatigue cumulative damage is determined based on the number of cycles and corresponding fatigue life of the curtain wall at different stress levels; The historical accumulation items of wind pressure, the gradient of temperature field at the spatial position and the fatigue accumulation damage are integrated to determine the static related items; Determine the maximum displacement of the curtain wall junction under dynamic load, that is, the peak value of dynamic flexural displacement, and the maximum vibration displacement of the curtain wall corner in high-frequency wind vibration after frequency band filtering, that is, the peak value of corner vibration displacement after filtering, and obtain the comprehensive parameters of the degree of mutual influence between deformation and vibration of the curtain wall structure under dynamic load; According to the statics-related items and the comprehensive parameters of the mutual influence between the deformation and vibration of the curtain wall structure under dynamic loads, the relationship between the statics risk and the deformation-vibration sensitivity of the curtain wall structure is determined; according to the peak value of the corner vibration displacement after filtering and the peak value of the dynamic displacement of the curtain wall junction after filtering in the same frequency band, the relationship between the vibration displacement synergy between the corner and the junction is obtained, that is, the dynamics-related items; According to the static related items and the dynamic related items, the contributions of the static factors and the dynamic factors to the target correction value are determined respectively to obtain the final target correction value.

6. The curtain wall deformation safety detection method based on laser point cloud according to claim 5 is characterized in that: The target correction value is superimposed on the initial safety coefficient to generate a corrected curtain wall structure safety coefficient, realize dynamic adjustment of the initial safety coefficient, and output a three-dimensional inspection report containing risk level, distribution of corrected safety coefficients and warning marks associated with the target correction value, including: The dynamic adjustment information reflected by the target correction value, including local deformation and dynamic displacement components, is integrated into the initial safety factor to generate a corrected curtain wall structure safety factor; According to the revised curtain wall structure safety factor and the preset risk level classification standard, the risk level of the curtain wall is determined; The corrected safety factor is mapped to the digital twin model of the curtain wall surface, and the distribution of the corrected safety factor in various parts of the curtain wall is processed to obtain the distribution of the corrected safety factor in various parts of the curtain wall; According to the size and change trend of the target correction value, as well as its correlation with the distribution of the safety factor in various parts of the curtain wall, corresponding warning marks are added to the digital twin model; The risk level of the curtain wall, the distribution of the safety factor in various parts of the curtain wall, the early warning marks, the basic information of the curtain wall and the detection process description, including the detection time, detection equipment and detection method, are integrated to generate a three-dimensional detection report.

7. The curtain wall deformation safety detection method based on laser point cloud according to claim 6 is characterized in that: The first detection target point is located at the four corner areas of the curtain wall panel to monitor the stress concentration effect at the corners; the second detection target point is located at the junction of the center and the edge of the curtain wall panel to monitor the temperature difference deformation gradient change.

8. A curtain wall deformation safety detection device based on laser point cloud, the device implements the method according to any one of claims 1 to 7, characterized in that: include: A data processing module is used to normalize the coordinate system of the three-dimensional point cloud data set covering the entire surface of the curtain wall, generate a digital twin model of the curtain wall surface, and define a first detection target point and a second detection target point in the digital twin model according to the geometric features and material bonding characteristics of the curtain wall panels; The model building module is used to extract the surface normal vector distribution characteristics based on the digital twin model, integrate the measured wind pressure load distribution data and temperature field parameters, build a dynamic flexural deformation coupling calculation model, and output global deformation parameters and local deformation parameters; A feature extraction module is used to extract the vibration displacement components of the first detection target point and the second detection target point according to the local deformation parameters through phase correlation analysis and wavelet transform to construct the curtain wall vibration modal feature matrix, and at the same time, analyze the displacement-time series data to identify the over-limit displacement events whose displacement amplitude is greater than or equal to the preset threshold; Safety assessment module, used to determine the initial safety coefficient of the curtain wall structure based on the global deformation parameters and the frequency of occurrence of over-limit displacement events, and calculate the target correction value based on the local deformation parameters and dynamic displacement components; The report generation module is used to superimpose the target correction value with the initial safety coefficient to generate a corrected curtain wall structure safety coefficient, realize dynamic adjustment of the initial safety coefficient, and output a three-dimensional detection report, including risk level, corrected safety coefficient distribution and warning marks associated with the target correction value.

9. A computing device, characterized in that include: one or more processors; A storage device for storing one or more programs, when the one or more programs are executed by the one or more processors, the one or more processors implement the method as claimed in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a program, which, when executed by a processor, implements the method according to any one of claims 1 to 7.

Images