Engineering structure simulation method, device, equipment and storage medium

Abstract

The present application provides an engineering structure simulation method, device, equipment and storage medium, which relate to the field of numerical simulation technology. The method includes: obtaining a first calculation result and a second result by calculation; determining the coupling characteristic parameters of the coupling area according to the first calculation result; adjusting the coupling characteristic parameters according to the second calculation result to obtain the iterated coupling characteristic parameters; then repeatedly adjusting the coupling characteristic parameters to obtain the iterated target coupling characteristic parameters; and regenerating the engineering structure simulation result. When performing simulation, the calculation methods of different simulation areas can be dynamically adjusted according to the calculation results, the coupling characteristic parameters of the contact area can be adjusted in a sequential iterative manner, and the calculation results obtained by different calculation methods can be automatically coupled, which can improve the continuity and coordination of the contact area data of the engineering simulation model, thereby improving the efficiency and accuracy of the engineering structure simulation.

Description

Technical Field

[0001] The present application relates to the field of numerical simulation technologies, and particularly relates to an engineering structure simulation method, device, equipment, and storage medium. Background Art

[0002] There are a wide variety of existing numerical simulation software, such as Ansys, Abaqus, and Flac3D, etc. Different software deal with different engineering problems and calculation methods. Moreover, when an engineering problem cannot be simulated by a single calculation method for an engineering structure model, the simulation results calculated are prone to errors. When using multiple calculation methods to calculate an engineering structure, it is difficult to integrate multiple calculation methods on a single software platform and achieve dynamic coupling. This leads to the numerical simulation of many complex engineering projects having to rely on a single calculation method, resulting in the problem that it is difficult to achieve coupling and the accuracy of the simulation results is relatively low. For some engineering problems, although it is possible to achieve coupled calculation by manually constructing a coupling transfer mechanism for the contact area by using multiple software platforms simultaneously, the pre-processing workload of this method is large, the efficiency is low, and it has quite high requirements for the professional knowledge of users, forming a relatively high theoretical threshold, which restricts its wide application and efficient implementation in actual engineering. Summary of the Invention

[0003] The present application provides an engineering structure simulation method, device, equipment, and storage medium, which can automatically adopt different calculation methods for different simulation areas in an engineering simulation model to perform calculations and obtain corresponding calculation results, thereby improving the efficiency of engineering structure simulation calculations. When performing simulation, the coupling characteristic parameters of the contact area can be adjusted in an iterative manner, and by automatically calculating the coupling characteristic parameters of the calculation results of adjacent areas, dynamic coupling processing of the calculation results of the contact area between the simulation areas using different calculation methods for parameter calculation can be achieved, thereby improving the efficiency and accuracy of engineering structure simulation.

[0004] In a first aspect, an embodiment of the present application provides an engineering structure simulation method, which may include:

[0005] For a first simulation area of an engineering simulation model, use a first simulation calculation method to calculate the simulation parameters in the first simulation area to obtain a first calculation result; wherein, the engineering simulation model includes multiple simulation areas; the first simulation area is any one of the simulation areas; the first calculation result represents the stress state information and displacement state information of the first simulation area;

[0006] Determine the coupling characteristic parameters of the coupling region based on the first calculation result, the position information of the known points in the coupling region, and the position information of the coupling points in the coupling region; wherein, in the calculation method, the coupling region is the contact region between the first simulation region and the second simulation region, and the second simulation region is the simulation region adjacent to the first simulation region; the coupling characteristic parameters characterize the displacement information and force information of the coupling points in the coupling region;

[0007] Calculate the simulation parameters of the second simulation region according to the second simulation calculation method to obtain a second calculation result, and adjust the coupling characteristic parameters according to the second calculation result to obtain the iterated coupling characteristic parameters;

[0008] Take the first simulation region as the second simulation region, and repeat the steps of calculating the simulation parameters of the second simulation region according to the second simulation calculation method to obtain a second calculation result, and adjusting the coupling characteristic parameters according to the second calculation result to obtain the iterated coupling characteristic parameters. When the error between the coupling characteristic parameters before iteration and the coupling characteristic parameters after iteration is less than the preset error threshold, determine the iterated coupling characteristic parameters as the target coupling characteristic parameters;

[0009] Generate an engineering structure simulation result based on one or more of the target coupling characteristic parameters.

[0010] In the above implementation manner, when simulating an engineering structure, different calculation methods can be automatically adopted for different simulation regions in the engineering simulation model to obtain corresponding calculation results. By automatically calculating the coupling characteristic parameters of the calculation results of adjacent regions and adjusting the coupling characteristic parameters of the contact region in an iterative manner, dynamic coupling processing of the calculation results of the contact region between simulation regions using different calculation methods is realized, thereby realizing dynamic adjustment of the simulation calculation methods of different simulation regions, and the efficiency and accuracy of engineering structure simulation calculation can be improved.

[0011] In some embodiments, the coupling characteristic parameters include a first coupling characteristic parameter, and the first coupling characteristic parameter characterizes the displacement information of the coupling points in the coupling region. The determining the coupling characteristic parameters of the coupling points in the coupling region based on the first calculation result, the position information of the known points in the coupling region, and the position information of the coupling points in the coupling region may include:

[0012] Determine the sum of the reciprocals of the distances from the known points in the first coupling region to all the coupling points in the first coupling region; wherein, the known points in the first coupling region are any known points in the coupling region; the coupling points in the first coupling region are any coupling points in the coupling region;

[0013] Determine a first weight coefficient based on the sum of reciprocals, the known points in the first coupling region, and the first distance between the coupling points in the first coupling region; the first weight coefficient characterizes the displacement distribution coefficient transmitted from the known points in the first coupling region to the coupling points in the first coupling region.

[0014] Determine the first coupling characteristic parameter based on the weight coefficient and the displacement information of the known points in the coupling region.

[0015] In the above implementation process, the weight coefficient determined by the distance between the known points in the coupling region and the coupling points in the coupling region characterizes the influence degree of each known point in the coupling region on the displacement of the coupling points in the coupling region. Therefore, the new displacement value of the coupling points in the coupling region can be determined based on the weight coefficient and the displacement information of the known points in the coupling region. By determining the first coupling characteristic parameter characterizing the coupling points in the coupling region in the above manner, the automatic calculation of the first coupling characteristic parameter between different calculation methods can be realized, so as to realize the dynamic coupling processing of the displacement information in the contact region between different simulation regions, and automatically realize the continuity and coordination of the displacements in the contact region between different regions, thereby improving the efficiency and accuracy of the engineering structure simulation calculation.

[0016] In some embodiments, the coupling characteristic parameter includes a second coupling characteristic parameter, and the second coupling characteristic parameter characterizes the force information of the coupling points in the coupling region.

