Method, system and equipment for post-processing simulation of force-bearing structures based on numerical calculation

Abstract

The present invention discloses a post-processing method, system and device for simulation of a force-bearing structure based on numerical calculation, and belongs to the technical field of structural deformation simulation. In the existing post-processing simulation scheme, the cloud map display effect will be disturbed by the overall displacement, resulting in the difficulty of the generated cloud map to accurately reflect the actual situation of the deformation of the structure itself. A post-processing method for simulation of a force-bearing structure based on numerical calculation of the present invention identifies the overall motion data by creating a deformation data acquisition unit, an interference motion identification unit, a local feature simulation unit, and a structural cloud map generation unit, and removes the overall motion data from the force motion results to obtain the self-deformation data, thereby being able to effectively distinguish between real deformation and overall motion phenomena, so that the self-deformation characteristics of the structural model caused by the load will be highlighted when the cloud map is displayed, which is helpful to quickly locate structural problems, improve the quality and efficiency of engineering design analysis, and meet the actual needs in engineering design.

Description

Technical Field

[0001] The invention relates to a force structure simulation post-processing method, system and equipment based on numerical calculation, and belongs to the technical field of structure deformation simulation. Background Art

[0002] Static analysis is an important tool in engineering design and scientific research, used to simulate the stress and deformation of complex structures under static loads. Through numerical calculations, engineers can accurately evaluate the strength, stiffness and stability of the structure, thereby optimizing the design and predicting potential failure risks. After the analysis is completed, the main task of the post-processing stage is to display the deformation and stress distribution of the structural model under load through visualization technology, among which cloud maps are the most commonly used form of expression.

[0003] Therefore, the quality of the cloud map directly affects the interpretation effect of the post-processing stage. Ideally, the cloud map should be able to clearly show the deformation characteristics and stress distribution of the structure itself, allowing engineers to quickly locate problem areas, identify potential design defects, and provide a basis for subsequent structural optimization.

[0004] However, in actual engineering, the structural model that requires static analysis is usually an assembly structure model containing multiple bodies. Several individuals in the assembly structure model are combined together through different connection methods (such as contact, welding or bolt connection), but due to the imperfect connection characteristics and degree of freedom settings, the overall stiffness of the assembly structure model may be low. In this case, even if sufficient degree of freedom restrictions are imposed on the entire assembly structure model, some individuals will still have a large overall displacement under the action of load, thereby masking the deformation characteristics of the model itself.

[0005] Furthermore, during the loading process, the structural model may exhibit an overall translation or rotation, making it difficult to clearly display local subtle deformations in the cloud map. Therefore, in actual engineering applications, the cloud map display effect will be disturbed by the overall displacement, resulting in the generated cloud map being difficult to accurately reflect the actual situation of the structure's own deformation and unable to meet the actual needs of engineering design.

[0006] The information disclosed in this Background Art is only for understanding the background of the inventive concept and therefore it may include information that does not constitute the prior art. Summary of the invention

[0007] In response to the above problem or one of the above problems, an object of the present invention is to provide a post-processing method for simulation of a force-bearing structure based on numerical calculation, by creating a deformation data acquisition unit, an interference motion identification unit, a local feature simulation unit, and a structural cloud map generation unit, identifying the overall motion data, and eliminating the overall motion data from the force motion results to obtain its own deformation data, and then converting its own deformation data into color scale distribution information to obtain a structural deformation cloud map that intuitively displays the results of static analysis, thereby being able to effectively distinguish between real deformation and overall motion phenomena, ensuring that the overall position movement in the post-processing result data is effectively shielded, so that when the cloud map is displayed, the deformation characteristics of the structural model itself caused by the load can be highlighted, rather than meaningless overall motion data, so that the local subtle deformation of the structural model can be clearly displayed from the cloud map, meeting the actual needs in engineering design.

[0008] In response to the above problem or one of the above problems, the second object of the present invention is to provide a post-processing method, system and equipment for simulation of a load-bearing structure based on numerical calculation, which can intuitively reflect the actual stress and deformation of each area of ​​the structural model, making the color gradation distribution of the cloud map more reasonable; at the same time, it also amplifies the effect of the model's own stress and deformation, avoids the problem of color gradation uniformity caused by overall large displacement, and greatly improves the data readability and utilization rate of the post-processing results. Engineers can more efficiently obtain and analyze the key features of the model deformation, which helps to quickly locate structural problems, improve the quality and efficiency of engineering design analysis, and meet the finite element analysis needs under complex working conditions.

