A coal yard steel structure simulation analysis method

Abstract

The present application relates to steel structure simulation technical field, disclose a kind of coal yard steel structure simulation analysis method, comprising: collecting coal mine field data constructs coal mine field steel structure geometric model;Grid division is carried out to steel structure geometric model, and material attribute in grid interior is set, and coal yard steel structure simulation model is generated;Reliability evaluation is carried out to coal yard steel structure simulation model based on historical data, and reliability evaluation result is generated;Coal yard steel structure simulation model is optimized based on reliability evaluation result.The present application carries out grid division to steel structure geometric model, calculates the complexity score C of first-level grid element, and different grid division strategies are executed based on this.This is helpful to reasonably divide grid according to the actual complexity of grid element, so that the generated coal yard steel structure simulation model has suitable grid density and accuracy in different regions.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of steel structure simulation, in particular to a coal yard steel structure simulation analysis method. BACKGROUND

[0002] The design and maintenance of coal yard steel structures are crucial for ensuring the safety of coal storage and transportation. With the development of mining technology and the increasing demand for environmental protection, the design of coal yard structures has become more complex, requiring consideration of factors such as structural stability, durability, and adaptability to environmental changes. Traditional design methods often rely on experience and simplified calculation models, which may not accurately predict the performance of structures under actual working conditions.

[0003] Collecting field data to construct geometric models, performing mesh division, and setting material properties are critical steps. These steps ensure that the simulation model accurately reflects the characteristics of the actual structure. Subsequently, reliability assessment of the model based on historical data can evaluate the performance of the structure under expected loads and environmental conditions. Finally, model optimization based on the evaluation results can significantly improve the safety and economy of the structure. SUMMARY

[0004] The purpose of the present application is to collect field data of coal mines to construct a geometric model of a coal mine steel structure, perform mesh division, and set material properties to generate a simulation model. Based on historical data, the reliability of the generated coal yard steel structure simulation model is evaluated, including structural reliability evaluation and response reliability evaluation. The weighted sum of the two results is obtained to obtain the reliability evaluation result, and the simulation model is optimized based on the result.

[0005] To achieve the above purpose, the present application provides a coal yard steel structure simulation analysis method, comprising:

[0006] Collecting field data of coal mines to construct a geometric model of a coal mine steel structure;

[0007] Mesh division is performed on the steel structure geometric model, and the material properties inside the mesh are set to generate a coal yard steel structure simulation model;

[0008] Based on historical data, the reliability of the coal yard steel structure simulation model is evaluated to generate a reliability evaluation result;

[0009] Based on the reliability evaluation result, the coal yard steel structure simulation model is optimized.

[0010] In some embodiments of the present application, when mesh division is performed on the steel structure geometric model, it comprises:

[0011] A steel structure geometric model is constructed by defining the size of the nodes, the size of the bars, and the shape of the bars;

[0012] performing preliminary stress analysis on the steel structure geometric model, performing preliminary mesh division on the steel structure geometric model according to the result of the preliminary stress analysis, and generating a plurality of first-level mesh units;

[0013] setting initial material properties of the first-level mesh units based on historical data;

[0014] calculating a complexity score C of the first-level mesh units;

[0015] and performing a mesh division strategy based on the complexity score C to generate second-level mesh units.

[0016] In some embodiments of the present application, when the plurality of first-level mesh units are generated, the method comprises:

[0017] performing simulation tests on the steel structure geometric model based on the obtained historical environmental data of the coal yard steel structure;

[0018] setting monitoring nodes at nodes of the steel structure geometric model and at centers of the members, and applying simulation forces to the steel structure geometric model to generate a stress data set D1 at the monitoring nodes, D1 = [d1, d2…di…dn];

[0019] wherein di represents the stress size of the i-th monitoring node, and n represents the number of monitoring nodes;

[0020] performing first-level mesh unit division on the stress data set D1 at the monitoring nodes to generate a stress data set D2 at the remaining monitoring nodes;

[0021] continuing to perform first-level mesh unit division on the stress data set D2 at the remaining monitoring nodes until all monitoring nodes are divided into first-level mesh units;

[0022] generating a plurality of first-level mesh units, including: M1, M2…Mi…Mm;

[0023] wherein Mi represents the i-th first-level mesh unit, m is the number of generated first-level mesh units, and m < n;

[0024] The first-level mesh unit division comprises:

[0025] setting a first deviation value a; and calculating a stress deviation value A of d1 and all adjacent monitoring nodes of d1;

[0026] selecting all adjacent monitoring nodes of A < a, then dividing d1 and the selected all adjacent monitoring nodes into a first-level mesh unit; and removing a plurality of monitoring nodes contained in the generated first-level mesh unit from the stress data set D1 at the monitoring nodes to generate a stress data set D2 at the remaining monitoring nodes.

[0027] In some embodiments of the present application, when the deviation value A of the calculation master monitoring node and the adjacent monitoring node is calculated, the following steps are included:

[0028] The comprehensive deviation value ΔFd in the calculation direction is calculated as follows:

[0029] The force borne by the i-th monitoring node and the force borne by the adjacent monitoring node are decomposed into x, y, and z direction components;

[0030] The difference values of the x, y, and z direction components are calculated respectively to generate ΔFx, ΔFy, and ΔFz;

[0031] Wherein, ΔFx represents the difference value between the x direction component of the i-th monitoring node and the x direction component of the adjacent monitoring node; ΔFy represents the difference value between the y direction component of the i-th monitoring node and the y direction component of the adjacent monitoring node; and ΔFz represents the difference value between the z direction component of the i-th monitoring node and the z direction component of the adjacent monitoring node.

