Twin power station collision visual simulation method and system based on physical simulation

Abstract

The invention provides a twin power station collision visual simulation method and system based on physical simulation, and relates to the technical field of data processing, and the method comprises the steps: obtaining a three-dimensional geometric model of a twin power station, carrying out space division, extracting geometric parameters, obtaining material data, constructing a material correlation model and a mechanical model, and carrying out fusion to obtain a physical simulation model; acquiring initial collision condition parameters as boundary conditions, performing local grid encryption, constructing a nonlinear kinetic equation of a collision process, performing time domain solution, determining equivalent stress distribution and equivalent plastic strain distribution, and generating a collision result; judging whether plastic deformation or fracture failure occurs, generating full collision process physical field data, generating a collision animation frame sequence, generating a stress distribution cloud picture and a plastic strain cloud picture, constructing a fracture special effect, and comprehensively generating a collision visual scene.

Description

Technical Field

[0001] The present invention relates to the technical field of data processing, and in particular to a twin power station collision visualization simulation method and system based on physical simulation. Background Art

[0002] As an important infrastructure for national economic and social development, the safe and stable operation of power stations is of vital importance. During the construction and operation of power stations, component collision accidents may occur due to natural disasters, human errors and other factors, which seriously threaten the safety and reliability of power stations. Traditional power station component collision analysis mainly relies on on-site investigation and experience judgment, which has problems such as low analysis accuracy and poor evaluation efficiency, and it is difficult to meet the increasing safety management needs of power stations. In the existing technology, physical simulation is disconnected from visualization, and the dynamic display of the entire collision process cannot be achieved. At the same time, there is a lack of high-fidelity digital twin models based on actual power stations, making it difficult to carry out accurate simulation analysis. Therefore, a solution is urgently needed to solve the problems existing in the existing technology. Summary of the invention

[0003] The embodiments of the present invention provide a twin power station collision visualization simulation method and system based on physical simulation, which can at least solve some of the problems existing in the prior art.

[0004] A first aspect of an embodiment of the present invention provides a twin power station collision visualization simulation method based on physical simulation, comprising: Acquire a three-dimensional geometric model of the twin power station, perform spatial division of the twin power station by using an octree space division algorithm, generate a voxelized model, and extract geometric parameters of power station components in the twin power station based on the topological relationship of the voxel units, acquire material data of the power station components based on the physical design specification and bill of materials of the twin power station, associate the geometric parameters with the material data, and construct a material association model, acquire a mechanical model of the internal frame structure of the physical power station according to the structural design drawing of the physical power station, calculate the force distribution in combination with mechanical principles, obtain a mechanical simulation model, and fuse it with the material association model to obtain a physical simulation model; Add the physical simulation model to a preset physical simulation engine to generate a basic collision simulation model, obtain the collision initial condition parameters of the two power station components that collide and add them to the basic collision simulation model as boundary conditions, perform local mesh encryption on the power station components through an adaptive mesh subdivision algorithm, define the material constitutive model through a preset finite element model in combination with the material data corresponding to the power station components, construct a nonlinear dynamic equation of the collision process based on the nonlinear properties of the power station components, solve the nonlinear dynamic equation in the time domain in combination with an explicit central difference algorithm, generate physical field information, calculate the strain field and stress field of the power station components through a conservative node average stress method, determine the equivalent stress distribution and the equivalent plastic strain distribution, and generate a collision result based on a yield criterion and a fracture failure criterion; Based on the collision result, it is judged whether the power station component has plastic deformation or fracture failure. If plastic deformation occurs, the finite element model is updated through a node displacement correction algorithm. If fracture failure occurs, the failed unit is removed from the finite element model through a unit deletion method and a free surface is generated. The simulation is repeated until the collision ends, and the physical field data of the entire collision process is generated. The displacement field data in the physical field data of the entire collision process is mapped to the three-dimensional geometric model. The position and shape of the power station component are updated through a vertex displacement difference algorithm to generate a collision animation frame sequence. Based on the equivalent stress distribution and the equivalent plastic strain distribution, a stress distribution cloud map and a plastic strain cloud map are generated in combination with a volume rendering algorithm. For the power station component that has fracture failure, according to the position and normal vector of the fracture unit, a mesh crushing algorithm is used to generate a fracture texture and a fracture morphology. A fracture special effect is constructed in combination with splashing fragments. A collision visualization scene is generated by integrating the collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effect.

[0005] In an optional embodiment, The three-dimensional geometric model of the twin power station is obtained, the twin power station is spatially divided by an octree space partitioning algorithm, a voxelized model is generated, and the geometric parameters of the power station components in the twin power station are extracted based on the topological relationship of the voxel units. The material data of the power station components are obtained based on the physical design specification and the bill of materials of the twin power station, the geometric parameters are associated with the material data, and a material association model is constructed. The mechanical model of the internal frame structure of the physical power station is obtained according to the structural design drawing of the physical power station, the force distribution is calculated in combination with the mechanical principle, and a mechanical simulation model is obtained and merged with the material association model to obtain a physical simulation model including: The point cloud data corresponding to the twin power station is acquired by three-dimensional laser scanning technology, and the three-dimensional geometric model corresponding to the twin power station is obtained by point cloud data preprocessing and registration and splicing operations. Based on the octree space partitioning algorithm, the root node of the octree is determined according to the bounding box of the twin power station, and the complexity metric is determined at each node according to the number of point clouds. The complexity threshold is combined to make a judgment. If the complexity of the current node is greater than the complexity threshold, the current node is recursively divided into eight equal parts, and the judgment is repeated until the complexity of all nodes is less than the complexity threshold to obtain a leaf node. For each leaf node, a corresponding voxel unit is generated and combined to obtain the voxelized model; Based on the topological relationship corresponding to the voxel units, different voxel sets corresponding to the power station components in the twin power station are identified, and the outer envelope box of the component is calculated based on the voxel set belonging to the same component to obtain the size information of the current component, and the shape feature of the current component is extracted by fitting the surface mesh of the voxel set, and the geometric center of the outer envelope box is used as the spatial position coordinate of the current component. The size information, shape feature and spatial position coordinate are combined to obtain the geometric parameters of the current component, and each component in the twin power station is traversed to obtain the corresponding geometric parameters; Based on the physical design specifications and bill of materials of the twin power stations, the material composition of the power station components in each power station is obtained. For metal materials, a material property database is constructed in combination with physical and mechanical property parameters, and the property parameters of the material composition of each component are retrieved from the material property database. For non-metallic materials and new materials, mechanical performance parameters are determined in combination with simulation experimental tests, and the material data is obtained by combining the property parameters and the mechanical performance parameters. The bill of materials is traversed and a first relationship table and a second relationship table corresponding to the power station components and the geometric parameters and the material data are respectively generated, the first relationship table and the second relationship table are connected, and material property parameters or mechanical performance parameters are retrieved, and the material association model is obtained by combining them; According to the structural design drawing of the physical power station, a mechanical model of the internal frame structure of the physical power station is obtained, and the key load mode is determined by combining mechanical principles and load condition analysis technology, wherein the key load mode is divided into static load and dynamic load. For the static load, the static force distribution is calculated by a finite element static analysis method, and for the dynamic load, the vibration mode and frequency response characteristics of the internal frame structure are obtained by a finite element dynamic analysis method. Based on the vibration mode and frequency response characteristics, combined with a transient mechanical analysis method, the dynamic force distribution is obtained by solving the dynamic equation by a direct integration method. The static force distribution and the dynamic force distribution are combined to obtain a mechanical simulation model, and the static force distribution and the dynamic force distribution in the mechanical simulation model are mapped to the material association model to obtain the physical simulation model.

[0006] In an optional embodiment, Based on the vibration mode and frequency response characteristics, combined with the transient mechanical analysis method, the dynamic force distribution is obtained by solving the dynamic equation through direct integration method as shown in the following formula: ; in, F ( t n+1 ) indicates that at the time step n +1 external force load, [M] represents the mass matrix, { u n+1} represents the time step n +1 displacement vector, [C] represents the damping matrix, Represents the time step n The velocity vector, △ t represents the time increment, Represents the time step n The acceleration vector, Represents the time step n +1 acceleration vector, [K] represents the stiffness matrix.

[0007] In an optional embodiment, The physical simulation model is added to a preset physical simulation engine to generate a basic collision simulation model, the collision initial condition parameters of the two power station components that collide are obtained and added to the basic collision simulation model as boundary conditions, the power station components are locally meshed by an adaptive mesh subdivision algorithm, the material constitutive model is defined by a preset finite element model in combination with the material data corresponding to the power station components, the nonlinear dynamic equation of the collision process is constructed based on the nonlinear properties of the power station components, the nonlinear dynamic equation is solved in the time domain in combination with an explicit central difference algorithm, the physical field information is generated, and the strain field and stress field of the power station components are calculated by the conservative node average stress method, the equivalent stress distribution and the equivalent plastic strain distribution are determined, and the collision results are generated based on the yield criterion and the fracture failure criterion, including: Building a physical simulation engine based on an explicit finite element algorithm, adding the physical simulation model to the physical simulation engine and generating the basic collision simulation model; Perform kinematic analysis on the two power station components that collided, track the trajectories and determine the center of mass position based on real-time monitoring data, determine the angular velocity and contact position at the time of collision in combination with the corresponding geometric shape and posture of each component, combine the center of mass position, angular velocity and contact position to obtain the collision initial condition parameters and add them to the basic collision simulation model as boundary conditions; For the power station components that collide, an initial grid is set by an adaptive grid subdivision algorithm, a stress analysis is performed on the initial grid, the equivalent stress of each unit is determined, and the units whose equivalent stress is greater than a preset grid stress threshold are marked, the marked areas are used as areas that need to be subdivided, and the areas that need to be subdivided are recursively subdivided by an adaptive space partitioning algorithm until the preset grid density requirement is met, and the subdivided grids are smoothed and coupled with the initial grid to obtain a second finite element model; Based on the second finite element model, in combination with the material data corresponding to the power plant component, through the pre-set finite element model, for isotropic materials, a bilinear elastoplastic constitutive model is constructed based on the elastic modulus and yield stress, and for anisotropic materials, a yield function is determined based on the anisotropic parameters corresponding to the current material, and a nonlinear dynamic equation corresponding to the collision process is constructed in combination with the nonlinear properties of the power plant component, wherein the nonlinear properties include large deformation, material nonlinearity and contact nonlinearity; The nonlinear dynamic equations are solved based on an explicit central difference algorithm, the stress, internal force and contact state of each component are updated, physical field information is generated, and the strain field and stress field of the power station components are calculated by the conservation node average stress method. The equivalent stress distribution and the equivalent plastic strain distribution are determined based on the strain field and the stress field. According to the material properties and current state of the power station components, combined with the corresponding yield criterion and fracture failure standard, the collision results are generated and output.

