Transient dynamic response calculation method and system of high-dimensional complex nonlinear rotor system

Abstract

The invention relates to a transient dynamic response calculation method and system for a high-dimensional complex nonlinear rotor system, and the method is characterized in that the method comprises the following steps: 1, obtaining a nonlinear dynamic equation of the rotor system, and carrying out the discretization of the nonlinear dynamic equation; 2, acquiring an algebraic equation which takes displacement response of a rotor system at a to-be-solved moment as a unique unknown quantity by utilizing two basic assumptions of a Newmak method, acquiring an accurate jacobian matrix J of the algebraic equation by utilizing an automatic differential function jacrev in PyTorch, and constructing a Newton-Raphson iterative format to perform iterative calculation; step 3, utilizing two basic assumptions of the Newmak method to further obtain the acceleration and the speed response of the rotor system at the moment to be solved; and step 4, repeating the step 2 to the step 3 to obtain the dynamic response of the rotor system at each moment until the set time is reached, and completing the calculation. According to the method, the problems of low efficiency and difficult solution of obtaining the transient dynamic response of the high-dimensional complex nonlinear rotor system at present are solved, the method can be widely applied to calculation of the transient dynamic response of the high-dimensional complex nonlinear rotor system, and the method has the characteristics of universality and modularization.

Description

TECHNICAL FIELD

[0001] The application belongs to the field of dynamic simulation algorithm, and relates to a transient dynamic response calculation method and system for a high-dimensional complex nonlinear rotor system. BACKGROUND

[0002] Due to complex nonlinearities of various connecting structures such as bearings, gear coupling, squeeze film dampers and other components in the rotor system, and the increasing demand for more detailed and accurate nonlinear transient dynamic response of the rotor system, the established rotor system dynamic model has a high degree of freedom, and the nonlinear expression is mostly complex mathematical forms such as power index and fraction, that is, the rotor system model presents the characteristics of high dimension and complex nonlinearity.

[0003] In the current research on the transient dynamic response of the rotor system, due to the high difficulty and low efficiency of solving the high-dimensional complex nonlinear rotor system, most of the complex nonlinear rotor models are greatly simplified, the rotor system model is low-dimensional and the nonlinear form is relatively simple, the model is seriously distorted, and the error of the transient dynamic response of the rotor system obtained by analysis and the actual system is large. Therefore, there is an urgent need for an effective transient dynamic response calculation method for a high-dimensional complex nonlinear rotor system. SUMMARY

[0004] To solve the technical problems in the prior art, the application provides a transient dynamic response calculation method and system for a high-dimensional complex nonlinear rotor system.

[0005] The transient dynamic response calculation method for a high-dimensional complex nonlinear rotor system provided by the application has the characteristics that it comprises the following steps:

[0006] Step 1, obtaining the nonlinear dynamic equation of the rotor system and discretizing it;

[0007] Step 2, using two basic assumptions of the Newmark method to obtain an algebraic equation with the displacement response of the rotor system at the time to be solved as the only unknown, using the automatic differentiation function jacrev in PyTorch to obtain the accurate Jacobian matrix J of the algebraic equation, and then constructing a Newton-Raphson iteration format for iterative calculation until the convergence condition is met, and outputting the displacement response of the rotor system at the time to be solved;

[0008] Step 3, using two basic assumptions of the Newmark method to further obtain the acceleration and speed response of the rotor system at the time to be solved;

[0009] Step 4, repeating steps 2-3 to obtain the transient dynamic response of the rotor system at each time until the set time is reached, and obtaining the transient dynamic response of the rotor system in the continuous time period.

[0010] Furthermore, the step 1 includes the following steps:

[0011] Step 1.1: By giving the linear mass matrix M, linear damping matrix C, linear stiffness matrix K, generalized nonlinear force matrix of the rotor system And the external excitation Q(t) of the system, the following nonlinear dynamic equation of the rotor system is established:

[0012]

[0013] Step 1.2: Discretize the nonlinear dynamic equations of the rotor system, given the known t n The dynamic response of the rotor system at time Satisfies the following equation:

[0014]

[0015] Step 1.3: Give the time t to be determined n+1 Dynamic response of the rotor system where t n+1 -t n =Δt, satisfying the following equation:

[0016]

[0017] Furthermore, the step 2 includes the following steps:

[0018] Step 2.1: Using the two basic assumptions of the Newmark method, obtain an algebraic equation with the displacement response of the rotor system at the time to be solved as the only unknown;

[0019] Step 2.2: Give the time t to be determined n+1 The unknown displacement response x of the rotor system n+1 The initial value of the iteration x 0 n+1 ;

[0020] Step 2.3: Use automatic differentiation to obtain the unknown displacement response x n+1 Jacobian matrix of the algebraic equation for the iterated values;

[0021] Step 2.4: Use the Jacobian matrix to construct the Newton-Raphson iterative format for iterative calculation and update the unknown displacement response x n+1 The value x of the kth iteration k n+1 , further update the Jacobian matrix;

[0022] Step 2.5: After iterating until the convergence condition is met, output the time t to be determined n+1 The displacement response of the system is x n+1 .

[0023] Furthermore, the step 2.1 includes the following steps:

[0024] Step 2.1.1: Using the two assumptions of the Newmark method, give the time t to be determined n+1 Unknown acceleration and velocity response of the rotor system and time t n Known dynamic response of the system Waiting time t n+1 The displacement response of the system is x n+1 The relationship is as follows:

[0025]

[0026] Where β and δ are parameters given by the Newmark method, β determines stability and accuracy, and δ controls numerical damping;

[0027] Step 2.1.2: Using Newmark’s hypothesis, the time t in step 1.3 is n+1 Dynamic response of the rotor system The equation that satisfies this is written as follows: n+1 The algebraic equation for :

[0028]

[0029] Furthermore, the step 2.3 includes the following steps:

[0030] Step 2.3.1: The system with unknown displacement response x n+1 The algebraic equation expression of the iterative value is compiled into Python language;

[0031] Step 2.3.2: Use automatic differentiation technology to accurately obtain the Jacobian matrix of the Newton-Raphson iterative format. That is, directly call the fully open source automatic differentiation function jacrev provided by PyTorch. According to the input parameter requirements of this function, input the Python language format containing only the unique unknown displacement response x n+1 The algebraic equation expression R of the iterative value can directly obtain the accurate displacement response x containing unknown n+1 The Jacobian matrix J of the algebraic equation of the iterative value is:

[0032]

[0033] Furthermore, the step 2.4 includes the following steps:

[0034] Step 2.4.1: Construct the Newton-Raphson iterative scheme using the Jacobian matrix J of the algebraic equation that changes with the iterative value of the unknown displacement response;

[0035] Step 2.4.2: based on the Newton-Raphson iteration formula, iteration is started from the initial value, when the convergence condition is not met, the Newton-Raphson iteration formula is updated, and the Newton-Raphson iteration is continued.

[0036] Further, the step 3 specifically comprises the following steps:

[0037] Step 3.1: save the to-be-solved time t n+1 The displacement response x of the rotor system n+1 ;

[0038] Step 3.2: based on the two assumptions of the Newmark method, the known displacement response x n The known dynamic response of the rotor system The to-be-solved time t n+1 The known displacement response x of the rotor system n+1 The unknown acceleration and velocity responses of the rotor system at the to-be-solved time t n+1 The unknown acceleration and velocity responses of the rotor system at the to-be-solved time t Can be given by the following formula:

[0039]

[0040] Further, the step 4 specifically comprises the following steps:

[0041] Step 4.1: obtain the to-be-solved time t n+1 The dynamic response of the system And save;

[0042] Step 4.2: take the system time t n+1 The known dynamic response of the rotor system As a new time, starting from the dynamic response of the rotor system at the time t n+1 , repeat steps 2-3 to obtain the unknown dynamic response of the rotor system at the next time, until the calculation time reaches the set time.

[0043] The application also relates to a high-dimensional complex nonlinear rotor system transient dynamic response calculation system, which comprises a computer module, and the computer module applies the high-dimensional complex nonlinear rotor system transient dynamic response calculation method.

[0044] Advantages

[0045] (1) The application solves the problems of low efficiency and difficult solution in obtaining the transient dynamic response of a high-dimensional complex nonlinear rotor system. Based on the Newmark method, the accurate Jacobian matrix in the Newton-Raphson iteration process is obtained by automatic differentiation, so that the transient dynamic response of the high-dimensional complex nonlinear rotor system can be obtained efficiently and accurately.

