Engineering structure simulation method and device, equipment and storage medium

Abstract

The invention provides an engineering structure simulation method and device, equipment and a storage medium, and relates to the technical field of numerical simulation. The method comprises the steps of obtaining a first calculation result and a second result through calculation; determining coupling characteristic parameters of the coupling region according to the first calculation result; adjusting the coupling characteristic parameter according to the second calculation result to obtain an iterated coupling characteristic parameter; then repeatedly adjusting the coupling characteristic parameters to obtain target coupling characteristic parameters after iteration; and generating an engineering structure simulation result. When simulation is carried out, different simulation area calculation methods can be dynamically adjusted according to calculation results, coupling characteristic parameters of a contact area are adjusted in a sequential iteration mode, and coupling processing is automatically carried out on the calculation results obtained through the different calculation methods. The continuity and the coordination of the data of the contact area of the engineering simulation model can be improved, so that the efficiency and the accuracy of engineering structure simulation are improved.

Description

Technical Field

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

[0002] There are many types of numerical simulation software, such as Ansys, Abaqus and Flac3D. Different software handles different engineering problems and calculation methods. When an engineering problem cannot be simulated by a single calculation method for the engineering structure model, the simulation results are prone to errors. When multiple calculation methods are used to calculate the engineering structure, it is difficult to integrate multiple calculation methods on a single software platform and realize dynamic coupling. This leads to the numerical simulation of many complex projects having to rely on a single calculation method, which is difficult to achieve coupling and leads to low accuracy of simulation results. Although some engineering problems can be achieved by manually constructing the coupling transfer mechanism of the contact area using multiple software platforms at the same time to achieve coupling calculation, this method has a large workload and low efficiency in pre-calculation processing, and requires high professional knowledge from users, forming a high theoretical threshold, which limits 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 use different calculation methods to calculate different simulation areas in the engineering simulation model to obtain corresponding calculation results, which can improve the efficiency of engineering structure simulation calculation. When performing simulation, the coupling characteristic parameters of the contact area can be adjusted in a sequential iterative manner, and the coupling characteristic parameters of the calculation results of each adjacent area can be automatically calculated to achieve dynamic coupling processing of the calculation results of the contact area between simulation areas that use different calculation methods to calculate parameters, 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: 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.

[0005] In the above implementation, when simulating an engineering structure, different calculation methods can be automatically used to calculate different simulation areas in the engineering simulation model to obtain corresponding calculation results. By automatically calculating the coupling characteristic parameters of the calculation results of each adjacent area, and adjusting the coupling characteristic parameters of the contact area in a sequential iterative manner, dynamic coupling processing of the calculation results of the contact area between simulation areas using different calculation methods for parameter calculation is achieved, thereby achieving dynamic adjustment of simulation calculation methods for different simulation areas, which can improve the efficiency and accuracy of engineering structure simulation calculation.

[0006] In some embodiments, the coupling characteristic parameter includes a first coupling characteristic parameter, the first coupling characteristic parameter characterizing the displacement information of the coupling point of the coupling region, and determining the coupling characteristic parameter of the coupling point of the coupling region according to the first calculation result, the position information of the known point of the coupling region, and the position information of the coupling point of the coupling region may include: 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.

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

[0008] In some embodiments, 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; Determining 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 may include: 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.

[0009] 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 distributed from the known points in the coupling region to the coupling points in the coupling region in one direction, so that the force information of the coupling points in the coupling region can be determined according to the proportion of the force distributed from the known points in the coupling region to the coupling points in each direction. The second coupling characteristic parameter characterizing the force information of the coupling points in the coupling region is determined by the above method, and the automatic calculation of the second coupling characteristic parameter between different calculation methods and the dynamic coupling processing of the stress information of the contact area between different simulation areas are realized, so that the dynamic adjustment of the simulation calculation methods of different simulation areas can be realized automatically to realize the continuity and coordination of the stress in the contact area between different areas, thereby improving the efficiency and accuracy of the simulation calculation of the engineering structure.

[0010] In some embodiments, after calculating the simulation parameters in the first simulation area of ​​the engineering simulation model using the first simulation calculation method and obtaining the first calculation result, the method may further include: 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.

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

[0012] In some embodiments, before calculating the simulation parameters in the first simulation area of ​​the engineering simulation model using the first simulation calculation method to obtain the first calculation result, the method may further include: 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.