[0032] The comprehensive deviation value ΔFd in the direction is generated by vector synthesis;

[0033]

[0034] The deviation value |di-d(i+j)| of the force magnitude is calculated;

[0035] Wherein, d(i+j) represents the force value of the adjacent monitoring node, and j represents the adjacent distance from the master monitoring node.

[0036] A=k1*ΔFd+k2*|di-d(i+j)|, (k1+k2=1);

[0037] k1 is the weight of the comprehensive deviation value in the direction, and k2 is the weight of the deviation value of the force magnitude.

[0038] In some embodiments of the present application, when the primary grid cell complexity score C is calculated, the following steps are included:

[0039] The number q of monitoring nodes contained in the current primary grid cell is obtained, and the aspect ratio r of the current primary grid cell and the force distribution degree σ of the current primary grid cell are calculated to generate the primary grid complexity score C;

[0040] C=w1×q+w2×r+w3×σ;

[0041] Wherein, w1 is the weight of the number q of monitoring nodes contained in the current primary grid cell, w2 is the weight of the aspect ratio r of the current primary grid cell, and w3 is the weight of the force distribution degree σ of the current primary grid cell, and w1+w2+w3=1, the values of w1, w2, and w3 are adjusted according to historical data.

[0042] In some embodiments of the present application, when the aspect ratio r in the primary grid unit and the stress distribution degree σ in the primary grid are calculated, the method comprises:

[0043] Obtaining each internal angle α in the primary grid unit, α=(α1, α2, …, αn);

[0044] Wherein, αn represents the nth internal angle;

[0045] Calculating the standard deviation σ of the internal angle in the primary grid unit α ;

[0046] Based on the standard deviation σ of the internal angle in the primary grid unit α , a plurality of internal angle standard deviation thresholds are set, and the aspect ratio r of the current primary grid unit is generated based on the internal angle standard deviation threshold;

[0047] Calculating the standard deviation σ of the stress f , based on the standard deviation σ of the stress f , a plurality of stress standard deviation thresholds are set, and the stress distribution degree σ of the current primary grid unit is generated based on the stress standard deviation threshold.

[0048] In some embodiments of the present application, when the grid division strategy is executed based on the complexity score C to generate the secondary grid unit, the method comprises:

[0049] Setting a complexity threshold C1 based on the complexity score C;

[0050] If C≤C1, the complexity level of the current primary grid unit is the first complexity,

[0051] Executing a primary grid division strategy;

[0052] If C1<C, the complexity level of the current primary grid unit is the second complexity, and a high-level grid division strategy is executed;

[0053] In some embodiments of the present application, the grid division strategy comprises:

[0054] The primary grid division strategy comprises: performing equal proportion division based on the shape of the current primary grid unit to generate the secondary grid unit;

[0055] The high-level grid division strategy comprises: calculating the stress distribution degree σ of the nodes in the current primary grid unit; and dividing the current primary grid unit based on the stress distribution degree σ inside the current primary grid unit to generate the secondary grid unit.

[0056] In some embodiments of the present application, the simulation model of the coal yard steel structure is optimized, which comprises:

[0057] Re-correct the initial material properties in the secondary grid unit, generate a coal yard steel structure simulation model;

[0058] Based on the obtained coal yard historical environment data, reliability evaluation is performed on the coal yard steel structure simulation model to generate a reliability evaluation result;

[0059] And based on the reliability evaluation result, the coal yard steel structure simulation model is optimized.

[0060] In some embodiments of the present application, when the reliability evaluation result is generated, it includes:

[0061] Structural reliability evaluation and response reliability evaluation;

[0062] The structural reliability evaluation includes:

[0063] Performing static analysis on the coal yard steel structure simulation model to determine the stability value of the coal yard steel structure simulation model under static load;

[0064] Determine the natural frequency and modal shape of the secondary grid unit through modal analysis to generate the resonance risk value of the coal yard steel structure simulation model;

[0065] Based on the stability value and the resonance risk value of the coal yard steel structure simulation model, a structural reliability evaluation value is generated;

[0066] The response reliability evaluation includes:

[0067] Based on historical data, seismic simulation is performed on the coal yard steel structure simulation model to calculate the residual safety reserve of the structure, and a first response reliability evaluation value is generated;

[0068] Based on historical data, the settlement process of the foundation is simulated by gradually applying load, the change of the stress field of the coal yard steel structure simulation model is calculated, and a second response reliability evaluation value is generated;

[0069] The first response reliability evaluation value and the second response reliability evaluation value are combined to generate a response reliability evaluation value;

[0070] The structural reliability evaluation value and the response reliability evaluation value are weighted and summed to generate a reliability evaluation result.

[0071] Compared with the prior art, the coal yard steel structure simulation analysis method according to the embodiments of the present application has the following advantages:

[0072] The steel structure geometry model is meshed by calculating the complexity score of the primary mesh unit, and different meshing strategies are executed, which is helpful to reasonably divide the mesh according to the actual complexity of the mesh unit, such as the number of nodes, the aspect ratio, the stress distribution degree and other factors, so that the generated simulation model of the coal yard steel structure has appropriate mesh density and accuracy in different areas.

[0073] The primary mesh unit is divided by setting the first deviation value and calculating the deviation value of the main monitoring node and the adjacent monitoring node, so that the size and shape of the mesh unit can be flexibly adjusted according to the stress change of the structure.

