Acceleration method, system, apparatus and medium for fuel cladding wear prediction

Abstract

The invention discloses an acceleration method for fuel cladding wear prediction, and relates to the technical field of fretting wear prediction of nuclear reactors, and the method comprises the following steps: establishing a part wear area contact local refinement three-dimensional finite element model; reducing the dimension of the three-dimensional model into a two-dimensional model, and converting the contact resultant force in the three-dimensional finite element model into equivalent force applied to the two-dimensional wear finite element model based on the principle that the maximum contact pressure of a contact area is equal; screening load fragments contributing to wear in the wear load, and splicing the load fragments to obtain an equivalent external load for wear calculation; and multiplying the node wear depth in the representative wear time period by the number of the representative wear time periods to obtain an accelerated wear result, and updating the shape of the part surface model and the finite element mesh according to the wear result. According to the method, fretting wear analysis can be carried out on the contact surface in any shape, and the surface wear morphology of the part at any moment can be obtained.

Description

Technical Field

[0001] This invention relates to the field of fretting wear prediction technology for nuclear reactors, and specifically to an accelerated method, system, device, and medium for predicting fuel cladding wear. Background Technology

[0002] Grid-to-rod fretting wear (GTRF) caused by flow-induced vibration of fuel rods is the most important fuel failure mechanism. During long-term service in nuclear reactors, the wear of fuel assemblies gradually intensifies, and when the wear reaches a certain depth, it will affect the safety of the structure.

[0003] Traditional methods obtain the interaction between the component and the support through flow-induced vibration analysis, and use the mathematical expression of wear volume versus wear power obtained from fretting wear tests to calculate the wear volume of the fuel rod during its service life. Then, based on an approximate formula for a specific support type, the wear volume is converted into an engineering design parameter—wear depth. This method is computationally efficient, but the calculations are based on ideal wear morphologies under a limited number of support types. The actual wear track morphology and depth deviate from the theoretical formula (e.g., ...). Figure 1 As shown in the figure, its applicability is not high and its accuracy needs to be improved.

[0004] Due to the special nature and importance of fuel rods, recent studies have begun to focus on refined numerical simulations based on the finite element method (element scale 10). -6 (m) Simulates the wear track morphology in the wear process from a microscopic perspective. By monitoring the wear power of each finite element between the fuel rod and the support during reciprocating relative motion, the finite element size of the contact surface is updated in a timely manner, thus achieving a more accurate simulation of the wear track morphology.

[0005] Fuel rod fretting wear has the following characteristics: Firstly, fretting wear is a long-cycle process with high vibration frequency, undergoing multiple changes in motion direction within a short period of time, and its service life is nearly 1×10⁻⁶ days. 9 In order to accurately reflect these fretting behaviors, the incremental time step size during the calculation and analysis cannot exceed the average turning time of the fretting, so as to avoid the large time increment step causing the fretting behavior to be swallowed up, thus leading to an underestimation of the relative sliding distance. Therefore, reflecting the dynamic evolution process of fuel assembly fretting wear through refined finite element numerical simulation requires massive computational resources. On the other hand, fretting wear has a small change in the surface morphology of the wear track in a short period of time, and the changes in vibration characteristics and wear behavior are extremely slow. The massive computational resources contain very low information density, resulting in low efficiency of finite element numerical simulation of fretting wear. Summary of the Invention

[0006] The purpose of this invention is to provide an accelerated method, system, device, and medium for predicting fuel cladding wear, thereby solving the technical problems existing in the prior art.

[0007] The present invention is achieved through the following technical solutions: In a first aspect, an accelerated method for predicting fuel cladding wear provided by embodiments of the present invention includes: establishing a refined three-dimensional finite element model of the contact area of ​​the wear region of a component; preserving geometric features, reducing the three-dimensional finite element model to a two-dimensional wear finite element model; based on the principle that the maximum contact pressure in the contact area is equal, converting the contact resultant force in the three-dimensional finite element model into an equivalent force applied to the two-dimensional wear finite element model; performing equivalent processing on the load, screening load segments that contribute to wear in the wear load, splicing the load segments to obtain the equivalent external load for wear calculation; controlling the movement of each part by applying tangential displacement boundary conditions to simulate the relative sliding displacement caused by the external load, and controlling the squeezing interaction of each part by applying normal force load conditions to simulate the contact surface squeezing caused by the external load to obtain the contact pressure; The nodal wear depth is calculated based on the wear depth formula, contact pressure, and relative sliding displacement for a representative wear period. Assuming that the wear power is the same as the wear power in the representative wear period within a wear morphology update cycle, the nodal wear depth in the representative wear period is multiplied by the number of representative wear periods to obtain the accelerated wear result. The surface model shape and finite element mesh of the component are updated based on the wear result.