[0013] In the above implementation process, the target calculation method is determined as 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 efficiency and accuracy when processing different types of materials. Therefore, according to the appropriate calculation method, the calculation efficiency can be improved while ensuring accuracy.

[0014] In some embodiments, generating an engineering structure simulation result according to one or more of the target coupling characteristic parameters may include: 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.

[0015] In some embodiments, the method may further include: generating a calculation report based on the engineering structure simulation results; the calculation report includes one or more of a model parameter summary, a simulation result summary, and a model stability assessment result.

[0016] In a second aspect, an embodiment of the present application provides an engineering structure simulation device, which may 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.

[0017] In a third aspect, an embodiment of the present application provides a computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method described in the above description when executing the computer program.

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

[0019] Compared with the prior art, the beneficial effects of the present application are: when simulating an engineering structure, different calculation methods can be automatically used 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, and the dynamic coupling processing of the calculation results can be realized by calculating the coupling characteristic parameters of the calculation results of each adjacent area, which can automatically realize the continuity and coordination of the contact area data of each simulation area, thereby improving the efficiency and accuracy of the engineering structure simulation calculation. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] Figure 1 A schematic diagram of the steps of the engineering structure simulation method provided in an embodiment of the present application.

[0021] Figure 2 A schematic diagram of a simulation model with multiple simulation areas provided in an embodiment of the present application.

[0022] Figure 3 A schematic diagram of a soil layer geometric model provided in an embodiment of the present application.

[0023] Figure 4 A stress cloud diagram of the generated soil layer geometric model provided in an embodiment of the present application.

[0024] Figure 5A schematic diagram of an engineering structure simulation device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0025] The present application is further described in detail below in conjunction with test examples and specific implementation methods. However, this should not be understood as the scope of the above subject matter of the present application being limited to the following embodiments, and all technologies implemented based on the content of the present application belong to the scope of protection of the present application.

[0026] Unless otherwise specified, in the description of the specific embodiments of the present application, the terms indicating the orientation or position relationship such as "up", "down", "left", "right", "center", "inside", "outside", "side", etc. are all expressions based on the orientation or position relationship shown in the drawings, or the orientation or position relationship when the product / equipment / device is usually used. These terms of orientation or position relationship are only for the convenience of describing the scheme of the present application or simplifying the description in the specific embodiments to facilitate the technicians to quickly understand the scheme, rather than indicating or implying that a specific device / component / element must have a specific orientation, or be constructed and operated in a specific position relationship, and therefore cannot be understood as a limitation on the present application.

[0027] In the description of the embodiments of the present application, the technical terms "first", "second", etc. 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 number, 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 "multiple" is two or more, unless otherwise clearly and specifically defined.

[0028] Reference to "embodiments" herein means that a particular feature, structure, or characteristic described in conjunction with the embodiments may be included in at least one embodiment of the present application. The appearance of the phrase in various locations in the specification does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment that is mutually exclusive with other embodiments. It is explicitly and implicitly understood by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0029] Example 1 During the research, the applicant found that when using simulation software such as Ansys, Abaqus and Flac3D to simulate engineering structures, if different calculation methods are needed to calculate the engineering structure model, cumbersome steps such as geometric modeling and meshing are required. The calculation method to be used for each area can only 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 actual engineering problems, in order to ensure the accuracy of simulation calculations, an engineering model may need to use more calculation methods to calculate different areas. At this time, the above software cannot realize the dynamic coupling of multiple calculation methods. Therefore, after studying this problem, the applicant proposed a method for simulating engineering structures. When simulating engineering structures, the coupling characteristic parameters of the calculation results of each adjacent area are automatically calculated to realize the dynamic coupling processing of the calculation results of the contact area between the simulation areas that use different calculation methods for parameter calculation, thereby improving the accuracy of engineering structure simulation.

[0030] Please see Figure 1 , Figure 1 A schematic diagram of the steps of the engineering structure simulation method provided in the embodiment of the present application. The engineering structure simulation method may include: S1. For a first simulation area of ​​an engineering simulation model, use a first simulation calculation method to calculate simulation parameters in the first simulation area to obtain a first calculation result.

[0031] The engineering simulation model includes multiple simulation areas; the first simulation area is any simulation area; and the first calculation result represents the stress state information and displacement state information of the first simulation area.

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

[0033] 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 point in the coupling region.