[0008] In an optional embodiment, For anisotropic materials, the yield function is determined based on the anisotropic parameters corresponding to the current material as shown in the following formula: ; in, F Indicates that the material y Axis and z Differential resistance to yield between shafts, G Indicates that the material z Axis and x Differential resistance to yield between shafts, H Indicates that the material x Axis and y Differential resistance to yield between shafts, s 11 Indicates that the material x The normal stress component in the axial direction, s 22 Indicates that the material y The normal stress component in the axial direction, s 33 Indicates that the material z The normal stress component in the axial direction,m represents the Horsford index, s y Represents the yield stress.

[0009] In an optional embodiment, Based on the collision result, it is judged whether the power station component has plastic deformation or fracture failure. If plastic deformation occurs, the finite element model is updated by a node displacement correction algorithm. If fracture failure occurs, the failed unit is removed from the finite element model by a unit deletion method and a free surface is generated. The simulation is repeated until the collision ends, and the physical field data of the entire collision process is generated. The displacement field data in the physical field data of the entire collision process is mapped to the three-dimensional geometric model. The position and shape of the power station component are updated by a vertex displacement difference algorithm, and a collision animation frame sequence is generated. Based on the equivalent stress distribution and the equivalent plastic strain distribution, a stress distribution cloud map and a plastic strain cloud map are generated in combination with a volume rendering algorithm. For the power station component that has fracture failure, according to the position and normal vector of the fracture unit, a mesh crushing algorithm is used to generate a fracture texture and a fracture morphology, and a fracture special effect is constructed in combination with splashing fragments. The collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effect are integrated to generate a collision visualization scene, including: Based on the collision result, determine whether the power station component has plastic deformation or fracture failure. If plastic deformation occurs, extract the displacement field data of the current component based on the node displacement correction algorithm. For each node, update the node coordinates according to the corresponding displacement vector, generate a finite element mesh based on the updated node coordinates, and update the material properties and plastic strain of the current component. Repeat the update until a finite element model sequence corresponding to the current component is obtained. If the collision result indicates that the current power plant component has fractured and failed, then based on the fracture failure criterion, the failed unit is identified and a failed unit set is generated, for each failed unit, the corresponding node set is extracted and the node set is removed, and nodes are inserted at the boundary of the fracture failure area and connected with the finite element mesh before the fracture to obtain the free surface, and the finite element mesh is regenerated based on the free surface and the material properties and plastic strain of the current component are updated, and the update is repeated until the finite element model sequence corresponding to the current component is obtained; Extracting the physical field data of the entire collision process based on the finite element model sequence, mapping the displacement field in the physical field data of the entire collision process to the three-dimensional geometric model to generate a vertex displacement field, updating the vertex coordinates of the three-dimensional geometric model through a vertex displacement difference algorithm and calculating the corresponding normal vector and texture coordinates, and updating the position and shape of the power station component in combination with the physical field data of the entire collision process to obtain the collision animation frame sequence; For the power station structure, based on the collision animation frame sequence and the equivalent stress distribution and equivalent plastic strain distribution, a color mapping table for each voxel is determined by a volume rendering algorithm, and a search is performed in the color mapping table based on the stress value and the strain value. A ray is emitted to each pixel point by a ray casting algorithm, and the intersection with different voxels is determined, and the color corresponding to each pixel point is determined, and the operation is repeated to obtain the stress distribution cloud map and the plastic strain cloud map; For power station components that have failed due to fracture, a random crack distribution is generated on the boundary surface corresponding to the fracture unit according to the position and normal vector of the fracture unit through a mesh breaking algorithm according to the fracture direction and the material data, a polygonal mesh is generated with each crack line in the random crack distribution as a starting point, the polygonal mesh is triangulated and a normal map and a height map are attached to generate an initial fracture morphology and random disturbances are added to obtain fracture texture and fracture morphology, fragment particles are randomly generated in the fracture unit, the kinematic parameters of each fragment particle are initialized, and the scattering trajectory corresponding to each fragment particle is determined by explicit integration, a fracture special effect is constructed based on the scattering trajectory and the fragment particles, and a collision visualization scene is generated by combining the collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effect.

[0010] In an optional embodiment, The ray casting algorithm is used to emit light to each pixel and determine the intersection points with different voxels. The color corresponding to each pixel is determined as shown in the following formula: ; in, CP ( u , v ) represents pixel points ( u , v ), CV ( x , y , z ) represents a voxel ( x , y , z ), S ( x , y , z ) represents a voxel ( x , y , z ) scattering coefficient, t ( x , y , z ) represents a voxel ( x , y , z )’s absorption coefficient, lis the path of light propagation, indicating the propagation distance of light in the voxel, dλ is the light element, used to calculate the path of the light l The points on .

[0011] A second aspect of an embodiment of the present invention provides a twin power station collision visualization simulation system based on physical simulation, including: The first unit is used to obtain a three-dimensional geometric model of a twin power station, perform spatial division of the twin power station by an octree space division algorithm, generate a voxelized model, and extract geometric parameters of power station components in the twin power station based on the topological relationship of the voxel units, obtain material data of the power station components based on the physical design specification and the bill of materials of the twin power station, associate the geometric parameters with the material data, construct a material association model, obtain a mechanical model of the internal frame structure of the physical power station according to the structural design drawing of the physical power station, calculate the force distribution in combination with the mechanical principle, obtain a mechanical simulation model, and fuse it with the material association model to obtain a physical simulation model; The second unit is used to add the physical simulation model to a preset physical simulation engine to generate a basic collision simulation model, obtain the collision initial condition parameters of the two power station components that collide and add them to the basic collision simulation model as boundary conditions, perform local mesh encryption on the power station components through an adaptive mesh subdivision algorithm, define the material constitutive model through a preset finite element model in combination with the material data corresponding to the power station components, construct a nonlinear dynamic equation of the collision process based on the nonlinear properties of the power station components, solve the nonlinear dynamic equation in the time domain in combination with an explicit central difference algorithm, generate physical field information, calculate the strain field and stress field of the power station components through a conservative node average stress method, determine the equivalent stress distribution and the equivalent plastic strain distribution, and generate a collision result based on the yield criterion and the fracture failure criterion; The third unit is used to judge whether the power station component has plastic deformation or fracture failure based on the collision result. If plastic deformation occurs, the finite element model is updated through a node displacement correction algorithm. If fracture failure occurs, the failed unit is removed from the finite element model through a unit deletion method and a free surface is generated. The simulation is repeated until the collision ends, and the physical field data of the entire collision process is generated. The displacement field data in the physical field data of the entire collision process is mapped to the three-dimensional geometric model, and the position and shape of the power station component are updated through a vertex displacement difference algorithm to generate a collision animation frame sequence. Based on the equivalent stress distribution and the equivalent plastic strain distribution, a stress distribution cloud map and a plastic strain cloud map are generated in combination with a volume rendering algorithm. For the power station component that has fracture failure, according to the position and normal vector of the fracture unit, a mesh crushing algorithm is used to generate a fracture texture and a fracture morphology, and a fracture special effect is constructed in combination with splashing fragments. The collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effect are integrated to generate a collision visualization scene.

[0012] According to a third aspect of the embodiments of the present invention, An electronic device is provided, comprising: processor; a memory for storing processor-executable instructions; The processor is configured to call the instructions stored in the memory to execute the aforementioned method.

[0013] According to a fourth aspect of the embodiments of the present invention, A computer-readable storage medium is provided, on which computer program instructions are stored. When the computer program instructions are executed by a processor, the aforementioned method is implemented.

[0014] In the present invention, by acquiring the three-dimensional geometric model and material data of the twin power station and combining it with the mechanical model of the physical power station, a sophisticated physical simulation model is constructed. Adaptive mesh subdivision and nonlinear dynamic solution techniques are used to accurately simulate the complex mechanical behaviors of components such as deformation and fracture, and high-precision physical field data of the entire collision process are obtained, which provides a reliable basis for accurately evaluating the impact of collision. Based on the physical field data obtained from the collision simulation, through visualization algorithms such as displacement mapping, volume rendering, and mesh crushing, various visualization expressions such as animation, cloud map, and fracture special effects of the component collision process are generated, which comprehensively, three-dimensionally, and realistically show the deformation and fracture of components under collision loads. , crushing and other complex mechanical behaviors, providing an immersive visualization experience, intuitively revealing the collision mechanism, and through the seamless connection of geometric modeling, physical simulation, visual drawing and other links, realizing the full process closed loop from real scenes to virtual simulation and then to visual results, greatly improving the pertinence, accuracy and persuasiveness of collision analysis. In summary, the present invention realizes the refined simulation and realistic display of the collision process through the deep integration of physical simulation and visualization technology, establishes a multidisciplinary collision analysis framework, greatly improves the efficiency of power plant safety assessment and accident investigation, provides strong technical support for the full life cycle management of power plants, and has good versatility and expansibility. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] Figure 1 It is a flow chart of a twin power station collision visualization simulation method based on physical simulation according to an embodiment of the present invention; Figure 2 It is a structural schematic diagram of a twin power station collision visualization simulation system based on physical simulation according to an embodiment of the present invention. DETAILED DESCRIPTION

