Bridge crane finite element calculation resource optimization method and system based on dynamic similarity discrimination

Abstract

The application discloses a bridge hoisting equipment finite element calculation resource optimization method based on dynamic similarity discrimination, which comprises the following steps: S10, setting a working condition and a finite element model, receiving a secondary input working condition, updating the finite element model, performing finite element calculation, storing information, outputting the calculation result, and outputting the secondary input working condition; S20, judging the source of the hoisting equipment of the input working condition and calling a library table, outputting if the input working condition exists, and otherwise, proceeding to S30; S30, if the model geometric parameters of the input working condition exist in the database table, proceeding to S40, otherwise, proceeding to S10; S40, finding the working condition closest to the input working condition according to the load information of the input working condition, updating the load information and proceeding to S10 if the closest working condition cannot be found, and otherwise, proceeding to S50; S50, if the input working condition and the working condition load information in S10 meet the conditions, performing interpolation calculation, outputting the calculation result, and storing the information to the corresponding database table; otherwise, updating the load difference information and proceeding to S10; and the corresponding system is also disclosed.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of bridge construction, more particularly, to a bridge hoisting equipment finite element calculation resource optimization method and system based on dynamic similarity discrimination. BACKGROUND

[0002] In the bridge construction site, it is crucial to ensure the safe operation of the hoisting equipment and effectively monitor its structural mechanical response. The traditional monitoring method relies on real-time measurement by sensors and real-time online calculation by finite element programs deployed on cloud servers, which can meet the monitoring needs to a certain extent. However, as the scale of the construction site expands and the number of hoisting equipment increases, the drawbacks of the traditional method become more apparent. The finite element online calculation program has high requirements for server performance, and when faced with multiple devices and high concurrency tasks, it is prone to insufficient computing resources such as high CPU occupancy and insufficient memory, which can cause real-time calculation monitoring to be interrupted for some devices, posing risks to construction safety. In addition, each calculation data of each hoisting equipment is stored in a database table, and the data volume grows explosively over time, which reduces the speed of writing cloud computing results and querying and reading edge data, seriously affecting the real-time performance of safety monitoring and failing to meet the efficient and stable monitoring needs.

[0003] To solve these problems, some existing technology patents disclose a large piece transportation bridge safety evaluation system based on distributed cloud computing, which includes multiple units such as evaluation and calculation, and related database tables. This system solves the problems of long evaluation and calculation time, non-standard reports prone to errors, and difficulty in modifying standards, but it mainly focuses on the architecture and process of safety evaluation and cannot reduce the demand for computing resources and storage capacity. Some other existing technology proposes a bridge response rapid calculation method, which involves determining bridge information, establishing a model, conducting experiments, and closed-loop calculation, and has the advantages of good real-time calculation and high precision. However, it focuses on bridge response calculation and vehicle-bridge interaction simulation and does not design a complete online calculation and storage logic, which cannot solve the problem of explosive growth of data volume, slow writing of cloud computing results, and slow querying and reading of edge data. SUMMARY

[0004] In view of the above defects or improvement needs of the prior art, the present application provides a bridge hoisting equipment finite element calculation resource optimization method based on dynamic similarity discrimination, which first stores multiple working condition finite element calculations, then judges whether the input working condition already exists in the database table through a similarity function, and outputs the result if it exists, avoiding repeated calculation. Based on the similarity of geometry and load information, the existing data and interpolation are used to avoid recalculation and reduce resource consumption, effectively solving the problem of insufficient computing resources in high concurrency.

[0005] To achieve the above purpose, according to the first aspect of the present application, a bridge hoisting equipment finite element calculation resource optimization method based on dynamic similarity discrimination is provided, comprising the following steps.

[0006] S10, database table setting and updating, setting multiple working conditions and finite element models according to practical data and receiving input working conditions and updating finite element models, performing finite element calculation on the working conditions, and storing the working conditions, corresponding finite element models and finite element calculation results into the database table corresponding to the hoisting equipment, if the input working conditions are secondary input, the calculation results also need to be output;

[0007] S20, equipment working condition searching, judging from which hoisting equipment the input working condition comes from, calling the corresponding database table, if the input working condition exists in the database table, directly outputting the result, otherwise entering S30;

[0008] S30, model geometric parameter checking and processing, if the model geometric parameter of the input working condition exists in the database table, entering S40, otherwise entering S10;

[0009] S40, load processing, finding the working condition in S10 closest to the input working condition according to the load information of the input working condition, if it cannot be found, entering S10 based on only updating the load information, otherwise entering S50;

[0010] S50, simplification processing, if the input working condition and the working condition in S10 meet the set conditions, performing interpolation calculation, outputting the interpolation calculation result and storing the input working condition and the result into the database table corresponding to the hoisting equipment; otherwise, entering S10 based on only updating the load difference information.

[0011] Further, the method for setting multiple working conditions and finite element models in S10 and receiving secondary input working conditions and updating finite element models, and performing finite element calculation on the working conditions is:

[0012] The total number of working conditions is N, i is the working condition number, The i-th group of working conditions includes a model geometric parameter set , a load information set and a boundary condition set Bi, the model geometric parameter set includes M model geometric parameters, the total number of model geometric parameters M is set according to experiments, m is the model geometric parameter number, The load information set includes load position and load size The boundary condition set Bi includes L boundary condition parameters, the total number of boundary condition parameters L is set according to experiments, l is the boundary condition parameter number, l=1, 2, …, L; for dynamic

