A large deformation finite element model construction method, device, terminal and medium

Abstract

The present invention relates to the technical field of large deformation model construction, and in particular to a large deformation finite element model construction method, device, terminal and medium, the method comprising: after assigning material parameters to the initial model, assembling the model and meshing, obtaining an initial boundary point array; performing the first pile penetration loading on the initial model, updating the initial boundary point array, reconstructing the initial model, and obtaining a first deformation model; after assigning material parameters to the first deformation model, assembling the model and meshing, obtaining a second boundary point array; performing the second pile penetration loading on the first deformation model, updating the second boundary point array, reconstructing the first deformation model, and obtaining a second deformation model; and repeating this cycle until a target deformation model with a predetermined deformation degree is obtained. The present application records the initial boundary point array when the model is not deformed, and by continuously updating the initial boundary point array, accurately records the changes in the model boundary, thereby improving the accuracy of the reconstructed model.

Description

Technical Field

[0001] The present invention relates to the field of large deformation model construction, and in particular to a large deformation finite element model construction method, device, terminal and medium. Background Art

[0002] The finite element method is a commonly used calculation method in scientific research and engineering, which can effectively analyze the stress and deformation of the structure. The classic finite element method mainly includes seven parts: model establishment, material parameter assignment, model assembly, analysis step setting, boundary condition setting, mesh division and calculation analysis. Taking pile penetration analysis as an example, the analysis process is as follows: Figure 1 However, when performing pile penetration analysis, when the pile penetration depth is small, the conventional finite element method can obtain better analysis results, such as Figure 2 As shown in Figure 2, when the pile penetration depth is large, the mesh deformation at the pile bottom and pile side is too large, such as Figure 3 As shown in the figure, conventional finite element methods cannot accurately analyze soil stress and deformation, and cannot perform pile penetration analysis. The actual pile construction penetration depth can reach tens of meters, and the mesh distortion at the pile bottom is more serious. Conventional finite element methods cannot complete pile penetration analysis.

[0003] Large deformation finite element analysis methods solve the problem of large deformation calculation, such as the arbitrary Lagrang-Euler coupling method, large deformation finite element analysis based on mesh remeshing and stress interpolation (RITSS) or particle finite element (PFEM) method, etc. In the RITSS method, large deformation is divided into several small deformations, and large deformation analysis is achieved by accumulating several small deformations. The analysis process of the RITSS method is as follows: Figure 4 As shown, the key to implementing this method is to reconstruct and update the deformed model after an analysis. Therefore, an accurate and efficient model reconstruction method is particularly important for analysis and calculation.

[0004] In the current model reconstruction method, when determining the position coordinates of the boundary points, an approximate fitting process is used, which results in an inaccurate reconstructed model.

[0005] Therefore, the prior art has defects and needs to be improved and developed. Summary of the invention

[0006] The present application provides a large deformation finite element model construction method, device, terminal and medium to solve the technical problem in the related art that when determining the position coordinates of the boundary points, it is an approximate fitting process, resulting in an inaccurate reconstructed model.

[0007] To achieve the above objectives, this application adopts the following technical solutions:

[0008] A large deformation finite element model construction method, wherein the method comprises:

[0009] Establishing an initial model, assigning material parameters to the initial model, assembling the model and dividing the mesh, and obtaining an initial boundary point array;

[0010] Performing a first pile penetration loading on the initial model, updating the initial boundary point array, and reconstructing the initial model according to the updated initial boundary point array to obtain a first deformation model;

[0011] After assigning material parameters, assembling the model and dividing the mesh on the first deformation model, a second boundary point array is obtained;

[0012] Performing a second pile penetration loading on the first deformation model, updating the second boundary point array, and reconstructing the first deformation model according to the updated second boundary point array to obtain a second deformation model;

[0013] This cycle is repeated until the target deformation model with a predetermined deformation degree is obtained.

[0014] In one embodiment of the present application, an initial model is established, and after assigning material parameters, model assembly and meshing the initial model, an initial boundary point array is obtained, including:

[0015] Establishing an initial model, assigning material parameters to the initial model, assembling the model and dividing the mesh, and obtaining all initial mesh nodes;

[0016] Searching all initial mesh nodes for an initial mesh node on a model boundary of the initial model;

[0017] An initial boundary point array is formed based on all initial mesh nodes on the model boundary of the initial model.

[0018] In one embodiment of the present application, searching among all initial grid nodes for an initial grid node on a model boundary of the initial model includes:

[0019] Using each corner point of the initial model as an initial boundary point of the initial model, and using the model boundary formed between each initial boundary point as a target straight line;

[0020] Screening the points to be searched from all the initial grid nodes, and determining the positional relationship between the points to be searched and the target straight line;

[0021] If the point to be searched is on the target straight line, it is determined that the initial grid node is an initial grid node on the model boundary of the initial model.

