¿Existe alguna forma no experimental de obtener CL, CD y CM en función del ángulo de ataque y el ángulo de desviación de la superficie de control?

Estoy trabajando en un simulador de vuelo para pruebas de algoritmos de control. El problema al que me enfrento es encontrar una relación entre CL/CD/CM y el ángulo de ataque ( α ) y el ángulo de desviación de la superficie de control ( d ).

Dado que estoy desarrollando un algoritmo de control para un ala voladora, es crucial saber cómo no solo α pero también d afectar las fuerzas y momentos en mi vehículo de vuelo. Y no puedo darme el lujo de medir las fuerzas y los momentos experimentalmente.

En el simulador Gazebo, el enfoque de LiftDragPlugin es asumir que la curva de elevación se desplaza hacia arriba y hacia abajo por d veces alguna constante. Sin embargo, ese no es el caso, como se puede ver en los resultados numéricos de XFLR5 a continuación.

CL contra ⍺

CD contra ⍺

CM contra ⍺

Leyenda: "Nombre del avión ±δ"

Donde “Nombre del avión ± d ” describe la curva para cualquier d .

Podemos ver que la curva se desplaza en ambas direcciones desde cero- d posición. Por ejemplo, el ángulo de pérdida disminuye a medida que d aumenta

Mi enfoque actual es construir una tabla de búsqueda tridimensional usando las curvas en d y + d como puntos finales e interpolar linealmente las dos curvas a lo largo del eje z. Cuando necesito el parámetro aerodinámico (CL, CD o CM) en algún α 0 y d 0 , solo puedo ubicar el punto en este espacio tridimensional desde α 0 y d 0 .

¿Es esta una forma válida de determinar CL/CD/CM en función de α y d ¿numéricamente? Si no, ¿cómo debo abordar este problema?

¿Dónde ves estancamiento en los datos que presentaste? No veo ningún puesto. El uso de XFLR5 sería computacional (es decir, numérico), no experimental.
Le sugiero que eche un vistazo a JSBsim ( jsbsim.sourceforge.net ) Si bien no es una respuesta directa a su pregunta, hay muchos ejemplos de aviones JSBSim que podrían incluir algunos buenos datos para usted.
No estoy seguro de si XFLR5 puede capturar el estancamiento para un análisis VLM.

Respuestas (3)

Los gráficos en el OP muestran resultados bastante típicos que uno esperaría de los análisis lineales (por ejemplo, VLM).

1. Ascensor

En el rango lineal, el coeficiente de sustentación de la superficie de sustentación ( C L ) se puede expresar como:

C L = a 0 α + a 1 d

dónde a 0 = C L α es la pendiente de la curva de sustentación y a 1 = C L d es la pendiente de elevación por deflexión de aleta plana, que son constantes en el rango lineal y se pueden ver fácilmente en su gráfico de elevación.

Conclusión : puede simplificar sus tablas de búsqueda a las ecuaciones anteriores (todas constantes) si solo está interesado en capturar el efecto lineal (que es todo lo que VLM puede capturar de todos modos).

2. Momento de lanzamiento

En el rango lineal, puede expresar el coeficiente de momento de cabeceo ( C metro ), que supongo que se calcula en un cuarto de cuerda, como:

C metro = C metro a C ( d ) + C L yo a C ( d ) C ¯

dónde C metro a C es el momento de cabeceo en el centro aerodinámico (AC) de la superficie y es una función de d ; yo a C es la distancia entre el 1/4c y AC y también es una función de d ; C ¯ es la longitud de la cuerda de referencia.

El C metro Las parcelas que ha mostrado parecen alinearse muy cerca del centro aerodinámico de la superficie de elevación (por el hecho de que hay muy poco cambio en C metro con respecto a AOA). La deflexión del flap solo cambia significativamente el desplazamiento ( C metro a C ). Las variaciones en las pendientes se deben a que no tiene un flap plano de tramo completo (tiene elevaciones que son tramos parciales), lo que cambia la distribución de sustentación en el tramo y el centro aerodinámico resultante.

Conclusión : puede reemplazar su tabla de búsqueda multidimensional por dos tablas de 1 dimensión para C metro a C y yo a C , y suplántela con la ecuación anterior. Incluso puede reemplazar estas tablas de búsqueda con relaciones lineales constantes si sus gráficos retienen agua.

3. Arrastre

El coeficiente de arrastre ( C D ) es el más interesante. Parece que hay un análisis de capa límite combinado con el VLM, ya que C D 0 cambios por colgajo. Si el número de Reynolds no cambia significativamente, puede intentar la siguiente simplificación:

C D = C D 0 ( d ) + k ( d ) C L 2

donde C_{D_0} es el formulario de arrastre y k es el factor de arrastre inducido, los cuales son una función de d .

Conclusión : puede volver a simplificar sus tablas de búsqueda de múltiples dimensiones a dos tablas de 1 dimensión.

4. Por último...

Definitivamente comenzará a ver efectos no lineales a 30 grados de deflexión de control en la vida real, tal vez incluso a 20 grados. Esto obviamente no es capturado por VLM.

¡Gracias! El ala que analicé en XFLR5 presenta un flap simple de envergadura completa. Todavía no he modelado el ala real y los alerones que quiero simular.

En última instancia, CFD (Dinámica de fluidos computacional), es decir, la simulación por computadora del flujo de aire y sus efectos en el ala, hará el trabajo. Muchos de estos programas están disponibles, incluidos los gratuitos.

Para este tipo de trabajos, VLM es una buena herramienta o existe alguna versión mejorada de xfoil. Creo que CFD será demasiado avanzado. estas en el camino correcto