79 results about "Differential algebraic equation" patented technology

The invention discloses a parallel iteration Runge-Kutta method for simulating a second-order nonlinear differential equation. The parallel iteration Runge-Kutta method comprises the steps that a simulated ordinary differential equation is converted into a target nonlinear implicit differential equation, and a system differential algebraic equation is established; judging whether a system matrix and a decomposition iteration matrix in the system differential algebraic equation are updated or not; a four-order extrapolation method is adopted to obtain a predicted differential variable derivative, and if a prediction result is not converged, parallel iterative operation is carried out through a four-stage four-order Runge-Kutta method; after the iteration process is converged, judging whether an operation result accords with a preset error range, if not, updating a system matrix according to a step size reduction strategy and re-decomposing an iteration matrix, and if yes, updating the operation result; and storing and displaying an operation result when judging that simulation is finished. The algorithm stability and the calculation efficiency of the rigidity problem can be considered.

The invention relates to a fault current analysis method based on first-order Taylor expansion, and the method comprises the following steps: 1, determining a system control strategy which is a network construction type control strategy and comprises a power loop, a virtual impedance link and a voltage and current loop; the power loop comprises an active power control loop and a reactive power control loop; 2, constructing a differential algebraic equation set model capable of representing the dynamic state of the system; 3, according to the characteristic that the integral value does not suddenly change at the moment before and after the fault occurs, the steady state before the fault is utilized to solve and obtain an instant state variable value after the fault, and then an instant dq-axis current value after the fault and an instant current amplitude value after the fault are obtained; 4, performing first-order Taylor expansion on the constructed differential algebraic equation set at the solved post-fault instant state variable value and the dq-axis current value, and converting a nonlinear differential algebraic equation set model into a linear differential algebraic equation set; and 5, analyzing the fault current.

The invention relates to the technical field of water supply pipe networks, and discloses a multi-pump-room collaborative intelligent water supply frequency conversion regulation and control method, which comprises the following steps: initializing sensor calibration and data acquisition, generating a standardized data packet by using Kalman filtering, constructing a differential algebraic equation model, and online updating a pipe network topological characteristic matrix in combination with a mixed basis function observer. Layered optimization control and frequency conversion cooperative regulation and control are driven, 30-100ms-level real-time solution of a differential equation is realized by means of a hardware acceleration module, a three-dimensional pressure cloud picture is finally generated, and early warning is triggered based on a current harmonic threshold value. The transient fault capturing capability is enhanced by dynamically adjusting the data acquisition frequency, and the problem of signal distortion caused by traditional fixed sampling is solved; resource allocation is optimized and calculated in combination with a pipe network topology model dimension reduction technology, and calculation power waste in large-scale monitoring is reduced; and dynamic reconstruction of the model is driven by abnormal data to realize system self-adaption, so that manual parameter adjustment deviation is avoided, and the prediction precision and the diagnosis reliability are remarkably improved.

The invention discloses a copper recovery processing system and method for flameless combustion of copper-containing sludge, and relates to the technical field of intelligent control. A distributed sensor network is deployed, multi-source combustion data in a combustion furnace are acquired in real time, a hybrid convolution-long and short-term memory network model is constructed, cross-modal feature fusion is performed on the multi-source combustion data, and the copper recovery processing system and method for flameless combustion of the copper-containing sludge are obtained. The method comprises the following steps: extracting a multi-dimensional combustion feature vector, correcting system parameters of a differential-algebraic equation on line based on the multi-dimensional combustion feature vector, solving a real-time system state variable after the system parameters are corrected, and constructing a constraint condition by using the system state variable. A multi-agent reinforcement learning model is constructed based on the multi-source combustion data and the multi-dimensional combustion feature vectors, dynamic combustion parameters are optimized in parallel by adopting a multi-agent depth deterministic strategy gradient algorithm, and the dynamic combustion parameters are sent to a combustion execution mechanism for combustion parameter updating; the stability and the environmental protection performance in the combustion process are obviously improved.