[0022] In one embodiment of the present application, screening the to-be-searched points in all initial grid nodes and determining the positional relationship between the to-be-searched points and the target straight line includes:

[0023] Obtaining a first horizontal coordinate, a first vertical coordinate, a second horizontal coordinate, and a second vertical coordinate corresponding to two initial boundary points of the target straight line, and determining a boundary straight line equation corresponding to the target straight line;

[0024] Filter out points to be searched whose abscissas are between the first abscissa and the second abscissa and whose ordinates are between the first ordinate and the second ordinate from all initial grid nodes;

[0025] Calculate the distance between the point to be searched and the boundary line equation according to the horizontal coordinate and the vertical coordinate of the point to be searched;

[0026] If the distance between the point to be searched and the boundary line equation is less than a preset value, the positional relationship between the point to be searched and the target line is that the point to be searched is on the target line;

[0027] If the distance from the point to be searched to the boundary line equation is greater than or equal to a preset value, the positional relationship between the point to be searched and the target line is that the point to be searched is not on the target line.

[0028] In one embodiment of the present application, the initial model is subjected to a first pile penetration load, the initial boundary point array is updated, and the initial model is reconstructed according to the updated initial boundary point array to obtain a first deformation model, including:

[0029] After the initial model is subjected to the first pile penetration loading, a first analysis result is obtained;

[0030] updating the initial boundary point array according to the first analysis result to obtain an updated initial boundary point array;

[0031] The initial model is reconstructed according to the updated initial boundary point array to obtain a first deformed model.

[0032] In one embodiment of the present application, after assigning material parameters, assembling the model and meshing the first deformation model, a second boundary point array is obtained, including:

[0033] After assigning material parameters, assembling the model and meshing the first deformation model, all second mesh nodes are obtained;

[0034] Searching for a second mesh node on a model boundary of the first deformed model among all second mesh nodes;

[0035] A second boundary point array is formed based on all second mesh nodes on the model boundary of the first deformation model.

[0036] In one embodiment of the present application, a second pile penetration loading is performed on the first deformation model, the second boundary point array is updated, and the first deformation model is reconstructed according to the updated second boundary point array to obtain a second deformation model, including:

[0037] After the first deformation model is subjected to a second pile penetration loading, the first deformation model is deformed, and a second analysis result is obtained;

[0038] updating the second boundary point array according to the second analysis result to obtain an updated second boundary point array;

[0039] The first deformation model is reconstructed according to the updated second boundary point array to obtain a second deformation model.

[0040] The present application also provides a large deformation finite element model construction device, wherein the device comprises:

[0041] The first model processing module is used to establish an initial model, assign material parameters to the initial model, assemble the model and divide the mesh to obtain an initial boundary point array;

[0042] A first updating module, configured to perform a first pile penetration loading on the initial model, update the initial boundary point array, and reconstruct the initial model according to the updated initial boundary point array to obtain a first deformation model;

[0043] A second model processing module, configured to obtain a second boundary point array after assigning material parameters, model assembly and meshing the first deformation model;

[0044] A second updating module, configured to perform a second pile penetration loading on the first deformation model, update the second boundary point array, and reconstruct the first deformation model according to the updated second boundary point array to obtain a second deformation model;

[0045] The loop module is used to loop in this way until a target deformation model with a predetermined deformation degree is obtained.

[0046] The present application also provides a terminal, which includes: a memory, a processor, and a large deformation finite element model construction program stored in the memory and executable on the processor, wherein the large deformation finite element model construction program implements the steps of the large deformation finite element model construction method as described above when executed by the processor.

[0047] The present application also provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and the computer program can be executed to implement the steps of the large deformation finite element model construction method as described above.

[0048] Beneficial effects of the present invention: The method of the embodiment of the present invention establishes an initial model, assigns material parameters to the initial model, assembles the model and divides the mesh to obtain an initial boundary point array; performs the first pile penetration loading on the initial model, updates the initial boundary point array, and reconstructs the initial model according to the updated initial boundary point array to obtain a first deformation model; performs the first deformation model assignment of material parameters, assembles the model and divides the mesh to obtain a second boundary point array; performs the second pile penetration loading on the first deformation model, updates the second boundary point array, and reconstructs the first deformation model according to the updated second boundary point array to obtain a second deformation model; and repeats this cycle until a target deformation model with a predetermined deformation degree is obtained. The present application records the initial boundary point array when the model is not deformed, and by continuously updating the initial boundary point array, accurately records the changes in the model boundary, and there is no approximate fitting process, thereby improving the accuracy of the reconstructed model. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] Figure 1 It is a classic finite element analysis flow chart.

[0050] Figure 2 It is a schematic diagram of pile penetration in a normal grid.