[0013] ,

[0014] ,

[0015] ,

[0016] wherein, is the mass matrix of the i-th group of working conditions, is the moving speed of the equipment, is the moving acceleration of the equipment, is the damping matrix, is the stiffness matrix, is the displacement of the equipment, and t is time, is the external load vector varying with time t;

[0017] is the density of the equipment, and the function is a function for calculating the mass matrix, and the function g is a load calculation function of the working condition;

[0018] For the static said hoisting equipment:

[0019] ,

[0020] is the equivalent node force vector of the i-th group of working conditions.

[0021] Further, the method of storing the working conditions, corresponding finite element models and finite element calculation results in the database table of the corresponding hoisting equipment in S10 is:

[0022] The database table includes working condition, finite element model and finite element calculation result fields, and a database table row is added when storing new data; using working condition summary set , finite element calculation result summary set The working conditions and finite element calculation results are summarized respectively, and the two summary sets are stored in the cells of the corresponding fields of the new row of the database table in a set format;

[0023] ,

[0024] ,

[0025] In the formula, , , are the geometric parameter set, the load information set and the boundary condition set of the i-th group of working conditions respectively, , , , are the node displacement set, the element axial stress set and the element combined stress set of the i-th group of finite element calculation results respectively, and the element safety factor set is;

[0026] The finite element model set in S10 and the updated finite element model generate a unique code when calculating, add the unique code to the model name and store separately, and store the finite element model address and name to the cell of the corresponding field of the new row of the database table.

[0027] Further, the method for determining in S20 which hoisting equipment the input working condition comes from is:

[0028] The input working condition is represented by input set I:

[0029] ,

[0030] In the formula, , , respectively represent the geometric parameter set, the load information set, and the boundary condition set of the input working condition;

[0031] The hoisting equipment is F, and one set of working conditions is taken from each hoisting equipment database table, a total of F sets of working conditions, f is the working condition serial number corresponding to the hoisting equipment, and the input working condition is compared with the F sets of working conditions in similarity:

[0032] ,

[0033] ,

[0034] ,

[0035] In the formula, is the similarity measure function of the input working condition and the fth set of working conditions, is the geometric parameter weight coefficient, is the load information weight coefficient, is the boundary condition weight coefficient; is the geometric similarity function of the input working condition and the fth set of working conditions, is the load similarity function of the input working condition and the fth set of working conditions, is the boundary condition similarity function of the input working condition and the fth set of working conditions, and the similarity threshold value of the fth set of working conditions is set according to experiments If , the input working condition comes from the hoisting equipment corresponding to the fth set of working conditions.

[0036] Further, the method for checking and processing the model geometric parameters in S30 is:

[0037] ,

[0038] ,

[0039] wherein, is the mth model geometry parameter of the input working condition, is the mth model geometry parameter of the ith group of working conditions in S10, is the geometric similarity measure function; the geometric similarity threshold value is obtained according to experiments ;

[0040] If there is i=k, such that , then the model geometry parameter of the input working condition has a sufficient similar match in the database table, and S40 is entered; if all satisfy , then the model geometry parameter of the input working condition is not found in the database table, and otherwise the historical finite element model is called, and S10 is entered based on the updated model geometry parameter and load information.

[0041] Further, the method for finding the closest working condition in S10 according to the load information of the input working condition in S40 is as follows:

[0042] S41, the load information of the input working condition is input as , ), and are the load position and load size of the input working condition, and the load information of the working condition in S10 is , );

[0043] S42, a flag function foundsimilar is defined and initialized to false, and when = but ≠ or ≠ but = , the flag function foundsimilar becomes true;

[0044] S43, the load information of the input working condition , is compared with that of each working condition in S10 , , and when the flag function foundsimilar is true and are simultaneously satisfied, the corresponding working condition in S10 is collected, and S44 is entered, and otherwise the historical finite element model is called, and S10 is entered based on the updated load information.

[0045] S44, screening the working conditions collected in S43 to find two working conditions closest to the input working condition; for = But ≠ , the method for finding two working conditions closest to the input working condition is:

[0046] Find , so that ,

[0047] Find , so that ,

[0048] Working conditions a, b are the two working conditions closest to the input working condition, , The load size of working conditions a, b respectively;

[0049] For ≠ But = , the method for finding two working conditions closest to the input working condition is:

[0050] ,

[0051] ,

[0052] ,

[0053] ,

[0054] ,

[0055] ,

[0056] Working conditions c, d are the two working conditions closest to the input working condition; Pc, Pd are the load positions of working conditions c, d respectively, is the load position displacement of working condition d to working condition c, is the displacement of working condition d to the input working condition, wherein , , are the horizontal, vertical and vertical coordinates of the input working condition position respectively, , , are the horizontal, vertical and vertical coordinates of the working condition c position respectively, , , are the horizontal, vertical and vertical coordinates of the working condition d position respectively;

[0057] t is any value, if t is not present, call the history finite element model, only update the load information under the basis of S10; otherwise, continue to calculate the working conditions c and d through the index variables e and f.

[0058] Further, the condition set in S50 is:

[0059] S51, if the load position of the working condition output in S40 is the same as that of the input working condition but the load size is different, enter S52, otherwise enter the load size process;

[0060] S52, the working conditions output in S40 are respectively recorded as working condition a and working condition b, and the absolute value of the difference is calculated and , wherein is the load size of working condition a, is the load size of working condition b, , respectively, the absolute value of the difference between the load sizes of working conditions a and b and the input working condition;

[0061] S53, set the load size error tolerance according to the experiment , if and , then interpolate the data of working conditions a and b to obtain the calculation result of the input working condition.

[0062] Further, the load size process is:

[0063] S1, for the working condition output in S40 and the load size of the input working condition is the same but the load position is different, the working condition output in S40 is recorded as working condition c and working condition d respectively;

[0064] S2, calculate the Euclidean distance between the load positions of the two working conditions and the load position of the input working condition,

[0065] ,

[0066] ,

[0067] In the formula, , respectively, the distance between working conditions c and d and the input working condition;

[0068] S3, set the load position error tolerance according to the experiment , if and , then interpolate the data of working conditions c and d to obtain the calculation result of the input working condition.

[0069] ​According to a second aspect of the present application, a bridge hoisting equipment finite element calculation resource optimization system based on dynamic similarity discrimination is provided, comprising:

[0070] A database table setting and updating module sets various working conditions and finite element models according to practical data, receives input working conditions and updates the finite element models, performs finite element calculation on the working conditions, and stores the working conditions, corresponding finite element models and finite element calculation results into a database table corresponding to the hoisting equipment. If the input working conditions are secondary input, the calculation results also need to be output.

[0071] An equipment working condition searching module judges from which hoisting equipment the input working conditions come from, calls the corresponding database table, and directly outputs the results if the input working conditions exist in the database table, otherwise, enters S30.

[0072] A model geometric parameter checking and processing module enters S40 if the model geometric parameters of the input working conditions exist in the database table, otherwise, enters S10.

[0073] A load processing module finds the working condition in S10 closest to the input working conditions according to the load information of the input working conditions, enters S10 on the basis of only updating the load information if it cannot be found, otherwise, enters S50.

[0074] A simplification processing module performs interpolation calculation if the input working conditions and the working condition load information in S10 meet the set conditions, outputs the interpolation calculation results and stores the input working conditions and the results into the database table corresponding to the hoisting equipment, otherwise, enters S10 on the basis of only updating the load difference information.

[0075] Overall, compared with the prior art, the above technical solutions conceived by the present application can achieve the following beneficial effects:

[0076] 1. The calculation method of the present application avoids repeated calculation by first performing finite element calculation on multiple working conditions and storing them, and then judging whether the input working conditions exist in the database table through a similarity function. If they exist, the results are output. Based on the similarity of geometric and load information, the existing data and interpolation are used to avoid recalculation, reduce resource consumption, and effectively solve the problem of insufficient calculation resources during high concurrency.

[0077] 2. The calculation method of the present application stores working conditions in a special database table for different equipment, and sets similarity metrics and threshold values according to equipment characteristics, so that the system is suitable for various construction hoisting equipment. BRIEF DESCRIPTION OF DRAWINGS

[0078] Figure 1 The calculation method provided by the preferred embodiment of the present application is shown in the figure;

[0079] Figure 2 A bridge hoisting equipment finite element calculation resource optimization system structure schematic diagram based on dynamic similarity discrimination is provided for an embodiment of the present application. DETAILED DESCRIPTION

[0080] In order to make the purpose, technical solutions and advantages of the present application clearer and more apparent, the present application is further described in detail below in combination with the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and do not limit the present application. In addition, the technical features involved in each embodiment of the present application described below can be combined with each other as long as they do not conflict with each other.

[0081] Based on the above problems, the present application designs a bridge hoisting equipment finite element calculation resource optimization method and data storage logic based on dynamic similarity discrimination according to whether the model geometric parameters and the load application size and position of the finite element model change. Through the following steps, the demand for computing resources of the finite element online calculation program can be significantly reduced when high-concurrency tasks are processed in the cloud, and the data storage amount of the database table can be reduced, so as to achieve the purpose of guaranteeing the real-time and continuity of the safety monitoring of multiple hoisting equipment in the bridge construction site. By setting practical data as multiple working conditions and storing the finite element calculation results of these working conditions in the database table corresponding to the hoisting equipment, a rich calculation history information library is constructed. Using this database table, the efficiency of subsequent similar working condition calculation can be improved, and repeated calculation can be avoided. For different working conditions, different situations of dynamic and static hoisting equipment are also considered. For dynamic hoisting equipment, the calculation of its working condition involves the mutual relationship of multiple complex physical quantities such as mass matrix, equipment moving speed, acceleration, damping matrix, stiffness matrix, equipment displacement, time and external load vector; for static hoisting equipment, there is the calculation of equivalent node force vector, etc.

[0082] Through equipment working condition searching, the input working condition set is used to represent (including geometric parameters, load information, boundary condition set), combined with the similarity measurement function, weight coefficient and similarity function, compared with the database table working condition. According to the similarity threshold value, it is judged whether the input working condition exists or not to decide the subsequent operation, and a comprehensive multi-factor similarity judgment mechanism is constructed.