[0074] By comparing the deviation value, the number of mesh units can be reduced in the area with relatively uniform stress, and more detailed division can be performed in the area with larger stress change, so as to improve the calculation efficiency, ensure the relative consistency of the stress in the mesh unit, and improve the rationality of the mesh division.

[0075] The deviation value A is calculated by combining the direction deviation and the size deviation, so that the calculation of the deviation value is more comprehensive, the boundary of the mesh division can be more accurately determined, and the primary mesh unit divided has better consistency in stress characteristics.

[0076] By comprehensively evaluating the mesh complexity from multiple aspects, the complexity of the primary mesh unit can be comprehensively and accurately described, and the one-sidedness of evaluating the mesh complexity by relying on a single factor is avoided.

[0077] By dividing the secondary mesh unit, reasonable mesh division can be performed according to the actual situation of each part, so that the simulation model can more truly reflect the actual working state of the coal yard steel structure, and provide more reliable basis for engineering design, analysis and optimization.

[0078] By re-calibrating the material properties, the material characteristics in the model can be more matched with the actual structure, the accuracy of the simulation model is improved, and the mechanical behavior of the coal yard steel structure under actual working conditions can be more truly simulated.

[0079] The structure reliability evaluation starts from two aspects of determining the stability value from static analysis and determining the resonance risk value from modal analysis, considers the stability under static load and the self-vibration characteristics of the structure, and comprehensively evaluates the reliability of the structure itself.

[0080] The response reliability evaluation calculates the residual safety reserve through seismic simulation and calculates the stress field change of the structure through simulation of foundation settlement, and considers the influence of external dynamic load (earthquake) and foundation settlement on the structure. BRIEF DESCRIPTION OF DRAWINGS

[0081] Figure 1 It is a coal yard steel structure simulation analysis method flowchart provided by the embodiment of the application. DETAILED DESCRIPTION

[0082] The specific embodiments of the present application will be further described in conjunction with the drawings and examples. The following examples are used to illustrate the present application, but are not used to limit the scope of the present application.

[0083] In the description of the present application, it needs to be understood that the terms "center", "upper", "lower", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer" and the like indicate the orientation or positional relationship based on the orientation or positional relationship shown in the drawings, and are only for the purpose of facilitating the description of the present application and simplifying the description, and do not indicate or imply that the device or element referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as a limitation of the present application.

[0084] The terms "first", "second" are only for the purpose of description, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the technical features indicated. Therefore, the features defined with "first", "second" can explicitly or implicitly include one or more of the features. In the description of the present application, unless otherwise specified, the meaning of "a plurality of" is two or more.

[0085] In the description of the present application, it needs to be explained that, unless otherwise explicitly specified and limited, the terms "mounting", "connecting", "connection" should be understood in a broad sense, for example, it can be fixed connection, or detachable connection, or integral connection; it can be mechanical connection, or electrical connection; it can be direct connection, or indirect connection through intermediate medium, or internal communication of two elements. For those skilled in the art, the specific meaning of the above terms in the present application can be understood according to the specific circumstances.

[0086] Example 1:

[0087] A coal yard steel structure simulation analysis method according to an embodiment of the present application, as shown in Figure 1 , comprises:

[0088] Collecting coal mine site data to construct a coal mine steel structure geometric model;

[0089] Grid division is performed on the steel structure geometric model, and material properties inside the grid are set to generate a coal yard steel structure simulation model;

[0090] Based on historical data, the reliability of the coal yard steel structure simulation model is evaluated to generate a reliability evaluation result;

[0091] Based on the reliability evaluation result, the coal yard steel structure simulation model is optimized.

[0092] In the process of meshing the steel structure geometric model, the following steps are included:

[0093] A steel structure geometric model is constructed by defining the size of nodes, the size of members, and the shape of members.

[0094] A preliminary stress analysis is performed on the steel structure geometric model, and based on the results of the preliminary stress analysis, a preliminary meshing is performed on the steel structure geometric model to generate a plurality of first-level mesh elements.

[0095] The initial material properties of the first-level mesh elements are set based on historical data.

[0096] The complexity score C of the first-level mesh elements is calculated.

[0097] Based on the complexity score C, a meshing strategy is executed to generate second-level mesh elements.

[0098] In this embodiment, the node size is defined as

[0099] First, the size of the node is determined according to the actual design requirements and the expected stress condition of the steel structure. For the nodes of the main bearing structure, considering the stress concentration and force transmission requirements, the node size should ensure sufficient connection strength. For example, at the beam-column connection node, according to the section size of the beam and column, the length, width and height of the node are set to 1.5 times the width of the beam, 1.2 times the width of the column, and the thickness value determined according to the connection structure.

[0100] At the same time, refer to relevant steel structure design specifications, such as "Steel Structure Design Standard" (GB 50017-2017), to ensure that the node size is within a reasonable range of values to meet the safety and reliability requirements of the structure.

[0101] Member size determination

[0102] The stress type of the members in the steel structure is analyzed, such as axial stress, bending shear combined stress, etc. For members subjected to axial compression, the cross-sectional size of the member is determined according to the Euler formula and the stability requirements of the structure. For example, the critical pressure of the member is calculated, and then according to the actual axial pressure value, a certain safety factor (such as 1.5-2.0) is reserved to determine the diameter or side length of the member.

[0103] For members subjected to bending shear combined stress, the combined action of bending moment and shear force is considered. According to the bending stress and shear stress calculation formula in material mechanics, combined with the allowable stress of steel, the cross-sectional shape and size of the member are determined. For example, when using a I-shaped cross section, the width and thickness of the flange plate are determined according to the bending moment, and the height and thickness of the web are determined according to the shear force.

