108 results about "Analytical expressions" patented technology

The ELM-based frequency safety constraint linearization optimization scheduling method comprises the following steps: establishing a fire-wind-light multi-energy complementary frequency dynamic response model, and deducing a time domain analytical expression of a frequency safety index by taking a frequency maximum change rate and a frequency maximum deviation as the frequency safety index of a system; considering a frequency lowest point security constraint and a frequency maximum change rate security constraint, and establishing a frequency security constraint-containing optimization scheduling model; a maximum bearable disturbance power linearization mapping method based on ELM is provided, nonlinear frequency security constraints are converted into a series of linear constraint sets, and linearization processing of the frequency security constraints is achieved; and piecewise linearization is carried out on quadratic terms in the established objective function containing the frequency security constraint optimization scheduling model, the overall model is converted into a mixed integer linear model, and a Gurobi solver is used to solve the linearly processed model. The scheduling method can ensure that the system frequency is safe and stable, so that the established optimal scheduling model gives consideration to both economy and frequency safety.

The invention belongs to the technical field of aircraft structure design, and particularly relates to an aircraft arbitrary space four-bar mechanism parameter optimization method, which comprises the following steps: firstly, determining input conditions or preposed requirements of a four-bar mechanism, secondly, establishing a mathematical model of a four-bar mechanism optimization problem, and then, carrying out derivation and analysis to obtain a four-bar mechanism parameter optimization result. The method comprises the following steps: establishing an analytic expression of parameter variables of the four-bar mechanism and compiling a corresponding calculation program, then carrying out iterative optimization on the parameters of the four-bar mechanism by utilizing an optimization algorithm, and finally outputting an optimization result to obtain an optimal design parameter.

The application discloses a kind of analytic expression probabilistic power flow calculation method, device and storage medium.The application is determined after the mapping relationship between the injection power of node and the operating state of power system by linear power flow model, using first-order polynomial fitting correction method to revise the operating state of power system, make the operating state of the power system after revision be closer to real operating state, according to the mapping relationship between the injection power of node and the operating state of power system after revision, update linear power flow model, carry out analytic expression probabilistic power flow calculation using updated linear power flow model, to effectively improve the accuracy of probabilistic power flow calculation.

The invention relates to a rapid calculation method and system for a space magnetic field of a circular coil. The method comprises the following steps: converting an integral expression of a Biot-Savart law into a closed analytical solution containing first and second types of complete elliptic integrals; on the basis of an analytical expression of elliptic integration, the magnetic field intensity and direction generated by the coil at any point in the space can be efficiently and precisely calculated, and the defects that traditional finite element simulation calculation is complex, consumed time is too long, and efficiency and precision of a direct numerical integration method are difficult to consider at the same time are overcome. The method is particularly suitable for application scenes such as WPT system design, multi-transmitting-coil amplitude-phase control optimization and real-time simulation which need to carry out rapid and frequent calculation on the magnetic field. According to the method, full-automatic rapid calculation from physical parameter input to magnetic field result output is achieved, the core advantage of the method is that an analytical formula is utilized, the tedious process of numerical simulation and integration is avoided, and the calculation efficiency is remarkably improved while high precision is guaranteed.

The application relates to the technical field of estimation of maximum frequency deviation of power grid nodes and is used for solving the problem that the existing estimation method of maximum frequency deviation of power grid nodes cannot reflect the time-space difference characteristics of frequency response of an interconnected power system. First, a two-region frequency response model TAFR considering the time-space difference characteristics of frequency of power grid nodes is established, a time-domain analytical expression describing the node frequency difference mechanism is obtained by adopting the mode analysis and average system frequency model aggregation idea, and a maximum frequency deviation estimation method reflecting the time-space difference is proposed on the basis of the analytical expression. The established TAFR model considers the distribution of tie-line parameters, system inertia and damping, and equivalently processes the new energy frequency modulation parameters, compared with the existing frequency dynamic response model, the model has improved accuracy and applicability while ensuring the time-space difference.

The invention discloses a turbulence model construction method for aerodynamic simulation of a hybrid three-compartment vehicle, and relates to the field of model simulation, and the method comprises the steps: S1, employing a flow field inversion method, and obtaining the dense space distribution of a turbulence model correction term; s2, screening local flow field characteristics related to the correction distribution characteristics; s3, training by adopting a symbol regression method to obtain an analytical expression of a turbulence model correction term; s4, through mathematical derivation, an analytical expression of a correction term is migrated into the standard model, so that the correction term is added in front of a damage term of a transport equation of the dissipation rate of the standard model; s5, selecting an initial value of a to-be-calibrated coefficient in the correction item according to the smoothness degree of the chassis of the hybrid three-compartment vehicle; and S6, calibrating the to-be-calibrated coefficient in the expression in the step S4 on the basis of real vehicle experimental data of the hybrid three-compartment vehicle to obtain the turbulence model for aerodynamic simulation of the hybrid three-compartment vehicle.