[0083] Model geometric parameter check: a special geometric similarity measurement function is used to compare the input and the existing geometric parameters in the database. According to the geometric similarity threshold value determined by the experiment, it is accurately judged whether the input geometric parameters can find similar matching in the library or not, and the fine consideration of the model geometric similarity is realized.

[0084] Load processing: distinguish the load position and size, set the corresponding processing flow for different similar cases. Use the flag function to mark the collection, calculate the absolute value and Euclidean distance of the difference, and combine the load size and position error tolerance set by the experiment to judge whether to interpolate calculation. This way can accurately find the closest working condition, take into account the subtle differences, ensure the accuracy of the calculation results, and ensure reasonable processing and storage of calculation results according to specific circumstances.

[0085] Please refer to Figure 1 The application relates to a bridge hoisting equipment finite element calculation resource optimization method based on dynamic similarity discrimination, comprising the following steps:

[0086] S10, database table establishment and update, a plurality of working conditions and finite element models are set according to practical data, and input working conditions and updated finite element models are received, the working conditions are subjected to finite element calculation, and the working conditions, corresponding finite element models and finite element calculation results are stored in the database table of the corresponding hoisting equipment, and the calculation results are also output if the input working conditions are secondary; that is, according to experience, common calculation working conditions 1, working condition 2, …, working condition N are set for each hoisting equipment, each calculation working condition contains finite element model geometric parameters and load size and position, and finite element calculation is performed in advance for each working condition, and the calculation working condition and result are stored in the database table list of the corresponding hoisting equipment.

[0087] In terms of data storage, an efficient database table management system is adopted, such as a relational database table (such as MySQL) or a non-relational database table (such as MongoDB), different tables or collections are used to store corresponding data according to different types of hoisting equipment. For each device, a special table or collection is created to ensure data isolation and manageability. Each table or collection assigns a unique identifier to each working condition to facilitate subsequent data retrieval and update operations. At the same time, the stored data not only includes the finite element calculation results, but also stores some key metadata, such as the creation time, last update time, and usage frequency of the working condition, so as to effectively manage and optimize the data.

[0088] The method for setting a plurality of working conditions and finite element models in S10 and receiving secondary input working conditions and updating finite element models, and performing finite element calculation on the working conditions is:

[0089] The total number of working conditions is N, i is the working condition number, The i-th group of working conditions includes a model geometric parameter set , a load information set and a boundary condition set Bi, the model geometric parameter set includes M model geometric parameters, the total number of model geometric parameters M is set according to the experiment, and m is the model geometric parameter number, ; the load information set includes load position and load size ; the boundary condition set Bi includes L boundary condition parameters, the total number of boundary condition parameters L is set according to experiments, l is the boundary condition parameter sequence number, l=1, 2, …, L; for dynamic hoisting equipment: (in some preferred embodiments, the boundary condition parameter is constant)

[0090] ,

[0091] ,

[0092] ,

[0093] wherein, is the mass matrix of the i-th group of working conditions, is the equipment moving speed, is the equipment moving acceleration, is the damping matrix, is the stiffness matrix, is the equipment displacement, t is time, is the external load vector varying with time t;

[0094] is the density of the equipment, the function is a function for calculating the mass matrix, and the function g is a load calculation function of the equipment;

[0095] For static hoisting equipment:

[0096] ,

[0097] is the equivalent node force vector of the i-th group of working conditions.

[0098] The method of storing the working conditions, corresponding finite element models, and finite element calculation results in the S10 to the database table of the corresponding hoisting equipment is:

[0099] The database table includes working condition, finite element model, and finite element calculation result fields, and a database table row is added when storing new data; using working condition summary set , finite element calculation result summary set The working conditions and finite element calculation results are summarized respectively, and the two summary sets are stored in the cells of the corresponding fields of the new row of the database table in a set format, and in some preferred embodiments, the json format is used;

[0100] ,

[0101] ,

[0102] wherein, , , are respectively the geometric parameter set, the load information set, and the boundary condition set of the i-th group of working conditions, , , , are respectively the node displacement set, the element axial stress set, the element combined stress set, and the element safety factor set of the i-th group of finite element calculation results;

[0103] The finite element model set in S10 and the updated finite element model generate a unique code when performing calculation, add the unique code to the model name and store separately, and store the finite element model address and name to the cell of the corresponding field of the newly added row of the database table. In some preferred embodiments, the method for separately storing the finite element model set in S10 and the updated finite element model is: in the local folder of the cloud server, the naming method is device name + UUID, UUID is a unique code generated simultaneously each time the working condition is calculated, and the cell under the "finite element model" field of the newly added row in the database table stores the finite element model address information and the finite element model name.

[0104] S20, device working condition searching, judging from which hoisting device the input working condition comes, calling the corresponding database table, if the input working condition exists in the database table, directly outputting the result, otherwise entering S30; that is, whenever the cloud receives a finite element calculation request, first judging which device it is, and then judging whether the input calculation working condition at this time exists in the database table, if it exists, directly outputting the calculation result under the corresponding working condition, saving the step of calling the finite element solver, if the input calculation working condition at this time does not exist in the database table, performing S30.

