Roadbed stability dynamic analysis method based on building information model and finite difference

Abstract

The invention relates to the field of building construction, in particular to a roadbed stability dynamic analysis method based on a building information model and finite difference, through dynamic updating of a BIM model and visual analysis result feedback, constructors can be guided to accurately control the compactness and layering thickness of filling materials, so that the construction quality is improved, and the construction period is shortened. And the quality problem caused by uneven filling or weak foundation is reduced. Through bidirectional association of BIM and FLAC3D, a closed-loop process of construction data, numerical analysis and result feedback is formed, and the time from monitoring to analysis is greatly shortened. In the construction process, the construction technology can be rapidly adjusted according to field data, reworking and delay are avoided, and the construction efficiency is improved. FLAC3D analysis results, such as a settlement curve, a safety coefficient distribution diagram and a stress distribution diagram, are visually displayed by adopting a three-dimensional cloud picture, so that constructors and managers can quickly understand the analysis results, the technical threshold is reduced, and meanwhile, a more scientific basis is provided for construction decision making.

Description

Technical Field

[0001] The present invention relates to the technical field of construction engineering, and particularly to a dynamic analysis method for subgrade stability based on building information modeling and finite difference method. Background Art

[0002] High fill subgrade is a common structure in transportation infrastructure construction (such as highways, railways, and large-scale site leveling projects), and its construction quality directly affects the service life and safety of the road. However, due to the large filling height, diverse material types, and complex foundation conditions during the construction of high fill subgrade, it often faces a series of technical problems, such as uneven settlement, sliding instability, insufficient compaction degree, etc. The occurrence of these problems is often closely related to the emergence of adverse working conditions during the construction process and the insufficient response to the dynamic changes in the construction process.

[0003] Currently, traditional subgrade stability analysis methods are mainly based on static models and manual calculations, which cannot reflect the dynamic changes during the construction process in real time and are difficult to cope with sudden adverse working conditions during construction (such as excessive settlement of the filling layer, foundation deformation, rainwater infiltration, etc.); at the same time, the existing building information modeling technology is mainly applied to construction design and progress management in subgrade construction, lacking the real-time interaction ability with the actual on-site situation. Summary of the Invention

[0004] In view of the above problems, an embodiment of the present invention provides a dynamic analysis method for subgrade stability based on building information modeling and finite difference method.

[0005] The dynamic analysis method for subgrade stability based on building information modeling and finite difference method provided by an embodiment of the present invention includes:

[0006] S1. Construct a BIM model of the high fill subgrade, where the BIM model includes the geometric attributes of the high fill subgrade and non-geometric attributes representing physical and mechanical parameters;

[0007] S2. Collect the geological change information at the construction site during the construction process, and input the geological change information into the BIM model. Through dynamic update, the real-time synchronization of the BIM model is realized, so that the BIM model can accurately reflect the changes in the construction state and geological conditions;

[0008] S3. Convert the real-time synchronized BIM model into a geometric structure file suitable for FLAC3D, and convert the geometric structure file into a FLAC3D numerical analysis model through geometric simplification and mesh generation, and define the material properties, boundary conditions, and loading conditions in the FLAC3D numerical analysis model;

[0009] S4. Conduct stability numerical analysis on the FLAC3D numerical analysis model to obtain analysis results. The content of the stability numerical analysis includes: settlement analysis, slip analysis, and stress distribution analysis. The analysis results include settlement data, safety factor distribution data, and stress distribution data.

[0010] S5. Feed the analysis results back to the BIM model after real-time synchronization, and visually display the settlement contour map, safety factor distribution contour map, and stress distribution contour map in the BIM model after real-time synchronization in the form of contour maps.

[0011] S6. Optimize the construction process according to the settlement contour map, safety factor distribution contour map, and stress distribution contour map displayed in the BIM model, and update the BIM model according to the optimized construction process.

[0012] Optionally,

[0013] The geometric properties include: the cross-sectional shape of the high fill subgrade, the thickness of the fill layer, the slope boundary, and the construction boundary conditions.

[0014] The non-geometric properties include: the physical and mechanical parameters of the fill material: density, elastic modulus, cohesion, internal friction angle, Poisson's ratio, and the physical and mechanical parameters of the foundation soil layer: shear strength, compression modulus.

[0015] Optionally, the process of constructing the BIM model of the high fill subgrade in S1 includes:

[0016] S11. Establish a layered three-dimensional model of the high fill subgrade, including the layered geometric information of different fill layers and the foundation soil.

[0017] S12. Define the physical and mechanical parameters of the fill material and the physical and mechanical parameters of the foundation soil layer.