[0008] Furthermore, the specific method for establishing a refined three-dimensional finite element model of the contact area in the wear region of the component includes: The geometric model is meshed using a mesh generation tool, and the mesh is refined locally in the contact area. Contact pairs and adaptive meshes are provided on the surface of the pre-wear area, which includes areas where contact and fretting wear may occur.

[0009] Furthermore, the specific method for splicing the load segments that contribute to wear in the screening of wear loads includes: Filter the time intervals in the contact force time history information of wear loads to find the time intervals in which the contact force appears, and extract the contact force and relative displacement information within the time interval; The time-history intervals of contact force and relative displacement information are spliced ​​together to form an equivalent external load for wear calculation.

[0010] Furthermore, the method for calculating the nodal wear depth of a representative wear period based on the wear depth formula, contact pressure, and relative sliding displacement, assuming that the wear power within a wear morphology update cycle is the same as the wear power within the representative wear period, multiplying the nodal wear depth within the representative wear period by the number of representative wear periods to obtain the accelerated wear result, and updating the component surface model shape and finite element mesh based on the wear result includes: Based on engineering experience and wear conditions, the design life is divided into N wear morphology update cycles. In each wear morphology update cycle, the shape of the component surface model and the finite element mesh remain unchanged, and the average wear power is used to replace the instantaneous wear power. One wear morphology update cycle is divided into n representative wear periods, where n is an integer greater than 1. The wear power in one wear morphology update cycle is set to be the same as the wear power in the representative wear periods. The wear depth of the nodes in the representative wear periods is multiplied by n to obtain the wear result in one wear morphology update cycle after acceleration. Then, the surface model shape of the component and the finite element mesh are updated according to the total wear depth in a single wear morphology update cycle, and the wear calculation for the next wear morphology update cycle begins.

[0011] Furthermore, the specific method for reducing the three-dimensional finite element model to a two-dimensional wear finite element model includes: Extract the geometric features of the three-dimensional model, reduce the three-dimensional geometric model to a two-dimensional geometric model, perform finite element meshing on the two-dimensional geometric model, refine the local mesh in the contact area, and form a two-dimensional wear finite element model. Based on the principle that the maximum contact pressure in the contact area is equal, the resultant contact force in the three-dimensional finite element model is transformed into an equivalent force applied to the two-dimensional wear finite element model.

[0012] Secondly, another embodiment of the present invention provides an acceleration system for predicting fuel cladding wear, comprising: a model building module, a model dimensionality reduction module, a load equivalence processing module, an external load application module, and a wear analysis acceleration module; The model building module is used to build a refined three-dimensional finite element model of the contact area of ​​the wear region of the component. The model dimensionality reduction module preserves geometric features and reduces the three-dimensional finite element model to a two-dimensional wear finite element model. Based on the principle that the maximum contact pressure in the contact area is equal, the contact resultant force in the three-dimensional finite element model is transformed into an equivalent force applied to the two-dimensional wear finite element model. The load equivalence processing module performs equivalent processing on the load, filters out load segments that contribute to wear in the wear load, splices the load segments, and obtains the equivalent external load for wear calculation. The external load application module controls the movement of each part by applying tangential displacement boundary conditions to simulate the relative sliding displacement caused by external loads. It controls the squeezing interaction of each part by applying normal force load conditions to simulate the contact surface squeezing caused by external loads and obtain the contact pressure. The wear analysis acceleration module calculates the node wear depth of a representative wear period based on the wear depth formula, contact pressure, and relative sliding displacement. It assumes that the wear power within a wear morphology update cycle is the same as the wear power within the representative wear period. The node wear depth within the representative wear period is multiplied by the number of representative wear periods to obtain the accelerated wear result. The surface model shape and finite element mesh of the component are updated based on the wear result.

[0013] Furthermore, the model building module includes a mesh generation unit and a setting unit. The mesh generation unit uses a mesh generation tool to perform finite element mesh generation on the geometric model and performs local mesh refinement in the contact area. The setting unit is used to set contact pairs and adaptive meshes on the surface of the pre-wear area, which includes areas where contact and fretting wear may occur.

[0014] Furthermore, the external load application module includes a screening unit and a splicing unit. The screening unit is used to screen the time interval in the contact force time history information of the wear load where the contact force appears, and extract the contact force and relative displacement information within the time interval. The splicing unit is used to splice together the time history intervals of contact force and relative displacement information to form an external load for analysis.

