Analysis Method for Electromagnetic Fluid Simulation Modeling of AC Contactor

Abstract

The present invention relates to the technical field of simulation modeling analysis, specifically an electromagnetic fluid simulation modeling analysis method applied to AC contactors. In the present invention, the displacement data of the armature is collected in real time through a displacement sensor, and a dynamic mapping relationship table containing a hysteresis compensation coefficient is constructed; based on the displacement data and the compensation coefficient, the air-gap magnetic permeability is dynamically corrected to generate a dynamic electromagnetic field coupling boundary condition, realizing the real-time matching of magnetic parameters and mechanical motion; in the fluid field simulation, a deformable sub-grid nesting technology is adopted, and the elastic deformation control of the arc region grid is implemented in combination with the displacement change rate, and the influence of the magnetic field distortion on the grid is synchronously corrected by associating the hysteresis compensation coefficient, solving the coupling problem of the dynamic deformation of the arc and the hysteresis non-linearity; at the same time, the main grid topology preservation module is used to fix the non-arc region, and the deformation-sensitive region is screened to perform local grid reorganization, realizing high-precision dynamic simulation of the severely deformed arc region on the premise of ensuring calculation stability.

Description

Technical Field

[0001] The present invention relates to the technical field of simulation modeling analysis, and specifically to an electromagnetic fluid simulation modeling analysis method applied to AC contactors. Background Art

[0002] As a core device for realizing circuit on-off control in the power system, the dynamic characteristics of an AC contactor directly affect the reliability of the device and the arc suppression performance. Traditional electromagnetic fluid simulation modeling methods usually adopt static mesh division and fixed magnetic permeability parameters, and it is difficult to accurately describe the influence of the dynamic change of the air gap on the magnetic field distribution during the high-speed movement of the armature. Especially in the transient process of closing and opening, due to the strong non-linear relationship between the armature displacement and the air gap magnetic permeability, the existing technologies rely on empirical formulas or off-line calibration data, resulting in the accumulation of calculation deviation of the magnetic field strength with the simulation time, insufficient compensation for the hysteresis effect, and the inability to accurately reflect the residual magnetism and eddy current loss under alternating current.

[0003] In addition, when the arc plasma is coupled with the flow field, the traditional fixed mesh generates local distortion due to the drastic change of the arc shape, and frequent global mesh reconstruction is required, resulting in low calculation efficiency and easy numerical oscillation; although the existing dynamic mesh technology can be adjusted locally, it lacks the coordinated control with the armature movement state and hysteresis compensation, and it is difficult to balance accuracy and stability. At the data acquisition level, conventional displacement sensors are limited by the physical installation position and fixed sampling frequency, and it is difficult to capture the instantaneous displacement details during the high-speed movement of the armature, resulting in the distortion of the displacement-magnetic permeability mapping relationship;

[0004] The above problems make the existing simulation models have significant errors in predicting the dynamic characteristics of contactors, restricting the optimal design and life assessment of high-performance contactors. Summary of the Invention

[0005] The purpose of the present invention is to provide an electromagnetic fluid simulation modeling analysis method applied to AC contactors to solve the problems raised in the above background art. The core problems to be solved include how to achieve high-precision acquisition of armature displacement data and dynamic air gap magnetic permeability correction to solve the problems of insufficient hysteresis effect compensation, poor mesh adaptability, and low calculation efficiency in the electromagnetic-fluid field coupling simulation; and how to solve the problem of simulation instability caused by grid distortion in the arc region through dynamic mesh reorganization and elastic deformation control.

[0006] To achieve the above purpose, the present invention provides the following technical solutions: An electromagnetic fluid simulation modeling analysis method applied to AC contactors, and the method steps include:

[0007] S1. By setting virtual nodes of displacement sensors at the closed critical points, maximum opening displacement points, and turning points of the armature movement trajectory, configuring independent data acquisition channels and overlapping detection areas, measurement errors are eliminated to ensure the accuracy of displacement data acquisition. Further, the sliding average filtering and cubic spline interpolation algorithms are used to preprocess the original armature displacement data with non-uniform sampling to generate a displacement data sequence with uniform intervals, solving the problem of insufficient capture of high-speed movement details;