The invention discloses a Bernstein polynomial-based multi-region frequency safety unit scheduling method and device, and the method comprises the steps: constructing a multi-region frequency response model based on the dynamic primary frequency modulation response of a synchronous generator, and the optimal distribution of a speed regulator dead zone and frequency modulation reserve; carrying out approximation on an integral-differential algebraic equation of the multi-region frequency response model through a Bernstein polynomial to obtain a dynamic algebraic expression of the multi-region frequency, and constructing a multi-region frequency security constraint according to the dynamic algebraic expression of the multi-region frequency; embedding the multi-region frequency security constraint into a unit commitment problem, and constructing a multi-region frequency security constraint unit commitment model; decomposing a mixed integer programming problem of the multi-region frequency security constraint unit commitment model into a main problem and a sub-problem through an adaptive double-layer decomposition algorithm; and solving to obtain a multi-region frequency safety unit scheduling scheme. According to the method, the frequency dynamic difference between the regions is accurately captured, and the frequency safety of the high-proportion new energy power system is improved.

Compiler systems and methods that solve the DAE transformation problem by simulating DAEs using general-purpose compilers are disclosed herein. The compiler systems and methods extend an existing general-purpose intermediate representation (IR) with a minimal set of intrinsics to represent the model function, which allows leveraging existing compiler technology. A first compiler input is a model of a physical system from source modeling languages. A second compiler input is a first IR for each subroutine specific to the model. The model is compiled to a second IR that is universal across multiple source modeling languages. The IR maps each primitive of the source modeling language to a subroutine call specific to the source modeling language from which the primitive is derived. The first and second IRs are combined into a third IR that describes the model in terms of intrinsics of the third IR.

The invention belongs to the technical field of photoelectricity, and relates to a detector moving laser irradiation damage assessment method and system. The evaluation method comprises the following steps: determining mobile laser parameters and detector structure material attribute parameters; constructing a detector solid heat transfer transient model; setting physical field boundary conditions of the solid heat transfer transient model and performing grid division; setting solver parameters, specifically, inputting the moving laser pulse width and the moving laser frequency into a pre-trained linear neural network model, and calculating by applying an expression of the model to obtain the solver parameters; and based on the solver parameters, solving the solid heat transfer transient model in the solver, judging the thermal damage condition of the detector by calling domain ordinary differential and differential algebraic equation modules, and if it is determined that the detector is damaged, calculating the damage area and volume of the detector to obtain the spatial distribution condition of the damage area of the detector.

The invention discloses a multi-converter parameter collaborative optimization method for improving the small interference stability of an active power distribution network, and belongs to the technical field of power systems, and the method comprises the steps: building a differential algebraic equation model, and carrying out the linearization of a steady-state operation point, and obtaining a state space model; and a control parameter optimization problem taking the minimization of the spectrum intercept construction as a target is defined, and control parameters are a GFL converter phase-locked loop parameter and a GFM converter virtual damping parameter. First-order sensitivity and second-order sensitivity of eigenvalues to control parameters are calculated based on a matrix perturbation theory and are converted into local quadratic constraint sub-problems. And performing convexity solution by adopting a positive semidefinite relaxation method, adjusting the step length in combination with a line search mechanism, iteratively updating control parameters, and outputting an optimization result after a convergence condition is met. According to the method, the spectral intercept is reduced through multi-converter parameter collaborative optimization, and the weak stability characteristic is improved; the convergence and solvability are improved by utilizing matrix perturbation and positive semidefinite relaxation; and the stability is enhanced and optimized by adopting line search step length adjustment.