[0009] To achieve one of the above purposes, the first technical solution of the present invention is: A post-processing method for simulation of a force-bearing structure based on numerical calculation comprises the following steps: Step 1: Obtain the force motion results of the structural model through the deformation data acquisition unit created in advance; Step 2: Using the previously created interference motion identification unit, based on the node state change, the force motion results are processed to obtain the overall motion data; Step 3: Using the local feature simulation unit created in advance, the overall motion data is removed from the force motion results to obtain the self-deformation data, which is used to characterize the stress distribution and deformation of the structural model itself; Step 4: Use the previously created structural cloud map generation unit to convert its own deformation data into color scale distribution information to obtain the structural deformation cloud map, thereby realizing the post-processing of the force-bearing structure simulation based on numerical calculation.

[0010] The present invention identifies the overall motion data by creating a deformation data acquisition unit, an interference motion identification unit, a local feature simulation unit, and a structural cloud map generation unit, and eliminates the overall motion data from the force motion results to obtain the self-deformation data, and then converts the self-deformation data into color scale distribution information to obtain a structural deformation cloud map that intuitively displays the static analysis results, thereby being able to effectively distinguish between real deformation and overall motion phenomena, and ensure that the overall position movement in the post-processing result data is effectively shielded. Therefore, when displaying the cloud map, the deformation characteristics of the structural model itself caused by the load can be highlighted, rather than the meaningless overall motion data, so that the local subtle deformation of the structural model can be clearly displayed from the cloud map, meeting the actual needs in engineering design.

[0011] Furthermore, the structural deformation cloud map obtained by the present invention can intuitively reflect the actual stress and deformation of each area of ​​the structural model, making the color scale distribution more reasonable, while also amplifying the effect of the model's own stress deformation, avoiding the problem of color scale uniformity caused by overall large displacement. Therefore, through the data processing and optimization method of the present invention, the data readability and utilization rate of the post-processing results are greatly improved, and engineers can more efficiently obtain and analyze the key features of model deformation, which helps to quickly locate structural problems, improve the quality and efficiency of engineering design analysis, and meet the needs of finite element analysis under complex working conditions.

[0012] As the preferred technical measures: Step 1: The method for obtaining the force motion results of the structural model through the deformation data acquisition unit created in advance is as follows: obtaining the concentrated load to be applied; Based on the structural model, set fixed boundary information; Divide the structural model into grids to obtain several grids; The vertices of the mesh are taken as nodes and numbered to obtain node information; Apply concentrated loads to the structural model, and perform static analysis on the nodes of the grid based on the fixed boundary information. Track the state of the nodes at each time step to obtain the state data of each node for monitoring the state changes of the nodes. The status data and node information of several nodes are aggregated to obtain the force and motion results of the structural model.

[0013] As the preferred technical measures: Step 2: Using the previously created interference motion identification unit, based on the node state change, the force motion results are processed to obtain the overall motion data as follows: Classify the results of force motion and get the motion type. The types of motion include at least overall translation, rotation about a point or rotation about a line; Based on the type of movement, obtain the corresponding discriminant parameters, The discrimination parameters at least include node stress value and / or contact relationship value; Based on the force motion results, obtain the state data of the node at each time step; According to the motion type, the node stress value and / or the contact relationship value are selected, and the state data at each time step is judged to obtain the overall motion position at the end time of the overall motion; Subtract the overall motion position from the initial position of the corresponding node to obtain the overall motion data of one or more nodes.

[0014] As the preferred technical measures: According to the motion type, the node stress value is selected, and the state data at each time step is judged to obtain the overall motion position at the end of the overall motion as follows: The motion type is overall translation, and the state data at each time step includes node stress values ​​and displacement values; The node stress value at each time step is judged. When the node stress value is greater than zero, it indicates that the overall movement is finished, and the time step at this time is recorded; According to the time step at this time, the displacement value of the corresponding node, that is, the overall motion position, is obtained.

[0015] As the preferred technical measures: According to the motion type, the contact relationship value is selected, and the state data at each time step is judged to obtain the overall motion position at the end of the overall motion as follows: The motion type is rotation around a point or rotation around a line; the contact relationship value is zero or the gap threshold, which is used to characterize the contact state of the node, which is divided into real contact and pseudo contact; real contact means that the node is in clear contact with a neighboring individual within a certain time step and satisfies the contact constraint condition; pseudo contact means that the node does not directly contact the neighboring individual, but when the distance between the two is less than the gap threshold, it is assumed that there is a contact relationship; Get the state data at each time step, which includes the distance value between the node and the adjacent individuals; According to the distance value at each time step, the contact state between the node and the adjacent individuals is judged to obtain the contact judgment result, which includes no contact, pseudo contact or real contact; When the node changes from no contact to pseudo contact or real contact, it indicates that the overall motion ends, and the end time step at this time is recorded; According to the end time step, the position change value of the corresponding node, that is, the overall motion position, is calculated.