The invention discloses a level set field reconstruction method and system based on a three-dimensional finite element grid and a medium. Relates to the technical field of geometric reconstruction. Positioning a space Cartesian grid area intersected with the boundary surface grid based on a separation axis theory; calculating a geometric space point-surface position relation, calculating a directed projection distance from a space Cartesian grid center to a two-dimensional boundary surface grid, and correcting the positive and negative characteristics of a Cartesian grid level set value based on a three-dimensional eight-direction rapid scanning method, so that the grid level set value is consistent with a region boundary value; according to the method, the geometric reconstruction problem of the deformable body which is complex in appearance and free of explicit analytic expression can be efficiently and accurately solved, a coupling approach between a finite element grid and a Cartesian grid is established under a non-conformal grid system, and the flexibility and universality of explicit calculation of a level set function under the Cartesian grid are improved.

The application provides a kind of multi-data stream security analysis method and device with dynamic index.The method comprises: obtaining corresponding dynamic indexable expression chain based on the correlation analysis expression of multi-data stream;Dynamic indexable expression graph is obtained based on the dynamic indexable expression chain;Determine the dynamic expression path in the dynamic indexable expression graph based on the graph node in the dynamic indexable expression graph and the edge between graph nodes;Determine the index corresponding to the flow table based on the dynamic expression path, and perform network security analysis on the multi-data stream based on the index of the flow table.The multi-data stream security analysis method with dynamic index provided by the application can perform dynamic index analysis on multi-data stream correlation analysis expression, improve real-time multi-data stream correlation analysis performance, and improve the efficiency of multi-data stream security analysis and detection.

The invention belongs to the technical field of ship intelligent digital optimization design, and particularly relates to a ship global deformation and optimization method, program and equipment considering drainage volume precision control and deformation rationality and a storage medium. The analytical expression of the drainage volume variation about the global deformation parameter is obtained based on the numerical integration and the recursion formula, the reasonable deformation parameter range of the ship body is determined based on the global deformation modification function and the regular design space, the ship body does not need to be actually deformed, the corresponding volume variation can be directly predicted according to the design parameter, and the calculation accuracy is improved. Accurate feedforward control of volume constraint and direct judgment of rationality of the deformed ship are realized. On the basis, rapid optimization of any ship type structure under the constraint of reasonable deformation and accurate volume control is realized in combination with an agent model. According to the method, a traditional'design-verification-adjustment 'experience iteration process is converted into deterministic one-step calculation, and subsequent optimization can be directly carried out in a pre-verified reasonable design space.

The invention discloses an approximate iteration analytic calculation method for a flight trajectory of a hypersonic flight vehicle, which comprises the following steps of: firstly, establishing a gliding flight dynamics mode of the hypersonic flight vehicle, and ignoring lateral motion to obtain a longitudinal plane dynamics model in a linear matrix differential equation form; secondly, according to the characteristics of different types of flight trajectories, performing targeted simplification on a linear matrix differential equation, and analyzing the characteristics of a jump gliding trajectory and a balanced gliding trajectory through matrix eigenvalues; then, on the basis of the quasi-equilibrium gliding condition, lift-drag ratio characteristics of the equilibrium gliding track are analyzed; then, an approximate dynamic model in a short-time period is established, and a flight path analytical expression in the short-time period is obtained; and finally, based on the flight path analytical calculation expression in the short-time period, proposing a flight path approximate iteration analytical calculation process in the long-time period. The method is simple in structure and easy for engineering realization, and has a relatively large popularization and application space.

The invention discloses an integral transformation fracture phase field model implementation method of a random-order polynomial cohesion law, and belongs to the technical field of phase field model finite element analysis. The method comprises the following steps: inputting a parameter equation of a polynomial cohesion law of any order and a target order into a subprogram; deducing a cohesion eigenfunction of a random-order polynomial cohesion law by using a subprogram to determine an analytical expression of a random-order polynomial energy density eigenfunction, solving a target-order polynomial energy density eigenfunction expression and first two-order derivatives thereof, and determining an expansion expression at a zero point by using Maclaurin series; and inputting the target-order polynomial energy density eigenfunction expression, the first two-order derivatives thereof and the expansion expression at the zero point into a phase field model solver to obtain a calculation result of the target-order polynomial integral transformation fracture phase field model. According to the method, the simulation calculation of the integral transformation fracture phase field model of the target-order polynomial cohesion law can be realized under a universal finite element framework, and the accuracy and the efficiency are improved.

