Method of simulation of unsteady aerodynamic loads on an external aircraft structure
원문보기
IPC분류정보
국가/구분
United States(US) Patent
등록
국제특허분류(IPC7판)
G01M-009/06
G06F-017/50
출원번호
US-0023981
(2013-09-11)
등록번호
US-9767229
(2017-09-19)
우선권정보
FR-12 58516 (2012-09-11)
발명자
/ 주소
Calderon, Raul
David, Christophe
출원인 / 주소
Airbus Operations (S.A.S.)
대리인 / 주소
Nixon & Vanderhye P.C.
인용정보
피인용 횟수 :
0인용 특허 :
10
초록▼
The invention relates to a method for simulating the unsteady aerodynamic loads being exerted on the external structure of an aircraft, notably in the context of a simulation of the buffeting of a wing surface in an airflow. The method includes a step of measurement of pressure (410) at different po
The invention relates to a method for simulating the unsteady aerodynamic loads being exerted on the external structure of an aircraft, notably in the context of a simulation of the buffeting of a wing surface in an airflow. The method includes a step of measurement of pressure (410) at different points of a grid, a step of calculation of the spectral density at these points (420) followed by extrapolation/interpolation operations to calculate the missing measurements (430), a step of estimation of the pressure coherence (440) for each pair of points of the grid, a step of estimation of the interspectral pressure density (450) for these same pairs of points, from the coherence thus estimated, and a step of calculation of the aerodynamic loads (460) by summing the real part of the interspectral density for the area of the wing surface having a separation of the boundary layer.
대표청구항▼
1. A method of simulation by computer of unsteady aerodynamic loads being exerted on an external structure of an aircraft, positioned in an airflow, the method comprising: (a) a step of measurement of pressure by multiple sensors located on the surface of the external structure, and positioned at fi
1. A method of simulation by computer of unsteady aerodynamic loads being exerted on an external structure of an aircraft, positioned in an airflow, the method comprising: (a) a step of measurement of pressure by multiple sensors located on the surface of the external structure, and positioned at first plurality of points of a grid having a first direction in a direction of the airflow, and a second direction, separate from the first direction, where the pressure signals of the sensors are stored in a memory;(b) a step of calculation of a spectral density of the pressure signals at the first plurality of points of the grid as a spectral pressure density at the first plurality of points;(c) a step of extrapolation and/or interpolation of the spectral pressure density at the first plurality of points in order to obtain a spectral pressure density at second plurality of points of the grid, the second plurality of points having no sensor;(d) a step of estimation of coherence of pressure for multiple pairs of points of the grid, from a predetermined model of a coherence function;(e) a step of calculation of an interspectral pressure density for the multiple pairs of points of the grid, from the coherence of the pressure estimated for these pairs of points, and a spectral pressure density being calculated for the points of these pairs; and(f) a step of calculation of the unsteady aerodynamic loads by summing the real part of the interspectral density in the area of the external structure having a separation of the boundary layer of the airflow. 2. The method of simulation according to claim 1, wherein in step (c) the spectral pressure density of the first plurality of points is extrapolated in linear fashion in the second direction and the spectral pressure density of the first plurality of points is interpolated using an interpolation polynomial of degree greater than or equal to three in the first direction. 3. The method of simulation according to claim 1, wherein in step (d) the model of the coherence function γnm(f), for a pair (n,m) of points of the grid, is given by: γnm(f)=exp(-2πf(xm-xnαx(f)+ym-ynαy(f)))·exp(jφ(f,xm-xn,ym-yn))φ(f,xm-xn,ym-yn)=2πf(xm-xnVx+ym-ynVy) where xn, yn are the respective coordinates of the first point of the pair of points in the first and second directions; xn, yn are the respective coordinates of the second point of the pair of points in the first and second directions; Vx is a speed of propagation of instability in the first direction; Vy is a speed of propagation of the instability in the second direction, f is a frequency; αx (f) is a decoherence factor in the first direction and αy (f) is a decoherence factor in the second direction. 4. The method of simulation according to claim 3, wherein an estimation {circumflex over (α)}x(f) of the decoherence factor in the first direction is calculated by the following average: α^x(f)=-2πfΔxK∑k=1K1ln(γkΔx(f)) where γkΔx(f) is coherence at frequency f between pressure signals measured respectively by two sensors separated by a distance Δx in the first direction, and in that an estimation {circumflex over (α)}y(f) of the decoherence factor in the second direction is calculated by the following average: α^y(f)=-2πfΔyK∑k=1K1ln(γkΔy(f)) where γkΔy(f) is coherence, at frequency f, of the pressure signals measured respectively by two sensors separated by a distance Δy in the second direction, and where K is the number of measurements considered for the average. 5. The method of simulation according to claim 3, wherein an estimation {circumflex over (V)}, of the speed of propagation of the instability in the first direction is calculated from the following average: V^x=2πCard(Ωx)∑n,m∈ΩxΔxnmλnmwhere λnm is the gradient of a coherence phase, as a function of frequency, of pressure signals measured respectively by two sensors n,m separated by a relative distance Δxnm=xm−xn in the first direction, Ωx is a first set of pairs of sensors, and Card (Ωx) is a cardinal of Ωx;and in that an estimation {circumflex over (V)}y of the speed of propagation of the instability in the second direction is calculated from an average of the following gradient: V^y=2πCard(Ωy)∑n,m∈ΩyΔynmμnm Where μnm is the gradient of a coherence phase, as a function of frequency, of the pressure signals measured respectively by two sensors n,m separated by a relative distance Δynm=ym−yn in the second direction, and Ωy is a second set of pairs of sensors. 6. The method of simulation according to claim 5, wherein the speed of propagation of the instability in the first direction and/or in the second direction is obtained for pairs of sensors the measured pressure signals of these pairs of sensors having a coherence module greater than a predetermined threshold (χTh). 7. The method of simulation according to claim 1, wherein in step (e) an interspectral pressure density Snm(f) for a pair of points of the grid is obtained by: Snm(f)=γnm(f)√{square root over (Snm(f))}√{square root over (Snm(f))} where γnm(f) is a pressure coherence between the first and second plurality of points of the pair of points; Snm(f) is the spectral pressure density at the first point of the pair of points and Snm(f) is the spectral pressure density at the second point of the pair of points. 8. The method of simulation according to claim 1, wherein in step (f) the aerodynamic loads on the external structure are obtained from: Fexcitation(f)=∑n,m∈ΩRe(Snm(f))σ→nσ→m where Snm(f) is an interspectral pressure density between a first point n and a second point m of the grid, f is a frequency, Re(.) designates the real part, where the surface of the external structure is divided into elementary surfaces associated with the different points of the grid, {right arrow over (σ)}n is the vector with norm equal to the elementary surface associated with a first point n and orthogonal to the external structure at this point, {right arrow over (σ)}m is the vector with a norm equal to the elementary surface associated with second point m and orthogonal to the external structure at this point, where the summation of the interspectral density is accomplished for a set Ω of pairs of points of the grid belonging to the area of the external structure having a separation of the boundary layer of said flow. 9. A method of simulation by computer of the buffeting of an external structure of an aircraft in an airflow, wherein a simulation is made of unsteady aerodynamic loads, excitation Fexcitation(f), being exerted on this external structure according to the simulation method of claim 1, and in that the equation is resolved to the following vibrations: M{umlaut over (z)}+C′ż+K′z=Fexcitation where M is the mass matrix of the external structure, C′ and K′ are respectively the damping and rigidity matrices of the external structure in the wind, and z is the coordinate in the direction orthogonal to the plane of the external structure.
연구과제 타임라인
LOADING...
LOADING...
LOADING...
LOADING...
LOADING...
이 특허에 인용된 특허 (10)
McIntyre, Melville Duncan Walter; Houck, Andrew William; Bridgewater, Russell Tanner; Freeman, Robert E.; Salo, Paul; Wilson, Douglas L.; Moore, Jonathan K., Alternative method to determine the air mass state of an aircraft and to validate and augment the primary method.
Holmes Bruce J. (Newport News VA) Carraway Debra L. (Norfolk VA) Holmes Harlan K. (Newport News VA) Moore Thomas C. (Poquoson VA), Crossflow vorticity sensor.
Velázquez López, Angel Gerardo; Alonso Fernández, Diego; Vega de Prada, José Manuel; Lorente Manzanares, Luis Santiago; De Pablo Fouce, Valentin, Method and system for a quick calculation of aerodynamic forces on an aircraft in transonic conditions.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.