[0015] Furthermore, the wear analysis acceleration module includes a grouping unit and a calculation unit. The grouping unit is used to divide the design life into N wear morphology update cycles. In each wear morphology update cycle, the shape of the component surface model and the finite element mesh remain unchanged, and the instantaneous wear power is replaced by the average wear power. The calculation unit is used to divide one wear morphology update cycle into n representative wear periods, where n is an integer greater than 1. It is set that the wear power in one wear morphology update cycle is the same as the wear power in the representative wear period. The wear depth of the nodes in the representative wear period is multiplied by n to obtain the wear result in one wear morphology update cycle after acceleration. Then, the surface model shape of the component and the finite element mesh are updated according to the total wear depth in a single wear morphology update cycle, and the wear calculation of the next wear morphology update cycle begins.

[0016] Furthermore, the model dimensionality reduction module includes a geometric model dimensionality reduction unit and a load equivalent unit. The geometric model dimensionality reduction unit extracts the geometric features of the three-dimensional geometric model, reduces the three-dimensional geometric model to a two-dimensional geometric model, performs finite element meshing on the two-dimensional geometric model, and refines the local mesh in the contact area to form a two-dimensional wear finite element model. The load equivalent element is based on the principle that the maximum contact pressure in the contact area is equal, and transforms the resultant contact force in the three-dimensional finite element model into an equivalent force applied to the two-dimensional wear finite element model.

[0017] Thirdly, another embodiment of the present invention provides an electronic device including a processor, an input device, an output device, and a memory. The processor is connected to the input device, the output device, and the memory, respectively. The memory is used to store a computer program, the computer program including program instructions, and the processor is configured to call the program instructions to execute the method described in the first embodiment.

[0018] Fourthly, another embodiment of the present invention provides a computer-readable storage medium storing a computer program, the computer program including program instructions that, when executed by a processor, cause the processor to perform the method described in the first embodiment.

[0019] Compared with the prior art, the present invention has the following advantages and beneficial effects: The present invention provides an accelerated method for predicting fuel cladding wear, offering an efficient and simple method to address the engineering problem of the dynamic evolution of long-cycle fretting wear of fuel cladding. This invention proposes model equivalence and dimensionality reduction techniques. Based on the principle of identical contact pressure, it can reduce the dimensionality of a three-dimensional finite element model to a two-dimensional wear finite element model while ensuring analytical accuracy, significantly reducing computational resource consumption. This invention combines short-cycle fretting wear data with nonlinear impact simulation, dividing the design life into multiple groups of equal duration. Within each group, the average wear power replaces the instantaneous wear power, significantly improving the update efficiency of the contact surface geometric model in wear numerical simulation. It balances computational efficiency and analytical accuracy, achieving accurate simulation of fretting wear. This invention identifies, filters, and splices external loads, eliminating invalid times where no contact occurs while retaining statistical information of vibration loads during contact. It performs equivalent processing on the service condition vibration loads that cause wear, significantly improving the information density in the computational data. The wear prediction acceleration method provided by this invention is not limited by structure or support form, and can perform fretting wear analysis on contact surfaces of arbitrary shapes, obtaining the surface wear morphology of components at any given time. The acceleration system, device and medium for predicting fuel cladding wear provided in this embodiment of the invention are based on the same inventive concept and have the same beneficial effects as the acceleration method for predicting fuel cladding wear, and will not be described again here. Attached Figure Description

[0020] To more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the accompanying drawings used in the embodiments will be briefly described below. It should be understood that the following drawings only show some embodiments of the present invention and should not be considered as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort. In the drawings: Figure 1 This is a comparison chart of actual wear marks and ideal wear marks. Figure 2 A flowchart of an accelerated method for predicting fuel cladding wear, provided in the first embodiment of the present invention; Figure 3 A schematic diagram of a classic wear case from existing literature; Figure 4 A schematic diagram of external loads that cause wear; Figure 5 This is a three-dimensional finite element model diagram in the first embodiment of the present invention; Figure 6 This is a diagram showing the correspondence between the contact resultant force of the three-dimensional model and the equivalent force of the two-dimensional model in the first embodiment of the present invention. Figure 7This is a schematic diagram of the load before and after splicing in the first embodiment of the present invention; Figure 8 This is a schematic diagram of time grouping in the first embodiment of the present invention; Figure 9 This is a comparison chart of the numerical simulation results of the wear surface morphology obtained by the method of the first embodiment of the present invention and the results of classic wear case studies; Figure 10 This is a structural block diagram of an acceleration system for predicting fuel cladding wear, provided as another embodiment of the present invention. Detailed Implementation

[0021] To make the objectives, technical solutions, and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the embodiments and accompanying drawings. The illustrative embodiments and descriptions of the present invention are only used to explain the present invention and are not intended to limit the present invention.

