This procedure proposes the appropriate equations that determine the rotor movement in order to obtain the aerodynamic actions for each blade in each iteration (the forces that are transmitted to the helicopter and the moments with respect to the articulations), and they are resolved by discretizati
This procedure proposes the appropriate equations that determine the rotor movement in order to obtain the aerodynamic actions for each blade in each iteration (the forces that are transmitted to the helicopter and the moments with respect to the articulations), and they are resolved by discretizations made in the rotor disc, such that that the blade elements of the partition are connected to the disc, which does not rotate, instead of to each blade. The blades pass through an azimuth range in each iteration, making a division of each blade longitudinally into ne elements, and the rotor disc into na sectors, corresponding to divisions of the range of azimuth angles. Therefore, partitions are made on ne·na blade elements, identifying each blade element by means of its distance to the rotor axis and its azimuth angle.
대표청구항▼
1. A real-time performance simulation method for a rotary wing aircraft having at least one rotor including a first rotor, the method comprising: determining geometric data of one blade of the first rotor, on the basis of which movement of the first rotor is determined;determining external forces an
1. A real-time performance simulation method for a rotary wing aircraft having at least one rotor including a first rotor, the method comprising: determining geometric data of one blade of the first rotor, on the basis of which movement of the first rotor is determined;determining external forces and moments that affect each blade of the first rotor in each iteration of numerical real-time integration of equations of motion of the rotary wing aircraft;discretizing the first rotor into virtual blade elements on which calculation to identify dynamic balance of the first rotor is to be performed, wherein the identifying dynamic balance includes:determining, by a processor, the movement of the first rotor to obtain real-time aerodynamic actions that include forces transmitted to the rotary wing aircraft and moments with respect to articulations for each blade in each iteration,wherein the discretization is done on a first rotor disc instead of on each blade itself generating the virtual blade elements on an area covering the first rotor disc and each of the virtual blade elements is connected to the first rotor disc, which does not rotate, instead of to each blade, andwherein the discretization divides the first rotor disc circumferentially into na sectors and radially into ne sectors, which results in na×ne virtual blade elements. 2. The method of real-time performance simulation according to claim 1, wherein each virtual blade element is identified by means of its distance r to an axis of the first rotor and its azimuth angle ψ. 3. The method of real-time performance simulation according to claim 1 or claim 2, wherein the na sectors of the first rotor disc correspond to divisions of the range of azimuth angles. 4. The method according to claim 1 or 2, wherein the geometric data that are determined for the blade are basically torque distribution along the blade θ0(r) and chord distribution along the blade c(r). 5. The method according to claim 1 or 2, wherein the external forces and moments that act on the blade include blade's local angle of attack αr, local Mach number Mr, and each blade's local lift, drag and moment coefficients cl, cd and cm, respectively. 6. The method according to claim 1 or 2, wherein the step of determining the movement further comprises: expressing flapping β and drag ζ as a function of azimuth ψ, using a Fourier series development: f(ψ)=a0+∑i=1n[ai·cos(iψ)+bi·sen(iψ)],minimizing errors associated with the dynamic balance calculation by calculating virtual works with coefficients corresponding to Lagrange multipliers for both cases of the flapping β and the drag ζ, expressing these multipliers as a function of coefficients a0, ai, and bi, andreiterating the calculation process in order to minimize total moments within a predetermined convergence time. 7. The method according to claim 6, wherein the predetermined characteristic convergence time is of at least four times of an integration time step. 8. The method according to claim 1 or 2, wherein the step of determining the movement further comprises calculating, at the end of each iteration, induced velocity for each virtual blade element, which will be used in a next iteration to obtain aerodynamic velocity of each virtual blade element. 9. The method according to claim 8, wherein, to calculate the induced velocity, local induced velocity vi(s) and the induced velocity resulting from aerodynamic force of the virtual blade element vi+(s) are related via a relationship vi(s)=G(s)vi+(s), where G(s) is a transfer function that can be approximated by the following expression as a function of the characteristic time Δt: G(s)=(2Δt)2(s+2ΔT)2. 10. The method according to claim 8, wherein, when there is aerodynamic interference between the first rotor and a second rotor, the step of determining the movement further comprises: including, in the calculation of the aerodynamic velocity of each virtual blade element of each rotor, not only the induced velocity of the first rotor but also the induced velocity of the second rotor, to obtain a reference induced velocity {right arrow over (vip)}, resulting from aerodynamic action of the second rotor at a specific point of the first rotor disc, defined by its coordinates (rip, ψip), andfiltering the reference induced velocity {right arrow over (vip)} such that the reference induced velocity is attenuated in time. 11. The method according to claim 1 or 2, wherein the step of determining the movement further comprises accounting for variation in lift due to ground effect, which is modeled by means of reduction in induced velocity, such that the induced velocity with ground effect (vi)IGE is related to the induced velocity without ground effect (vi)OGE by means of ground effect factor, kG, which is a ratio therebetween. 12. The method according to claim 1 or 2, wherein the step of determining the movement further comprises accounting for variation in lift due to ground effect by means of ground effect factor, kG, as a function of the following parameters: rotor radius R, rotor distance to the ground zG, and distance to the ground at which it is considered that the first rotor is outside the ground effect zG max, where kG kG=1-(R4·zG)21-(R4·zGmax)2. 13. A real-time rotor performance simulator for a rotary wing aircraft having at least one rotor including a first rotor, the simulator comprising: means for determining geometric data of one blade of the first rotor, on the basis of which movement of the first rotor is determined;means for determining external forces and moments that affect each blade of the first rotor in each iteration of numerical real-time integration of equations of motion of the rotary wing aircraft;means for discretizing the first rotor into virtual blade elements on which calculation to identify dynamic balance of the first rotor is to be performed, wherein the identifying dynamic balance includes:means for determining, by a processor, the movement of the first rotor to obtain real-time aerodynamic actions that include forces transmitted to the rotary wing aircraft and moments with respect to articulations for each blade in each iteration,wherein the discretization is done on a first rotor disc instead of on each blade itself generating the virtual blade elements on an area covering the first rotor disc and each of the virtual blade elements is connected to the first rotor disc, which does not rotate, instead of to each blade, andwherein the discretization divides the first rotor disc circumferentially into na sectors and radially into ne sectors, which results in na×ne virtual blade elements.
연구과제 타임라인
LOADING...
LOADING...
LOADING...
LOADING...
LOADING...
이 특허에 인용된 특허 (15)
Kaul, Upender K., Finite-difference simulation and visualization of elastodynamics in time-evolving generalized curvilinear coordinates.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.