[0104] Member shape design

[0105] According to the overall layout and force transmission path of the steel structure, the appropriate shape of the rod is selected. In the space force system, for the rod that needs to transmit multi-directional force, box section can be used to improve the torsional stiffness of the rod.

[0106] In the plane force system, such as the rod in the truss structure, according to the force direction, round or rectangular section can be used. For the rod mainly bearing tension, round section is more reasonable in material utilization; for the rod bearing pressure and bending moment, rectangular section can better meet the stability requirement.

[0107] Primary grid cell initial material property setting

[0108] Historical data sources

[0109] From the design documents of the steel structure, material test reports and the database of similar steel structure engineering in the past, the historical data of the steel material, including yield strength, elastic modulus, Poisson's ratio and other material properties, are obtained.

[0110] For special steel or new type of steel, special material test may be needed to obtain accurate material property data.

[0111] Property setting

[0112] According to the obtained historical data, the material properties of the primary grid cells are assigned to the corresponding cells. For example, for ordinary Q235 steel, the yield strength is set to 235 MPa, the elastic modulus is set to 206 GPa, and the Poisson's ratio is set to 0.3.

[0113] Considering the actual dispersion of the material, when setting the material properties, a certain fluctuation range can be introduced according to the statistical data. For example, the yield strength can fluctuate within ±10% to more truly simulate the stress performance of the structure.

[0114] In embodiment 3, when the plurality of primary grid cells are generated, the method comprises:

[0115] Based on the obtained historical environmental data of the coal yard steel structure, a simulation test is performed on the steel structure geometric model;

[0116] Monitoring nodes are set at the nodes of the steel structure geometric model and the centers of the rods, and a simulation force is applied to the steel structure geometric model to generate a stress data set D1 at the monitoring nodes, D1 = [d1, d2…di…dn];

[0117] Wherein, di represents the stress size of the i th monitoring node, and n represents the number of monitoring nodes;

[0118] The stress data set D1 at the monitoring nodes is divided into primary grid cells to generate a stress data set D2 at the remaining monitoring nodes;

[0119] Continue the first-level grid cell division on the stress data set D2 at the remaining monitoring nodes until all monitoring nodes are divided into first-level grid cells;

[0120] Generate a plurality of first-level grid cells, including: M1, M2…Mi…Mm;

[0121] Wherein, Mi represents the i-th first-level grid cell, m is the number of generated first-level grid cells, m < n;

[0122] The first-level grid cell division includes:

[0123] Set the first deviation value a; Calculate the stress deviation value A of d1 and all adjacent monitoring nodes of d1;

[0124] Select all adjacent monitoring nodes with A < a, then divide d1 and the selected all adjacent monitoring nodes into a first-level grid cell; and remove the plurality of monitoring nodes contained in the generated first-level grid cell in the stress data set D1 at the monitoring nodes; generate the stress data set D2 at the remaining monitoring nodes.

[0125] In this embodiment,

[0126] Obtain historical stress data of the steel structure from the operation records of the coal yard. These data include the pressure, lateral force and other forces borne by each part of the steel structure under different coal pile heights. For example, under different working conditions such as coal pile height of 5 meters, 10 meters, 15 meters, etc., the stress of the steel structure at the corresponding time is recorded.

[0127] At the same time, consider the stress data of the coal yard under special circumstances (such as strong wind weather, local collapse of coal pile, etc.). Extract real-time stress data under these special working conditions from the monitoring system (if any) to ensure the completeness and diversity of the data.

[0128] Data preprocessing

[0129] Screen the collected historical stress data to remove obviously incorrect or abnormal data points. For example, due to sensor failure, the stress value is too large or too small.

[0130] Classify and organize the data according to different stress types (such as vertical pressure, horizontal thrust, etc.) and structure parts (such as columns, beams, supports, etc.) for subsequent analysis and use.

[0131] Finite element simulation settings

[0132] Use professional finite element analysis software (such as ABAQUS) to simulate the stress test of the steel structure geometric model. Apply the sorted historical stress data as external load to the steel structure geometric model.

[0133] According to the actual connection mode and constraint conditions of the steel structure, the constraints of the nodes are accurately set in the finite element model. For example, for the column connected with the foundation at the bottom, it is set as a fixed end constraint; for the connection nodes of beams and columns, the corresponding constraint settings are made according to the actual hinge or connection conditions.

[0134] Test process and result analysis

[0135] The load corresponding to the historical stress data is gradually applied to simulate the stress process of the steel structure in the actual coal yard operation process. The deformation of the steel structure in the loading process is observed, including displacement, angle, etc.

[0136] The stress distribution of each part of the steel structure is analyzed. The stress concentration area is intuitively viewed through the stress nephogram, such as the stress size and distribution law of the beam-column connection node, the support and the main structure connection part, etc., which provides a basis for the subsequent setting of monitoring nodes.

[0137] Stress data collection

[0138] After setting the monitoring nodes, the stress data of each monitoring node is collected during the stress test process through the finite element analysis software. For each monitoring node, the stress size under different load steps is recorded.

[0139] According to the numbering order of the nodes, the collected stress data is arranged into a data set D1=[d1,d2…di…dn]. Wherein, di represents the stress size of the i th monitoring node, and n represents the number of monitoring nodes.

[0140] Data accuracy verification

[0141] The accuracy of the generated stress data set is verified. By comparing the stress relationship of different monitoring nodes under the same load condition, it is checked whether the data conforms to the basic principles of structural mechanics.

