Method and device for measuring residual deformation rate of cylinder volume

Abstract

The application provides a kind of cylinder volume residual deformation rate measurement method and device, comprising: obtaining the first point cloud data before the pressure test of the measured cylinder, the second point cloud data in pressure maintaining stage and the third point cloud data after pressure relief;According to the first point cloud data, the second point cloud data and the third point cloud data, respectively through the ring surface point set of layered extraction along point cloud main shaft direction and the cross-sectional area of generation initial layered ring sequence are calculated;According to the initial layered ring sequence, the interval of adjacent layered ring axial spacing greater than the preset threshold is filled by interpolation to generate the final layered ring sequence;According to the final layered ring sequence, the product of each ring cross-sectional area and corresponding axial spacing is added to obtain the volume before pressure test, the volume in pressure maintaining stage and the volume after pressure relief in corresponding state;According to the volume before pressure test, the volume in pressure maintaining stage and the volume after pressure relief, the residual deformation rate of cylinder volume is calculated.

Description

Technical Field

[0001] This application relates to the field of pressure vessel testing technology, and more particularly to a method and apparatus for measuring the residual deformation rate of gas cylinder volume. Background Technology

[0002] Gas cylinders, as important mobile pressure vessels, are widely used in industrial and civil fields. During repeated filling and discharging, they are prone to fatigue and plastic deformation, which may eventually lead to failure. Therefore, accurately measuring the residual volumetric deformation rate of gas cylinders during periodic inspections is crucial to ensuring their safe operation.

[0003] Currently, the most widely adopted methods for measuring the residual deformation rate of gas cylinder volume in relevant domestic and international standards are the external and internal methods in hydrostatic testing. The external method requires placing the gas cylinder under test in a sealed water jacket filled with water. By measuring the amount of water discharged from the jacket during pressurization and the amount of water returned after depressurization, the total deformation value and elastic deformation value are obtained, and then the residual deformation rate is calculated. The internal method calculates the volume deformation by directly measuring the change in the amount of water injected into the gas cylinder. However, these traditional methods have significant limitations: they require fixed facilities such as large water jackets, hydrostatic pump stations, or precise water volume measuring devices, resulting in large equipment investments and large land areas; during testing, the gas cylinder must be disassembled from the usage site and transported to a dedicated testing station, making in-situ testing impossible. For fixed or large gas cylinder groups such as those using long-tube trailers, disassembly and transportation costs are high and the cycle is long; furthermore, the measurement results are only a single value characterizing the overall volume change of the gas cylinder, and cannot reflect the local plastic deformation distribution of different parts of the cylinder (such as the cylinder body and end caps), making it difficult to accurately locate weak points and assess failure risks.

[0004] In recent years, some studies have attempted to use 3D vision or laser scanning technology to acquire surface data of gas cylinders to measure residual deformation rates. Other studies have explored methods such as strain measurement and optical imaging to measure residual deformation of gas cylinders. Among these, optical imaging methods involve placing multiple cameras around the gas cylinder to capture images of the cylinder surface at different pressure stages. Feature point extraction and matching analysis of feature point displacement are then used to calculate diameter or volume changes. However, these methods rely on image resolution and feature point matching accuracy. For the minute elastic and plastic deformations generated during gas cylinder hydrostatic testing, measurement accuracy is difficult to guarantee. Furthermore, they are easily affected by environmental factors such as ambient light, gas cylinder surface reflection, and oil contamination, limiting their applicability and reliability.

[0005] Therefore, there is an urgent need to provide a gas cylinder residual deformation rate detection scheme that does not rely on geometric model reconstruction or regular shape fitting and can achieve high-precision and high-robust volume measurement. Summary of the Invention

[0006] The purpose of this application is to provide a method and device for measuring the residual deformation rate of gas cylinder volume, so as to solve the problems of computational complexity, insufficient robustness or poor accuracy caused by relying on geometric model reconstruction or regular shape fitting in the prior art, and to achieve high-precision and high-robust measurement of gas cylinder volume.

[0007] To achieve the above objectives, the method for measuring the residual deformation rate of gas cylinder volume provided in this application specifically includes: acquiring first point cloud data before the pressure resistance test of the gas cylinder to be tested, second point cloud data during the pressure holding stage, and third point cloud data after depressurization; generating an initial layered ring sequence by extracting toroidal point sets layer by layer along the main axis of the point cloud and calculating the cross-sectional area based on the first, second, and third point cloud data; generating a final layered ring sequence by interpolating and filling the intervals where the axial spacing between adjacent layered rings is greater than a preset threshold based on the initial layered ring sequence; obtaining the volume before the pressure resistance test, the volume during the pressure holding stage, and the volume after depressurization in the corresponding state by multiplying and summing the cross-sectional areas of each ring with the corresponding axial spacing based on the final layered ring sequence; and calculating the residual deformation rate of the gas cylinder volume based on the volume before the pressure resistance test, the volume during the pressure holding stage, and the volume after depressurization.

[0008] Optionally, the residual deformation rate of the gas cylinder volume can be calculated based on the volume before the pressure test, the volume during the pressure holding stage, and the volume after depressurization. This includes: calculating the residual deformation rate of the gas cylinder volume by dividing the difference between the volume after depressurization and the volume before the pressure test by the difference between the volume during the pressure holding stage and the volume before the pressure test.

[0009] Optionally, in the above method, the process of extracting toroidal point sets and calculating cross-sectional areas along the principal axis of the point cloud to generate an initial layered toroidal sequence includes: calculating the covariance matrix of the point cloud data using principal component analysis, determining the direction of the eigenvector corresponding to the largest eigenvalue as the principal axis of the point cloud; dividing the point cloud data into multiple initial intervals along this direction, and selecting points whose normal vectors satisfy preset conditions to form toroidal point sets; projecting the toroidal point sets onto a plane perpendicular to the principal axis of the point cloud and performing triangulation, and obtaining the layer cross-sectional area by accumulating the areas of each triangular element obtained from the triangulation to form the initial layered toroidal sequence.

[0010] Optionally, in the above method, the triangulation is a Delaunay triangulation.