[0046] (2) The present invention can be widely used in the calculation of transient dynamic responses of high-dimensional complex nonlinear rotor systems. It has the characteristics of versatility and modularity, that is, it is not limited by the complex nonlinear forms of high-dimensional rotor systems, such as rotor systems considering thermal-solid coupling, rotor-casing high-dimensional whole machine systems considering complex nonlinear connection structures, and rotor systems considering complex nonlinear damping. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] Figure 1 1 is a flow chart of the nonlinear dynamic response calculation method of the present invention;

[0048] Figure 2 is a schematic diagram of a model of a complex nonlinear rotor system embodiment of the present invention;

[0049] Figure 3 It is a schematic diagram of the nonlinear dynamic response results of the present invention and other methods for comparison. DETAILED DESCRIPTION

[0050] The following will describe in detail the implementation of the present invention in conjunction with the accompanying drawings and specific implementation cases.

[0051] like Figure 1 As shown, a method for calculating the transient dynamic response of a high-dimensional complex nonlinear rotor system of the present invention includes the following steps:

[0052] Step 1: Obtain the nonlinear dynamic equations of the rotor system and discretize them:

[0053] Step 1.1: By giving the linear mass matrix M, linear damping matrix C, linear stiffness matrix K, generalized nonlinear force matrix of the rotor system And the external excitation Q(t) of the system, the following nonlinear dynamic equation of the rotor system is established:

[0054]

[0055] Step 1.2: Discretize the nonlinear dynamic equations of the rotor system, given the known t n The dynamic response of the rotor system at time Satisfies the following equation:

[0056]

[0057] Step 1.3: Give the time t to be determined n+1 Dynamic response of the rotor system where t n+1 -t n =Δt, satisfying the following equation:

[0058]

[0059] Step 2, using the two basic assumptions of Newmark method, obtain the algebraic equation with the displacement response of the rotor system at the time to be solved as the only unknown, use the automatic differentiation function jacrev in PyTorch to obtain the accurate Jacobian matrix J of the algebraic equation, and then construct the Newton-Raphson iteration format for iterative calculation until the convergence condition is met, output the displacement response of the rotor system at the time to be solved, the specific operation is as follows:

[0060] Step 2.1: Using the two basic assumptions of Newmark method, obtain the algebraic equation with the displacement response of the rotor system at the time to be solved as the only unknown;

[0061] Step 2.1.1: Using the two assumptions of Newmark method, give the time t n+1 Unknown acceleration and velocity response of the rotor system and the time t n Known dynamic response of the rotor system The relationship between the time t n+1 Displacement response x of the rotor system n+1 of the rotor system is as follows:

[0062]

[0063] In the formula, β and δ are parameters given by Newmark method, β determines the stability and accuracy, and δ controls the numerical damping;

[0064] Step 2.1.2: Through the assumption of Newmark, the time t n+1 Dynamic response of the rotor system satisfies the algebraic equation containing only the unknown displacement response x n+1 as shown below:

[0065]

[0066] Step 2.2: Give the initial value x n+1 of the unknown displacement response x n+1 of the rotor system at the time t 0 n+1 ;

[0067] Step 2.3: Use automatic differentiation technology to obtain the Jacobian matrix of the algebraic equation containing the iterative value of the unknown displacement response x n+1 ;

[0068] Step 2.3.1: Compile the algebraic equation expression of the rotor system containing the iterative value of the unknown displacement response x n+1 into Python language;

[0069] Step 2.3.2: Use automatic differentiation technology to accurately obtain the Jacobian matrix of the Newton-Raphson iterative format. That is, directly call the fully open source automatic differentiation function jacrev provided by PyTorch. According to the input parameter requirements of this function, input the Python language format containing only the unique unknown displacement response x n+1 The algebraic equation expression R of the iterative value can directly obtain the accurate displacement response x containing unknown n+1 The Jacobian matrix J of the algebraic equation of the iterative value is:

[0070]

[0071] Step 2.4: Use the Jacobian matrix to construct the Newton-Raphson iterative format for iterative calculation and update the unknown displacement response x n+1 The value x of the kth iteration k n+1 , further update the Jacobian matrix;

[0072] Step 2.4.1: Construct the Newton-Raphson iterative scheme using the Jacobian matrix J of the algebraic equation that changes with the iterative value of the unknown displacement response;

[0073] Step 2.4.2: Based on the Newton-Raphson iterative formula, start the iteration from the initial value. If the convergence condition is not met, update the Newton-Raphson iterative formula and continue the Newton-Raphson iteration.