[0016] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0017] The technical solution of the present invention is described in detail with specific embodiments below. The following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described in detail in some embodiments.

[0018] Figure 1 FIG. 1 is a flow chart of a twin power station collision visualization simulation method based on physical simulation according to an embodiment of the present invention. Figure 1 As shown, the method includes: S1. Obtain a three-dimensional geometric model of the twin power station, perform spatial division of the twin power station by using an octree space division algorithm, generate a voxelized model, and extract geometric parameters of power station components in the twin power station based on the topological relationship of the voxel units, obtain material data of the power station components based on the physical design specification and bill of materials of the twin power station, associate the geometric parameters with the material data, and construct a material association model, obtain a mechanical model of the internal frame structure of the physical power station according to the structural design drawing of the physical power station, calculate the force distribution in combination with mechanical principles, obtain a mechanical simulation model, and merge it with the material association model to obtain a physical simulation model; The octree is a data structure for segmenting a three-dimensional space, which recursively divides the space into eight sub-areas to form a tree structure. The voxel is the abbreviation of "volume pixel" and is the smallest unit in a three-dimensional space. The voxelized model is the process of decomposing a three-dimensional model into a voxel grid, and each voxel represents a volume part of the model. The voxel unit is the basic unit in the voxelized model, similar to a pixel in a two-dimensional image, and each voxel unit represents a small cube in the space. The solid design specification refers to the standards and rules used to create and describe three-dimensional models. The material association model mechanical model takes into account the physical properties of the material, such as elastic modulus, Poisson's ratio, etc., and is used to analyze the mechanical behavior of the material under different conditions. The mechanical principle is the basic theory for studying the behavior and design of mechanical systems, including mechanics, dynamics, kinematics, etc. The mechanical simulation model is a computational model used to simulate and analyze the behavior of objects under various loads and conditions.

[0019] In an optional embodiment, The three-dimensional geometric model of the twin power station is obtained, the twin power station is spatially divided by an octree space partitioning algorithm, a voxelized model is generated, and the geometric parameters of the power station components in the twin power station are extracted based on the topological relationship of the voxel units. The material data of the power station components are obtained based on the physical design specification and the bill of materials of the twin power station, the geometric parameters are associated with the material data, and a material association model is constructed. The mechanical model of the internal frame structure of the physical power station is obtained according to the structural design drawing of the physical power station, the force distribution is calculated in combination with the mechanical principle, and a mechanical simulation model is obtained and merged with the material association model to obtain a physical simulation model including: The point cloud data corresponding to the twin power station is acquired by three-dimensional laser scanning technology, and the three-dimensional geometric model corresponding to the twin power station is obtained by point cloud data preprocessing and registration and splicing operations. Based on the octree space partitioning algorithm, the root node of the octree is determined according to the bounding box of the twin power station, and the complexity metric is determined at each node according to the number of point clouds. The complexity threshold is combined to make a judgment. If the complexity of the current node is greater than the complexity threshold, the current node is recursively divided into eight equal parts, and the judgment is repeated until the complexity of all nodes is less than the complexity threshold to obtain a leaf node. For each leaf node, a corresponding voxel unit is generated and combined to obtain the voxelized model; Based on the topological relationship corresponding to the voxel units, different voxel sets corresponding to the power station components in the twin power station are identified, and the outer envelope box of the component is calculated based on the voxel set belonging to the same component to obtain the size information of the current component, and the shape feature of the current component is extracted by fitting the surface mesh of the voxel set, and the geometric center of the outer envelope box is used as the spatial position coordinate of the current component. The size information, shape feature and spatial position coordinate are combined to obtain the geometric parameters of the current component, and each component in the twin power station is traversed to obtain the corresponding geometric parameters; Based on the physical design specifications and bill of materials of the twin power stations, the material composition of the power station components in each power station is obtained. For metal materials, a material property database is constructed in combination with physical and mechanical property parameters, and the property parameters of the material composition of each component are retrieved from the material property database. For non-metallic materials and new materials, mechanical performance parameters are determined in combination with simulation experimental tests, and the material data is obtained by combining the property parameters and the mechanical performance parameters. The bill of materials is traversed and a first relationship table and a second relationship table corresponding to the power station components and the geometric parameters and the material data are respectively generated, the first relationship table and the second relationship table are connected, and material property parameters or mechanical performance parameters are retrieved, and the material association model is obtained by combining them; According to the structural design drawing of the physical power station, a mechanical model of the internal frame structure of the physical power station is obtained, and the key load mode is determined by combining mechanical principles and load condition analysis technology, wherein the key load mode is divided into static load and dynamic load. For the static load, the static force distribution is calculated by a finite element static analysis method, and for the dynamic load, the vibration mode and frequency response characteristics of the internal frame structure are obtained by a finite element dynamic analysis method. Based on the vibration mode and frequency response characteristics, combined with a transient mechanical analysis method, the dynamic force distribution is obtained by solving the dynamic equation by a direct integration method. The static force distribution and the dynamic force distribution are combined to obtain a mechanical simulation model, and the static force distribution and the dynamic force distribution in the mechanical simulation model are mapped to the material association model to obtain the physical simulation model.

[0020] The point cloud data preprocessing includes operations such as cleaning, denoising, and downsampling of the point cloud data to improve data quality and analysis accuracy. The point cloud registration is the process of aligning multiple point cloud data sets to the same coordinate system. The bounding box is a minimum rectangle or cube that is used to enclose and represent the spatial range of a three-dimensional object. The complexity metric is used to evaluate the complexity of a model or algorithm, including time complexity and space complexity. The outer envelope is a minimum convex polyhedron that contains all point cloud data. The surface mesh is a mesh composed of polygons (such as triangles) that is used to represent the surface of a three-dimensional object. The physical and mechanical property parameters describe the physical and mechanical properties of a material or structure, such as density. , elastic modulus, hardness, etc., the load condition analysis technology is used to analyze and predict the behavior and performance of the structure under different load conditions, the finite element static analysis method is used to calculate and analyze the stress, strain and displacement of the structure under static load, the finite element dynamic analysis method is used to analyze the behavior of the structure under dynamic load, including vibration and impact, the dynamic mode is the characteristic that describes the different deformation modes of the structure when vibrating, the frequency response characteristic describes the response behavior of the structure at different frequencies, the transient mechanics analysis method is used to analyze and predict the response of the structure under transient loads (such as impact), and the direct integration method is a numerical method used to solve the equations of the dynamic system and directly calculate the response of the system.

[0021] Use a 3D laser scanner to conduct a comprehensive scan of the twin power stations to obtain point cloud data. Preprocess the obtained point cloud data, including denoising, filtering, and downsampling. Use a point cloud registration algorithm (such as the ICP algorithm) to stitch point cloud data from different perspectives to obtain a complete 3D point cloud model. Use the octree space partitioning algorithm to voxelize the point cloud model, determine the bounding box of the point cloud model as the root node of the octree, calculate the complexity measure of each node based on the number of point clouds, and compare it with the preset complexity threshold. If it is greater than the threshold, divide the current node into eight equal parts, and recursively divide the child nodes until the complexity of all nodes is less than the threshold. Generate voxel units corresponding to leaf nodes and combine them to obtain a complete voxelized model. Based on the topological relationship of voxel units in the voxelized model, the voxel sets corresponding to different power station components are identified. For each power station component, the outer envelope of the component voxel set is calculated to obtain the size information of the component, and the surface mesh of the component voxel set is fitted to extract the shape features of the component. The geometric center of the component outer envelope is used as the spatial position coordinate of the component. The size information, shape features and spatial position coordinates are integrated to obtain the geometric parameters of the component. All components in the twin power stations are traversed to obtain their respective geometric parameters. Obtain the physical design specifications and bill of materials of the twin power station, determine the material composition of each power station component, for metal materials, build a material property database in combination with physical and mechanical property parameters, retrieve the property parameters of each component material composition in the material property database, for non-metallic materials and new materials, determine the mechanical performance parameters of the materials through simulation experiments, combine the property parameters and mechanical performance parameters to obtain complete material data, traverse the bill of materials, generate a first relationship table between power station components and geometric parameters, and a second relationship table between power station components and material data, connect the first relationship table and the second relationship table, retrieve material property parameters or mechanical performance parameters, and combine to obtain a material association model; According to the structural design drawing of the physical power station, the mechanical model of the internal frame structure of the physical power station is obtained. In combination with mechanical principles and load condition analysis technology, the key load mode is determined, which is divided into static load and dynamic load. For static load, the finite element static analysis method is used to calculate the static force distribution. For dynamic load, the finite element dynamic analysis method is used to obtain the vibration mode and frequency response characteristics of the internal frame structure. Based on the vibration mode and frequency response characteristics, combined with the transient mechanical analysis method, the dynamic equation is solved using the direct integration method to obtain the dynamic force distribution. The static force distribution and the dynamic force distribution are combined to obtain the mechanical simulation model. The static force distribution and the dynamic force distribution in the mechanical simulation model are mapped to the material association model to obtain the final physical simulation model. For example, taking a simplified twin power station as an example, assuming that it consists of three main components: a main power plant, a cooling tower and a transmission line, a 3D laser scanner is used to comprehensively scan the power station to obtain point cloud data. After preprocessing, the point cloud is spliced ​​into a complete 3D model using the ICP algorithm, and the point cloud model is voxelized using the octree algorithm to obtain a voxelized model. Based on the voxelized model, the voxel sets corresponding to the main power plant, cooling tower and transmission line are identified. For each component, its outer envelope box size is calculated, the surface mesh is fitted to extract shape features, and the geometric center is used as the spatial position coordinate. Finally, the geometric parameters of the three components are obtained, and the design specifications and material list of the power station are obtained. Assume Assume that the main plant is made of reinforced concrete, the cooling tower is made of steel structure and concrete, and the transmission line is made of steel and aluminum alloy. Retrieve the property parameters of these materials from the material property database, and determine their mechanical performance parameters through simulation experiments. Generate a relationship table between components and geometric parameters and material data, and combine them to obtain a material association model. According to the structural design drawing of the power station, obtain the mechanical model of its internal frame structure, determine the static load (such as dead weight) and dynamic load (such as wind load), use the finite element method to calculate the static force distribution and dynamic force distribution respectively, and obtain a mechanical simulation model. Map the force distribution in the mechanical simulation model to the material association model to obtain a physical simulation model.