[0105] In the data retrieval process, first, according to the device type of the input working condition, locate the corresponding database table or set. Then, using the characteristic information of the input working condition, such as geometric parameters, load information and boundary conditions, construct a query statement. For the calculation of the similarity function, the required data will be retrieved from the database table, through multi-table connection or set query operation, to obtain the historical working conditions related to the input working condition. For example, for a relational database table, using SQL statements, according to the geometric parameters, load information and boundary conditions of the input working condition, using JOIN operation and WHERE clause to filter out potentially similar working conditions; for non-relational database tables, using document query or graph query method to find similar documents or nodes.

[0106] The method for judging from which hoisting device the input working condition comes in S20 is:

[0107] The input working condition is represented by an input set I:

[0108] ,

[0109] wherein, 、 、 respectively represent a set of geometric parameters, a set of load information, and a set of boundary conditions of the input working condition;

[0110] The F lifting devices each take one set of working conditions in each of the database tables of the lifting devices, and a total of F sets of working conditions are taken, f is the working condition serial number corresponding to the lifting device, and the input working condition is compared with the F sets of working conditions in terms of similarity:

[0111] ,

[0112] ,

[0113] ,

[0114] wherein, is a similarity measurement function of the input working condition and the fth set of working conditions, is a geometric parameter weight coefficient, is a load information weight coefficient, is a boundary condition weight coefficient; is a geometric similarity function of the input working condition and the fth set of working conditions, is a load similarity function of the input working condition and the fth set of working conditions, is a boundary condition similarity function of the input working condition and the fth set of working conditions, and a similarity threshold value is set for the fth set of working conditions according to experiments If , the input working condition is from the lifting device corresponding to the fth set of working conditions.

[0115] S30, model geometric parameter inspection and processing, if the model geometric parameter of the input working condition exists in the database table, S40 is entered, otherwise S10 is entered; that is, it is judged whether the model geometric parameter in the working condition at this time exists in the database table, if it exists, it is judged whether the load size or load position exists in the database table.

[0116] The method for model geometric parameter inspection and processing in S30 is:

[0117] ,

[0118] ,

[0119] wherein, is the mth model geometric parameter of the input working condition, for the mth model geometry parameter of the ith group of working conditions in S10, is a geometric similarity measure function; the geometric similarity threshold value is obtained according to experiments ;

[0120] If there is i=k, such that , the model geometry parameter of the input working condition has a similar match in the database table, then enter S40; if all satisfy , the model geometry parameter of the input working condition has no similar match in the database table, otherwise, the historical finite element model is called, and the model geometry parameter and the load information are updated to enter S10.

[0121] If the similarity threshold value is set to 1, even if the geometric change of the device is very small, the system will consider it a new working condition and need to be recalculated. However, these calculations may not be necessary because the geometric change is too small to have a significant impact on the final calculation results. Recalculation will consume a large amount of computing resources, including calculation time, performance of the computing device, etc. For example, in finite element analysis, each calculation may need to perform complex matrix operations and model solving, which will occupy a large amount of CPU time and memory resources. Therefore, the similarity threshold value is not directly set to 1 in the present application.

[0122] S40, load processing, according to the load information of the input working condition, find the working condition in S10 closest to the input working condition according to the set method, if it cannot be found, enter S10 on the basis of only updating the load information, otherwise enter S50;

[0123] The method for finding the closest working condition in S10 according to the load information of the input working condition in S40 is:

[0124] S41, input the load information of the input working condition as , , and are the load position and load size of the input working condition, and the load information of the working condition in S10 is , ;

[0125] S42, define a flag function foundsimilar and initialize it to false, when = but ≠ or ≠ but = When the flag function foundsimilar becomes true;

[0126] S43, compare the load information of the input working condition (S10) with each of the working conditions (S10), when the flag function foundsimilar is true and , , the corresponding working condition in S10 is collected and S44 is entered, otherwise the historical finite element model is called and only the load information is updated to enter S10; , When the flag function foundsimilar becomes true;

[0127] S44, the working conditions collected in S43 are screened to find the two working conditions closest to the input working condition; for But ≠ , the method for finding the two working conditions closest to the input working condition is:

[0128] Find such that ,

[0129] Find such that ,

[0130] The working conditions a and b are the two working conditions closest to the input working condition, , are the load sizes of working conditions a and b, respectively;

[0131] For ≠ But = , the method for finding the two working conditions closest to the input working condition is:

[0132] ,

[0133] ,

[0134] ,

[0135] ,

[0136] ,

[0137] ,

[0138] ​​The working conditions c and d are two working conditions closest to the input working condition, and Pc and Pd are load positions of the working conditions c and d, respectively, is the load position displacement of the working condition d to the working condition c, is the displacement of the working condition d to the input working condition, wherein , , are horizontal, vertical and vertical coordinates of the input working condition position, respectively, , , are horizontal, vertical and vertical coordinates of the working condition c position, respectively, , , are horizontal, vertical and vertical coordinates of the working condition d position, respectively;

[0139] t is an arbitrary value, if t does not exist, the historical finite element model is called, and only the load information is updated to enter S10; otherwise, the working conditions c and d are obtained by continuing to calculate through the index variables e and f

[0140] S50, simplification processing, if the input working condition and the working condition load information in S10 meet the set condition, interpolation calculation is performed, the interpolation calculation result is output, and the input working condition and the result are stored in the database table of the corresponding hoisting equipment; otherwise, only the load difference information is updated to enter S10.

