Intelligent modeling and analysis method for shield tunnel distribution tenon type segment lining structure

Abstract

The invention belongs to the technical field of modeling and analysis of a shield tunnel segment lining structure, particularly relates to an intelligent modeling and analysis method for a shield tunnel distributed tenon type segment lining structure, and aims to solve the problems of high manual dependence, difficulty in batch operation, low result processing efficiency and the like in the modeling process of the shield tunnel segment structure. A whole set of shield segment modeling and simulation full-process automation method which takes a Python script as a core and is based on an Abaqus secondary development platform is provided. According to the method, a traditional modeling process of gradual interaction of a graphical user interface depending on Abaqus is converted into an automatic modeling mode driven by a full-process script, and a complete closed loop of automatic input of structural parameters, automatic construction of a geometric structure, automatic execution of Boolean cutting and automatic deployment of grid division is achieved. Meanwhile, the technology of the patent supports parameterized modification, can easily adapt to various structures with different radiuses, concave-convex tenon arrangements and bolt angles, and greatly enhances the reusability and engineering universality of the model.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of modeling and analysis of shield tunnel segment lining structure, and particularly relates to an intelligent modeling and analysis method for a shield tunnel distributed dovetail segment lining structure. BACKGROUND

[0002] With the development of shield technology, the demand for large shield tunnel engineering design and construction is increasing. In large-diameter shield tunnels, the distributed dovetail is one of the most commonly used inter-ring connection structures. Therefore, the demand for the design of inter-ring segment lining structures with distributed dovetail is increasing.

[0003] In the prior art, a user first manually and interactively establishes a three-dimensional entity model in the graphical user interface of Abaqus software through the menu bar, toolbar and model tree, including key steps such as geometry creation, material property definition, mesh division, and application of loads and boundary conditions. Then, the user needs to manually select the analysis step and solver options in the software interface and submit the calculation task. After the calculation is successfully completed, the user still needs to manually operate, for example, selecting the corresponding output variables such as displacement and internal force in the result visualization module, manually extracting internal force and deformation by creating a chart report or extracting numerical values on a specific path.

[0004] The prior art needs to use the graphical user interface of Abaqus for manual modeling, and its core defect is that the entire process highly depends on repetitive manual operation and lacks automation, resulting in extremely low modeling efficiency. Any minor modification of model size, load or boundary condition requires the user to spend a lot of time manually adjusting each item in the interface, especially when conducting series modeling or parameterized research. This manual reconstruction method takes a long time. At the same time, the complex interaction process is prone to human error, making it difficult to ensure the consistency and accuracy of model definition such as geometric relationship, material property, mesh division rule, and load boundary application position, thereby affecting the reliability of the simulation results. In the result extraction stage, the user must manually enter the post-processing module to obtain internal forces such as cross-section axial force, shear force and bending moment, and deformations such as key node displacement by repeatedly selecting, querying and positioning. This process is not only time-consuming and labor-intensive, but also has a very low degree of automation. Moreover, it is prone to errors and inefficient when extracting complex data such as path integral results or processing a large number of results. The existing manual mode has an insurmountable bottleneck in processing speed and scale, which greatly limits the deep application value of simulation analysis.

[0005] The core defect of the prior art lies in the lack of parameterized modeling and process automation capability at the bottom. Specifically, the modeling process completely relies on the interactive operation of the graphical user interface, such as menu clicking, object picking, parameter input, such as geometric size, material parameter, load value, boundary position, mesh seed, analysis step setting, and other model elements are isolated from each other and cannot be globally associated and dynamically updated through explicitly defined parameter variables. This makes any design change require the user to re-execute lengthy and error-prone operations to rebuild the model, rather than modifying the key parameters to achieve model modeling. This leads to a series of problems such as serious lack of adaptability to complex engineering requirements such as multi-parameter optimization.

[0006] In summary, in the design process of the distributed dowel type segment lining structure of the shield tunnel, a three-dimensional refined numerical model calculation analysis is performed to obtain the internal force and deformation of the structure. The traditional method is manual modeling and post-processing data extraction, which is time-consuming, laborious, low in efficiency and prone to errors, resulting in inaccurate or even incorrect structure analysis results. The present application provides an intelligent modeling and analysis method for shield tunnel distributed dowel type segment lining structure, which is a parameterized method for establishing distributed dowel type segment lining structure and automatically extracting structural stress and deformation data. Only a small amount of design parameters need to be input to intelligently obtain the calculation and analysis results, which greatly improves the efficiency of modeling and analysis compared with the traditional method. SUMMARY