[0034] S3. Calculate 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 iterated coupling characteristic parameters.

[0035] S4. Taking the first simulation area as the second simulation area, repeatedly 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 coupling characteristic parameters after iteration. When the error between the coupling characteristic parameters before iteration and the coupling characteristic parameters after iteration is less than a preset error threshold, determining the coupling characteristic parameters after iteration as the target coupling characteristic parameters.

[0036] S5. Generate engineering structure simulation results according to one or more target coupling characteristic parameters.

[0037] 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; automobile collision safety analysis, simulation of vehicle behavior 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.

[0038] The engineering simulation model can be constructed by the user inputting simulation values ​​in the relevant interface. Figure 2 , Figure 2 A schematic diagram of a simulation model with multiple simulation areas provided for an embodiment 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, a different calculation method 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.

[0039] 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 using different calculation methods for simulation calculations, the calculation results need to be coupled between regions of different calculation methods, and channels for data transmission between different regions are established by creating contacts between regions. In the contact area, data transmission rules and conversion mechanisms need to be defined, which includes 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 that the calculation results between different regions remain numerically consistent.

[0040] When performing coupling processing, the displacement information and force information of the nodes in the coupling area are characterized by determining the coupling characteristic parameters of the coupling area, and then the coupling characteristic parameters are iterated. When the iteration is completed, the target coupling characteristic parameters are obtained to accurately characterize the displacement information and force information in the contact area of ​​each adjacent simulation area, so that the engineering structure simulation results can be generated according to the target coupling characteristic parameters, which can improve the continuity and coordination between the contacts in different simulation areas. The specific implementation of the coupling processing can be referred to the description in the following embodiment, which will not be expanded here.

[0041] The process of iterating the coupling characteristic parameters is specifically to correct the parameters such as force and displacement in the coupling characteristic parameters according to the second calculation result to reflect the influence of the second simulation area on the coupling area. Then the adjusted coupling area characteristic parameters are input into the calculation model of the first simulation area, so that the influence of the second simulation area is transferred to the first simulation area through the adjustment of the coupling area characteristic parameters. Then, based on the adjusted coupling area characteristic parameters, the model parameters of the first simulation area and the first calculation method, the first simulation area is simulated and calculated. By repeatedly executing step S3, the coupling calculation is performed in a sequential iterative manner until the error of the two iterations is less than the preset error threshold, then the iteration is stopped to obtain the target coupling characteristic parameters that accurately characterize the displacement information and force information in the contact area of ​​each adjacent simulation area.

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

[0043] In the above implementation, when simulating an engineering structure, different calculation methods can be automatically used for different simulation areas in the engineering simulation model to obtain corresponding calculation results, which can improve the efficiency of engineering simulation calculations. When simulating an engineering structure, the calculation methods for different simulation areas in the engineering simulation model can be dynamically adjusted according to the calculation results, and the dynamic coupling processing of the calculation results of the contact areas between the simulation areas using different calculation methods can be realized by automatically calculating the coupling characteristic parameters of the calculation results of each adjacent area, thereby improving the efficiency and accuracy of the engineering structure simulation calculation.

[0044] Example 2 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-mentioned embodiment 1.

[0045] In the embodiment of the present application, the coupling characteristic parameter includes a first coupling characteristic parameter, and the first coupling characteristic parameter represents the displacement information of the coupling point in the coupling region.

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

[0047] The embodiment of the present application uses the finite element-material point method coupling as an example for explanation. When other calculation methods are used for coupling, this method can also be referred to for data coupling. Among them, the finite element-material point method coupling is to use the finite element method to discretize small deformation objects, and the material point method to discretize large deformation objects, and the different discrete areas are coupled to each other through the contact algorithm.

[0048] The step of determining the coupling characteristic parameters of the coupling region according to the first calculation result and the second calculation result may include: S21. Determine the sum of the reciprocals of distances from the known point in the first coupling region to all coupling points in the first coupling region.

[0049] The first coupling region known point is any coupling region known point; the first coupling region coupling point is any coupling region coupling point.

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

[0051] The first weight coefficient represents the influence of the first coupling region known point on the displacement of the first coupling region coupling point. The first coupling region known point is any coupling region known point; the first coupling region coupling point is any coupling region coupling point; the first weight coefficient represents the displacement distribution coefficient transmitted from the first coupling region known point to the first coupling region coupling point.