[0022] In this embodiment, the octree space partitioning algorithm can adaptively generate a voxelized model according to the complexity of the model, which can not only ensure geometric details but also improve data processing efficiency. By calculating the outer envelope box, fitting the surface mesh and other operations, it can automatically extract the geometric parameters such as the size, shape, and position of the component, greatly reducing the workload of manual modeling, and improving the efficiency and accuracy of parametric modeling. A relationship table between the component and the geometric and material data is generated, and connected and combined, which can achieve seamless integration of data in different fields and formats, provide consistent data support, map the mechanical simulation results with the material association model, and achieve deep coupling of multiple physical fields (geometry, materials, mechanics), greatly improving the accuracy and reliability of the simulation. In summary, this embodiment constructs a high-fidelity, multi-dimensional, full-life cycle power station digital twin model through systematic data acquisition, processing, integration and analysis, which provides powerful digital tools and decision-making support for various stages such as design optimization, construction management, operation and maintenance of the power station, and is expected to greatly improve the efficiency, quality and safety of the power station project, and promote the development of the power station towards intelligence and digitalization.

[0023] In an optional embodiment, Based on the vibration mode and frequency response characteristics, combined with the transient mechanical analysis method, the dynamic force distribution is obtained by solving the dynamic equation through direct integration method as shown in the following formula: ; in, F ( t n+1 ) indicates that at the time step n +1 external force load, [M] represents the mass matrix, { u n+1} represents the time step n +1 displacement vector, [C] represents the damping matrix, Represents the time step n The velocity vector, △ t represents the time increment, Represents the time step n The acceleration vector, Represents the time step n +1 acceleration vector, [K] represents the stiffness matrix.

[0024] In this embodiment, by reasonably setting the mass matrix, damping matrix and stiffness matrix, the dynamic behavior of the structure can be accurately described, and the accuracy of dynamic response calculation can be improved. By discretizing the time domain, the continuous dynamic problem is converted into a discrete algebraic equation group, which has the characteristics of stable calculation and high efficiency, and is suitable for large-scale, nonlinear dynamic problems. By applying the corresponding external force load at each time step, the time-varying characteristics of the dynamic load, such as wind load, seismic load, etc., can be simulated. By solving the dynamic equation, the displacement vector, velocity vector and acceleration vector of the structure at each time step can be obtained, which can be used to evaluate the deformation, vibration and stress state of the structure under the action of the dynamic load, and provide a basis for the analysis of the dynamic performance of the structure such as earthquake resistance and wind resistance. In summary, this embodiment fully considers the dynamic characteristics of the structure and the time-varying nature of the external force load, and obtains the dynamic response of the structure through numerical solution, which provides important data support for subsequent dynamic performance analysis, fatigue assessment, etc., and realizes the comprehensive consideration of static and dynamic loads through coupling with static analysis, which improves the accuracy of the simulation results, and is of great significance for analyzing and optimizing the dynamic performance of the power station structure such as earthquake resistance and wind resistance.

[0025] S2. Add the physical simulation model to a pre-set physical simulation engine to generate a basic collision simulation model, obtain the collision initial condition parameters of the two power station components that collide and add them to the basic collision simulation model as boundary conditions, perform local mesh encryption on the power station components through an adaptive mesh subdivision algorithm, define the material constitutive model through a pre-set finite element model in combination with the material data corresponding to the power station components, construct a nonlinear dynamic equation of the collision process based on the nonlinear properties of the power station components, solve the nonlinear dynamic equation in the time domain in combination with an explicit central difference algorithm, generate physical field information, and calculate the strain field and stress field of the power station components through the conservative node average stress method, determine the equivalent stress distribution and the equivalent plastic strain distribution, and generate collision results based on the yield criterion and the fracture failure criterion; The physical simulation engine is used to simulate physical phenomena in the real world, such as mechanics, collision, fluid, etc. The basic collision simulation model is used to simulate collisions between objects, including collision detection and response. The collision initial condition parameters define the initial state of the collision simulation, including the position, speed, direction and collision angle of the object. The adaptive grid subdivision algorithm dynamically adjusts the fineness of the grid according to calculation needs to improve the simulation accuracy. The local grid encryption refers to increasing the grid density in a specific area to more accurately simulate complex physical phenomena. The finite element model approximates the solution of physical problems by dividing complex structures into small discrete units. The material constitutive model describes the response of the material under various loads, including the stress-strain relationship. The nonlinear property refers to the relationship between stress and strain of the material or structure. The system is not linear and may include effects such as plasticity and creep. The nonlinear dynamic equation is used to describe the dynamic behavior of nonlinear mechanical effects in the system. The explicit central difference algorithm is a numerical method for solving time-stepping dynamic problems based on the central difference formula. The time domain solution is a method for analyzing the system behavior that changes with time, including time stepping and dynamic response analysis. The conservative node average stress method is used to calculate and average node stresses in finite element analysis to ensure the conservation of stress. The equivalent stress distribution and equivalent plastic strain distribution are used to describe the internal state of the structure under stress and deformation. The yield criterion is a standard for predicting the beginning of plastic deformation of a material under stress. The fracture failure criterion is used to predict the fracture behavior of a material under the influence of stress and cracks, such as the stress intensity factor criterion.

[0026] In an optional embodiment, The physical simulation model is added to a preset physical simulation engine to generate a basic collision simulation model, the collision initial condition parameters of the two power station components that collide are obtained and added to the basic collision simulation model as boundary conditions, the power station components are locally meshed by an adaptive mesh subdivision algorithm, the material constitutive model is defined by a preset finite element model in combination with the material data corresponding to the power station components, the nonlinear dynamic equation of the collision process is constructed based on the nonlinear properties of the power station components, the nonlinear dynamic equation is solved in the time domain in combination with an explicit central difference algorithm, the physical field information is generated, and the strain field and stress field of the power station components are calculated by the conservative node average stress method, the equivalent stress distribution and the equivalent plastic strain distribution are determined, and the collision results are generated based on the yield criterion and the fracture failure criterion, including: Building a physical simulation engine based on an explicit finite element algorithm, adding the physical simulation model to the physical simulation engine and generating the basic collision simulation model; Perform kinematic analysis on the two power station components that collided, track the trajectories and determine the center of mass position based on real-time monitoring data, determine the angular velocity and contact position at the time of collision in combination with the corresponding geometric shape and posture of each component, combine the center of mass position, angular velocity and contact position to obtain the collision initial condition parameters and add them to the basic collision simulation model as boundary conditions; For the power station components that collide, an initial grid is set by an adaptive grid subdivision algorithm, a stress analysis is performed on the initial grid, the equivalent stress of each unit is determined, and the units whose equivalent stress is greater than a preset grid stress threshold are marked, the marked areas are used as areas that need to be subdivided, and the areas that need to be subdivided are recursively subdivided by an adaptive space partitioning algorithm until the preset grid density requirement is met, and the subdivided grids are smoothed and coupled with the initial grid to obtain a second finite element model; Based on the second finite element model, in combination with the material data corresponding to the power plant component, through the pre-set finite element model, for isotropic materials, a bilinear elastoplastic constitutive model is constructed based on the elastic modulus and yield stress, and for anisotropic materials, a yield function is determined based on the anisotropic parameters corresponding to the current material, and a nonlinear dynamic equation corresponding to the collision process is constructed in combination with the nonlinear properties of the power plant component, wherein the nonlinear properties include large deformation, material nonlinearity and contact nonlinearity; The nonlinear dynamic equations are solved based on an explicit central difference algorithm, the stress, internal force and contact state of each component are updated, physical field information is generated, and the strain field and stress field of the power station components are calculated by the conservation node average stress method. The equivalent stress distribution and the equivalent plastic strain distribution are determined based on the strain field and the stress field. According to the material properties and current state of the power station components, combined with the corresponding yield criterion and fracture failure standard, the collision results are generated and output.

[0027] The mesh density requirement refers to the fineness of the model mesh in finite element analysis or other numerical simulations. The smoothing process refers to improving the surface quality of the mesh model by smoothing, removing sharp corners and irregularities in the model, and making the surface look smoother and more natural. The isotropic material refers to a material that has the same physical properties in all directions, and the anisotropic material has different physical properties in different directions. The elastic modulus is a measure of the material's ability to resist elastic deformation, and the yield stress is the stress value at which the material begins to undergo plastic deformation.