[0011] Optionally, in the above method, generating the final layered ring sequence by interpolating and filling intervals where the axial spacing between adjacent layered rings is greater than a preset threshold based on the initial layered ring sequence includes: calculating the axial spacing between adjacent layered rings along the principal axis of the point cloud; when the axial spacing is greater than a preset spacing threshold, inserting a filling ring in the corresponding interval; wherein, the cross-sectional area of ​​the filling ring is determined by linear interpolation based on the cross-sectional areas of two adjacent layered rings, and the axial coordinate of the filling ring is set within the corresponding interval by linear interpolation.

[0012] Optionally, in the above method, obtaining the volume before the pressure test, the volume during the pressure holding stage, and the volume after pressure release in the corresponding state by multiplying and summing the cross-sectional areas of each ring with their corresponding axial spacing according to the final layered ring sequence includes: for each layered ring in the final layered ring sequence, multiplying the cross-sectional area of ​​the layered ring by the axial spacing between the layered ring and the next adjacent layered ring along the principal axis of the point cloud as the micro-element volume, and summing the micro-element volumes corresponding to all layered rings in the final layered ring sequence to obtain the corresponding volume.

[0013] Optionally, the above method may also include: based on the first point cloud data and the third point cloud data, comparing the changes in the cross-sectional area of ​​each layer along the main axis of the point cloud, generating and outputting a deformation distribution map to characterize the axial local plastic deformation of the gas cylinder.

[0014] Optionally, in the above method, the first point cloud data, the second point cloud data, and the third point cloud data are all acquired by a 3D laser scanner.

[0015] This application also provides a device for measuring the residual deformation rate of a gas cylinder volume. The device includes a point cloud acquisition module, an initial layered ring sequence generation module, a final layered ring sequence generation module, a volume calculation module, and a residual deformation rate calculation module. The point cloud acquisition module is used to acquire first point cloud data before the pressure resistance test of the gas cylinder, second point cloud data during the pressure holding stage, and third point cloud data after pressure release. The initial layered ring sequence generation module is used to generate corresponding state-specific point cloud data by extracting ring surface point sets layer by layer along the main axis of the point cloud and calculating the cross-sectional area based on the first point cloud data, the second point cloud data, and the third point cloud data. The initial layered ring sequence is generated; the final layered ring sequence generation module is used to generate the final layered ring sequence in the corresponding state by interpolating and filling the intervals where the axial spacing of adjacent layered rings is greater than a preset threshold, based on the initial layered ring sequence; the volume calculation module is used to obtain the volume before the pressure test, the volume during the pressure holding stage, and the volume after depressurization in the corresponding state by multiplying and summing the cross-sectional areas of each ring with the corresponding axial spacing, based on the final layered ring sequence; the residual deformation rate calculation module is used to calculate the residual deformation rate of the gas cylinder volume based on the volume before the pressure test, the volume during the pressure holding stage, and the volume after depressurization.

[0016] This application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the above-described method.

[0017] This application also provides a computer-readable storage medium storing a computer program that performs the above-described methods.

[0018] This application also provides a computer program product, including a computer program / instructions that, when executed by a processor, implement the steps of the above-described method.

[0019] The beneficial technical effects of this application are as follows: It directly calculates the gas cylinder volume from point cloud data by using a method of extracting toroidal point sets layer by layer along the main axis of the point cloud, adaptive interpolation filling, and integral accumulation. This eliminates the need to pre-reconstruct the point cloud data into a closed 3D mesh model or fit it to a regular geometric body. This fundamentally avoids the problems of volume calculation failure or error amplification caused by local point cloud omissions, noise, or topological ambiguity during mesh reconstruction. Even if there are missing data in local areas, as long as the remaining point set of that layer is sufficient to fit the cross-sectional profile, the cross-sectional area of ​​that layer can be accurately calculated, significantly improving the robustness of volume calculation and the tolerance for point cloud data quality. Simultaneously, the aforementioned layered processing naturally produces cross-sectional area data at various positions along the gas cylinder's axial direction. By comparing the changes in the cross-sectional area of ​​each layer before and after the pressure test, a deformation distribution map of the cylinder along the axis can be generated. This allows for intuitive and quantitative location of the area where the cylinder undergoes the maximum plastic deformation, providing crucial diagnostic information that traditional methods cannot obtain for the refined assessment of the cylinder's safety status and prediction of its remaining life. Furthermore, by reducing the three-dimensional volume calculation to the integral accumulation of two-dimensional cross-sectional areas, the cross-sectional area calculation of each layer can be processed independently and in parallel, significantly reducing computational complexity and memory usage. Moreover, the entire process is based entirely on the two-dimensional geometric relationship of the local point cloud, without relying on the selection of external reference points. This avoids the floating-point numerical errors introduced by improper reference points in the traditional tetrahedral directional volume method, resulting in superior numerical stability. Attached Figure Description

[0020] The accompanying drawings, which are included to provide a further understanding of this application and form part of this application, do not constitute a limitation thereof. In the drawings: Figure 1 This is a schematic flowchart of a method for measuring the residual deformation rate of gas cylinder volume provided in an embodiment of this application; Figure 2 This is a schematic diagram illustrating the process of generating an initial hierarchical ring sequence according to an embodiment of this application; Figure 3 This is a schematic diagram illustrating the process of generating the final hierarchical ring sequence according to an embodiment of this application; Figure 4A schematic diagram of the gas cylinder volume residual deformation rate measuring device provided in an embodiment of this application; Figure 5 This is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. Detailed Implementation

[0021] Various exemplary embodiments, features, and aspects of this application will now be described in detail with reference to the accompanying drawings. The same reference numerals in the drawings denote elements that have the same or similar functions. Although various aspects of the embodiments are shown in the drawings, they are not necessarily drawn to scale unless specifically indicated otherwise.

[0022] The term “exemplary” as used herein means “serving as an example, embodiment, or illustration.” Any embodiment illustrated herein as “exemplary” is not necessarily to be construed as superior to or better than other embodiments.

[0023] In this document, the term "and / or" is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent three cases: A alone, A and B simultaneously, and B alone. Furthermore, the term "at least one" in this document means any combination of at least two of any one or more elements. For example, including at least one of A, B, and C can mean including any one or more elements selected from the set consisting of A, B, and C.

