21 agosto 2011

Empezando a analizar resultados

Bueno, pues más o menos el modelo está encauzado y con una linea clara de desarrollo. Con él ya se puede simular el comportamiento de una moto en cualquiera de las tres situaciones fundamentales de estudio, que son, la aceleración, la frenada y el paso por curva.

Como ya he comentado en anteriores entradas, en el paso por curva existen ciertos mecanismos de comportamiento del neumático difíciles de modelar sobre todo cuando no se dispone de los datos necesarios del mismo o no se tienen medios para ensayar. El modelo que quería implementar desde un principio es el más usado en el mundo de la competición que es la llamada "Fórmula Mágica" o modelo de Placejka.
Finalmente, teniendo en cuenta los datos de entrada utilizados en mi modelo, no ha sido necesario utilizar esto, simplemente he hecho un equilibrio de fuerzas en todo momento sobre cada neumático para que fueran capaces de "copiar" las aceleraciones medidas con la telemetría. En mi opinión esto funciona perfectamente según he considerado mi modelo pero es una lástima no disponer de estos datos para plantar el modelo de otra forma y conocer el deslizamiento y la deriva de los mismos.

A continuación muestro un esquema básico del modelo completo, para que de un simple vistazo se entienda su mecánica. Si lo escribo en texto como he hecho en anteriores ocasiones, me da la sensación de que aburre y la gente no lo lee o quizás quede menos claro. Una imagen vale más que mil palabras:


He empezado a lanzar un primer test del modelo que simula el comportamiento de la Ducati 1098 en 35 segundos.

Los 5 primeros segundos comienzan con la moto colgada, con las suspensiones completamente extendidas. Posteriormente aplico la aceleración de la gravedad de una forma más suave y no comenzando en el instante cero, ya que al hacer esto me da problemas de convergencia, así que decidí darle un poco de tiempo para que resolviera las primeras ecuaciones de equilibrio para estabilizar y luego aplico la gravedad.
En el segundo 5, empezamos en condiciones de marcha que se dan en la recta de meta del circuito de Alcañiz y las siguientes tres curvas que la siguen, izquierda, derecha y derecha. Me gustaría simular una vuelta entera, pero supone bastante tiempo de cálculo y prefiero ir piano piano.

Teniendo en cuenta las hipotesis de cálculo obligadas por no disponer de datos, los medios y experiencia en este campo, me sorprende muy gratamente los resultados obtenidos puesto que aunque existen ligeras diferencias entre los resultados y las medidas de la telemetría, creo que tengo aun algo de margen para mejorar y creo que con un poco de trabajo y mejores medios podría llegar a clavar los resultados.

A continuación muestro las entradas introducidas al modelo:


Estos son los datos de entrada necesarios para como he modelado mi sistema, pero existe la posibilidad de dar otros como características del circuito, relación de marchas del cambio, grado de apertura del gas.... todo en función de como se quiera estudiar.

Seguidamente se puede apreciar una comparativa entre los resultados obtenidos matemáticamente y los medidos realmente con la telemetría en esta zona del circuito de Alcañiz. La primera muestra la compresión de la horquilla y la segunda la distensión en el colín (compresión del amortiguador).


Resulta muy gratificante ver estos resultados, en los que practicamente sin ajustes, haya llegado a unos resultados tan buenos. Ahora es cuando pienso que merece la pena tanto esfuerzo y tiempo de trabajo, el placer de llegar a esto.

Sigo mostrando más resultados. Estos son geométricos, y corresponden a como varia el avance y la batalla en estas condiciones. Finalmente unos datos, solamente orientativos del par necesario a aplicar sobre los semimanillares para salir de la condición de equilibrio. Es un cálculo sencillo función del coeficiente de rozamiento, carga normal sobre la rueda delantera y el avance medido anteriormente. Esto da una idea de la maniobrabilidad de la moto al variar el avance. Digo que es orientativo porque no he entrado en detalle de otros factores como el par autoalineante del neumático, efecto giroscópico... así que lo muestro solamente por interés, no hay que tener en cuenta sus valores.

Seguimos sacando resultados. Ahora muestro el llamado circulo de adherencia de los neumáticos, donde se aprecia la fuerza generada en el neumático, lateral en X y longitudinal en Y.

Se sacan conclusiones muy interesantes de como trabajan los neumáticos, invito a los visitantes que las estudien un poco y verán cosas muy curiosas.

Seguimos con el reparto de carga sobre cada rueda. La segunda muestra el porcentaje de carga que soporta la rueda delantera, la trasera soportará lo restante.

Finalmente, y por no cargar mucho esto, podemos ver la fuerza total y par al que se ve sometida la pipa de dirección.

Un siguiente estudio consistirá en hacer un análisis estadístico de los resultados y sacar algunas coclusiones de las tendencias del modelo e intentar mejorar poco a poco.

Mañana ya retomo el trabajo en mi empresa, la que me da de comer, asi que dedicaré mucho menos tiempo a esto pero no lo dejaré aparcado, continuaré con ello, realizando mejoras y aprendiendo un poco más sobre como funciona una moto.... es muy interesante y llegas a entender un poco mejor esas máquinas que tanto nos gustan pilotar.

Espero que os guste... ciao a tutti.

Escrito por Iván Martinez

04 agosto 2011

Seguimos trabajando con el modelo

Aunque con ciertas complicaciones en el tema personal ya que me vuelven a operar de la clavícula, he dedicado bastante tiempo a "pulir" el modelo matemático de predicción de la dinámica del chasis de la Ducati 1098.