[0008] The hysteresis compensation coefficient is generated by combining the air-gap permeability differences at the same displacement points during opening and closing processes and is dynamically updated through direction marking. The amplitude is corrected by superimposing the oscillation suppression factor at the closed critical point to achieve real-time compensation of the hysteresis effect. A dynamic mapping relationship table including the reference air-gap permeability, hysteresis compensation coefficient, and motion state is constructed to provide an accurate data basis for dynamic correction of the air-gap permeability.

[0009] S2. Based on the reference air-gap permeability and hysteresis compensation coefficient in the dynamic mapping relationship table, the permeability is corrected in real time. The air-gap region is divided into high-density grids through adaptive grid meshing technology. Combining the finite-difference time-domain method to solve the three-dimensional transient Maxwell's equations, a dynamic coupling boundary condition including the magnetic field strength gradient distribution, electromagnetic force density distribution, and dynamically updated region coordinate set is generated. The boundary element set to be refreshed in real time is screened through the magnetic field gradient threshold to drive the local grid deformation calculation of the fluid field simulation module, solving the problem that traditional fixed grids cannot adapt to dynamic air-gap changes;

[0010] S3. In the fluid field simulation module, according to the magnetic field strength gradient distribution and electromagnetic force density distribution of the dynamic coupling boundary condition, the arc region grid is separated from the main grid to generate a deformed sub-grid. The non-matching interface interpolation algorithm is used to establish a boundary data transfer channel between the deformed sub-grid and the main grid to ensure seamless coupling of electromagnetic-fluid field parameters;

[0011] The deformation sensitivity is dynamically adjusted in combination with the change rate of the armature displacement data. When the change rate exceeds the set change threshold, the deformation margin amplification mode is activated to increase the maximum allowable deformation of the nodes. Based on the change trend, the arc expansion direction is predicted to adjust the deformation parameters, and the deformation freedom of the nodes is optimized by synchronously associating the hysteresis compensation coefficient, solving the problem of simulation instability caused by grid distortion in the arc region.

[0012] S4. Extract the gradient value of the reference air-gap permeability and the change rate of the armature displacement data through the dynamic mapping relationship table. When both exceed the corresponding thresholds, mark them as deformation-sensitive regions. Based on the distortion rate data of the deformable sub-grids, use the local grid recombination algorithm to re-divide the sub-grid units whose distortion rates exceed the grid distortion rate threshold, synchronously update the node coordinates and the boundary connection relationships, and re-couple with the main grid through the non-matching interface interpolation algorithm. On the premise of keeping the topology structure of the main grid fixed, solve the problem of the decrease in simulation accuracy caused by the over-limit of the global grid distortion rate.

[0013] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0014] Through the multi-position deployment and data consistency verification of the virtual nodes of the displacement sensor, combined with the moving average filtering and cubic spline interpolation algorithms, eliminate the displacement acquisition noise during the high-speed movement of the armature, and ensure the data reliability of the reference air-gap permeability and the hysteresis compensation coefficient in the dynamic mapping relationship table. Based on the reference air-gap permeability driven by real-time correction, drive the adaptive grid meshing of the electromagnetic field simulation module, and accurately transfer the magnetic field intensity gradient and the electromagnetic force density distribution to the fluid field through the dynamic coupling boundary conditions, reducing the magnetic field-fluid field coupling phase deviation caused by the lagging update of the air-gap permeability in the traditional method.

[0015] Adopt the deformable sub-grid nesting technology to dynamically adjust the deformation sensitivity in combination with the displacement change rate, synchronously correlate the hysteresis compensation coefficient to compensate for the influence of the residual magnetism of the material, suppress the local distortion caused by the contact bounce, and maintain the convergence stability of the arc plasma temperature gradient and the fluid field continuity equation. Through the main grid topology preservation and the threshold screening mechanism of the deformation-sensitive region, on the premise of keeping the non-arc region grid fixed, quickly execute local recombination calculations for the sub-grid units with over-limit distortion, avoid the consumption of computing resources for global grid reconstruction, and improve the iteration efficiency and dynamic condition adaptability of the electromagnetic-fluid two-way coupling simulation. Description of the Drawings

[0016] Figure 1 It is a schematic diagram of the method steps of the present invention. Detailed Embodiments

[0017] Next, the technical solutions in the embodiments of the present invention will be clearly and completely described in conjunction with the drawings in the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments in the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