[0017] The determining the coupling characteristic parameter of the coupling points in the coupling region according to the first calculation result, the position information of the known points in the coupling region, and the position information of the coupling points in the coupling region may include:

[0018] Determine the sum of the distances from the known points in the second coupling region to all the coupling points in the second coupling region; wherein, the known points in the second coupling region are any known points in the coupling region, and the coupling points in the second coupling region are any coupling points in the coupling region.

[0019] Determine a first distribution coefficient according to the sum of the distances; wherein, the first distribution coefficient characterizes the proportion of the force distributed by the known points in the second coupling region to the coupling points in the second coupling region in the first direction.

[0020] Determine the second coupling characteristic parameter according to the distribution matrix; wherein, the distribution matrix is obtained according to the distribution coefficients in each direction.

[0021] In the above implementation process, the distribution coefficient determined by the sum of the distances between the known points in the coupling region and the coupling points in the coupling region characterizes the proportion of the force that the known points in the coupling region are distributed to the coupling points in the coupling region in one direction. Thus, the force information of the coupling points in the coupling region can be determined according to the proportion of the force that the known points in the coupling region are distributed to the coupling points in the coupling region in each direction. By determining the second coupling characteristic parameter characterizing the force information of the coupling points in the coupling region in the above manner, the automatic calculation of the second coupling characteristic parameter between different calculation methods and the dynamic coupling processing of the stress information in the contact region between different simulation regions are realized. Thus, the dynamic adjustment of the simulation calculation methods in different simulation regions can be automatically realized, and the continuity and coordination of the stress in the contact region between different regions can be achieved, thereby improving the efficiency and accuracy of the engineering structure simulation calculation.

[0022] In some embodiments, after using the first simulation calculation method to calculate the simulation parameters in the first simulation region of the engineering simulation model and obtaining the first calculation result, the method may further include:

[0023] In the case where the first calculation result obtained by calculating the first simulation region based on the first simulation calculation method cannot converge, it is determined to use the third calculation method to calculate the first simulation region.

[0024] In the above implementation process, the calculation methods of each simulation region can be dynamically adjusted according to the simulation calculation results of each simulation region, and a more suitable calculation method can be determined for each simulation region, thereby improving the calculation convergence and calculation accuracy.

[0025] In some embodiments, before using the first simulation calculation method to calculate the simulation parameters in the first simulation region of the engineering simulation model and obtaining the first calculation result, the method may further include:

[0026] Receive engineering parameters and construct the engineering simulation model based on the engineering parameters;

[0027] Determine the target calculation method as the first simulation calculation method according to the material type of the first simulation region, and initialize the coupling region node coupling characteristic parameters between all simulation regions in the engineering simulation model, so as to calculate the simulation parameters in the first simulation region based on the first simulation calculation method and obtain the first calculation result.

[0028] In the above implementation process, the target calculation method is determined to be the first simulation calculation method according to the material type of the first simulation area. This targeted selection can ensure the matching of the calculation method with the material characteristics and improve the accuracy of the simulation. Selecting a calculation method that matches the material type of the target simulation area can optimize the calculation process and reduce unnecessary calculation overhead. Different calculation methods may have different calculation efficiencies and accuracies when dealing with different types of materials. Therefore, using the appropriate calculation method can improve the calculation efficiency while ensuring the accuracy.

[0029] In some embodiments, generating the engineering structure simulation result according to one or more of the target coupling characteristic parameters may include:

[0030] Determining the displacement information and force information of each simulation area in the engineering simulation model according to one or more of the target coupling characteristic parameters;

[0031] Generating one or more of the stress nephogram, displacement distribution diagram, and parameter table of the engineering simulation model according to the displacement information and the force information.

[0032] In some embodiments, the method may further include: generating a calculation report according to the engineering structure simulation result; the calculation report includes one or more of a model parameter summary, a simulation result summary, and a model stability evaluation result.

[0033] In a second aspect, an embodiment of the present application provides an engineering structure simulation device, which may include:

[0034] A calculation module, configured to calculate the simulation parameters in the first simulation area of the engineering simulation model using the first simulation calculation method to obtain a first calculation result; wherein, the engineering simulation model includes multiple simulation areas; the first simulation area is any one of the simulation areas; the first calculation result represents the stress state information and displacement state information of the first simulation area;

[0035] A first determination module, configured to determine the coupling characteristic parameters of the coupling area according to the first calculation result, the position information of the known points in the coupling area, and the position information of the coupling points in the coupling area; wherein, the coupling area is the contact area between the first simulation area and the second simulation area, and the second simulation area is a simulation area adjacent to the first simulation area; the coupling characteristic parameters represent the displacement information and force information of the coupling points in the coupling area;

[0036] An iterative module, configured to calculate simulation parameters of the second simulation area according to a second simulation calculation method, obtain a second calculation result, and adjust the coupling characteristic parameters according to the second calculation result to obtain the iterated coupling characteristic parameters;

[0037] The first determination module is further configured to use the first simulation area as the second simulation area, and repeatedly execute the steps of calculating the simulation parameters of the second simulation area according to the second simulation calculation method to obtain a second calculation result, and adjusting the coupling characteristic parameters according to the second calculation result to obtain the iterated coupling characteristic parameters. When the error between the coupling characteristic parameters before iteration and the coupling characteristic parameters after iteration is less than a preset error threshold, determine the iterated coupling characteristic parameters as the target coupling characteristic parameters;

[0038] A generation module, configured to generate an engineering structure simulation result according to one or more of the target coupling characteristic parameters.

[0039] In a third aspect, an embodiment of the present application provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the method described in the above description is implemented.

[0040] In a fourth aspect, an embodiment of the present application further provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the method described in the above description is implemented.

[0041] Compared with the prior art, the beneficial effects of the present application are as follows: When simulating an engineering structure, different calculation methods can be automatically adopted for different simulation areas in the engineering simulation model to obtain corresponding calculation results, which can improve the efficiency of engineering structure simulation. When simulating an engineering structure, the calculation methods of different simulation areas in the engineering simulation model can be dynamically adjusted according to the calculation results. By calculating the coupling characteristic parameters of the calculation results of adjacent areas, dynamic coupling processing of the calculation results can be realized, and the continuity and coordination of the contact area data of each simulation area can be automatically achieved, thereby improving the efficiency and accuracy of engineering structure simulation calculation. Description of the Drawings

[0042] Figure 1 It is a schematic diagram of the steps of the engineering structure simulation method provided by an embodiment of the present application.

[0043] Figure 2 It is a schematic diagram of a simulation model with multiple simulation areas provided by an embodiment of the present application.

[0044] Figure 3Schematic diagram of the soil layer geometric model provided by the embodiments of the present application.

[0045] Figure 4 Stress nephogram of the generated soil layer geometric model provided by the embodiments of the present application.

[0046] Figure 5 Schematic diagram of the engineering structure simulation device provided by the embodiments of the present application. Detailed implementation manners

