Simulation modeling analysis method applied to electromagnetic fluid of alternating current contactor

Abstract

The invention relates to the technical field of simulation modeling analysis, in particular to a simulation modeling analysis method applied to electromagnetic fluid of an alternating current contactor, which comprises the following steps: acquiring armature displacement data in real time through a displacement sensor, and constructing a dynamic mapping relation table containing hysteresis compensation coefficients; dynamically correcting the air gap permeability based on the displacement data and the compensation coefficient, generating an electromagnetic field dynamic coupling boundary condition, and realizing real-time matching of magnetic parameters and mechanical motion; in fluid field simulation, a deformable sub-grid nesting technology is adopted, elastic deformation control is carried out on grids in an electric arc area in combination with a displacement change rate, the influence of magnetic field distortion on the grids is synchronously corrected by associating hysteresis compensation coefficients, and the coupling problem of electric arc dynamic deformation and hysteresis nonlinearity is solved; and meanwhile, a non-arc region is fixed by utilizing a main grid topology keeping module, a deformation sensitive region is screened to execute local grid recombination, and high-precision dynamic simulation of an arc severe deformation region is realized on the premise of ensuring the 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, making it 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 inability to accurately reflect the remanence 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 in the arc shape, requiring frequent global mesh reconstruction, with low calculation efficiency and prone to numerical oscillation; while the existing dynamic mesh technologies can only make local adjustments, lacking 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, making it 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. 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

[0004] The purpose of the present invention is to provide an electromagnetic fluid simulation modeling analysis method applied to AC contactors to solve the problems proposed 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 compensation for the hysteresis effect, 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 grid reorganization and elastic deformation control.

[0005] To achieve the above purpose, the present invention provides the following technical solution: An electromagnetic fluid simulation modeling analysis method applied to AC contactors, and the method steps include: S1. By setting virtual nodes of displacement sensors at the closed critical point, maximum opening displacement point, and turning point of the moving direction of the armature, configuring independent data acquisition channels and overlapping detection areas, the measurement error is 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, generating a displacement data sequence with uniform intervals to solve the problem of insufficient capture of high-speed movement details. The hysteresis compensation coefficient is generated by combining the difference in air-gap permeability 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.

[0006] 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 the adaptive grid meshing technology, and the three-dimensional transient Maxwell's equations are solved by combining the finite-difference time-domain method to generate dynamic coupling boundary conditions including the magnetic field strength gradient distribution, electromagnetic force density distribution, and dynamically updated region coordinate set. The boundary element set to be refreshed in real time is screened through the magnetic field gradient threshold, driving the local grid deformation calculation of the fluid field simulation module to solve the problem that the traditional fixed grid cannot adapt to the dynamic air-gap change. 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 conditions, 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. 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, and the deformation parameters are adjusted based on the predicted arc expansion direction according to the change trend. The hysteresis compensation coefficient is synchronously associated to optimize the deformation degree of freedom of the nodes, solving the problem of simulation instability caused by grid distortion in the arc region.

[0007] S4. The gradient value of the reference air-gap permeability and the change rate of the armature displacement data are extracted through the dynamic mapping relationship table. When both exceed the corresponding thresholds, they are marked as deformation-sensitive regions. Based on the distortion rate data of the deformed sub-grid, the local grid reorganization algorithm is used to re-divide the sub-grid cells with distortion rates exceeding the grid distortion rate threshold, synchronously updating the node coordinates and boundary connection relationships, and re-coupling with the main grid through the non-matching interface interpolation algorithm. On the premise of keeping the main grid topology structure fixed, the problem of reduced simulation accuracy caused by excessive global grid distortion rate is solved.