[0016] As the preferred technical measures: According to the end time step, the method to calculate the position change value of the corresponding node is as follows: Obtain the rotation center of the structural model and the three-dimensional coordinates of the nodes when the initial time is zero; According to the end time step, determine the real-time coordinate value of the node; Calculate the rotation difference between the real-time coordinates of the node and the rotation center, and establish the unit vector of the line direction; Construct a rotation matrix based on the unit vector; Multiply the rotation matrix by the rotation difference to get the position change value of the corresponding node.

[0017] As the preferred technical measures: The method for judging the contact state between a node and an adjacent object is as follows: When the distance value is greater than the gap threshold, there is no contact, indicating that the nodes have neither real contact nor pseudo contact, and the position change of the nodes mainly comes from the overall movement; When the distance value is greater than zero and less than or equal to the gap threshold, it is a pseudo contact; When the distance value equals zero, it is a true contact.

[0018] As the preferred technical measures: Step 4: Use the previously created structural cloud map generation unit to convert the deformation data into color distribution information. The method for obtaining the structural deformation cloud map is as follows: Get the self-deformation data of all nodes; Map the deformation data to the standardized interval [0,1] to obtain the deformation standard data; Set a color gradient, which includes at least a first color, a second color, and a third color, and is used to express continuous changes in the deformation size; the first color is used to represent the low value in the deformation data itself, the second color is a transition color, and the third color is used to represent the high value in the deformation data itself; According to the color gradient, the deformed standard data is processed using the linear interpolation method to obtain the corresponding color value; Assign color values ​​to corresponding nodes to obtain color scale distribution information; According to the color distribution information, the structural deformation cloud map is drawn.

[0019] To achieve one of the above purposes, the second technical solution of the present invention is: A stress structure simulation post-processing system based on numerical calculation, comprising a deformation data acquisition module, an interference motion identification module, a local feature simulation module and a structure cloud map generation module; Deformation data acquisition module, used to obtain the force and motion results of the structural model; Interference motion identification module, used to process the force motion results based on the node state changes to obtain the overall motion data; The local feature simulation module is used to remove the overall motion data from the force motion results and obtain the self-deformation data to characterize the stress distribution and deformation of the structural model itself; The structural cloud map generation module is used to convert the self-deformation data into color-scale distribution information to obtain the structural deformation cloud map.

[0020] The present invention can intuitively reflect the actual stress and deformation conditions of each area of ​​the structural model by setting a deformation data acquisition module, an interference motion recognition module, a local feature simulation module, and a structural cloud map generation module, so that the color scale distribution of the cloud map is more reasonable, and the effect of the model's own force deformation is also amplified, avoiding the problem of color scale uniformity caused by overall large displacement.

[0021] Furthermore, through the data processing and optimization method of the present invention, the data readability and utilization rate of the post-processing results are greatly improved, and engineers can more efficiently obtain and analyze the key features of the model deformation, which helps to quickly locate structural problems, improve the quality and efficiency of engineering design analysis, and meet the finite element analysis needs under complex working conditions.

[0022] To achieve one of the above purposes, the third technical solution of the present invention is: An electronic device comprising: one or more processors; A storage device for storing one or more programs; When the one or more programs are executed by the one or more processors, the one or more processors implement the above-mentioned force-bearing structure simulation post-processing method based on numerical calculation.

[0023] Compared with the prior art solutions, the present invention has the following beneficial effects: The present invention identifies the overall motion data by creating a deformation data acquisition unit, an interference motion identification unit, a local feature simulation unit, and a structural cloud map generation unit, and eliminates the overall motion data from the force motion results to obtain the self-deformation data, and then converts the self-deformation data into color scale distribution information to obtain a structural deformation cloud map that intuitively displays the static analysis results, thereby being able to effectively distinguish between real deformation and overall motion phenomena, and ensure that the overall position movement in the post-processing result data is effectively shielded. Therefore, when displaying the cloud map, the deformation characteristics of the structural model itself caused by the load can be highlighted, rather than the meaningless overall motion data, so that the local subtle deformation of the structural model can be clearly displayed from the cloud map, meeting the actual needs in engineering design.

[0024] The present invention can intuitively reflect the actual stress and deformation conditions of each area of ​​the structural model by setting a deformation data acquisition module, an interference motion recognition module, a local feature simulation module, and a structural cloud map generation module, so that the color scale distribution of the cloud map is more reasonable, and the effect of the model's own force deformation is also amplified, avoiding the problem of color scale uniformity caused by overall large displacement.

[0025] Furthermore, the structural deformation cloud map obtained by the present invention can intuitively reflect the actual stress and deformation of each area of ​​the structural model, making the color scale distribution more reasonable, while also amplifying the effect of the model's own stress deformation, avoiding the problem of color scale uniformity caused by overall large displacement. Therefore, through the data processing and optimization method of the present invention, the data readability and utilization rate of the post-processing results are greatly improved, and engineers can more efficiently obtain and analyze the key features of model deformation, which helps to quickly locate structural problems, improve the quality and efficiency of engineering design analysis, and meet the needs of finite element analysis under complex working conditions. BRIEF DESCRIPTION OF THE DRAWINGS