[0142] For example, for two adjacent monitoring nodes, if they are on the same member, according to the force balance equation of the member, the difference between their stresses should be within a reasonable range. If unreasonable data is found, recheck the setting of the monitoring nodes and the data collection process to ensure the accuracy of the data.

[0143] Consideration of structural stress characteristics

[0144] According to the stress characteristics of the coal yard steel structure, the first deviation value a is set. If the steel structure mainly bears relatively uniform pressure load, such as the vertical pressure of the coal pile on the structure, a can be set relatively large. For example, for this case, a can be set to 10 kN (assuming the unit of stress is kN).

[0145] If the steel structure is also subjected to large uneven loads or dynamic loads (such as lateral forces caused by strong winds or loads caused by coal yard equipment vibration), in order to more finely divide the grid to accurately capture the stress changes, a should be set smaller. For example, in this case, a can be set to 5 kN.

[0146] Structural accuracy requirements

[0147] Adjust the value of a considering the accuracy requirements for steel structure analysis. If high-precision analysis results are required, such as stress analysis of critical parts of the structure (such as support structures, areas near nodes), a should be set smaller to ensure that the grid division can reflect smaller stress deviations.

[0148] If it is a preliminary overall analysis of the entire steel structure, the accuracy requirement is relatively low, a can be appropriately increased to reduce the amount of calculation.

[0149] In the embodiment 4, the deviation value A of the main monitoring node and the adjacent monitoring node is calculated, comprising:

[0150] The comprehensive deviation value ΔFd in the calculation direction is calculated, comprising:

[0151] The force borne by the i-th monitoring node and the force borne by the adjacent monitoring node are decomposed into x, y, z direction components;

[0152] The difference values of the x, y, z direction components are calculated respectively to generate ΔFx, ΔFy and ΔFz;

[0153] Wherein, ΔFx represents the difference value between the x direction component of the i-th monitoring node and the x direction component of the adjacent monitoring node; ΔFy represents the difference value between the y direction component of the i-th monitoring node and the y direction component of the adjacent monitoring node; ΔFz represents the difference value between the z direction component of the i-th monitoring node and the z direction component of the adjacent monitoring node;

[0154] The comprehensive deviation value ΔFd in the direction is generated by vector synthesis;

[0155]

[0156] The deviation value of the force size |di-d(i+j)| is calculated;

[0157] Wherein, d(i+j) represents the force value of the adjacent monitoring node, and j represents the adjacent distance from the main monitoring node.

[0158] A=k1*ΔFd+k2*|di-d(i+j)|,(k1+k2=1);

[0159] k1 is the weight of the comprehensive deviation value in the direction, and k2 is the weight of the deviation value of the force size.

[0160] Example 5: The calculation of the first-level grid cell complexity score C includes:

[0161] Obtain the number of monitoring nodes q contained in the current first-level grid cell, calculate the aspect ratio r and the stress distribution degree σ of the current first-level grid cell, and generate the complexity score C of the first-level grid cell;

[0162] C = w1 × q + w2 × r + w3 × σ;

[0163] Where w1 is the weight of the number of monitoring nodes q in the current first-level grid cell, w2 is the weight of the aspect ratio r of the current first-level grid cell, and w3 is the weight of the force distribution degree σ of the current first-level grid cell, and w1+w2+w3=1, the values ​​of w1, w2 and w3 are adjusted according to historical data.

[0164] In this embodiment, the values ​​of w1, w2, and w3 are adjusted based on historical data. This characteristic allows the scoring method to adapt to different types of coal yard steel structures. For example, for some coal yard steel structures with relatively regular shapes and relatively uniform stress distribution, the weight of the stress distribution degree σ, w3, may be appropriately reduced, while for structures with complex shapes and large stress variations, the value of w3 can be increased. This flexible adjustment of weights based on historical data allows the complexity score C of the first-level grid unit to better reflect the actual situation of different structures, improving the accuracy and applicability of complexity assessment.

[0165] Example 6: When calculating the aspect ratio r and the force distribution σ in the first-level grid cell, the following steps are included:

[0166] Obtain the interior angles α in the first-level grid cell, where α = (α1, α2, ..., αn);

[0167] Where αn represents the nth interior angle;

[0168] Calculate the standard deviation σ of the interior angles in a first-level mesh cell. α ;

[0169] Based on the standard deviation σ of the interior angles in the first-level grid cell α Set multiple interior angle standard deviation thresholds, and generate the aspect ratio r of the current first-level grid cell based on the interior angle standard deviation thresholds;

[0170] Calculate the standard deviation σ of the force. f Based on the standard deviation of the force σ f Multiple stress standard deviation thresholds are set, and the stress distribution degree σ of the current first-level grid cell is generated based on the stress standard deviation thresholds.

[0171] In the embodiment 7, the grid division strategy based on the complexity score C includes:

[0172] Setting a complexity threshold C1 based on the complexity score C;

[0173] If C≤C1, the complexity level of the current primary grid cell is primary complexity, and the primary grid division strategy is executed;

[0174] If C1<C, the complexity level of the current primary grid cell is secondary complexity, and the advanced grid division strategy is executed;

[0175] In the embodiment, the complexity threshold C1 is set according to historical data. When the complexity score C is lower than the complexity threshold C1, the grid only needs to be divided into secondary grid cells in equal proportions, and then the stress analysis of the secondary grid cells is performed. At this time, it is primary complexity.