[0008] Compared with the prior art, the beneficial effects of the present invention are as follows: Through the multi-position deployment of virtual nodes of displacement sensors and data consistency verification, combined with the moving average filtering and cubic spline interpolation algorithms, the displacement acquisition noise during the high-speed movement of the armature is eliminated, ensuring the data reliability of the reference air-gap permeability and the hysteresis compensation coefficient in the dynamic mapping table; based on the reference air-gap permeability corrected in real time, the adaptive mesh generation of the electromagnetic field simulation module is driven, and the magnetic field strength gradient and electromagnetic force density distribution are accurately transmitted to the fluid field through 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; The deformable sub-grid nesting technology is adopted to dynamically adjust the deformation sensitivity in combination with the displacement change rate, and the hysteresis compensation coefficient is synchronously correlated to compensate for the influence of the residual magnetism of the material, suppressing the local distortion caused by contact bounce, and maintaining 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 for the deformation-sensitive area, on the premise that the grid in the non-arc area is fixed, local recombination calculations are quickly performed on the sub-grid units with excessive distortion, avoiding the consumption of computing resources for global grid reconstruction, and improving the iterative efficiency and dynamic working condition adaptability of the electromagnetic-fluid two-way coupling simulation. BRIEF DESCRIPTION OF THE DRAWINGS

[0009] Figure 1 It is a schematic diagram of the method steps of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0010] The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to 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. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0011] 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: S1. Virtual nodes are set at key positions of the armature movement trajectory, the sampling frequency is dynamically adjusted according to the armature acceleration to capture the details of displacement changes, and the real-time displacement data is preprocessed through the moving average filtering and cubic spline interpolation algorithms; 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 to ensure that the acquisition accuracy error of the original armature displacement data is less than 0.05 mm; the original armature displacement data collected by the displacement sensor virtual nodes is subjected to moving average filtering, and the width of the filtering window is dynamically adjusted according to the displacement data sampling frequency; the cubic spline interpolation algorithm is used to resample the unevenly sampled displacement data sequence at equal intervals to generate a displacement data sequence with uniform time intervals as the armature displacement data.

[0012] 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; the reference displacement data sampling frequency is restored in the uniform motion stage, and the reference displacement data sampling frequency is set to 1 kHz. Calculate the hysteresis compensation coefficient based on the air gap permeability difference at the same displacement point in historical switching 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, a dynamically updated hysteresis compensation coefficient is generated; an oscillation suppression factor is superimposed at the closing critical point to correct the amplitude of the hysteresis compensation coefficient. 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. Construct a four-dimensional dynamic mapping relationship table including armature displacement data, reference air gap permeability, hysteresis compensation coefficient, and motion state; after each switching operation, update the dynamic mapping relationship table with the new armature displacement data and hysteresis compensation coefficient according to a 7:3 weight ratio.

[0013] 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 moving 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 are 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 are updated according to the weight ratio to solve the problems of displacement-permeability mapping distortion and insufficient hysteresis effect compensation.

[0014] 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: , where represents the corrected reference air-gap permeability; represents the reference air-gap permeability; represents the hysteresis compensation coefficient; Input the corrected reference air-gap permeability into the finite element solver, use the adaptive mesh refinement technology to perform high-density mesh division on the air-gap area, solve the three-dimensional transient Maxwell's equations by the finite-difference time-domain method, and 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. Synchronously extract the temperature gradient and ion concentration distribution data of the arc plasma, and define the threshold condition of the flow field energy exchange in combination with the flow field continuity equation. Finally, a dynamic coupling boundary condition including the magnetic field strength gradient, electromagnetic force density distribution, and dynamically updated region coordinate set is generated. The dynamically updated region coordinate set is screened by the magnetic field gradient threshold and marked as a set of boundary elements that need 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 bidirectional coupling. At the same time, verify the accuracy of the boundary conditions by comparing historical simulation data. When it is detected that the magnetic field strength deviation exceeds 5%, trigger the resampling verification of the armature displacement data and reverse-correct the weight ratio of the reference air-gap permeability in the dynamic mapping relationship table to form a closed-loop feedback mechanism.