[0052] In an embodiment of the calculation using the finite element-material point method, in the first iteration, the material point method area is calculated first, and then the finite element method area is calculated. The known points in the coupling area are the material points, and the coupling points in the coupling area are the finite element points.

[0053] At the coupling interface, in order to make the displacement continuous, the displacement information of the material point method area needs to be transferred to the finite element method area. By setting shared nodes or interpolation points on the coupling interface, the displacement of the material point in the material point method area close to the coupling interface is distributed to the nodes in the finite element method area by interpolation. For example, assuming there are material point i and finite element point j on the coupling interface, the weight coefficient w can be established according to the distance between them. ij :

[0054] Among them, r ij is the distance from material point i to finite element point j, n is any finite element point (with values ​​ranging from 1, ..., k), k is the number of all finite element points in the coupling region, r in represents the distance from material point i to finite element point n, Represents the sum of the inverses of the distances from material point i to all finite element points. ij When the distance r from the material point to all finite element points is in The reciprocal of the distance between the material point and the finite element point is summed up, and the ratio of the reciprocal of the distance between the material point and the finite element point to the sum is used as the weight w of the displacement of the material point i to the finite element point j. ij In this way, the displacement of the finite element point can be updated according to the displacement of the material point. Based on the weights assigned to the displacement of all material points to the finite element point j, the displacement of the finite element point can be obtained .in, is the displacement of material point i.

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

[0056] On the coupling interface, the displacement should be continuous. If the displacement in the material point method region is u MPM , the displacement of the finite element method area is u FEM , that is u MPM =u FEM , if the displacement of the material point method region is obtained by the motion of the material point, for the material point i close to the coupling interface, its displacement , where x i is the current position of material point i, is the initial position of material point i.

[0057] In the finite element method area, the displacement u in the finite element method area FEM It can be obtained by finite element point displacement interpolation, that is:

[0058] Among them, N jIt is the unit shape function, which is defined based on the unit's geometric shape and node position. It can describe the distribution of displacement inside the unit. The unit is the basic geometric body or element that constitutes the entire finite element method area; ε, μ, δ are the natural coordinates within the unit, d j is the first coupling characteristic parameter, i.e. the displacement of the finite element point.

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

[0060] Since different calculation methods use different data transmission methods, in order to ensure the continuity and harmony of the data in the contact areas between different areas, it is necessary to determine the data transmission mechanism between different areas, so as to couple the calculation results of different areas. Among them, data transmission can also include force transmission, that is, determining the force information of the coupling points in the coupling area on the coupling interface. In the embodiment of the present application, the finite element-material point method coupling is continued as an example to provide a simulation method for coupling the forces of the calculation results of different areas of the engineering structure model.

[0061] The step of determining the coupling characteristic parameters of the coupling region according to the first calculation result and the second calculation result may further include: S24: Determine the sum of the distances from the known point in the second coupling region to all the coupling points in the second coupling region.

[0062] 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.

[0063] The transmission of force is the key to ensure the conservation of momentum. When the finite element method region is deformed, the force generated needs to be transmitted to the material point method region, and vice versa. On the coupling interface, the force in the material point method region 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.

[0064] There are n material points and m finite element points on the coupling interface. Assume that the force on material point i is:

[0065] The equivalent nodal force at finite element point j is:

[0066] We can assume that there exists an allocation matrix T whose elements t ij It represents the proportionality coefficient of the force distribution of material point i to node j.

[0067] S25. Determine a first allocation coefficient according to the sum of the distances.

[0068] 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 known point of the second coupling region is any known point of the coupling region, and the coupling point of the second coupling region is any coupling point of the coupling region.

[0069] The first direction may be the x direction of the spatial coordinate system. In the following embodiments, the first distribution coefficient t ijx Take the ratio of the force distributed to the current finite element point in the x direction by the distribution matrix T as an example to illustrate: At 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 is MPM-FEM The force F acting on the material point method area by the finite element method area FEM-MPM Equal in size and opposite in direction, that is:

[0070] In the material point method, the force can be expressed by the stress tensor σ MPM and the area A of the material point MPM To calculate, that is:

[0071] In the finite element method, the force is expressed through the element stiffness matrix K FEM and the node displacement d FEM To calculate, that is:

[0072] The transmission of force is the key to ensure the conservation of momentum. When the finite element method area is deformed, the force generated needs to be transmitted to the material point method area, and vice versa.

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