[0026] Figure 1 A schematic diagram of a flow chart of a stress-bearing structure simulation post-processing method of the present invention; Figure 2 A cloud map generated by an existing post-processing method; Figure 3 A cloud map generated by applying the post-processing method of the present invention. DETAILED DESCRIPTION

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

[0028] On the contrary, the present invention covers any substitution, modification, equivalent method and scheme made on the essence and scope of the present invention as defined by the claims. Further, in order to make the public have a better understanding of the present invention, some specific details are described in detail in the following detailed description of the present invention. Those skilled in the art can fully understand the present invention without the description of these details.

[0029] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art to which the present invention pertains. The terms used herein are only for the purpose of describing specific embodiments and are not intended to limit the present invention. The term "or / and" as used herein includes any and all combinations of one or more of the associated listed items.

[0030] like Figure 1As shown, the first specific embodiment of the post-processing method of the force-bearing structure simulation based on numerical calculation of the present invention is: A post-processing method for simulation of a force-bearing structure based on numerical calculation comprises the following steps: Step 1: Obtain the force motion results of the structural model through the deformation data acquisition unit created in advance; Step 2: Using the previously created interference motion identification unit, based on the node state change, the force motion results are processed to obtain the overall motion data; Step 3: Using the local feature simulation unit created in advance, the overall motion data is removed from the force motion results to obtain the self-deformation data, which is used to characterize the stress distribution and deformation of the structural model itself; Step 4: Use the previously created structural cloud map generation unit to convert its own deformation data into color scale distribution information to obtain the structural deformation cloud map, thereby realizing the post-processing of the force-bearing structure simulation based on numerical calculation.

[0031] A second specific embodiment of the force-bearing structure simulation post-processing method based on numerical calculation of the present invention: A post-processing method for simulation of a force-bearing structure based on numerical calculation includes the following contents: In order to solve the problem of large overall displacement interfering with the display effect of cloud charts in finite element static analysis, the large overall displacement situations that may occur in the analysis can be classified and solved one by one. The results of static analysis need to be in a steady state, so the large overall displacement is mainly manifested in three situations: overall translation, rotation around a point, and rotation around a line. For these situations, it is necessary to classify and process them in combination with the complexity of the structural model and the assembly characteristics, which specifically include the following: Taking individuals as classification units, each component in the assembly structure model is considered separately. In the post-processing of the cloud map, by calculating the difference between the self-deformation of each body and its overall movement separately, the interference of the overall displacement on the result can be eliminated, thereby highlighting the self-deformation characteristics of the model. At the same time, it is also necessary to combine the classification of large displacements, such as overall translation, rotation around a point, and rotation around a line, and design corresponding processing algorithms respectively. For example, in the case of overall rotation around a point or around a line, it is necessary to introduce coordinate system transformation technology to separate the overall movement of the model and provide more intuitive deformation results for cloud map display.

[0032] In this embodiment, the displacement may be a translation, a rotation, or a movement.

[0033] In this embodiment, when the structural model has an overall translation, the initial stage of the calculation usually manifests itself as a consistent overall displacement of all nodes of a single model in a stress-free state. This situation usually occurs at the initial stage of loading, when the model individual is not in contact with the other individuals, and its degree of freedom allows overall translation. After the model individual contacts other individuals and establishes a constraint relationship, the node gradually begins to bear the load, the stress begins to distribute and shows local deformation. In the post-processing cloud map results at this stage, the overall translation often obscures the local deformation characteristics of the structural model, so a special algorithm mechanism is required to track the node displacement changes and eliminate the overall translation, thereby highlighting the stress distribution and deformation characteristics of the model itself.

[0034] From the perspective of numerical calculation, this embodiment separates the overall displacement and self-deformation data by improving the post-processing algorithm, and enhances the deformation feature effect of the model itself. By redefining the generation process of the cloud map, this embodiment can effectively highlight the local deformation characteristics of the model, enhance the color scale contrast of the cloud map, and provide engineers with more intuitive and accurate post-processing results. Therefore, the present invention can effectively solve the problem of poor cloud map display effect. By improving the numerical decomposition technology and dynamic color scale mapping algorithm, the present invention can automatically eliminate the interference of the overall displacement, only retain the data of the structure's own deformation, and generate a cloud map with higher color scale contrast.

[0035] A third specific embodiment of the force-bearing structure simulation post-processing method based on numerical calculation of the present invention: A post-processing method for simulation of a force-bearing structure based on numerical calculation includes the following contents: First, it is necessary to analyze the displacement changes of the structural model nodes in the result data at each time step and identify the critical moment when the stress begins to change.