[0024] Furthermore, to better illustrate this application, numerous specific details are provided in the following detailed description. Those skilled in the art should understand that this application can be implemented without certain specific details. In some instances, methods, means, components, and circuits well-known to those skilled in the art have not been described in detail in order to highlight the main points of this application.

[0025] To facilitate a clearer explanation of the specific implementation details of this application, the following provides a unified interpretation of the key terms involved in the specific implementation process of this application: Layered ring: Along the main axis of the point cloud, at a specific axial position, a ring-shaped profile is formed by a set of points extracted from the point cloud data that are approximately located in the same cross section. Each layered ring corresponds to the cross section of the gas cylinder at that axial position.

[0026] Toroidal point set: A subset of point clouds that meets preset screening conditions (the angle between the normal vector and the principal axis of the point cloud is within a preset range) and is used to form a certain layered toroidal ring.

[0027] Initial layered ring sequence: An ordered set of layered rings obtained after extracting point cloud data in layers along the main axis of the point cloud according to a preset initial spacing.

[0028] Final layered ring sequence: After interpolating and filling the intervals where the distance between adjacent layered rings in the initial layered ring sequence is greater than a preset threshold, the resulting ordered set of a series of layered rings in which the distance between adjacent layered rings does not exceed the preset threshold is used for volume integral calculation.

[0029] Infinite element volume: For a small axial segment defined by two adjacent layered rings in the final layered ring sequence, the infinitesimal volume obtained by multiplying the cross-sectional area of ​​one of the layered rings by the axial spacing of the segment is used to calculate the total volume by integration and summation.

[0030] Principal axis direction of point cloud: After calculating the point cloud data using principal component analysis, the direction is determined by the eigenvector corresponding to the largest eigenvalue. This direction represents the main extension direction of the spatial distribution of the point cloud data.

[0031] Based on the above explanation, please refer to Figure 1 As shown in one embodiment of this application, the method for measuring the residual deformation rate of gas cylinder volume provided in this application specifically includes: S101 acquires the first point cloud data before the pressure resistance test of the gas cylinder under test, the second point cloud data during the pressure holding stage, and the third point cloud data after the pressure is released. In practical applications, the three-dimensional geometric information of the outer surface of the gas cylinder under test is obtained at three different stages of the pressure test, providing raw data for subsequent volume calculations. The acquisition of point cloud data is fundamental to the method described in this application, and its quality and completeness directly affect the accuracy of subsequent volume calculations.

[0032] In this application, the first point cloud data, the second point cloud data, and the third point cloud data can all be acquired using a 3D laser scanner. A 3D laser scanner is a precision measuring device that uses a laser beam to scan the surface of an object and obtains the three-dimensional coordinates of a large number of points on the object's surface by measuring laser reflection time or using triangulation principles. Alternatively, structured light scanners or multi-view photogrammetry can also be used to acquire point cloud data; these technologies can also achieve the purpose of obtaining the object's three-dimensional geometric information.

[0033] S102 Based on the first point cloud data, the second point cloud data, and the third point cloud data, an initial layered ring sequence is generated by extracting toroidal point sets layer by layer along the main axis of the point cloud and calculating the cross-sectional area. In practice, this step transforms the problem of calculating the volume of a three-dimensional gas cylinder into a problem of accumulating two-dimensional cross-sectional integrals along the main axis of the cylinder. The advantages of this "dimensionality reduction" are: first, it avoids the complex process of reconstructing a three-dimensional surface; second, each layer can be calculated independently, facilitating parallel processing and improving computational efficiency; and third, the calculation process naturally generates cross-sectional area information at each axial position, providing a data foundation for subsequent local deformation analysis.

[0034] The principal axis direction of the point cloud is the main extension direction of the point cloud spatial distribution, determined through statistical analysis of the point cloud data. Ideally, this direction coincides with the geometric axis of the gas cylinder. By layering along this direction, each cross-section can be made as perpendicular as possible to the gas cylinder surface, thus making the calculation of the cross-sectional area more accurate.

[0035] Extracting the toroidal point set requires satisfying certain constraints on the angle between the normal vectors. The purpose is to filter out the points that are truly located on the circumferential contour of the outer surface of the gas cylinder and exclude the deviation points caused by scanning angle deviation or surface defects.

[0036] The triangulation is a Delaunay triangulation. Delaunay triangulation has good empty circle properties and maximizes the minimum angle property, enabling the generation of high-quality triangular meshes, which is beneficial for improving the accuracy of area calculation.

[0037] S103 Based on the initial layered ring sequence, the final layered ring sequence is generated by interpolating and filling the intervals where the axial spacing between adjacent layered rings is greater than a preset threshold. Because the initial layering is performed at a fixed interval, and the point cloud density on the surface of the gas cylinder may be uneven, or the toroidal point sets in some areas may be of poor quality and discarded, there are significant differences in the axial spacing between adjacent layered rings. If the integration is performed directly using a sequence of layered rings with uneven spacing, it will result in a large error in volume calculation.

[0038] In practical applications, the adaptive filling mechanism is designed to address this issue. By inserting filling rings into intervals with spacing greater than a preset threshold, the spacing between adjacent rings in the final layered ring sequence is ensured to remain within the preset threshold, thus guaranteeing the accuracy of the volume integral. The cross-sectional area and axial coordinates of the filling rings are determined through linear interpolation, based on the reasonable assumption that the change in cylinder profile between adjacent cross-sections is typically gradual.

[0039] S104 obtains the volume before the pressure test, the volume during the pressure holding stage, and the volume after pressure release in the corresponding state by summing the product of the cross-sectional area of ​​each ring and the corresponding axial spacing according to the final layered ring sequence. In practice, this step mainly applies the concept of integration to volume calculation. The gas cylinder is considered along its axial direction as being composed of countless tiny cylindrical elements stacked together. The volume of each element is approximately the cross-sectional area of ​​one end face multiplied by its height (i.e., the axial spacing). Summing the volumes of all the elements yields an approximate value for the total volume of the gas cylinder. When the spacing between the layers is sufficiently small, this approximation can approach the true volume with arbitrarily high precision.

