A life evaluation method and device of a high-voltage circuit breaker based on field-circuit coupling, a terminal device and a computer readable storage medium

Abstract

The application discloses a kind of life evaluation method, device, terminal equipment and computer readable storage medium of high-voltage circuit breaker based on field-circuit coupling, belong to high-voltage switchgear simulation and reliability evaluation technical field.Method includes: fluid simulation is carried out to high-voltage circuit breaker geometric model, and gas pressure distribution data is calculated;Identify the force interface of transmission component not parallel to the direction of motion, extract radial profile coordinate data to calculate annular effective projection area, and map the pressure data to the area to calculate dynamic load reaction data;The reaction data is imported into the multi-body dynamics model to calculate kinematic state data;Based on kinematic state and reaction data, calculate mechanism output power data and component stress distribution law, to evaluate mechanical life attenuation condition.The application solves the problem that existing method is inaccurate in evaluating mechanical life of high-voltage circuit breaker due to distortion in extracting dynamic load reaction.

Description

Technical Field

[0001] This invention relates to the field of high-voltage switchgear simulation and reliability assessment technology, and in particular to a method, apparatus, terminal equipment, and computer-readable storage medium for assessing the lifespan of a high-voltage circuit breaker based on field-circuit coupling. Background Technology

[0002] Co-simulation technology based on field-circuit coupling has been widely used in the mechanical characteristic analysis of the opening and closing process of high-voltage circuit breakers; the accurate extraction of the dynamic load reaction force generated by the gas in the arc-extinguishing chamber of the circuit breaker on the solid transmission components will directly affect the accurate calculation of the output power of the hydraulic operating mechanism and the stress distribution law of key components.

[0003] Existing force analysis methods, when dealing with complex geometric models of transmission components, often directly integrate pressure across the entire surface of the component, failing to identify and distinguish the effectiveness of different force-bearing surfaces. This leads to distorted calculated load reaction force data, failing to accurately reflect the actual reaction force of gas hindering the movement of the transmission component. Consequently, the calculated results of kinematic states, output work, and stress distribution in multibody dynamics simulations are inaccurate. Therefore, improving the accuracy of dynamic load reaction force extraction and output work calculation for high-voltage circuit breakers has become a pressing technical problem to be solved in this field. Summary of the Invention

[0004] This invention provides a life assessment method for high-voltage circuit breakers based on field-circuit coupling, which can solve the problem of inaccurate mechanical life assessment of high-voltage circuit breakers caused by distortion in the extraction of dynamic load reaction force in existing methods.

[0005] One embodiment of the present invention provides a life assessment method for high-voltage circuit breakers based on field-circuit coupling, comprising: Geometric model data of a high-voltage circuit breaker is acquired, and fluid dynamics simulation is performed on the geometric model data to calculate the gas pressure distribution data of the gas inside the arc-extinguishing chamber acting on the solid-state transmission component during the opening and closing process of the high-voltage circuit breaker. The force interface on the surface of the solid-state transmission component that is not parallel to the direction of movement of the solid-state transmission component is identified, and the radial contour coordinate data of the force interface is extracted. Based on the radial contour coordinate data, the annular effective projected area of ​​the force interface in the direction of movement is calculated. Then, the gas pressure distribution data is mapped onto the annular effective projected area to calculate the dynamic load reaction force data of the gas hindering the movement of the solid-state transmission component. The dynamic load reaction force data is used as an external load condition and imported into a preset multibody dynamics simulation model of the transmission system. The kinematic state data of the solid transmission component under the condition of overcoming the dynamic load reaction force is calculated through the multibody dynamics simulation model of the transmission system. Based on the kinematic state data and the dynamic load reaction force data, the output power data of the hydraulic operating mechanism and the stress distribution law of the solid transmission component are calculated respectively. Based on the stress distribution law of the solid transmission component, the mechanical life decay of the high-voltage circuit breaker is evaluated.

[0006] Furthermore, fluid dynamics simulation is performed on the geometric model data to calculate the gas pressure distribution data of the gas inside the arc-extinguishing chamber acting on the solid-state transmission component during the opening and closing process of the high-voltage circuit breaker. This includes: constructing a two-dimensional axisymmetric magnetohydrodynamic model based on the geometric model data; acquiring the short-circuit current data of the high-voltage circuit breaker during the breaking process, and calculating the Joule heating and radiative energy loss of the arc plasma based on the short-circuit current data; substituting the Joule heating and radiative energy loss as energy source terms into the energy conservation equation of the two-dimensional axisymmetric magnetohydrodynamic model; and solving the energy conservation equation to obtain the gas pressure distribution data of the gas inside the arc-extinguishing chamber acting on the solid-state transmission component.

[0007] Further, the radial contour coordinate data of the force-bearing interface is extracted, and the annular effective projected area of ​​the force-bearing interface in the direction of motion is calculated based on the radial contour coordinate data, including: extracting the radial coordinates of the outer edge and the inner edge of the force-bearing interface in the radial direction; calculating the difference between the square of the outer edge radial coordinate and the square of the inner edge radial coordinate; and multiplying the difference by pi as the annular effective projected area of ​​the force-bearing interface in the direction of motion.

[0008] Further, mapping the gas pressure distribution data onto the annular effective projected area and calculating the dynamic load reaction force data of the gas hindering the movement of the solid transmission component includes: extracting the target node gas pressure data corresponding to the force interface from the gas pressure distribution data; multiplying the target node gas pressure data by the annular effective projected area as the partial reaction force data borne by the force interface in the direction of movement; and summing the partial reaction force data of all the force interfaces to obtain the dynamic load reaction force data of the gas hindering the movement of the solid transmission component.

[0009] Furthermore, the kinematic state data of the solid transmission component under the condition of overcoming the dynamic load reaction force are calculated using the multibody dynamics simulation model of the transmission system, including: In the preset simulation platform, the connection pair attributes, friction attributes, and contact constraint attributes between the solid transmission components are configured, and a multibody dynamics simulation model of the transmission system is constructed based on the connection pair attributes, friction attributes, and contact constraint attributes between the solid transmission components. Based on the lumped parameter method, a mathematical model of the hydraulic cylinder is constructed that interacts with the multibody dynamics simulation model of the transmission system; The dynamic load reaction force data is input into the multibody dynamics simulation model of the transmission system, and combined with the system oil pressure data provided by the hydraulic cylinder mathematical model for joint solution to obtain the displacement and velocity data of the solid transmission component during the motion process, which are used as the kinematic state data.

[0010] Further, based on the kinematic state data and the dynamic load reaction force data, the output work data of the hydraulic operating mechanism and the stress distribution law of the solid transmission component are calculated respectively, including: calculating the kinetic energy data of the solid transmission component during the motion process based on the component velocity data and the mass data of the solid transmission component in the kinematic state data; calculating the load work data required for the solid transmission component to overcome gas resistance based on the component displacement data and the corresponding dynamic load reaction force data in the kinematic state data; and combining the kinetic energy data and the load work data to obtain the output work data of the hydraulic operating mechanism driving the solid transmission component. The dynamic load reaction data and the kinematic state data are used as force boundary conditions and input into a pre-constructed transmission system stress analysis model for solution, thereby obtaining the stress distribution law of the solid transmission component.

[0011] Furthermore, based on the stress distribution pattern of the solid-state transmission component, the mechanical life degradation of the high-voltage circuit breaker is evaluated, including: extracting stress extreme value data of key connecting components in the solid-state transmission component from the stress distribution pattern, wherein the key connecting components include at least one of connecting rod, insulating tie rod, and pin; and evaluating the mechanical life degradation of the high-voltage circuit breaker based on preset material fatigue life parameters and the stress extreme value data.