The invention provides a three-dimensional geological modeling method, system and equipment based on a particle dynamics system. The method comprises the following steps: carrying out data preprocessing on a drilling data set of a to-be-built geological model region to obtain a stratum attribute known point set; discretizing the to-be-built geologic model area, and performing stratum parameter initialization processing on each spatial discrete coordinate point in the to-be-built geologic model area according to the stratum attribute known point set to obtain initial stratum parameters of a to-be-built geologic model; based on the initial stratum parameters and the drilling data set, a weight function expressed by spatial distribution information is established; establishing a link relation information matrix according to grid division and section discontinuity information of the to-be-established geological model area; based on a particle dynamics theory, establishing a system dynamics equation in a differential algebraic equation form by applying the initial stratum parameters, the weighting function and the link relation information matrix; and solving the system kinetic equation to establish a three-dimensional geological model of the to-be-established geological model area.

The invention discloses a deep and far sea wind power time-varying impedance model grid connection analysis method based on a time domain residual error method, relates to the technical field of power system modeling simulation, and can efficiently and accurately analyze deep and far sea wind power grid connection. The method comprises the following steps: establishing a differential algebraic equation set containing high-capacity submarine cable impedance and a plurality of nonlinear links, and constructing a recursive coefficient matrix; dividing a preset global time domain into a plurality of time domain subintervals, and solving the recursion coefficient matrix in each time domain subinterval by referring to a universal analytic solution represented by a power series explicit operator to generate a global time domain explicit solution corresponding to the preset global time domain; performing discrete Fourier transform on the global time domain explicit solution to obtain a time-frequency domain spectrum result; and constructing a time-varying impedance model according to the differential algebraic equation, setting the time-varying impedance model by adopting a time-frequency domain spectrum result, and analyzing the deep and far sea wind power grid-connected control system by utilizing the set time-varying impedance model.

The invention discloses a neural network construction method and system for solving a high-index differential algebraic equation, and the method achieves the efficient and precise solving of the high-index differential algebraic equation through the preprocessing of an equation, the construction of a physical information neural network model, the definition of a loss function, the setting of training parameters and data, and the iterative training through an L-BFGS algorithm. And finally, the trained neural network model is used for solving a differential algebraic equation of the constrained mechanical arm system, variable values are uniformly selected in a definition domain and input into the trained model, so that a solution of the system is obtained, in addition, a predicted solution and an actual solution of the neural network are compared and analyzed, and the optimal solution is obtained. By means of the mode, the nonlinear high-index differential algebraic equation problem in the constrained mechanical arm system can be efficiently and accurately solved.

The invention relates to the field of electric power, discloses a thermal power generating unit APS automatic start-stop full-process control method and system, and is used for solving the contradiction among safety, economy and rapidity of a traditional APS system. The method comprises the steps of establishing a differential algebraic equation system based on unit physical parameters to generate a unified model constraint set; generating a reference trajectory representing energy conversion efficiency and power in a differential flat space according to the start-stop target curve; predicting a thermal stress change rate through a fractional order differential operator based on the metal temperature to generate a safety margin instruction; generating an energy coordination correction instruction according to the Hamiltonian energy network balance point of the boiler heat storage and steam turbine kinetic energy solving port; and mapping the reference trajectory into a control instruction, integrating thermal stress constraint and energy coordination correction, and then distributing and executing. Safe and efficient start-stop control under full-process physical constraints is achieved, start-stop safety and economy are improved, and the service life of equipment is prolonged.

The invention relates to the technical field of power system simulation, and discloses a power system time sequence simulation system considering security constraints. The system comprises a power grid topology analysis unit, a dynamic security domain calculation unit, a time sequence constraint modeling unit, a simulation track generation unit and a check feedback unit. The power grid topology analysis unit obtains element connection relations and electrical parameters, and generates a topology description file containing a node admittance matrix; the dynamic safety domain calculation unit obtains the voltage stability limit and the power angle stability boundary of each working condition, and forms a safety domain boundary set; a time sequence constraint modeling unit combines the set with a predefined time interval to build a security constraint time sequence model containing a transient process; the simulation trajectory generation unit solves a differential algebraic equation set by using an implicit trapezoidal integral method, and outputs a state variable trajectory with a timestamp; and the check feedback unit compares the track with the security domain boundary, marks an abnormal segment and generates a correction instruction. The system realizes fusion of security constraint and time sequence simulation.