[0051] Figure 3 It is a schematic diagram of pile penetration using a distorted grid.

[0052] Figure 4 It is a flow chart of large deformation finite element analysis.

[0053] Figure 5 It is a schematic diagram of slit retrieval.

[0054] Figure 6 It is a flow chart of a preferred embodiment of the method for constructing a large deformation finite element model in the present invention.

[0055] Figure 7 It is a schematic diagram of the steps of a preferred embodiment of the method for constructing a large deformation finite element model in the present invention.

[0056] Figure 8 It is a schematic diagram of boundary point modeling of a preferred embodiment of the large deformation finite element model construction method in the present invention.

[0057] Fig. 9 It is a schematic diagram of the model reconstructed three times in the present invention.

[0058] Fig.10 It is a schematic diagram of the model reconstructed 10 times in the present invention.

[0059] Fig.11 It is a schematic diagram of the model reconstructed 20 times in the present invention.

[0060] Fig.12 It is a functional principle block diagram of a preferred embodiment of the large deformation finite element model building device in the present invention.

[0061] Fig.13 It is a functional principle block diagram of a preferred embodiment of the terminal in the present invention. DETAILED DESCRIPTION

[0062] In order to make the purpose, technical solution and advantages of the present invention clearer and more specific, the present invention is further described in detail below with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0063] At present, the model reconstruction methods include spline curve method. For example, the three-dimensional model of each analysis is constructed by fitting the deformed free surface according to the B-spline curve. Alternatively, the alpha shape method is used to search the outer boundary of the node and then reconstruct the model. The B-spline curve fitting and alpha shape methods need to calculate all points during the search to determine the boundary points. As the calculation progresses, the grid is continuously encrypted, and the calculation time of this search method will increase significantly. For models with slits, such as Figure 5 As shown in FIG. 1 , when the rolling circle diameter in the alpha shape method is large, the slit will be ignored, resulting in the constructed model not being consistent with the actual model and causing calculation errors. Therefore, the existing model reconstruction method still has the problems of low efficiency and poor effect.

[0064] The following describes the large deformation finite element model construction method, device, terminal and medium of the embodiment of the present application with reference to the accompanying drawings. In view of the problem that in the related art mentioned in the above background technology, when determining the position coordinates of the boundary points, it is an approximate fitting process, which leads to the problem that the reconstructed model is not accurate, the present application provides a large deformation finite element model construction method, in which an initial model is established, and the initial model is assigned material parameters, model assembly and mesh division to obtain an initial boundary point array; the initial model is subjected to the first pile penetration loading, the initial boundary point array is updated, and the initial model is reconstructed according to the updated initial boundary point array to obtain a first deformation model; the first deformation model is subjected to the second pile penetration loading, the second boundary point array is updated, and the first deformation model is reconstructed according to the updated second boundary point array to obtain a second deformation model; this cycle is repeated until a target deformation model with a predetermined deformation degree is obtained. The present application records the initial boundary point array when the model is not deformed, and accurately records the changes in the model boundary by continuously updating the initial boundary point array, and there is no approximate fitting process, thereby improving the accuracy of the reconstructed model.

[0065] See also Figure 6 The large deformation finite element model construction method described in the embodiment of the present invention comprises the following steps:

[0066] Step S100, establishing an initial model, assigning material parameters to the initial model, assembling the model and dividing the mesh, and obtaining an initial boundary point array.

[0067] like Figure 7 As shown in (a), the initial model is established and the initial model boundary is obtained [B] initinal Specifically, in the embodiment of the present application, since the rigidity of the pile is much greater than that of the soil, the pile is regarded as a rigid body without deformation, and only the deformation of the soil is considered. The present application uses finite element software to process each step.

[0068] In the embodiment of the present application, the step S100 specifically includes:

[0069] Step S110, establishing an initial model, assigning material parameters to the initial model, assembling the model and meshing the model to obtain all initial mesh nodes;

[0070] Step S120, searching among all initial grid nodes for an initial grid node on the model boundary of the initial model;

[0071] Step S130: forming an initial boundary point array based on all initial grid nodes on the model boundary of the initial model.

[0072] like Figure 7 As shown in (b), grid division, search for boundary points, and form an initial boundary point array [B] 1 Specifically, since the model boundary coordinates are known when the model is established, the model boundary can be directly obtained according to the model boundary coordinates, which is denoted as [B] initinal After the material parameters are assigned and the model is assembled, the model enters the meshing step. During meshing, grid points are formed inside the model and on the boundary (internal grid points do not affect the model construction and can be ignored). The grid points on the boundary are searched to form an initial boundary point array [B] 1 .

[0073] The embodiment of the present application records the initial boundary point array when the model is not deformed, and then continuously updates the initial boundary point array to accurately record the changes in the model boundary and improve the accuracy of the reconstructed model.