The invention provides a fault rate calculation method for a transformer cooling system considering a leakage fault, which comprises the following steps of: firstly, analyzing a dynamic transfer process among a normal operation state, a fracture fault state and a leakage fault state of a cooler pipeline, and constructing a three-state Markov transfer model of the cooler pipeline; then, sensor information such as main transformer oil temperature, winding temperature, cooler outlet oil pressure, water pressure, oil temperature and water temperature is collected, and a knowledge and data driving-based fracture fault related transfer rate calculation method is provided; thirdly, based on the frequency and the transfer rate which are easy to count or calculate, an analytical expression of the transfer rate and the frequency is established by combining a Markov model and a frequency duration method, and the transfer rate related to the leakage fault in the model is calculated; and finally, integrating the transfer rate, constructing a fault rate calculation model of the transformer cooling system, and proposing to predict the future fault rate of the cooling system based on a time convolutional network and a Transform hybrid model. According to the method provided by the invention, fault rate monitoring, analysis and early warning of the transformer cooling system can be realized.

The invention relates to the technical field of power electronics, and discloses a universal optimization design method of a multi-pulse rectifier direct current side harmonic suppression circuit, which comprises the following steps: S1, establishing a mathematical model of N-step wave input current of a multi-pulse rectifier; s2, based on the structure and rule of the main circuit, extracting a characteristic parameter n, and deducing a general analytical expression of the output current; and S3, mapping the relation between the optimal N step wave input current of the alternating current side and the reference sine wave to the direct current side, and establishing an associated equation of the step wave injection current and the zero-harmonic injection current of the harmonic suppression circuit in combination with the topological law and the ampere-turn balance theorem of the harmonic suppression circuit. According to the method, the optimization standard of the N-step wave input current is constructed, the optimal design of the direct-current side harmonic suppression circuit of any multi-pulse rectifier can be realized by adjusting the turn ratio to meet the standard, and the universality is high.

The invention discloses a deflection calculation method of a steel-concrete composite beam under uniformly distributed loads, and relates to the technical field of high-speed railway bridge construction. The method comprises the following steps: establishing a differential equation of the continuous steel-concrete composite beam and carrying out Laplace transformation to obtain an expression of a non-zero real solution; inverse transformation is carried out on the equation obtained after Laplace transformation, an analytical expression of the deflection of the composite beam is obtained, and the expression comprises a plurality of unknown quantities. Substituting the structure parameters and the material parameters of the composite beam into the non-zero real number solution expression to solve Laplace transformation parameters; the parameters and the boundary conditions are substituted into the deflection analytical expression, and matrix elements used for solving unknown quantities are obtained; and constructing a linear equation set, solving an unknown quantity, substituting the unknown quantity into a deflection analytical expression, and finally obtaining a complete expression for solving the deflection of the continuous composite beam. The complete expression of the deflection of the continuous composite beam provided by the invention is suitable for analyzing the deflection of a single-span or continuous steel-concrete composite beam, and has relatively high calculation precision.

This invention discloses a design method and power supply system for a non-clamped vertical magnetic field induction energy harvesting device. The method includes analyzing the constraints between the output characteristics of the energy harvesting device and its parameters, transmission line current, and connected load parameters using a block-superimposed magnetic flux calculation method, deriving an analytical expression for the maximum output power of the energy harvesting device; employing a dimension-reduced distributed performance parameter optimization method to analyze the maximum output power of the vertical magnetic field induction energy harvesting device when it is located at the optimal energy harvesting position, and designing the energy harvesting device structure; and introducing X-axis and Y-axis rotation matrices and Z-axis and Y-axis offset parameters to analyze the fluctuation of the maximum output power of the energy harvesting device under offset and rotation conditions. This invention improves the installation freedom, scenario applicability, and anti-saturation capability of the magnetic field induction energy harvesting device, and can replace the empirical columnar energy harvesting device in its prototype development. It fills the research gap in the mechanism and performance optimization theory of non-clamped vertical magnetic field induction energy harvesting, and has practical value and guiding significance.