[0012] Another embodiment of the present invention provides a life assessment device for a high-voltage circuit breaker based on field-circuit coupling, comprising: a fluid field simulation module, a load reaction force mapping module, a multibody dynamics solution module, and a comprehensive life assessment module; The fluid field simulation module is used to acquire the geometric model data of the high-voltage circuit breaker and perform fluid dynamics simulation on the geometric model data to calculate the gas pressure distribution data of the gas inside the arc extinguishing chamber acting on the solid transmission component during the opening and closing process of the high-voltage circuit breaker. The load reaction force mapping module is used to identify the force interface on the surface of the solid transmission component that is not parallel to the direction of motion of the solid transmission component, extract the radial contour coordinate data of the force interface, and calculate the annular effective projection area of ​​the force interface in the direction of motion based on the radial contour coordinate data. Then, the gas pressure distribution data is mapped onto the annular effective projection area to calculate the dynamic load reaction force data of the gas that hinders the movement of the solid transmission component. The multibody dynamics solution module is used to import the dynamic load reaction force data as an external load condition into a preset transmission system multibody dynamics simulation model, and calculate the kinematic state data of the solid transmission component under the condition of overcoming the dynamic load reaction force through the transmission system multibody dynamics simulation model. The lifespan comprehensive assessment module is used to calculate the output power data of the hydraulic operating mechanism and the stress distribution law of the solid transmission component based on the kinematic state data and the dynamic load reaction force data, and to assess the mechanical lifespan degradation of the high-voltage circuit breaker based on the stress distribution law of the solid transmission component.

[0013] Another embodiment of the present invention provides a terminal device, including: a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor. When the processor executes the computer program, it implements the life assessment method for a high-voltage circuit breaker based on field-circuit coupling of the present invention.

[0014] Another embodiment of the present invention also provides a computer-readable storage medium item, including: a stored computer program, wherein, when the computer program is running, it controls the device where the computer-readable storage medium is located to execute the life assessment method for a high-voltage circuit breaker based on field-circuit coupling of the present invention.

[0015] The embodiments of the present invention have the following beneficial effects: This invention provides a life assessment method for high-voltage circuit breakers based on field-circuit coupling. By acquiring the geometric model data of the high-voltage circuit breaker and performing fluid dynamics simulation on the geometric model data, the gas pressure distribution data of the gas inside the arc-extinguishing chamber acting on the solid-state transmission component during the opening and closing process of the high-voltage circuit breaker is calculated. The method identifies the force-bearing interfaces on the surface of the solid-state transmission component that are not parallel to the direction of motion of the solid-state transmission component, extracts the radial contour coordinate data of the force-bearing interfaces, and calculates the annular effective projected area of ​​the force-bearing interfaces in the direction of motion based on the radial contour coordinate data. Finally, the gas pressure distribution data is mapped to the annular effective projected area. On the shadow area, the dynamic load reaction force data of the gas hindering the movement of the solid transmission component is calculated; the dynamic load reaction force data is used as an external load condition and imported into a preset multibody dynamics simulation model of the transmission system; and the kinematic state data of the solid transmission component under the condition of overcoming the dynamic load reaction force is calculated through the multibody dynamics simulation model of the transmission system; based on the kinematic state data and the dynamic load reaction force data, the output power data of the hydraulic operating mechanism and the stress distribution law of the solid transmission component are calculated respectively; and the mechanical life decay of the high-voltage circuit breaker is evaluated according to the stress distribution law of the solid transmission component.

[0016] By employing the above technical solution, this invention addresses the problem of inaccurate mechanical life assessment of high-voltage circuit breakers caused by distortion in the extraction of dynamic load reaction forces in existing methods. It identifies force interfaces that are not parallel to the direction of motion and extracts radial contour coordinate data, thus reconstructing the complex surface force features into an effective annular projected area along the direction of motion at the algorithmic geometric level. This removes mutually canceling radial invalid components from a physical computational mechanism perspective. Simultaneously, based on the aforementioned geometric reconstruction features, this invention directionally maps gas pressure distribution data onto this effective annular projected area for calculation, avoiding the error accumulation introduced by traditional spatial vector synthesis calculations, and achieving… High-precision quantitative extraction of dynamic load reaction force data that truly hinders the movement of transmission components is achieved. Based on this, the extracted high-precision dynamic load reaction force data is imported into the multibody dynamics simulation model of the transmission system as an external load condition. This effectively constructs a low-level data closed loop between fluid field resistance and mechanical field response, ensuring that the final calculated kinematic state data, hydraulic operating mechanism output power data, and stress distribution law of solid transmission components can accurately reflect the coupling effect under real physical conditions. This overcomes the interference of reaction force extraction distortion on subsequent evaluation and significantly improves the accuracy and reliability of the assessment of the mechanical life decay of high-voltage circuit breakers. Attached Figure Description

[0017] To more clearly illustrate the technical solution of this application, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0018] Figure 1 This is a schematic flowchart of a life assessment method for a high-voltage circuit breaker based on field-circuit coupling provided by an embodiment of the present invention.

[0019] Figure 2 This is a schematic diagram of the device for a life assessment method for high-voltage circuit breakers based on field-circuit coupling, provided in an embodiment of the present invention. Detailed Implementation

[0020] To make the objectives, technical solutions, and advantages of this application clearer, the technical solutions of this application will be clearly and completely described below with reference to the accompanying drawings of the embodiments. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0021] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains; the terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the application; the terms “comprising” and “having”, and any variations thereof, in the specification, claims, and foregoing description of the drawings are intended to cover non-exclusive inclusion.

[0022] In the description of the embodiments of this application, technical terms such as "first" and "second" are used only to distinguish different objects and should not be construed as indicating or implying relative importance or implicitly specifying the number, specific order, or primary and secondary relationship of the indicated technical features. In the description of the embodiments of this application, "multiple" means two or more, unless otherwise explicitly defined.

[0023] In this document, the term "embodiment" means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of this application. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.

[0024] In the description of the embodiments in this application, the term "and / or" is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, and B existing alone. Additionally, the character " / " in this document generally indicates that the preceding and following related objects have an "or" relationship.

[0025] In the description of the embodiments of this application, the term "multiple" refers to two or more (including two), similarly, "multiple sets" refers to two or more (including two sets), and "multiple pieces" refers to two or more (including two pieces).

[0026] In the description of the embodiments of this application, unless otherwise expressly specified and limited, technical terms such as "installation," "connection," "joining," and "fixing" should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral part; they can refer to a mechanical connection or an electrical connection; they can refer to a direct connection or an indirect connection through an intermediate medium; they can refer to the internal communication of two components or the interaction between two components. For those skilled in the art, the specific meaning of the above terms in the embodiments of this application can be understood according to the specific circumstances.

[0027] It is particularly important to clarify the following basic and explanatory notes regarding the hardware and software environment and execution entity involved in the embodiments of this application: The life assessment method for high-voltage circuit breakers based on field-circuit coupling provided in the embodiments of this application is essentially a data processing and virtual simulation calculation method executed by a computer, server, or related electronic computing device. In the description of the specific embodiments below, for the sake of logical coherence and conciseness, the term "system" (e.g., "system extraction," "system calculation," "system call," etc.) is used in many places as the subject of various instructions, mathematical operations, and data flow actions. This "system" is only a general term for the "execution entity" that executes the algorithm steps of this application. This "system" can refer to any electronic computing device with data processing capabilities (such as a personal computer, high-performance workstation, distributed server cluster, etc.), or it can refer to a CAE co-simulation software platform running on the above-mentioned hardware platform (such as a comprehensive computing platform containing fluid mechanics, multibody dynamics, and finite element analysis modules), and does not refer to a specific physical hardware system, nor does it constitute any limitation on the physical hardware architecture of the scope of protection of this application. Those skilled in the art should clearly understand that all data acquisition, identification, mapping, solving, and evaluation actions performed by the "system" in the following embodiments are essentially implemented by the processor (CPU / GPU) in the electronic computing device reading and running the corresponding computer program code or instructions in the memory.