[0040] S105 calculates the residual deformation rate of the gas cylinder volume based on the volume before the pressure test, the volume during the pressure holding stage, and the volume after the pressure is released.

[0041] Residual deformation rate is a core indicator for evaluating the safety of gas cylinders. According to relevant standards, gas cylinders with a residual deformation rate exceeding a certain limit should be deemed unqualified and scrapped or downgraded for use.

[0042] In one embodiment of this application, the residual deformation rate of the gas cylinder volume is calculated based on the volume before the pressure test, the volume during the pressure holding stage, and the volume after depressurization. This includes: calculating the residual deformation rate of the gas cylinder volume by dividing the difference between the volume after depressurization and the volume before the pressure test by the difference between the volume during the pressure holding stage and the volume before the pressure test.

[0043] Specifically, in practical work, the volume before the pressure test is set. The above calculations were performed on the point cloud model after the experiment to obtain the post-experiment volume. .

[0044] The difference in volume before and after the test is the residual deformation of the gas cylinder. : ; Obtain the total deformation value measured during the pressure test. The total deformation value can be obtained by scanning the gas cylinder and calculating its volume during the pressure holding stage, or by estimating the elastic deformation value using traditional methods (such as external measurement).

[0045] residual deformation Divide by the total deformation value Multiply by 100% to get the residual deformation rate. : ×100%; The physical meaning of residual deformation rate is: what proportion of the total deformation of a gas cylinder under test pressure is irreversible permanent deformation. The higher this ratio, the more severe the plastic damage to the gas cylinder material, and the higher the safety risk of continued use.

[0046] Please refer to Figure 2 As shown, in one embodiment of this application, generating an initial layered ring sequence by extracting toroidal point sets layer by layer along the principal axis of the point cloud and calculating the cross-sectional area includes: S201 calculates the covariance matrix of point cloud data using principal component analysis and determines the direction of the eigenvector corresponding to the largest eigenvalue as the direction of the principal axis of the point cloud. Principal Component Analysis (PCA) is a classic method for data dimensionality reduction and feature extraction. Its core idea is to extract the principal feature components of the data through linear transformation. In this application, the logic of PCA is as follows: the gas cylinder has a long and slender shape, and its point cloud data exhibits the greatest dispersion in the axial direction and relatively less dispersion in the radial direction. PCA, through eigenvalue decomposition of the covariance matrix, can find the direction of maximum data dispersion, which is the principal axis direction of the point cloud. Therefore, the principal axis direction is entirely determined by the statistical characteristics of the point cloud data itself, without the need to pre-assume a standard axis for the gas cylinder, thus enabling adaptive adaptation to any placement of the gas cylinder in the scanning coordinate system.

[0047] S202 divides the point cloud data into multiple initial intervals along this direction, and selects points whose normal vectors and the angle between them and the principal axis of the point cloud satisfy a preset condition to form a toroidal point set; The point cloud data is divided into multiple initial intervals along the principal axis of the point cloud at a preset initial spacing (e.g., 0.1 mm). Each interval corresponds to a thin slice at an axial position. For each point within an initial interval, the angle between its normal vector and the principal axis of the point cloud is calculated. The wind vector can be obtained by performing local PCA calculation on the set of points in the neighborhood of that point. The purpose of setting angle constraints (e.g., within 0.5°) is to filter out points in the point cloud whose normal vectors are approximately perpendicular to the principal axis of the point cloud, i.e., points that are truly located on the circumferential contour of the outer surface of the gas cylinder. Points with large deviations in normal vector direction may come from scanning noise, surface impurities, or attachments on the gas cylinder; excluding them can improve the accuracy of cross-sectional area calculation.

[0048] S203 projects the toroidal point set onto a plane perpendicular to the principal axis of the point cloud and performs triangulation. The layer cross-sectional area is obtained by accumulating the areas of each triangular element obtained by the triangulation to form the initial layered toroidal sequence.

[0049] The selected toroidal point set is projected onto a plane perpendicular to the principal axis of the point cloud, resulting in a two-dimensional planar point set. This projection operation transforms the toroidal contour in three-dimensional space into a closed curve on a two-dimensional plane, facilitating area calculation. Delaunay triangulation is then performed on the two-dimensional point set, dividing the interior region of the closed curve into several triangles. The cross-sectional area of ​​this layer is obtained by summing the areas of each triangular element. The layered rings and their cross-sectional areas obtained along all intervals along the principal axis of the point cloud are arranged in axial order to form the initial layered ring sequence.

[0050] Please refer to Figure 3 As shown, in one embodiment of this application, generating a final layered ring sequence by interpolating and filling intervals where the axial spacing between adjacent layered rings is greater than a preset threshold based on the initial layered ring sequence includes: S301 calculates the axial spacing between adjacent layered rings along the principal axis of the point cloud; Traverse each pair of adjacent layered rings in the initial layered ring sequence and calculate their axial spacing along the principal axis of the point cloud.

[0051] S302 When the axial spacing is greater than the preset spacing threshold, insert a filler ring in the corresponding interval; The cross-sectional area of ​​the filling ring is determined by linear interpolation based on the cross-sectional areas of two adjacent layered rings, and the axial coordinate of the filling ring is set within the corresponding interval by linear interpolation.

[0052] In practical work, a preset spacing threshold is used. These parameters are set according to the accuracy requirements of volume calculation. In this invention, It can be set to 0.025mm to match the single-point accuracy of the 3D laser scanner used. The logic behind this setting is that the accuracy of volume integration should not be limited by the uniformity of the layer spacing, but should remain on the same order of magnitude as the accuracy of the original data.

[0053] When d ≤ In this case, the interval does not need to be filled, and the original layered ring is directly retained.

[0054] When d > When k fill rings are inserted in this interval, where: ; The cross-sectional area of ​​the filling ring is determined by linear interpolation based on the cross-sectional areas of two adjacent layered rings. The axial coordinates of the filling ring are set within corresponding intervals using linear interpolation. The axial coordinates of the filling ring are evenly distributed between adjacent layered rings, ensuring that the distance between each sub-interval does not exceed [a certain value]. The cross-sectional area of ​​the filling ring is determined by linear interpolation, based on the reasonable assumption that the cross-sectional area of ​​the gas cylinder varies linearly along the axial direction within a sufficiently small axial spacing. The filling ring is merged with the original layered rings and sorted according to their axial coordinates to form the final layered ring sequence. After adaptive filling, the axial spacing between adjacent layered rings is no greater than [value missing]. This provides a uniform, high-density layered basis for subsequent volume integral calculations.