[0022] Example 1 like Figure 2 As shown, the first embodiment of the present invention provides an accelerated method for predicting fuel cladding wear, comprising: establishing a refined three-dimensional finite element model of the contact area of ​​the component wear region; preserving geometric features, reducing the three-dimensional finite element model to a two-dimensional wear finite element model; based on the principle that the maximum contact pressure in the contact region is equal, converting the contact resultant force in the three-dimensional finite element model into an equivalent force applied to the two-dimensional wear finite element model; performing equivalent processing on the load, screening the load segments that contribute to wear in the wear load, splicing the load segments to obtain the equivalent external load for wear calculation; controlling the movement of each part by applying tangential displacement boundary conditions to simulate the relative sliding displacement caused by the external load, and controlling the squeezing interaction of each part by applying normal force load conditions to simulate the contact surface squeezing caused by the external load to obtain the contact pressure; The nodal wear depth is calculated based on the wear depth formula, contact pressure, and relative sliding displacement for a representative wear period. Assuming that the wear power is the same as the wear power in the representative wear period within a wear morphology update cycle, the nodal wear depth in the representative wear period is multiplied by the number of representative wear periods to obtain the accelerated wear result. The surface model shape and finite element mesh of the component are updated based on the wear result.

[0023] Specifically, this invention embodiment is illustrated using a simulation of the classic wear case in the literature "Mccoll IR, Ding J, Leen S B. Finite element simulation and experimental validation of fretting wear[J]. Wear, 2004, 256(11):1114-1127" using the commercial finite element software Abaqus, specifically the fretting wear caused by the sliding between the arc-shaped cylindrical indenter and the planar slider. Figure 3 As shown, a semi-cylindrical indenter and a square slider are pressed together to create a contact pair of curved and flat surfaces. Then, the two parts are allowed to slide relative to each other. The pressing force initially increases linearly from zero, reaching 185N and then remaining constant. At this point, the two parts begin to slide periodically relative to each other. Load data and sliding data are... Figure 4 The study presents the interaction between components under loads of this characteristic, simulates the wear process, and investigates the dynamic evolution of wear depth in the contact area over a long period of time.

[0024] In order to better understand the above technical solution, the above technical solution will be described in detail below with reference to the accompanying drawings and specific implementation methods.