[0007] The purpose of the present application is to provide an intelligent modeling and analysis method for shield tunnel distributed dowel type segment lining structure, which can solve the problem of "high dependence on Abaqus graphical user interface for modeling process and post-processing extraction process, repetitive and tedious modeling process and prone to errors".

[0008] The technical scheme adopted by the present application is as follows:

[0009] The intelligent modeling and analysis method for shield tunnel distributed dowel type segment lining structure is based on the automatic modeling method of Abaqus secondary development. The whole modeling process and post-processing extraction process are programmed and controlled through Python script to realize parameterization and automatic modeling of tunnel segment structure, and internal force and deformation data are automatically extracted after calculation to form a complete reusable technical process; including the following steps:

[0010] Step 1: parameter definition and rationality verification: define the structure parameters through the script input, and set the parameter legality judgment mechanism. If the structure closure or arrangement rule is not met, it will be directly exited. The structure parameters include segment outer diameter, block number, dowel slot position, bolt parameters, etc.

[0011] Step 2: Automatic generation of module structure: Generate standard blocks, abutting blocks, and capping blocks three-dimensional entity models according to parameters, and combine to form a ring structure, then automatically copy and cut to form front and rear ring structures;

[0012] Step 3: Arrangement of mortise and tenon and bolt structure: Automatically build mortise body, tenon body and bolt structure, generate front, middle and rear ring mortise and tenon combination and anchoring structure through geometric boolean and rotation arrangement operation, and complete automatic combination of connecting components;

[0013] Step 4: Generation and arrangement of loading components: Generate loading beams and place them on the top and side of the structure to facilitate subsequent load application;

[0014] Step 5: Structure grid division: Automatically generate grids according to preset block grid seed density parameters, and realize automatic grid division; Block types include standard blocks, abutting blocks, capping blocks, and mortise and tenon;

[0015] Step 6: Boundary condition and load definition: Automatically set material properties, such as interaction of contact pair definition, constraints and boundary conditions; Boundary conditions such as segment connection and fixed boundary, and define radial pressure load and tension load;

[0016] Step 7: Output complete model structure: Obtain a complete three-dimensional entity model for subsequent simulation analysis;

[0017] Step 8: Submit tasks and perform batch calculation: Generate multiple different working condition combinations, submit tasks for calculation, and obtain multiple ODB result files;

[0018] Step 9: Read ODB file: Automatically traverse all ODB result files in the specified directory and read the finite element simulation information therein;

[0019] Step 10: Read joint opening amount: Automatically identify joint node pairs according to model structure, obtain U1 and U2 displacement components of these nodes, calculate the Euclidean distance change between each pair of nodes, and output as joint opening amount data;

[0020] Step 11: Read structure deformation amount: Automatically select feature nodes, extract their displacement information, and output as structure deformation amount data combined with original radius data;

[0021] Step 12: Generate and add path: Automatically generate a circumferential path in the radial direction at the key section of the structure to facilitate subsequent extraction of internal force data on the free body cut surface;

[0022] Step 13: Extract internal force: Extract internal force data report on the path, and parse node position, time, and component force value in the report through regular expression, then extract bending moment, axial force and shear force data;

[0023] Step 14: result output and archiving: all extracted data is unified in CSV file format and named by structure name and variable type, and automatically classified and archived in a folder named by the ODB output file name.

[0024] The technical effects obtained by the present application are:

[0025] The patent proposes a complete set of automatic modeling and simulation method for shield segment based on Python script and Abaqus secondary development platform, aiming at the problems of high artificial dependence, batch operation difficulty and low result processing efficiency in the process of shield tunnel segment structure modeling. The patent changes the traditional modeling process which depends on Abaqus graphical user interface and gradually interacts into automatic modeling mode driven by the whole process script, realizes the complete closed loop of automatic input of structure parameters, automatic construction of geometric structure, automatic execution of Boolean cutting and automatic deployment of meshing. The test data shows that the traditional operation modeling of a single three-ring structure takes more than 10 hours, while the complete model generated by the patent method only takes 3-5 minutes, and the efficiency is greatly improved. At the same time, the patent technology supports parameterized modification, which can easily adapt to various structures with different radii, concave-convex tenon arrangement and bolt angles, greatly enhancing the reusability and engineering generality of the model. BRIEF DESCRIPTION OF DRAWINGS