[0055] In one embodiment of this application, the volume before the pressure test, the volume during the pressure holding stage, and the volume after the pressure release are obtained by multiplying and summing the cross-sectional areas of each ring with their corresponding axial spacing according to the final layered ring sequence. This includes: for each layered ring in the final layered ring sequence, multiplying the cross-sectional area of ​​the layered ring by the axial spacing between the layered ring and the next adjacent layered ring along the principal axis of the point cloud as the micro-element volume, and summing the micro-element volumes corresponding to all layered rings in the final layered ring sequence to obtain the corresponding volume.

[0056] The physical meaning of this calculation method is: the gas cylinder is regarded along the axial direction as being composed of countless tiny cylindrical micro-elements stacked together, and the volume of each micro-element is approximately the cross-sectional area of ​​one of its end faces multiplied by the height of the micro-element; since the spacing between adjacent layered rings does not exceed the preset spacing threshold, the error introduced by this approximation is extremely small.

[0057] In one embodiment of this application, the method further includes: comparing the changes in the cross-sectional area of ​​each layer along the main axis of the point cloud based on the first point cloud data and the third point cloud data, and generating and outputting a deformation distribution map to characterize the axial local plastic deformation of the gas cylinder.

[0058] The specific implementation method is as follows: For the first point cloud data before the pressure test and the third point cloud data after depressurization, the final layered ring sequence is generated according to the aforementioned method, obtaining the cross-sectional area data of the layered rings at each axial position along the principal axis of the point cloud. The same principal axis of the point cloud is used for both processing to ensure the comparability of the cross-sectional area data in the axial position. For each axial position, the change in cross-sectional area after depressurization relative to the cross-sectional area before the test is calculated, and then the radial deformation is deduced.

[0059] Color mapping technology is used to map radial deformation onto a 3D model of the gas cylinder surface. A color scale is set: red indicates large deformation (greater than 0.5 mm), yellow indicates medium deformation (0.2 mm to 0.5 mm), and green indicates small deformation (less than 0.2 mm). Through the deformation distribution map, inspectors can visually locate the area of ​​maximum plastic deformation in the cylinder, providing crucial diagnostic information that traditional methods cannot obtain for refined assessment of the cylinder's safety status and prediction of its remaining life. The ability to generate the deformation distribution map in this application is naturally obtained during the volume calculation process. This is because the layered ring integral method already generates cross-sectional area data for each axial position during the calculation; comparing the cross-sectional areas at different stages only requires simple subtraction, which is not available in the traditional tetrahedral directional volume method.

[0060] Please refer to Figure 4 As shown, this application also provides a gas cylinder volume residual deformation rate measuring device, the device including a point cloud acquisition module, an initial layered ring sequence generation module, a final layered ring sequence generation module, a volume calculation module and a residual deformation rate calculation module; The point cloud acquisition module is used to acquire the first point cloud data before the pressure resistance test of the gas cylinder under test, the second point cloud data during the pressure holding stage, and the third point cloud data after the pressure is released; the module can communicate with a 3D laser scanner to receive the point cloud data of each stage collected by the scanner.

[0061] The initial layered ring sequence generation module is used to generate an initial layered ring sequence in the corresponding state by extracting toroidal point sets and calculating the cross-sectional area along the principal axis direction of the point cloud based on the first point cloud data, the second point cloud data, and the third point cloud data. The module determines the principal axis direction of the point cloud by principal component analysis, divides the interval along the direction, selects toroidal point sets that meet the constraint conditions of the angle between the normal vectors, projects and triangulates them, and calculates the cross-sectional area of ​​each layer to form the initial layered ring sequence.

[0062] The final layered ring sequence generation module is used to generate the final layered ring sequence in the corresponding state by interpolating and filling the intervals where the axial spacing between adjacent layered rings is greater than a preset threshold, based on the initial layered ring sequence; the cross-sectional area and axial coordinates of the filled rings are determined by linear interpolation.

[0063] The volume calculation module is used to obtain the volume before the pressure test, the volume during the pressure holding stage, and the volume after the pressure release in the corresponding state by summing the product of the cross-sectional area of ​​each ring and the corresponding axial spacing according to the final layered ring sequence. The residual deformation rate calculation module is used to calculate the residual deformation rate of the gas cylinder volume based on the volume before the pressure test, the volume during the pressure holding stage, and the volume after the pressure release.

[0064] The specific implementation logic of each of the above modules has been described in detail in the foregoing embodiments, and will not be described in detail here.

[0065] To facilitate a clearer understanding of the gas cylinder volume residual deformation rate measurement method and apparatus provided in this application, a large-volume seamless steel gas cylinder with a nominal volume of 3755L is used as the test object to illustrate the complete implementation process and technical effects of the method and apparatus of this application in detail; those skilled in the art will understand that this example does not limit this application in any way.

[0066] I. Test Subjects and Test Equipment

[0067] A batch of large-capacity seamless steel gas cylinders with a nominal volume of 3755L was selected as the test object. 1. 3D laser scanner: The selected device is the TSCAN hawk2, which has a single-point accuracy of 0.025mm, and can meet the requirements for high-precision surface data acquisition.

[0068] 2. Hydraulic pressure test system: As a comparative verification and to provide a reference for the full deformation value, the traditional external measurement method hydraulic pressure test device is used, equipped with a high-precision water meter and pressure sensor.

[0069] II. Experimental Procedure

[0070] Step 1: Point Cloud Data Acquisition

[0071] 1. Marker Placement: Evenly attach positioning markers to the outer surface of the gas cylinder. The markers are highly reflective circular targets with a diameter of 6 mm. The spacing between the markers is controlled within 200 mm, and they are evenly distributed along the axial and circumferential directions of the cylinder to ensure that multiple scan data can be stitched together with high precision through the markers. The markers are evenly distributed in different parts of the cylinder.