The invention discloses a unit commitment optimization method and system containing transient voltage security constraints, and belongs to the field of transient security evaluation and control of a power system. The method comprises the following steps: constructing a unit commitment model according to transient voltage safety constraint conditions; simplifying an embedded differential-algebraic equation set constraint in the unit commitment model through a trajectory sensitivity algorithm to obtain an optimization model under transient voltage under multiple time scales; decomposing the optimization model according to a generalized Benders decomposition algorithm to obtain a mixed integer linear programming problem which does not contain alternating current power flow calculation and a nonlinear programming problem which meets transient voltage constraint after iteration; and respectively solving the mixed integer linear programming problem and the nonlinear programming problem through a solver to obtain an optimal solution of the unit commitment model. Therefore, by simplifying the constraint of the embedded differential-algebraic equation set, the model is decomposed into the main problem and the sub-problem with different emphasis, and the main problem and the sub-problem are solved respectively, so that the efficiency of solving the original problem is greatly improved.

The invention provides a dynamic grouping method and device for a heterogeneous direct current delivery system, and relates to the technical field of dynamic analysis and control of a power system, and the method comprises the steps: deducing a linear time-varying variation equation which is satisfied by the trajectory sensitivity of a system state trajectory to a system parameter vector based on a parameterized differential-algebraic equation model; calculating the track sensitivity of each device to a system parameter vector, projecting the track sensitivity to a common connection point corresponding to each device, and constructing a unified sensitivity feature vector of each device; calculating the dynamic similarity between the devices based on the unified sensitivity feature vector, and constructing a mixed similarity matrix in combination with the electrical distance similarity between the devices; and feature decomposition and feature vector clustering are carried out based on the mixed similarity matrix, and dynamic clustering of heterogeneous equipment is realized. According to the dynamic grouping method provided by the invention, unified grouping of heterogeneous equipment such as a synchronous machine and a converter can be realized in a heterogeneous power system containing high-proportion new energy and direct-current power transmission.

The application relates to the technical field of power electronics, and discloses a press-contact type IGBT device packaging structure optimization method and system. The method comprises the following steps: constructing an equivalent model of a press-contact type IGBT device; constructing a differential algebraic equation for uniformly describing the equivalent model of the press-contact type IGBT device and circuit port variables; converting the differential algebraic equation into a nonlinear algebraic equation group through time domain discretization, and constructing an iteration matrix with a global state vector as a derivation object according to a solving mode of the nonlinear algebraic equation group; performing low-dimensional subspace projection on the iteration matrix based on a Krylov subspace, and combining residual minimization to solve a linear subproblem of the iteration matrix, so that current uniformity characteristic results of the press-contact type IGBT device are obtained; and the packaging structure of the press-contact type IGBT device is optimized based on current density distribution reflected by the current uniformity characteristic results. The application can improve the current uniformity and operation reliability of the system.

The invention discloses a multi-stage separation unit initialization method and system based on a virtual transient model and application, and relates to the technical field of chemical process simulation, and the method comprises the steps: constructing a virtual transient differential equation of a multi-stage separation unit; constructing a virtual differential variable associated equation of the multi-stage separation unit; constructing a virtual transient model of the whole tower of the multistage separation unit; solving the virtual transient model; and outputting the initial estimation value. According to the multi-stage separation unit initialization method and system based on the virtual transient model and the application, when a dynamic process differential algebraic equation model is constructed, the use of a complex hydraulics association formula and structural information of the separation unit is avoided, the time-varying item is virtual, the method does not need to be based on a real physical process, and the time-varying item is saved. However, the integral is completely equivalent to the algebraic equation model when reaching the steady state, so that the initialization process of the multi-stage separation unit based on the virtual transient model is easier to solve, and the method has the advantages of high calculation efficiency and higher robustness.