[0074] There are n material points and m finite element points on the coupling interface. Assume that the force on material point i is:

[0075] The equivalent nodal force at finite element point j is:

[0076] Assume that there is a distribution matrix T whose element t ij It represents the proportionality factor of the force distribution of material point i to finite element point j.

[0077] According to the force equivalence principle, the equivalent nodal force at finite element point j is is the sum of all material point forces distributed to the node. The equivalent nodal force at 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, which indicates the proportion of the x-direction force of material point i distributed to node j, that is, the proportion of the force of all second material points distributed to the current finite element point in the first direction; F ix is the force on material point i in the x direction.

[0078] S26. Determine a second coupling characteristic parameter according to the allocation matrix.

[0079] The allocation matrix is ​​obtained based on the allocation coefficients in each direction. The allocation matrix T can be composed of the first allocation coefficient in the x direction of the spatial coordinate system, the second allocation coefficient in the y direction, and the ijy and the third distribution coefficient t in the z direction ijz Sure.

[0080] The following is an example of constructing a distribution matrix based on geometric relationships. The first distribution coefficient can be determined according to the distance between the material point and the finite element point in the first direction. Suppose the vector from material point i to node j is , 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 embodiment of the present application, the first direction is the x direction of the spatial coordinate system, the second direction is the y direction of the spatial coordinate system, and the third direction is the z direction of the spatial coordinate system. The template is:

[0081] 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 separately and the third distribution coefficient After determining the distribution coefficient in each direction, we can determine the ijx ,t ijy and t ijzThe constructed matrix is ​​the allocation matrix T.

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

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

[0084] In the above implementation process, the distribution coefficient determined by the sum of the distances between the known point in the coupling region and the coupling point in the coupling region characterizes the proportion of the force distributed from the known point in the coupling region to the coupling point in the coupling region in one direction, so that the force information of the coupling point in the coupling region can be determined according to the proportion of the force distributed from the known point in the coupling region to the coupling point in each direction. By determining the second coupling characteristic parameter characterizing the force information of the coupling point in the coupling region in the above manner, the continuity and coordination of the data transmission of the contact area between different areas can be automatically realized, thereby improving the efficiency and accuracy of the engineering structure simulation calculation.

[0085] Example 3 This embodiment is an implementation method of dynamically adjusting the calculation method of different simulation areas according to the simulation results in the above embodiment 1.

[0086] In a first simulation area of ​​the engineering simulation model, after using a first simulation calculation method to calculate simulation parameters in the first simulation area and obtaining a first calculation result, the method may further include: When a first calculation result obtained by calculating the first simulation region based on the first simulation calculation method fails to converge, it is determined to use a third calculation method to calculate the first simulation region.

[0087] During the calculation process, the system can automatically modify the calculation methods of different areas according to the stress and deformation conditions obtained during the simulation calculation process. For example, if the engineering simulation model includes two simulation areas A and B, the finite element method is used to calculate area A, but the system finds that the calculation cannot converge during the calculation process, then the stress or deformation representing the first simulation area exceeds the set tolerance range. The system can determine the conversion of the first calculation method to a third calculation method that is more suitable for the area based on the engineering parameters of this area. For example, the calculation method of area A can be adjusted from the finite element method to the material point method, thereby recalculating the entire model.

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

[0089] Example 4 This embodiment is a specific implementation plan of receiving engineering simulation parameters and constructing an engineering simulation model in the above-mentioned embodiment 1.

[0090] 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 provided in the embodiment of the present application may further include: receiving engineering parameters and constructing an 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 a first calculation result.

[0091] In an embodiment 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 engineering parameters input by a user through a web page, and automatically create an engineering simulation model according to the engineering parameters, automatically select a calculation engine, automatically divide the grid, and automatically perform calculations based on the selected calculation engine.

[0092] Among them, the engineering structure simulation system includes multiple computing engines, each of which can support one or more computing methods, such as finite element, finite difference, discrete element, material point, etc. Different computing engines can extract simulation parameters from the database through the data interface, and automatically create geometric models, assign material parameters, add boundary conditions, and divide computing units according to material areas; each computing engine can calculate different areas of the model at the same time, and can use distributed parallel computing methods to speed up the simulation.