The invention discloses a turbulence model construction method and device, electronic equipment, a storage medium and a vehicle, and relates to the field of sound production control, and the method comprises the steps: obtaining the experiment data of a target vehicle type separation flow, and determining the dense space distribution characteristics of a turbulence model correction term according to the experiment data; determining a first flow field feature related to the dense spatial distribution feature in the flow field features of the separation flow; according to the first flow field feature and the dense space distribution feature, determining an analytical expression of a turbulence model correction term; and constructing a turbulence model for calculating the air resistance of the target vehicle model based on the analytical expression. The precision of the turbulence model for calculating the air resistance of the target vehicle model can be improved, so that the research and development efficiency is improved and the research and development cost is reduced.

The invention provides a new energy networking flexible direct current system fault current sequence component characteristic characterization method, which comprises the following steps of: considering a fault ride-through control strategy of alternating current transmission line side equipment in a new energy networking flexible direct current system, and establishing an equivalent sequence network diagram of the system; solving analytic expressions of fault current sequence components on two sides of the line under different fault types based on the equivalent sequence network diagram; quantizing a fault current sequence component corresponding to each fault type through an analytical expression; and extracting influence parameters of the fault current sequence component, and representing the characteristics of the fault current sequence component in combination with the analytical expression. According to the method, the fault current sequence components under different fault types are systematically quantified, the association rule of the influence parameters and the fault current sequence components is represented, a theoretical basis is provided for formulation, optimization and parameter adjustment of a conventional control strategy, the adaptability and the operation reliability of the control strategy to different fault scenes are remarkably improved, and the method is suitable for large-scale popularization and application. And technical support can be directly provided for safe and stable operation of the new energy network construction flexible direct current system.

The invention provides an uplift pile bearing capacity prediction method and system, and belongs to the field of uplift pile bearing capacity prediction. The technical problem to be solved is that the method comprises the following steps: acquiring parameters including pile parameters, soil body parameters and interface parameters; numerical value calibration and sensitivity identification: establishing a three-dimensional / axisymmetric pile-soil finite element example library, obtaining a load-displacement curve and displacement / plastic zone distribution for mechanism identification and evaluation of the influence of the obtained parameters on the net uplift bearing capacity, and identifying a key control quantity; taking test pile / historical engineering data as a true value, and performing regression setting of model parameters according to an objective function RMSE; analyzing, modeling and solving, including assumed conditions of the model, modeling of a failure surface and expression of an equilibrium equation and bearing capacity, and finally obtaining a net uplift bearing capacity analytic expression; checking and outputting the structure; according to the method, the pile-soil coupling effect and the multi-parameter sensitivity are comprehensively considered, so that the prediction precision and the engineering applicability are improved while the safety is ensured.

The invention discloses a milling chatter stability domain prediction method and system based on an approximate analysis method, and the method comprises the steps: employing a Fourier approximation method to represent a milling steady-state response based on a milling direction coefficient and a time delay characteristic in a milling dynamic model; applying small disturbance to the steady-state response, substituting the steady-state response into the milling dynamic model, and obtaining a milling stability analysis expression by applying a Galerkin process; substituting the bifurcation condition into the stability analysis expression based on the characteristics of the stability domain boundary corresponding to the milling stability bifurcation condition, and deducing an approximate implicit analytical expression of the milling stability boundary; solving the implicit analytical expression by adopting a numerical continuation method to obtain a stability boundary expressed by rotating speed-cutting depth; according to the method, the stability lobe graph is obtained, the problems of initial continuation point selection and transition between lobes with different stability in the solving process are solved, complete and accurate construction of the stability lobe graph is achieved, theoretical guidance of flutter-free working conditions is provided for milling, and the method has important engineering application value.

The invention discloses a multi-scene power system frequency safety index self-adaptive partition reconstruction method, which comprises the following steps: acquiring multi-type frequency modulation resource parameters of a system, and establishing an analytical expression of nonlinear frequency safety index constraints based on a system multi-machine frequency response dynamic model considering multi-type frequency modulation resources; acquiring an original feasible range of system frequency modulation parameters, and sampling to generate an original frequency modulation parameter region and a frequency safety index sample curved surface set in the feasible range of the system frequency modulation parameters; acquiring all sample point subsets on an inertia boundary in the frequency modulation parameter area by traversing sample points meeting a specific equivalent inertia value of the system, acquiring intersection points of frequency safety indexes corresponding to the sample points and constraint limit values of the sample points on an original frequency modulation parameter area boundary, and re-partitioning the original frequency modulation parameter area into a plurality of subareas; and linearizing the linear sub-region which needs to be linearized, outputting a linearization parameter which meets a linear error threshold value requirement, and reconstructing a frequency safety index constraint after partition linearization.