[0072] 2. Pre-test scanning: Using a TSCAN hawk2 high-precision 3D laser scanner, the gas cylinder was scanned from different angles to obtain a 3D point cloud dataset of the outer surface of the gas cylinder before the test. The scanning parameters are set as follows: Scan resolution: 0.025 mm The scanning path covers all parts of the cylinder, including the body, end caps, and neck. The overlap of each scan is no less than 30% to ensure the accuracy of subsequent stitching. After scanning, a 3D point cloud dataset of the cylinder's outer surface before the experiment is obtained. ).

[0073] 3. Pressure Resistance Test: The pressure test is conducted according to GB / T 9251 "Method for Hydrostatic Testing of Gas Cylinders": The gas cylinder is placed in a water jacket, the pressurization pipeline and measuring instrument are connected, and water is injected into the water jacket to the specified water level to remove air from the system; Pressurization stage: The pressure is slowly increased to the test pressure of 33.4 MPa at a rate not exceeding 0.5 MPa / s; Pressure holding stage: The pressure is held at 33.4 MPa for 120 s, and the gas cylinder is scanned again to obtain a three-dimensional point cloud dataset of the outer surface of the gas cylinder after the test. The measured total deformation value was subsequently calculated and recorded. (); Depressurization stage: Slowly depressurize to normal pressure and record the residual deformation.

[0074] 4. Post-test scanning: After the water pressure test is depressurized, the gas cylinder is removed from the water jacket. After the cylinder surface is dry, the gas cylinder is scanned again using the same scanning parameters and marker point positions as in step 1.2 to obtain a 3D point cloud dataset of the outer surface of the gas cylinder after the test. By using pasted positioning markers, it is ensured that the scan data before and after the experiment have a unified coordinate system reference.

[0075] Step 2: Point cloud data preprocessing

