Nomenclature A s= cross-sectional area of suction line, m 2 A th = nozzle throat area, m 2 C f = thrust coefficient, dimensionless g = standard gravitational acceleration, m=s 2 h s = height of propellant in suction line, m h t = height of propellant in tank, m I s = inertance of propellant in suction line, s 2 =m 2 K = effective stiffness of accumulator which is defined as K a K p =K a K p , 1=m 2 K a = stiffness of the accumulator, 1=m 2 K p = stiffness of pump cavitation, 1=m 2 k = the ratio between effective stiffness and stiffness of accumulator defined as K =K a , dimensionless L = length, m M s = modal mass of structure mode, kg i, j, n = index, dimensionless m 1 = pump gain, dimensionless N = number of thrust chambers, dimensionless p = oscillatory pressure, Pa R c = linearalized combustion resistance for propulsion flow perturbation, s=m 2 w = weight displacement of flow, N = mass flow gain factor of a pump, s = the propulsion modal damping coefficient, dimensionless = dimensionless mass flow gain factor of a pump defined as , dimensionless = couple strength between the structure mode and the propulsion system, dimensionless = frequency ratio between the propulsion mode and the structure mode, dimensionless = the structure damping ratio, dimensionless = density of the propellant, kg=m 3 ' = structural mode gain, dimensionless = circular natural frequency of the structure mode, 1=s ! = circular natural frequency of the propulsion mode, 1=sSubscripts c = combustion chamber p = pump s = suction line t = tank
A method for deriving a non-singular state-space formulation based on Rubin's model for Pogo analysis is presented in this study. Because of the non-singularity, the state-space model can be directly used in frequency-domain analyses and time-domain simulations. To describe the assembly method concisely, the dynamic equations of the nine types of independent elements are described in a standard manner. The nine types of elements are divided into two classes according to characteristics of the dynamic equations. The mapping relationship between the local and global numbers of elements and nodes is obtained by numbering all of the elements and nodes. By integrating the element stiffness matrixes to obtain the total stiffness matrix used in the finite element method, the coefficient matrixes of the improved Rubin's model are assembled from the coefficient matrixes of all of the elements according to the mapping relationship. Based on the non-singular model, the time-varying simulation with the nonlinear property of the accumulator can be done conveniently by revising the state-space model. The successful application of this method to a Pogo analysis of a certain type of CZ rocket used in China verifies the correctness and efficiency of the method of this study.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
customersupport@researchsolutions.com
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
This site is protected by reCAPTCHA and the Google Privacy Policy and Terms of Service apply.
Copyright © 2025 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.