[0047] The present application will be further described in detail below in conjunction with test examples and detailed implementation manners. However, it should not be understood that the scope of the above-mentioned subject matter of the present application is limited to the following embodiments. All technologies implemented based on the content of the present application fall within the scope of protection of the present application.

[0048] In the description of the specific embodiments of the present application, unless otherwise specified, the expression terms indicating the orientation or positional relationship such as "upper", "lower", "left", "right", "center", "inner", "outer", "side", etc. are all based on the orientation or positional relationship shown in the drawings, or the orientation or positional relationship when the product / device / device is usually used. These terms of orientation or positional relationship are only for the convenience of describing the solution of the present application or simplifying the description in the specific embodiments, so as to facilitate technicians to quickly understand the solution, rather than indicating or implying that a specific device / component / element must have a specific orientation, or be constructed and operated in a specific positional relationship. Therefore, it should not be construed as a limitation to the present application.

[0049] In the description of the embodiments of the present application, technical terms such as "first" and "second" only distinguish one entity or operation from another entity or operation, and cannot be understood as indicating or implying relative importance or implicitly indicating the quantity, specific order or primary and secondary relationship of the indicated technical features. In the description of the embodiments of the present application, the meaning of "a plurality" is two or more, unless otherwise specifically defined.

[0050] Referring to "embodiments" herein means that the specific features, structures or characteristics described in connection with the embodiments may be included in at least one embodiment of the present application. The phrase appears in various places in the specification and does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment mutually exclusive with other embodiments. Those skilled in the art explicitly and implicitly understand that the embodiments described herein may be combined with other embodiments.

[0051] Embodiment 1

[0052] During the research process, the applicant found that when using simulation software such as Ansys, Abaqus, and Flac3D to perform simulation calculations on engineering structures, if different calculation methods are required to calculate the engineering structure model, cumbersome steps such as geometric modeling and mesh generation are needed. Only the calculation method to be used for each region can be manually specified before the calculation, and only the calculation results of two fixed calculation methods can be coupled on multiple software platforms. When solving practical engineering problems, for the accuracy of simulation calculations, an engineering model may need to use more calculation methods to calculate different regions. At this time, the above-mentioned software cannot achieve the dynamic coupling of multiple calculation methods. Therefore, after researching this problem, the applicant proposed an engineering structure simulation method. When simulating an engineering structure, by automatically calculating the coupling characteristic parameters of the calculation results of adjacent regions, the dynamic coupling processing of the calculation results of the contact regions between the simulation regions using different calculation methods for parameter calculation can be realized, thereby improving the accuracy of engineering structure simulation.

[0053] Please refer to Figure 1 , Figure 1 which is a schematic diagram of the steps of the engineering structure simulation method provided by the embodiment of the present application. The engineering structure simulation method may include:

[0054] S1. For the first simulation region of the engineering simulation model, use the first simulation calculation method to calculate the simulation parameters in the first simulation region to obtain the first calculation result.

[0055] Among them, the engineering simulation model includes multiple simulation regions; the first simulation region is any simulation region; the first calculation result represents the stress state information and displacement state information of the first simulation region.

[0056] S2. Determine the coupling characteristic parameters of the coupling region according to the first calculation result, the position information of the known points in the coupling region, and the position information of the coupling points in the coupling region.

[0057] Among them, the coupling region is the contact region between the first simulation region and the second simulation region, and the second simulation region is the simulation region adjacent to the first simulation region; the coupling characteristic parameters represent the displacement information and force information of the coupling points in the coupling region.

[0058] S3. Calculate the simulation parameters of the second simulation region according to the second simulation calculation method to obtain the second calculation result, and adjust the coupling characteristic parameters according to the second calculation result to obtain the iterated coupling characteristic parameters.

[0059] S4. Take the first simulation area as the second simulation area, and repeatedly execute the steps of calculating the simulation parameters of the second simulation area according to the second simulation calculation method to obtain a second calculation result, and adjusting the coupling characteristic parameters according to the second calculation result to obtain the iterated coupling characteristic parameters. When the error between the coupling characteristic parameters before iteration and the iterated coupling characteristic parameters is less than the preset error threshold, determine the iterated coupling characteristic parameters as the target coupling characteristic parameters.

[0060] S5. Generate an engineering structure simulation result according to one or more target coupling characteristic parameters.

[0061] The engineering structure simulation method provided in the embodiments of the present application can be applied to many technical fields, such as the structural design and optimization of aircraft, rockets and other aircraft, including stress analysis and fatigue life prediction; automotive crash safety analysis, simulating the behavior of vehicles under different collision conditions, and body structure strength, stiffness and durability analysis; stability and safety assessment of building structures and geotechnical engineering; strength, fatigue and durability simulation analysis of mechanical parts, etc.

[0062] The engineering simulation model can be constructed by the user inputting simulation values in the relevant interface. Please refer to Figure 2 , Figure 2 which is a schematic diagram of a simulation model with multiple simulation areas provided in the embodiments of the present application. The engineering simulation model may include multiple first simulation areas composed of different materials, such as simulation area 1, simulation area 2, simulation area 3, simulation area 4, simulation area 5, and simulation area 6. Since each simulation area is composed of different materials, different calculation methods can be used for calculation. For example, simulation areas 1 and 4 are calculated using calculation method A, simulation areas 2 and 3 are calculated using calculation method B, and simulation areas 5 and 6 are calculated using calculation method C. The first calculation method and the second calculation method are different calculation methods, so it is necessary to couple the calculation results of adjacent simulation areas using different calculation methods.

[0063] In the physical world, physical quantities between different regions are continuous, such as displacement, stress, etc. Different calculation methods may use different mathematical models and numerical algorithms, resulting in numerical differences in the calculation results. Therefore, when performing simulation calculations using different calculation methods, it is necessary to couple the calculation results between the regions using different calculation methods. By creating contacts between the regions to establish a channel for transferring calculation data between different regions, in the contact area, it is necessary to define data transfer rules and conversion mechanisms, which include determining which calculation parameters need to be transferred between different regions, and how to convert the results calculated in one region into a form that can be understood and used by another region, so as to make the calculation results between different regions numerically consistent.

[0064] When performing coupling processing, the displacement information and force information of the nodes in the coupling region are characterized by determining the coupling characteristic parameters of the coupling region. Then, by iterating the coupling characteristic parameters, the target coupling characteristic parameters are obtained when the iteration is completed, so as to accurately characterize the displacement information and force information in the contact region of each adjacent simulation region. Thus, the engineering structure simulation results can be generated according to the target coupling characteristic parameters, and the continuity and coordination between the contacts of different simulation regions can be improved. The specific implementation manner of the coupling processing can refer to the description content in the following embodiments and will not be elaborated here for the time being.