[0074] In one embodiment of the present application, step S120 specifically includes:

[0075] Step S121, taking each corner point of the initial model as an initial boundary point of the initial model, and taking a model boundary formed between each initial boundary point as a target straight line;

[0076] Step S122, screening the points to be searched in all the initial grid nodes, and determining the positional relationship between the points to be searched and the target straight line;

[0077] Step S123: If the point to be searched is on the target straight line, it is determined that the initial grid node is an initial grid node on the model boundary of the initial model.

[0078] Specifically, the boundary point search method of the present application is: first search for boundary points based on the initial model boundary. The coordinates of the initial boundary points are known. After the first meshing, some grid points will be on the boundary. These grid points on the boundary need to be searched. The method of searching boundary points is to determine which points are on the boundary. Unlike other methods, the present application can accurately determine the boundary points by introducing an array of initial boundary points with known positions to retrieve the deformed grid. Figure 8 It can be seen that the initial boundary point is K 1 , K 2 , K 3 and K 4 To determine which of the added grid points are boundary points, we only need to determine which nodes are on the straight line formed by the initial boundary points. 2 With K 3 is the initial boundary point. To determine whether a grid point is a boundary point, it is necessary to determine whether the grid point is within K2 and K 3 To reduce the number of points that need to be judged, Figure 8 It can be seen that for the possible K 2 and K 3 The points on the straight line will fall within the narrow band shown by the dotted line, so we can use K 2 and K 3 The coordinates of determine the points in the narrow band, such as Figure 8 G 1 The point is obviously not in the narrow band, and according to K 2 and K 3 The coordinates are directly excluded. This step can exclude most points and reduce subsequent judgment work.

[0079] The embodiment of the present application excludes most of the points by screening the points to be searched, that is, when retrieving boundary points, only the points near the boundary are searched instead of the global points. Fewer points are calculated, which saves time and improves search efficiency.

[0080] In the embodiment of the present application, the step S122 specifically includes:

[0081] Step S1221, obtaining a first horizontal coordinate, a first vertical coordinate, a second horizontal coordinate, and a second vertical coordinate corresponding to two initial boundary points of the target straight line, and determining a boundary straight line equation corresponding to the target straight line;

[0082] Step S1222, selecting points to be searched whose abscissas are between the first abscissa and the second abscissa and whose ordinates are between the first ordinate and the second ordinate from all initial grid nodes;

[0083] Step S1223, calculating the distance from the point to be searched to the boundary line equation according to the horizontal coordinate and the vertical coordinate of the point to be searched;

[0084] Step S1224: if the distance from the point to be searched to the boundary line equation is less than a preset value, the positional relationship between the point to be searched and the target line is that the point to be searched is on the target line;

[0085] Step S1225: If the distance from the point to be searched to the boundary line equation is greater than or equal to a preset value, the positional relationship between the point to be searched and the target line is that the point to be searched is not on the target line.

[0086] Specifically, firstly, a simple judgment is made based on the coordinates of the initial grid nodes, such as Figure 8As shown in the figure, if the coordinates to be searched are between the two initial boundary points, then its horizontal and vertical coordinates should also be between the two initial boundary points. Therefore, firstly, delete the points that do not meet this condition and filter out the points to be searched to reduce the number of initial grid nodes that need to be judged. Then, the boundary line equation is calculated based on the coordinates of the two initial boundary points, and then the distance between the point to be searched and the line is calculated. If the line distance is 0, it means that this point is on the boundary and belongs to the boundary point. After judging all the points, the initial boundary point array is formed [B] 1 It should be pointed out that due to the existence of errors, the distance between the search point and the straight line is often difficult to be accurately zero, so a distance error is introduced. When the distance error is less than the preset value, the point is determined to be a boundary point. The introduction of this error can also adjust the search accuracy. The size of the distance error determines the width of the narrowband. When the distance error is large, the narrowband is wider, and when the distance error is small, the narrowband is also smaller.

[0087] ;

[0088] ;

[0089] in, are the coordinates of the point to be searched, is the coordinate of one of the initial boundary points, For The coordinates of another initial boundary point on the same straight line. A, B, and C are the parameters of the straight line formed by the initial boundary points, and d represents the distance between the point to be searched and the straight line.

[0090] The embodiment of the present application determines whether the point to be searched is on the model boundary by judging the difference between the distance from the point to be searched to the boundary line equation and a preset value, thereby improving the judgment efficiency.

[0091] like Figure 6 As shown, the large deformation finite element model construction method also includes the following steps:

[0092] Step S200: performing a first pile penetration loading on the initial model, updating the initial boundary point array, and reconstructing the initial model according to the updated initial boundary point array to obtain a first deformation model.

[0093] In the embodiment of the present application, the step S200 specifically includes:

[0094] Step S210, after performing the first pile penetration loading on the initial model, a first analysis result is obtained;