[0025] S10: Use a meshing tool to perform finite element meshing on the geometric model, refine the mesh locally in the contact area, and establish a locally refined 3D finite element model of the contact area. The 3D finite element model is shown in the figure below. Figure 5 As shown. The area where contact and fretting wear may occur is called the pre-wear zone (in Figure 5 (Marked in red in the image) On the surface of the pre-wear zone, the tangential mesh size is approximately 0.01 mm to ensure that the mesh size condition for accurate calculation of contact pressure is met. Contact pairs are set for the surface of the pre-wear zone, and adaptive mesh settings are applied to the pre-wear zone.

[0026] After step S10, the method of this embodiment of the invention proceeds to step S20: model dimensionality reduction and load equivalence. The geometric features of the three-dimensional model are preserved: the radius of the arc-shaped cylindrical indenter. A two-dimensional geometric model is established with a semicircle and a rectangular plane in contact. A mesh generation tool is used to perform finite element mesh generation on the two-dimensional geometric model, and local mesh refinement is performed in the contact area to form a two-dimensional finite element model.

[0027] Based on the principle that the maximum contact pressure in the contact area is equal, the contact resultant force in the three-dimensional finite element model is transformed into an equivalent force applied to the two-dimensional wear finite element model: a force load is applied to the three-dimensional model of the contact area to cause mutual compression, the contact pressure distribution in the contact area is analyzed, and the maximum contact pressure at the center of the area under different compression force loads is extracted. A compression load is applied to the two-dimensional model, and the maximum contact pressure is extracted. When the maximum contact pressure in the three-dimensional model and the two-dimensional model are equal, the equivalent force load in the two-dimensional model is considered to have been obtained. The equivalent forces in the two-dimensional model corresponding to different contact resultant forces in the three-dimensional model are plotted on the same graph, and these correspondences are fitted to obtain the correspondence between the equivalent force load in the two-dimensional model and the actual contact resultant force in the three-dimensional model. Figure 6 As shown.

[0028] After step S20, the method of this embodiment of the invention proceeds to step S30: load and displacement splicing. A total duration of T is defined as the representative wear period, where T is an integer greater than 1. The contact force and relative displacement information within this period are considered representative of structural wear. Time intervals showing significant contact force in the contact force time history information of the wear load are filtered out and considered as time intervals where contact and fretting wear processes occur. Contact force and relative displacement information within these time intervals are extracted, and the time history intervals of these contact force and relative displacement information are spliced ​​together. This operation filters out "blank times" in the representative wear period of total duration T that do not involve contact and contribute no to wear, shortening the calculation time of the representative wear period to t, where t is an integer greater than 1. A schematic diagram of the external load causing wear in this embodiment of the invention is shown below. Figure 4 As shown, Figure 7 This is a schematic diagram of the load before and after splicing.

[0029] After step S30, the method of this embodiment proceeds to step S40: applying an external load that causes wear. In the finite element software, the motion of each part is controlled by applying tangential displacement boundary conditions to simulate the relative sliding caused by the external load; the squeezing interaction of each part is controlled by applying normal force load conditions to simulate the contact surface squeezing caused by the external load, thus obtaining the contact pressure.

[0030] After step S40, the method of this embodiment proceeds to step S50: accelerated fretting wear analysis. The wear depth of the node within a representative wear period is calculated based on the wear depth formula, contact pressure, and relative sliding displacement. Since the structural dynamic response and energy dissipation that cause wear are highly similar within a certain time period, and the support clearance and wear power remain essentially constant, the design life is divided into N wear morphology update cycles, where N is an integer greater than 1. One wear morphology update cycle is divided into n representative wear periods, where n is an integer greater than 1. The average wear power within each wear morphology update cycle replaces the instantaneous wear power. That is, within a certain time range, the wear depth value caused by fretting reciprocation multiplied by n can represent the total wear depth caused by n similar fretting reciprocations. In other words, the total wear depth within a single wear morphology update cycle can be obtained by multiplying n by the calculated wear depth of the representative wear period. Then, the component surface model shape and finite element mesh are updated based on the total wear depth within a single wear morphology update cycle, and wear calculation for the next wear morphology update cycle begins.

[0031] Figure 8 The design life is divided into N wear morphology update cycles. Only wear calculations of duration t are required in each wear morphology update cycle. Therefore, the wear calculation time for the entire design life is reduced to N×t.

[0032] The wear depth formula is derived from the Archard wear formula:

[0033] In the formula: W This represents the local wear depth, expressed in meters (m). s The sliding distance is in meters (m). K The wear coefficient is expressed in Pa. -1 Determined by experiments; p This is the local contact pressure, expressed in Pa. t The time unit is seconds (s).

[0034] After step S50, the method of this embodiment of the invention proceeds to step S60: generating an Abaqus computation inp file and running an Abaqus computation job associated with the user subroutine UMESHMOTION.

[0035] The numerical simulation results of the wear surface morphology obtained by the method used in the embodiments of the present invention are compared with the results of the classic wear case in the literature "MccollI R, Ding J, Leen S B. Finite element simulation and experimental validation of fretting wear[J]. Wear, 2004, 256(11):1114-1127". The comparison figure is shown in the figure below. Figure 9 As shown in the figure, the method used in this embodiment of the invention yielded simulation results that are very close to the simulation data in the literature. The numerical simulation results demonstrate that the method of this embodiment can accelerate the calculation of long-term, deep-wear numerical simulations.

[0036] Therefore, while ensuring the authenticity and validity of the numerical simulation results, the method of this invention can perform long-term, deep wear numerical simulations on surfaces with arbitrary morphology. Furthermore, it can be modularly combined with commercial finite element software, making it more applicable and easier to promote and use.

[0037] This invention provides an accelerated method for predicting fuel cladding wear. It merges the long-cycle wear problem of structures operating under similar statistical load conditions into a short-term, continuous, multi-rate wear problem. A numerical simulation acceleration strategy is proposed to accelerate the dynamic evolution of wear morphology in long-term wear problems, enabling the acquisition of the dynamic wear evolution process of the structure and significantly saving computation time and resources used for fretting wear. Compared to methods that estimate wear volume through dynamic response based on the Archard wear volume formula and wear work, the method provided in this invention effectively obtains more critical wear depth characteristics for engineering applications. It proposes an equivalent dimensionality reduction method, transforming the three-dimensional problem into a two-dimensional wear problem, significantly reducing computational resource consumption. It divides the design life into multiple equal-duration groups, replacing instantaneous wear power with average wear power within each group, greatly improving the update efficiency of the contact surface geometry model. It performs equivalent processing on the vibration loads of the service conditions that cause wear, significantly increasing the information density in the computational data. It enables wear simulation under long-term service conditions, significantly improving the analytical accuracy of wear damage numerical simulation in engineering practice, and is more suitable for practical engineering applications.

[0038] Example 2 like Figure 10 As shown, another embodiment of the present invention provides an acceleration system for predicting fuel cladding wear, which corresponds one-to-one with the acceleration method for predicting fuel cladding wear in Embodiment 1. The system includes: a model building module, a model dimensionality reduction module, a load equivalence processing module, an external load application module, and a wear analysis acceleration module. The model building module is used to build a refined three-dimensional finite element model of the contact area of ​​the wear region of the component. The model dimensionality reduction module preserves geometric features and reduces the three-dimensional finite element model to a two-dimensional wear finite element model. Based on the principle that the maximum contact pressure in the contact area is equal, the contact resultant force in the three-dimensional finite element model is transformed into an equivalent force applied to the two-dimensional wear finite element model. The load equivalence processing module performs equivalent processing on the load, filters out load segments that contribute to wear in the wear load, splices the load segments, and obtains the equivalent external load for wear calculation. The external load application module controls the movement of each part by applying tangential displacement boundary conditions to simulate the relative sliding displacement caused by external loads. It controls the squeezing interaction of each part by applying normal force load conditions to simulate the contact surface squeezing caused by external loads and obtain the contact pressure. The wear analysis acceleration module calculates the node wear depth of a representative wear period based on the wear depth formula, contact pressure, and relative sliding displacement. It assumes that the wear power within a wear morphology update cycle is the same as the wear power within the representative wear period. The node wear depth within the representative wear period is multiplied by the number of representative wear periods to obtain the accelerated wear result. The surface model shape and finite element mesh of the component are updated based on the wear result.

[0039] In this embodiment, the model building module includes a mesh generation unit and a setting unit. The mesh generation unit uses a mesh generation tool to perform finite element mesh generation on the geometric model and performs local mesh refinement in the contact area. The setting unit is used to set contact pairs and adaptive meshes on the surface of the pre-wear area. The pre-wear area includes areas where contact and fretting wear may occur.

[0040] In this embodiment, the external load application module includes a screening unit and a splicing unit. The screening unit is used to screen the time intervals in the contact force time history information of the wear load where the contact force appears, and extract the contact force and relative displacement information within the time interval. The splicing unit is used to splice the time history intervals of the contact force and relative displacement information together to form an external load for analysis.

[0041] In this embodiment, the wear analysis acceleration module includes a grouping unit and a calculation unit. The grouping unit is used to divide the design life into N wear morphology update cycles. In each wear morphology update cycle, the shape of the component surface model and the finite element mesh remain unchanged, and the average wear power is used to replace the instantaneous wear power. The calculation unit is used to divide one wear morphology update cycle into n representative wear periods, where n is an integer greater than 1. It is assumed that the wear power in one wear morphology update cycle is the same as the wear power in the representative wear period. The wear depth of the nodes in the representative wear period is multiplied by n to obtain the wear result in one accelerated wear morphology update cycle. Then, the shape of the component surface model and the finite element mesh are updated according to the total wear depth in a single wear morphology update cycle, and the wear calculation for the next wear morphology update cycle begins.

[0042] In this embodiment, the model dimensionality reduction module includes a geometric model dimensionality reduction unit and a load equivalent unit. The geometric model dimensionality reduction unit extracts the geometric features of the three-dimensional geometric model, reduces the three-dimensional geometric model to a two-dimensional geometric model, performs finite element meshing on the two-dimensional geometric model, and refines the mesh locally in the contact area to form a two-dimensional wear finite element model. The load equivalent unit, based on the principle that the maximum contact pressure in the contact area is equal, transforms the resultant contact force in the three-dimensional finite element model into an equivalent force applied to the two-dimensional wear finite element model.

[0043] The execution process of each module can be carried out according to the steps of the acceleration method for fuel cladding wear prediction in Example 1, and will not be described in detail in this example.

[0044] The accelerated system for fuel cladding wear prediction provided in this invention can effectively obtain more critical wear depth characteristics for engineering applications; it proposes an equivalent dimensionality reduction method to transform the three-dimensional problem into a two-dimensional wear problem, significantly reducing computational resource consumption; it divides the design life into multiple groups of equal duration, and replaces the instantaneous wear power with the average wear power in each group, greatly improving the update efficiency of the contact surface geometric model; it performs equivalent processing on the service condition vibration loads that cause wear, significantly improving the information density in the computational data; it can realize wear simulation under long-term service conditions, significantly improving the analytical accuracy of wear damage numerical simulation in engineering practice, and is more suitable for practical engineering use.

[0045] Example 3 Another embodiment of the present invention provides a structural block diagram of an electronic device, the device including a processor, an input device, an output device, and a memory. The processor is connected to the input device, the output device, and the memory, respectively. The memory is used to store a computer program, the computer program including program instructions, and the processor is configured to call the program instructions to execute the method described in the first embodiment above.

[0046] It should be understood that, in the embodiments of the present invention, the processor may be a Central Processing Unit (CPU), but it may also be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor.

[0047] Input devices may include touchpads, fingerprint sensors (for collecting the user's fingerprint information and fingerprint orientation information), microphones, etc., while output devices may include displays (LCDs, etc.), speakers, etc.

[0048] The memory may include read-only memory and random access memory, and provides instructions and data to the processor. A portion of the memory may also include non-volatile random access memory. For example, the memory may also store information about the device type.

[0049] In specific implementations, the processor, input device, and output device described in the embodiments of the present invention can execute the implementation methods described in the method embodiments of the present invention, or they can execute the implementation methods described in the system embodiments of the present invention, which will not be repeated here.

[0050] Example 4 In another embodiment of the present invention, a computer-readable storage medium is also provided, which stores a computer program, the computer program including program instructions that, when executed by a processor, cause the processor to perform the method described in the first embodiment above.

[0051] The computer-readable storage medium can be an internal storage unit of the terminal described in the foregoing embodiments, such as the terminal's hard drive or memory. The computer-readable storage medium can also be an external storage device of the terminal, such as a plug-in hard drive, Smart Media Card (SMC), Secure Digital (SD) card, or Flash Card equipped on the terminal. Furthermore, the computer-readable storage medium can include both internal storage units and external storage devices of the terminal. The computer-readable storage medium is used to store the computer program and other programs and data required by the terminal. The computer-readable storage medium can also be used to temporarily store data that has been output or will be output.

[0052] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0053] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the terminals and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0054] In the several embodiments provided in this application, it should be understood that the disclosed systems and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed may be indirect coupling or communication connection through some interfaces, devices or units, or may be electrical, mechanical or other forms of connection.