Ya tengo terminadas las condiciones de estudio plano, es decir, en situación de aceleración y frenada, ya que el paso por curva llevará más tiempo por la dificultad que tiene modelar el neumático en esas condiciones. Se necesitan muchos parámetros y características que solamente con datos del fabricante o con ensayos podría conseguir. Otra solución es ir ajustando esos parametros hasta llevarlos a darnos soluciones coherentes y adaptadas a las medidas realizadas con la telemetría en pista, cuando lo tenga más avanzado daré más detalles.

De momento, las simulaciones lanzadas, se adaptan muy bien a los valores recogidos en la telemetría, permitiendome conocer en todo instante (en aceleración y frenada por el momento) toda la física de la moto, geometría, cinemática, esfuerzos soportados por todos los componentes ... TODO.
De momento el modelo es algo "rigido" en el sentido de que no existe interfaz para la entrada de datos, si no un código de programación en el que resulta un poco laborioso modificar por ejemplo la precarga del muelle, clicks de compresión de horquilla, avances.... para ver como impactan en el comportamiento final de la moto. No es que sea largo, si no que un trabajo final podría ser hacer un interfaz sencillo y gráfico donde cualquier usuario pudiese explotar el modelo y hacer su puesta a punto o analisis del chasis.

Voy a intentar listar los parametros que entran en el modelo y que dan una idea del nivel de detalle hasta donde, por el momento, he llegado.

Se han discretizado las masas, inercias y sus centros de gravedad en todos estos elementos:

- Rueda trasera y delantera con discos de freno.
- Buje trasero, pinza de freno trasero y corona de transmisión.
- Baslulante.
- Amortiguador, bieleta de regulación de altura y palanca de progresividad.
- Chasis, cableado, radiador y periféricos unidos al chasis.
- Motor.
- Depósito de combustible con la variante de cantidad de combustible en su interior.
- Tijas, rodamientos...
- Botellas y barras de horquilla con pinzas de freno.
- Subchasis con colas de escape.
- Piloto.

Existe la posibilidad de variar cualquier geometría de la moto, longitud de basculante, regulación de la bieleta de regulación de altura, batalla, posición de cualquier elemento antes mencionado con el fin de variar el c.d.g. global del conjunto, offset de las horquillas, coeficientes aerodinámicos de carenados o capacidad de acople del piloto, altura eje de basculante, diferentes desarrollos con coronas/piñones diferentes para ver el impacto que esto produce en el squat por el "tiro" de la cadena en aceleraciones.... Esto solo es una pincelada de variables que tengo en cuenta, si entramos en suspensiones, tengo modelado:

- La rigidez del muelle amortiguador y horquillas.
- Los clicks de compresión y exensión de los hidraulicos de horquillas y amortiguador (a pulir con datos del fabricante)
- La precarga de los muelles.
- El nivel de aceite en las horquillas.
- La rigidez y presión de inflado de los neumáticos

... y más...

- Rigidez torsional del chasis.
- Coeficientes aerodinámicos.
- Densidad del aire dependiendo donde se ruede, altura, temperatura....



Con las primeras simulaciones lanzadas, se ha discretizado el chasis como elemento aislado en elementos finitos, metiendo como condiciones de contorno los datos obtenidos de Matlab como son:

- Todas las masas, inercias impactadas en el chasis.
- La aceleración longitudinal (aceleración y frenada) importada de la telemetría.
- Las fuerzas resistentes aerodinámicas.
- Las fuerzas y par en la pipa de dirección.
- La fuerza de tensión de la cadena de transmisión.
- Las fuerzas generadas por el sistema de progresividad de suspensión trasera.

Finalmente mido reacciones en el eje del basculante y verifico que coinciden con la fuerza que actua sobre el chasis, y efectivamente, así ha sido.

Ahora mismo puedo predecir como se "retuerce" el chasis en una aceleración y frenada, los niveles de deformaciónes, tensiones, las barras que mejor trabajan, donde se puede optimizar para tener la mejor relación peso-rigidez.
En mi opinión esta es la vía de desarrollo que se debería tomar en el diseño de un chasis y que las empresas más punteras de la industria utilizan para la optimización de sus productos y llegar a ser los mejores y destacar del resto.
Espero equivocarme pero tengo la sensación de que la competición se basa más en la experiencia que en el uso de estas herramientas y la ingeniería.

Volviendo al modelo, me encantaría haberlo podido validar con unos ensayos instrumentando el chasis con galgas extensométricas y dando unas vueltas en pista. Sería muy gratificante el certificar que el modelo matemático es correcto o al menos tomar las medidas oportunas para corregirlo.

Hace unos años trabajé en la instrumentación de los comandos de alerones para los primeros ensayos de vuelo de un avión con galgas extensométricas y me ayudó a coger algo de experiencia en este campo, tipos de montaje, configuración eléctrica, puente de Wheatstone a implimentar... y a priori, creo que resultaría algo complicado en la moto, ya que el nivel de tensiones encontrados son bajos, muy lejos del límite elástico del acero, y en estas condiciones con deformaciones tan bajas es dificil medir.
En estas condiciones lo que hacia es crear un amplificador mecánico para aumentar el nivel de deformaciones en la zona a medir y mejorar los resultados y resolución del ensayo. En el caso del chasis.... habría que pensarlo un poco, pero seguro que existe una solución.

Me encantaría que algún constructor de chasis estuviese interesado en seguir esta linea de trabajo e invirtiese un poco de dinero en investigar todo esto, ya que se necesitan medios que un pelamanillas como yo no tiene.

Aunque hacer este tipo de cosas me es gratificante, sobre todo en estos momentos de inutilidad a causa de la puñetera clavícula, a la vez es un poco frustrante el que queden en saco roto y no puedan ser explotadas... en fin.

Escrito por Iván Martinez