Nonlinear power flow feedback control for improved stability and performance of airfoil sections
원문보기
IPC분류정보
국가/구분
United States(US) Patent
등록
국제특허분류(IPC7판)
G06F-017/50
G06F-007/60
G06G-007/48
G06G-007/50
B63H-003/00
출원번호
US-0633045
(2009-12-08)
등록번호
US-8527247
(2013-09-03)
발명자
/ 주소
Wilson, David G.
Robinett, III, Rush D.
출원인 / 주소
Sandia Corporation
대리인 / 주소
Peacock Myers P.C.
인용정보
피인용 횟수 :
3인용 특허 :
4
초록▼
A computer-implemented method of determining the pitch stability of an airfoil system, comprising using a computer to numerically integrate a differential equation of motion that includes terms describing PID controller action. In one model, the differential equation characterizes the time-dependent
A computer-implemented method of determining the pitch stability of an airfoil system, comprising using a computer to numerically integrate a differential equation of motion that includes terms describing PID controller action. In one model, the differential equation characterizes the time-dependent response of the airfoil's pitch angle, α. The computer model calculates limit-cycles of the model, which represent the stability boundaries of the airfoil system. Once the stability boundary is known, feedback control can be implemented, by using, for example, a PID controller to control a feedback actuator. The method allows the PID controller gain constants, KI, Kp, and Kd, to be optimized. This permits operation closer to the stability boundaries, while preventing the physical apparatus from unintentionally crossing the stability boundaries. Operating closer to the stability boundaries permits greater power efficiencies to be extracted from the airfoil system.
대표청구항▼
1. A method of determining pitch stability characteristics of an airfoil, comprising using a computer to numerically integrate a differential equation of motion characterizing a time-dependent response of an airfoil's pitch angle, α, measured relative to direction of an incoming fluid flowing toward
1. A method of determining pitch stability characteristics of an airfoil, comprising using a computer to numerically integrate a differential equation of motion characterizing a time-dependent response of an airfoil's pitch angle, α, measured relative to direction of an incoming fluid flowing towards the airfoil; wherein the method comprises steps of: a) inputting airfoil design parameters and system initial conditions to a memory unit of a computer;b) programming into the memory unit a differential equation of motion that models a response of the airfoil's pitch angle as a function of time;c) the computer calculating the time-dependent response of the airfoil's pitch angle, energy generation (Egen), and energy dissipation (Ediss), by using a microprocessor in the computer to numerically integrate the differential equation of motion over a specified period of time;d) storing results of the numerical integration step in the computer's memory, wherein the results comprise: α, Egen, and Ediss as a function of time;e) inspecting the results, and then categorizing an airfoil's transient response as being: (1) Passive if the response is damped and dissipative, (2) Neutral Limit-Cycle if the response is balanced and stable, or (3) Unstable if the response is generative and divergent. 2. The method of claim 1, further comprising: f) after step e), changing one or more airfoil design parameters, or initial conditions, or both;g) then repeating another analysis run comprising steps a) through e); andh) iterating steps a) through g) until a Neutral Limit-Cycle response has been achieved. 3. The method of claim 2, wherein if the airfoil's transient response in step e) is Passive; then step f) comprises changing a design of the airfoil in such a way that increases Egen, or decreases Ediss, or that does both, in the next analysis run, when compared to a previous analysis run. 4. The method of claim 2, wherein if the airfoil's transient response in step e) is Unstable; then step f) comprises changing a design of the airfoil in a way that decreases Egen, or increases Ediss, or that does both, in the next analysis run, when compared to a previous analysis run. 5. The method of claim 1, wherein the differential equation of motion comprises a single independent variable, α; and the equation comprises: I{umlaut over (α)}+C{dot over (α)}+CNL sin ({dot over (α)})+Kα+KNLα3=−Mα(α)+M{dot over (α)}({dot over (α)},α)wherein α is the pitch angle, K is torsional stiffness, KNL is nonlinear torsional stiffness, C is torsional damping, CNL is nonlinear torsional damping, I is an airfoil's rotational moment of inertia, and Mα, M{dot over (α)} are applied aerodynamic moments. 6. The method of claim 1, further comprising: using the computer to calculate the energy terms Egen and Ediss, where: Egen=∮τ[Mα(α)+Mα.(α.,α)]α.ⅆtEdiss=∮τ[C(α.)+CNLsin(α.)]α.ⅆtwherein C is torsional damping, CNL is nonlinear torsional damping, and Mα, M{dot over (α)} are applied aerodynamic moments;comparing Egen to Ediss; andthen deciding, if Egen is equal to Ediss, that a Neutral Limit-Cycle response has been produced. 7. The method of claim 1, wherein the numerical integration step uses a fixed time-step, 4th order Runge-Kutta solution technique. 8. The method of claim 1, wherein the torsional stiffness, K, is selected to represent a torsional modulus of elasticity of a long blade or wing having said airfoil shape. 9. A method of determining pitch stability of an airfoil, comprising using a computer to numerically integrate a differential equation of motion that includes terms describing PID controller action; wherein the differential equation characterizes a time-dependent response of an airfoil's pitch angle, α, measured relative to direction of an incoming fluid flowing towards the airfoil; wherein the method comprises steps of a) inputting airfoil design parameters, system initial conditions, and feedback control reference values to a memory unit of a computer;b) programming into the memory unit a differential equation of motion that has terms describing PID controller action, which models a response of the airfoil's pitch angle as a function of time;c) the computer calculating the time-dependent response of the airfoil's pitch angle, energy generation (Egen), and energy dissipation (Ediss), by using a microprocessor in the computer to numerically integrate the differential equation of motion over a specified period of time;d) storing results of the numerical integration step in the computer's memory, wherein the results comprise: α, Egen, and Ediss as a function of time;e) inspecting the results, and then categorizing an airfoil's transient response as being: (1) Passive if the response is damped and dissipative, (2) Neutral Limit-Cycle if the response is balanced and stable, or (3) Unstable if the response is generative and divergent. 10. The method of claim 9, further comprising: f) after step e), adjusting one or more PID controller gain values, integral gain (KI), proportional gain (KP), or derivative gain (KD);g) then repeating another analysis run comprising steps a) through e), and then:h) iterating steps a) through g) until a Neutral Limit-Cycle response has been achieved. 11. The method of claim 10, wherein if the airfoil's transient response in step e) is Passive; then step f) comprises changing KI or KD in such a way that increases Egen, or decreases Ediss, or that does both, in the next analysis run, when compared to a previous analysis run. 12. The method of claim 10, wherein if the airfoil's transient response in step e) is Unstable; then step f) comprises changing KI or KD in such a way that decreases Egen, or increases Ediss, or that does both, in the next analysis run, when compared to a previous analysis run. 13. The method of claim 9, wherein the differential equation of motion comprises a single independent variable, α; and the equation comprises: I{umlaut over (α)}+[K+KP]α+KNLα3=−[C+KD]{dot over (α)}−CNL sin ({dot over (α)})+Mα(α)+M{dot over (α)}({dot over (α)},α)−KI∫0tαdτwherein α is the pitch angle, K is torsional stiffness, KNL is nonlinear torsional stiffness, C is torsional damping, CNL is nonlinear torsional damping, I is airfoil's rotational moment of inertia, and Mα, M{dot over (α)} are applied aerodynamic moments. 14. The method of claim 9, further comprising: using the computer to calculate the energy terms Egen and Ediss, where: Egen=∮τ[Mα(α)+Mα.(α.,α)-KI∫0tαⅆτ]α.ⅆtEdiss=∮τ[(C+KD)α.+CNLsin(α.)]α.ⅆtwherein KI is integral gain, KD is derivative gain, C is torsional damping, CNL is nonlinear torsional damping, and Mα, M{dot over (α)} are applied aerodynamic moments;comparing Egen to Ediss; andthen deciding, if Egen is equal to Ediss, that a Neutral Limit-Cycle response has been produced. 15. The method of claim 14, further comprising performing steps a) through k) in real-time during operation of a system comprising said airfoil; and then providing updated values of the PID controller gain constants, KI, KP, and KD to a PID controller unit that supplies signals in real-time to an actuator coupled to the airfoil, that generates opposing feedback control forces, μ, on the airfoil according to the following equation: I{umlaut over (α)}+Kα+KNLα3=−C{dot over (α)}−CNL sin ({dot over (α)})+μ+Mα(α)+M{dot over (α)}({dot over (α)},α)whereμ=−KPα−KI∫0tαdτ−KD{dot over (α)}. 16. The method of claim 9, wherein the numerical integration step uses a fixed time-step, 4th order Runge-Kutta solution technique. 17. The method of claim 9, wherein the torsional stiffness, K, is selected to represent the torsional modulus of elasticity of a long blade or wing having said airfoil shape.
연구과제 타임라인
LOADING...
LOADING...
LOADING...
LOADING...
LOADING...
이 특허에 인용된 특허 (4)
Steven Y. Goldsmith, Alpha-beta coordination method for collective search.
Raymond H. Byrne ; John J. Harrington ; Steven E. Eskridge ; John E. Hurtado, Cooperative system and method using mobile robots for testing a cooperative search controller.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.