[0176] When the complexity score C is higher than the complexity threshold C1, the accuracy of stress analysis of the secondary grid cells generated by direct equal proportion division will be reduced due to the increase of complexity, so the primary grid after equal proportion division needs to be analyzed again. At this time, it is secondary complexity.

[0177] Through statistical analysis of historical data, a threshold that can distinguish different complexity levels can be determined, or through a series of experiments, the threshold is adjusted until a balance point that can guarantee analysis quality and maintain efficiency is found.

[0178] In the embodiment 8, the grid division strategy includes:

[0179] The primary grid division strategy includes: dividing the current primary grid cell in equal proportions based on the shape of the current primary grid cell to generate secondary grid cells;

[0180] The advanced grid division strategy includes: calculating the stress distribution degree σ of the nodes in the current primary grid cell; and dividing the current primary grid cell based on the stress distribution degree σ of the current primary grid cell to generate secondary grid cells.

[0181] In the embodiment, the shape of the primary grid cell is determined

[0182] Firstly, the shape of the primary grid cell is identified. If it is a triangular cell, the length proportion relationship of its three sides needs to be determined; if it is a quadrilateral cell, the length and internal angle relationship needs to be determined, etc. For example, for a quadrilateral primary grid cell, the lengths of its four sides are a=5 units, b=5 units, c=5 units, and d=5 units (assuming the unit is centimeter), and the four internal angles are close to 90 ° , which can be regarded as a square-shaped cell.

[0183] According to the shape of the unit, an equal proportion is divided. For the above-mentioned square shape of the first-level grid unit, if it is to be divided into smaller second-level grid units, it can be divided according to the same edge length proportion. For example, each edge is divided into n = 2 parts (here n can be set according to actual needs).

[0184] (I) Primary grid division strategy execution

[0185] The primary grid division strategy is executed on the current first-level grid unit. According to the primary grid division strategy described above, an equal proportion is divided according to the shape of the first-level grid unit. For example, for a hexagonal first-level grid unit, each edge is divided into n = 3 parts, obtaining a plurality of small units after preliminary division, which constitute the first-level grid unit after preliminary processing.

[0186] (II) Advanced grid division strategy execution

[0187] The advanced grid division strategy is executed on the first-level grid unit after preliminary processing. First, the stress distribution degree σ of the nodes in each first-level grid unit after preliminary processing is calculated, and then the stress distribution level (A level, B level or C level) is determined according to the value of σ.

[0188] Finally, according to the corresponding precision grid unit division instruction, the division is carried out to generate the second-level grid unit. For example, a quadrilateral first-level grid unit after preliminary processing is calculated to have a stress distribution degree of B level, and then it is divided according to the second-level precision grid unit division instruction to obtain the final second-level grid unit.

[0189] Embodiment 9: The coal yard steel structure simulation model is optimized, comprising:

[0190] Based on the second-level grid unit, the initial material properties in the set first-level grid unit are re-corrected to generate a coal yard steel structure simulation model;

[0191] Based on the obtained coal yard historical environment data, the reliability of the coal yard steel structure simulation model is evaluated to generate a reliability evaluation result;

[0192] And based on the reliability evaluation result, the coal yard steel structure simulation model is optimized.

[0193] In this embodiment, the initial material property correction is based on

[0194] The initial material properties are corrected based on the stress conditions and geometric characteristics of the secondary grid elements. Different stress and geometric shapes can affect the stress distribution within the material, which in turn affects the actual performance of the material. For example, if a secondary grid element is in a high stress area (calculated based on node stress) and has a relatively long shape (such as a long and narrow four-sided element), it may cause the material in that area to be more prone to fatigue or deformation, requiring adjustments to the material's elastic modulus, yield strength, and other properties.

[0195] Generating a simulation model of the coal yard steel structure

[0196] Using the corrected material properties, combined with the overall geometric shape and boundary conditions of the coal yard steel structure, a simulation model of the coal yard steel structure is constructed. For example, in a finite element analysis software, input the corrected material properties, the beam and column sizes of the steel structure, the connection method, and the support conditions of the steel structure (such as fixed at the bottom, free at the top, etc.), to generate a simulation model of the coal yard steel structure for reliability assessment.

[0197] (I) Determining the optimization target

[0198] Determining the target based on the evaluation results

[0199] If the reliability evaluation result shows that the strength reliability is low, the optimization target is to improve the strength of the structure. For example, when the strength reliability Rs = 0.8, measures need to be taken to improve the carrying capacity of the structure, so that the strength reliability reaches 1.2 or above.

[0200] If the stiffness reliability is insufficient, the optimization target is to reduce the deformation of the structure, such as increasing the stiffness reliability from 1.5 to 2 or above.

[0201] If the stability reliability is low, the optimization target is to enhance the stability of the structure to ensure that the stability reliability is greater than 1.3.

[0202] (II) Optimization operation

[0203] Adjusting the geometric shape of the structure

[0204] If the strength or stability needs to be improved, the geometric shape of the steel structure can be adjusted. For example, increase the cross-sectional size of the beam and column. For columns, the diameter of the circular cross-section can be increased from d = 0.5 meters to d = 0.6 meters; for beams, the height of the I-beam can be increased from h = 0.3 meters to h = 0.35 meters. Or change the layout of the structure, such as adding support structures, and adding diagonal bracing in areas with large spans to improve the overall stability of the structure.

[0205] Improving material performance

[0206] If the intensity reliability is low, consider replacing the material with higher strength. For example, replace the original Q345 steel with Q420 steel, or perform surface strengthening treatment on the existing material, such as shot blasting to improve the hardness and strength of the material surface, thereby improving the overall strength of the structure.