[0018] Please refer to Figure 1 , the present invention provides a technical solution: an electromagnetic-fluid simulation modeling and analysis method applied to an AC contactor, including the following method steps:

[0019] S1. Set virtual nodes at key positions of the armature movement trajectory, dynamically adjust the sampling frequency according to the armature acceleration to capture the details of displacement changes, and preprocess the real-time displacement data through moving average filtering and cubic spline interpolation algorithms;

[0020] Set displacement sensor virtual nodes at key positions of the armature movement trajectory, including the closing critical point, the maximum opening displacement point, and the movement direction turning point; each displacement sensor virtual node is configured with an independent displacement data acquisition channel, and an overlapping detection area is set at the junction of the opening acceleration section and the deceleration section. By comparing the displacement data consistency of adjacent displacement sensor virtual nodes, measurement errors are eliminated, and the original armature displacement data is generated, ensuring that the acquisition accuracy error of the original armature displacement data is less than 0.05 mm; perform moving average filtering on the original armature displacement data collected by the displacement sensor virtual nodes, and the width of the filtering window is dynamically adjusted according to the displacement data sampling frequency; use the cubic spline interpolation algorithm to resample the non-uniformly sampled displacement data sequence at equal intervals to generate a displacement data sequence with uniform time intervals as the armature displacement data.

[0021] Dynamically adjust the sampling frequency of the armature displacement data according to the acceleration of the armature displacement data. When the absolute value of the acceleration exceeds the set acceleration threshold, the sampling frequency of the armature displacement data is increased to 3 times the reference value; restore the reference displacement data sampling frequency during the uniform motion stage, and the reference displacement data sampling frequency is set to 1 kHz;

[0022] Calculate the hysteresis compensation coefficient based on the air-gap permeability difference at the same displacement point in historical opening and closing operations. By extracting the air-gap permeability difference at the same displacement point during opening and closing processes, combined with the direction mark (opening / closing) of the armature displacement data, generate a dynamically updated hysteresis compensation coefficient; superimpose an oscillation suppression factor at the closing critical point to correct the amplitude of the hysteresis compensation coefficient;

[0023] Fix the armature position through a high-precision displacement control platform, traverse the full stroke at intervals of 0.1 mm, and combine the measured magnetic induction intensity by a Hall effect magnetic flux detector to generate initial reference data on the relationship between armature displacement data and air-gap permeability as the reference air-gap permeability, and store it in the dynamic mapping relationship table;

[0024] Construct a four-dimensional dynamic mapping relationship table including armature displacement data, reference air-gap permeability, hysteresis compensation coefficient, and motion state; after each opening and closing operation, update the dynamic mapping relationship table with the new armature displacement data and hysteresis compensation coefficient according to a 7:3 weight ratio.

[0025] In step S1, aiming at the problems of large displacement data acquisition error and hysteresis compensation lag in traditional methods, in this step, virtual displacement sensor nodes are set at the closed critical point, the maximum opening displacement point, and the turning point of the moving direction of the armature. Independent data acquisition channels and overlapping detection areas are configured, and measurement errors are eliminated through the consistency verification of adjacent node data. The sliding average filter (the window width is dynamically adjusted according to the sampling frequency) and the cubic spline interpolation algorithm are used to resample the non-uniformly sampled data at equal intervals to generate armature displacement data with an error less than 0.05 mm. The hysteresis compensation coefficient is generated based on the air-gap permeability difference at the same displacement point in historical opening and closing operations, dynamically updated in combination with the direction mark, and an oscillation suppression factor is superimposed at the closed critical point to correct the amplitude. Initial reference data is established through a high-precision displacement platform and Hall magnetic flux detection, and a four-dimensional dynamic mapping relationship table including displacement data, reference air-gap permeability, hysteresis compensation coefficient, and motion state is constructed. The data is updated according to the weight ratio to solve the problems of displacement-permeability mapping distortion and insufficient hysteresis effect compensation.

[0026] S2. Query the dynamic mapping relationship table through the armature displacement data to obtain the reference air-gap permeability at the current displacement point, and superimpose the hysteresis compensation coefficient to correct the reference air-gap permeability in real time for dynamically updating the reference air-gap permeability. The correction formula is as follows:

[0027] , where represents the corrected reference air-gap permeability; represents the reference air-gap permeability; represents the hysteresis compensation coefficient;

[0028] The corrected reference air-gap permeability is input into the finite element solver. The adaptive mesh refinement technology is used to perform high-density mesh division on the air-gap region, and the three-dimensional transient Maxwell's equations are solved by the finite-difference time-domain method to calculate the magnetic field strength gradient distribution. Based on the Maxwell stress tensor formula, the spatial distribution of the electromagnetic force density of the magnetic field-arc force field is derived. The temperature gradient and ion concentration distribution data of the arc plasma are synchronously extracted, and the flow field energy exchange threshold condition is defined in combination with the flow field continuity equation. Finally, a dynamic coupling boundary condition including the magnetic field strength gradient, the electromagnetic force density distribution, and the dynamically updated region coordinate set is generated. The dynamically updated region coordinate set is screened by the magnetic field gradient threshold and marked as the boundary element set that needs to be refreshed in real time, which is used to drive the local mesh deformation calculation of the fluid field simulation module to ensure the synchronization and calculation efficiency of the electromagnetic-fluid two-way coupling.

[0029] At the same time, the accuracy of the boundary condition is verified by comparing historical simulation data. When the detected magnetic field strength deviation exceeds 5%, the armature displacement data resampling verification is triggered and the weight ratio of the reference air-gap permeability in the dynamic mapping relationship table is corrected in the reverse direction to form a closed-loop feedback mechanism.

[0030] In step S2, aiming at the problem of cumulative deviation in magnetic field calculation caused by the update delay of air-gap magnetic permeability, the reference air-gap magnetic permeability is corrected in real time through a dynamic mapping relationship table. After being input into the finite element solver, the adaptive mesh refinement technology is adopted to densely divide the air-gap region, and the three-dimensional transient Maxwell's equations are solved by the finite-difference time-domain method to calculate the magnetic field intensity gradient distribution. Based on the Maxwell stress tensor, the electromagnetic force density distribution is derived, and the energy exchange threshold is defined in combination with the fluid flow continuity equation to generate a dynamic coupling boundary condition that includes the magnetic field gradient, electromagnetic force density, and dynamically updated region coordinate set. The real-time refresh region is screened through the magnetic field gradient threshold, and the displacement resampling verification and reverse correction of magnetic permeability weight are triggered based on the over-limit deviation of the magnetic field intensity, forming a closed-loop feedback mechanism to solve the problem of cumulative magnetic field calculation error caused by traditional empirical formulas.

[0031] S3. Based on the dynamic coupling boundary condition in step S2, the deformable sub-grid nesting technology is used to generate deformed sub-grids in the fluid field simulation module. Among them, according to the magnetic field intensity gradient and electromagnetic force density distribution in the dynamic coupling boundary condition, the arc region grid is separated from the main grid, and based on the magnetic field intensity gradient distribution and electromagnetic force density distribution in the dynamic coupling boundary condition, the high-gradient plasma region within the arc region grid is identified and marked as the deformed sub-grid; the non-matching interface interpolation algorithm is used to establish the boundary data transfer channel between the deformed sub-grid and the main grid to ensure the seamless coupling of electromagnetic-fluid field parameters; the coverage range of the deformed sub-grid is dynamically expanded or contracted according to the magnetic field gradient threshold. For example, when the magnetic field gradient exceeds 50 T / m, the adjacent high-gradient cells are automatically included in the deformed sub-grid to form a deformed sub-grid that changes in real time with the magnetic field.

[0032] Combined with the change rate of the armature displacement data in step S1, the deformation amplitude of the deformed sub-grid is dynamically controlled. According to the real-time amplitude of the change rate of the armature displacement data, the deformation sensitivity of the nodes in the deformed sub-grid is adjusted. When the change rate of the armature displacement data exceeds the set change threshold, the deformation margin amplification mode is activated to increase the maximum allowable deformation amount of the nodes in the deformed sub-grid. For example, the maximum allowable deformation amount of the nodes in the deformed sub-grid is increased to twice the reference value; at the same time, based on the change trend of the change rate of the armature displacement data, the expansion direction of the arc shape is predicted, and the deformation direction and rate of the deformed sub-grid are adjusted in advance to ensure that the deformed sub-grid is strictly synchronized with the armature motion state.