[0028] Select a suitable explicit finite element algorithm, such as the explicit central difference algorithm, build a numerical solver for the physical simulation engine, input the geometric model, material properties, boundary conditions and other information of the power station components into the physical simulation engine, generate the corresponding finite element model, set the collision detection parameters in the physical simulation engine, such as the contact surface search algorithm, contact stiffness, etc., and build a basic collision simulation model; The displacement, velocity and other motion state data of the two power station components that collided are collected in real time by sensors. Based on the kinematic model and numerical integration method, the collected data are tracked to obtain the center of mass positions of the two components at the moment of collision. According to the geometric shape and posture data of the components, the angular velocity and contact position at the time of collision are calculated by analytical or numerical methods. The center of mass position, angular velocity, contact position and other parameters are used as initial collision conditions and added to the basic collision simulation model. According to the geometric characteristics and material properties of the components, the initial finite element mesh is generated by an adaptive meshing algorithm. Stress analysis is performed on the initial mesh, and the equivalent stress of each unit is calculated. The units with equivalent stress greater than the preset threshold are marked as areas that need to be subdivided. Adaptive space partitioning algorithms such as octrees and quadtrees are used to recursively subdivide the areas that need to be subdivided until the preset mesh density requirements are met. The subdivided mesh is smoothed to eliminate mesh distortion and improve mesh quality. The smoothed subdivided mesh is coupled with the initial mesh to form a new finite element model (the second finite element model). According to the material type (isotropic / anisotropic) of the power station components, the corresponding constitutive model is selected. For isotropic materials, the elastic modulus and yield stress of the material are used to construct a bilinear elastoplastic constitutive relationship. For anisotropic materials, anisotropic yield functions are established based on anisotropic parameters (such as the six parameters of the Hill yield criterion). Considering the large deformation, material nonlinearity and contact nonlinearity of the power station components, nonlinear dynamic equations are established on the basis of the second finite element model. The explicit central difference algorithm is used to discretize and solve the nonlinear dynamic equations in the time domain. In each time increment, the node coordinates and node velocity are updated, and the unit stress and node internal force are calculated. The contact point is found according to the contact search algorithm, and the contact force is calculated. The strain field and stress field of the component are obtained based on the extrapolation of the node stress using the conservative node average stress method, and the distribution of equivalent stress and equivalent plastic strain is calculated. Combined with the yield criterion of the material and the fracture failure standard (such as the maximum tensile stress criterion), it is judged whether the component has yielded or fractured, and the collision result data is output.

[0029] For example, suppose that in a nuclear power plant, a steam generator collides with an in-core component. An explicit dynamic finite element solver is constructed through a software package, a finite element model including the steam generator and the in-core components is established, collision detection parameters are set, and a basic collision simulation model is formed. The motion state data of the two components are collected through high-speed video and an inertial measurement unit, and the center of mass position, angular velocity and collision point position when the collision occurs are calculated, and applied to the finite element model as the initial conditions of the collision. The mesh of the collision area is encrypted, and the stress concentration area is subdivided three times by an octree algorithm. The mesh quality is improved by the Laplace smoothing method to form a second finite element model. For the steam generator, considering its elastic-plastic behavior, a bilinear elastic-plastic constitutive model is adopted. The in-core components are alloys, and the corresponding anisotropy is considered. The Hill-48 yield criterion is selected. At the same time, the large deformation, material nonlinearity and contact nonlinearity of the two components during the collision process are considered. The nonlinear dynamic equation is established, and the explicit central difference algorithm is used to solve the dynamic equation with a time step of 1×10-6s to obtain the stress field and strain field of the entire collision process. The Von The Mises criterion and maximum tensile stress criterion are used to evaluate the yield and fracture conditions of components and generate analysis reports.

[0030] In this embodiment, nonlinear factors in the collision process, such as large deformation, material nonlinearity, contact nonlinearity, etc., are fully considered, so that the simulation results are closer to the actual situation. Mesh encryption can better capture the gradient changes of stress and strain, avoid the occurrence of stress singularities, and improve the convergence and stability of the calculation. By adopting different material constitutive models (such as bilinear elastic-plastic model, anisotropic yield criterion, etc.), the mechanical behavior of different materials in the collision process can be accurately described, and the universality of the simulation results can be improved. By explicitly solving nonlinear dynamic equations, various physical field information such as stress field, strain field, displacement field, etc. of power station components in the collision process can be obtained, which can not only intuitively display the impact range and damage degree of the collision, but also provide important data support for subsequent damage assessment, life prediction and other analyses. By introducing the fracture failure criterion, the crack initiation and propagation of the component can be predicted, and the integrity of the component can be quantitatively evaluated. In summary, this embodiment generates rich physical field information and evaluates the impact of the collision on the component, which provides important technical support for accident emergency decision-making and structural optimization design, and realizes high-fidelity dynamic simulation of the collision process of power station components.

[0031] In an optional embodiment, For anisotropic materials, the yield function is determined based on the anisotropic parameters corresponding to the current material as shown in the following formula: ; in, F Indicates that the material y Axis andz Differential resistance to yield between shafts, G Indicates that the material z Axis and x Differential resistance to yield between shafts, H Indicates that the material x Axis and y Differential resistance to yield between shafts, s 11 Indicates that the material x The normal stress component in the axial direction, s 22 Indicates that the material y The normal stress component in the axial direction, s 33 Indicates that the material z The normal stress component in the axial direction, m represents the Horsford index, s y Represents the yield stress.

[0032] In this embodiment, the Horsford criterion takes into account the anisotropy of the material, so that the yield function is closer to the actual mechanical behavior of the material and the accuracy of the yield judgment is improved. By comparing the difference between the three normal stress components, the yield of the material under different stress states can be judged, which improves the comprehensiveness of the yield judgment. By using the yield stress as the judgment standard, it is possible to quantitatively evaluate whether the material is in a yield state. By adjusting the Horsford index, the yield function can better fit the actual yield behavior of the material and improve the applicability of the yield judgment. By embedding the Horsford yield function into the finite element algorithm, the accurate simulation of the nonlinear behavior of the material can be achieved. The explicit yield function is also convenient for structural optimization design. By optimizing the coupling of the objective function and the yield function, structural optimization based on the yield criterion can be achieved. In summary, this embodiment fully considers the influence of the anisotropic characteristics and stress state of the material, improves the accuracy and comprehensiveness of the yield judgment, facilitates numerical solution and optimization design, is widely applicable to various engineering materials, and has important application value in material nonlinear analysis, structural strength assessment, and impact resistance design.

[0033] S3. Based on the collision result, determine whether the power station component has undergone plastic deformation or fracture failure. If plastic deformation has occurred, the finite element model is updated through a node displacement correction algorithm. If fracture failure has occurred, the failed unit is removed from the finite element model through a unit deletion method and a free surface is generated. The simulation is repeated until the collision ends, and the physical field data of the entire collision process is generated. The displacement field data in the physical field data of the entire collision process is mapped to the three-dimensional geometric model. The position and shape of the power station component are updated through a vertex displacement difference algorithm to generate a collision animation frame sequence. Based on the equivalent stress distribution and the equivalent plastic strain distribution, a stress distribution cloud map and a plastic strain cloud map are generated in combination with a volume rendering algorithm. For power station components that have undergone fracture failure, fracture texture and fracture morphology are generated according to the position and normal vector of the fracture unit through a grid breaking algorithm. Fracture special effects are constructed in combination with splashing fragments. A collision visualization scene is generated by integrating the collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effects.

[0034] The node displacement correction algorithm is used to adjust the node displacement in the finite element model, which is usually a correction step applied in the calculation process to improve the accuracy of the calculation results or solve numerical instability. The unit deletion method is a technology for dealing with fracture or failure in the finite element model. The fracture process of the material or structure is simulated by deleting the failed unit from the model. The free surface refers to the surface that is not constrained or affected by external forces, which is usually used to describe the surface of fluids, particles or debris. The full collision process physical field data includes all physical data collected in the collision event, such as force, stress, deformation and temperature change. The vertex displacement difference algorithm is used to calculate and interpolate the position of the vertex in the three-dimensional mesh model. The collision animation frame sequence refers to the animation sequence recorded and displayed frame by frame in the collision simulation, which is used to show the dynamic process of the collision event. The volume rendering algorithm is used to convert voxel data (such as medical scan data) into a visual image. The mesh breaking algorithm is used to simulate and generate the breaking effect in the mesh model, and the mesh is usually decomposed into fragments or broken parts through an algorithm. The fracture texture refers to the pattern or texture used to describe and render the fracture surface. The fracture morphology describes the surface characteristics and shape of the material or structure after fracture, including the shape, direction and roughness of the crack. The splashing fragments refer to the fragments flying from the surface of the object during the collision or explosion.

[0035] In an optional embodiment, Based on the collision result, it is judged whether the power station component has plastic deformation or fracture failure. If plastic deformation occurs, the finite element model is updated by a node displacement correction algorithm. If fracture failure occurs, the failed unit is removed from the finite element model by a unit deletion method and a free surface is generated. The simulation is repeated until the collision ends, and the physical field data of the entire collision process is generated. The displacement field data in the physical field data of the entire collision process is mapped to the three-dimensional geometric model. The position and shape of the power station component are updated by a vertex displacement difference algorithm, and a collision animation frame sequence is generated. Based on the equivalent stress distribution and the equivalent plastic strain distribution, a stress distribution cloud map and a plastic strain cloud map are generated in combination with a volume rendering algorithm. For the power station component that has fracture failure, according to the position and normal vector of the fracture unit, a mesh crushing algorithm is used to generate a fracture texture and a fracture morphology, and a fracture special effect is constructed in combination with splashing fragments. The collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effect are integrated to generate a collision visualization scene, including: Based on the collision result, determine whether the power station component has plastic deformation or fracture failure. If plastic deformation occurs, extract the displacement field data of the current component based on the node displacement correction algorithm. For each node, update the node coordinates according to the corresponding displacement vector, generate a finite element mesh based on the updated node coordinates, and update the material properties and plastic strain of the current component. Repeat the update until a finite element model sequence corresponding to the current component is obtained. If the collision result indicates that the current power plant component has fractured and failed, then based on the fracture failure criterion, the failed unit is identified and a failed unit set is generated, for each failed unit, the corresponding node set is extracted and the node set is removed, and nodes are inserted at the boundary of the fracture failure area and connected with the finite element mesh before the fracture to obtain the free surface, and the finite element mesh is regenerated based on the free surface and the material properties and plastic strain of the current component are updated, and the update is repeated until the finite element model sequence corresponding to the current component is obtained; Extracting the physical field data of the entire collision process based on the finite element model sequence, mapping the displacement field in the physical field data of the entire collision process to the three-dimensional geometric model to generate a vertex displacement field, updating the vertex coordinates of the three-dimensional geometric model through a vertex displacement difference algorithm and calculating the corresponding normal vector and texture coordinates, and updating the position and shape of the power station component in combination with the physical field data of the entire collision process to obtain the collision animation frame sequence; For the power station structure, based on the collision animation frame sequence and the equivalent stress distribution and equivalent plastic strain distribution, a color mapping table for each voxel is determined by a volume rendering algorithm, and a search is performed in the color mapping table based on the stress value and the strain value. A ray is emitted to each pixel point by a ray casting algorithm, and the intersection with different voxels is determined, and the color corresponding to each pixel point is determined, and the operation is repeated to obtain the stress distribution cloud map and the plastic strain cloud map; For power station components that have failed due to fracture, a random crack distribution is generated on the boundary surface corresponding to the fracture unit according to the position and normal vector of the fracture unit through a mesh breaking algorithm according to the fracture direction and the material data, a polygonal mesh is generated with each crack line in the random crack distribution as a starting point, the polygonal mesh is triangulated and a normal map and a height map are attached to generate an initial fracture morphology and random disturbances are added to obtain fracture texture and fracture morphology, fragment particles are randomly generated in the fracture unit, the kinematic parameters of each fragment particle are initialized, and the scattering trajectory corresponding to each fragment particle is determined by explicit integration, a fracture special effect is constructed based on the scattering trajectory and the fragment particles, and a collision visualization scene is generated by combining the collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effect.