[0141] The set condition in S50 is:

[0142] S51, if the load position of the working condition output in S40 is the same as that of the input working condition but the load size is different, S52 is entered, otherwise, the load size flow is entered;

[0143] S52, the working condition output in S40 is recorded as working condition a and working condition b, respectively, and the difference absolute value and are calculated, wherein is the load size of working condition a, is the load size of working condition b, , are the load size difference absolute values of working conditions a and b and the input working condition, respectively;

[0144] S53, according to the experiment, the load size error tolerance is set, if and , interpolation calculation is performed on the working condition a and b data to obtain the input working condition calculation result. In some preferred embodiments, node displacement and element stress calculation are taken as examples, and other results are the same.

[0145] ,

[0146] ,

[0147] where q, j are node and element serial numbers, Q, J are node and element numbers.

[0148] The load size procedure is as follows:

[0149] S1, for the working condition output in S40 and the working condition with the same load size but different load position, the working condition output in S40 is recorded as working condition c and working condition d respectively;

[0150] S2, the Euclidean distance of the load position of the two working conditions and the load position of the input working condition is calculated,

[0151] ,

[0152] ,

[0153] where, , are the distances of working condition c and d and the input working condition respectively;

[0154] S3, according to the experiment, the load position error tolerance is set , if and , the input working condition calculation result is obtained by interpolating the data of working condition c and d.

[0155] In some preferred embodiments, the interpolation calculation formula is as follows, taking node displacement and element stress calculation as an example, and the same applies to other results.

[0156] ,

[0157] ,

[0158] ,

[0159] where q, j are node and element serial numbers, Q, J are node and element numbers.

[0160] The setting of the load position error tolerance will comprehensively consider factors such as the structural characteristics and working environment of the equipment to ensure the applicability of the results in different situations. When storing and retrieving these data, for the stored position information, appropriate data structures such as spatial index (for geographic spatial data) will be used according to its spatial distribution characteristics to speed up data retrieval. For the calculated distance information, it will be stored in a temporary table or cache for subsequent comparison and judgment operations. If and Then, the input working condition is obtained by interpolating the data of working conditions c and d.

[0161] As shown in Figure 2 Another aspect of the present application relates to a bridge crane finite element calculation resource optimization system based on dynamic similarity discrimination, which comprises:

[0162] The database table setting and updating module sets various working conditions and finite element models according to practical data, receives input working conditions and updates the finite element models, performs finite element calculation on the working conditions, stores the working conditions, corresponding finite element models and finite element calculation results in the database table of the corresponding crane, and outputs the calculation results if the input working conditions are secondary input.

[0163] (1) Practical data application: Based on the data collected and sorted from actual engineering or experiments, various working conditions covering all types of working states of the bridge crane are set.

[0164] (2) Finite element calculation: Finite element calculation is performed on each working condition.

[0165] (3) Data storage: The detailed data of each working condition and the corresponding finite element calculation results are stored in the database table of the crane, providing reference for subsequent calculation, improving the calculation efficiency and performance of the system, and avoiding repeated calculation.

[0166] The equipment working condition searching module determines from which crane the input working condition comes from, calls the corresponding database table, and directly outputs the results if the input working condition exists in the database table, otherwise, enters S30.

[0167] (1) Equipment identification: By analyzing the input working condition, it is determined to which crane the working condition belongs. Different cranes have different structural characteristics and working conditions.

[0168] (2) Database table calling and result outputting: Once the equipment category is determined, the corresponding database table is called. If the data of the input working condition already exists in the database table, it means that the same or similar working condition has been calculated before, and the calculation results stored in the database table can be directly outputted to quickly obtain the required information and save calculation resources. If the input working condition is new and has not appeared in the database table, the next processing flow, i.e. the model geometric parameter checking and processing module (S30), is entered for further analysis and calculation.

[0169] The model geometric parameter checking and processing module enters S40 if the model geometric parameters of the input working condition exist in the database table, otherwise, enters S10.

[0170] The module mainly focuses on the model geometry parameters of the input working condition. It compares the geometry parameters of the input working condition with the geometry parameters of various working conditions stored in the database table. Through certain comparison methods, it is determined whether the geometry parameters of the input working condition are similar to the existing data in the database table. If similar geometry parameters can be found, it indicates that the existing calculation results can be used or a similar processing method can be used, thereby entering the load processing module (S40); if similar geometry parameters cannot be found, it indicates that the input working condition may involve new structural forms or sizes, and it is necessary to return to the database table setting module (S10) for new finite element calculation to obtain more accurate analysis results.

[0171] The load processing module finds the working condition in S10 closest to the input working condition according to the load information of the input working condition according to the setting method, and if it cannot be found, it enters S10 on the basis of only updating the load information, otherwise it enters S50;

[0172] According to the load information of the input working condition, the closest working condition in the S10 working condition in the database table is found according to the preset method. The load information can include various types. Here, "closest" is according to certain comparison logic, and factors such as the position, size, and distribution law of the load are considered. If a working condition close enough cannot be found in the database table, it indicates that the load situation of the input working condition is special, at which time only the load information is updated and the full calculation is returned to S10; if a close working condition is found, it enters the simplification processing module (S50) to use the existing information for subsequent processing.