[0055] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. An accelerated method for predicting fuel cladding wear, characterized in that, include: Establish a detailed three-dimensional finite element model of the contact area in the wear region of the component; Preserving geometric features, the three-dimensional finite element model is reduced to a two-dimensional wear finite element model. Based on the principle that the maximum contact pressure in the contact area is equal, the contact resultant force in the three-dimensional finite element model is transformed into an equivalent force applied to the two-dimensional wear finite element model. The load is processed by equivalent treatment, and the load segments that contribute to wear are selected. The load segments are spliced ​​together to obtain the equivalent external load for wear calculation. By applying tangential displacement boundary conditions to control the motion of each part, the relative sliding displacement caused by external loads is simulated. By applying normal force load conditions to control the squeezing interaction of each part, the contact surface squeezing caused by external loads is simulated, and the contact pressure is obtained. The nodal wear depth is calculated based on the wear depth formula, contact pressure, and relative sliding displacement for a representative wear period. Assuming that the wear power is the same as the wear power in the representative wear period within a wear morphology update cycle, the nodal wear depth in the representative wear period is multiplied by the number of representative wear periods to obtain the accelerated wear result. The surface model shape and finite element mesh of the component are updated based on the wear result.

2. The method according to claim 1, characterized in that, The specific method for establishing a refined three-dimensional finite element model of the contact area of ​​the component wear region includes: The geometric model is meshed using a mesh generation tool, and the mesh is refined locally in the contact area. Contact pairs and adaptive meshes are provided on the surface of the pre-wear area, which includes areas where contact and fretting wear may occur.