[0018] S13. Update the geometric structure and physical properties of the layered three-dimensional model according to the geological boundary, soil layer distribution, and thickness change obtained from drilling and testing.

[0019] Optionally, the dynamic update in S2 includes:

[0020] Update of the foundation geometric shape: Dynamically adjust the geometric shape attributes of the BIM model foundation according to the change of the foundation contour during the construction excavation or filling process.

[0021] Addition of new soft soil layers: Input the thickness and distribution range of the soft soil layer found by on-site investigation into the BIM model, and adjust the shear strength and compression modulus parameters of the corresponding soil layer.

[0022] Change of the groundwater level: Monitor the dynamic change of the groundwater level during the construction period and use it as a parameter for updating the BIM model.

[0023] Optionally, the process of performing stability numerical analysis in S4 includes:

[0024] S41. Define the material properties of the high fill subgrade, including: the density, Poisson's ratio, internal friction angle, and cohesion of the fill material, and the elastic model, shear strength, and compression modulus of the foundation soil;

[0025] S42. Set the dynamic load changes during the construction of layered filling and simulate the interaction between the filling layer and the foundation;

[0026] S43. Apply the strength reduction method to analyze the stability of the high fill subgrade, and calculate the safety factor and the position of the most dangerous slip surface;

[0027] S44. Simulate the dynamic changes of the settlement of the filling layer during the construction process, and generate the settlement curve and settlement distribution map;

[0028] S45. Analyze the stress concentration areas of the foundation and the filling layer, and predict the potential instability risks.

[0029] Optionally, S5 specifically includes:

[0030] S51. Use a Python script to read the analysis results and generate a structured output file in CSV format;

[0031] S52. Use the Dynamo tool to read the structured output file and reconstruct the 3D Mesh grid based on the BIM model after real-time synchronization;

[0032] S53. Generate the following cloud maps in the reconstructed BIM model environment:

[0033] Settlement cloud map: Different color areas correspond to different settlement amounts, stress distribution cloud map: Intuitively display the stress distribution of the filling layer and the foundation, safety factor distribution cloud map: Identify the position of the most dangerous slip surface and the corresponding safety factor.

[0034] Optionally, the process of optimizing the construction technology in S6 includes:

[0035] Adjust the compaction technology of the fill material, including increasing the frequency of compaction equipment or adjusting the parameters of compaction equipment;

[0036] Improve the mix ratio of the fill material, including adjusting the proportion of sand or stone;

[0037] Increase the foundation reinforcement measures in local areas, including using replacement soil, reinforcement piles or geogrids;

[0038] Modify the construction sequence or the thickness of layered filling to reduce the slip risk during the process of layered filling;

[0039] Optimize the construction drainage measures to reduce the impact of rain erosion on the stability of the fill subgrade.

[0040] Optionally, the dynamic analysis method for subgrade stability is applied to the construction of high-fill earthwork subgrades under the following complex working conditions:

[0041] Areas where the distribution of foundation soil layers is uneven, there are soft interlayers, or the groundwater level is relatively high;

[0042] Construction environments that require real-time monitoring of geological changes and rapid adjustment of construction techniques;

[0043] Construction environments with extreme working conditions such as heavy rain, seismic loads, or sudden changes in construction loads;

[0044] Highways, railways, airports, and other transportation infrastructure projects with high requirements for foundation stability.

[0045] Compared with the prior art, the beneficial effects of the present invention are as follows: Through the dynamic update of the BIM model and the feedback of the visualized analysis results, it is possible to guide construction personnel to accurately control the compaction degree and layer thickness of the filling materials, thereby improving the construction quality and reducing quality problems caused by uneven filling or weak foundations. Through the bidirectional association between BIM and FLAC3D, a closed-loop process of construction data, numerical analysis, and result feedback is formed, significantly shortening the time from monitoring to analysis. During the construction process, the construction techniques can be quickly adjusted according to the on-site data, avoiding rework and delays, and improving the construction efficiency.

[0046] The analysis results of FLAC3D (such as settlement curves, safety factor distribution maps, and stress distribution maps) are intuitively displayed using three-dimensional cloud maps, enabling construction personnel and management personnel to quickly understand the analysis results, reducing the technical threshold, and providing a more scientific basis for construction decisions.

[0047] The present invention is particularly suitable for construction scenarios with a large filling height, complex foundation conditions (such as sites with soft soil layers or significant hydrogeological effects), and the need for real-time monitoring and dynamic response, and has significant application effects in large-scale projects such as highways, railways, and airports. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] The drawings described herein are used to provide a further understanding of the present invention, form a part of this application, and do not limit the present invention. In the drawings:

[0049] Figure 1 is a flowchart of a dynamic analysis method for subgrade stability based on building information model and finite difference provided by an embodiment of the present invention;

[0050] Figure 2 is a principle block diagram of a dynamic analysis method for subgrade stability based on building information model and finite difference provided by an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0051] To make the objectives, technical solutions, and advantages of the present invention more clear and understandable, the present invention will be further described in detail below in conjunction with the embodiments and the accompanying drawings. Herein, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, but do not limit the present invention.

[0052] See Figure 1 and Figure 2 , a method for dynamic analysis of subgrade stability based on building information modeling and finite difference provided by an embodiment of the present invention includes:

[0053] S1. Construct a BIM model of a high embankment subgrade, where the BIM model includes the geometric attributes of the high embankment subgrade and non-geometric attributes representing physical and mechanical parameters;

[0054] S2. Collect the geological change information at the construction site during the construction process, and use BIM secondary development technology to input the geological change information into the BIM model. Through dynamic update, real-time synchronization of the BIM model is achieved, so that the BIM model can accurately reflect the changes in the construction state and geological conditions;

[0055] S3. Convert the real-time synchronized BIM model into a geometric structure file suitable for FLAC3D. Through geometric simplification and mesh generation for the geometric structure file, a FLAC3D numerical analysis model is obtained, and the material properties, boundary conditions, and loading conditions in the FLAC3D numerical analysis model are defined;

[0056] S4. Conduct stability numerical analysis on the FLAC3D numerical analysis model to obtain analysis results. The content of the stability numerical analysis includes: settlement analysis for calculating the settlement amount distribution and settlement curve at each stage during the construction of the embankment subgrade; slip analysis using the strength reduction method to calculate the safety factor and the position of the most dangerous slip surface of the high embankment subgrade, and stress distribution analysis to simulate the stress distribution in the embankment materials and the foundation, and predict the instability area; the analysis results include settlement data, safety factor distribution data, and stress distribution data;

[0057] S5. Feed back the analysis results to the real-time synchronized BIM model, and visually display the settlement contour map, safety factor distribution contour map, and stress distribution contour map in the real-time synchronized BIM model in the form of contour maps;

[0058] S6. Optimize the construction process according to the settlement contour map, safety factor distribution contour map, and stress distribution contour map shown in the BIM model, and update the BIM model according to the optimized construction process.

[0059] After updating the BIM model, repeat the above analysis process to form a dynamic response closed-loop process for construction stability analysis.

[0060] In implementation, the geometric properties include: the cross-section shape of the high embankment subgrade, the thickness of the filling layer, the slope boundary, and the construction boundary conditions;

[0061] The non-geometric properties include: the physical and mechanical parameters of the filling material: density, elastic modulus, cohesion, internal friction angle, Poisson's ratio, and the physical and mechanical parameters of the foundation soil layer: shear strength, compression modulus.

[0062] In implementation, the process of constructing the BIM model of the high embankment subgrade in S1 includes:

[0063] S11. A layered 3D model of the high embankment subgrade can be established through Revit, including the layered geometric information of different filling layers and the foundation soil;