[0173] The simplification processing module performs interpolation calculation if the load information of the input working condition and the working condition in S10 meets the set condition, outputs the interpolation calculation result, and stores the input working condition and the result to the database table corresponding to the hoisting equipment; otherwise, it enters S10 on the basis of only updating the load difference information.

[0174] Interpolation calculation: when the load information of the input working condition and the working condition in S10 meets the set condition, interpolation calculation is performed. Interpolation is a method of estimating unknown data points between known data points, and in the system, it is to estimate the result of the newly input working condition according to the existing working condition and its calculation result.

[0175] Those skilled in the art will readily understand that the above description is only a preferred embodiment of the present application and is not intended to limit the present application, and any modifications, equivalent replacements, and improvements made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A bridge crane finite element calculation resource optimization method based on dynamic similarity discrimination, characterized in that, The method comprises the following steps: S10, database table setting and updating, setting multiple working conditions and finite element models according to practical data, receiving secondary input working conditions and updating the finite element models, performing finite element calculation on the working conditions, and storing the working conditions, corresponding finite element models and finite element calculation results into the database table of the corresponding hoisting equipment, and if the input working conditions are secondary input, the calculation results are also outputted; S20, equipment working condition searching, judging from which hoisting equipment the input working conditions come from, calling the corresponding database table, if the input working conditions exist in the database table, the results are directly outputted, otherwise, S30 is entered; S30, model geometric parameter checking and processing, if the model geometric parameters of the input working conditions exist in the database table, S40 is entered, otherwise, S10 is entered; S40, load processing, finding the working conditions in S10 closest to the input working conditions according to the load information of the input working conditions, if the working conditions cannot be found, S10 is entered on the basis of only updating the load information, otherwise, S50 is entered; S50, simplification processing, if the input working conditions and the working conditions in S10 meet the set conditions, interpolation calculation is performed, the interpolation calculation results are outputted, and the input working conditions and the results are stored into the database table of the corresponding hoisting equipment; otherwise, S10 is entered on the basis of only updating the load difference information.

2. The dynamic similarity-based discrimination method for optimizing finite element calculation resources of bridge hoisting equipment according to claim 1, characterized in that, The method for setting multiple working conditions and finite element models in S10 and receiving secondary input working conditions and updating the finite element models, and performing finite element calculation on the working conditions is as follows: The total number of the working conditions is N, i is the working condition serial number, The i-th group of the working conditions comprises a model geometry parameter set , a load information set and a boundary condition set Bi, the model geometry parameter set comprises M model geometry parameters, the total number of the model geometry parameters M is set according to experiments, m is the model geometry parameter serial number, The load information set comprises a load position and a load size The boundary condition set Bi comprises L boundary condition parameters, the total number of the boundary condition parameters L is set according to experiments, l is the boundary condition parameter serial number, l = 1, 2, …, L; for the dynamic hoisting equipment: ， ， ， wherein, is a mass matrix for the i-th group of conditions, is a device movement velocity, is a device movement acceleration, is a damping matrix, is a stiffness matrix, is a device displacement, t is time, is an external load vector as a function of time t; For the density of the device, the function is a function for calculating the mass matrix, the function g is a load calculation function; For static hoisting equipment: ， Equivalent nodal force vector for the i-th group of conditions.

3. The dynamic similarity-based discrimination method for optimizing finite element calculation resources of bridge hoisting equipment according to claim 2, characterized in that, The method for storing the working conditions, corresponding finite element models and finite element calculation results into the database table of the corresponding hoisting equipment in S10 is as follows: The database table includes working conditions, finite element models and finite element calculation result fields, a database table row is added when storing new data, and working condition summary collections are used , finite element calculation result summary collections The working conditions and finite element calculation results are summarized respectively, and the two summary collections are stored in cells corresponding to the fields of the new row of the database table in a set format. , , wherein, , , are respectively the geometry parameter set, the load information set, the boundary condition set of the i-th group of working conditions, , , , are respectively the node displacement set, the element axial stress set, the element combined stress set, the element safety factor set of the i-th group of finite element calculation results; The finite element models set in S10 and the updated finite element models generate unique codes when performing calculation, the unique codes are added to the model names, and the finite element models are stored separately, the addresses and names of the finite element models are stored into the cells of the corresponding fields of the newly added rows of the database table.

4. The dynamic similarity-based discrimination method for optimizing finite element calculation resources of bridge hoisting equipment according to claim 3, characterized in that, The method for judging from which hoisting equipment the input working conditions come from in S20 is as follows: The input working conditions are represented by an input set I: , In the formula, , , respectively represent the input working condition geometry parameter set, load information set, boundary condition set; There are F hoisting equipments, each hoisting equipment database table takes a group of working conditions, a total of F groups of working conditions, f is the working condition serial number of the hoisting equipment, the input working conditions are compared with the F groups of working conditions in similarity: ， , , In the formula, is a similarity measurement function of the input working condition and the fth group of working conditions, is a geometric parameter weight coefficient, is a load information weight coefficient, is a boundary condition weight coefficient; is a geometric similarity function of the input working condition and the fth group of working conditions, is a load similarity function of the input working condition and the fth group of working conditions, is a boundary condition similarity function of the input working condition and the fth group of working conditions, and a similarity threshold value of the fth group of working conditions is set according to an experiment , if , the input working condition comes from the fth group of working conditions corresponding to the hoisting equipment.