The invention discloses a distributed photovoltaic aggregation frequency modulation standby optimal configuration method considering node frequency security constraints under transmission and distribution cooperation. The method comprises the following steps: describing uncertainty of photovoltaic aggregation frequency modulation capability by using a Gaussian mixture model, and establishing a power transmission network node frequency response model containing frequency space distribution characteristics; taking minimization of system power generation and frequency modulation reserve cost as a target, and combining opportunity constraint planning to construct an optimal configuration model; a response model is reduced based on a frequency divider theory, a Bernstein polynomial approximates frequency dynamics, and a differential algebraic equation is converted into an algebraic equation; then, binary expansion and Gaussian mixture model linear invariance are utilized to convert the model into a mixed integer linear programming model, and Bernstein polynomial convex hull property tightening constraint is combined; and finally, solving the model to obtain the optimal photovoltaic aggregation frequency modulation reserve capacity and control parameters. The frequency space distribution is accurately described, the problem of insufficient frequency safety risk identification is solved, and the frequency safety level of a power system is improved.

The application discloses a multi-region frequency safety unit scheduling method and device based on Bernstein polynomials, and the method is based on the optimal distribution of the dynamic primary frequency modulation response of a synchronous generator, a speed governor dead zone and a frequency modulation backup, and a multi-region frequency response model is constructed; the integral-differential algebraic equation of the multi-region frequency response model is approximated through Bernstein polynomials to obtain an algebraic expression of the multi-region frequency dynamics, and a multi-region frequency safety constraint is constructed according to the algebraic expression; the multi-region frequency safety constraint is embedded into a unit combination problem to construct a multi-region frequency safety constraint unit combination model; through a self-adaptive double-layer decomposition algorithm, a mixed integer programming problem of the multi-region frequency safety constraint unit combination model is decomposed into a main problem and a sub-problem; and through solving, a multi-region frequency safety unit scheduling scheme is obtained. The application accurately captures the frequency dynamic difference between regions and improves the frequency safety of a high-proportion new energy power system.

The application provides a multi-inverter hybrid system transient stability boundary delineation method based on reverse integration, and belongs to the field of electrical engineering. The delineation method combines a circuit model of a grid-constructed inverter, a grid-following inverter and an infinite grid hybrid system and a virtual synchronous machine mathematical model of the grid-constructed inverter, derives a second-order differential-algebraic equation set of the hybrid system, further obtains system steady-state equilibrium points by solving the equation set and classifies, and reversely integrates from the type I unstable equilibrium point to obtain a transient stability boundary of the hybrid system. The application is simple and intuitive, has a solid theoretical basis, can quickly and accurately identify the transient stability boundary of the grid-connected inverter system, and has important significance for analyzing the transient stability of the multi-grid-connected inverter hybrid system.

The invention provides an electric power system step backspacing simulation method and system considering switching events, and belongs to the technical field of electric power system electromechanical transient simulation, and the method comprises the steps: executing transient simulation step by step, carrying out the alternate iteration of a differential algebraic equation set, scanning the switching events of a DC system in the iteration process, scanning other switching events after the iteration convergence, and carrying out the step backspacing simulation. Whether step backspacing needs to be carried out or not is checked according to a backspacing criterion of the step backspacing strategy; if so, carrying out a rollback operation and adjusting the simulation step length according to a rollback simulation strategy; switching the scanned switching event, and keeping the state variable of the system unchanged after the switching is completed; according to the scheme, whether the next step of simulation continues to be executed or not is checked, if simulation continues to be executed, the step length of the next step needs to be adjusted according to the criterion of step length increase and step length decrease, and the scheme has good technical compatibility and also has excellent expansibility.