[0065] The process of iterating the coupling characteristic parameters is specifically to correct the parameters related to force, displacement, etc. in the coupling characteristic parameters according to the second calculation result to reflect the influence of the second simulation region on the coupling region. Then, the adjusted coupling region characteristic parameters are input into the calculation model of the first simulation region, so that the influence of the second simulation region is transmitted to the first simulation region through the adjustment of the coupling region characteristic parameters. Then, based on the adjusted coupling region characteristic parameters, the model parameters of the first simulation region, and the first calculation method, the first simulation region is simulated and calculated. By repeating step S3, the coupling calculation is performed based on the sequential iteration method until the error between two iterations is less than the preset error threshold, and then the iteration is stopped to obtain the target coupling characteristic parameters that accurately characterize the displacement information and force information in the contact region of each adjacent simulation region.

[0066] If the engineering simulation model has only two simulation regions, the engineering structure simulation results can be generated according to one target coupling characteristic parameter. If the engineering simulation model has more simulation regions, the target coupling characteristic parameters of each adjacent region are determined respectively, and then the engineering structure simulation results are generated according to all the target coupling characteristic parameters.

[0067] In the above implementation manner, when simulating the engineering structure, different calculation methods can be automatically adopted for different simulation regions in the engineering simulation model to obtain the corresponding calculation results, which can improve the efficiency of the engineering simulation calculation. When simulating the engineering structure, the calculation methods of different simulation regions in the engineering simulation model can be dynamically adjusted according to the calculation results. By automatically calculating the coupling characteristic parameters of the calculation results of adjacent regions, the dynamic coupling processing of the calculation results in the contact region between the simulation regions using different calculation methods for simulation calculation is realized, so that the efficiency and accuracy of the engineering structure simulation calculation can be improved.

[0068] Embodiment 2

[0069] This embodiment is a specific implementation scheme for determining the coupling characteristic parameters of the coupling region according to the first calculation result and the second calculation result in the above Embodiment 1.

[0070] In the embodiments of the present application, the coupling characteristic parameter includes a first coupling characteristic parameter, and the first coupling characteristic parameter characterizes the displacement information of the coupling points in the coupling region.

[0071] There are various calculation methods in the embodiments of the present application, such as the finite element method, the material point method, the finite difference method, etc. Since the data transfer methods are different under different calculation methods, in order to improve the continuity and coordination of the data in the contact area between different regions, it is necessary to determine the data transfer between different regions, so as to couple the calculation results of different regions. Among them, the data transfer may include displacement transfer, and the new displacement value of the coupling region on the coupling interface is determined based on the simulation calculation results.

[0072] The embodiments of the present application are described by taking the coupling of the finite element-material point method as an example. When other calculation methods are used for coupling, this method can also be referred to for data coupling. Among them, the coupling of the finite element-material point method is to discretize small deformation objects by the finite element method and large deformation objects by the material point method, and couple them with each other through a contact algorithm between different discrete regions.

[0073] The steps of determining the coupling characteristic parameter of the coupling region according to the first calculation result and the second calculation result may include:

[0074] S21. Determine the sum of the reciprocals of the distances from the known points in the first coupling region to all the coupling points in the first coupling region.

[0075] Among them, the known point in the first coupling region is any known point in the coupling region; the coupling point in the first coupling region is any coupling point in the coupling region.

[0076] S22. Determine the first weight coefficient according to the sum of the reciprocals, the known point in the first coupling region, and the first distance between the known point in the first coupling region and the coupling point in the first coupling region.

[0077] Among them, the first weight coefficient characterizes the influence degree of the known point in the first coupling region on the displacement of the coupling point in the first coupling region. The known point in the first coupling region is any known point in the coupling region; the coupling point in the first coupling region is any coupling point in the coupling region; the first weight coefficient characterizes the displacement distribution coefficient transmitted from the known point in the first coupling region to the coupling point in the first coupling region.

[0078] In the embodiments calculated by using the finite element-material point method, in the first iteration, the material point method region is calculated first, and then the finite element method region is calculated. The known points in the coupling region are the material points, and the coupling points in the coupling region are the finite element points.

[0079] At the coupling interface, to make the displacement continuous, it is necessary to transfer the displacement information of the material point method region to the finite element method region. By setting shared nodes or interpolation points on the coupling interface, for the displacement of the material points near the coupling interface in the material point method region, it is distributed to the nodes in the finite element method region by interpolation. For example, assuming there is a material point i and a finite element point j on the coupling interface, according to the distance between them, the weight coefficient w can be established ij :

[0080]

[0081] where r ij is the distance from the material point i to the finite element point j, n is any finite element point (taking values 1, …, k), k is the number of all finite element points in the coupling region, r in represents the distance from the material point i to the finite element point n, represents the sum of the reciprocals of the distances from the material point i to all finite element points. When calculating the weight coefficient w ij , the reciprocals of the distances r in from the material point to all finite element points are summed, and the ratio of the reciprocal of the distance between the material point and the finite element point to the total sum is used as the allocation weight w ij of the displacement of the material point i to the finite element point j. In this way, the displacement of the finite element point can be updated according to the displacement of the material point. Based on the displacement allocation weights of all material points to the finite element point j, the displacement of the finite element point can be obtained. Among them, is the displacement of the material point i.

[0082] S23. Determine the first coupling characteristic parameter according to the weight coefficient and the displacement information of the known points in the coupling region.

[0083] On the coupling interface, the displacement should be continuous. If the displacement in the material point method region is u MPM , and the displacement in the finite element method region is u FEM , that is, u MPM = u FEM . If the displacement in the material point method region is obtained through the movement of the material points, for the material point i near the coupling interface, its displacement , where x i is the current position of the material point i, is the initial position of the material point i.

[0084] In the finite element method region, the displacement u FEM of the finite element method region can be obtained by interpolating the displacements of the finite element points, that is:

[0085]

[0086] where Nj is the element shape function, which is defined based on the geometric shape of the element and the node positions, and can describe the distribution of the displacement within the element. The element is the basic geometric body or element that constitutes the entire finite element method region; ε, μ, δ are the natural coordinates within the element, d j is the first coupling characteristic parameter, that is, the displacement of the finite element point.

[0087] In the above implementation process, the weight coefficient determined by the distance between the known points in the coupling region and the coupling points in the coupling region characterizes the influence degree of each known point in the coupling region on the displacement of the coupling points in the coupling region. Thus, the new displacement value of the coupling points in the finite element method region where the coupling points in the coupling region are located can be determined according to the weight coefficient and the displacement information of the known points in the coupling region. By determining the first coupling characteristic parameter representing the displacement of the finite element method region in the above manner, the continuity and coordination of the data in the contact region between different regions can be automatically realized, thereby improving the efficiency and accuracy of the simulation calculation of the engineering structure.