The invention relates to a snowflake-type power distribution network load space transfer optimization method considering N-1 safety and reconstruction. The method comprises the steps of basic data input, load transfer capability analysis and expression, reconstruction time period division, model establishment, conversion, solution and result output. According to the method, the N-1 safety criterion is considered, and the snowflake type power distribution network load transfer capacity is analyzed and expressed based on the topological incidence matrix. Then, taking the load balance degree and the minimum active power loss as targets, and establishing a snowflake type power distribution network multi-period load space transfer model considering N-1 security constraint and switch combination reconstruction; and then, converting the model through a mathematical programming method, solving the model by using commercial solvers such as CPLEX and the like, and outputting a related result. Based on the advantages of a snowflake type power distribution network frame structure, the load space transfer capacity of a feeder line cluster can be effectively utilized through sectional interconnection switch combination reconstruction, so that light and heavy load main transformers and line loads are balanced, large-range flexible transfer of the load space is achieved, and a larger safety transfer margin is provided for power grid operation.

This invention provides an adaptive admittance method and system based on the optimal control LQR theory, including constructing an environmental dynamics model and a robotic arm admittance control model, and deriving their interactive dynamics model; defining state variables and control inputs, and establishing a standard form state-space equation; constructing a quadratic performance index function, transforming admittance parameter optimization into an LQR problem; substituting the linear quadratic regulator framework into the quadratic performance index function, transforming the LQR problem into solving the Riccati matrix differential equation, and deriving analytical expressions for virtual optimal stiffness and virtual optimal damping; within the robotic arm control sampling frequency, inputting real-time collected environmental contact forces and the actual pose of the robotic arm end effector, along with set weighting coefficients, into the analytical expressions to calculate optimal virtual stiffness and damping; and then calculating pose corrections based on the admittance control model to correct the robotic arm's motion trajectory.

The application provides a system and method for generating an index topology graph based on an analytical expression, and relates to the technical field of computer processing, comprising an index topology configuration module, a task scheduling module, an index calculation module and a data display module; the index topology configuration module is used for performing entity recognition analysis on a plurality of user input expression texts respectively to obtain corresponding index topology structures, and sending the index topology structures to the task scheduling module; the task scheduling module is used for creating tasks according to a preset scheduling rule to obtain corresponding index calculation tasks; the index calculation module is used for calculating task priorities according to all the index calculation tasks to obtain index priority calculation data; and the data display module is used for receiving the index priority calculation data and obtaining a target topology graph through the index priority calculation data. The application realizes improved topology graph generation efficiency.

This invention discloses a design method for a Type I ternary coprime array based on coprime sparse arrangement, relating to the field of array structure design. Two coprime positive integers M and N are selected to construct the basic structure of a prototype coprime array. A subarray with n elements and n spacing is added in reverse to the prototype coprime array to form a Type I ternary coprime array. Based on the distribution law of holes in the prototype coprime array, the lower bound expression of the optimal parameters of the subarray can be obtained, ensuring that the holes in the prototype coprime array are completely filled. Analytical expressions for the distinct and continuous degrees of freedom of the Type I ternary coprime array are derived. The method is scalable, allowing for the expansion of array degrees of freedom without destroying the prototype coprime array structure by adding array elements in reverse. The Type I ternary coprime array designed using this method exhibits weak mutual coupling, scalability, and high distinct degrees of freedom.

The application discloses a kind of pure angle target trajectory estimation methods, this method will typical target motion model be converted into the autoregressive moving average model of position component to realize state decoupling, by adjusting the relationship between pseudo-linear least square cost function and process noise energy to obtain regularized pseudo-linear least square model.To reduce the estimation deviation generated by the correlation between pseudo-linear measurement matrix and noise, the regularized instrumental variable least square position estimation is obtained by using instrumental variable estimation method.The method of the application obtains the analytical expression of batch estimation of azimuth position component in the framework of regularized least square estimation.The method has linear time calculation complexity, and its estimation performance is irrelevant to initial state.