[0026] Figure 1 is the intelligent modeling and analysis method for shield tunnel distributed tenon type segment lining structure of the present application;

[0027] Figure 2 is the overall model structure diagram of the distributed tenon type segment lining generated by the present application;

[0028] Figure 3 is the schematic diagram of automatic selection of joint opening measurement points of the present application;

[0029] Figure 4 is the standard block segment structure diagram when the number of concave-convex tenon of single block segment is 2;

[0030] Figure 5 is the standard block segment structure diagram when the number of concave-convex tenon of single block segment is 3. DETAILED DESCRIPTION

[0031] In order to make the purpose and advantages of the present application clearer and more apparent, the present application will be specifically described below in combination with examples. It should be understood that the following text is only used to describe one or several specific embodiments of the present application, and does not strictly limit the specific protection scope requested by the present application.

[0032] As Figures 1-5The shield tunnel distribution mortise type segment lining structure intelligent modeling and analysis method is shown, an automatic modeling method based on Abaqus secondary development, process programming control of the whole modeling process and post-processing extraction process is realized through a Python script, parametric and automatic modeling of the tunnel segment structure is realized, and internal force and deformation data are automatically extracted after calculation to form a complete reusable technical process; including the following steps:

[0033] Step 1: parameter definition and rationality check: define the structure parameters through the script input, and set the parameter legality judgment mechanism, if the structure closure or arrangement rule is not satisfied, it will be directly exited; the structure parameters are, for example, segment outer diameter, block number, mortise position, bolt parameters, etc.;

[0034] Step 2: automatic generation of module structure: generate standard block, adjacent block and capping block three-dimensional entity model according to the parameters, and combine to form the middle ring structure, then automatically copy and cut to form the front ring and rear ring structure;

[0035] Step 3: concave-convex mortise and bolt structure arrangement: automatically build concave groove body, convex mortise body and bolt structure, generate front, middle and rear ring concave-convex mortise combination and anchoring structure through geometric boolean and rotation arrangement operation, and complete the automatic combination of connecting components;

[0036] Step 4: loading component generation and arrangement: generate loading beam and place it on the top and side of the structure, which is convenient for subsequent load application;

[0037] Step 5: structure mesh division: automatically generate mesh according to the preset block mesh seed density parameter, realize automatic mesh division; the block types include standard block, adjacent block, capping block and mortise;

[0038] Step 6: boundary condition and load definition: automatically set material properties, such as interaction of contact pair definition, constraint and boundary condition; boundary conditions such as segment connection and fixed boundary, and define radial pressure load and tension load;

[0039] Step 7: output complete model structure: get the complete three-dimensional entity model, which can be directly used for subsequent simulation analysis;

[0040] Step 8: submit task and perform batch calculation: generate multiple different working condition combinations, submit tasks for calculation, and get multiple ODB result files;

[0041] Step 9: read ODB file: automatically traverse all ODB result files in the specified directory, and read the finite element simulation information in them;

[0042] Step 10: Read the joint opening amount: automatically identify the joint node pairs according to the model structure, and obtain the U1 and U2 displacement components of the nodes. Calculate the Euclidean distance change between each pair of nodes, and output as joint opening amount data;

[0043] Step 11: Read the structure deformation amount: automatically select feature nodes and extract their displacement information. Combine with the original radius data to output as structure deformation amount data;

[0044] Step 12: Generate and add paths: automatically generate circumferential paths in the radial direction of the key sections of the structure to facilitate the subsequent extraction of internal force data on the free body cut surface;

[0045] Step 13: Extract internal forces: extract internal force data reports on the path, and parse node positions, time, and component force values in the report through regular expressions, then extract bending moment, axial force, and shear force data;