The invention provides a reverse solution method, device and equipment for a vehicle motion model, relates to the technical field of dynamics and control, and aims to solve the problems that in the prior art, a forward solution model and a reverse solution model cannot be shared, and reverse solution universality is poor. The method comprises the following steps: acquiring dynamic data of a vehicle motion model, and establishing a vehicle motion model equation set according to the dynamic data; based on a to-be-solved target, target parameters are added into the vehicle motion model equation set, and an implicit differential algebraic equation is added for the vehicle motion model equation set according to the target parameters and the dynamic data; and performing reverse solution on the vehicle motion model by using the vehicle motion model equation set after the target parameter and the implicit differential algebraic equation are added.

The invention belongs to the technical field of automatic driving, and particularly relates to a position and posture derivation method and device of a tractor-trailer system, a computer readable storage medium and electronic equipment. The method comprises the following steps: determining an initial pose state of the tractor-trailer system; acquiring a real-time position of a tractor, and determining the real-time position of the tractor as a target position; based on a preset trajectory tracking controller, the initial pose state and the target position are processed, and the linear speed and the steering angle of the tractor are obtained; and based on a preset differential algebraic equation solver, the initial pose state, the linear speed and the steering rudder angle are processed, and the target pose state of the tractor-trailer system is obtained. According to the invention, accurate pose derivation can be carried out on the tractor-trailer system, so that the automatic driving technology can be conveniently applied to the tractor-trailer system.

A system for solving nonlinear problems in implicit differential equations comprises hardware processors that transform differential equation systems into a unified canonical form enabling efficient numerical solving across multiple integration methods. The processors represent differential equations in a reduced canonical form G(x; v1, v2, γ, c)=v1+Mx−γ(h) f(x+v2, tn+ch) incorporating solution vectors, integrator parameters, mass matrices, and time parameters. The system generates model-specific mapping functions N(x)=y and N−1(y, γ, v1, v2)=x that transform between differential algebraic equation states and compressed nonlinear solver states, eliminating algebraically redundant variables through symbolic analysis. This compressed representation reduces computational complexity from O(n3) to O(m3) where m<n, while integrator-specific functions v1, v2, and γ(ρh) encode method-specific discretization parameters.

The present invention belongs to the field of optoelectronic technology and relates to a method and system for evaluating the damage of a detector caused by moving laser irradiation. The evaluation method includes: determining the parameters of the moving laser and the parameters of the detector's structural material properties; constructing a transient solid heat transfer model of the detector; setting the physical field boundary conditions of the transient solid heat transfer model and performing mesh division; setting the solver parameters, specifically: inputting the moving laser pulse width and the moving laser frequency into a pre-trained linear neural network model and using the expression of the model to calculate the solver parameters; based on the solver parameters, solving the transient solid heat transfer model in the solver, and determining the thermal damage condition of the detector by calling the domain ordinary differential and differential algebraic equation module. If it is determined that the detector has been damaged, calculate the damage area and volume of the detector to obtain the spatial distribution of the damaged area of the detector.

The invention discloses an electrical integrated energy system fault analysis method based on semi-implicit simulation, and belongs to the field of integrated energy system simulation analysis. Specifically, in order to quantify the propagation velocity of the pipe explosion and leakage faults of the natural gas system and analyze the propagation mechanism of the faults in an electrical bidirectional coupling system, the method comprises the following steps: firstly, establishing a natural gas system differential algebraic equation model for describing the pipe explosion and leakage faults based on a characteristic line method; then, in consideration of the high-dimensional, nonlinear and rigid characteristics of the electrical bidirectional coupling system, a semi-implicit dynamic simulation method based on a rigid precise Rosenbrock format is provided; and finally, a key event positioning strategy based on a continuous Runge-Kutta algorithm is used to accurately position a key simulation event. According to the invention, high efficiency and robustness of simulation can be ensured. In combination with a key event positioning strategy, the fault analysis of the electrical integrated energy system can be efficiently and accurately realized.