[0036] After the structural model is divided by the finite element mesh, it will contain multiple nodes and units. Each node is a vertex in the mesh, and the calculated result data includes its node number, displacement and other data.

[0037] Stress is the equivalent stress, which is calculated from the principal stress. The calculation formula is as follows:

[0038] in, is the equivalent stress, , , are the three stress components of the principal stress, representing the stress components of the node in the three principal axis directions.

[0039] Since the displacement changes of each node are the same during the overall translation phase, it is only necessary to record a single overall displacement value without processing the data of each node separately. Assume that the positions of all nodes at the initial moment are , retrieve the stress of all time steps. When the stress changes at the nodes in the calculation result data, that is, when the stress value is greater than zero, record this time step The location is , at this time, its position change can be expressed as the overall displacement vector , and its calculation formula is as follows:

[0040]

[0041] in, is the node number, It is the total set of all node data.

[0042] When the stress changes at the node begin to occur, the node position changes from overall translation to and local deformation Together, when the final time step is t, its position is , then the expression of local deformation is as follows, and its calculation formula is as follows:

[0043] In order to highlight the local deformation, the overall translation should be eliminated from the cloud map results. ,reserve As the final display data.

[0044] In this embodiment, the specific processing flow is as follows: Step 1: Record the initial positions of all nodes , and initialize it.

[0045] Step 2: Perform overall displacement tracking and calculate the displacement changes of all nodes at each time step. , if there is no stress change, then record the overall translation vector .

[0046] Step 3: Perform stress change detection and monitor the stress value of each node ,when , marking the end moment of the overall translation.

[0047] Step 4: Perform data elimination and display, eliminate the overall displacement from the cloud map data, and display the local deformation results.

[0048] A fourth specific embodiment of the force-bearing structure simulation post-processing method based on numerical calculation of the present invention: A post-processing method for simulation of a force-bearing structure based on numerical calculation includes the following contents: Based on the node contact relationship, the overall translation of the structural model is analyzed. In the early stage of calculation, the phenomenon that all nodes of a single structural model are displaced as a whole actually reflects that there is no effective contact relationship between the structural model and other structural models. By analyzing the contact situation of the nodes in the structural model, the states of the nodes of real contact and pseudo contact (i.e., nodes that are not in direct contact but are considered to be in contact within the gap range set by the pre-processing) can be determined, so as to accurately distinguish the staged changes of overall displacement and local deformation.

[0049] Based on the monitoring of the node contact state, the displacement data is recorded and adjusted, and the cloud map display effect is finally optimized, so that the influence of the overall translation can be eliminated and the true characteristics of the local deformation and stress distribution are revealed. The specific steps include: The first step is to monitor the node contact status, which includes the following: During the calculation process, the contact status between each node in the structural model and other bodies is checked in real time. The contact status is divided into two categories: real contact (the node is clearly in contact with other bodies during calculation and meets the contact constraint conditions) and pseudo contact (the node is not directly in contact with other bodies, but within the gap range set in the finite element pre-processing stage, the default contact relationship).

[0050] The second step is to determine the contact status, which includes the following: For each node , the closest distance to other bodies is defined as ,if , then it is a real contact; if ,in is the gap threshold, then it is pseudo contact. When the node has neither real contact nor pseudo contact, the displacement change of the node mainly comes from the overall translation, and its expression is as follows:

[0051] The third step is to record and adjust the node displacement, which includes the following: In the overall translation stage, the displacement changes of each node are the same, so the overall displacement value is recorded , uniformly adjust the displacement changes of all nodes to eliminate the overall translation effect. The adjusted node displacement changes are:

[0052] The fourth step is to optimize the cloud map results, and use the node deformation data after removing the overall displacement to generate the cloud map, so as to highlight the local deformation and stress distribution characteristics. The specific steps include the following: Step 11: The cloud map rendering logic is reflected through the color gradation changes to reflect the size and distribution of the node data.

[0053] Step 12, normalize the node deformation data and map the deformation values ​​of all nodes to the standardized interval [0,1].

[0054] Step 13: The color gradient is from blue (low value) to green (transition color) to red (high value), showing the continuous change of deformation size.

[0055] Step 14, using linear interpolation to calculate the color value corresponding to the node, and then assigning the calculated color value to the corresponding node.

[0056] In this embodiment, the structural model nodes are divided into contact nodes and non-contact nodes by judging the contact state of the nodes. The displacement vector of the overall translation is recorded for the non-contact nodes. , and uniformly remove the overall displacement from the cloud map data, so that the stress distribution and deformation characteristics can be intuitively displayed. Compared with the traditional method, this embodiment judges the overall displacement based on the node contact relationship, which further improves the processing accuracy and result credibility.