[0095] Step S220: updating the initial boundary point array according to the first analysis result to obtain an updated initial boundary point array;

[0096] Step S230: reconstruct the initial model according to the updated initial boundary point array to obtain a first deformed model.

[0097] like Figure 7 As shown in (c), the numerical simulation is loaded, the model shape changes, and the updated initial boundary point array is obtained [B] 1 new . Specifically, after the grid division is completed, the embodiment of the present application performs pile penetration loading (calculation and analysis step). At this time, the pile will move downward to squeeze the soil, and the soil will deform under the squeezing of the pile. Due to the compaction and shear expansion characteristics of the soil, the soil at the bottom of the pile moves downward after being compressed by the pile, while the soil on the side of the pile will bulge upward. At this time, the deformation of the soil is unpredictable. After the calculation and analysis is completed, the grid is deformed, and the grid nodes will move accordingly, and the coordinates of the moved grid nodes are recorded. That is, the first analysis result includes the coordinates of each grid node after the initial model is deformed; the initial boundary point array [B] is calculated according to the coordinates of each grid node. 1 The node coordinates within are updated to form an updated initial boundary point array [B] 1 new . Based on the updated initial boundary point array [B] 1 new Reconstruct the model and get the first deformation model, such as Figure 7 As shown in (d) in .

[0098] The present application continuously updates the initial boundary point array, accurately records the changes in the model boundary, and does not have an approximate fitting process, thereby improving the accuracy of the reconstructed model.

[0099] like Figure 6 As shown, the large deformation finite element model construction method also includes the following steps:

[0100] Step S300: After assigning material parameters, assembling the model and dividing the mesh, the first deformation model is subjected to a second boundary point array.

[0101] In one embodiment of the present application, step S300 specifically includes:

[0102] Step S310, after assigning material parameters, model assembly and meshing the first deformation model, all second mesh nodes are obtained;

[0103] Step S320: searching for a second mesh node on the model boundary of the first deformed model among all second mesh nodes;

[0104] Step S330: forming a second boundary point array based on all second mesh nodes on the model boundary of the first deformation model.

[0105] like Figure 7As shown in (e), the first deformed model is meshed and new boundary points are searched to form a second boundary point array [B] 2 Specifically, the first deformation model is meshed again after material parameter assignment, model assembly and other steps. At this time, the mesh points on the boundary will change, and the boundary mesh points need to be identified. Based on array [B] 1 new Node coordinates, identify the re-divided grid nodes, and determine the second boundary point array [B] 2 .

[0106] The present application accurately records the changes in the model boundary by continuously updating the boundary point array, and there is no approximate fitting process, thereby improving the accuracy of the reconstructed model.

[0107] In one embodiment, the step S320 specifically includes:

[0108] Step S321, taking the line between adjacent corner points of the first deformation model as the target straight line;

[0109] Step S322, screening the points to be searched in all the second grid nodes, and determining the positional relationship between the points to be searched and the target straight line;

[0110] Step S323: If the point to be searched is on the target straight line, determine that the second mesh node is a second mesh node on the model boundary of the first deformation model.

[0111] The step S322 specifically includes:

[0112] Saving all the second grid nodes to the array to be determined, and comparing all the second grid nodes in the array to be determined with the updated initial boundary point array;

[0113] If the coordinates of the current second mesh node are the same as the node coordinates in the updated initial boundary point array, determining that the current second mesh node is a target boundary point on the model boundary of the first deformation model;

[0114] Using the second grid nodes in the array to be determined except the target boundary point as the first screening second grid nodes;

[0115] Obtaining the third abscissa, the third ordinate, the fourth abscissa and the fourth ordinate of the two corner points corresponding to the target straight line, and determining the boundary straight line equation corresponding to the target straight line;

[0116] Screen out the points to be searched whose abscissas are between the third abscissa and the fourth abscissa and whose ordinates are between the third ordinate and the fourth ordinate from all the preliminarily screened second grid nodes;

[0117] Calculate the distance between the point to be searched and the boundary line equation according to the horizontal coordinate and the vertical coordinate of the point to be searched;

[0118] If the distance between the point to be searched and the boundary line equation is less than a preset value, the positional relationship between the point to be searched and the target line is that the point to be searched is on the target line;

[0119] If the distance from the point to be searched to the boundary line equation is greater than or equal to a preset value, the positional relationship between the point to be searched and the target line is that the point to be searched is not on the target line.

[0120] Specifically, when determining the boundary points of the first deformation model, if the second grid node has the same node coordinates as those in the updated initial boundary point array, it can be directly determined that the point is the second grid node on the model boundary of the first deformation model, because the grid nodes in the updated initial boundary point array are all boundary points, thereby improving the determination efficiency.

[0121] like Figure 6 As shown, the large deformation finite element model construction method also includes the following steps:

[0122] Step S400: performing a second pile penetration loading on the first deformation model, updating the second boundary point array, and reconstructing the first deformation model according to the updated second boundary point array to obtain a second deformation model.

[0123] In the embodiment of the present application, the step S400 specifically includes:

[0124] Step S410: After the first deformation model is subjected to a second pile penetration loading, the first deformation model is deformed to obtain a second analysis result;

[0125] Step S420: updating the second boundary point array according to the second analysis result to obtain an updated second boundary point array;

[0126] Step S430: reconstruct the first deformation model according to the updated second boundary point array to obtain a second deformation model.

[0127] like Figure 7 As shown in (f), the numerical simulation is loaded, the model shape changes, and the updated second boundary point array is obtained [B] 2 new Specifically, the first deformation model is subjected to a second pile penetration load, the pile continues to be inserted into the soil, the soil deforms, and the boundary points continue to move. After the calculation is completed, the second boundary point array is updated according to the second analysis result to obtain an updated second boundary point array [B] 2 new .

[0128] like Figure 7As shown in (g), the model is reconstructed according to the updated second boundary point array to obtain the second deformed model.

[0129] The present application adopts a node tracking method to retrieve the model boundary points according to the initial model boundary of the initial model, form a model boundary array, and record the model boundary points. When the model boundary needs to be reconstructed, the model boundary is retrieved according to the model boundary array. At this time, only the points near the boundary need to be searched, which greatly reduces the search points and can effectively improve the search and modeling speed. In addition, the present application accurately records the changes in the model boundary by continuously updating the boundary point array, and the model boundary with tiny slits can still be accurately identified and constructed.

[0130] like Figure 6 As shown, the large deformation finite element model construction method also includes the following steps:

[0131] Step S500: repeat this process until a target deformation model with a predetermined deformation degree is obtained.

[0132] like Figure 7 As shown in (h), the second deformed model is meshed and new edge connections are searched to form an updated third boundary point array [B] 3 Specifically, the second deformation model is subjected to a third penetration loading and meshing, and the search process is repeated to obtain an updated third boundary point array [B] 3 As the pile is penetrated, the pile goes deeper and deeper into the soil. At this time, the traditional method has serious grid distortion and cannot complete the analysis. The embodiment of the present application repeats the above process cyclically to reconstruct the model and realize large deformation simulation. The method of the present application can effectively reconstruct the model, such as Fig. 9 , Fig.10 and Fig.11 As shown in the figure, the constructed soil model increases with the penetration depth, and the subsidence arc formed on the surface due to the pile penetration and extrusion can also be well displayed.

[0133] The method proposed in the present invention is applicable to the model reconstruction process of RITSS or PFEM method in large deformation analysis. In the past, the method used in the model was relatively cumbersome. The method of the present application records the boundary point information when the model is not deformed, tracks the movement of the point during deformation, updates the position coordinates, and realizes the accurate position tracking of the boundary node. The position coordinates obtained in this process are error-free, and there is no approximate fitting process. This method also improves the accuracy of the stress transfer process. Moreover, when searching for boundary points, only the points near the boundary are searched, not the global point judgment, and fewer points are calculated, saving time.

[0134] In one embodiment, if Fig.12As shown, based on the above large deformation finite element model construction method, the present invention also provides a large deformation finite element model construction device, including:

[0135] The first model processing module 100 is used to establish an initial model, assign material parameters to the initial model, assemble the model and divide the mesh to obtain an initial boundary point array;

[0136] A first updating module 200 is used to perform a first pile penetration loading on the initial model, update the initial boundary point array, and reconstruct the initial model according to the updated initial boundary point array to obtain a first deformation model;

[0137] A second model processing module 300 is used to assign material parameters, assemble the model and divide the mesh into meshes for the first deformation model to obtain a second boundary point array;

[0138] A second updating module 400 is used to perform a second pile penetration loading on the first deformation model, update the second boundary point array, and reconstruct the first deformation model according to the updated second boundary point array to obtain a second deformation model;

[0139] The loop module 500 is used to loop in this way until a target deformation model with a predetermined deformation degree is obtained.

[0140] It should be noted that the aforementioned explanation of the embodiment of the large deformation finite element model construction method is also applicable to the large deformation finite element model construction device of this embodiment, and will not be repeated here.

[0141] The present invention discloses a large deformation finite element model construction device, which establishes an initial model, assigns material parameters to the initial model, assembles the model and divides the mesh to obtain an initial boundary point array; performs the first pile penetration loading on the initial model, updates the initial boundary point array, and reconstructs the initial model according to the updated initial boundary point array to obtain a first deformation model; performs the first deformation model assignment of material parameters, assembles the model and divides the mesh to obtain a second boundary point array; performs the second pile penetration loading on the first deformation model, updates the second boundary point array, and reconstructs the first deformation model according to the updated second boundary point array to obtain a second deformation model; and repeats this cycle until a target deformation model with a predetermined deformation degree is obtained. The present application records the initial boundary point array when the model is not deformed, and by continuously updating the initial boundary point array, accurately records the changes in the model boundary, and there is no approximate fitting process, thereby improving the accuracy of the reconstructed model.