[0076] The collected raw point cloud data ( ) and( Perform the following preprocessing operations respectively: 1. Noise Reduction: A statistical filtering method is used to calculate the average distance between each point and other points in its neighborhood, eliminating outlier noise points that do not conform to a Gaussian distribution. Specific parameter settings are as follows: number of neighborhood points: 50; standard deviation factor: 2.0; effectively eliminating isolated noise points and flying points generated during the scanning process.

[0077] 2. Point Cloud Stitching: Based on the positioning markers, the Iterative Closest Point (ICP) algorithm is used to accurately stitch together multiple point clouds to form a complete point cloud model of the gas cylinder's outer surface. The specific steps are as follows: Initial registration is performed using the marker coordinates to obtain a coarse transformation matrix; the ICP algorithm is used to iteratively optimize the transformation parameters to minimize the distance error between point clouds.

[0078] 3. Volume Calculation: The cylinder volume is calculated using an adaptive hierarchical loop integral method, which includes the following sub-steps: Reference axis determination: The principal axis direction of the gas cylinder point cloud data was calculated using Principal Component Analysis (PCA). Calculate the centroid coordinates of the point cloud data: ; Construct the covariance matrix: ; Find the eigenvalues ​​and eigenvectors of the covariance matrix; The direction of the eigenvector corresponding to the largest eigenvalue is taken as the reference axis. This axis minimizes the sum of the squares of the perpendicular distances from all points to the axis, ensuring that the slice direction matches the bottle's geometry.

[0079] Layered loop extraction: Based on the point cloud data and the normal vector constraint of the reference axis, layered loops are extracted along the reference axis direction. The point cloud data is sorted along the reference axis to obtain the axial coordinate range. In this embodiment ( ), ( ).

[0080] Set the initial layer spacing ( The gas cylinder is divided into several initial intervals along its axial direction, for a total of: ; For each initial interval, point clouds that satisfy the condition that the angle between the normal vector and the reference axis is within (0.5°) are selected as the torus point set for that interval. The normal vector is calculated through principal component analysis of the point cloud neighborhood.

[0081] For each torus point set, the geometric center of the layered torus is obtained. The torus point set is projected onto a plane perpendicular to the reference axis to obtain a two-dimensional point set. This two-dimensional point set is then decoupled using Delaunay triangulation to form a mesh composed of two-dimensional triangular elements. All triangular elements within the layer are traversed, and the area of ​​each triangular element is calculated using Heron's formula. The areas of all triangular elements are then summed to obtain the cross-sectional area of ​​the layer. ): ; in( ) represents the total number of triangular fractals within this layer. ) is the ( The area of ​​each triangular piece.

[0082] Taking a typical sheet as an example, the calculated cross-sectional area ( The corresponding diameter is approximately 1200 mm.

[0083] If the number of toroidal point sets in a certain interval is insufficient or the fitting error exceeds the threshold (the threshold is set to the number of toroidal point sets being less than 10 or the relative error of the cross-sectional area being greater than 5%), then linear interpolation of the toroidal point sets in adjacent intervals is used to supplement the layered toroidal ring.

[0084] Adaptive padding: The spacing between adjacent layer rings is adaptively padded to ensure that the distance between adjacent rings along the reference axis does not exceed a set accuracy threshold. ): Calculate adjacent layered rings ( ) and( Distance along the reference axis: ; like( If the original layered ring is not filled, then no filling is needed; the original layered ring can be directly preserved.

[0085] like( If the two layered rings are inserted, a filling ring is inserted between them, as follows: 1) Calculate the number of rings that need to be inserted: ; 2) Determine the axial position of each filling ring: ; 3) Determine the area of ​​each filling ring using a linear interpolation method: ; The filled rings are merged with the original layered rings to form a new sequence of layered rings, ensuring that the spacing between adjacent rings is no greater than ( ). ).

[0086] Calculations show that the initial number of layered rings was 9000. After adaptive filling, the total number of layered rings expanded to: ; The spacing between adjacent rings is controlled within 0.025 mm.

[0087] Volume integral calculation: Calculate the volume of each infinitesimal element from left to right, and sum all the infinitesimal element volumes to obtain the total volume of the gas cylinder. For two adjacent layered rings (including the original ring and the filling ring), let their axial distance be ( The area of ​​the left-side layered ring is ( ). If ), then the volume of the infinitesimal element is: ; By summing the volumes of all the infinitesimal elements, we obtain the total volume of the gas cylinder: ; in( ) represents the total number of layered rings after filling. ) is the ( ) layered ring and the ( The axial distance between the layered rings.

[0088] Perform the above calculations on the point cloud model before the experiment to obtain the volume before the experiment ( Perform the above calculations on the point cloud model after the experiment to obtain the volume after the experiment ( ).

[0089] 4. Calculation and evaluation of residual deformation

[0090] Calculation of residual deformation: The difference in volume before and after the test is the residual deformation of the gas cylinder. ): ; Calculation of residual deformation rate: Obtaining the total deformation value during the pressure test ( In this embodiment, the total deformation value is measured during the pressure holding stage using the traditional external measurement method. The residual deformation amount ( Divide by the total deformation value ( Multiply by 100% to get the residual deformation rate. ): ; Local deformation visualization analysis: By comparing and analyzing the two 3D models before and after the experiment, deformation distribution maps of different parts of the gas cylinder are generated. The distance difference between corresponding points in the point cloud model before and after the experiment was calculated. For each point cloud location, its radial deformation was calculated. Color mapping technology was used to map the radial deformation onto the surface of the gas cylinder. A color scale was set: red indicates large deformation (>0.5 mm), yellow indicates medium deformation (0.2~0.5 mm), and green indicates small deformation (<0.2 mm). The area where the cylinder experienced the maximum plastic deformation was located intuitively.

[0091] The beneficial technical effects of this application are as follows: It directly calculates the gas cylinder volume from point cloud data by using a method of extracting toroidal point sets layer by layer along the main axis of the point cloud, adaptive interpolation filling, and integral accumulation. This eliminates the need to pre-reconstruct the point cloud data into a closed 3D mesh model or fit it to a regular geometric body. This fundamentally avoids the problems of volume calculation failure or error amplification caused by local point cloud omissions, noise, or topological ambiguity during mesh reconstruction. Even if there are missing data in local areas, as long as the remaining point set of that layer is sufficient to fit the cross-sectional profile, the cross-sectional area of ​​that layer can be accurately calculated, significantly improving the robustness of volume calculation and the tolerance for point cloud data quality. Simultaneously, the aforementioned layered processing naturally produces cross-sectional area data at various positions along the gas cylinder's axial direction. By comparing the changes in the cross-sectional area of ​​each layer before and after the pressure test, a deformation distribution map of the cylinder along the axis can be generated. This allows for intuitive and quantitative location of the area where the cylinder undergoes the maximum plastic deformation, providing crucial diagnostic information that traditional methods cannot obtain for the refined assessment of the cylinder's safety status and prediction of its remaining life. Furthermore, by reducing the three-dimensional volume calculation to the integral accumulation of two-dimensional cross-sectional areas, the cross-sectional area calculation of each layer can be processed independently and in parallel, significantly reducing computational complexity and memory usage. Moreover, the entire process is based entirely on the two-dimensional geometric relationship of the local point cloud, without relying on the selection of external reference points. This avoids the floating-point numerical errors introduced by improper reference points in the traditional tetrahedral directional volume method, resulting in superior numerical stability.

[0092] This application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the above-described method.

[0093] This application also provides a computer-readable storage medium storing a computer program that performs the above-described methods.

[0094] This application also provides a computer program product, including a computer program / instructions that, when executed by a processor, implement the steps of the above-described method.

[0095] like Figure 5 As shown, the electronic device 600 may also include: a communication module 110, an input unit 120, an audio processor 130, a display 160, and a power supply 170. It is worth noting that the electronic device 600 does not necessarily need to include these components. Figure 5 All components shown; in addition, the electronic device 600 may also include Figure 5 For components not shown, please refer to existing technologies.

[0096] like Figure 5 As shown, the central processing unit 100, sometimes also referred to as a controller or operating control, may include a microprocessor or other processor device and / or logic device. The central processing unit 100 receives inputs and controls the operation of various components of the electronic device 600.

[0097] The memory 140 may be, for example, one or more of a cache, flash memory, hard drive, removable media, volatile memory, non-volatile memory, or other suitable devices. It may store the aforementioned failure-related information, and also store a program for executing that information. The central processing unit 100 may execute the program stored in the memory 140 to perform information storage or processing, etc.

[0098] Input unit 120 provides input to central processing unit 100. Input unit 120 may be, for example, a keypad or touch input device. Power supply 170 provides power to electronic device 600. Display 160 displays images and text. Display may be, for example, an LCD display, but is not limited thereto.

[0099] The memory 140 can be a solid-state memory, such as a read-only memory (ROM), random access memory (RAM), a SIM card, etc. It can also be a memory that retains information even when power is off, can be selectively erased, and contains more data; examples of this type of memory are sometimes referred to as EPROMs. The memory 140 can also be some other type of device. The memory 140 includes a buffer memory 141 (sometimes referred to as a buffer). The memory 140 may include an application / function storage unit 142 for storing application programs and function programs or processes for executing the operation of the electronic device 600 via the central processing unit 100.

[0100] The memory 140 may also include a data storage unit (data 143) for storing data, such as contacts, digital data, pictures, sounds, and / or any other data used by the electronic device. The driver storage unit (driver 144) of the memory 140 may include various drivers for the electronic device's communication functions and / or for performing other functions of the electronic device (such as messaging applications, address book applications, etc.).

[0101] The communication module 110 is a transmitter / receiver 110 that transmits and receives signals via antenna 111. The communication module (transmitter / receiver) 110 is coupled to the central processing unit 100 to provide input signals and receive output signals, which can be the same as in a conventional mobile communication terminal.

[0102] Based on different communication technologies, multiple communication modules 110 can be configured in the same electronic device, such as cellular network modules, Bluetooth modules, and / or wireless LAN modules. The communication module (transmitter / receiver) 110 is also coupled to a speaker 131 and a microphone 132 via an audio processor 130 to provide audio output via the speaker 131 and receive audio input from the microphone 132, thereby enabling typical telecommunications functions. The audio processor 130 may include any suitable buffer, decoder, amplifier, etc. Additionally, the audio processor 130 is coupled to a central processing unit 100, enabling on-device recording via the microphone 132 and on-device playback of stored audio via the speaker 131.

[0103] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0104] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0105] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1The function specified in one or more boxes.

[0106] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0107] The description of the various embodiments above tends to emphasize the differences between the various embodiments. The similarities or similarities between them can be referred to, and for the sake of brevity, they will not be repeated here.