[0036] The displacement vector represents the displacement of an object in a certain direction, and usually includes three components: displacement in the X, Y and Z directions. The finite element mesh is a mesh that divides the calculation area into small discrete units (such as triangles or tetrahedrons), which is used to approximately solve problems in finite element analysis. The failure unit set refers to a set of mesh units that are identified as failed or broken in finite element analysis. The fracture failure area refers to an area in a material or structure where fracture or failure occurs, usually manifested as a crack or fracture. The vertex displacement field describes the spatial distribution of the displacement of each vertex in the mesh model, including the actual displacement of each vertex. The color mapping table is used to map numerical data (such as temperature, stress) to a color range for visualization in a graphical display. The ray casting algorithm is used to sample the scene along the direction of the light from the viewpoint to generate an image or perform volume rendering. The random crack distribution refers to generating cracks in a random pattern on the surface or structure of the material to simulate the distribution of natural cracks. The normal map is a texture mapping technology used to simulate details and bump effects on the surface, and enhance surface details by storing the normal direction. The height map is an image used to represent the change in surface height, and the terrain or surface bumps are defined by grayscale values. The debris particles refer to small broken pieces or particles generated during the rupture or collision process, which are used to simulate the scattering effect of the material. The scattering trajectory describes the path and movement trajectory of the debris or splashes during the collision or explosion.

[0037] Extract the equivalent stress and equivalent plastic strain data from the collision results, compare the equivalent stress with the yield stress of the material, determine whether plastic deformation occurs, if the equivalent stress is greater than the yield stress, plastic deformation occurs, extract the displacement field data of the component, update the node coordinates for each node according to the displacement vector, generate a new finite element mesh based on the updated node coordinates, update the material properties and plastic strain of the component, repeat the update until the finite element model sequence of the current component is obtained, if the equivalent stress is less than the yield stress, no plastic deformation occurs; Substitute the equivalent stress and equivalent plastic strain into the fracture failure criterion to determine whether fracture failure occurs. If the fracture failure criterion is met, fracture failure occurs. Failure units are identified based on the fracture failure criterion, and a set of failure units is generated. For each failure unit, the corresponding node set is extracted and removed. New nodes are inserted at the boundary of the fracture failure area and connected with the mesh before fracture to form a free surface. The finite element mesh is regenerated based on the free surface, and the material properties and plastic strain of the component are updated. The update is repeated until the finite element model sequence of the component is obtained. If the fracture failure criterion is not met, fracture failure does not occur. Extract the physical field data of the whole collision process based on the finite element model sequence, map the displacement field data to the three-dimensional geometric model, generate the vertex displacement field, update the vertex coordinates of the three-dimensional model through the vertex displacement difference algorithm, calculate the updated vertex normal vector and texture coordinates, update the position and shape of the component in combination with the physical field data, and repeat the update until the collision animation frame sequence is obtained. Based on the collision animation frame sequence and the equivalent stress and plastic strain distribution, determine the color mapping table of each voxel through the volume rendering algorithm, for each pixel point, find the color in the color mapping table based on the stress value and strain value, emit light to each pixel point through the ray casting algorithm, determine the intersection with different voxels, determine the color of each pixel point based on the intersection information, and repeat the determination until the stress distribution cloud map and the plastic strain cloud map are obtained; Generate random crack distribution on the boundary surface of the fracture unit through the mesh breaking algorithm, generate polygonal mesh with each crack line as the starting point, triangulate the polygonal mesh, and add normal map and height map to generate the initial fracture morphology, add random perturbation to obtain the final fracture texture and fracture morphology, randomly generate fragment particles in the fracture unit, initialize the kinematic parameters of each fragment particle, determine the scattering trajectory of each fragment particle through explicit integration, build fracture special effects based on the scattering trajectory and fragment particles, and generate the final collision visualization scene by integrating collision animation frame sequence, stress distribution cloud map, plastic strain cloud map and fracture special effects; For example, a pump body of a power station was hit by an external object during operation, and a visual analysis of the collision process was required. Based on the collision simulation results, it was determined whether the pump body had plastic deformation or fracture failure. By comparing the equivalent stress and yield stress, it was found that the maximum equivalent stress was 400MPa, which was greater than the yield strength of the pump body material of 350MPa. Therefore, the pump body had plastic deformation. The equivalent stress and plastic strain were substituted into the fracture failure criterion for judgment. It was found that some units met the failure conditions, so the pump body also had local fracture failure. For the area where plastic deformation occurs, the displacement field is extracted and the node coordinates are updated, a new finite element mesh is generated, the material properties and plastic strain are updated, and the process is repeated many times to obtain a sequence of finite element models of the deformed pump body. For the area where fracture failure occurs, the nodes of the failed unit are extracted and removed from the mesh, new nodes are inserted at the fracture boundary and the mesh is regenerated, the material properties and plastic strain are updated, and the process is repeated many times to obtain a sequence of finite element models of the pump body after fracture. Based on the finite element model sequence, the physical field data of the entire collision process is extracted, the displacement field is mapped to the three-dimensional model of the pump body to generate a vertex displacement field, the vertex coordinates, normal vectors and texture coordinates are updated, and the position and shape of the pump body are updated in combination with the physical field data to generate a collision animation frame sequence. Based on the collision animation frame sequence and the equivalent stress / strain distribution, the stress cloud map and plastic strain cloud map of the pump body are generated through volume rendering and ray casting algorithms, which intuitively display the stress and deformation of the pump body. For the fracture failure area, a random crack distribution is generated through a mesh breaking algorithm, a polygonal mesh is constructed and triangulated, and the fracture morphology and fracture texture are generated. Fragment particles are generated in the fracture unit, the scattering trajectory is calculated, and the fracture special effects are constructed. Finally, the collision animation, stress cloud map, strain cloud map, and fracture special effects are integrated to generate a visualization scene of the entire pump body collision process, which intuitively displays the complete process of deformation, fracture, and crushing of the pump body under the action of collision.

[0038] In this embodiment, by extracting the equivalent stress and equivalent plastic strain data in the collision results and comparing them with the material yield strength and fracture failure criterion, it is accurately judged whether the component has plastic deformation and fracture failure, which provides a basis for subsequent visualization processing. The node coordinates are updated by the node displacement correction algorithm, and the finite element mesh is generated based on the new coordinates. The material properties and plastic strain are updated at the same time. The deformed component model sequence is obtained through multiple iterations, which truly reproduces the progressive deformation process of the component under the collision load. By identifying the failed units and removing the corresponding nodes, new nodes are inserted at the fracture boundary and the mesh is regenerated. The material properties and plastic strain are updated at the same time. The component model sequence after fracture is obtained through multiple iterations, which truly reproduces the component fracture. The complex process of failure is analyzed. By extracting the physical field data corresponding to the finite element model sequence, rich data information such as the displacement field, stress field and strain field of the component in the whole collision process is obtained, which provides a data basis for visual rendering. By mapping the displacement field to the three-dimensional model of the component and updating the vertex coordinates, normal vectors and texture coordinates, the position and shape of the component are updated in real time in combination with the physical field data, and a continuous, smooth and fluent collision animation frame sequence is generated, which intuitively displays the movement changes of the component during the collision process. In summary, this embodiment not only intuitively displays the deformation and destruction process of the component, but also reveals the stress and strain distribution law inside the component, which provides intuitive and reliable technical support for collision consequence analysis, accident investigation, safety assessment and other work.

[0039] In an optional embodiment, The ray casting algorithm is used to emit light to each pixel and determine the intersection points with different voxels. The color corresponding to each pixel is determined as shown in the following formula: ; in, CP ( u , v ) represents pixel points ( u , v ), CV ( x , y , z ) represents a voxel ( x , y , z ), S ( x , y , z ) represents a voxel ( x , y , z ) scattering coefficient, t ( x , y , z ) represents a voxel ( x ,y , z )’s absorption coefficient, l is the path of light propagation, indicating the propagation distance of light in the voxel, dλ is the light element, used to calculate the path of the light l The points on .