[0046] Step 14: Result output and archiving: all extracted data is output in CSV file format, named according to structure name and variable type, and automatically archived in a folder named after the ODB output file.

[0047] Preferably, in the prior art, the angle division of the structural unit depends on manual specification, which is not accurate and efficient. In the existing modeling process of shield tunnel segments, the angle division of the capping block, the standard block and the adjacent block usually relies on experience formula or drawing assistance, and is completed by manual calculation and specification. In the modeling scene with multiple parameters and high precision requirements, manual specification of angles is prone to errors and cannot meet the automation needs of rapid iteration and batch modeling;

[0048] Therefore, the present application automatically deduces the angle of the structural unit based on the number of concave-convex tenons. The present application provides a parameterized modeling scheme for automatically deriving the angle of the structure, which automatically calculates the circumferential angle of the capping block, the standard block and the adjacent block by inputting the number of concave-convex tenons on each type of segment block, thereby realizing intelligent angle configuration of the structural blocks. Specifically, the user inputs the number of concave-convex tenons on each type of block (the number of concave-convex tenons on the standard block, the adjacent block and the capping block) and the total number of standard blocks in the ring, and the angles of each type of block can be automatically calculated. The number of concave-convex tenons includes the number of concave-convex tenons on the standard block, the adjacent block and the capping block. The formula is as follows:

[0049]

[0050]

[0051] ;

[0052] wherein, represents the angle corresponding to the standard block, representing the angle corresponding to the capping block, representing the angle corresponding to the adjacent block, representing the number of concave-convex tenons of the standard block, representing the number of concave-convex tenons of the capping block, representing the number of concave-convex tenons of the adjacent block, representing the total number of standard blocks;

[0053] automatically calculated , and will be used as control parameters in the geometric modeling function for geometric sketching and part arrangement, realizing the automatic layout of pipe piece blocking. As shown in Figure 2 .

[0054] Preferably, a structural rationality checking mechanism is introduced in the automatic modeling system of the pipe piece, which judges whether the closure condition is met by summing up the structural circumferential angle, and further judges whether the current blocking and concave-convex tenon arrangement scheme can constitute a complete structure; at the same time, it is verified whether the preset angle distribution of all concave-convex tenons conflicts with the key joint position, so that the concave-convex tenon elements reasonably avoid in the ring body; in the parameter setting stage, a variable automatic expansion method is adopted, according to the key parameters provided by the user, the key parameters including the number of standard blocks, the number of concave-convex tenons, and the width of the block, the total number of concave-convex tenons, the angle occupied by each block, the interval angle of the concave-convex tenons and other core variables are automatically derived, and multiple logical judgments and constraint rules are embedded, to ensure that the geometric logic is self-consistent and the mechanical parameters are effective, greatly reducing the modeling failure or simulation distortion phenomenon caused by input errors, and providing a key guarantee for the stable operation of automatic model generation and structural analysis.

[0055] Preferably, in the prior art, the joint surface node information is complex and difficult to manually specify, and it is difficult to realize automatic extraction and simulation. In the simulation analysis of shield tunnel structure, the opening amount (cracking amount), relative deformation and normal displacement of the joint surface are the key indicators for judging the structural integrity and joint performance. However, in the current mainstream modeling process, these joint surface nodes need to be selected or manually specified one by one in the graphical user interface of Abaqus, which is tedious and prone to error, and is difficult to popularize in complex or large-scale models. At the same time, due to the ring-like distribution of the pipe piece structure, the joint angles and coordinates of each block are different and do not have a unified rule, so it is difficult to establish an automatic extraction template by traditional methods.

[0056] Therefore, the present application proposes an automatic joint node extraction and displacement analysis method based on spatial distribution logic and angle symmetry matching, which can automatically filter out all effective joint nodes by inputting a small number of control parameters and generate the corresponding opening amount change data file; the control parameters are, for example, ring layer, radius, and concave-convex tenon arrangement; first, the joint node filtering conditions need to be defined, including:

[0057] (A) The coordinates satisfy the tolerance judgment at a certain Z level;

[0058] (B) The coordinate radius is close to the inner / outer ring radius; a very small deviation is allowed;

[0059] (C) The ring angle satisfies the closed interval of the actual concave-convex tenon structure arrangement;