[0033] By correlating the hysteresis compensation coefficient in step S1, the generation accuracy of the deformed sub-grid is optimized. According to the magnitude of the hysteresis compensation coefficient, the deformation degrees of freedom of the nodes in the deformed sub-grid are proportionally enlarged to compensate for the arc shape hysteresis caused by the remanence of the material. Near the closing critical point, an oscillation suppression mechanism is superimposed. By reducing the deformation sensitivity, the distortion of the sub-grid caused by contact bounce is suppressed, so that the generated deformed sub-grid maintains a stable shape under dynamic working conditions, and the maximum distortion rate is controlled within 15%.

[0034] In step S3, aiming at the problem of simulation instability caused by arc grid distortion, based on the dynamic coupling boundary conditions, the deformable sub-grid nesting technology is used to separate the arc region from the main grid, and a parameter transfer channel between the main and sub-grids is established through the non-matching interface interpolation algorithm. Combining the change rate of the armature displacement data, the deformation sensitivity is dynamically adjusted. When the change rate exceeds the threshold, the deformation margin amplification mode is activated (such as the maximum deformation amount is increased to 2 times the benchmark), and the deformation parameters are adjusted by predicting the arc expansion direction. Synchronously correlate the hysteresis compensation coefficient to optimize the node deformation degrees of freedom and compensate for the arc shape hysteresis caused by the material remanence. Near the closing critical point, an oscillation suppression mechanism is superimposed to reduce the deformation sensitivity and suppress the distortion caused by contact bounce, and the maximum distortion rate of the sub-grid is controlled within 15%, solving the problem of asynchronous deformation in the existing dynamic grid technology.

[0035] S4. Fix the topological structure of the non-arc region grid through the main grid topology preservation module to ensure that the node coordinates and element connection relationships of the non-arc region grid remain unchanged during the simulation process. Based on the reference air-gap permeability gradient and the change rate of the armature displacement data recorded in the dynamic mapping relationship table in step S1, the deformation-sensitive regions are screened. The screening process of the deformation-sensitive regions is as follows:

[0036] Extract the gradient values of the reference air-gap permeability and the change rate of the armature displacement data at each displacement point from the dynamic mapping relationship table, set the permeability gradient threshold and the displacement change rate threshold. When the permeability gradient or displacement change rate of a certain region exceeds the corresponding threshold, it is marked as a deformation-sensitive region; and the distortion rate data of the corresponding region in the deformed sub-grid generated in step S3 is extracted.

[0037] For the screened deformation-sensitive regions, the local grid reorganization algorithm is used to re-divide the deformed sub-grid whose distortion rate exceeds the grid distortion rate threshold, synchronously update the node coordinates and boundary connection relationships of the deformed sub-grid, and at the same time maintain the topological consistency between the main grid and the non-arc region grid. The reorganized deformed sub-grid is recoupled with the main grid through the non-matching interface interpolation algorithm to ensure the seamless transfer of electromagnetic-fluid field parameters. Finally, an updated simulation model that meets the grid quality constraints is output, the local reorganization calculation time is controlled within 20 ms / frame, the global grid distortion rate is stable below 5%, and the change characteristics of the permeability and displacement data in the reorganization region are recorded in real time through the dynamic mapping relationship table to form a closed-loop optimization link.

[0038] In step S4, aiming at the problems of low global reconstruction efficiency and local distortion diffusion, the main grid topology preservation module is used to fix the grid structure of the non-arc region. Based on the dynamic mapping relation table, the deformation-sensitive regions with excessive permeability gradient and displacement change rate are screened out. For the deformed sub-grids with distortion rate exceeding the grid distortion rate threshold, a local recombination algorithm is executed to update the node coordinates and boundary connections within 20 milliseconds / frame, and re-couple with the main grid through the non-matching interface interpolation algorithm to ensure that the global grid distortion rate is stabilized below 5%. The data of the recombination region is recorded in real time into the dynamic mapping relation table to form a closed-loop optimization link, solving the problem of low global reconstruction efficiency of traditional methods.

[0039] As can be seen from the above description, the electromagnetic fluid simulation modeling analysis method provided in this embodiment has the following technical effects:

[0040] Through the multi-channel deployment of the displacement sensor virtual nodes at the closing critical point, the maximum opening displacement point and the turning point of the movement direction, combined with the data consistency verification of the overlapping detection region in the opening acceleration section and the deceleration section, the displacement measurement error is eliminated. And by dynamically adjusting the width of the sliding average filter window according to the sampling frequency and resampling the non-uniform sampling data by the cubic spline interpolation algorithm, the time interval of the armature displacement data is made uniform, and the error is controlled within 0.05 mm. Based on the acceleration threshold, the sampling frequency is dynamically increased to 3 times the reference value to accurately capture the details of the high-speed movement of the armature. The hysteresis compensation coefficient is generated by the air-gap permeability difference at the same displacement point in the historical opening and closing operations, and the amplitude is corrected by combining the oscillation suppression factor at the closing critical point. A four-dimensional dynamic mapping relation table is constructed and a 7:3 weight ratio update mechanism is adopted to realize the closed-loop compensation of the permeability hysteresis effect.

[0041] In the electromagnetic field simulation, the three-dimensional transient Maxwell's equations are solved by adaptive grid meshing and the finite-difference time-domain method. By combining the magnetic field gradient threshold, the coordinate set of the dynamically updated region is screened, and the dynamic coupling boundary conditions including the electromagnetic force density distribution and the energy exchange threshold conditions are generated. And the displacement data resampling verification is triggered by the overlimit of the magnetic field strength deviation, and the permeability weight ratio is corrected in reverse to improve the accuracy of the boundary conditions.

[0042] In the fluid field simulation, the high-gradient plasma region is identified based on the magnetic field strength gradient and the electromagnetic force density distribution to generate deformed sub-grids. The seamless coupling of the main grid and the sub-grids is realized through the non-matching interface interpolation algorithm. By combining the armature displacement change rate, the deformation sensitivity and the deformation margin amplification mode are dynamically adjusted to predict the arc expansion direction and correlate with the hysteresis compensation coefficient to optimize the node deformation degree of freedom, and the maximum distortion rate is controlled within 15%. A deformation sensitivity suppression mechanism is superimposed at the closing critical point to reduce the grid distortion caused by contact bounce.

[0043] The grid structure of the non-arc region is fixed through the main grid topology. The deformation-sensitive region is screened based on the magnetic permeability gradient threshold and the displacement change rate threshold. The local recombination algorithm is executed for the sub-grids with the distortion rate exceeding the limit. The node coordinates and boundary connections are updated within 20 milliseconds per frame. The global grid distortion rate is stabilized below 5%. At the same time, the data characteristics of the recombination region are recorded in real time through the dynamic mapping relation table, forming a closed-loop link from data acquisition, magnetic permeability correction, grid deformation to local optimization, significantly improving the convergence speed and stability of the electromagnetic-fluid coupling simulation under dynamic conditions.

[0044] The above shows and describes the basic principles, main features and advantages of the present invention. Those skilled in the art should understand that the present invention is not limited by the above embodiments. The above embodiments and the descriptions in the specification are only preferred examples of the present invention and are not used to limit the present invention. Without departing from the spirit and scope of the present invention, the present invention will have various changes and improvements, and these changes and improvements all fall within the scope of the present invention claimed. The scope of the present invention claimed is defined by the appended claims and their equivalents.

Claims

1. An electromagnetic fluid simulation modeling and analysis method applied to AC contactors, characterized in that, The method steps are as follows: S1. Real-time collect armature displacement data through the displacement sensor virtual nodes, construct a dynamic mapping relationship table between the armature displacement data and the reference air-gap permeability, and the dynamic mapping relationship table includes a hysteresis compensation coefficient; S2. Dynamically update the reference air-gap permeability based on the armature displacement data and the hysteresis compensation coefficient in step S1, and generate dynamic coupling boundary conditions in the electromagnetic field simulation module; S3. According to the dynamic coupling boundary conditions in step S2, adopt the deformable sub-grid nesting technology in the fluid field simulation module, combine with the change rate of the armature displacement data in step S1, and implement elastic deformation control on the arc region grid to generate deformed sub-grids; Among them, the hysteresis compensation coefficient is synchronously correlated during the elastic deformation control process; S4. Fix the non-arc region grid through the main grid topology preservation module, screen the deformation-sensitive regions based on the dynamic mapping relationship table in step S1, and perform local grid reorganization calculations on the deformed sub-grids generated in step S3.