This invention discloses a rapid voltage calculation method for a radial virtual power plant based on the branch U-S theorem, belonging to the field of power system analysis and calculation technology. The method includes: acquiring the power network topology of the radial virtual power plant; numbering each node in the network layer by layer; collecting the voltage and power of each node and the impedance data of each branch; calculating the voltage of each node layer by layer based on the branch U-S theorem, according to the known head-end voltage and the power flowing into the tail end; and outputting the voltage calculation results for each node. The branch U-S theorem provides an analytical expression for the tail voltage of a single power branch. This invention avoids the complex iterative process required by traditional voltage calculation methods, enabling rapid and accurate calculation of the voltage of each node in a radial virtual power plant, significantly improving computational efficiency, and is suitable for scenarios such as real-time monitoring, rapid dispatching, and market response of virtual power plants.

This invention belongs to the field of rotating machinery condition monitoring and fault diagnosis technology. It proposes a bearing failure assessment method and system. The method involves acquiring bearing vibration signals and performing mean and absolute value removal processing to obtain the signal to be processed. Next, a lower threshold for the tail is determined, and a set of tail samples with amplitudes greater than this threshold is selected. Then, a soft-truncated power-law distribution model is constructed based on the tail sample set, and the distribution parameters, including the power-law tail exponent, truncation scale parameter, and truncation smoothing exponent, are solved using maximum likelihood estimation. Furthermore, an analytical expression for the improved Gini index is constructed based on the distribution parameters, and the bearing degradation degree index is calculated. Finally, a bearing degradation curve is constructed based on the degradation degree index for different time periods. This invention accurately characterizes the tail impact features of vibration signals using a soft-truncated power-law model and quantifies the degradation degree using an improved Gini index, effectively assessing the health status and failure trend of rolling bearings.

This invention belongs to the field of seismic design technology for railway bridges, and relates to a seismic analysis method for railway bridges using a variable friction coefficient friction pendulum bearing for seismic isolation and vibration reduction. This method is used to determine the design displacement D of the variable friction coefficient friction pendulum bearing; specifically, it includes the following steps: establishing an equivalent period T based on the design parameters of the variable friction coefficient friction pendulum bearing. eff The relationship between the coefficient of friction and the variable friction coefficient; based on the bridge design conditions, the relationship between the equivalent period and the acceleration response spectrum value is established, and S(T) is obtained. eff Determine the damping adjustment coefficient η; equivalent damping ratio ξ; establish support displacement and equivalent period T. eff S(T) eff By considering the relationship between η and η, and combining it with Newton's iteration process, the final support displacement D is determined. This invention saves time by using analytical expressions for calculation, avoiding tedious modeling, long computation times, and convergence issues compared to the finite element model nonlinear time history analysis method. This allows for rapid evaluation of the seismic performance of different support parameters and structural configurations during the scheme design, parameter optimization, and preliminary design stages.

The invention relates to a VSG-based network construction flexible direct current system sending-out line fault current analysis method and device, and the method comprises the steps: firstly constructing an equivalent network model based on a low-pass control strategy of converters at two sides; and then sequence components of output currents of the fan side and the flexible direct current side are calculated based on the equivalent network model, so that a complete analytical expression of fault currents on the two sides is obtained. Compared with an existing method, the method has the advantages that the calculation value can be obtained only through one-time calculation while the calculation precision is guaranteed, analytic calculation of the short-circuit current is achieved, the calculation speed is obviously increased, the defect that analytic calculation of the fault current of the VSG-controlled flexible direct current system transmission line is relatively deficient in existing research is overcome, and the method is suitable for large-scale popularization and application. And a theoretical support is provided for subsequent traditional protection adaptability analysis and criterion optimization.

The application relates to the technical field of bridge engineering and provides a method for determining the spatial dynamic response of a four-corner elastic support bridge, which comprises the following steps: introducing a virtual rigid arm, converting a complex four-corner elastic support bridge system into a simplified mechanical model with clear physical meaning, and obtaining a parameter-adjustable analytical simplified theoretical model; based on the parameter-adjustable analytical simplified theoretical model, solving bridge vertical and torsional dynamic response equations by using Duhamel integral, and outputting dynamic response analytical solutions. The application converts the traditional bridge dynamic response problem depending on high-dimensional finite element simulation into an analytical expression which can be obtained by using Duhamel integral, greatly improving the solution efficiency; the application can be used for bridge design, evaluation and health monitoring: the dynamic response analytical solution of the simplified theoretical model can be embedded into a bridge structure response analysis tool chain, realizing bridge performance prediction and diagnosis under the working conditions of extreme load action, support degradation and the like, and having good engineering popularization value.