[0057] A fifth specific embodiment of the force-bearing structure simulation post-processing method based on numerical calculation of the present invention: A post-processing method for simulation of a force-bearing structure based on numerical calculation includes the following contents: For the case where the structural model rotates around a point, the node movement at the initial stage of the calculation shows a trend of overall rotation around a fixed point. In this case, it is not applicable to judge simply by the change of node stress, because when a node in the target structural model comes into contact with other structural models, the contact stress will be transmitted to a large number of adjacent nodes in a short time through the connection relationship between the nodes, resulting in the stress values ​​of multiple nodes being non-zero. This stress transfer phenomenon will mask the true characteristics of the overall rotation around the point, which is easy to cause misjudgment and cannot accurately separate the overall rotation from the local deformation. Therefore, other judgment mechanisms must be used to eliminate the influence of rotation around the point and highlight the effective local deformation information.

[0058] In this case, the core of this embodiment is to determine the contact state between all nodes in the target structure model and the surface features of other structure models, including real contact and pseudo contact. By judging the contact relationship between nodes and surface features, it can be ensured that there is no real or pseudo contact relationship between the target structure model and other structure models, thereby eliminating the interference caused by node stress transfer. When there is no contact, the overall displacement of the target structure model can be obtained by rotating the center coordinate O and the rotation matrix Perform unified elimination and finally retain the local effective deformation displacement , and its calculation formula is as follows:

[0059] By eliminating the influence of the overall rotation around a point, the cloud map can accurately reflect the local deformation distribution and avoid the phenomenon that the overall displacement masks the local features. The specific judgment steps are as follows: The first step is to determine the contact relationship between nodes and surface features, which includes the following: For all node coordinates in the target structure model , it is necessary to compare with the surface feature set S of other structural models one by one to determine the minimum distance between the node and the surface feature :

[0060] in, The first Face features, Indicates the Euclidean distance between the node and the face feature, that is, the node With face features The minimum distance.

[0061] The method of obtaining surface features is as follows: By traversing the geometric feature data of the model, each surface in the model will have a corresponding face number, which can be used to determine (lock) the face. Then, according to the locked face, the corresponding face feature can be obtained.

[0062] Furthermore, according to the gap threshold Determine the contact status between the node and the surface feature, which includes the following: like , then the node With face features There is real contact.

[0063] like , although the node There is no direct contact, but the nodes are considered to be in contact with the faces due to the contact tolerance introduced by the finite element mesh discretization.

[0064] like , then the node No contact occurred.

[0065] In this embodiment, the gap threshold is 0.2 mm.

[0066] The second step is to determine the overall winding rotation, which includes the following: If all nodes in the target structure model If only one line feature has real contact or pseudo contact with other structural models, it can be determined that the target structural model is currently in a state of overall rotation around the fixed line. At this time, the overall node displacement It is mainly caused by the rotation transformation around the rotation center O. The influence of the overall rotation is eliminated by the following formula:

[0067] in, Is a node The total displacement vector, is a line feature, is the rotation angle around L The rotation matrix of is the initial position of the node, To eliminate the local deformation displacement after the whole body rotates around the point.

[0068] The total displacement vector is obtained by subtracting the three-dimensional coordinates of each node at the initial time 0 from the three-dimensional coordinates of each node after the solution is completed. The initial position of the node is the three-dimensional coordinate value of each node in the global coordinate system when the initial time is 0.

[0069] In this embodiment, the method for constructing the rotation matrix is ​​as follows: Unit vector establishing the line direction , which is expressed as follows:

[0070] in, The rotation axis is Directional weight.

[0071] According to the unit vector, the rotation matrix is ​​constructed, and its expression is as follows:

[0072] In this embodiment, the overall winding rotation judgment is mainly manifested in that all nodes in a single structural model in the initial stage of calculation revolve around a certain line feature of itself. This winding rotation will cause the overall rotation phenomenon of the node displacement. At the same time, due to the force transfer effect, the nodes near the line feature are easy to contact with other structural models, resulting in the contact stress being quickly transferred to the surrounding nodes, thereby causing stress value changes in the local area. In this case, it is unreliable to simply rely on the change of node stress to determine the overall winding phenomenon, because the real-time transmission of stress will mask the characteristics of the overall winding rotation and easily lead to misjudgment. Therefore, it is necessary to identify and eliminate the influence of the overall winding rotation by judging the contact relationship between the line features of the target structural model and the surface features of other structural models.

[0073] A specific embodiment of applying the present invention to simulate a mechanical assembly model: The post-processing method of the present invention is applied to simulate a mechanical assembly model, which includes the following contents: In a static analysis, a mechanical assembly model includes a large plate-like structure and a small block structure on it. The loading condition is that the block structure is subjected to a concentrated downward load, and the fixed boundaries are set at the edges of the plate. In the analysis, the block structure initially shifted downward as a whole due to insufficient freedom restrictions, with a displacement of 10 mm, while the small block structure only slightly deformed (less than 1 mm) after contacting the plate-like structure. In the post-processing results, the color gradation distribution of the cloud map is dominated by the huge value range of the overall displacement, resulting in the entire cloud map showing a single dark color. See Figure 2 The red block structure in the figure cannot directly reflect the actual deformation of the plate itself. This makes it difficult for engineers to locate the problem area through result analysis, significantly reducing the analysis efficiency.