[0028] To address the problem of inaccurate mechanical life assessment of high-voltage circuit breakers due to distortion in dynamic load reaction force extraction in existing methods, an embodiment of the present invention provides a life assessment method for high-voltage circuit breakers based on field-circuit coupling, comprising: Step S1: Obtain the geometric model data of the high-voltage circuit breaker and perform fluid dynamics simulation on the geometric model data to calculate the gas pressure distribution data of the gas inside the arc extinguishing chamber acting on the solid transmission component during the opening and closing process of the high-voltage circuit breaker. In a preferred embodiment, hydrodynamic simulation is performed on the geometric model data to calculate the gas pressure distribution data of the gas inside the arc-extinguishing chamber acting on the solid-state transmission component during the opening and closing process of the high-voltage circuit breaker. This includes: constructing a two-dimensional axisymmetric magnetohydrodynamic model based on the geometric model data; acquiring the short-circuit current data of the high-voltage circuit breaker during the breaking process, and calculating the Joule heating and radiative energy loss of the arc plasma based on the short-circuit current data; substituting the Joule heating and radiative energy loss as energy source terms into the energy conservation equation of the two-dimensional axisymmetric magnetohydrodynamic model; and solving the energy conservation equation to obtain the gas pressure distribution data of the gas inside the arc-extinguishing chamber acting on the solid-state transmission component.

[0029] Preferably, a two-dimensional axisymmetric magnetohydrodynamic model is constructed based on the geometric model data through the following steps: First, the three-dimensional geometric model data of the arc-extinguishing chamber of the actual high-voltage circuit breaker prototype is obtained. Since the arc-extinguishing chamber structure of the real high-voltage circuit breaker has many non-core geometric features, in order to balance computational accuracy and solution efficiency, this embodiment performs structural simplification preprocessing on the three-dimensional geometric model data: specifically, non-fluid-critical structures such as threads, screw holes, fillets, and fixing bolts are removed from the model; simultaneously, for non-axisymmetric structures such as the air outlet of the compressor cylinder and the pressure relief valve, the flow equivalence principle is used to equate them to axisymmetric annular channels. Based on this, utilizing the central axisymmetric property of the overall structure of the high-voltage circuit breaker, the three-dimensional geometric model is reduced in dimension and extracted into two-dimensional axisymmetric geometric model data. Further, the attributes of each physical region are defined in the two-dimensional axisymmetric geometric model data. Specifically, this includes: defining the moving region containing movable parts (such as the moving arc contactor and piston rod) as the dynamic mesh region; defining the lateral relative sliding interface of the gas domain as the interface (fluid interface), the longitudinal interface of adjacent gas domains as the interior (internal surface), the pressure outlet as the outlet, and the solid surface as the wall (wall boundary). Subsequently, the processed two-dimensional axisymmetric geometric model is meshed, and the dynamic mesh region is set as a quadrilateral mesh. It is particularly noteworthy that, due to the severe temperature and pressure gradients in the arc region (i.e., the core region where the arc is generated), the mesh quality requirements are extremely high. This embodiment pre-limits the upper limit of the mesh size for the key arc-focused region, for example, limiting its maximum mesh size to no more than 2 mm, to ensure the convergence of turbulence and arc energy transfer calculations. Finally, based on the aforementioned mesh model, magnetohydrodynamics (MHD) theory is introduced. The arc plasma generated by the ionization of gas (such as SF6) in the arc-extinguishing chamber at high temperature is considered as a compressible and conductive fluid medium, thereby constructing a two-dimensional axisymmetric magnetohydrodynamic model that includes mass conservation equations, momentum conservation equations (Navier-Stokes equations), and a standard k-epsilon turbulence model. The standard k-epsilon turbulence model considers momentum exchange and viscosity between fluid layers, effectively simulating the high-speed unsteady vortex flow generated by the arc.

[0030] Preferably, the short-circuit current data of the high-voltage circuit breaker during the breaking process is obtained through the following steps, and the Joule heating and radiant energy loss of the arc plasma are calculated based on the short-circuit current data: the "short-circuit current data" refers to pre-configured large current characteristic waveform data (e.g., a 50kA power frequency sinusoidal current waveform containing a DC bias component) flowing through the contact system of the high-voltage circuit breaker when a short-circuit fault occurs in the power grid. In this embodiment, the short-circuit current data is input into the two-dimensional axisymmetric magnetohydrodynamic model in the form of a transient time function using a user-defined function (UDF, written in C language and defined using the DEFINE macro). Regarding the calculation of Joule heating: since the arc plasma is a conductive fluid, a huge Joule heat is generated when the short-circuit current flows through the arc. The model internally solves a simplified Ohm's law equation (…). ), Current continuity equation ( ) and the potential gradient equation ( The spatial current density distribution and electric field intensity within the arc region were calculated (where...). Plasma conductivity, (where is the electric potential, j is the current density, and E is the electric field strength). Further calculations using Ampere's circuital law are performed to determine the azimuth magnetic field distribution, and finally, Joule's law is used to calculate the Joule heat generation at each node in space. .