5. The dynamic similarity-based discrimination method for optimizing finite element calculation resources of bridge hoisting equipment according to claim 4, characterized in that, The method for model geometric parameter checking and processing in S30 is as follows: , ， In the formula, is the mth model geometry parameter of the input working condition, is the mth model geometry parameter of the i th group of working conditions in S10, is a geometry similarity degree measurement function; a geometry similarity degree threshold value is obtained according to experiments ; If there exists i=k such that , the model geometry parameters of the input working condition have a similar match in the database table, then go to S40; if all satisfy , the model geometry parameters of the input working condition have no similar match in the database table, otherwise, the historical finite element model is called, the model geometry parameters and the load information are updated, and then go to S10.

6. The dynamic similarity-based discrimination method for optimizing finite element calculation resources of bridge hoisting equipment according to claim 5, characterized in that, The method for finding the working conditions in S10 closest to the input working conditions according to the load information of the input working conditions in S40 is as follows: S41, input the load information of the input working condition as , ), and are the load position and load size of the input working condition respectively, the working condition load information in the S10 is , ); S42, define a flag function foundsimilar and initialize it to false, when a = But ≠ Or ≠ But = the flag function foundsimilar becomes true; S43, input the load information of the working condition ( , ) and the working condition of each of the S10 ( , The comparison is performed when both the flag function `foundsimilar` and `foundsimilar` are true. When the condition is not in the specified state, the corresponding working condition in S10 is collected and the process proceeds to S44; otherwise, the historical finite element model is retrieved and the process proceeds to S10 with only the load information updated. S44, screening the working conditions collected in S43 to find two groups of working conditions closest to the input working condition; for = But ≠ The method for finding two groups of working conditions closest to the input working condition is: find such that , find such that , The working conditions a and b are two groups of working conditions closest to the input working condition, , are the load sizes of the working conditions a and b, respectively. For ≠ But = The method for finding two groups of working conditions closest to the input working condition is: ， ， ， , ， , The working conditions c and d are two groups of working conditions closest to the input working condition; Pc and Pd are load positions of the working conditions c and d, respectively, is the displacement of the load position from the working condition d to the working condition c, is the displacement from the working condition d to the input working condition, wherein , , are the horizontal, vertical and vertical coordinates of the input working condition position, respectively, , , are the horizontal, vertical and vertical coordinates of the working condition c position, respectively, , , are the horizontal, vertical and vertical coordinates of the working condition d position, respectively. t is an arbitrary value, if t does not exist, the historical finite element models are called, S10 is entered on the basis of only updating the load information; otherwise, the working conditions c and d are obtained through the index variables e and f.

7. The dynamic similarity-based discrimination method for optimizing finite element calculation resources of bridge hoisting equipment according to claim 6, characterized in that, The set conditions in S50 are as follows: S51, if the working conditions outputted in S40 and the input working conditions have the same load position but different load sizes, S52 is entered, otherwise, the load size flow is entered; S52, the S40 output respectively as the working condition a and working condition b, calculate the difference absolute value and Wherein The load size of working condition a, The load size of working condition b, , The absolute value of the difference between the load size of working condition a and b and the input working condition, respectively; S53, according to the experiment, the load size error tolerance is set , if and , the interpolation calculation of working condition a and b data is carried out to obtain the input working condition calculation result.

8. The dynamic similarity-based discrimination method for optimizing finite element calculation resources of bridge hoisting equipment according to claim 7, characterized in that, The load size flow is as follows: S1, for the S40 output conditions and the load size of the input conditions of the different load position of the working condition, the S40 output conditions are recorded as working condition c and working condition d respectively; S2, calculate the Euclidean distance of the load position of the two working conditions and the input working condition load position, ， ， wherein , are the distances of the operating conditions c and d, respectively, from the input operating condition. S3, according to the experimental load position error tolerance setting , if and then the input operating condition calculation result is obtained by interpolating the operating condition c and d data.

9. A bridge crane finite element calculation resource optimization system based on dynamic similarity discrimination, which adopts the method according to claim 1, characterized in that, Comprise: Database table setting and updating module, according to the practice data setting multiple working conditions and finite element model and receiving secondary input working condition and updating finite element model, the working condition is calculated by finite element, and the working condition, corresponding finite element model and finite element calculation result are stored to the database table of corresponding hoisting equipment, if it is secondary input working condition, then it also needs to output calculation result; Equipment working condition searching module, judge input working condition from which hoisting equipment, call corresponding database table, if the input working condition exists in the database table, then directly output result, otherwise enter S30; Model geometric parameter inspection and processing module, if the model geometric parameter of the input working condition exists in the database table, then enter S40, otherwise enter S10; Load processing module, according to the load information of the input working condition, find out the working condition in S10 closest to the input working condition according to the setting method, if it cannot be found, then enter S10 on the basis of only updating load information, otherwise enter S50; Simplified processing module, if the input working condition and the working condition in S10 load information meet the setting condition, then interpolation calculation is carried out, interpolation calculation result is output, and the input working condition and result are stored to the database table of corresponding hoisting equipment, otherwise enter S10 on the basis of only updating load difference information.

Images