[0093] During the calculation process, the system can automatically select the calculation engine for each area 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 material does not contain a calculation method, the pre-set unified calculation method can be used for calculation. If the calculation results of a material area show that the stress or deformation exceeds the set tolerance range, and the system has been pre-set to allow the calculation method to be changed, the system will automatically convert the calculation method for this area to a calculation method that is more suitable for this area.

[0094] The engineering parameters input by the user may involve aspects such as material parameter creation, geometric model parameter input, mesh division, boundary condition entry, and calculation method setting. For example, the user can enter the geometric model and material parameters of the project on the web page, and complete the setting of boundary conditions, mesh refinement, and calculation methods through dragging and selecting the interface.

[0095] Among them, in terms of material parameter creation, users can set parameters for different materials, such as elastic modulus, Poisson's ratio, etc., and specify 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 areas, and allow users to set conversion rules for calculation methods. In terms of meshing, the system can automatically divide the mesh based on the engineering geometric model parameters entered by the user in the web interface, and receive the size, shape, mesh accuracy and other data of the geometric model defined by the user according to the actual engineering needs. In terms of boundary condition entry, the system can receive boundary conditions specified by the user for different areas of the simulation model, such as concentrated loads, distributed loads, fixed boundaries, free boundaries, etc., to ensure that the calculation boundary of the model meets the actual engineering needs. In terms of calculation method setting, the system can provide users with the option of selecting the initial calculation method of the overall model, the unified calculation method of the simulation calculation, and receive settings for whether to allow the calculation method to be dynamically transformed during the simulation process to improve the simulation accuracy.

[0096] The system can realize real-time visualization of the model through WebGL technology. Users can directly interact with the model through dragging, zooming, etc., and modify the input simulation parameters. All parameters are automatically saved in the server-side database after confirmation.

[0097] In the above implementation process, the target calculation method is determined as 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 efficiency and accuracy when processing different types of materials. Therefore, according to the appropriate calculation method, the calculation efficiency can be improved while ensuring accuracy.

[0098] Example 5 This embodiment is a specific implementation scheme for determining the engineering structure simulation result according to one or more coupling characteristic parameters in Embodiment 1.

[0099] Determining engineering structure simulation results based on one or more coupled characteristic parameters may include: 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; and generate one or more of the stress cloud map, displacement distribution map, and parameter table of the engineering simulation model according to the displacement information and force information.

[0100] 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 results in the form of a three-dimensional cloud map on the web page.

[0101] In some embodiments, the engineering simulation results include stress cloud maps, displacement distribution maps, parameter tables, etc. Among them, in the stress cloud map, the stress distribution of each area can be drawn on the model according to the predefined color interval based on the unit stress value. The displacement values ​​of each unit of the model in the displacement distribution map are also presented in color intervals, and the user can view the overall deformation of the model through the interface. In the parameter table, the system can display the detailed values ​​of the unit's displacement, stress, strain, etc. 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.

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

[0103] 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 position in the model to help users understand the calculation results of specific areas. The model stability evaluation result can evaluate the overall stability of the model, and can include the analysis results of key data such as plastic zone, stress concentration zone and maximum displacement.

[0104] After obtaining the simulation results, users can also choose different simulation schemes. By systematically comparing the key results (such as stress, deformation, stability analysis, etc.) in different simulation calculation reports, engineers and technicians can quickly find the optimal design scheme. They can also automatically optimize the model based on the comparative analysis results, such as strengthening the mesh refinement or adjusting the material properties in areas with larger deformation, thereby optimizing the simulation effect.

[0105] Example 6 The following is an embodiment of the present application for analyzing the stability of a retaining wall on a roadbed slope.

[0106] In the case of performing a retaining wall stability analysis on a roadbed slope, the material parameters of the model may be first input from the web page, and the material parameters may specifically include material type, material name, failure criterion, gravity, elastic modulus, Poisson's ratio, etc. In some embodiments, a variety of model material parameters may be pre-stored in the system material library, which may be directly called when used, or newly defined when used.

[0107] After defining the material parameters, you can enter the geometric model parameters. The toe vertex of the retaining wall is the base point of the model, which is also the origin of the calculated coordinates. Figure 3 , Figure 3 The schematic diagram of the soil layer geometry model provided in the embodiment of the present application is a two-dimensional analysis, and the soil layer geometry model is a polyline enclosed by multiple vertex coordinates. The soil layer geometry model includes multiple material areas such as backfill soil, silty clay, and retaining walls. The stress conditions of the soil layer geometry model are shown by the arrows in the figure. Then enter the boundary condition parameters, and the distributed load parameters include the starting position of 8.1m, the end position of 14.1m, the starting pressure of 150kPa, the end pressure of 150kPa, the angle type of radial, the angle of 90°, etc.