[0074] Step 2.5: After iterating until the convergence condition is met, output the time t to be determined n+1 The displacement response x of the rotor system n+1 .

[0075] Step 3: Using the two basic assumptions of the Newmark method, further obtain the acceleration and velocity response of the rotor system at the time to be solved. The specific operations are as follows:

[0076] Step 3.1: Save the time t to be solved that meets the convergence conditions n+1 The displacement response of the rotor system is x n+1 ;

[0077] Step 3.2: Using the two assumptions of the Newmark method, take time t n Known dynamic response of the rotor system Waiting time t n+1 The known displacement response x of the rotor system n+1 Based on, the time t n+1 Unknown acceleration and velocity response of the rotor system It can be given by the following formula:

[0078]

[0079] Step 4: Repeat steps 2 to 3 to obtain the transient dynamic response of the rotor system at each moment until the set time is reached, and obtain the transient dynamic response of the rotor system in a continuous time period:

[0080] Step 4.1: Obtain the desired time t of the rotor system n+1 The dynamic response of and save;

[0081] Step 4.2: Set the rotor system at time t n+1 The known dynamic response of For the new moment, from t n+1 Starting from the system dynamic response at the moment, steps 2 to 3 are repeated to obtain the unknown dynamic response of the system at the next moment until the calculation time reaches the set time.

[0082] Example

[0083] Squeeze film dampers are often used in rotor systems to provide vibration reduction, but their stiffness and damping contain nonlinear characteristics, making the rotor system containing a squeeze film damper a complex nonlinear dynamic system with more complex dynamic response. Figure 2 The nonlinear dynamic response of the rotor system with a squeeze film damper is calculated. The nonlinear dynamic equations of the system are as follows:

[0084]

[0085] In the above formula, x, y, θ x ,θ y is the displacement response of the rotor system, is the speed response of the rotor system, is the acceleration response of the rotor system, F x , F y The nonlinear oil film force provided by the squeeze film damper to the rotor structure satisfies the following equation:

[0086]

[0087] in, F r , F t Satisfies the following equation:

[0088]

[0089] in, r=e / C, I satisfies the following equation:

[0090]

[0091] wherein, θ2= θ1+ π. It can be found that F x , F y is a nonlinear support force, and its mathematical form is extremely complex, which is a typical complex nonlinear force.

[0092] The parameters in the above rotor system nonlinear dynamics model are as follows: m = 37.62, J d = 0.8, J P = 1.6, k = 5.4 x 10 6 , c = 265, δ = 6.508 x 10 -4 , C = 2.5 x 10 -4 , R = 3.915 x 10 -2 , L = 0.015, μ = 6.76 x 10 -3 , and w = 1400.

[0093] According to the above step 1.1, the following rotor system nonlinear dynamics equation can be obtained:

[0094]

[0095] wherein, the linear mass matrix M of the system is the linear damping matrix C is the linear stiffness matrix K is the generalized nonlinear force is that is, the nonlinear support force provided by the squeeze film damper, and the external excitation Q(t) suffered by the rotor system is wherein, x is that is, the displacement response of the system, is the first order derivative of x with respect to time t, that is, the velocity response of the system, is the second order derivative of x with respect to time t, that is, the acceleration response of the system.

[0096] According to step 1.2, the velocity response and the displacement response of the rotor system at time t0 (t = 0) are given as then the acceleration response is which can be obtained by the following equation:

[0097]

[0098] According to step 1.3, a time interval is given, which satisfies t1-t0= Δt, then the unknown dynamics response of the rotor system at time t1 to be solved is which satisfies the following equation:

[0099]

[0100] According to step 2.1.1, the unknown acceleration response and the velocity response of the rotor system at the time t1 to be solved are given by using two assumptions of the Newmark method and the known dynamic response of the rotor system at the time t0 The relationship of the unknown displacement response x1 of the rotor system at the time t1 to be solved is as follows:

[0101]

[0102] In the formula, β is 0.25, and δ is 0.5.

[0103] According to step 2.1.2, the obtained unknown acceleration and velocity response of the rotor system at the time t1 to be solved are expressed as the known dynamic response of the system at the time t0 and the relationship of the unknown displacement response x1 of the rotor system at the time t1 to be solved is substituted into the unknown dynamic response of the rotor system at the time t1 to be solved obtained in step 1.3 The equation is satisfied, and after arrangement, the algebraic equation containing only the unknown displacement response x1 is as follows:

[0104]

[0105] According to step 2.3.1, the algebraic equation expression of the value x k 1 of the kth iteration of the system containing the unknown displacement response x1 is compiled into Python language.

[0106] According to step 2.3.2, the fully open-source automatic differentiation function jacrev provided by PyTorch is directly called, and the algebraic equation expression of the value x k 1 of the kth iteration of the unknown displacement response x1 obtained in step 2.3.1 is input, without additional preprocessing, so that the accurate Jacobian matrix J of different iteration values x k 1 can be directly obtained, that is:

[0107]

[0108] According to step 2.4.1, the Newton-Raphson iteration format is constructed by the Jacobian matrix J of the algebraic equation changing with the value x k 1 of the kth iteration of the unknown displacement response x1:

[0109]

[0110] In the above formula, the value x k 1 of the kth iteration of the unknown displacement response x1 can be updated every time the iteration is performed;

[0111] According to step 2.4.2, based on the Newton-Raphson iteration in step 2.4.1, from the initial value x 0 1 starts the iteration. If the convergence condition is not met, it is iterated to obtain the value x of the kth iteration of the unknown displacement response x1. k 1. Update different iteration values ​​x k 1 accurate Jacobian matrix J, further update the Newton-Raphson iteration, perform a new round of Newton-Raphson iteration, and save the displacement response x1 when the output meets the convergence condition and the displacement response x1 of the system at time t1 is to be calculated.

[0112] According to step 3.2, the known dynamic response of the system at time t0 is obtained After finding the known displacement response x1 of the system at time t1, the unknown acceleration and velocity response of the system at time t1 are obtained according to the Newmark hypothesis. As shown below:

[0113]

[0114] According to step 4.1, the known displacement response x1 of the system at time t1 obtained by iteration, the known acceleration and velocity response of the system at time t1 obtained in step 3.2 Arrange into known form and save it.

[0115] According to step 4.2, the dynamic response of the system The known time t1 is set as the new time, and the system dynamic response at time t1 is known. Start from step 2 and repeat step 3 to obtain the unknown dynamic response of the system at time t2 and save it until the calculation time reaches the set time, completing the calculation of the nonlinear dynamic response of the system.

[0116] The complex nonlinear rotor system θ obtained by solving the case according to the present invention y The results are plotted on Figure 3 The solution results are consistent with those of other methods. The present invention can accurately solve the complex nonlinear dynamic equations of the rotor system.

[0117] The above contents of the present invention are only preferred embodiments of the present invention and are not intended to limit the implementation scheme of the present invention. Ordinary technicians in this field can easily make corresponding changes or modifications based on the main concepts and spirit of the present invention. Therefore, the scope of protection of the present invention shall be based on the scope of protection required by the claims.

Claims

1. A method for calculating the transient dynamic response of a high-dimensional complex nonlinear rotor system, characterized in that: The following steps are involved: Step 1: Obtain the nonlinear dynamic equations of the high-dimensional complex nonlinear rotor system and discretize them; Step 2: Using the two basic assumptions of the Newmark method, an algebraic equation with the displacement response of the rotor system at the time to be solved as the only unknown is obtained. The automatic differentiation function jacrev in PyTorch is used to obtain the accurate Jacobian matrix J of the algebraic equation. Then, a Newton-Raphson iterative scheme is constructed to perform iterative calculations until the convergence conditions are met, and the displacement response of the rotor system at the time to be solved is output. Step 3: Using the two basic assumptions of the Newmark method, further obtain the acceleration and velocity response of the rotor system at the time to be solved; Step 4: Repeat steps 2 to 3 to obtain the transient dynamic response of the rotor system at each moment until the set time is reached, and obtain the transient dynamic response of the rotor system in a continuous time period.

2. The method for calculating the transient dynamic response of a high-dimensional complex nonlinear rotor system according to claim 1, characterized in that: The step 1 comprises the following steps: Step 1.1: By giving the linear mass matrix M, linear damping matrix C, linear stiffness matrix K, generalized nonlinear force matrix of the rotor system And the external excitation Q(t) of the system, the following nonlinear dynamic equation of the rotor system is established: Step 1.2: Discretize the nonlinear dynamic equations of the rotor system, given the known t n The dynamic response of the rotor system at time Satisfies the following equation: Step 1.3: Give the time t to be determined n+1 Dynamic response of the rotor system where t n+1 -t n =Δt, satisfying the following equation:

3. The method for calculating the transient dynamic response of a high-dimensional complex nonlinear rotor system according to claim 1, characterized in that: The step 2 comprises the following steps: Step 2.1: Using the two basic assumptions of the Newmark method, obtain an algebraic equation with the displacement response of the rotor system at the time to be solved as the only unknown; Step 2.2: Give the time t to be determined n+1 The unknown displacement response x of the rotor system n+1 The initial value of the iteration x 0 n+1 ; Step 2.3: Use automatic differentiation to obtain the unknown displacement response x n+1 Jacobian matrix of the algebraic equation for the iterative value; Step 2.4: Use the Jacobian matrix to construct the Newton-Raphson iterative format for iterative calculation and update the unknown displacement response x n+1 The value x of the kth iteration k n+1 , further update the Jacobian matrix; Step 2.5: After iterating until the convergence condition is met, output the time t to be determined n+1 The displacement response of the rotor system is x n+1 .

4. The method for calculating the transient dynamic response of a high-dimensional complex nonlinear rotor system according to claim 3, characterized in that: The step 2.1 includes the following steps: Step 2.1.1: Using the two assumptions of the Newmark method, give the time t to be determined n+1 Unknown acceleration and velocity response of the rotor system and time t n Known dynamic response of the rotor system Waiting time t n+1 The displacement response of the rotor system is x n+1 The relationship is as follows: Where β and δ are parameters given by the Newmark method, β determines stability and accuracy, and δ controls numerical damping; Step 2.1.2: Using Newmark’s hypothesis, the time t in step 1.3 is n+1 Dynamic response of the rotor system The equation that satisfies this is written as follows: n+1 The algebraic equation for :

5. The method for calculating the transient dynamic response of a high-dimensional complex nonlinear rotor system according to claim 3, characterized in that: The step 2.3 includes the following steps: Step 2.3.1: The system with the unknown displacement response x n+1 The algebraic equation expression of the iterative value is compiled into Python language; Step 2.3.2: Use automatic differentiation technology to accurately obtain the Jacobian matrix of the Newton-Raphson iterative format. That is, directly call the fully open source automatic differentiation function jacrev provided by PyTorch. According to the input parameter requirements of the function, input the displacement response x in Python language containing only the unique unknown n+1 The algebraic equation expression R of the iterative value can directly obtain the accurate displacement response x containing unknowns. n+1 The Jacobian matrix J of the algebraic equation of the k-th iteration value is:

6. The method for calculating the transient dynamic response of a high-dimensional complex nonlinear rotor system according to claim 3, characterized in that: The step 2.4 includes the following steps: Step 2.4.1: Construct the Newton-Raphson iterative scheme using the Jacobian matrix J of the algebraic equation that varies with the iterative value of the unknown displacement response; Step 2.4.2: Based on the Newton-Raphson iterative formula, start the iteration from the initial value. If the convergence condition is not met, update the Newton-Raphson iterative formula and continue the Newton-Raphson iteration.

7. The method for calculating the transient dynamic response of a high-dimensional complex nonlinear rotor system according to claim 1, characterized in that: The step 3 specifically includes the following steps: Step 3.1: Save the time t to be solved that meets the convergence conditions n+1 The displacement response of the rotor system is x n+1 ; Step 3.2: Using the two assumptions of the Newmark method, take time t n Known dynamic response of the rotor system Waiting time t n+1 The known displacement response x of the rotor system n+1 Based on, the time t n+1 Unknown acceleration and velocity response of the rotor system It can be given by the following formula:

8. The method for calculating the transient dynamic response of a high-dimensional complex nonlinear rotor system according to claim 1, characterized in that: The step 4 specifically includes the following steps: Step 4.1: Get the time t to be found n+1 Dynamic response of the rotor system and save; Step 4.2: Set time t n+1 Known dynamic response of the rotor system For the new moment, from t n+1 Starting from the dynamic response of the rotor system at the moment, steps 2 to 3 are repeated to obtain the unknown dynamic response of the rotor system at the next moment until the calculation time reaches the set time.

9. A transient dynamic response calculation system for a high-dimensional complex nonlinear rotor system, characterized in that: The system includes a computer module, which applies the transient dynamic response calculation method of the high-dimensional complex nonlinear rotor system according to any one of claims 1 to 8.

Images