[0040] In this embodiment, a ray casting algorithm is used to emit light to each pixel, and the intersection of the light and the voxel is calculated. The color of the pixel is calculated according to the voxel properties (color, scattering coefficient, absorption coefficient) at the intersection, so that realistic rendering of three-dimensional volume data is achieved, and realistic visualization results with translucent, lighting, shadow and other effects are generated. The absorption, scattering and transmission effects of voxels on light are comprehensively considered, so that the rendering results are more realistic and can accurately represent the optical properties of different materials. By integrating the absorption coefficient on the light propagation path, the energy attenuation process of light propagating in the voxel is simulated, so that the voxels farther away from the light source have darker colors, and the voxels farther away from the light source have darker colors. Voxels closer to the light source are brighter in color, producing a realistic light attenuation effect. By introducing light elements in the integral operation and discretizing the light propagation path, the efficiency of the rendering calculation is improved. By setting appropriate sampling steps and termination conditions, unnecessary calculations are reduced while ensuring the rendering quality, thereby improving the execution efficiency of the algorithm. In summary, this embodiment achieves high-quality, highly realistic volume data rendering, generates visualization results with effects such as translucency, illumination, and shadow, restores the physical mechanism of the interaction between light and objects in the real world, makes the rendering results more realistic and natural, and provides strong technical support for the visualization analysis of complex volume data.

[0041] Figure 2 FIG. 1 is a schematic diagram of a twin power station collision visualization simulation system based on physical simulation according to an embodiment of the present invention. Figure 2 As shown, the system comprises: The second unit is used to add the physical simulation model to a preset physical simulation engine to generate a basic collision simulation model, obtain the collision initial condition parameters of the two power station components that collide and add them to the basic collision simulation model as boundary conditions, perform local mesh encryption on the power station components through an adaptive mesh subdivision algorithm, define the material constitutive model through a preset finite element model in combination with the material data corresponding to the power station components, construct a nonlinear dynamic equation of the collision process based on the nonlinear properties of the power station components, solve the nonlinear dynamic equation in the time domain in combination with an explicit central difference algorithm, generate physical field information, calculate the strain field and stress field of the power station components through a conservative node average stress method, determine the equivalent stress distribution and the equivalent plastic strain distribution, and generate a collision result based on the yield criterion and the fracture failure criterion; The third unit is used to judge whether the power station component has plastic deformation or fracture failure based on the collision result. If plastic deformation occurs, the finite element model is updated through a node displacement correction algorithm. If fracture failure occurs, the failed unit is removed from the finite element model through a unit deletion method and a free surface is generated. The simulation is repeated until the collision ends, and the physical field data of the entire collision process is generated. The displacement field data in the physical field data of the entire collision process is mapped to the three-dimensional geometric model, and the position and shape of the power station component are updated through a vertex displacement difference algorithm to generate a collision animation frame sequence. Based on the equivalent stress distribution and the equivalent plastic strain distribution, a stress distribution cloud map and a plastic strain cloud map are generated in combination with a volume rendering algorithm. For the power station component that has fracture failure, according to the position and normal vector of the fracture unit, a mesh crushing algorithm is used to generate a fracture texture and a fracture morphology, and a fracture special effect is constructed in combination with splashing fragments. The collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effect are integrated to generate a collision visualization scene.

[0042] According to a third aspect of the embodiments of the present invention, An electronic device is provided, comprising: processor; a memory for storing processor-executable instructions; The processor is configured to call the instructions stored in the memory to execute the aforementioned method.

[0043] According to a fourth aspect of the embodiments of the present invention, A computer-readable storage medium is provided, on which computer program instructions are stored. When the computer program instructions are executed by a processor, the aforementioned method is implemented.

[0044] The present invention may be a method, an apparatus, a system and / or a computer program product. The computer program product may include a computer-readable storage medium carrying computer-readable program instructions for executing various aspects of the present invention.

[0045] 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 aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A twin power station collision visualization simulation method based on physical simulation, characterized in that: include: Acquire a three-dimensional geometric model of the twin power station, perform spatial division of the twin power station by using an octree space division algorithm, generate a voxelized model, and extract geometric parameters of power station components in the twin power station based on the topological relationship of the voxel units, acquire material data of the power station components based on the physical design specification and bill of materials of the twin power station, associate the geometric parameters with the material data, and construct a material association model, acquire a mechanical model of the internal frame structure of the physical power station according to the structural design drawing of the physical power station, calculate the force distribution in combination with mechanical principles, obtain a mechanical simulation model, and fuse it with the material association model to obtain a physical simulation model; Add the physical simulation model to a preset physical simulation engine to generate a basic collision simulation model, obtain the collision initial condition parameters of the two power station components that collide and add them to the basic collision simulation model as boundary conditions, perform local mesh encryption on the power station components through an adaptive mesh subdivision algorithm, define the material constitutive model through a preset finite element model in combination with the material data corresponding to the power station components, construct a nonlinear dynamic equation of the collision process based on the nonlinear properties of the power station components, solve the nonlinear dynamic equation in the time domain in combination with an explicit central difference algorithm, generate physical field information, calculate the strain field and stress field of the power station components through a conservative node average stress method, determine the equivalent stress distribution and the equivalent plastic strain distribution, and generate a collision result based on a yield criterion and a fracture failure criterion; Based on the collision result, it is judged whether the power station component has plastic deformation or fracture failure. If plastic deformation occurs, the finite element model is updated through a node displacement correction algorithm. If fracture failure occurs, the failed unit is removed from the finite element model through a unit deletion method and a free surface is generated. The simulation is repeated until the collision ends, and the physical field data of the entire collision process is generated. The displacement field data in the physical field data of the entire collision process is mapped to the three-dimensional geometric model. The position and shape of the power station component are updated through a vertex displacement difference algorithm to generate a collision animation frame sequence. Based on the equivalent stress distribution and the equivalent plastic strain distribution, a stress distribution cloud map and a plastic strain cloud map are generated in combination with a volume rendering algorithm. For the power station component that has fracture failure, according to the position and normal vector of the fracture unit, a mesh crushing algorithm is used to generate a fracture texture and a fracture morphology. A fracture special effect is constructed in combination with splashing fragments. A collision visualization scene is generated by integrating the collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effect.

2. The method according to claim 1, characterized in that The three-dimensional geometric model of the twin power station is obtained, the twin power station is spatially divided by an octree space partitioning algorithm, a voxelized model is generated, and the geometric parameters of the power station components in the twin power station are extracted based on the topological relationship of the voxel units. The material data of the power station components are obtained based on the physical design specification and the bill of materials of the twin power station, the geometric parameters are associated with the material data, and a material association model is constructed. The mechanical model of the internal frame structure of the physical power station is obtained according to the structural design drawing of the physical power station, the force distribution is calculated in combination with the mechanical principle, and a mechanical simulation model is obtained and merged with the material association model to obtain a physical simulation model including: The point cloud data corresponding to the twin power station is acquired by three-dimensional laser scanning technology, and the three-dimensional geometric model corresponding to the twin power station is obtained by point cloud data preprocessing and registration and splicing operations. Based on the octree space partitioning algorithm, the root node of the octree is determined according to the bounding box of the twin power station, and the complexity metric is determined at each node according to the number of point clouds. The complexity threshold is combined to make a judgment. If the complexity of the current node is greater than the complexity threshold, the current node is recursively divided into eight equal parts, and the judgment is repeated until the complexity of all nodes is less than the complexity threshold to obtain a leaf node. For each leaf node, a corresponding voxel unit is generated and combined to obtain the voxelized model; Based on the topological relationship corresponding to the voxel units, different voxel sets corresponding to the power station components in the twin power station are identified, and the outer envelope box of the component is calculated based on the voxel set belonging to the same component to obtain the size information of the current component, and the shape feature of the current component is extracted by fitting the surface mesh of the voxel set, and the geometric center of the outer envelope box is used as the spatial position coordinate of the current component. The size information, shape feature and spatial position coordinate are combined to obtain the geometric parameters of the current component, and each component in the twin power station is traversed to obtain the corresponding geometric parameters; Based on the physical design specifications and bill of materials of the twin power stations, the material composition of the power station components in each power station is obtained. For metal materials, a material property database is constructed in combination with physical and mechanical property parameters, and the property parameters of the material composition of each component are retrieved from the material property database. For non-metallic materials and new materials, mechanical performance parameters are determined in combination with simulation experimental tests, and the material data is obtained by combining the property parameters and the mechanical performance parameters. The bill of materials is traversed and a first relationship table and a second relationship table corresponding to the power station components and the geometric parameters and the material data are respectively generated, the first relationship table and the second relationship table are connected, and material property parameters or mechanical performance parameters are retrieved, and the material association model is obtained by combining them; According to the structural design drawing of the physical power station, a mechanical model of the internal frame structure of the physical power station is obtained, and the key load mode is determined by combining mechanical principles and load condition analysis technology, wherein the key load mode is divided into static load and dynamic load. For the static load, the static force distribution is calculated by a finite element static analysis method, and for the dynamic load, the vibration mode and frequency response characteristics of the internal frame structure are obtained by a finite element dynamic analysis method. Based on the vibration mode and frequency response characteristics, combined with a transient mechanical analysis method, the dynamic force distribution is obtained by solving the dynamic equation by a direct integration method. The static force distribution and the dynamic force distribution are combined to obtain a mechanical simulation model, and the static force distribution and the dynamic force distribution in the mechanical simulation model are mapped to the material association model to obtain the physical simulation model.

3. The method according to claim 2, characterized in that Based on the vibration mode and frequency response characteristics, combined with the transient mechanical analysis method, the dynamic force distribution is obtained by solving the dynamic equation through direct integration method as shown in the following formula: ； in, F ( t n+1 ) indicates that at the time step n +1 external force load, [M] represents the mass matrix, { u n+1 } represents the time step n +1 displacement vector, [C] represents the damping matrix, Represents the time step n The velocity vector, △ t represents the time increment, Represents the time step n The acceleration vector, Represents the time step n +1 acceleration vector, [K] represents the stiffness matrix.