[0142] Fig.13 A schematic diagram of the structure of a terminal provided in an embodiment of the present application. The terminal may include:

[0143] A memory 501 , a processor 502 , and a computer program stored in the memory 501 and executable on the processor 502 .

[0144] When the processor 502 executes the program, the large deformation finite element model building method provided in the above embodiment is implemented.

[0145] Furthermore, the terminal further includes:

[0146] The communication interface 503 is used for communication between the memory 501 and the processor 502 .

[0147] The memory 501 is used to store computer programs that can be executed on the processor 502 .

[0148] The memory 501 may include a high-speed RAM memory, and may also include a non-volatile memory, such as at least one disk memory.

[0149] If the memory 501, the processor 502 and the communication interface 503 are implemented independently, the communication interface 503, the memory 501 and the processor 502 can be connected to each other through a bus and communicate with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component (PCI) bus or an Extended Industry Standard Architecture (EISA) bus. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, only one line is used in the figure, but it does not mean that there is only one bus or one type of bus.

[0150] Optionally, in a specific implementation, if the memory 501, the processor 502 and the communication interface 503 are integrated on a chip, the memory 501, the processor 502 and the communication interface 503 can communicate with each other through an internal interface.

[0151] The processor 502 may be a central processing unit (CPU), or an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application.

[0152] This embodiment also provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the above-mentioned large deformation finite element model construction method is implemented.

[0153] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or N embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

[0154] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Therefore, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the features. In the description of this application, "N" means at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0155] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or N executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present application includes additional implementations, in which the order shown or discussed may not be followed, including performing functions in a substantially simultaneous manner or in reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present application belong.

[0156] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can read instructions from an instruction execution system, device or apparatus and execute instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or N wirings (electronic devices), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), a fiber optic device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically by optically scanning the paper or other medium and then editing, interpreting or otherwise processing in a suitable manner if necessary and then storing it in a computer memory.

[0157] It should be understood that the various parts of the present application can be implemented by hardware, software, firmware or a combination thereof. In the above embodiment, the N steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. If implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0158] A person skilled in the art may understand that all or part of the steps in the above-mentioned embodiment method may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

[0159] In addition, each functional unit in each embodiment of the present application may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0160] The storage medium mentioned above may be a read-only memory, a magnetic disk or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limiting the present application. A person of ordinary skill in the art may change, modify, replace and modify the above embodiments within the scope of the present application.

Claims

1. A large deformation finite element model construction method, characterized in that: The method comprises: Establishing an initial model, assigning material parameters to the initial model, assembling the model and dividing the mesh, and obtaining an initial boundary point array; Performing a first pile penetration loading on the initial model, updating the initial boundary point array, and reconstructing the initial model according to the updated initial boundary point array to obtain a first deformation model; After assigning material parameters, assembling the model and dividing the mesh on the first deformation model, a second boundary point array is obtained; Performing a second pile penetration loading on the first deformation model, updating the second boundary point array, and reconstructing the first deformation model according to the updated second boundary point array to obtain a second deformation model; This cycle is repeated until a target deformation model with a predetermined deformation degree is obtained; An initial model is established, and after assigning material parameters, model assembly and meshing the initial model, an initial boundary point array is obtained, including: Establishing an initial model, assigning material parameters to the initial model, assembling the model and dividing the mesh, and obtaining all initial mesh nodes; Searching all initial mesh nodes for an initial mesh node on a model boundary of the initial model; forming an initial boundary point array based on all initial grid nodes on the model boundary of the initial model; Searching for an initial mesh node on a model boundary of the initial model among all initial mesh nodes includes: Using each corner point of the initial model as an initial boundary point of the initial model, and using the model boundary formed between each initial boundary point as a target straight line; Screening the points to be searched from all the initial grid nodes, and determining the positional relationship between the points to be searched and the target straight line; If the point to be searched is on the target straight line, determining that the initial grid node is an initial grid node on a model boundary of the initial model; Screening the points to be searched from all the initial grid nodes and determining the positional relationship between the points to be searched and the target straight line includes: Obtaining a first horizontal coordinate, a first vertical coordinate, a second horizontal coordinate, and a second vertical coordinate corresponding to two initial boundary points of the target straight line, and determining a boundary straight line equation corresponding to the target straight line; Filter out points to be searched whose abscissas are between the first abscissa and the second abscissa and whose ordinates are between the first ordinate and the second ordinate from all initial grid nodes; Calculate the distance between the point to be searched and the boundary line equation according to the horizontal coordinate and the vertical coordinate of the point to be searched; If the distance between the point to be searched and the boundary line equation is less than a preset value, the positional relationship between the point to be searched and the target line is that the point to be searched is on the target line; If the distance from the point to be searched to the boundary line equation is greater than or equal to a preset value, the positional relationship between the point to be searched and the target line is that the point to be searched is not on the target line.