[0015] In step S2, to address the problem of cumulative errors 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 table. After being input into the finite element solver, the adaptive mesh refinement technique is used to densely divide the air-gap region. The three-dimensional transient Maxwell's equations are solved by the finite-difference time-domain method to calculate the distribution of magnetic field strength gradient. Based on the Maxwell stress tensor, the electromagnetic force density distribution is derived. Combining with the continuity equation of the flow field, the energy exchange threshold is defined to generate a dynamic coupling boundary condition that includes the magnetic field gradient, electromagnetic force density, and the coordinate set of the dynamically updated region. The real-time refresh region is screened through the magnetic field gradient threshold, and the displacement resampling verification and the reverse correction of the magnetic permeability weight are triggered based on the over-limit deviation of the magnetic field strength, forming a closed-loop feedback mechanism to solve the problem of cumulative errors in magnetic field calculation caused by traditional empirical formulas.

[0016] S3. Based on the dynamic coupling boundary condition in step S2, the deformable sub-grid nesting technique is used in the fluid field simulation module to generate deformable sub-grids. Among them, according to the magnetic field strength gradient and electromagnetic force density distribution in the dynamic coupling boundary condition, the arc region grid is separated from the main grid. Based on the magnetic field strength 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 deformable sub-grid. The non-matching interface interpolation algorithm is used to establish the boundary data transfer channel between the deformable sub-grid and the main grid to ensure the seamless coupling of electro-magnetic and fluid field parameters. The coverage range of the deformable 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 deformable sub-grid to form a deformable sub-grid that changes in real time with the magnetic field. Combined with the change rate of the armature displacement data in step S1, the deformation amplitude of the deformable 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 deformable 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 deformable sub-grid. For example, the maximum allowable deformation amount of the nodes in the deformable sub-grid is increased to 2 times 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 deformable sub-grid are adjusted in advance to ensure that the deformable sub-grid is strictly synchronized with the movement state of the armature. By correlating the hysteresis compensation coefficient in step S1, the generation accuracy of the deformable sub-grid is optimized. According to the magnitude of the hysteresis compensation coefficient, the deformation freedom degree of the nodes in the deformable sub-grid is proportionally increased 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 to suppress the sub-grid distortion caused by contact bounce by reducing the deformation sensitivity, so that the generated deformable sub-grid maintains a stable shape under dynamic conditions, and the maximum distortion rate is controlled within 15%.

[0017] In step S3, for the problem of simulation instability caused by arc grid distortion, based on 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. The deformation sensitivity is dynamically adjusted in combination with the change rate of armature displacement data. When the change rate exceeds the threshold, the deformation margin amplification mode is activated (for example, the maximum deformation amount is increased to 2 times the benchmark), and the deformation parameters are adjusted by predicting the arc expansion direction. The node deformation degrees of freedom are optimized by synchronously correlating the hysteresis compensation coefficient to compensate for the arc shape hysteresis caused by the remanence of the material. An oscillation suppression mechanism is superimposed at the closing critical point 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 of the existing dynamic grid technology.

[0018] S4. The topological structure of the non-arc region grid is fixed by 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 benchmark air-gap permeability gradient and the change rate of 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: The gradient values of the benchmark air-gap permeability and the change rate of armature displacement data at each displacement point are extracted from the dynamic mapping relationship table. The permeability gradient threshold and the displacement change rate threshold are set. 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. For the screened deformation-sensitive regions, the local grid recombination algorithm is used to re-divide the deformed sub-grids whose distortion rate exceeds the grid distortion rate threshold, and the node coordinates and boundary connection relationships of the deformed sub-grids are updated synchronously, while maintaining the topological consistency between the main grid and the non-arc region grid. The recombined deformed sub-grids are re-coupled 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 recombination calculation time is controlled within 20 ms / frame, the global grid distortion rate is stabilized below 5%, and the change characteristics of the permeability and displacement data in the recombination region are recorded in real time through the dynamic mapping relationship table to form a closed-loop optimization link.