[0207] Model re-evaluation

[0208] After adjusting the structure geometry or improving the material performance, rebuild the coal yard steel structure simulation model and perform reliability evaluation based on the historical environmental data of the coal yard again to check whether the optimized model meets the predetermined reliability target. If not, continue to adjust and optimize until a satisfactory reliability level is reached.

[0209] In the generation of the reliability evaluation result, the embodiment 10 includes:

[0210] Structural reliability evaluation and response reliability evaluation;

[0211] The structural reliability evaluation includes:

[0212] Perform static analysis on the coal yard steel structure simulation model to determine the stability value of the coal yard steel structure simulation model under static load;

[0213] Determine the natural frequency and modal shape of the secondary grid element through modal analysis to generate the resonance risk value of the coal yard steel structure simulation model;

[0214] Based on the stability value and resonance risk value of the coal yard steel structure simulation model, generate a structural reliability evaluation value;

[0215] The response reliability evaluation includes:

[0216] Perform seismic simulation on the coal yard steel structure simulation model based on historical data to calculate the residual safety reserve of the structure and generate a first response reliability evaluation value;

[0217] Based on historical data, simulate the settlement process of the foundation by gradually applying load, calculate the change of the stress field of the coal yard steel structure simulation model, and generate a second response reliability evaluation value;

[0218] Combine the first response reliability evaluation value and the second response reliability evaluation value to generate a response reliability evaluation value;

[0219] Weighted sum of the structural reliability evaluation value and the response reliability evaluation value to generate a reliability evaluation result.

[0220] In this embodiment, the static analysis determines the stability value

[0221] When performing static analysis on the simulation model of the coal yard steel structure, various static loads that the coal yard steel structure bears in the normal working state need to be considered, such as self-weight, weight of coal pile, weight of equipment, etc. These loads act on the corresponding nodes of the structure in the form of force.

[0222] According to the principle of structural mechanics, the internal force and deformation of the structure under these static loads are solved by establishing a balance equation. For example, for a simple beam structure, according to the force condition of the beam (such as concentrated force, uniformly distributed force, etc.), the bending moment, shear force and deflection of the beam are solved by using the bending moment-curvature relationship, force balance equation, etc.

[0223] Determination of stability value

[0224] After the internal force and deformation are solved, the stability value of the simulation model of the coal yard steel structure is determined by a specific stability criterion. For example, for a frame structure, the critical load coefficient can be used to measure the stability. Assuming that the critical load coefficient of the coal yard steel structure under static load is k = 1.5 (when k > 1, it indicates that the structure is stable under the current load) through static analysis calculation, this k value is a kind of representation of the stability value.

[0225] Foundation settlement simulation process