[0060] Only the nodes that satisfy the above conditions will be used as candidate joint nodes. In addition, the present patent also adopts the polar angle sorting and difference matching method to sort the nodes in all different blocks according to the polar angle from small to large, and automatically identifies pairs of nodes with close angles by setting a tolerance of ±1°, to ensure that each pair of joint surface nodes has a corresponding relationship in the geometric sense. This algorithm can be used for both the butt joint between different blocks such as standard blocks and adjacent blocks, and the bolt position butt joint between different ring layers of the same block such as the front ring and the middle ring, as shown in Figure 3 .

[0061] In step 1, parameter definition and reasonableness verification: define the structure parameters through scripts; structure parameters such as pipe piece outer diameter, block number, tenon and groove position, and bolt parameters, and set a parameter legality judgment mechanism. If the structure closure or arrangement rule is not satisfied, the modeling process will be directly exited.

[0062] When detecting the structure closure, the angle sum is calculated by the formula as shown below :

[0063]

[0064] wherein, represents the angle corresponding to the standard block, represents the angle corresponding to the capping block, represents the total number of standard blocks. If , the current structure is considered invalid, and the modeling process will be automatically exited.

[0065] When detecting the arrangement rule of concave-convex tenon, the formula as shown below will be used for calculation.

[0066] For the th concave-convex tenon ( ), its angle position can be calculated by the formula:

[0067]

[0068] wherein, is the total number of concave-convex tenons, represents the reference angle of the first concave-convex tenon, represents the distribution interval of the concave-convex tenon.

[0069] Set of joint angles ; wherein is the total number of joints, in order to avoid the overlap between the position of the tongue-and-groove and the joint, it is required that for each angle position of the tongue-and-groove and each joint angle satisfies:

[0070]

[0071] for all tongue-and-groove and all joints , the above inequality must be established; if there is a pair ( ) such that the inequality does not hold, the current structure is considered invalid and the modeling process will automatically exit.

[0072] In step 5, structure meshing: through the preset block mesh seed density parameter, the mesh is automatically generated according to the block type, realizing automatic meshing;

[0073] When dividing the mesh, according to the judgment of different line segment types, the mesh is divided according to the preset block mesh seed density; for the standard block, the following piecewise constant function is used to distribute the seed number:

[0074]

[0075] wherein, represents a very small constant, , represents the length of the line segment, represents the length of the upper and lower long edges of the tongue-and-groove, represents the radius of the arc edge of the tongue-and-groove, represents the thickness of the tongue-and-groove, represents the radius shrinkage amount of the top surface of the tongue-and-groove, represents the width of the concrete segment, represents the number of mesh seeds corresponding to the upper and lower long edges of the tongue-and-groove, represents the number of mesh seeds corresponding to the arc edge of the tongue-and-groove, represents the number of mesh seeds corresponding to the inclined edge of the tongue-and-groove, represents the number of mesh seeds corresponding to the edge of the concrete segment parallel to the tunneling direction.

[0076] For the edges that cannot satisfy the above conditions, the preset distance set is also searched, if there is satisfies:

[0077]

[0078] the number of mesh seeds is distributed ; is the value of the corresponding index in the preset seed number set, otherwise the grid seed number is assigned .

[0079] In the step 13, the internal force is extracted: the internal force data report is extracted on the path, and the node position, time, and component force value in the report are parsed through a regular expression to further extract the data of the bending moment, axial force, and shear force.

[0080] In order to accurately extract the boundary nodes and internal nodes of the pipe joint position, first, the specified nodes need to be screened through the formula as follows:

[0081]

[0082]

[0083] In the formula, represents the coordinate of the node, represents the scaling factor of the coordinate, represents the width of the pipe, represents a small constant , represents the coordinate of the node, represents the coordinate of the node, represents the outer radius of the pipe, represents the radius offset, represents a small constant .

[0084] For the nodes meeting the above two conditions, their angles are calculated according to the following formula:

[0085]

[0086] Subsequently, two boundary nodes can be defined:

[0087]

[0088]

[0089] In the formula, represents the node set meeting the above two screening conditions.

[0090] Finally, the internal node set is obtained:

[0091]

[0092] ​Boundary node set:

[0093]