[0088] Since the data transfer methods are different under different calculation methods, in order to ensure the continuity and harmony of the data in the contact region between different regions, it is necessary to determine the data transfer mechanism between different regions, so as to couple the calculation results of different regions. Among them, the data transfer can also include the transfer of force, that is, to determine the force information of the coupling points in the coupling region on the coupling interface. In the embodiments of the present application, the coupling of the finite element - material point method is continued as an example for illustration, and a simulation method for coupling the forces of the calculation results of different regions of the engineering structure model is provided.

[0089] The step of determining the coupling characteristic parameter of the coupling region according to the first calculation result and the second calculation result may further include:

[0090] S24. Determine the sum of the distances from the known points in the second coupling region to all the coupling points in the second coupling region.

[0091] Among them, the known points in the second coupling region are any known points in the coupling region, and the coupling points in the second coupling region are any coupling points in the coupling region.

[0092] The transfer of force is the key to ensuring the conservation of momentum. When the finite element method region deforms, the force generated by it needs to be transferred to the material point method region, and vice versa. On the coupling interface, for the force in the material point method region, it needs to be distributed to the nodes in the finite element method region in an appropriate manner. According to the principle of force equivalence, the concentrated force in the material point method region can be equivalent to the equivalent nodal force on the nodes in the finite element method region.

[0093] There are n material points and m finite element points on the coupling interface. Let the force received by material point i be:

[0094]

[0095] The equivalent nodal force of the finite element point j is as follows:

[0096]

[0097] It can be assumed that there exists a distribution matrix T, whose element t ij represents the proportionality coefficient of the force of material point i distributed to node j.

[0098] S25. Determine the first distribution coefficient according to the sum of distances.

[0099] Among them, the first distribution coefficient characterizes the proportion of the force of the known points in the second coupling region distributed to the coupling points in the second coupling region in the first direction; the known points in the second coupling region are any known points in the coupling region, and the coupling points in the second coupling region are any coupling points in the coupling region.

[0100] The first direction can be the x direction of the space coordinate system. The following embodiments will be described by taking the first distribution coefficient t ijx characterizing the proportion of the force distributed to the current finite element point in the x direction of the distribution matrix T as an example:

[0101] On the coupling interface, according to Newton's third law, the force F exerted by the material point method region on the finite element method region MPM-FEM and the force F exerted by the finite element method region on the material point method region FEM-MPM are equal in magnitude and opposite in direction, that is:

[0102]

[0103] In the material point method, the force can be calculated through the stress tensor σ MPM and the area A of the material point MPM , that is:

[0104]

[0105] In the finite element method, the force is calculated through the element stiffness matrix K FEM and the nodal displacement d FEM , that is:

[0106]

[0107] The transfer of force is the key to ensuring the conservation of momentum. When the finite element method region deforms, the force generated needs to be transferred to the material point method region, and vice versa.

[0108] Therefore, on the coupling interface, for the force of the material point method region, it can be distributed to the nodes of the finite element method region in an appropriate manner. According to the principle of force equivalence, the concentrated force of the material point method region can be equivalent to the equivalent nodal force on the nodes of the finite element method region.

[0109] There are n material points and m finite element points on the coupling interface. Let the force acting on material point i be:

[0110]

[0111] The equivalent nodal force of finite element point j is:

[0112]

[0113] Assume there exists a distribution matrix T, whose element t ij represents the proportion coefficient of the force of material point i distributed to finite element point j.

[0114] According to the principle of force equivalence, the equivalent nodal force of finite element point j is the sum of the forces of all material points distributed to this node. The equivalent nodal force of finite element point j in the x direction , where, t ijx is the first distribution coefficient of the distribution matrix T in the x direction, representing the proportion of the x-direction force of material point i distributed to node j, that is, the proportion of the forces of all second material points in the first direction distributed to the current finite element point; F ix is the force acting on material point i in the x direction.

[0115] S26. Determine the second coupling characteristic parameter according to the distribution matrix.

[0116] Among them, the distribution matrix is obtained according to the distribution coefficients in each direction. The distribution matrix T can be determined by the first distribution coefficient in the x direction of the space coordinate system, the second distribution coefficient t ijy in the y direction, and the third distribution coefficient t ijz in the z direction.

[0117] The following takes constructing the distribution matrix based on geometric relationships as an example for illustration. The first distribution coefficient can be determined according to the distance between the material point and the finite element point in the first direction. Let the vector from material point i to node j be , where, r ijx , r ijy and r ijz are the distances between the material point and the finite element point in the first direction, the second direction, and the third direction respectively. In the embodiments of the present application, the first direction is the x direction of the space coordinate system, the second direction is the y direction of the space coordinate system, and the third direction is the z direction of the space coordinate system. The modulus of is:

[0118]

[0119] By summing the distances from material point i to all finite element points j Then the first distribution coefficient can be determined. Similarly, the second distribution coefficient can be determined respectively. And the third distribution coefficient After determining the distribution coefficient in each direction, the matrix composed of t ijx , t ijy and t ijz can be determined as the distribution matrix T.

[0120] At this time, the equivalent nodal force of the finite element point j in the x direction , the equivalent nodal force of the finite element point j in the y direction, , and the equivalent nodal force of the finite element point j in the z direction can be determined.

[0121] Finally, the second coupling characteristic parameter can be determined, that is, the equivalent nodal force of the finite element point j in the coupling region .

[0122] In the above implementation process, the distribution coefficient determined by the sum of the distances between the known points in the coupling region and the coupling points in the coupling region represents the proportion of the force distributed by the known points in the coupling region to the coupling points in the coupling region in one direction. Therefore, the force information of the coupling points in the coupling region can be determined according to the proportion of the force distributed by the known points in the coupling region to the coupling points in the coupling region in each direction. By determining the second coupling characteristic parameter representing the force information of the coupling points in the coupling region in the above manner, the continuity and coordination of data transfer in the contact region between different regions can be automatically realized, thereby improving the efficiency and accuracy of engineering structure simulation calculation.

[0123] Embodiment 3

[0124] This embodiment is an implementation manner of dynamically adjusting the calculation methods of different simulation regions according to the simulation results in Embodiment 1 above.

[0125] After using the first simulation calculation method to calculate the simulation parameters in the first simulation region of the engineering simulation model and obtaining the first calculation result, the method may further include:

[0126] In the case where the first calculation result obtained by calculating the first simulation region based on the first simulation calculation method cannot converge, it is determined to use the third calculation method to calculate the first simulation region.

[0127] During the calculation process, the system can automatically modify the calculation methods for different regions according to the forces and deformations obtained during the simulation calculation. For example, if the engineering simulation model includes two simulation regions, A and B, and the finite element method is used for region A during the calculation, but the system finds that the calculation cannot converge during the calculation, it indicates that the stress or deformation in the first simulation region exceeds the set tolerance range. The system can determine to convert the first calculation method to a third calculation method more suitable for this region based on the engineering parameters of this region. For example, the calculation method for region A can be adjusted from the finite element method to the material point method, and then the entire model can be recalculated.