[0019] In step S4, to address the problems of low global reconstruction efficiency and local distortion diffusion, the main grid topology preservation module is used to fix the grid structure in the non-arc region. Based on the dynamic mapping relation table, the deformation-sensitive regions with excessive magnetic permeability gradient and displacement change rate are screened. 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 per frame, and they are recoupled 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 recombined 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 in traditional methods.

[0020] As can be seen from the above description, the electromagnetic fluid simulation modeling analysis method provided in this embodiment has the following technical effects: 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. By dynamically adjusting the width of the sliding average filter window with the sampling frequency and resampling the non-uniform sampling data through the cubic spline interpolation algorithm, the time interval of the armature displacement data is made uniform, and the error is controlled within 0.05 millimeters. 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 based on the air gap magnetic 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 achieve the closed-loop compensation of the magnetic permeability hysteresis effect.

[0021] 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 to generate a dynamic coupling boundary condition including the electromagnetic force density distribution and the energy exchange threshold condition. When the magnetic field strength deviation exceeds the limit, the displacement data resampling verification is triggered to reversely correct the magnetic permeability weight ratio and improve the accuracy of the boundary condition.

[0022] 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 achieved through the non-matching interface interpolation algorithm. By dynamically adjusting the deformation sensitivity and the deformation margin amplification mode in combination with the armature displacement change rate, the arc expansion direction is predicted and the node deformation degree of freedom is optimized by associating with the hysteresis compensation coefficient, controlling the maximum distortion rate within 15%. At the closing critical point, a deformation sensitivity suppression mechanism is superimposed to reduce the grid distortion caused by contact bounce.

[0023] The grid structure of the non-arc region is fixed by 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 on the sub-grid with an excessive distortion rate. 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.

[0024] 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 the armature displacement data through the displacement sensor virtual node, 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 a dynamic coupling boundary condition in the electromagnetic field simulation module; S3. According to the dynamic coupling boundary condition in step S2, adopt the deformable sub-grid nesting technology in the fluid field simulation module, and combine the change rate of the armature displacement data in step S1 to implement elastic deformation control on the arc region grid and generate a deformed sub-grid; Among them, the hysteresis compensation coefficient is synchronously associated during the elastic deformation control process; S4. Fix the non-arc region grid through the main grid topology preservation module, screen the deformation-sensitive region based on the dynamic mapping relationship table in step S1, and perform local grid recombination calculation on the deformed sub-grid generated in step S3.

2. The electromagnetic fluid simulation modeling analysis method for an AC contactor according to claim 1, wherein The setting positions of the displacement sensor virtual nodes in step S1 include the closed critical point, the maximum opening displacement point and the turning point 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 and analysis method for an AC contactor according to claim 1, wherein The specific process of collecting the armature displacement data includes: The displacement sensor virtual node first collects the original armature displacement data, and processes the original armature displacement data by 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 uniform time intervals as the armature displacement data.

4. The electromagnetic fluid simulation modeling and analysis method for AC contactors according to claim 1, characterized in that The generation of the hysteresis compensation coefficient in step S1 includes extracting the air-gap permeability difference at the same displacement point 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, wherein 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 a dynamic coupling boundary condition.

6. The electromagnetic fluid simulation modeling and analysis method for AC contactors according to claim 5, characterized in that, The dynamic coupling boundary condition includes 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-grid 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-grid and the main grid by using a non-matching interface interpolation algorithm, and dynamically expand or contract the coverage range of the deformed sub-grid according to the magnetic field gradient threshold.

8. The electromagnetic fluid simulation modeling analysis method for an AC contactor according to claim 7, wherein In step S3, the deformation sensitivity of the deformed 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 expansion direction of the arc shape is predicted based on the change trend of the change rate of the armature displacement data to adjust the deformation parameters of the deformed sub - grid.

9. The electromagnetic fluid simulation modeling and 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 relation 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 deformed 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 deformed sub - grid whose distortion rate exceeds the set grid distortion rate threshold, synchronously updates the node coordinates and boundary connection relationships of the deformed sub - grid, and recouples with the main grid through the non - matching interface interpolation algorithm.

Images