[0094] After obtaining the internal node and boundary node set, the corresponding mechanical data is extracted by using the session.xyDataListFromField() function of Abaqus.

[0095] Preferably, in the prior art, different working condition simulations need to manually modify parameters and submit calculations one by one, which is inefficient and prone to errors. Currently, in the analysis of shield tunnel segments, in the face of water pressure, longitudinal pressure and main and secondary tension forces, etc. Multiple working conditions, it is usually necessary to manually modify the load parameters one by one and save them as new INP calculation files, and then manually click to submit the job before calculation. This working method is extremely tedious, and the workload increases exponentially when performing parameter sensitivity analysis or multi-working condition simulation tasks, which is prone to cause parameter configuration errors, naming confusion, file omission, etc., which seriously restricts the simulation efficiency and modeling controllability.

[0096] Therefore, the present application is based on a parameterized INP generation and automatic calculation scheduling system; the present patent proposes an automatic simulation control scheme, which can batch generate INP input files of different working condition combinations and automatically complete task submission and calculation scheduling. In the model building stage, the modeling parameters such as water pressure, longitudinal load, main and secondary tension forces are abstracted as variables through Python scripts, and users only need to specify the minimum value, maximum value and step size to realize the parameterized generation of multiple working conditions, and generate the corresponding INP calculation file and then submit the task for calculation.

[0097] The present patent proposes a complete set of automatic modeling and simulation whole-process automation method for shield segment structure modeling process, which is based on Python script as the core and Abaqus secondary development platform. The present patent changes the traditional modeling process which depends on Abaqus graphical user interface and gradually interacts into automatic modeling mode driven by whole-process script, realizes complete closed loop of automatic input of structure parameters, automatic construction of geometric structure, automatic execution of Boolean cutting, automatic deployment of meshing. The test data shows that the traditional operation modeling of a single three-ring structure takes more than 10 hours, while the complete model generated by the method of the present patent only takes 3-5 minutes, and the efficiency is significantly improved. At the same time, the present patent technology supports parameterized modification, and can easily adapt to various structures with different radii, concave-convex tenon arrangement and bolt angles, greatly enhancing the reusability and engineering universality of the model.

[0098] In the present application, only the middle ring structure modeling mode: by setting the parameter middle_ring_only as True, only the middle ring structure composed of standard blocks, adjacent blocks and capping blocks can be automatically generated, and the cutting and Boolean operation process of the front ring and the rear ring is skipped. It is suitable for single-ring rapid modeling, internal force analysis and parameter sensitivity research of structures that do not need to consider the interference of three-ring mortise and tenon; as shown in Figure 2 In the present application, automatic angle adaptation under various mortise and tenon arrangement schemes: under the premise of keeping the total number of mortise and tenon arrangements the same, the user can flexibly set the number of mortise and tenon on the standard block, the adjacent block and the capping block; as shown in Figure 4 and Figure 5 For example, the capping block is changed from 2 mortise and tenon to 3 mortise and tenon, which will automatically re-derive the circumferential angle relationship, improving the engineering adaptability.

[0099] The above is only the preferred embodiment of the present application, it should be noted that for those skilled in the art, without departing from the principles of the present application, can make a number of improvements and refinements, these improvements and refinements should also be considered within the scope of the present application. The structures, devices and operation methods not specifically described and explained in the present application, if not specifically described and limited, are implemented according to the conventional means in the art.

Claims