[0128] In the above implementation process, the calculation methods for each simulation region can be dynamically adjusted according to the simulation calculation results of each simulation region, and a more suitable calculation method can be determined for each simulation region, thereby improving the calculation convergence and calculation accuracy.

[0129] Embodiment 4

[0130] This embodiment is a specific implementation scheme for receiving engineering simulation parameters and constructing an engineering simulation model in the above Embodiment 1.

[0131] Before calculating the simulation parameters in the first simulation region using the first simulation calculation method to obtain the first calculation result, the method provided in the embodiments of the present application may further include:

[0132] Receiving engineering parameters and constructing an engineering simulation model based on the engineering parameters;

[0133] Determining that the target calculation method is the first simulation calculation method according to the material type of the first simulation region, and initializing the coupling region node coupling characteristic parameters between all simulation regions in the engineering simulation model, so as to calculate the simulation parameters in the first simulation region based on the first simulation calculation method to obtain the first calculation result.

[0134] In the embodiments of the present application, the engineering structure simulation method can be applied to an engineering structure simulation system. The engineering structure simulation system can receive the engineering parameters input by the user through the web page, and automatically create an engineering simulation model, automatically select a calculation engine, automatically divide the grid, and automatically perform calculations based on the selected calculation engine.

[0135] Among them, the engineering structure simulation system includes multiple calculation engines, and each calculation engine can support one or more calculation methods, such as finite element, finite difference, discrete element, material point and other calculation methods. Different calculation engines can extract simulation parameters from the database through a data interface, and automatically create a geometric model, allocate material parameters, add boundary conditions, and divide calculation units according to the material region; each calculation engine can calculate different regions of the model simultaneously, and can use a distributed parallel calculation method to accelerate the simulation speed.

[0136] During the calculation process, the system can automatically select the calculation engine for each region according to the calculation method corresponding to the material, and carry out numerical simulation calculations on the same simulation model based on different calculation engines. When the calculation method is not included in the material, the pre-set unified calculation method can be used for calculation. If the calculation results of a certain material region show that the stress or deformation exceeds the set tolerance range, and the system has been pre-set to allow the transformation of the calculation method, the system will automatically convert the calculation method of this region to a more suitable calculation method for this region.

[0137] The engineering parameters input by the user can involve aspects such as material parameter creation, geometric model parameter input, mesh generation, boundary condition entry, and calculation method setting. Exemplarily, the user can enter the geometric model and material parameters of the project on the web side, and complete the setting of boundary conditions, mesh refinement, and calculation methods through the drag-and-drop and selection interfaces.

[0138] Among them, in terms of material parameter creation, the user can set parameters for different materials, such as elastic modulus, Poisson's ratio, etc., and specify the corresponding numerical calculation methods (such as finite element method, finite difference method, etc.). In terms of geometric model parameter input, the system can support specifying different calculation methods for different material regions and allow the user to set the conversion rules of the calculation methods. In terms of mesh generation, the system can automatically generate meshes based on the engineering geometric model parameters input by the user on the web side interface and receive data such as the size, shape, and mesh accuracy of the geometric model defined by the user according to the actual engineering needs. In terms of boundary condition entry, the system can receive the boundary conditions specified by the user for different regions of the simulation model, such as concentrated loads, distributed loads, fixed boundaries, free boundaries, etc., to ensure that the calculation boundaries of the model meet the actual engineering requirements. In terms of calculation method setting, the system can provide the user with options to select the initial calculation method for the overall model, the unified calculation method for simulation calculation, and receive the setting of whether to allow the calculation method to change dynamically during the simulation process to improve the simulation accuracy.

[0139] The system can realize the real-time visual display of the model through WebGL technology. The user can directly interact with the model through operations such as dragging and zooming, and modify the input simulation parameters. All parameters are automatically saved in the database on the server side after confirmation.

[0140] In the above implementation process, the target calculation method is determined to be the first simulation calculation method according to the material type of the first simulation area. This targeted selection can ensure the matching of the calculation method with the material characteristics and improve the accuracy of the simulation. Selecting a calculation method that matches the material type of the target simulation area can optimize the calculation process and reduce unnecessary calculation overhead. Different calculation methods may have different calculation efficiencies and accuracies when dealing with different types of materials. Therefore, using an appropriate calculation method can improve the calculation efficiency while ensuring the accuracy.

[0141] Embodiment 5

[0142] This embodiment is a specific implementation of determining the simulation result of an engineering structure according to one or more coupling characteristic parameters in Embodiment 1.

[0143] Determining the simulation result of an engineering structure according to one or more coupling characteristic parameters may include:

[0144] Determining the displacement information and force information of each simulation area in the engineering simulation model according to one or more target coupling characteristic parameters; and generating one or more of a stress nephogram, a displacement distribution diagram, and a parameter table of the engineering simulation model according to the displacement information and the force information.

[0145] In the embodiment of the present application, after the simulation calculation is completed, the system draws the target simulation model through WebGL technology according to the calculated coupling characteristic parameters and displays the calculation result in the form of a three-dimensional nephogram on the web page.

[0146] In some embodiments, the engineering simulation result includes a stress nephogram, a displacement distribution diagram, a parameter table, etc. Among them, in the stress nephogram, the stress distribution of each area can be drawn on the model according to the unit stress value in a predefined color interval. The displacement value of each unit in the model is also presented in a color interval in the displacement distribution diagram, and the user can view the overall deformation of the model through the interface. In the parameter table, the system can display the detailed numerical values of the displacement, stress, strain, etc. of the unit in the form of a list, and the user can compare the specific calculation results of each calculation unit to assist in analyzing the simulation effect.

[0147] In other embodiments, the engineering structure simulation method may further include: generating a calculation report according to the engineering structure simulation result; the calculation report includes one or more of a model parameter summary, a simulation result summary, and a model stability evaluation result.

[0148] Among them, the model parameter summary can summarize the geometric model, material parameters, mesh division, and boundary conditions input by the user. The simulation result summary can include the stress, strain, and displacement data of each calculation unit, and associate these data with the geometric positions in the model to help the user understand the calculation results of specific regions. The model stability evaluation result can evaluate the overall stability of the model and can include the analysis results of key data such as the plastic zone, stress concentration zone, and maximum displacement.

[0149] After obtaining the simulation results, the user can also select different simulation schemes. By comparing the key results (such as stress, deformation, stability analysis, etc.) in different simulation calculation reports through the system, it helps engineering and technical personnel quickly find the optimized design scheme. It can also further automatically optimize and adjust the model according to the comparative analysis results, such as strengthening the mesh refinement or adjusting the material properties in the regions with large deformations, so as to optimize the simulation effect.

[0150] Embodiment 6