A kind of electromechanical coupling dynamics modeling method of electric vertical take-off and landing aircraft based on differential algebraic equation belongs to electric vertical take-off and landing aircraft technical field.First, the motor dynamics equation in eVTOL electric propulsion system is established according to motor parameters;Second, according to the geometric configuration and inertia parameters of eVTOL, the dynamics equation of fuselage multibody system is established, and the aerodynamic model is applied;Third, the motion driving equation is used to couple based on the general mechanical load calculation scheme of Lagrange multiplier;Finally, the electromechanical coupling dynamics model of eVTOL is established, which is discretized by trapezoidal-Newmark hybrid discrete format, and efficient solution is realized based on Newton iteration method, and it is verified by eVTOL vertical take-off.The invention avoids the defects of traditional simplified model ignoring motor transient characteristics and non-coupling modeling without real-time load feedback, improves the motion response and energy consumption calculation accuracy, and provides reliable dynamics modeling basis for eVTOL energy consumption optimization and commercial design.

The invention discloses an electric power system load flow calculation method based on a rigid accurate Rosenbrock format, and relates to the field of electric power system load flow analysis. The method comprises the following steps: converting an alternating current load flow calculation model of a power system into a group of differential algebraic equations; according to a given voltage phase angle and an amplitude iteration initial value, calculating a differential algebraic equation state variable and an algebraic variable initial value; solving the differential algebraic equation set according to the obtained initial values of the state variable and the algebraic variable by using a rigid accurate Rosenbrock format until the error of the original power system load flow calculation model is smaller than a given error limit; and finally deriving a voltage phase angle and amplitude calculation result. According to the method, the characteristics of non-iteration, high efficiency and high numerical stability of the rigid accurate Rosenbrock format are combined, and the convergence performance of power flow calculation of the power system is improved. The method is suitable for secondary checking of a non-convergent scene of a traditional load flow calculation method, and fast robust convergence of load flow calculation in an ill-conditioned load flow initial value scene is guaranteed.

The invention belongs to the technical field of thermal hydraulics, and particularly relates to a Modelica-based two-fluid six-equation and differential algebraic equation coupling system compiling and solving device and a Modelica-based two-fluid six-equation and differential algebraic equation coupling system compiling and solving method. Comprising a compiling stage function system and a solving stage function system. The method has the beneficial effects that the problem of coupling, compiling and solving of a two-fluid six-equation system and a differential algebraic equation system under a unified modeling framework based on Modelica is solved, and a compiling and solving device with complete functions is provided; the whole process compiling solution from topological structure generation, coefficient matrix assembly, time step size control, source term function calculation, matrix equation solution to convergence judgment based on the quality error of the two-fluid six-equation system is realized; the whole-process compilation solution from grammar semantic check, planarization processing, compilation analysis to nonlinear solution of the differential algebraic equation system is realized, and the real-time data communication and interaction of the two-fluid six-equation system and the differential algebraic equation system under the unified solution sequence planning are realized.

The present invention discloses a method, system, terminal device and storage medium for transient stability assessment of a multi-machine system based on a hyper-tangent plane. The method constructs a differential algebraic equation group of the multi-machine system according to the state quantity of the multi-machine system, calculates the equilibrium point and Jacobian matrix of the multi-machine system according to the differential algebraic equation group of the multi-machine system, and then determines a class of unstable equilibrium points and Jacobian matrices thereof. The hyper-tangent plane equation is calculated according to the Jacobian matrix of the class of unstable equilibrium points. The nonlinear characteristics of the multi-machine system are analyzed by the hyper-tangent plane equation and the fault trajectory to calculate the critical fault removal time of the multi-machine system, and the transient stability of the multi-machine system is evaluated according to the critical fault removal time, thereby improving the accuracy of the transient stability assessment of the multi-machine system, ensuring the stable operation of the multi-machine system, and solving the problem that the transient stability assessment result of the multi-machine system is low due to the lack of analysis of the nonlinear characteristics of the multi-machine system in the time domain simulation method.