[0108] Those skilled in the art will understand that, in the above-described method of the specific implementation, the order in which each step is written does not imply a strict execution order and does not constitute any limitation on the implementation process. The specific execution order of each step should be determined by its function and possible internal logic.

[0109] The various embodiments of this disclosure have been described above. These descriptions are exemplary and not exhaustive, and are not limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments. The terminology used herein is chosen to best explain the principles, practical application, or improvement of the technology in the market, or to enable others skilled in the art to understand the embodiments disclosed herein.

Claims

1. A method for measuring the residual deformation rate of a gas cylinder volume, characterized in that, include: Acquire the first point cloud data before the pressure resistance test of the gas cylinder under test, the second point cloud data during the pressure holding stage, and the third point cloud data after the pressure is released; Based on the first point cloud data, the second point cloud data, and the third point cloud data, an initial layered ring sequence is generated by extracting toroidal point sets layer by layer along the main axis of the point cloud and calculating the cross-sectional area. The final layered ring sequence is generated by interpolating and filling the intervals where the axial spacing between adjacent layered rings is greater than a preset threshold, based on the initial layered ring sequence. Based on the final layered ring sequence, the volumes before the pressure test, the pressure holding stage, and the pressure release stage under the corresponding conditions are obtained by multiplying and summing the cross-sectional areas of each ring with the corresponding axial spacing. The residual deformation rate of the gas cylinder volume is calculated based on the volume before the pressure test, the volume during the pressure holding stage, and the volume after the pressure is released.

2. The method according to claim 1, characterized in that, The residual deformation rate of the gas cylinder volume, calculated based on the volume before the pressure test, the volume during the pressure holding stage, and the volume after depressurization, includes: Based on the volume before the pressure test, the volume during the pressure holding stage, and the volume after depressurization, the residual deformation rate of the gas cylinder volume is calculated by dividing the difference between the volume after depressurization and the volume before the pressure test by the difference between the volume during the pressure holding stage and the volume before the pressure test.

3. The method according to claim 1, characterized in that, The initial layered ring sequence is generated by extracting toroidal point sets layer by layer along the principal axis of the point cloud and calculating the cross-sectional area, including: The covariance matrix of the point cloud data is calculated by principal component analysis, and the direction of the eigenvector corresponding to the largest eigenvalue is determined as the direction of the principal axis of the point cloud. The point cloud data is divided into multiple initial intervals along this direction, and points whose normal vectors make an angle with the principal axis of the point cloud that satisfies a preset condition are selected to form a toroidal point set. The toroidal point set is projected onto a plane perpendicular to the principal axis of the point cloud and triangulated. The cross-sectional area of ​​each triangular element obtained by the triangulation is accumulated to form the initial layered toroidal sequence.

4. The method according to claim 3, characterized in that, The triangulation is the Delaunay triangulation.

5. The method according to claim 1, characterized in that, The final layered ring sequence is generated by interpolating and filling intervals where the axial spacing between adjacent layered rings is greater than a preset threshold, based on the initial layered ring sequence. Calculate the axial spacing between adjacent layered rings along the principal axis of the point cloud; When the axial spacing is greater than the preset spacing threshold, a filler ring is inserted in the corresponding interval; The cross-sectional area of ​​the filling ring is determined by linear interpolation based on the cross-sectional areas of two adjacent layered rings, and the axial coordinate of the filling ring is set within the corresponding interval by linear interpolation.

6. The method according to claim 1 or 5, characterized in that, Based on the final layered ring sequence, by summing the product of the cross-sectional area of ​​each ring and the corresponding axial spacing, the volumes before the pressure test, the pressure holding stage, and the pressure release stage under the corresponding conditions are obtained, including: For each layer ring in the final layer ring sequence, the cross-sectional area of ​​the layer ring is multiplied by the axial distance between the layer ring and the next adjacent layer ring along the principal axis of the point cloud as the micro-element volume. The corresponding micro-element volumes of all layer rings in the final layer ring sequence are accumulated and summed to obtain the corresponding volume.

7. The method according to claim 1, characterized in that, Also includes: Based on the first and third point cloud data, the changes in the cross-sectional area of ​​each layer along the main axis of the point cloud are compared to generate and output a deformation distribution map to characterize the axial local plastic deformation of the gas cylinder.

8. The method according to claim 1, characterized in that, The first point cloud data, the second point cloud data, and the third point cloud data were all acquired using a 3D laser scanner.

9. A device for measuring the residual deformation rate of gas cylinder volume, characterized in that, The device includes a point cloud acquisition module, an initial layered ring sequence generation module, a final layered ring sequence generation module, a volume calculation module, and a residual deformation rate calculation module. The point cloud acquisition module is used to acquire the first point cloud data before the pressure resistance test of the gas cylinder under test, the second point cloud data during the pressure holding stage, and the third point cloud data after the pressure is released. The initial layered ring sequence generation module is used to generate an initial layered ring sequence in the corresponding state by extracting ring surface point sets layer by layer along the main axis direction of the point cloud and calculating the cross-sectional area based on the first point cloud data, the second point cloud data and the third point cloud data. The final layered ring sequence generation module is used to generate the final layered ring sequence in the corresponding state by interpolating and filling the intervals where the axial spacing between adjacent layered rings is greater than a preset threshold, based on the initial layered ring sequence. The volume calculation module is used to obtain the volume before the pressure test, the volume during the pressure holding stage, and the volume after the pressure release in the corresponding state by summing the product of the cross-sectional area of ​​each ring and the corresponding axial spacing according to the final layered ring sequence. The residual deformation rate calculation module is used to calculate the residual deformation rate of the gas cylinder volume based on the volume before the pressure test, the volume during the pressure holding stage, and the volume after the pressure release.

10. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the method of any one of claims 1 to 8.

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

12. A computer program product comprising a computer program / instructions, characterized in that, When the computer program / instructions are executed by the processor, they implement the steps of the method according to any one of claims 1 to 8.

Images