[0064] S12. Define the physical and mechanical parameters of the filling material (such as density, elastic modulus, cohesion, internal friction angle, Poisson's ratio, etc.) and the physical and mechanical parameters of the foundation soil layer (such as shear strength, compression modulus, etc.);

[0065] S13. Update the geometric structure and physical properties of the layered 3D model according to the geological boundary, soil layer distribution, and thickness changes obtained from drilling and testing.

[0066] In implementation, the dynamic update in S2 includes:

[0067] Update of the foundation geometric shape: Dynamically adjust the geometric shape attributes of the BIM model foundation according to the changes in the foundation contour during construction excavation or filling;

[0068] Addition of new soft soil layers: Input the thickness and distribution range of the soft soil layers found through on-site investigation into the BIM model, and adjust the shear strength and compression modulus parameters of the corresponding soil layers;

[0069] Change in the groundwater level: Monitor the dynamic changes in the groundwater level during the construction period and use it as a parameter for BIM model update.

[0070] In implementation, the process of performing stability numerical analysis in S4 includes:

[0071] S41. Define the material properties of the high embankment subgrade, including: the density, Poisson's ratio, internal friction angle, and cohesion of the filling material, and the elastic model, shear strength, and compression modulus of the foundation soil;

[0072] S42. Set the dynamic load changes during construction layered filling and simulate the interaction between the filling layer and the foundation;

[0073] S43. Apply the strength reduction method to analyze the stability of the high embankment subgrade, and calculate the safety factor and the location of the most dangerous slip surface;

[0074] S44. Simulate the dynamic changes in the settlement of the filled layer during the construction process, and generate a settlement curve and a settlement distribution map;

[0075] S45. Analyze the stress concentration areas of the foundation and the filled layer, and predict potential instability risks.

[0076] During implementation, S5 specifically includes:

[0077] S51. Use a Python script to read the analysis results and generate a structured output file in CSV format;

[0078] S52. Use the Dynamo tool to read the structured output file and reconstruct a three-dimensional Mesh grid based on the BIM model after real-time synchronization;

[0079] S53. Generate the following cloud maps in the reconstructed BIM model environment:

[0080] Settlement cloud map: Different color areas correspond to different settlement amounts, stress distribution cloud map: Intuitively display the stress distribution of the filled layer and the foundation, safety factor distribution cloud map: Identify the location of the most dangerous slip surface and the corresponding safety factor.

[0081] During implementation, the process of optimizing the construction technology in S6 includes:

[0082] Adjust the compaction technology of the filled material, including increasing the frequency of the compaction equipment or adjusting the parameters of the compaction equipment;

[0083] Improve the mixture ratio of the filled material, including adjusting the proportion of sand or stone;

[0084] Add foundation reinforcement measures in local areas, including using replacement soil, reinforcement piles or geogrids;

[0085] Modify the construction sequence or the thickness of layered filling to reduce the slip risk during the layered filling process;

[0086] Optimize the construction drainage measures to reduce the impact of rain erosion on the stability of the filled subgrade.

[0087] During implementation, the dynamic analysis method for subgrade stability is applied to the construction of high-filled earth subgrades under the following complex working conditions:

[0088] Areas with uneven distribution of foundation soil layers, presence of soft interlayers or relatively high groundwater levels;

[0089] Construction environments that require real-time monitoring of geological changes and rapid adjustment of construction technology;

[0090] Construction environments with extreme working conditions such as heavy rain, seismic loads or sudden changes in construction loads;

[0091] Highway, railway, airport and other transportation infrastructure projects with high requirements for foundation stability.

[0092] During implementation, the dynamic update of the BIM model and the numerical analysis of FLAC3D are achieved through two-way interaction: the BIM model obtains construction data in real time and updates model attributes, which are used as input parameters for FLAC3D analysis; the analysis results of FLAC3D are fed back into the BIM model through Mesh reconstruction technology to visually present the analysis results for construction personnel to reference; the BIM model parameters are adjusted according to the analysis feedback of FLAC3D and numerical analysis is carried out again to ensure the reliability after the construction process is optimized.

[0093] The solution provided by the embodiment of the present invention can guide construction personnel to accurately control the compaction degree and layer thickness of filling materials through the dynamic update of the BIM model and the feedback of visual analysis results, thereby improving the construction quality and reducing quality problems caused by uneven filling or weak foundation. Through the two-way association between BIM and FLAC3D, a closed-loop process of construction data, numerical analysis and result feedback is formed, greatly shortening the time from monitoring to analysis. During the construction process, the construction process can be quickly adjusted according to on-site data, avoiding rework and delays and improving construction efficiency.

[0094] The analysis results of FLAC3D (such as settlement curves, safety factor distribution maps and stress distribution maps) are intuitively displayed using three-dimensional cloud maps, enabling construction personnel and management personnel to quickly understand the analysis results, lowering the technical threshold, and providing a more scientific basis for construction decisions.

[0095] The above are only the preferred embodiments of the present invention and are not intended to limit the protection scope of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention are included in the protection scope of the present invention.

Claims

1. A dynamic analysis method for roadbed stability based on building information model and finite difference, characterized in that: include: S1. Construct a BIM model of the high fill roadbed, which includes geometric properties of the high fill roadbed and non-geometric properties representing physical and mechanical parameters; S2. Collect geological change information at the construction site during the construction process, and input the geological change information into the BIM model, and achieve real-time synchronization of the BIM model through dynamic update, so that the BIM model can accurately reflect the changes in construction status and geological conditions; S3, converting the real-time synchronized BIM model into a geometric structure file suitable for FLAC3D, converting the geometric structure file into a FLAC3D numerical analysis model through geometric simplification and meshing, and defining the material properties, boundary conditions and loading conditions in the FLAC3D numerical analysis model; S4. Perform stability numerical analysis on the FLAC3D numerical analysis model to obtain analysis results. The stability numerical analysis includes: settlement analysis, slip analysis, and stress distribution analysis. The analysis results include settlement data, safety factor distribution data, and stress distribution data. S5. Feedback the analysis results to the real-time synchronized BIM model, and visualize the settlement cloud map, safety factor distribution cloud map and stress distribution cloud map in the real-time synchronized BIM model in the form of cloud maps; S6. Optimize the construction process according to the settlement cloud map, safety factor distribution cloud map and stress distribution cloud map displayed in the BIM model, and update the BIM model according to the optimized construction process.

2. The method for dynamic analysis of roadbed stability based on building information model and finite difference according to claim 1, characterized in that: Geometric attributes include: high fill roadbed section shape, fill layer thickness, slope boundary, and construction boundary conditions; Non-geometric properties include: physical and mechanical parameters of fill materials: density, elastic modulus, cohesion, internal friction angle, Poisson's ratio; physical and mechanical parameters of foundation soil layers: shear strength and compression modulus.

3. The method for dynamic analysis of roadbed stability based on building information model and finite difference as claimed in claim 2, characterized in that: The process of constructing the BIM model of the high fill roadbed in S1 includes: S11, establishing a layered three-dimensional model of the high fill roadbed, including layered geometric information of different fill layers and foundation soil; S12. Define the physical and mechanical parameters of the filling material and the physical and mechanical parameters of the foundation soil layer; S13. Update the geometric structure and physical properties of the layered 3D model based on the geological boundaries, soil layer distribution and thickness changes obtained through drilling and testing.

4. The method for dynamic analysis of roadbed stability based on building information model and finite difference according to claim 1, characterized in that: Dynamic updates in S2 include: Foundation geometry update: Dynamically adjust the geometric attributes of the BIM model foundation according to the changes in the foundation contour during construction excavation or filling; Add soft soil layer: input the thickness and distribution range of the soft soil layer found in the on-site survey into the BIM model, and adjust the shear strength and compression modulus parameters of the corresponding soil layer; Groundwater level changes: Monitor the dynamic changes of groundwater levels during the construction period as a parameter for updating the BIM model.

5. The method for dynamic analysis of roadbed stability based on building information model and finite difference according to claim 1, characterized in that: The process of numerical stability analysis in S4 includes: S41. Define the material properties of high fill roadbed, including density, Poisson's ratio, internal friction angle, cohesion of fill material, elastic model, shear strength, and compression modulus of foundation soil; S42, set the dynamic load change of the construction layer filling, and simulate the interaction between the filling layer and the foundation; S43. Use the strength reduction method to analyze the stability of high fill roadbed and calculate the safety factor and the location of the most dangerous sliding surface; S44, simulating the dynamic changes of the settlement of the filling layer during the construction process, and generating a settlement curve and a settlement distribution diagram; S45. Analyze the stress concentration areas of the foundation and fill layer and predict potential instability risks.

6. The method for dynamic analysis of roadbed stability based on building information model and finite difference according to claim 5, characterized in that: S5 specifically includes: S51, using a Python script to read the analysis results and generate a structured output file in CSV format; S52. Use Dynamo to read the structured output file and reconstruct the 3D Mesh based on the real-time synchronized BIM model; S53. Generate the following cloud map in the reconstructed BIM model environment: Settlement cloud map: different color areas correspond to different settlement amounts, stress distribution cloud map: intuitively displays the stress distribution of the fill layer and the foundation, safety factor distribution cloud map: identifies the location of the most dangerous sliding surface and the corresponding safety factor.

7. The method for dynamic analysis of roadbed stability based on building information model and finite difference according to claim 1, characterized in that: The process of optimizing the construction process in S6 includes: Adjust the compaction process of the fill material, including increasing the frequency of compaction equipment or adjusting compaction equipment parameters; Improve the proportion of fill materials, including adjusting the proportion of sand or stone; Adding foundation reinforcement measures in local areas, including using replacement fill, reinforcement piles or geogrids; Modify the construction sequence or layered filling thickness to reduce the risk of slippage during layered filling; Optimize construction drainage measures to reduce the impact of rainwater erosion on the stability of fill roadbed.

8. The method for dynamic analysis of roadbed stability based on building information model and finite difference according to any one of claims 1 to 7, characterized in that: The dynamic analysis method of roadbed stability is applied to the construction of high fill roadbed under the following complex working conditions: Areas with uneven foundation soil distribution, weak interlayers or high groundwater levels; Construction environments that require real-time monitoring of geological changes and rapid adjustment of construction techniques; Construction environments with extreme conditions such as heavy rain, earthquake loads or sudden changes in construction loads; Highways, railways, airports and other transportation infrastructure projects that require high foundation stability.

Images