[0074] The method of the present invention is used for the above-mentioned structural model. By analyzing the contact relationship between nodes and the deformation behavior of the block structure, the overall large displacement data is eliminated, and the focus is placed on the deformation characteristics of the block structure itself. In the result cloud map, the displacement range is optimized to its own deformation range, that is, the deformation value after eliminating invalid displacements, so that the deformation characteristics of the small block structure are clearly visible, and the deformation gradient color level in the center of the plate transitions from light to dark, accurately showing the stress distribution and deformation characteristics of different areas after loading. The optimized cloud map intuitively displays the local deformation characteristics, and engineers can quickly identify weak areas, which improves the efficiency of data analysis and the reliability of decision-making. Please refer to Figure 3 , Figure 3 The blue areas in the graph represent low values ​​in the self-deformation data, the green areas are transition areas, and the red areas represent high values ​​in the self-deformation data.

[0075] Therefore, the cloud map obtained by applying the present invention can avoid the interference of invalid displacement, so that the structural deformation results that are easy to observe with the naked eye can be intuitively displayed, which can help engineers identify areas with larger deformations and facilitate engineers to make accurate judgments.

[0076] An embodiment of a device applying the method of the present invention: An electronic device comprising: one or more processors; A storage device for storing one or more programs; When the one or more programs are executed by the one or more processors, the one or more processors implement the above-mentioned force-bearing structure simulation post-processing method based on numerical calculation.

[0077] A computer medium embodiment using the method of the present invention: A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the above-mentioned numerical calculation-based stress structure simulation post-processing method.

[0078] Those skilled in the art should understand that the embodiments of the present application can be provided as methods, systems, and computer program products. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, optical storage, etc.) containing computer-usable program codes.

[0079] The present application is described by flowcharts or / and block diagrams of the methods, devices (systems), and computer program products of the embodiments of the present application. It should be understood that each process or / and block in the flowchart or / and block diagram and the combination of the processes or / and blocks in the flowchart or / and block diagram can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart or / and block diagram. Figure 1 Process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0080] These computer program instructions may also be stored in a computer readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture including an instruction device, which implements the process Figure 1 Process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0081] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process in the computer or other programmable device. Figure 1 Process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0082] The unit in this application is an object that objectively describes the morphological structure with the help of physical or virtual expressions. The object is not equal to the object and is not limited to physical and virtual. It can be a data processing function, software program, processing mode, usage method, operation method, workflow, application process, electronic hardware, circuit module, processing system, system imitation or simulation object.

[0083] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the relevant field can still modify or replace the specific implementation methods of the present invention with equivalents. Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be included in the scope of protection of the claims of the present invention.

Claims

1. A post-processing method for simulation of a force-bearing structure based on numerical calculation, characterized in that: The following steps are involved: Step 1: Obtain the force motion results of the structural model through the deformation data acquisition unit created in advance; Step 2: Using the previously created interference motion identification unit, based on the node state change, the force motion results are processed to obtain the overall motion data; Step 3: Using the local feature simulation unit created in advance, the overall motion data is removed from the force motion results to obtain the self-deformation data, which is used to characterize the stress distribution and deformation of the structural model itself; Step 4: Use the previously created structural cloud map generation unit to convert its own deformation data into color scale distribution information to obtain the structural deformation cloud map, thereby realizing the post-processing of the force-bearing structure simulation based on numerical calculation.

2. A method for post-processing simulation of a force-bearing structure based on numerical calculation according to claim 1, characterized in that: Step 1: The method for obtaining the force motion results of the structural model through the deformation data acquisition unit created in advance is as follows: obtaining the concentrated load to be applied; Based on the structural model, set fixed boundary information; Divide the structural model into grids to obtain several grids; The vertices of the mesh are taken as nodes and numbered to obtain node information; Apply concentrated loads to the structural model, and perform static analysis on the nodes of the grid based on the fixed boundary information. Track the state of the nodes at each time step to obtain the state data of each node for monitoring the state changes of the nodes. The status data and node information of several nodes are aggregated to obtain the force and motion results of the structural model.