3. The method according to claim 1, characterized in that, The specific method for splicing the load segments that contribute to wear in the screening of wear loads includes: Filter the time intervals in the contact force time history information of wear loads to find the time intervals in which the contact force appears, and extract the contact force and relative displacement information within the time interval; The time-history intervals of contact force and relative displacement information are spliced ​​together to form an equivalent external load for wear calculation.

4. The method according to claim 1, characterized in that, The method for calculating the nodal wear depth of a representative wear period based on the wear depth formula, contact pressure, and relative sliding displacement, assuming that the wear power within a wear morphology update cycle is the same as the wear power within the representative wear period, multiplying the nodal wear depth within the representative wear period by the number of representative wear periods to obtain the accelerated wear result, and updating the component surface model shape and finite element mesh based on the wear result includes: Based on engineering experience and wear conditions, the design life is divided into N wear morphology update cycles. In each wear morphology update cycle, the shape of the component surface model and the finite element mesh remain unchanged, and the average wear power is used to replace the instantaneous wear power. One wear morphology update cycle is divided into n representative wear periods, where n is an integer greater than 1. The wear power in one wear morphology update cycle is set to be the same as the wear power in the representative wear periods. The wear depth of the nodes in the representative wear periods is multiplied by n to obtain the wear result in one wear morphology update cycle after acceleration. Then, the surface model shape of the component and the finite element mesh are updated according to the total wear depth in a single wear morphology update cycle, and the wear calculation for the next wear morphology update cycle begins.