[0151] The following is an embodiment provided by this application for analyzing the stability of a retaining wall for a subgrade slope.

[0152] When analyzing the stability of a retaining wall for a subgrade slope, the material parameters of the model can be input from the web page first. The material parameters can specifically include material type, material name, failure criterion, unit weight, elastic modulus, Poisson's ratio, etc. In some embodiments, multiple model material parameters can be pre-stored in the system material library and can be directly called during use, or new model material parameters can be defined during use.

[0153] After defining the material parameters, the geometric model parameters can be entered. The toe vertex of the retaining wall is the base point of the model, that is, the coordinate origin for calculation. Please refer to Figure 3 , Figure 3 which is a schematic diagram of the soil layer geometric model provided by this application embodiment. This application embodiment is a two-dimensional analysis, and the soil layer geometric model is a polyline enclosed by multiple vertex coordinates. The soil layer geometric model contains multiple material regions such as backfill soil, silty clay, and retaining wall. The force condition of the soil layer geometric model is shown by the arrows in the figure. Subsequently, the boundary condition parameters are entered. The distributed load parameters include the starting position of 8.1 m, the ending position of 14.1 m, the starting pressure of 150 kPa, the ending pressure of 150 kPa, the angle type of radial, and the angle of 90°.

[0154] Subsequently, calculations are performed on each region and simulation analysis is carried out. The calculation method can refer to the description content in the above embodiments and will not be elaborated here. After the calculation is completed, a stress nephogram of the soil layer geometric model can be generated according to the calculation results. Please refer to Figure 4 , Figure 4It is the stress nephogram of the generated soil layer geometric model provided by the embodiment of the present application. By calculating each of the above material regions, the stress distribution in the soil layer geometric model can be obtained, mainly including five stress levels: stress level 1, stress level 2, stress level 3, stress level 4, and stress level 5. Each stress level represents a stress range respectively. Among them, the point with the maximum stress is the MAX point, and the point with the minimum stress is the MIN point.

[0155] Embodiment 7

[0156] Please refer to Figure 5 , Figure 5 It is the schematic diagram of the engineering structure simulation device provided by the embodiment of the present application. The engineering structure simulation device 50 may include:

[0157] A calculation module 51, configured to calculate the simulation parameters in the first simulation area of the engineering simulation model using the first simulation calculation method for the first simulation area of the engineering simulation model, and obtain a first calculation result; wherein, the engineering simulation model includes multiple simulation areas; the first simulation area is any simulation area; the first calculation result characterizes the stress state information and displacement state information calculation method of the first simulation area.

[0158] A first determination module 52, configured to determine the coupling characteristic parameters of the coupling points in the coupling area according to the first calculation result, the position information of the known points in the coupling area, and the position information of the coupling points in the coupling area; the coupling area is the contact area between the first simulation area and the second simulation area, and the second simulation area is the simulation area adjacent to the first simulation area; the coupling characteristic parameters characterize the displacement information and force information of the coupling points in the coupling area.

[0159] An iteration module 53, configured to determine the second simulation calculation method for the second simulation area according to the coupling characteristic parameters, and calculate the simulation parameters of the second simulation area through the second simulation calculation method to obtain the iterated coupling characteristic parameters;

[0160] The first determination module 52 may also be configured to use the first simulation area as the second simulation area, repeat the steps of determining the second simulation calculation method for the second simulation area according to the coupling characteristic parameters and calculating the simulation parameters of the second simulation area through the second simulation calculation method, and determine the current coupling characteristic parameters as the target coupling characteristic parameters when the error of the iteration result is less than the preset error threshold;

[0161] A generation module 54, configured to generate an engineering structure simulation result according to one or more target coupling characteristic parameters.

[0162] Optionally, the first determination module 52 may specifically be configured to:

[0163] Determine the sum of the reciprocals of the distances from the known points in the first coupling area to all the coupling points in the coupling area;

[0164] Determine the first weight coefficient according to the sum of reciprocals, the known points in the first coupling region, and the first distance from the coupling points in the first coupling region; wherein, the known points in the first coupling region are any known points in the coupling region; the coupling points in the first coupling region are any coupling points in the coupling region; the first weight coefficient characterizes the displacement distribution coefficient transmitted from the known points in the first coupling region to the coupling points in the first coupling region.

[0165] Determine the first coupling characteristic parameter according to the weight coefficient and the displacement information of the known points in the coupling region.

[0166] Optionally, the first determination module 52 may be specifically configured to:

[0167] Determine the sum of the distances from the known points in the second coupling region to all the coupling points in the coupling region; the known points in the second coupling region are any known points in the coupling region.

[0168] Determine the first distribution coefficient according to the sum of the distances; the first distribution coefficient characterizes the proportion of the force distributed to the coupling points in the second coupling region by the known points in the second coupling region in the first direction; the known points in the second coupling region are any known points in the coupling region, and the coupling points in the second coupling region are any coupling points in the coupling region.

[0169] Determine the second coupling characteristic parameter according to the distribution matrix; wherein, the distribution matrix is obtained according to the distribution coefficients in each direction.

[0170] Optionally, the first determination module 52 may also be used to:

[0171] In the case where the first calculation result obtained by calculating the first simulation region based on the first simulation calculation method cannot converge, determine to use the third calculation method to calculate the first simulation region.

[0172] Optionally, the engineering structure simulation device 50 may further include a receiving module, configured to receive engineering parameters and construct an engineering simulation model based on the engineering parameters.

[0173] Determine that the target calculation method is the first simulation calculation method according to the material type of the first simulation region, and initialize the coupling region node coupling characteristic parameters between all the simulation regions in the engineering simulation model, so as to calculate the simulation parameters in the first simulation region based on the first simulation calculation method to obtain the first calculation result.

[0174] Optionally, the generation module 54 may be specifically configured to:

[0175] Determine the displacement information and force information of each simulation region in the engineering simulation model according to one or more target coupling characteristic parameters.

[0176] Generate one or more of the stress nephogram, displacement distribution diagram, and parameter table of the engineering simulation model based on the displacement information and force information.

[0177] It should be understood that when each module of the engineering structure simulation 30 provided in the above embodiments performs engineering structure simulation calculations, only the division of each functional module in the above description content is used as an example. In actual applications, the above functions can be allocated to different functional modules according to needs, that is, the internal structure of the device can be divided into different functional modules to complete all or part of the functions described above.

[0178] Each functional module in the above embodiments can be integrated in a processing unit, or each unit can exist physically alone, or two or more units can be integrated in one unit. The above integrated unit can be implemented in the form of hardware or in the form of a software functional unit. In addition, the specific names of each functional unit and module are only for the convenience of mutual distinction and do not limit the protection scope of the embodiments of the present application.

[0179] Based on the same application concept, the embodiments of the present application also provide a computer device, which may include a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the engineering structure simulation method as described in the above description content.