2. The large deformation finite element model construction method according to claim 1, characterized in that: Performing a first pile penetration loading on the initial model, updating the initial boundary point array, and reconstructing the initial model according to the updated initial boundary point array to obtain a first deformation model, including: After the initial model is subjected to the first pile penetration loading, a first analysis result is obtained; updating the initial boundary point array according to the first analysis result to obtain an updated initial boundary point array; The initial model is reconstructed according to the updated initial boundary point array to obtain a first deformed model.

3. The large deformation finite element model construction method according to claim 1, characterized in that: After assigning material parameters, assembling the model and meshing the first deformation model, a second boundary point array is obtained, including: After assigning material parameters, assembling the model and meshing the first deformation model, all second mesh nodes are obtained; Searching for a second mesh node on a model boundary of the first deformed model among all second mesh nodes; A second boundary point array is formed based on all second mesh nodes on the model boundary of the first deformation model.

4. The large deformation finite element model construction method according to claim 1, characterized in that: Performing a second pile penetration loading on the first deformation model, updating the second boundary point array, and reconstructing the first deformation model according to the updated second boundary point array to obtain a second deformation model, including: After the first deformation model is subjected to a second pile penetration loading, the first deformation model is deformed, and a second analysis result is obtained; updating the second boundary point array according to the second analysis result to obtain an updated second boundary point array; The first deformation model is reconstructed according to the updated second boundary point array to obtain a second deformation model.

5. A large deformation finite element model construction device, characterized in that: The device comprises: The first model processing module is used to establish an initial model, assign material parameters to the initial model, assemble the model and divide the mesh to obtain an initial boundary point array; A first updating module, configured to perform a first pile penetration loading on the initial model, update the initial boundary point array, and reconstruct the initial model according to the updated initial boundary point array to obtain a first deformation model; A second model processing module, configured to assign material parameters, assemble the model and divide the mesh into meshes for the first deformation model, thereby obtaining a second boundary point array; A second updating module, configured to perform a second pile penetration loading on the first deformation model, update the second boundary point array, and reconstruct the first deformation model according to the updated second boundary point array to obtain a second deformation model; A loop module, used for looping in this way until a target deformation model with a predetermined deformation degree is obtained; An initial model is established, and after assigning material parameters, model assembly and meshing the initial model, an initial boundary point array is obtained, including: Establishing an initial model, assigning material parameters to the initial model, assembling the model and dividing the mesh, and obtaining all initial mesh nodes; Searching all initial mesh nodes for an initial mesh node on a model boundary of the initial model; forming an initial boundary point array based on all initial grid nodes on the model boundary of the initial model; Searching for an initial mesh node on a model boundary of the initial model among all initial mesh nodes includes: Using each corner point of the initial model as an initial boundary point of the initial model, and using the model boundary formed between each initial boundary point as a target straight line; Screening the points to be searched from all the initial grid nodes, and determining the positional relationship between the points to be searched and the target straight line; If the point to be searched is on the target straight line, determining that the initial grid node is an initial grid node on a model boundary of the initial model; Screening the points to be searched from all the initial grid nodes and determining the positional relationship between the points to be searched and the target straight line includes: Obtaining a first horizontal coordinate, a first vertical coordinate, a second horizontal coordinate, and a second vertical coordinate corresponding to two initial boundary points of the target straight line, and determining a boundary straight line equation corresponding to the target straight line; Filter out points to be searched whose abscissas are between the first abscissa and the second abscissa and whose ordinates are between the first ordinate and the second ordinate from all initial grid nodes; Calculate the distance between the point to be searched and the boundary line equation according to the horizontal coordinate and the vertical coordinate of the point to be searched; If the distance between the point to be searched and the boundary line equation is less than a preset value, the positional relationship between the point to be searched and the target line is that the point to be searched is on the target line; If the distance from the point to be searched to the boundary line equation is greater than or equal to a preset value, the positional relationship between the point to be searched and the target line is that the point to be searched is not on the target line.

6. A terminal, characterized in that: include: A memory, a processor, and a large deformation finite element model construction program stored in the memory and executable on the processor, wherein the large deformation finite element model construction program, when executed by the processor, implements the steps of the large deformation finite element model construction method as described in any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and the computer program can be executed to implement the steps of the large deformation finite element model construction method according to any one of claims 1 to 4.

Images