1. A method for intelligent modeling and analysis of distributed tenon-type segment lining structures in shield tunnels, characterized by: Includes the following steps: Step 1: Parameter definition and validity check: Define structural parameters through script input and set up a parameter validity check mechanism. If the structural closure or layout rules are not met, the program will exit directly. Step 2: Automatic generation of module structure: Generate 3D solid models of standard blocks, adjacent blocks and capping blocks according to parameters, and combine them to form a middle ring structure. Then, automatically copy and cut them to form the front ring and rear ring structures. Step 3: Arrangement of tenon and bolt structure: Automatically construct groove body, tenon body and bolt structure, and generate front, middle and rear ring tenon and bolt joint and anchoring structure through geometric Boolean and rotation arrangement operations, and complete the automatic combination of connecting components; Step 4: Loading Component Generation and Placement: Generate loading beams and place them on the top and sides of the structure; Step 5: Structure mesh generation: Automatically generate meshes according to the block type using preset block mesh seed density parameters to achieve automated mesh generation; Step 6: Boundary conditions and load definition: Automatically set material properties, interactions, constraints and boundary conditions, and define radial compressive loads and tensile loads; Step 7: Output the complete model structure: Obtain the complete 3D solid model, which can be directly used for subsequent simulation analysis; Step 8: Submit the task and perform batch calculations: Generate multiple different working condition combinations, submit the task for calculation, and obtain multiple ODB result files; Step 9: Read ODB files: Automatically traverse all ODB result files in the specified directory and read the finite element simulation information from them; Step 10: Read the seam opening amount: Automatically identify the seam node pairs according to the model structure, obtain the U1 and U2 displacement components of these nodes, calculate the change in Euclidean distance between each pair of nodes, and output the seam opening amount data. Step 11: Read structural deformation: Automatically select feature nodes, extract their displacement information, combine with the original radius data, and output as structural deformation data; Step 12: Generate and add paths: Automatically generate circumferential paths in the radial direction at key structural sections; Step 13: Extract internal forces: Extract the internal force data report on the path, and parse the node positions, time, and force values ​​of each component in the report using regular expressions. Then extract the data of bending moment, axial force, and shear force. Step 14: Results Output and Archiving: All extracted data is output in CSV file format and named according to structure name and variable type, and automatically classified and archived in a folder named after the ODB output file name.

2. The intelligent modeling and analysis method for distributed tenon-type segment lining structure of shield tunnels according to claim 1, characterized in that: By inputting the number of tenons and mortises on each type of segment, the system automatically calculates the circumferential angles of the capping block, standard block, and adjacent blocks, achieving intelligent angle configuration for structural segmentation. Specifically, the user inputs the number of tenons and mortises on each type of block and the total number of standard blocks in the entire ring, which automatically calculates the angles for each type of block. The number of tenons and mortises includes those on the standard block, adjacent block, and capping block; the formula is shown below: ； ； ； in, This represents the angle corresponding to the standard block. This represents the angle corresponding to the capping block. This represents the angle corresponding to the adjacent block. The number of tenons and mortises in the standard block. This represents the number of tenons and mortises in the capping block. This represents the number of tenons and mortises in adjacent blocks. Represents the total number of standard blocks; Automatically calculated , and It will be used as a control parameter in the geometric modeling function for drawing geometric sketches and arranging parts, so as to realize the automatic layout of segment blocks.

3. The intelligent modeling and analysis method for distributed tenon-type segment lining structure of shield tunnels according to claim 2, characterized in that: A structural rationality verification mechanism is introduced into the automatic segment modeling system. This mechanism determines whether the total circumferential angles of the structure meet the closure condition, thereby determining whether the current block and tenon / concave-convex tenon arrangement scheme can form a complete structure. At the same time, it checks whether the preset angle distribution of all tenons / concave-convex tenons conflicts with the key joint positions, ensuring that the tenon / concave-convex tenon elements reasonably avoid each other in the ring. In the parameter setting stage, an automatic variable expansion method is adopted. Based on the key parameters provided by the user, including the number of standard blocks, the number of tenons / concave-convex tenons, and the block width, the system automatically derives the total number of tenons / concave-convex tenons, the angle occupied by each block, and the tenon / concave-convex tenon interval angle, and embeds multiple logical judgments and constraint rules.

4. The intelligent modeling and analysis method for distributed tenon-type segment lining structure of shield tunnels according to claim 3, characterized in that: By inputting a small number of control parameters, all valid joint nodes are automatically selected, and corresponding opening change data files are generated. Control parameters such as ring layer, radius, and tenon / parietal arrangement; First, it is necessary to define the screening criteria for seam nodes, including: (A) The coordinates satisfy the tolerance judgment at a specific Z-level; (B) The coordinate radius is close to the inner / outer annular radius; (C) The circumferential angle satisfies the closed interval of the actual tongue and groove tenon structure arrangement.