3. A method for post-processing simulation of a force-bearing structure based on numerical calculation as claimed in claim 2, characterized in that: Step 2: Using the previously created interference motion identification unit, based on the node state change, the force motion results are processed to obtain the overall motion data as follows: Classify the results of force motion and get the motion type. The types of motion include at least overall translation, rotation about a point or rotation about a line; Based on the type of movement, obtain the corresponding discriminant parameters, The discrimination parameters at least include node stress value and / or contact relationship value; Based on the force motion results, obtain the state data of the node at each time step; According to the motion type, the node stress value and / or the contact relationship value are selected, and the state data at each time step is judged to obtain the overall motion position at the end time of the overall motion; Subtract the overall motion position from the initial position of the corresponding node to obtain the overall motion data of one or more nodes.

4. A method for post-processing simulation of a force-bearing structure based on numerical calculation as claimed in claim 3, characterized in that: According to the motion type, the node stress value is selected, and the state data at each time step is judged to obtain the overall motion position at the end of the overall motion as follows: The motion type is overall translation, and the state data at each time step includes node stress values ​​and displacement values; The node stress value at each time step is judged. When the node stress value is greater than zero, it indicates that the overall movement is finished, and the time step at this time is recorded; According to the time step at this time, the displacement value of the corresponding node, that is, the overall motion position, is obtained.

5. A method for post-processing simulation of a force-bearing structure based on numerical calculation as claimed in claim 3, characterized in that: According to the motion type, the contact relationship value is selected, and the state data at each time step is judged to obtain the overall motion position at the end of the overall motion as follows: The motion type is rotation around a point or rotation around a line; the contact relationship value is zero or the gap threshold, which is used to characterize the contact state of the node, which is divided into real contact and pseudo contact; real contact means that the node is in clear contact with a neighboring individual within a certain time step and satisfies the contact constraint condition; pseudo contact means that the node does not directly contact the neighboring individual, but when the distance between the two is less than the gap threshold, it is assumed that there is a contact relationship; Get the state data at each time step, which includes the distance value between the node and the adjacent individuals; According to the distance value at each time step, the contact state between the node and the adjacent individuals is judged to obtain the contact judgment result, which includes no contact, pseudo contact or real contact; When the node changes from no contact to pseudo contact or real contact, it indicates that the overall motion ends, and the end time step at this time is recorded; According to the end time step, the position change value of the corresponding node, that is, the overall motion position, is calculated.

6. A method for post-processing simulation of a force-bearing structure based on numerical calculation according to claim 5, characterized in that: According to the end time step, the method to calculate the position change value of the corresponding node is as follows: Obtain the rotation center of the structural model and the three-dimensional coordinates of the nodes when the initial time is zero; According to the end time step, determine the real-time coordinate value of the node; Calculate the rotation difference between the real-time coordinates of the node and the rotation center, and establish the unit vector of the line direction; Construct a rotation matrix based on the unit vector; Multiply the rotation matrix by the rotation difference to get the position change value of the corresponding node.

7. A method for post-processing simulation of a force-bearing structure based on numerical calculation as claimed in claim 5, characterized in that: The method for judging the contact state between a node and an adjacent object is as follows: When the distance value is greater than the gap threshold, there is no contact, indicating that the nodes have neither real contact nor pseudo contact, and the position change of the nodes mainly comes from the overall movement; When the distance value is greater than zero and less than or equal to the gap threshold, it is a pseudo contact; When the distance value equals zero, it is a true contact.

8. A method for post-processing simulation of a force-bearing structure based on numerical calculation as claimed in claim 1, characterized in that: Step 4: Use the previously created structural cloud map generation unit to convert the deformation data into color distribution information. The method for obtaining the structural deformation cloud map is as follows: Get the self-deformation data of all nodes; Map the deformation data to the standardized interval [0,1] to obtain the deformation standard data; Set a color gradient, which includes at least a first color, a second color, and a third color, and is used to express continuous changes in the deformation size; the first color is used to represent the low value in the deformation data itself, the second color is a transition color, and the third color is used to represent the high value in the deformation data itself; According to the color gradient, the deformed standard data is processed using the linear interpolation method to obtain the corresponding color value; Assign color values ​​to corresponding nodes to obtain color scale distribution information; According to the color distribution information, the structural deformation cloud map is drawn.

9. A force structure simulation post-processing system based on numerical calculation, characterized in that: It includes a deformation data acquisition module, an interference motion recognition module, a local feature simulation module and a structure cloud map generation module; Deformation data acquisition module, used to obtain the force and motion results of the structural model; Interference motion identification module, used to process the force motion results based on the node state changes to obtain the overall motion data; The local feature simulation module is used to remove the overall motion data from the force motion results and obtain the self-deformation data to characterize the stress distribution and deformation of the structural model itself; The structural cloud map generation module is used to convert the self-deformation data into color-scale distribution information to obtain the structural deformation cloud map.

10. An electronic device, characterized in that: It includes: one or more processors; A storage device for storing one or more programs; When the one or more programs are executed by the one or more processors, the one or more processors implement a post-processing method for simulation of a force-bearing structure based on numerical calculation as described in any one of claims 1 to 8.

Images