[0180] Based on the same application concept, the embodiments of the present application also provide a computer-readable storage medium. The computer-readable storage medium stores a computer program, and when the computer program is executed by the processor, it implements the engineering structure simulation method as described in the above description content.

[0181] The above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them; although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements on some of the technical features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present application, and should all be included in the protection scope of the present application.

Claims

1. A method for simulating an engineering structure, characterized in that: include: For a first simulation area of ​​the engineering simulation model, a first simulation calculation method is used to calculate simulation parameters in the first simulation area to obtain a first calculation result; wherein the engineering simulation model includes multiple simulation areas; the first simulation area is any of the simulation areas; the first calculation result represents stress state information and displacement state information of the first simulation area; Determine the coupling characteristic parameters of the coupling region according to the first calculation result, the position information of the known points in the coupling region, and the position information of the coupling points in the coupling region; wherein the coupling region is a contact region between the first simulation region and the second simulation region, and the second simulation region is a simulation region adjacent to the first simulation region; and the coupling characteristic parameters represent the displacement information and force information of the coupling points in the coupling region; Calculating the simulation parameters of the second simulation area according to a second simulation calculation method to obtain a second calculation result, and adjusting the coupling characteristic parameters according to the second calculation result to obtain the iterated coupling characteristic parameters; The first simulation area is used as the second simulation area, and the steps of calculating the simulation parameters of the second simulation area according to the second simulation calculation method to obtain a second calculation result, and adjusting the coupling characteristic parameter according to the second calculation result to obtain the coupling characteristic parameter after iteration are repeated, and when the error between the coupling characteristic parameter before iteration and the coupling characteristic parameter after iteration is less than a preset error threshold, determining the coupling characteristic parameter after iteration as the target coupling characteristic parameter; An engineering structure simulation result is generated according to one or more of the target coupling characteristic parameters.

2. The method according to claim 1, characterized in that The coupling characteristic parameter includes a first coupling characteristic parameter, and the first coupling characteristic parameter represents the displacement information of the coupling point of the coupling region; The determining the coupling characteristic parameters of the coupling region according to the first calculation result, the position information of the known points in the coupling region, and the position information of the coupling points in the coupling region includes: Determine the sum of the reciprocals of the distances from the first coupling region known point to all first coupling region coupling points; wherein the first coupling region known point is any coupling region known point; and the first coupling region coupling point is any coupling region coupling point; Determine a first weight coefficient according to the sum of the reciprocals, a first distance between the first coupling region known point and the first coupling region coupling point; the first weight coefficient represents a displacement distribution coefficient transmitted from the first coupling region known point to the first coupling region coupling point; The first coupling characteristic parameter is determined according to the weight coefficient and the displacement information of the known point in the coupling area.

3. The method according to claim 1, characterized in that The coupling characteristic parameter includes a second coupling characteristic parameter, and the second coupling characteristic parameter represents the force information of the coupling point of the coupling region; The determining the coupling characteristic parameters of the coupling region according to the first calculation result, the position information of the known points in the coupling region, and the position information of the coupling points in the coupling region includes: Determine the sum of the distances from the second coupling region known point to all second coupling region coupling points; wherein the second coupling region known point is any coupling region known point, and the second coupling region coupling point is any coupling region coupling point; Determine a first distribution coefficient according to the sum of the distances; wherein the first distribution coefficient represents the proportion of the force distributed from the known point of the second coupling region to the coupling point of the second coupling region in the first direction; The second coupling characteristic parameter is determined according to a distribution matrix; wherein the distribution matrix is ​​obtained according to a distribution coefficient in each direction.

4. The method according to claim 1, characterized in that: In the first simulation area of ​​the engineering simulation model, after using a first simulation calculation method to calculate the simulation parameters in the first simulation area and obtaining a first calculation result, the method further includes: When the first simulation region is calculated based on the first simulation calculation method and the first calculation result obtained cannot converge, it is determined to use a third calculation method to calculate the first simulation region.

5. The method according to claim 1, characterized in that: In the first simulation area of ​​the engineering simulation model, before using the first simulation calculation method to calculate the simulation parameters in the first simulation area to obtain the first calculation result, the method further includes: Receiving engineering parameters and constructing the engineering simulation model based on the engineering parameters; According to the material type of the first simulation area, the target calculation method is determined to be the first simulation calculation method, and the coupling characteristic parameters of the coupling area nodes between all simulation areas in the engineering simulation model are initialized to calculate the simulation parameters in the first simulation area based on the first simulation calculation method to obtain the first calculation result.

6. The method according to any one of claims 1 to 5, characterized in that: Generating the engineering structure simulation result according to one or more of the target coupling characteristic parameters includes: Determine the displacement information and force information of each simulation area in the engineering simulation model according to one or more target coupling characteristic parameters; One or more of a stress cloud map, a displacement distribution map, and a parameter table of the engineering simulation model are generated according to the displacement information and the force information.

7. The method according to claim 6, characterized in that The method further includes: generating a calculation report according to the engineering structure simulation result; the calculation report includes one or more of a model parameter summary, a simulation result summary, and a model stability evaluation result.

8. An engineering structure simulation device, characterized in that: include: A calculation module, configured to calculate, for a first simulation area of ​​an engineering simulation model, simulation parameters in the first simulation area using a first simulation calculation method to obtain a first calculation result; wherein the engineering simulation model includes a plurality of simulation areas; the first simulation area is any of the simulation areas; and the first calculation result represents stress state information and displacement state information of the first simulation area; A first determination module is used to determine the coupling characteristic parameters of the coupling area according to the first calculation result, the position information of the known points in the coupling area, and the position information of the coupling points in the coupling area; wherein the coupling area in the calculation method is the contact area between the first simulation area and the second simulation area, and the second simulation area is the simulation area adjacent to the first simulation area; the coupling characteristic parameters represent the displacement information and force information of the coupling points in the coupling area; an iteration module, configured to calculate the simulation parameters of the second simulation area according to a second simulation calculation method to obtain a second calculation result, and adjust the coupling characteristic parameters according to the second calculation result to obtain the iterated coupling characteristic parameters; The first determination module is further used to use the first simulation area as the second simulation area, repeatedly perform the steps of calculating the simulation parameters of the second simulation area according to the second simulation calculation method to obtain a second calculation result, and adjust the coupling characteristic parameters according to the second calculation result to obtain the coupling characteristic parameters after iteration, and when the error between the coupling characteristic parameters before iteration and the coupling characteristic parameters after iteration is less than a preset error threshold, determine the coupling characteristic parameters after iteration as target coupling characteristic parameters; A generation module is used to generate engineering structure simulation results according to one or more of the target coupling characteristic parameters.

9. A computer device, characterized in that: The computer device comprises a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the method according to any one of claims 1 to 7 when executing the computer program.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images