5. The intelligent modeling and analysis method for distributed tenon-type segment lining structure of shield tunnels according to claim 1, characterized in that: In step 1, parameter definition and rationality verification: structural parameters are defined by inputting them through a script; structural parameters include the outer diameter of the segment, the number of segments, the position of the tenon and groove, and the bolt parameters, and a parameter validity judgment mechanism is set. If the structural closure or layout rules are not met, the process will exit directly. When checking the closure of a structure, the sum of angles is calculated using the formula shown below. : ； in, This represents the angle corresponding to the standard block. This represents the angle corresponding to the capping block. Represents the total number of standard blocks; if If the current structure is deemed invalid, the modeling process will automatically exit. When checking the arrangement rules of the tenon and mortise joints, the following formula is used for calculation; For the A tenon and a mortise ( ), its angular position It can be calculated using the formula: ； in, This refers to the total number of tenons and mortises. The reference angle representing the first tenon joint. Represents the spacing between the male and female tenons; Set the set of joint angles ;in Given the total number of joints, the angle position of each tenon and mortise is required. and the angle of each seam satisfy: ； For all tongue and groove joints and all seams For a pair of pairs to hold, the above inequality must hold; if there exists a pair ( If the inequality is invalid, the current structure is considered invalid and the modeling process will automatically exit.

6. The intelligent modeling and analysis method for distributed tenon-type segment lining structure of shield tunnels according to claim 1, characterized in that: In step 5, the structural mesh is generated automatically according to the block type by using the preset block mesh seed density parameters, thus achieving automated mesh generation. When dividing the grid, different line segment types are identified to achieve grid division based on a preset seed density. For standard blocks, the seed number is allocated using the piecewise constant function shown below: ； in, The constant representing the infinitesimal is taken as , Represents the length of the line segment. This represents the length of the top and bottom long sides of the tenon and mortise. Represents the radius of the curved edge of the tenon and mortise joint. This represents the thickness of the tenon and mortise joint. This represents the amount of shrinkage in the top radius of the tenon and mortise joint. Represents the width of the concrete tunnel segment. This represents the number of grid seeds corresponding to the upper and lower long sides of the tenon and mortise. This represents the number of mesh seeds corresponding to the curved edge of the tenon and mortise. This represents the number of grid seeds corresponding to the bevel of the tenon and mortise. This represents the number of grid seeds corresponding to the edges of the concrete segments parallel to the tunneling direction; For edges that do not meet the above conditions, it is necessary to further search the preset distance set. If it exists satisfy: ； Then allocate the number of grid seeds. ; The corresponding index in the preset seed set The value is used to determine the number of grid seeds; otherwise, the number of seeds is allocated. .

7. The intelligent modeling and analysis method for distributed tenon-type segment lining structure of shield tunnels according to claim 1, characterized in that: In step 13, internal forces are extracted: an internal force data report is extracted from the path, and the node positions, time, and force values ​​of each component in the report are parsed using regular expressions to further extract the data of bending moment, axial force, and shear force. First, you need to filter the specified nodes using the formula shown below: ； ； In the formula, Representing nodes coordinate, represent The scaling factor of the coordinates. Represents the width of the tunnel segments. The constant representing the infinitesimal is taken as , Representing nodes coordinate, Representing nodes coordinate, Represents the outer radius of the tunnel segment. This represents the radius offset. The smaller constant is taken ; For nodes that satisfy the above two conditions, their angles are calculated using the following formula: ； Then, two boundary nodes can be defined: ； ； In the formula, This represents the set of nodes that satisfy the two filtering conditions mentioned above. Finally, we obtain the set of internal nodes: ； Set of boundary nodes: ； After obtaining the set of internal nodes and boundary nodes, the corresponding mechanical data is extracted using the Abaqus session.xyDataListFromField() function.

8. The intelligent modeling and analysis method for distributed tenon-type segment lining structure of shield tunnel according to claim 1, characterized in that: A parametric INP generation and automatic calculation scheduling system; batch generation of INP input files for different working condition combinations, and automatic completion of task submission and calculation scheduling; during the model building stage, modeling parameters such as water pressure, longitudinal load, main tension force and secondary tension force are abstracted into variables through Python scripts. Users only need to specify the minimum value, maximum value and step size to achieve parametric generation of multiple working conditions, generate corresponding INP calculation files and then submit the task for calculation.

Images