2. The electromagnetic fluid simulation modeling analysis method for AC contactors according to claim 1, characterized in that The setting positions of the displacement sensor virtual nodes in step S1 include the closed critical points, the maximum opening displacement points, and the turning points of the movement direction of the armature; each displacement sensor virtual node is configured with an independent data acquisition channel, and an overlapping detection region is set at the junction of the opening acceleration section and the deceleration section, and the measurement error is eliminated through the data consistency verification of adjacent displacement sensor virtual nodes.

3. The electromagnetic fluid simulation modeling analysis method for an AC contactor according to claim 1, characterized in that The specific process of collecting the armature displacement data includes: The displacement sensor virtual nodes first collect the original armature displacement data, and process the original armature displacement data using a sliding average filter. The width of the filter window is dynamically adjusted according to the sampling frequency of the original armature displacement data, and the unevenly sampled original armature displacement data is resampled at equal intervals through a cubic spline interpolation algorithm to generate a displacement data sequence with a uniform time interval as the armature displacement data.

4. The electromagnetic fluid simulation modeling and analysis method for an AC contactor according to claim 1, wherein The generation of the hysteresis compensation coefficient in step S1 includes extracting the difference in air-gap permeability at the same displacement points during the opening and closing processes, dynamically updating it in combination with the direction mark of the armature displacement data, and superimposing an oscillation suppression factor at the closed critical point to correct the amplitude of the hysteresis compensation coefficient.

5. The electromagnetic fluid simulation modeling and analysis method for AC contactors according to claim 1, characterized in that When dynamically updating the reference air-gap permeability in step S2, input the corrected reference air-gap permeability into the finite element solver, use the adaptive grid meshing technology to perform high-density grid division on the air-gap region, and solve the three-dimensional transient Maxwell's equations through the finite difference time domain method to generate dynamic coupling boundary conditions.

6. The electromagnetic fluid simulation modeling and analysis method for AC contactors according to claim 5, characterized in that, The dynamic coupling boundary conditions include the magnetic field strength gradient distribution, the electromagnetic force density distribution, and the dynamically updated region coordinate set; the dynamically updated region coordinate set is generated by screening through the magnetic field gradient threshold and is used to drive the local grid deformation calculation of the fluid field simulation module.

7. The electromagnetic fluid simulation modeling and analysis method for AC contactors according to claim 1, characterized in that When generating the deformed sub-grids in step S3, separate the arc region grid from the main grid according to the magnetic field strength gradient distribution and the electromagnetic force density distribution, establish a boundary data transfer channel between the deformed sub-grids and the main grid using the non-matching interface interpolation algorithm, and dynamically expand or contract the coverage range of the deformed sub-grids according to the magnetic field gradient threshold.

8. The application to the electromagnetic fluid simulation modeling and analysis method of an AC contactor according to claim 7, characterized in that, In step S3, the deformation sensitivity of the deformation sub-grid is adjusted in combination with the change rate of the armature displacement data. When the change rate of the armature displacement data exceeds the set change threshold, the deformation margin amplification mode is activated, and the deformation parameters of the deformation sub-grid are adjusted based on the change trend of the change rate of the armature displacement data to predict the expansion direction of the arc shape.

9. The electromagnetic fluid simulation modeling analysis method for AC contactors according to claim 7, characterized in that, In step S4, the screening of the deformation-sensitive area includes extracting the gradient value of the reference air-gap permeability and the change rate of the armature displacement data from the dynamic mapping relationship table. When the gradient value of the reference air-gap permeability and the change rate of the armature displacement data exceed the corresponding thresholds, they are marked as deformation-sensitive areas, and the local grid reorganization algorithm is executed based on the distortion rate data of the deformation sub-grid.

10. The electromagnetic fluid simulation modeling analysis method for AC contactors according to claim 9, characterized in that, The local grid reorganization algorithm re-divides the deformation sub-grid whose distortion rate exceeds the set grid distortion rate threshold, synchronously updates the node coordinates and boundary connection relationships of the deformation sub-grid, and recouples with the main grid through the non-matching interface interpolation algorithm.

Images