4. The method according to claim 1, characterized in that: The physical simulation model is added to a preset physical simulation engine to generate a basic collision simulation model, the collision initial condition parameters of the two power station components that collide are obtained and added to the basic collision simulation model as boundary conditions, the power station components are locally meshed by an adaptive mesh subdivision algorithm, the material constitutive model is defined by a preset finite element model in combination with the material data corresponding to the power station components, the nonlinear dynamic equation of the collision process is constructed based on the nonlinear properties of the power station components, the nonlinear dynamic equation is solved in the time domain in combination with an explicit central difference algorithm, the physical field information is generated, and the strain field and stress field of the power station components are calculated by the conservative node average stress method, the equivalent stress distribution and the equivalent plastic strain distribution are determined, and the collision results are generated based on the yield criterion and the fracture failure criterion, including: Building a physical simulation engine based on an explicit finite element algorithm, adding the physical simulation model to the physical simulation engine and generating the basic collision simulation model; Perform kinematic analysis on the two power station components that collided, track the trajectories and determine the center of mass position based on real-time monitoring data, determine the angular velocity and contact position at the time of collision in combination with the corresponding geometric shape and posture of each component, combine the center of mass position, angular velocity and contact position to obtain the collision initial condition parameters and add them to the basic collision simulation model as boundary conditions; For the power station components that collide, an initial grid is set by an adaptive grid subdivision algorithm, a stress analysis is performed on the initial grid, the equivalent stress of each unit is determined, and the units whose equivalent stress is greater than a preset grid stress threshold are marked, the marked areas are used as areas that need to be subdivided, and the areas that need to be subdivided are recursively subdivided by an adaptive space partitioning algorithm until the preset grid density requirement is met, and the subdivided grids are smoothed and coupled with the initial grid to obtain a second finite element model; Based on the second finite element model, in combination with the material data corresponding to the power plant component, through the pre-set finite element model, for isotropic materials, a bilinear elastoplastic constitutive model is constructed based on the elastic modulus and yield stress, and for anisotropic materials, a yield function is determined based on the anisotropic parameters corresponding to the current material, and a nonlinear dynamic equation corresponding to the collision process is constructed in combination with the nonlinear properties of the power plant component, wherein the nonlinear properties include large deformation, material nonlinearity and contact nonlinearity; The nonlinear dynamic equations are solved based on an explicit central difference algorithm, the stress, internal force and contact state of each component are updated, physical field information is generated, and the strain field and stress field of the power station components are calculated by the conservation node average stress method. The equivalent stress distribution and the equivalent plastic strain distribution are determined based on the strain field and the stress field. According to the material properties and current state of the power station components, combined with the corresponding yield criterion and fracture failure standard, the collision results are generated and output.

5. The method according to claim 4, characterized in that For anisotropic materials, the yield function is determined based on the anisotropic parameters corresponding to the current material as shown in the following formula: ； in, F Indicates that the material y Axis and z Differential resistance to yield between shafts, G Indicates that the material z Axis and x Differential resistance to yield between shafts, H Indicates that the material x Axis and y Differential resistance to yield between shafts, σ 11 Indicates that the material x The normal stress component in the axial direction, σ 22 Indicates that the material y The normal stress component in the axial direction, σ 33 Indicates that the material z The normal stress component in the axial direction, m represents the Horsford index, σ y Represents the yield stress.

6. The method according to claim 1, characterized in that Based on the collision result, it is judged whether the power station component has plastic deformation or fracture failure. If plastic deformation occurs, the finite element model is updated by a node displacement correction algorithm. If fracture failure occurs, the failed unit is removed from the finite element model by a unit deletion method and a free surface is generated. The simulation is repeated until the collision ends, and the physical field data of the entire collision process is generated. The displacement field data in the physical field data of the entire collision process is mapped to the three-dimensional geometric model. The position and shape of the power station component are updated by a vertex displacement difference algorithm, and a collision animation frame sequence is generated. Based on the equivalent stress distribution and the equivalent plastic strain distribution, a stress distribution cloud map and a plastic strain cloud map are generated in combination with a volume rendering algorithm. For the power station component that has fracture failure, according to the position and normal vector of the fracture unit, a mesh crushing algorithm is used to generate a fracture texture and a fracture morphology, and a fracture special effect is constructed in combination with splashing fragments. The collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effect are integrated to generate a collision visualization scene, including: Based on the collision result, determine whether the power station component has plastic deformation or fracture failure. If plastic deformation occurs, extract the displacement field data of the current component based on the node displacement correction algorithm. For each node, update the node coordinates according to the corresponding displacement vector, generate a finite element mesh based on the updated node coordinates, and update the material properties and plastic strain of the current component. Repeat the update until a finite element model sequence corresponding to the current component is obtained. If the collision result indicates that the current power plant component has fractured and failed, then based on the fracture failure criterion, the failed unit is identified and a failed unit set is generated, for each failed unit, the corresponding node set is extracted and the node set is removed, and nodes are inserted at the boundary of the fracture failure area and connected with the finite element mesh before the fracture to obtain the free surface, and the finite element mesh is regenerated based on the free surface and the material properties and plastic strain of the current component are updated, and the update is repeated until the finite element model sequence corresponding to the current component is obtained; Extracting the physical field data of the entire collision process based on the finite element model sequence, mapping the displacement field in the physical field data of the entire collision process to the three-dimensional geometric model to generate a vertex displacement field, updating the vertex coordinates of the three-dimensional geometric model through a vertex displacement difference algorithm and calculating the corresponding normal vector and texture coordinates, and updating the position and shape of the power station component in combination with the physical field data of the entire collision process to obtain the collision animation frame sequence; For the power station structure, based on the collision animation frame sequence and the equivalent stress distribution and equivalent plastic strain distribution, a color mapping table for each voxel is determined by a volume rendering algorithm, and a search is performed in the color mapping table based on the stress value and the strain value. A ray is emitted to each pixel point by a ray casting algorithm, and the intersection with different voxels is determined, and the color corresponding to each pixel point is determined, and the operation is repeated to obtain the stress distribution cloud map and the plastic strain cloud map; For power station components that have failed due to fracture, a random crack distribution is generated on the boundary surface corresponding to the fracture unit according to the position and normal vector of the fracture unit through a mesh breaking algorithm according to the fracture direction and the material data, a polygonal mesh is generated with each crack line in the random crack distribution as a starting point, the polygonal mesh is triangulated and a normal map and a height map are attached to generate an initial fracture morphology and random disturbances are added to obtain fracture texture and fracture morphology, fragment particles are randomly generated in the fracture unit, the kinematic parameters of each fragment particle are initialized, and the scattering trajectory corresponding to each fragment particle is determined by explicit integration, a fracture special effect is constructed based on the scattering trajectory and the fragment particles, and a collision visualization scene is generated by combining the collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effect.

7. The method according to claim 6, characterized in that The ray casting algorithm is used to emit light to each pixel and determine the intersection points with different voxels. The color corresponding to each pixel is determined as shown in the following formula: ； in, CP ( u , v ) represents pixel points ( u , v ), CV ( x , y , z ) represents a voxel ( x , y , z ), S ( x , y , z ) represents a voxel ( x , y , z ) scattering coefficient, τ ( x , y , z ) represents a voxel ( x , y , z )’s absorption coefficient, λ is the path of light propagation, indicating the propagation distance of light in the voxel, dλ is the light element, used to calculate the path of the light λ The points on .

8. A twin power plant collision visualization simulation system based on physical simulation, used to implement the method described in any one of claims 1 to 7, characterized in that: include: The first unit is used to obtain a three-dimensional geometric model of a twin power station, perform spatial division of the twin power station by an octree space division algorithm, generate a voxelized model, and extract geometric parameters of power station components in the twin power station based on the topological relationship of the voxel units, obtain material data of the power station components based on the physical design specification and the bill of materials of the twin power station, associate the geometric parameters with the material data, construct a material association model, obtain a mechanical model of the internal frame structure of the physical power station according to the structural design drawing of the physical power station, calculate the force distribution in combination with the mechanical principle, obtain a mechanical simulation model, and fuse it with the material association model to obtain a physical simulation model; The second unit is used to add the physical simulation model to a preset physical simulation engine to generate a basic collision simulation model, obtain the collision initial condition parameters of the two power station components that collide and add them to the basic collision simulation model as boundary conditions, perform local mesh encryption on the power station components through an adaptive mesh subdivision algorithm, define the material constitutive model through a preset finite element model in combination with the material data corresponding to the power station components, construct a nonlinear dynamic equation of the collision process based on the nonlinear properties of the power station components, solve the nonlinear dynamic equation in the time domain in combination with an explicit central difference algorithm, generate physical field information, calculate the strain field and stress field of the power station components through a conservative node average stress method, determine the equivalent stress distribution and the equivalent plastic strain distribution, and generate a collision result based on the yield criterion and the fracture failure criterion; The third unit is used to judge whether the power station component has plastic deformation or fracture failure based on the collision result. If plastic deformation occurs, the finite element model is updated through a node displacement correction algorithm. If fracture failure occurs, the failed unit is removed from the finite element model through a unit deletion method and a free surface is generated. The simulation is repeated until the collision ends, and the physical field data of the entire collision process is generated. The displacement field data in the physical field data of the entire collision process is mapped to the three-dimensional geometric model, and the position and shape of the power station component are updated through a vertex displacement difference algorithm to generate a collision animation frame sequence. Based on the equivalent stress distribution and the equivalent plastic strain distribution, a stress distribution cloud map and a plastic strain cloud map are generated in combination with a volume rendering algorithm. For the power station component that has fracture failure, according to the position and normal vector of the fracture unit, a mesh crushing algorithm is used to generate a fracture texture and a fracture morphology, and a fracture special effect is constructed in combination with splashing fragments. The collision animation frame sequence, the stress distribution cloud map, the plastic strain cloud map and the fracture special effect are integrated to generate a collision visualization scene.

9. An electronic device, characterized in that: include: processor; a memory for storing processor-executable instructions; The processor is configured to call the instructions stored in the memory to execute the method described in any one of claims 1 to 7.

10. A computer-readable storage medium having computer program instructions stored thereon, characterized in that: When the computer program instructions are executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images