[0031] Specifically, regarding Joule heating The refined calculation steps are as follows: Since the arc column has strict macroscopic axisymmetric properties in the two-dimensional axisymmetric geometric model, the main component of the electromagnetic field originates from the azimuth direction. To accurately quantify the influence of this magnetic field on the self-magnetic compression effect and energy distribution of the arc, this embodiment assumes that the permeability of the arc medium is a preset parameter. (That is, the product of vacuum permeability and plasma relative permeability), and the magnetic field in the two-dimensional computational domain is solved numerically by means of Ampere's circuital law. Specifically, for any fluid mesh node with radial coordinate r in space, its azimuth magnetic field strength is... By measuring the axial current density It is obtained by integrating along the radial section; its mathematical expression logic is: calculate the axial current density Jz and the area of ​​the infinitesimal element (where the radial coordinate is 0) within the integration range from the central axis of symmetry (i.e., the radial coordinate is 0) to the current node radius r. ,in The integral sum of the product of the radial virtual space variables (in the integration process) is multiplied by the magnetic permeability. Dividing by the circumference of the circle containing that node (2πr), the azimuth magnetic field distribution of the entire flow field is dynamically obtained. The current density vector of each grid node in space is then acquired. (Including axial and radial current density components) and the corresponding electric field intensity vector Then, according to the differential form of Joule's law, the current density vector is... With the electric field intensity vector Perform dot product (i.e.) ), or by dividing the square of the nodal current density amplitude by the current plasma conductivity σ of that node (i.e. This allows for the precise calculation of the volumetric heat source power density of each grid node in space at the current time step, i.e., the Joule heat generation. .

[0032] Regarding the amount of radiation energy loss The detailed calculation steps are as follows: When a high-voltage circuit breaker interrupts a large short-circuit current, the core temperature of the arc plasma can rise sharply to tens of thousands of Kelvin. Thermal radiation becomes the most crucial physical pathway for it to dissipate energy to the surrounding cold air (such as unionized SF6 gas). Considering the large characteristic scale and extremely high plasma concentration of the arc region, its optical thickness is much greater than 1 (i.e., it belongs to a strongly absorbing and strongly emitting participating radiation medium). In this embodiment, the P-1 radiation model is preferred for solution. Using the P-1 radiation model can significantly reduce the computational load of solving the multidimensional integral-differential radiation transfer equation by using a first-order truncation expansion of the complex directional radiation intensity using spherical harmonic functions, while ensuring the accuracy of the radiation energy transfer calculation. Specifically, based on the node thermodynamic temperature T and gas pressure data of the fluid field at the current time step, the pre-built-in arc-extinguishing medium property data mapping table is called to dynamically match and obtain the absorption coefficient a (used to characterize the local medium's ability to absorb radiation energy) and scattering coefficient of each spatial calculation node. (Used to characterize the ability of a medium to change the direction of radiation propagation). Subsequently, the P-1 radiation model constructs a partial differential diffusion equation for the incident spatial radiation energy G. The solver incorporates the temperature gradient of the spatial grid nodes, the absorption coefficient a, and the scattering coefficient. The scalar distribution of incident radiation energy G in the entire computational domain is obtained by matrix iteration.

[0033] Based on this, and using the assumption of local thermodynamic equilibrium, the spontaneous emission power density emitted outward from each node (i.e., ,in Let T be the Stefan-Boltzmann constant and T be the absolute temperature of the node, and let aG be the difference between the actual incident radiation power density absorbed by the node. This difference represents the net radiative heat released (or absorbed) by each node in space to the surrounding environment per unit volume. In this embodiment, it is defined as the radiative energy loss. (Its mathematical expression is:) ).

[0034] Preferably, the energy conservation equation is solved through the following steps to obtain the gas pressure distribution data of the gas inside the arc-extinguishing chamber acting on the solid transmission component: First, the basic material properties of the arc-extinguishing medium (such as SF6 gas) are input through the aforementioned UDF custom function. These properties include density, specific heat at constant pressure, thermal conductivity, viscosity, and velocity of sound, which vary nonlinearly with temperature and pressure. The gas state equation ρ=f(P,T) is defined. Next, the solver parameters are configured: a pressure-based solver is used, the time type is set to transient, the velocity equation is specified as absolute velocity, and global initial conditions are set (e.g., initial temperature is 300K, and initial operating pressure is set to the preset circuit breaker working charging pressure). After the simulation calculation starts, the solver performs transient iterations at preset time steps. Within each time step, the moving mesh region is reconstructed and its boundaries are updated according to the pre-input motion travel curve of the movable component; simultaneously, the solver solves the mass conservation, momentum conservation, k-epsilon turbulence equation, and energy conservation equation that introduces the energy source term S. When the equation residuals at the current time step reach the convergence criterion, the flow field variables of all grid nodes in the computational domain of the arc-extinguishing chamber at the current moment can be extracted. From this, the set of wall pressures that dynamically change over time at the fluid node pressure values ​​located on the wall boundary of the solid transmission component is extracted, which is the required gas pressure distribution data.

[0035] In this embodiment, the magnetohydrodynamic modeling and solution processing for the characteristics of the arc-extinguishing chamber of a high-voltage circuit breaker, as described above, achieves the following significant technical effects: On the one hand, by using flow equivalence and two-dimensional axisymmetric dimensionality reduction, the massive computational load of the three-dimensional mesh is greatly reduced, while satisfying the feasibility of transient calculations in simulation software, and the core characteristics of fluid volume and axial flow are preserved. On the other hand, it breaks through the technical bottleneck of traditional cold-state fluid dynamics simulation being unable to simulate the force of large-current interruption. Through the UDF interface, an arc energy source term including Joule heating and radiation loss is innovatively introduced, and solved under the constraints of a highly accurate moving mesh and a k-epsilon turbulence model, realistically reproducing the physical phenomenon of explosive pressure surge caused by the high temperature of the arc. The gas pressure distribution data obtained from this calculation eliminates pressure distortion caused by neglecting the thermal-fluid coupling effect, providing a data benchmark input with extremely high physical confidence for subsequent accurate separation of radial forces and calculation of the real axial dynamic load reaction force of the transmission system.

[0036] Step S2: Identify the force interface on the surface of the solid transmission component that is not parallel to the direction of motion of the solid transmission component, extract the radial contour coordinate data of the force interface, and calculate the annular effective projected area of ​​the force interface in the direction of motion based on the radial contour coordinate data. Then, map the gas pressure distribution data onto the annular effective projected area to calculate the dynamic load reaction force data of the gas that hinders the movement of the solid transmission component. In a preferred embodiment, the radial contour coordinate data of the force-bearing interface is extracted, and the annular effective projected area of ​​the force-bearing interface in the direction of motion is calculated based on the radial contour coordinate data. This includes: extracting the radial coordinates of the outer edge and the inner edge of the force-bearing interface in the radial direction; calculating the difference between the square of the outer edge radial coordinate and the square of the inner edge radial coordinate; and multiplying the difference by pi as the annular effective projected area of ​​the force-bearing interface in the direction of motion.

[0037] Preferably, mapping the gas pressure distribution data onto the annular effective projected area and calculating the dynamic load reaction force data of the gas hindering the movement of the solid transmission component includes: extracting the target node gas pressure data corresponding to the force interface from the gas pressure distribution data; multiplying the target node gas pressure data by the annular effective projected area as the partial reaction force data borne by the force interface in the direction of movement; and summing the partial reaction force data of all the force interfaces to obtain the dynamic load reaction force data of the gas hindering the movement of the solid transmission component.

[0038] Preferably, the following computer graphics and mesh topology traversal steps are used to identify the force interfaces on the surface of the solid transmission component that are not parallel to the direction of motion of the solid transmission component: First, in the post-processing module of the fluid field solver, a preset direction of motion is defined. This preset direction of motion refers to the direction of macroscopic mechanical displacement (usually defined as the axial direction in a two-dimensional axisymmetric coordinate system, i.e., the Z-axis or X-axis direction) that occurs along the central axis of symmetry during the opening and closing operations of the high-voltage circuit breaker piston rod and the moving arc contact. Subsequently, all wall mesh elements marked as fluid-structure interaction boundaries in the two-dimensional axisymmetric geometric model are extracted. For each wall mesh element, its surface normal vector is calculated. The surface normal vector is multiplied by the direction vector corresponding to the preset direction of motion to obtain the angle between the two. If the included angle is equal to 90 degrees (i.e., the surface normal vector is perpendicular to the direction of motion), then the wall mesh element is determined to be completely parallel to the direction of motion (e.g., the pure cylindrical side surface of a piston rod). Ideally, this surface only experiences radial compression and does not generate axial resistance, so it is removed from the force calculation set. If the included angle is not equal to 90 degrees (i.e., the surface normal vector has a non-zero projection component in the direction of motion), then the wall surface is an inclined surface, curved surface, or vertical step surface. These geometric surfaces will generate axial forces that directly hinder the movement of the component when subjected to high-pressure gas compression. The system clusters and extracts all wall mesh elements that satisfy this included angle condition, marking them as force interfaces that are not parallel to the direction of motion.

[0039] Preferably, the radial contour coordinate data of the stress interface is extracted and the annular effective projected area is calculated through the following geometric mapping and dimensionality reduction steps: For any stress interface identified above (taking a certain inclined section as an example), the system extracts the coordinates of the two endpoints of the section in a two-dimensional axisymmetric coordinate system. The radial coordinate of the endpoint farther from the central axis of symmetry is defined as the radial coordinate of the outer edge (denoted as...). The radial coordinates of the endpoints closer to the central axis of symmetry are defined as the radial coordinates of the inner edge (denoted as ). In three-dimensional physical space, this force-bearing interface is the side surface of a frustum or a complex surface of revolution formed by rotating a line segment on a two-dimensional plane 360 ​​degrees around its central axis of symmetry. Due to the high axisymmetry of the arc-extinguishing chamber structure of a high-voltage circuit breaker, the radial component of the force (perpendicular to the direction of motion) acting on this surface by the high-pressure gas is equal in magnitude and opposite in direction along the 360-degree circumference, macroscopically canceling each other out and only acting as a compression component. Therefore, to calculate the axial resistance that truly affects the output power of the hydraulic operating mechanism, a geometrically reduced-dimensional projection of the area is necessary. Specifically, according to the formula for the area of ​​a circle, the orthographic projection of this force-bearing interface onto a plane parallel to the direction of motion is actually a surface formed by the outer diameter... and inner diameter The system forms a two-dimensional planar annulus. Therefore, the system directly calculates the square of the radial coordinates of the outer edge ( ) and the square of the radial coordinate of the inner edge ( The difference between the two is calculated by multiplying this difference by pi (π), i.e., by using the mathematical formula... The effective annular projected area S of the force-bearing interface in the direction of motion is obtained.

[0040] Preferably, the dynamic load reaction force data of the gas hindering the movement of the solid transmission component is calculated through the following discrete numerical mapping and summation steps: Considering the severe pressure gradient within the arc-extinguishing chamber under real high-voltage electric arc impact, the pressure distribution at a single force interface often exhibits high non-uniformity. To ensure the accuracy of the counterforce quantification, this embodiment abandons the coarse average pressure product method in its algorithm implementation, and instead adopts a high-precision discrete numerical integration based on grid nodes and a vector direction determination rule.

[0041] Specifically, at any time step in the current fluid dynamics simulation Within the system, assume that a total of M force interfaces (j=1,2,...,M) not parallel to the direction of motion are identified. For the j-th force interface, the system discretizes it radially according to the topology of its surface mesh. Each infinitesimal toroidal surface (i=1,2,..., For the i-th infinitesimal annular surface on the j-th force interface, extract the transient scalar pressure at the corresponding grid center position from the already converged fluid field results, and define it as the gas pressure data of the target node. Simultaneously, the radial coordinates of the outer edge of the i-th infinitesimal annular surface are extracted. radial coordinates of the inner edge Based on the principle of area projection dimensionality reduction, the effective projected area of ​​this infinitesimal annular surface in the direction of motion is... The calculation formula is: ; The gas pressure data of the target node With the effective projected area of ​​the micro-element Perform scalar product, and at all forces at the j-th interface. The summation is performed on each infinitesimal element to obtain the value of the j-th force interface. Data on the partial reaction forces at all times Its mathematical expression is: ; Subsequently, to obtain the macroscopic aerodynamic drag experienced by the entire solid-state transmission system at the current moment, a global algebraic summation of the partial reaction force data at all M force-interfaces is necessary. Since inclined surfaces at different locations (e.g., the front and back surfaces of a piston) may generate resistance opposite to the direction of motion or a driving force in the same direction as the motion when subjected to gas compression, a direction determination mechanism must be introduced. The system extracts the surface unit normal vector of the j-th force-interface. (Defining the direction as perpendicular from the fluid domain to the solid domain), and obtaining the current macroscopic motion direction vector of the solid transmission component. Define the direction determination coefficient. : ; In the formula, sgn() is the symbolic function. When When (i.e., the direction of gas pressure is the same as the direction of component movement). This indicates that the force generated by that surface is a driving force (negative resistance); when That is, the direction of gas pressure is opposite to the direction of component movement. This indicates that the force generated by this surface is pure resistance.

[0042] Combined with the aforementioned direction determination coefficient The partial reaction force data at all force interfaces are globally summed to obtain the solid transmission component's... Total dynamic load reaction force at all times Its mathematical expression is: ; The system will simulate all time steps throughout the entire opening and closing simulation cycle. The calculation obtained below The discrete time series data are recombined according to the time series and the final output set of discrete time series data is the dynamic load reaction force data of the gas hindering the movement of the solid transmission component.

[0043] In this embodiment, the geometric feature extraction, discrete projection dimensionality reduction, and numerical summation with direction determination described above achieve the following technical effects: Traditional three-dimensional force simulation typically employs the full-surface infinitesimal vector integration method. This method is not only extremely computationally intensive but also prone to introducing significant numerical calculation errors when dealing with complex meshes due to mesh distortion, which can lead to incomplete cancellation of radial forces. This embodiment innovatively proposes an algebraic topological algorithm based on infinitesimal radial contour coordinate extraction and discrete accumulation of the effective projected area of ​​the infinitesimal annular ring. Furthermore, it combines a normal vector dot product mechanism to rigorously distinguish between driving force and resistance boundaries. This method establishes a rigid physical filtering mechanism at the computer algorithm level. By utilizing the geometric and mathematical rules of area projection, it directly removes all invalid radial mechanical components, reducing the three-dimensional spatial force calculation of complex curved surfaces to a simplified two-dimensional scalar mapping and one-dimensional algebraic summation. This not only improves the computational efficiency of reaction force extraction but also prevents the introduction of error accumulation during spatial vector synthesis. It ensures that the final extracted dynamic load reaction force data can purely reflect the real resistance effect of high-pressure gas on the axial movement of the hydraulic mechanism, providing a high-precision external load input boundary for subsequent multibody dynamics models.

[0044] Step S3: The dynamic load reaction force data is used as an external load condition and imported into a preset transmission system multibody dynamics simulation model. The kinematic state data of the solid transmission component under the condition of overcoming the dynamic load reaction force is calculated through the transmission system multibody dynamics simulation model. In a preferred embodiment, the kinematic state data of the solid transmission component under the condition of overcoming the dynamic load reaction force are calculated using the multibody dynamics simulation model of the transmission system, including: In the preset simulation platform, the connection pair attributes, friction attributes, and contact constraint attributes between the solid transmission components are configured, and a multibody dynamics simulation model of the transmission system is constructed based on the connection pair attributes, friction attributes, and contact constraint attributes between the solid transmission components. Based on the lumped parameter method, a mathematical model of the hydraulic cylinder is constructed that interacts with the multibody dynamics simulation model of the transmission system; The dynamic load reaction force data is input into the multibody dynamics simulation model of the transmission system, and combined with the system oil pressure data provided by the hydraulic cylinder mathematical model for joint solution to obtain the displacement and velocity data of the solid transmission component during the motion process, which are used as the kinematic state data.

[0045] Preferably, the following multibody system topology definition steps are used to configure the connection pair attributes, friction attributes, and contact constraint attributes between the solid transmission components in a preset simulation platform to construct a multibody dynamics simulation model of the transmission system: First, import the three-dimensional solid assembly model of the high-voltage circuit breaker transmission system into a preset mechanical system dynamics simulation platform (such as ADAMS or other multibody dynamics software). For key components in the transmission system (including but not limited to connecting rods, insulating tie rods, pins, moving arc contacts, and piston rods), assign inertial parameters such as mass, center of mass position, and moment of inertia based on their actual material properties, defining them as rigid or flexible body components. Subsequently, based on the actual relative motion relationships between the components, the connection properties of the geometric topology are configured. Specifically, a revolute joint is configured at the connection between the pin and the connecting rod, restricting all translational and rotational degrees of freedom except for single-axis rotation. A prismatic joint is configured for the fit between the piston rod and the hydraulic cylinder, and between the moving arc contact and the guide sleeve, allowing them to only perform linear translation along a preset direction of motion. Furthermore, to realistically reproduce the energy loss of the mechanical mechanism during high-speed motion, friction properties are configured. Based on the theories of Coulomb friction and viscous friction, preset static and dynamic friction coefficients (e.g., considering the parameters of the Stribek friction model under hydraulic oil lubrication conditions) are input into the prismatic and revolute joints to calculate the dynamic frictional resistance during the mechanism's movement. Finally, contact constraint attributes are configured to address potential mechanical limit collisions that may occur at the end of the opening and closing phases. A nonlinear spring-damped contact model based on the penalty function method (such as Hertz contact theory) is employed to define the contact stiffness, force exponent, damping coefficient, and penetration depth between colliding components, accurately simulating the microscopic collision bounce and peak contact force of the transmission components during the buffering phase. Combining the defined connection pairs, friction, and contact constraints, the system automatically generates a system of multibody dynamic differential-algebraic equations (DAEs) described by Lagrange multipliers, completing the construction of the multibody dynamics simulation model of the transmission system.

[0046] Preferably, the following fluid network modeling steps are used to construct a hydraulic cylinder mathematical model that interacts with the multibody dynamics simulation model of the transmission system based on the lumped parameter method: Since the hydraulic operating mechanism is the power source driving the entire high-voltage circuit breaker, the transient changes in its internal oil pressure directly determine the magnitude of the driving force. To avoid the inefficiency of joint solution caused by establishing a complex three-dimensional hydraulic flow field model, this embodiment preferably uses the lumped parameter method to discretize the continuous hydraulic physical system into a fluid network topology composed of pressure nodes and flow channels. Specifically, a hydraulic cylinder mathematical model including the working cylinder, buffer chamber, high-pressure oil circuit, and control main valve is established. For the rodless and rod chambers of the hydraulic cylinder, based on the principles of volume conservation and fluid compressibility, the nodal pressure differential equations are established: ; Where p represents the current system oil pressure data of the chamber. The effective bulk modulus of hydraulic oil. Let x be the instantaneous chamber volume as a function of piston displacement x. and These represent the volumetric flow rates flowing into and out of the chamber, respectively, and A is the effective working area of ​​the hydraulic cylinder. The instantaneous velocity of the piston.

[0047] Preferably, the dynamic load reaction force data is input and combined with the system hydraulic pressure data for joint solution through the following unidirectional load mapping and Newton-Euler iteration steps: Before starting the joint solution, the system extracts the dynamic load reaction force data (i.e., discrete time series) of the gas hindering the movement of the solid transmission component from step S2. The function is fitted to an external resistance function that varies continuously with time through spline interpolation. Within the iterative time step of the joint solution, the system will use the external resistance function. Configured as an independent temporal load decoupled from the fluid dynamics simulation mesh, it is applied unidirectionally to the force-bearing nodes representing the moving arc contact or piston rod in the multibody dynamics simulation model, with the direction strictly opposite to the direction of motion; simultaneously, the hydraulic cylinder mathematical model applies an active driving force to the piston rod node based on the current system oil pressure data. At this point, the multibody dynamics solver, based on the Newton-Euler equations, performs dynamic equilibrium calculations for the entire transmission system. For the core transmission component that translates along the direction of motion, its dynamic equilibrium equations can be simplified as follows: ; Where M is the system's equivalent mass referred to the direction of motion. For the aforementioned frictional resistance, This is the reaction force of the contact spring or the trip spring. This refers to the instantaneous acceleration of the system.

[0048] The solver employs a high-order numerical integration algorithm (such as the Gear method or the implicit backward differential formula BDF) to solve the aforementioned nonlinear differential-algebraic equations through time integration. At the end of each integration step, it outputs the actual spatial position and velocity variables of the solid transmission component after overcoming the resistance of the external high-pressure gas. These time-varying displacement data x(t) and velocity data v(t) represent the desired kinematic state data.

[0049] In this embodiment, the multibody dynamics topology construction, hydraulic lumped parameter modeling, and unidirectional load joint solution processing described above achieve the following technical effects: Traditional high-voltage circuit breaker co-simulation often attempts to use two-way fluid-structure interaction (2-Way FSI) technology, that is, continuously exchanging forces and displacements in each time step. This method is not only prone to fluid calculation divergence and collapse due to severe mesh deformation, but also has an extremely long calculation cycle. This embodiment constructs a unidirectional joint solution architecture based on pre-solved resistance spline mapping, transforming the complex arc-extinguishing chamber fluid field pressure into a clean, independent external time-series resistance curve, and combining it with a lumped parameter hydraulic model and a refined friction and contact rigid body topology for solution. This architecture fundamentally decouples the nonlinear oscillation conflict between fluid mesh deformation and mechanical rigid body displacement. While ensuring nanosecond-level high-precision dynamic balance of hydraulic driving force, gas resistance, spring force, and friction force, it also ensures that the entire multibody dynamics model has extremely high numerical convergence stability. This allows for accurate and efficient output of component displacement and velocity states under real physical conditions, laying a reliable kinematic foundation for subsequent accurate stripping of output power and stress fatigue assessment.

[0050] Step S4: Based on the kinematic state data and the dynamic load reaction force data, calculate the output power data of the hydraulic operating mechanism and the stress distribution law of the solid transmission component, and evaluate the mechanical life decay of the high-voltage circuit breaker according to the stress distribution law of the solid transmission component.

[0051] In a preferred embodiment, based on the kinematic state data and the dynamic load reaction force data, the output work data of the hydraulic operating mechanism and the stress distribution law of the solid transmission component are calculated, including: calculating the kinetic energy data of the solid transmission component during its motion based on the component velocity data and the mass data of the solid transmission component in the kinematic state data; calculating the load work data required for the solid transmission component to overcome gas resistance based on the component displacement data and the corresponding dynamic load reaction force data in the kinematic state data; and combining the kinetic energy data and the load work data to obtain the output work data of the hydraulic operating mechanism driving the solid transmission component. The dynamic load reaction data and the kinematic state data are used as force boundary conditions and input into a pre-constructed transmission system stress analysis model for solution, thereby obtaining the stress distribution law of the solid transmission component.

[0052] Preferably, assessing the mechanical life degradation of the high-voltage circuit breaker based on the stress distribution pattern of the solid-state transmission component includes: extracting stress extreme value data of key connecting components in the solid-state transmission component from the stress distribution pattern, wherein the key connecting components include at least one of connecting rods, insulating tie rods, and pins; and assessing the mechanical life degradation of the high-voltage circuit breaker based on preset material fatigue life parameters and the stress extreme value data.

[0053] Preferably, the output power data of the hydraulic operating mechanism driving the solid transmission component is calculated through the following energy decoupling and numerical integration steps: Since the energy consumption of the circuit breaker during opening and closing is highly dynamic and nonlinear, this embodiment, based on the law of conservation of energy, strictly decomposes the total output power of the hydraulic operating mechanism into two parts in a physical mechanism: the first part is the kinetic energy necessary to maintain the accelerated movement of each transmission component of the mechanism, and the second part is the work consumed to overcome the resistance of high-pressure gas. Regarding the kinetic energy data (denoted as...),... The system extracts the instantaneous velocity component v(t) of the solid transmission component at the current time t from the kinematic state data and obtains the preset equivalent mass m of the component. Based on the classical kinetic energy theorem, the formula... The kinetic energy sequence at each time point throughout the entire opening and closing stroke is continuously calculated. The work done by the load (denoted as...) is then analyzed. The system extracts the instantaneous displacement data x(t) of the solid transmission component along the direction of motion, and the dynamic load reaction force data extracted in step S2. Since both force and displacement are discrete sequences that change with time, the system employs numerical integration algorithms (such as the trapezoidal rule or Simpson's rule) to calculate and accumulate the work done within the infinitesimal displacement segment. Its discrete integral mathematical expression is: ; Where n represents the total number of discrete time steps from the start of the component's movement to the current time t, and k is the cyclic index variable of the time steps. This refers to the specific time corresponding to the k-th time step. The displacement travel within the k-th time step precisely quantifies the mechanical energy consumed by the transmission system to overcome the pressure explosion effect of the airflow field in the arc-extinguishing chamber. Finally, the kinetic energy data at the same moment... With the load work data Perform algebraic addition and combination (i.e.) This outputs a complete time-energy curve, which represents the output power data of the hydraulic operating mechanism.

[0054] Preferably, the dynamic load reaction data and the kinematic state data are input into a pre-built transmission system stress analysis model for solution through the following finite element boundary mapping and transient structure solution steps: The pre-built transmission system stress analysis model refers to a three-dimensional structural mechanics mathematical model formed by discretizing the solid geometric model of the transmission system with a high-quality mesh in professional finite element analysis (FEA) software and then assigning constitutive relations (such as Young's modulus, Poisson's ratio, yield strength, etc.) to linear elastic or elastoplastic materials. During the solution process, the system performs rigorous boundary condition mapping: on the one hand, it extracts the instantaneous acceleration time series from the kinematic state data and, based on d'Alembert's principle, transforms the instantaneous acceleration time series into inertial force loads, which are then applied to the centroid nodes of each component in the stress analysis model to realistically reproduce the structural inertial impact caused by high-speed motion; on the other hand, it treats the dynamic load reaction data as external concentrated forces or equivalent surface loads, accurately mapping them onto the force-bearing end faces of the piston rod or moving arc contact; simultaneously, it applies nonlinear contact constraint boundaries (such as face-to-face contact pairs with friction coefficients) at the relative hinges of the connecting rod, insulating tie rod, and pin. After the boundary conditions are applied, the finite element solver uses an implicit time integration algorithm (such as the Newmark-β method) to iteratively solve the transient dynamic equations of the system, which include the mass matrix, damping matrix, and stiffness matrix. After the solution converges, the system outputs the stress tensor of each grid node inside the transmission component during the entire opening and closing cycle, and calculates the equivalent stress spatiotemporal distribution cloud map under the fourth strength theory. The cloud map and the set of stress time series of the underlying nodes contained therein are the stress distribution law.

[0055] Preferably, the mechanical life decay of the high-voltage circuit breaker is evaluated based on preset material fatigue life parameters and the stress extreme value data through the following extreme value extraction and fatigue damage accumulation steps: First, the system activates a stress hotspot monitoring mechanism to scan and lock stress concentration areas (such as the shear section of the pin and the variable cross-section fillet of the connecting rod) on key connecting components such as connecting rods, insulating tie rods, and pins from the massive stress distribution data. Using the rainflow counting method, the complex alternating stress time series of these hotspot areas is extracted cyclically to isolate the effective stress cycles causing fatigue damage, thereby accurately obtaining the maximum stress amplitude and average stress of the key connecting components in a single opening and closing operation, which are used as the stress extreme value data. Subsequently, the system calls the preset material fatigue life parameters, which refer to the SN curve (stress-life curve) data matrix pre-determined through standard material fatigue tests and stored in the database, as well as empirical correction factors used to correct surface roughness, size effect, and stress concentration coefficient. During the evaluation phase, the system uses the Goodman mean stress correction equation to perform an equivalent transformation on the extracted stress extremum data. Then, the transformed stress extrema are substituted into the corresponding SN curve model, and the limit number of cycles that the component can withstand at that stress level is calculated by looking up a table or interpolation. .

[0056] Finally, the Miner linear fatigue damage accumulation rule is introduced ( ),in, This represents the number of cycles actually experienced by the component at the i-th equivalent alternating stress amplitude level, obtained through the rainflow counting method within the analyzed opening and closing operation cycle. This number is used to calculate the local fatigue damage degree caused to the critical component by a single high-capacity short-circuit breaking operation. By comparing this damage degree with a warning threshold, or converting it into the remaining operable cycles, a quantitative assessment report of the high-voltage circuit breaker's mechanical life degradation is finally output.

[0057] In this embodiment, the above-mentioned energy decomposition, transient finite element mapping, and fatigue damage accumulation assessment process achieves the following technical effects: Traditional high-voltage circuit breaker life assessment often relies on expensive and time-consuming destructive testing of physical prototypes, or simply uses macroscopic motion speed to roughly infer the mechanism's state, which is a black-box empirical assessment that cannot perceive the microscopic damage of internal parts. This invention precisely decomposes hydraulic output power into kinetic energy and work done to overcome air resistance using physical formulas, providing direct numerical indicators for the refined design of hydraulic cylinder volume and operating pressure. On the other hand, by using high-fidelity finite element transient mapping and rainflow counting method, the abstract aerodynamic resistance and inertial impact are perfectly transformed into intuitive stress extremes and fatigue damage degrees of weak links such as pins and connecting rods. This white-box quantitative assessment architecture based on full-link physical field coupling not only eliminates the excessive reliance on physical experiments, but also achieves accurate prediction of hidden mechanical fatigue of transmission systems under complex load conditions, providing reliable theoretical and data support for preventive maintenance, lightweight structural design, and full life-cycle reliability management of high-voltage switchgear.

[0058] like Figure 2 As shown, another embodiment of the present invention also provides a life assessment device for high-voltage circuit breakers based on field-circuit coupling, including: a fluid field simulation module, a load reaction force mapping module, a multibody dynamics solution module, and a comprehensive life assessment module; The fluid field simulation module is used to acquire the geometric model data of the high-voltage circuit breaker and perform fluid dynamics simulation on the geometric model data to calculate the gas pressure distribution data of the gas inside the arc extinguishing chamber acting on the solid transmission component during the opening and closing process of the high-voltage circuit breaker. The load reaction force mapping module is used to identify the force interface on the surface of the solid transmission component that is not parallel to the direction of motion of the solid transmission component, extract the radial contour coordinate data of the force interface, and calculate the annular effective projection area of ​​the force interface in the direction of motion based on the radial contour coordinate data. Then, the gas pressure distribution data is mapped onto the annular effective projection area to calculate the dynamic load reaction force data of the gas that hinders the movement of the solid transmission component. The multibody dynamics solution module is used to import the dynamic load reaction force data as an external load condition into a preset transmission system multibody dynamics simulation model, and calculate the kinematic state data of the solid transmission component under the condition of overcoming the dynamic load reaction force through the transmission system multibody dynamics simulation model. The lifespan comprehensive assessment module is used to calculate the output power data of the hydraulic operating mechanism and the stress distribution law of the solid transmission component based on the kinematic state data and the dynamic load reaction force data, and to assess the mechanical lifespan degradation of the high-voltage circuit breaker based on the stress distribution law of the solid transmission component.

[0059] It is understood that the above-described device embodiments correspond to the method embodiments of the present invention, and can realize the life assessment method for high-voltage circuit breakers based on field-circuit coupling provided by any of the above-described method embodiments of the present invention.

[0060] It should be noted that the device embodiments described above are merely illustrative, and some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Furthermore, in the accompanying drawings of the device embodiments provided by this invention, the connection relationships between modules indicate that they have communication connections, which can specifically be implemented as one or more communication buses or signal lines. Those skilled in the art can understand and implement this without any creative effort.

[0061] Based on the above embodiment of the life assessment method for high-voltage circuit breakers based on field-circuit coupling, another embodiment of the present invention provides a terminal device, which includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor. When the processor executes the computer program, it implements any one of the life assessment methods for high-voltage circuit breakers based on field-circuit coupling of the present invention.

[0062] For example, in this embodiment, the computer program can be divided into one or more modules, which are stored in the memory and executed by the processor to complete the present invention. The one or more modules may be a series of computer program instruction segments capable of performing a specific function, which describe the execution process of the computer program in the terminal device.

[0063] The terminal device may be a desktop computer, laptop, handheld computer, or cloud server, etc. The terminal device may include, but is not limited to, a processor and a memory.

[0064] The processor can be a Central Processing Unit (CPU), or 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 can be a microprocessor or any conventional processor. The processor is the control center of the terminal device, connecting all parts of the terminal device via various interfaces and lines.

[0065] Based on the above-described method embodiments, another embodiment of the present invention provides a computer-readable storage medium including a stored computer program, wherein, when the computer program is executed, it controls the device where the computer-readable storage medium is located to execute the life assessment method for a high-voltage circuit breaker based on field-circuit coupling as described in any of the above-described method embodiments of the present invention.

[0066] The modules / units integrated in the device / terminal equipment, if implemented as software functional units and sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, all or part of the processes in the above embodiments of the present invention can also be implemented by a computer program instructing related hardware. The computer program can be stored in a computer-readable storage medium, and when executed by a processor, it can implement the steps of the various method embodiments described above. The computer program includes computer program code, which can be in the form of source code, object code, executable files, or certain intermediate forms. The computer-readable medium can include: any entity or device capable of carrying the computer program code, recording media, USB flash drives, portable hard drives, magnetic disks, optical disks, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signals, telecommunication signals, and software distribution media, etc.

[0067] The above description represents the preferred embodiments of the present invention. It should be noted that those skilled in the art can make various improvements and modifications without departing from the principles of the present invention, and these improvements and modifications are also considered to be within the scope of protection of the present invention.

Claims

1. A life assessment method for high-voltage circuit breakers based on field-circuit coupling, characterized in that, include: Geometric model data of a high-voltage circuit breaker is acquired, and fluid dynamics simulation is performed on the geometric model data to calculate the gas pressure distribution data of the gas inside the arc-extinguishing chamber acting on the solid-state transmission component during the opening and closing process of the high-voltage circuit breaker. The force interface on the surface of the solid-state transmission component that is not parallel to the direction of movement of the solid-state transmission component is identified, and the radial contour coordinate data of the force interface is extracted. Based on the radial contour coordinate data, the annular effective projected area of ​​the force interface in the direction of movement is calculated. Then, the gas pressure distribution data is mapped onto the annular effective projected area to calculate the dynamic load reaction force data of the gas hindering the movement of the solid-state transmission component. The dynamic load reaction force data is used as an external load condition and imported into a preset multibody dynamics simulation model of the transmission system. The kinematic state data of the solid transmission component under the condition of overcoming the dynamic load reaction force is calculated through the multibody dynamics simulation model of the transmission system. Based on the kinematic state data and the dynamic load reaction force data, the output power data of the hydraulic operating mechanism and the stress distribution law of the solid transmission component are calculated respectively. Based on the stress distribution law of the solid transmission component, the mechanical life decay of the high-voltage circuit breaker is evaluated.

2. The life assessment method for high-voltage circuit breakers based on field-circuit coupling as described in claim 1, characterized in that, A fluid dynamics simulation is performed on the geometric model data to calculate the gas pressure distribution data of the gas inside the arc-extinguishing chamber acting on the solid-state transmission component during the opening and closing process of the high-voltage circuit breaker. This includes: constructing a two-dimensional axisymmetric magnetohydrodynamic model based on the geometric model data; acquiring the short-circuit current data of the high-voltage circuit breaker during the breaking process, and calculating the Joule heating and radiative energy loss of the arc plasma based on the short-circuit current data; substituting the Joule heating and radiative energy loss as energy source terms into the energy conservation equation of the two-dimensional axisymmetric magnetohydrodynamic model; solving the energy conservation equation to obtain the gas pressure distribution data of the gas inside the arc-extinguishing chamber acting on the solid-state transmission component.

3. The life assessment method for high-voltage circuit breakers based on field-circuit coupling as described in claim 2, characterized in that, Extracting the radial contour coordinate data of the force-bearing interface, and calculating the annular effective projected area of ​​the force-bearing interface in the direction of motion based on the radial contour coordinate data, includes: extracting the radial coordinates of the outer edge and the inner edge of the force-bearing interface in the radial direction; calculating the difference between the square of the outer edge radial coordinate and the square of the inner edge radial coordinate; and multiplying the difference by pi as the annular effective projected area of ​​the force-bearing interface in the direction of motion.

4. The life assessment method for high-voltage circuit breakers based on field-circuit coupling as described in claim 3, characterized in that, Mapping the gas pressure distribution data onto the annular effective projected area and calculating the dynamic load reaction force data of the gas hindering the movement of the solid transmission component includes: extracting the target node gas pressure data corresponding to the force interface from the gas pressure distribution data; multiplying the target node gas pressure data by the annular effective projected area as the partial reaction force data borne by the force interface in the direction of movement; and summing the partial reaction force data of all the force interfaces to obtain the dynamic load reaction force data of the gas hindering the movement of the solid transmission component.

5. The life assessment method for high-voltage circuit breakers based on field-circuit coupling as described in claim 4, characterized in that, The kinematic state data of the solid transmission component under the condition of overcoming the dynamic load reaction force are calculated using the multibody dynamics simulation model of the transmission system, including: In the preset simulation platform, the connection pair attributes, friction attributes, and contact constraint attributes between the solid transmission components are configured, and a multibody dynamics simulation model of the transmission system is constructed based on the connection pair attributes, friction attributes, and contact constraint attributes between the solid transmission components. Based on the lumped parameter method, a mathematical model of the hydraulic cylinder is constructed that interacts with the multibody dynamics simulation model of the transmission system; The dynamic load reaction force data is input into the multibody dynamics simulation model of the transmission system, and combined with the system oil pressure data provided by the hydraulic cylinder mathematical model for joint solution to obtain the displacement and velocity data of the solid transmission component during the motion process, which are used as the kinematic state data.

6. The life assessment method for high-voltage circuit breakers based on field-circuit coupling as described in claim 5, characterized in that, Based on the kinematic state data and the dynamic load reaction force data, the output work data of the hydraulic operating mechanism and the stress distribution law of the solid transmission component are calculated, including: calculating the kinetic energy data of the solid transmission component during its motion based on the component velocity data and the mass data of the solid transmission component in the kinematic state data; calculating the load work data required for the solid transmission component to overcome gas resistance based on the component displacement data and the corresponding dynamic load reaction force data in the kinematic state data; and combining the kinetic energy data and the load work data to obtain the output work data of the hydraulic operating mechanism driving the solid transmission component. The dynamic load reaction data and the kinematic state data are used as force boundary conditions and input into a pre-constructed transmission system stress analysis model for solution, thereby obtaining the stress distribution law of the solid transmission component.

7. The life assessment method for high-voltage circuit breakers based on field-circuit coupling as described in claim 6, characterized in that, The mechanical life degradation of the high-voltage circuit breaker is evaluated based on the stress distribution pattern of the solid-state transmission component, including: extracting stress extreme value data of key connecting components in the solid-state transmission component from the stress distribution pattern, wherein the key connecting components include at least one of connecting rod, insulating tie rod and pin; and evaluating the mechanical life degradation of the high-voltage circuit breaker based on preset material fatigue life parameters and the stress extreme value data.

8. A life assessment device for high-voltage circuit breakers based on field-circuit coupling, characterized in that, include: The system includes a fluid field simulation module, a load reaction force mapping module, a multibody dynamics solution module, and a comprehensive life assessment module. The fluid field simulation module is used to acquire the geometric model data of the high-voltage circuit breaker and perform fluid dynamics simulation on the geometric model data to calculate the gas pressure distribution data of the gas inside the arc extinguishing chamber acting on the solid transmission component during the opening and closing process of the high-voltage circuit breaker. The load reaction force mapping module is used to identify the force interface on the surface of the solid transmission component that is not parallel to the direction of motion of the solid transmission component, extract the radial contour coordinate data of the force interface, and calculate the annular effective projection area of ​​the force interface in the direction of motion based on the radial contour coordinate data. Then, the gas pressure distribution data is mapped onto the annular effective projection area to calculate the dynamic load reaction force data of the gas that hinders the movement of the solid transmission component. The multibody dynamics solution module is used to import the dynamic load reaction force data as an external load condition into a preset transmission system multibody dynamics simulation model, and calculate the kinematic state data of the solid transmission component under the condition of overcoming the dynamic load reaction force through the transmission system multibody dynamics simulation model. The lifespan comprehensive assessment module is used to calculate the output power data of the hydraulic operating mechanism and the stress distribution law of the solid transmission component based on the kinematic state data and the dynamic load reaction force data, and to assess the mechanical lifespan degradation of the high-voltage circuit breaker based on the stress distribution law of the solid transmission component.

9. A terminal device, characterized in that, The device includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, wherein when the processor executes the computer program, it implements the life assessment method for a high-voltage circuit breaker based on field-circuit coupling as described in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, include: A stored computer program, wherein, when the computer program is executed, it controls the device containing the computer-readable storage medium to perform the life assessment method for a high-voltage circuit breaker based on field-circuit coupling as described in any one of claims 1-7.

Images