[0108] Then, calculations are performed on each area and simulation analysis is carried out. The calculation method can refer to the description of the above embodiment and will not be repeated here. After the calculation is completed, a stress cloud diagram of the soil layer geometric model can be generated according to the calculation results. Figure 4 , Figure 4 The stress cloud map of the generated soil layer geometry model provided in the embodiment of the present application. By calculating the above-mentioned various material areas, the stress distribution in the soil layer geometry model can be obtained, which mainly includes five stress levels, namely stress level 1, stress level 2, stress level 3, stress level 4 and stress level 5. Each stress level represents a stress range, wherein the maximum stress point is the MAX point and the minimum stress point is the MIN point.

[0109] Example 7 Please see Figure 5 , Figure 5 Schematic diagram of an engineering structure simulation device provided in an embodiment of the present application. The engineering structure simulation device 50 may include: The calculation module 51 is used 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 simulation area; the first calculation result represents the calculation method of the stress state information and displacement state information of the first simulation area.

[0110] The first determination module 52 is used 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.

[0111] An iteration module 53 is used to determine a second simulation calculation method for a second simulation region according to the coupling characteristic parameters, and calculate the simulation parameters of the second simulation region by the second simulation calculation method to obtain the iterated coupling characteristic parameters; The first determination module 52 may also be used to use the first simulation area as the second simulation area, repeatedly perform the second simulation calculation method for determining the second simulation area according to the coupling characteristic parameters, and calculate the simulation parameters of the second simulation area by 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 a preset error threshold; The generating module 54 is used to generate engineering structure simulation results according to one or more target coupling characteristic parameters.

[0112] Optionally, the first determining module 52 may be specifically configured to: Determine the sum of the reciprocals of the distances from the known point in the first coupling region to the coupling points in all coupling regions; Determine a first weight coefficient according to the sum of the reciprocals, the first distance between the first coupling region known point and the first coupling region coupling point; wherein the first coupling region known point is any coupling region known point; the first coupling region coupling point is any 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.

[0113] Optionally, the first determining module 52 may be specifically configured to: Determine the sum of the distances from the second coupling region known point to all coupling region coupling points; the second coupling region known point is any coupling region known point; The first distribution coefficient is determined according to the sum of the distances; 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 known point of the second coupling region is any known point of the coupling region, and the coupling point of the second coupling region is any coupling point of the coupling region; The second coupling characteristic parameter is determined according to the allocation matrix; wherein the allocation matrix is ​​obtained according to the allocation coefficients in each direction.

[0114] Optionally, the first determining module 52 may also be used to: When a first calculation result obtained by calculating the first simulation region based on the first simulation calculation method fails to converge, it is determined to use a third calculation method to calculate the first simulation region.

[0115] Optionally, the engineering structure simulation device 50 may further include a receiving module for receiving engineering parameters and constructing an 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 a first calculation result.

[0116] Optionally, the generating module 54 may be specifically configured to: Determine displacement information and force information of each simulation area in the engineering simulation model according to one or more target coupling characteristic parameters; Generate one or more of the following: a stress cloud diagram, a displacement distribution diagram, and a parameter table of the engineering simulation model according to the displacement information and the force information.

[0117] It should be understood that, when performing engineering structure simulation calculations, the various modules of the engineering structure simulation 30 provided in the above embodiment are only illustrated by the division of the various functional modules in the above description. In actual applications, the above functions can be assigned to different functional modules as needed, 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.

[0118] The functional modules in the above embodiments may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit, and the above integrated unit may be implemented in the form of hardware or in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the protection scope of the embodiments of the present application.

[0119] Based on the same application concept, an embodiment of the present application also provides 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, the engineering structure simulation method as described in the above description is implemented.

[0120] Based on the same application concept, 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 engineering structure simulation method described in the above description is implemented.

[0121] The embodiments described above 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 aforementioned embodiments, a person skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some of the technical features may be replaced by equivalents. Such modifications or replacements do not deviate the essence of the corresponding technical solutions 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