Metodo de Dugoff Caracterizacion de Neumaticos
Metodo de Dugoff Caracterizacion de Neumaticos
Metodo de Dugoff Caracterizacion de Neumaticos
ÍNDICE DE CONTENIDOS
1. Documento de aceptación
2. Nomenclatura
3. Análisis del problema
3.1. Contexto
3.2. Metodología
3.3. Modelo de análisis usado: modelo de Dugoff
4. Caracterización del neumático
4.1. Método analítico
4.2. Método empírico
4.3. Diferencias entre modelos
5. Modelado del neumático usando el modelo de Dugoff
5.1. Consideraciones iniciales: parámetros
5.2. Problema propuesto
5.2.1. Cálculo de rigidez radial y al giro
5.2.2. Variación de la fuerza de frenado en función del deslizamiento
longitudinal
5.2.3. Variación de la fuerza transversal variando el ángulo de deriva
5.2.4. Variación de la fuerza longitudinal en función de la fuerza transversal
variando el deslizamiento longitudinal
5.3. Conclusiones
6. Bibliografía
1. Documento de aceptación
Yo, Pedro José García Villanueva, manifiesto que he participado en la redacción de esta
práctica, con título "Caracterización de neumáticos", correspondiente a la unidad
didáctica 1 de la guía docente de la asignatura que lleva por nombre "Neumáticos"1.
Así mismo, manifiesto que he leído, corregido y estoy de acuerdo con todo lo expuesto
en las páginas que siguen a este documento. Porcentaje aproximado del trabajo
realizado: 60%
Yo, Iván Martínez Benito, alego que he sido participe de este documento con título
“Caracterización de neumáticos” correspondiente a la práctica 1 de la asignatura
Automóviles. También alego que he leído la totalidad de este documento y estoy de
acuerdo con todo lo expuesto en él. Porcentaje aproximado del trabajo realizado: 40%
1
Guía completa en: http://www.upv.es/titulaciones/GIM/menu_921106c.html
Modelado por Dugoff Automóviles 34O
4
Práctica 1- Caracterización de neumáticos
Automóviles
2. Nomenclatura
Símbolos genéricos
3.1 Contexto
En este mundo en el que el transporte por carretera dispone de más de mil millones de
vehículos, sin contar maquinaria de obra o de minería, el estudio del comportamiento
de los neumáticos ha sido, es y será una parte fundamental en los nuevos avances en
la automoción.
Ver, suponer y descubrir la forma en la que algo tan simple como es un disco que gira,
se ha convertido en objeto de innumerables estudio, análisis y modelos matemáticos
para acercarnos lo más posible a saberlo todo sobre los neumáticos. Si queremos
mejorar la seguridad, confort y prestaciones de nuestros vehículos, el avance en el la
fabricación y diseño es un camino que obligatoriamente debemos travesar.
3.2 Metodología
También asume que la carcasa del neumático es indeformable así como que el ángulo
de caída de este es nulo. Tampoco existe par de autoalineamiento y las rigideces a la
deriva y al deslizamiento son constantes. La ventaja de este modelo radica en la
sencillez final de las expresiones analíticas de la fuerza longitudinal y lateral, aunque
estas solo son válidas para deslizamientos puros.
2
*En nuestro caso, Ki será igual a la rigidez radial del neumático.
3
*En nuestro estudio, Kα será la rigidez al giro calculada multiplicada por la carga Fz.
De tal forma:
Siendo
Siendo el valor de
Siendo
Y con
Con todo lo anterior claro, debemos saber de qué tipo es el modelo que vamos a usar,
si es analítico o empírico. En los siguientes capítulos se analizaran ambos tipos,
retratando pros y contras de ellos.
Modelo de Fiala
Modelo de Bakker, Nyborg y Pacejka
Modelo de la "Fórmula mágica" de Bakker y Pacejka
Modelo 521-Tire de TNO Automotive
Modelo UA-Tire de la Universidad de Arizona
A parte de estos modelos, existen muchos más o variaciones de los ya creados, como
el "Dugoff modificado" de Nenggen Ding and Saied Taheri(2).
Entre estos modelos empíricos usados en el mundo del automóviles estarían el modelo
de deformación por galgas extensiométricas o el análisis de vibraciones por medio de
foto receptores, entre muchos otros.
Las diferencias entre el modelo analítico y el modelo empírico las hemos resumido en
los puntos más críticos entre ellos.
Modelo analítico
o Permite cambiar datos y parámetros rápidamente
o El tiempo de cálculo depende de la complejidad, no obstante suele ser
menor que el empírico
o El coste de estos modelos es mucho menor que en el empírico, no se
usa aparataje extra
o Los resultados obtenidos pueden no ser totalmente verídicos
o Pueden llegar a ser muy complejos
Modelo empírico
o Su uso es caro y se ha de disponer de herramientas e instalaciones
preparadas para tal fin
o El tiempo de ejecución suele ser mayor que el analítico ya que hay que
realizar varias pruebas para hacer el error de medición más pequeño
o Los resultados, si se realizan bien, muestran la que ocurre en la realidad
En este momento ya tenemos todos los datos para la obtención de lo que se nos pide,
la práctica consta de 4 apartados A, B, C y D los cuales serán estudiados a
continuación.
Hay dos tipos de rigidez radial, la estática que es la que tiene el neumático cuando
este está quieto y la dinámica que es la que tiene cuando está rodando. En nuestro
caso analizaremos la rigidez radial estática, Kt. Esta se determina por las pendientes de
la curva de carga estática-deformación
Para obtener la rigidez radial estática, solo necesitamos saber la carga y la presión de
inflado. En nuestro caso a analizar, tenemos un vehículo de pasajeros que monta un
neumático 165/60 R13 de estructura radial a una presión de 30 psi soportando una
carga vertical de 600Kg.
La relación Fyα-α es prácticamente lineal cuando los valores del ángulo de deriva son
pequeños (α < 4°). Esto corresponde a aceleraciones laterales pequeñas y es
despreciable el deslizamiento lateral. A medida que la fuerza lateral incrementa ya sea
por viento o por otros factores, el deslizamiento lateral aumenta decreciendo
rápidamente la pendiente de la curva Fyα-α. La fuerza transversal seguirá creciendo
junto con α hasta el límite de adherencia, produciendo que la rueda derrape ya sea en
aceleración o frenado.
La rigidez de deriva, Kα, se define por la derivada de la fuerza transversal con respecto
al ángulo de deriva, considerando nula la influencia del ángulo de caída de la rueda.
En este apartado B, se nos pedía graficar lo que ocurría con la fuerza de frenado según
fuésemos variando i. Para llevar a cabo esta tarea se han usado las ecuaciones que nos
aportaba el modelo de Dugoff además de tener en cuenta ciertos puntos.
Para el cálculo, como bien hemos dicho, usaremos las ecuaciones de Dugoff, estas
ecuaciones son las explicadas en el punto 3.3 y son:
Siendo
Entonces
Con todo lo anterior resuelto, la gráfica que nos queda después de operar es:
3600
3250
Fuerza longitudinal N
2900
50 km/h
2550
80 km/h
2200
1850
1500
0% 2% 4% 6% 8% 10% 12% 14% 16% 18% 20% 22% 24% 26% 28% 30% 32% 34% 36% 38% 40%
Deslizamiento longitudinal %
3000
2500
Fuerza longitudinal N
2000
1500
1000
150 km/h
500
0
0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100%
Deslizamiento longitudinal %
Vemos que ahora la pendiente es más acusada, y que el pico de máxima adherencia se
reduce a un zona más estrecha de deslizamiento. Esta zona corresponde con la zona de
uso del ABS, el cual, por medio de la electrónica, juega con la fuerza de frenado para
evitar que las ruedas patinen y se mantengan lo más posibles pegadas al asfalto.
Los datos en forma de tabla de las gráficas anteriores se pueden resumir en la
siguiente:
i vi1 vi2 μa v1 μa v2 Fx1 Fx2 i vi1 vi2 μa v1 μa v2 Fx1 Fx2
0% 0.000 0.000 0.600 0.600 0.0 0.0 21% 2.917 4.667 0.583 0.572 3362.3 3302.8
1% 0.139 0.222 0.599 0.599 1683.5 1683.5 22% 3.056 4.889 0.582 0.571 3361.4 3298.9
2% 0.278 0.444 0.598 0.597 2610.2 2607.3 23% 3.194 5.111 0.581 0.569 3360.1 3294.7
3% 0.417 0.667 0.598 0.596 2917.0 2911.2 24% 3.333 5.333 0.580 0.568 3358.5 3290.2
4% 0.556 0.889 0.597 0.595 3068.0 3059.2 25% 3.472 5.556 0.579 0.567 3356.7 3285.3
5% 0.694 1.111 0.596 0.593 3156.5 3144.8 26% 3.611 5.778 0.578 0.565 3354.6 3280.3
6% 0.833 1.333 0.595 0.592 3213.9 3199.2 27% 3.750 6.000 0.578 0.564 3352.3 3275.0
7% 0.972 1.556 0.594 0.591 3253.5 3235.8 28% 3.889 6.222 0.577 0.563 3349.8 3269.5
8% 1.111 1.778 0.593 0.589 3282.0 3261.3 29% 4.028 6.444 0.576 0.561 3347.2 3263.9
9% 1.250 2.000 0.593 0.588 3303.0 3279.3 30% 4.167 6.667 0.575 0.560 3344.4 3258.1
10% 1.389 2.222 0.592 0.587 3318.8 3292.1 31% 4.306 6.889 0.574 0.559 3341.4 3252.2
11% 1.528 2.444 0.591 0.585 3330.9 3301.2 32% 4.444 7.111 0.573 0.557 3338.3 3246.2
12% 1.667 2.667 0.590 0.584 3340.1 3307.4 33% 4.583 7.333 0.573 0.556 3335.2 3240.0
13% 1.806 2.889 0.589 0.583 3347.1 3311.5 34% 4.722 7.556 0.572 0.555 3331.9 3233.7
14% 1.944 3.111 0.588 0.581 3352.4 3313.8 35% 4.861 7.778 0.571 0.553 3328.5 3227.4
15% 2.083 3.333 0.588 0.580 3356.4 3314.8 36% 5.000 8.000 0.570 0.552 3325.0 3220.9
16% 2.222 3.556 0.587 0.579 3359.2 3314.7 37% 5.139 8.222 0.569 0.551 3321.5 3214.4
17% 2.361 3.778 0.586 0.577 3361.1 3313.6 38% 5.278 8.444 0.568 0.549 3317.8 3207.8
18% 2.500 4.000 0.585 0.576 3362.3 3311.8 39% 5.417 8.667 0.568 0.548 3314.1 3201.1
19% 2.639 4.222 0.584 0.575 3362.8 3309.3 40% 5.556 8.889 0.567 0.547 3310.4 3194.4
20% 2.778 4.444 0.583 0.573 3362.8 3306.3 - - - - - - -
El análisis se hizo con las dos velocidades de estudio, además de velocidad 0 km/h y se
compararon con gráficas de neumáticos reales. Los datos que se obtuvieron son los
que siguen:
A simple vista podemos ver poco de lo que nos está aportando cada valor, para poder
verlos mejor se han graficado las dos series de estudio y una serie más de referencia a
velocidad 0:
3600
3400
Fuerza transversal N
3200
3000
2800
2000
0 2 4 6 8 10 12 14
Ángulo de deriva º
En el gráfico se pueden observar 3 zonas, marcadas por las líneas verticales negras, a la
izquierda de la primera línea el comportamiento es lineal y conforme se aumenta el
ángulo lo hace la fuerza. Entre las líneas ocurre la zona de transición, en la que la rueda
esta en un estado crítico de deslizamiento. A la derecha de la segunda línea ocurre el
deslizamiento para cualquier valor de α (3).
3000
Fuerza lateral N
2000
50 km/h
1000
80 km/h
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Ángulo de deriva
Variación de la fuerza lateral en función del ángulo de deriva en ejes totales. Vemos
que la variación de 50 a 80 km/h es bastaste pequeña. Si las comparamos con curvas
teóricas vemos que los resultados obtenidos por el método de Dugoff son muy
acertados:
Podemos apreciar que para la zona estudiada de 0º a 14º de deriva se cumplen ambas
curvas. Teniendo en cuenta que las presiones y cargas de [1] son diferentes a las
estudiadas podemos ver una gran similitud. Empieza por una zona lineal hasta la de
deslizamiento pasando por la de transición.
Como punto de partida empezamos calculando las fuerzas horizontales, hecho ahora
fácil a que teníamos la experiencia de los apartados anteriores. A partir de ahí nos tocó
obtener los esfuerzos longitudinales, los cuales se obtuvieron de dos formas distintas.
3500
3000
2500
Esfuerzo longitudinal
2000
1500
1000
500
0
0 500 1000 1500 2000 2500 3000 3500
Esfuerzo transversal
También nos resultó interesante ver como influía el hecho de usar la fuerza de frenado
obtenida por el método de Dugoff o por el método del círculo de Kamm, esto nos
arrojaba que el modelo de Dugoff se alejaba del de Kamm una cantidad, no obstante
dicha cantidad se reducía conforme aumentaba i.
3500
3000
2500
Esfuerzo longitudinal
2000
1500
1000
500
0
0 500 1000 1500 2000 2500 3000 3500
Esfuerzo transversal
Fuerza longitudinal calculada por el método del círculo de Kamm y transversal por
Dugoff.
-1000
-1200
-1400
-1600 Diferencia
-1800
Deslizamiento longitudinal
Como punto final se quiso ver que ocurría si el coeficiente de adherencia se disminuía,
si se cumpliría con el círculo de Kamm y su radio descendería de valor. Se eligió como
nuevo μ un valor de 0.3, haciendo referencia a una calzada en malas condiciones o
unos neumáticos peores que los del caso anterior. Se obtuvo la siguiente gráfica:
3500
3000
2500
Esfuerzo longitudinal
2000
1500
1000
500
0
0 500 1000 1500 2000 2500 3000 3500
Esfuerzo transversal
Se ve que la nueva línea de fuerza máxima, línea roja, se reduce bastante respecto al
caso previamente estudiado, esto nos cumple con la realidad, de que malos
neumáticos o calzadas reducen las fuerzas máximas que se pueden transmitir. Los
datos del apartado primero con un deslizamiento hasta 40% son los siguientes:
3500
3000
Fuerza transversal N
500
0
0 1000 2000 3000 4000
Fuerza de frenado N
5.3 Conclusiones
6. Bibliografía y referencias
http://e-
archivo.uc3m.es/bitstream/handle/10016/10785/PFC_Carlos_Heras_Barras.pdf
2. A Modified Dugoff Tire Model for Combined-slip Forces: Nenggen Ding and Saied
Taheri- Tire Science and Technology 2010 38:3, 228-244
http://tiresciencetechnology.org/doi/abs/10.2346/1.3481696
http://e-
archivo.uc3m.es/bitstream/handle/10016/10818/PFC_Lorena_Hernandez_Penafiel.pd
f
http://arxiv.org/ftp/arxiv/papers/1008/1008.5041.pdf