5. The method according to claim 1, characterized in that, The specific methods for reducing a three-dimensional finite element model to a two-dimensional wear finite element model include: Extract the geometric features of the three-dimensional model, reduce the three-dimensional geometric model to a two-dimensional geometric model, perform finite element meshing on the two-dimensional geometric model, refine the local mesh in the contact area, and form a two-dimensional wear finite element model. Based on the principle that the maximum contact pressure in the contact area is equal, the resultant contact force in the three-dimensional finite element model is transformed into an equivalent force applied to the two-dimensional wear finite element model.

6. An acceleration system for predicting fuel cladding wear, characterized in that, include: The module includes a model building module, a model dimensionality reduction module, a load equivalence processing module, an external load application module, and a wear analysis acceleration module. The model building module is used to build a refined three-dimensional finite element model of the contact area of ​​the wear region of the component. The model dimensionality reduction module preserves geometric features and reduces the three-dimensional finite element model to a two-dimensional wear finite element model. Based on the principle that the maximum contact pressure in the contact area is equal, the contact resultant force in the three-dimensional finite element model is transformed into an equivalent force applied to the two-dimensional wear finite element model. The load equivalence processing module performs equivalent processing on the load, filters out load segments that contribute to wear in the wear load, splices the load segments, and obtains the equivalent external load for wear calculation. The external load application module controls the movement of each part by applying tangential displacement boundary conditions to simulate the relative sliding displacement caused by external loads. It controls the squeezing interaction of each part by applying normal force load conditions to simulate the contact surface squeezing caused by external loads and obtain the contact pressure. The wear analysis acceleration module calculates the node wear depth of a representative wear period based on the wear depth formula, contact pressure, and relative sliding displacement. It assumes that the wear power within a wear morphology update cycle is the same as the wear power within the representative wear period. The node wear depth within the representative wear period is multiplied by the number of representative wear periods to obtain the accelerated wear result. The surface model shape and finite element mesh of the component are updated based on the wear result.

7. The system as described in claim 6, characterized in that, The model building module includes a mesh generation unit and a setting unit. The mesh generation unit uses a mesh generation tool to perform finite element mesh generation on the geometric model and performs local mesh refinement in the contact area. The setting unit is used to set contact pairs and adaptive meshes on the surface of the pre-wear area, which includes areas where contact and fretting wear may occur.

8. The system as described in claim 6, characterized in that, The external load application module includes a screening unit and a splicing unit. The screening unit is used to screen the time interval in the contact force time history information of the wear load where the contact force appears, and extract the contact force and relative displacement information within the time interval. The splicing unit is used to splice together the time history intervals of contact force and relative displacement information to form an external load for analysis.

9. The system as described in claim 6, characterized in that, The wear analysis acceleration module includes a grouping unit and a calculation unit. The grouping unit is used to divide the design life into N wear morphology update cycles. In each wear morphology update cycle, the shape of the component surface model and the finite element mesh remain unchanged, and the instantaneous wear power is replaced by the average wear power. The calculation unit is used to divide one wear morphology update cycle into n representative wear periods, where n is an integer greater than 1. It is set that the wear power in one wear morphology update cycle is the same as the wear power in the representative wear period. The wear depth of the nodes in the representative wear period is multiplied by n to obtain the wear result in one wear morphology update cycle after acceleration. Then, the surface model shape of the component and the finite element mesh are updated according to the total wear depth in a single wear morphology update cycle, and the wear calculation of the next wear morphology update cycle begins.

10. The system as described in claim 7, characterized in that, The model dimensionality reduction module includes a geometric model dimensionality reduction unit and a load equivalent unit. The geometric model dimensionality reduction unit extracts the geometric features of the three-dimensional geometric model, reduces the three-dimensional geometric model to a two-dimensional geometric model, performs finite element meshing on the two-dimensional geometric model, and refines the local mesh in the contact area to form a two-dimensional wear finite element model. The load equivalent element is based on the principle that the maximum contact pressure in the contact area is equal, and transforms the resultant contact force in the three-dimensional finite element model into an equivalent force applied to the two-dimensional wear finite element model.

11. An electronic device comprising a processor, an input device, an output device, and a memory, wherein the processor is connected to the input device, the output device, and the memory, and the memory is used to store a computer program, the computer program comprising program instructions, characterized in that... The processor is configured to invoke the program instructions to perform the method as described in any one of claims 1-5.

12. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program, the computer program including program instructions that, when executed by a processor, cause the processor to perform the method as described in any one of claims 1-5.

Images