[0226] Based on historical data, the settlement process of the foundation is simulated by gradually applying load. According to the principle of soil mechanics, the settlement of the foundation is related to the properties of the foundation soil (such as compression modulus, Poisson's ratio, etc.), load distribution of the upper structure, etc.

[0227] In the simulation process, the load of the simulation model of the coal yard steel structure is gradually increased, while considering the deformation characteristics of the foundation soil, the settlement of the foundation is calculated. According to the change of the foundation settlement, the change of the stress field of the simulation model of the coal yard steel structure is calculated

[0228] Finally, it should be noted that those skilled in the art can make various modifications and variations to the present application without departing from the spirit and scope of the present application. Thus, if these modifications and variations of the present application fall within the scope of the claims and their equivalent technologies, the present application also intends to include these modifications and variations.

[0229] The above is only a preferred embodiment of the present application, and is not intended to limit the protection scope of the present application.

Claims

1. A coal yard steel structure simulation analysis method, characterized in that, The method comprises the following steps: collecting field data of a coal mine site to construct a geometric model of a steel structure of the coal mine site; dividing the geometric model of the steel structure into grids and setting material properties inside the grids to generate a simulation model of the steel structure of the coal mine site; performing reliability evaluation on the simulation model of the steel structure of the coal mine site based on historical data to generate a reliability evaluation result; optimizing the simulation model of the steel structure of the coal mine site based on the reliability evaluation result; when the geometric model of the steel structure is divided into grids, the method comprises the following steps: constructing the geometric model of the steel structure by defining the size of nodes, the size of members and the shape of members; performing preliminary stress analysis on the geometric model of the steel structure, and performing preliminary grid division on the geometric model of the steel structure according to the result of the preliminary stress analysis to generate a plurality of first-level grid units; setting initial material properties of the first-level grid units based on historical data; calculating the complexity score C of the first-level grid units; and executing a grid division strategy based on the complexity score C to generate second-level grid units; when the plurality of first-level grid units are generated, the method comprises the following steps: performing simulation tests on the geometric model of the steel structure based on the obtained historical environmental data of the steel structure of the coal mine site; setting monitoring nodes at the nodes of the geometric model of the steel structure and at the centers of the members, and applying simulation forces to the geometric model of the steel structure to generate a stress data set D1 at the monitoring nodes, D1 = [d1, d2…di…dn]; wherein di represents the stress size of the i-th monitoring node, and n represents the number of monitoring nodes; performing first-level grid unit division on the stress data set D1 at the monitoring nodes to generate a stress data set D2 at the remaining monitoring nodes; continuing to perform first-level grid unit division on the stress data set D2 at the remaining monitoring nodes until all monitoring nodes are divided into first-level grid units; generating a plurality of first-level grid units, including M1, M2…Mi…Mm; wherein Mi represents the i-th first-level grid unit, and m is the number of generated first-level grid units, m < n; the first-level grid unit division comprises: setting a first deviation value a; calculating the stress deviation value A of d1 and all adjacent monitoring nodes of d1; selecting all adjacent monitoring nodes of A < a, then dividing d1 and the selected all adjacent monitoring nodes into a first-level grid unit; and removing the plurality of monitoring nodes contained in the generated first-level grid unit from the stress data set D1 at the monitoring nodes; generating a stress data set D2 at the remaining monitoring nodes; when the stress deviation value A of d1 and all adjacent monitoring nodes of d1 is calculated, it comprises the following steps: calculating the comprehensive deviation value ΔFd in the direction: taking the 1st monitoring node d1 as the current monitoring node; decomposing the forces borne by the current monitoring node and the adjacent monitoring nodes into x, y and z direction components; respectively calculating the difference values of the x, y and z direction components to generate ΔFx, ΔFy and ΔFz; wherein ΔFx represents the difference value between the x direction component of the current monitoring node and the x direction component of the adjacent monitoring node; ΔFy represents the difference value between the y direction component of the current monitoring node and the y direction component of the adjacent monitoring node; and ΔFz represents the difference value between the z direction component of the current monitoring node and the z direction component of the adjacent monitoring node; Generate the comprehensive deviation value ΔFd in the direction by vector synthesis; Calculate the deviation value of the force size |d1-d(1+j)| Wherein, d(1+j) represents the force value of the adjacent monitoring node of the current monitoring node, j represents the adjacent distance of the current monitoring node; A=k1*ΔFd+k2*|d1-d(1+j)|,k1+k2=1; K1 is the weight of the comprehensive deviation value in the direction, and k2 is the weight of the deviation value of the force size.

2. The coal yard steel structure simulation analysis method according to claim 1, wherein, When calculating the complexity score C of the primary grid unit, the number q of monitoring nodes contained in the current primary grid unit is obtained, and the aspect ratio r of the current primary grid unit and the force distribution degree σ of the current primary grid unit are calculated to generate the primary grid complexity score C. C=w1×q+w2×r+w3×σ; Wherein, w1 is the weight of the number q of monitoring nodes contained in the current primary grid unit, w2 is the weight of the aspect ratio r of the current primary grid unit, and w3 is the weight of the force distribution degree σ of the current primary grid unit, and w1+w2+w3=1, the values of w1, w2 and w3 are adjusted according to historical data.

3. The coal yard steel structure simulation analysis method according to claim 2, wherein, When calculating the aspect ratio r in the primary grid unit and the force distribution degree σ in the primary grid, each internal angle α in the primary grid unit is obtained, α=(α1, α2, …, αn). Wherein, αn represents the nth internal angle. Calculate the standard deviation σa of the internal angle in the primary grid unit; Based on the standard deviation σa of the internal angle in the primary grid unit, a plurality of internal angle standard deviation thresholds are set, and the aspect ratio r of the current primary grid unit is generated based on the internal angle standard deviation threshold. Calculate and set a plurality of force standard deviation thresholds based on the standard deviation σf of the force; Based on the force standard deviation threshold, the force distribution degree σ of the current primary grid unit is generated.

4. The coal yard steel structure simulation analysis method according to claim 2, wherein, When the grid division strategy is executed based on the complexity score C to generate the secondary grid unit, the complexity threshold C1 is set based on the complexity score C. If C≤C1, the complexity level of the current primary grid unit is first-level complexity, and the primary grid division strategy is executed; If C1<c, the complexity level of the current primary grid unit is second-level complexity, and the advanced grid division strategy is executed.

5. The coal yard steel structure simulation analysis method according to claim 4, wherein, The grid division strategy includes: the primary grid division strategy includes: equal proportion division based on the shape of the current primary grid unit to generate the secondary grid unit; The advanced grid division strategy includes: calculating the force distribution degree σ of the nodes in the current primary grid unit; and dividing the current primary grid unit based on the force distribution degree σ inside the current primary grid unit to generate a plurality of secondary grid units.

6. The coal yard steel structure simulation analysis method according to claim 5, wherein, When the coal yard steel structure simulation model is optimized, the initial material properties in the secondary grid unit are corrected to generate the coal yard steel structure simulation model; Based on the obtained historical environmental data of the coal yard, the reliability of the coal yard steel structure simulation model is evaluated to generate a reliability evaluation result; And based on the reliability evaluation result, the coal yard steel structure simulation model is optimized.

7. The coal yard steel structure simulation analysis method according to claim 6, wherein When the reliability evaluation result is generated, structure reliability evaluation and response reliability evaluation are included. The structural reliability evaluation includes: Performing static analysis on the coal yard steel structure simulation model to determine the stability value of the coal yard steel structure simulation model under static load; Determine the natural frequency and modal shape of the secondary grid element through modal analysis to generate the resonance risk value of the coal yard steel structure simulation model; Based on the stability value and the resonance risk value of the coal yard steel structure simulation model, generate the structural reliability evaluation value; The response reliability evaluation includes: Performing seismic simulation on the coal yard steel structure simulation model based on historical data to calculate the residual safety reserve of the structure and generate the first response reliability evaluation value; Based on historical data, simulate the settlement process of the foundation by gradually applying load, calculate the change of the stress field of the coal yard steel structure simulation model, and generate the second response reliability evaluation value; Combine the first response reliability evaluation value and the second response reliability evaluation value to generate the response reliability evaluation value; Weighted sum of the structural reliability evaluation value and the response reliability evaluation value to generate the reliability evaluation result.

Images