Ups CT004599
Ups CT004599
Ups CT004599
SEDE CUENCA
_______________________________
de Ingeniero Electrónico.
_______________________________
Autores:
Director:
Ing. Néstor Diego Rivera Campoverde
Nosotros, Cristian Paúl Calle Jara y Darwin Gonzalo Vicuña Pacheco, declaramos bajo
juramento que el trabajo aquí descrito es de nuestra autoría; que no ha sido previamente
presentado para ningún grado o calificación profesional; y, que hemos consultado las
referencias bibliográficas que se incluyen en este documento.
Firma:
II
CERTIFICACIÓN
Firma:
III
AGRADECIMIENTO
Primero que nada agradezco a Dios, por haberme permitido llegar a cumplir esta meta,
por guiarme por el buen camino y haberme acompañado todos los días en mi dificultosa
labor.
motivación constante que me ha permitido plasmar como una persona de bien, pero
más que nada por acompañarme durante todo este arduo camino y compartir conmigo
alegrías y fracasos.
A mi madre Bertha Jara y a mi Padre Enrique Calle, por los ejemplos de perseverancia,
consejos, pero sobre todo por los valores que me han inculcado siempre, por haberme
de igual modo a los profesores por compartir sus conocimientos, en especial al Ing.
A mi hermano Freddy Calle por ser un gran amigo para mí, que junto a sus ideas hemos
Cristian Calle
IV
AGRADECIMIENTO
Agradezco a Dios, a mi familia y a mis amigos que siempre me han dado su apoyo
incondicional para realizar esta tesis, a la Universidad Politécnica Salesiana por darme
que con sus enseñanzas me formaron como un profesional con calidad humana, en
especial a el Ing. Néstor Rivera, que con su paciencia y conocimientos nos guió el
presente trabajo, y a todas las personas que de una u otra forma permitieron el
Darwin Vicuña
V
DEDICATORIA
fuerza necesaria. También dedicarles con mucho cariño a mis hijos Emily y
profesional. A mi hermano Freddy por su apoyo en los estudios, y a todos los familiares
en general, por el apoyo que siempre me brindaron día a día en el transcurso de cada
Cristian Calle
VI
DEDICATORIA
La elaboración de esta tesis la dedico a Dios quién es el motor de mi vida, quien siempre
estuvo a mi lado guiando mis pasos y dándome fuerzas para seguir adelante ante las
brindarme todo su amor, paciencia, comprensión y sus consejos que me han permitido
ser una persona de bien con principios, valores y perseverancia para conseguir mis
objetivos. A mis hermanas, Paola y Nelly; por la motivación constante que me brindan
para cumplir mis metas, a mi sobrino Matías por ser mi mayor inspiración, y a todas
las personas que estuvieron apoyándome y bridándome su ayuda para hacer realidad
este sueño.
Darwin Vicuña
VII
RESUMEN
Un motor de combustión interna, es una máquina que obtiene energía mecánica a partir de la
energía química contenida en un combustible que arde dentro de una cámara de combustión.
Estos motores se denominan así porque realizan su trabajo en el interior de la cámara de
combustión mediante la aportación del calor producido al quemarse el combustible.
Árbol de levas: llamado también eje de levas, es el que regula el movimiento de las válvulas
de admisión y de escape, determinando el tiempo en que la válvula de admisión este abierta
para el ingreso de la mezcla aire combustible, como el tiempo en que la válvula demora en
cerrarse en el proceso de compresión de la mezcla.
Válvula de Admisión: es aquella que permite que la mezcla aire combustible, entre en la
cámara de combustión del motor para que se efectúe el tiempo de admisión.
Válvula de Escape: es aquella que permite expulsar al medio ambiente los gases de escape
que se generan dentro del cilindro del motor después que se quema la mezcla aire combustible.
Biela: es un elemento mecánico que es sometido a esfuerzos, permite la unión del pistón al
cigüeñal.
1
Figura 1. 1 Partes principales de motor de combustión interna.
Fuente. (Aprendamos Tegnologia. (s.f.). Obtenido de Motores de COMBUSTIÓN).
A los motores de gasolina de 4 tiempos se les conoce como Ciclo Otto y es el proceso que
sucede dentro de la cámara de combustión cuando entra la mezcla aire combustible, y dicha
mezcla sale quemada. Para todo el ciclo intervienen el pistón, válvulas tanto de admisión como
de escape, bujías, biela y cigüeñal (Klever & Solano, 2010).
Este ciclo está compuesto de 4 tiempos, cada tiempo permanece una carrera por lo que se
requiere que el pistón suba dos veces y baje dos veces generando dos vueltas completas del
cigüeñal.
El Punto Muerto Superior (PMS) es el punto más alto al que llega el pistón. El Punto Muerto
Inferior (PMI) es el punto más bajo que llega el pistón, y finalmente tenemos a la carrera que
se define como el recorrido que realiza el pistón entre el PMS al PMI.
2
El ciclo ideal de un motor Otto es de cuatro tiempos.
En este tiempo se abre la válvula de admisión para que entre la mezcla aire combustible, el
pistón inicia la carrera de aspiración desde el PMS al PMI, llenando así al cilindro de la
conocida mezcla, una vez que el pistón llegue al PMI se cierra la válvula de admisión.
A todo este proceso se considera como un evento mecánico ya que el gas no modifica sus
propiedades termodinámicas. (Núñez Orozco & González Oropeza, 2004)
3
Figura 1. 5 Tiempo de Compresión.
Fuente. (Ayala, D. (s.f.). Ciclo Otto Cuatro Tiempos).
Es el primer proceso termodinámico, llamado compresión adiabática, por ser rápida y no hay
tiempo a que se establezca la transmisión de calor. (Núñez Orozco & González Oropeza, 2004)
Este tiempo se da justo cuando el pistón llega al PMS. La bujía origina una chispa que inflama
la mezcla, esta explosión originada empuja al pistón al PMI. Al inicio de la explosión de dicha
mezcla la presión en el cilindro es máxima y el volumen mínimo, pero una vez que el pistón se
desplaza hacia el PMI transmite toda la potencia al cigüeñal, y la presión dentro del cilindro
disminuye mientras el volumen aumenta (Klever & Solano, 2010).
4
Figura 1. 7 a) Encendido de la mezcla b) Explosión de la mezcla.
Fuente. Ayala, D. (s.f.). Ciclo Otto Cuatro Tiempos.
Este tiempo es la última carrera del ciclo, es decir comienza a subir el pistón del PMI al PMS,
accionándose la válvula de escape para que los gases quemados que arrastra el pistón salgan
del cilindro, en este tiempo el cilindro disminuye el volumen y la presión no aumenta.
5
Figura 1. 9 Tiempo de Escape.
Fuente. Ayala, D. (s.f.). Ciclo Otto Cuatro Tiempos.
6
1.1.4 Combustión a presión constante.
Es llamado así en honor del ingeniero Alemán Rudolf Diésel, en este tipo de motor la
combustión se realiza de manera diferente a lo que ocurre en un motor de gasolina. La
combustión no se produce por una chispa en el interior de una cámara, si no que se aprovecha
las propiedades químicas que posee el gasóleo, el encendido se produce por una
alta temperatura que posibilita la compresión del aire al interior del cilindro.
El ciclo ideal de un motor a Diésel puede ser de dos y cuatro tiempos el que se analizara será
de cuatro tiempos.
El pistón se desplaza desde el PMS (Punto muerto superior) al PMI (punto muerto inferior),
con la válvula de admisión abierta aspirando aire de la atmosfera, aumentando la cantidad de
aire en la cámara.
7
“Esto se modela como una expansión a presión constante (ya que al estar la válvula abierta la
presión es igual a la exterior). En el diagrama Presión, Volumen aparece como una recta
horizontal” (Motor de cuatro tiempos. Ciclo Diesel).
Este es el segundo tiempo, aquí las válvulas tanto de admisión como de escape se encuentran
cerradas y el pistón se mueve hacia arriba comprimiendo aire, según se vaya comprimiendo el
aire la temperatura va aumentando considerablemente.
“Dada la velocidad del proceso se supone que el aire no tiene posibilidad de intercambiar calor
con el ambiente, por lo que el proceso es adiabático. Se modela como la curva adiabática”
(Laplace, 2013).
8
Figura 1. 16 Tiempo de Compresión A-B.
Fuente: Autores.
1.1.4.3 Tiempo de Combustión.
En este tercer tiempo siguen cerradas las válvulas de admisión y escape, mientras tanto un poco
antes de que el pistón llegue al PMS y hasta un poco después que el pistón comienza a
descender, se inyecta el combustible en la cámara a una presión elevada.
9
1.1.4.4 Tiempo de Expansión.
Al finalizar el proceso de la inyección se produce una Expansión adiabática por ser un proceso
muy rápido, las válvulas tanto de admisión como de escape se mantienen cerradas.
En este tiempo se abre la válvula de escape y procede a expulsar los gases producidos por la
combustión, empujado por el pistón.
10
Figura 1.20 Tiempo de Escape.
Fuente. Autores.
Teniendo en cuenta que la cantidad de aire que sale y la cantidad de aire que entra en el cilindro
es idéntica podemos considerar que ha sufrido un proceso de enfriamiento, el cual se produce
en dos fases.
Cuando el pistón se encuentra en el punto más bajo ósea el PMI (Punto muerto
inferior), el volumen se mantiene aproximadamente constante y se representa en el
diagrama como la isócora D-A.
Cuando el pistón empuja el aire hacia el exterior con la válvula de escape abierta se
llama isóbara A-E. Con lo que se cierra el ciclo.
11
En la figura de a continuación se observa el ciclo Diésel ideal.
“Observando el ciclo Diésel ideal, podemos considerar despreciables los procesos de admisión
y de escape a presión constante A-E y E-A, puesto que son idénticos en la gráfica y de sentido
opuesto, por lo que el calor y el trabajo intercambiados entre ellos se anulan mutuamente”
(Laplace, 2013).
12
1.2.1 Admisión.
Las pérdidas se dan debido a las limitaciones que encuentra el flujo de aire hasta llegar al
interior del cilindro, en la cual presión no se mantiene invariable. La limitación más indicadora
es la apertura tardía de las válvulas de admisión.
1.2.2 Compresión.
Al momento del llenado defectuoso del aire en las fases de admisión, dentro del cilindro se
tendrá una masa menor de aire al finalizar la compresión, por lo cual hay menor presión y
menor temperatura, factor que retarda el encendido de la mezcla.
1.2.3 Explosión.
1.2.4 Escape.
Cuando el pistón comienza ascender del PMI al PMS hay una cierta resistencia porque que
hay una presión interna en el cilindro, debido a la apertura retardada de la válvula de escape.
Las diferencias que se dan durante el proceso de la combustión interna entre el ciclo teórico y
el ciclo real, tanto en los motores de ciclo Diésel, como ciclo Otto, están dadas por:
Este es muy importante en el ciclo real, ya que el cilindro se encuentra refrigerado para asegurar
un correcto funcionamiento del pistón, una cierta parte de calor del fluido que se transmite a
las paredes, lo cual hace que se produzca una pérdida de trabajo útil tanto en las líneas de
compresión y expansión.
En el ciclo teórico la apertura y cierre de válvulas ocurre instantáneamente. Así, en el ciclo real
es imposible, realizar esta acción debido a un tiempo relativamente largo, por lo que, para
mejorar el llenado y vaciado del cilindro, las válvulas tanto de admisión como de escape se
abren anticipadamente, lo que provoca una pérdida de trabajo útil (PorrasySoriano, 2014).
13
1.2.8 Pérdidas por bombeo.
Las cotas de reglaje ayudan a identificar la apertura y cierre de las válvulas de admisión
y escape, para analizar los diagramas de distribución y ver donde se dan los cruces de
válvulas.
14
1.3.1 Adelanto en la apertura de la admisión (AAA).
Es el momento en el cual la válvula de admisión se abre antes que el pistón llegue al PMS,
después de haber completado la carrera de escape (Castro, 2013).
1.3.2 Retraso en el cierre de la admisión (RCA).
Debido a la inercia de los gases al final de la admisión estos siguen entrando en el cilindro, por
lo que válvula de admisión se cierra un poco después de que el pistón llegue a su PMI.
Este momento se da para compensar el tiempo necesario para que al final de la combustión, el
movimiento del pistón en su fase de trabajo sea mínimo.
La válvula de escape se abre completamente antes de que el pistón comience hacer el barrido
de los gases, debido a que puede haber presión en los gases al momento que el pistón comienza
a desplazarse hacia el PMS, por lo que hay pérdidas de energía.
15
1.3.6 Sincronización de las válvulas.
Antes que comience el pistón a bajar del PMS al PMI en el tiempo de admisión la válvula de
admisión se abre y continúa abierta mucho después del P.M.I. es decir en pleno tiempo de
compresión, esto se da para rendir los gases entrantes. El punto de cierre de la válvula de
admisión nos muestra la relación de comprensión segura. La compresión del motor será menor
cuando la válvula de admisión se cierre tarde.
Antes que termine el tiempo de explosión la válvula de escape se abre para liberar los gases
que están en el cilindro. La potencia del motor no es afectada por las válvulas de escape cuando
se abren en ese punto, debido a que la mayor parte de la potencia de los gases en explosión ha
sido transmitida al pistón durante el tiempo de explosión. Cuando el pistón alcance su velocidad
máxima, la válvula de escape debe estar abierta en su totalidad para que no haya resistencia al
movimiento, la cual nos daría perdidas en el bombeo (Ugalde, 2007).
Los ángulos de adelanto y retraso de las válvulas de admisión y escape (cotas de reglaje), varían
según la fabricación del motor, para las cotas de reglaje del motor que se va hacer las pruebas
se describe en el capítulo 2.
16
CAPÍTULO 2
2.1.1 Sensores.
El sensor de posición de cigüeñal es también conocido como CKP por las siglas de CranKshaft
Position Sensor, el cual es un dispositivo de efecto Hall que detecta la posición de cigüeñal y
la velocidad (rpm) del motor, esta información es enviada a la computadora del automóvil para
calcular el tiempo de apertura de los inyectores y el tiempo de encendido.
17
2.1.3 Sensor de posición del árbol de levas (CMP).
El sensor de posición de árbol de levas es también conocido como CMP por las siglas de
CranKshaft Position Sensor, es aquel que permite conocer el momento exacto en que el pistón
de referencia (cilindro número uno) se encuentra en el PMS (Punto Muerto Superior),
convirtiéndole en voltaje para que la ECU (unidad de control de motor) pueda medir el tiempo
de la inyección secuencial del combustible. Este sensor proporciona una señal cada 720 grados
de rotación del cigüeñal (Loaiza, 2013).
La adquisición de datos o adquisición de señales, consiste en medir variables del mundo real
(sistema analógico) para convertirlas en formato digital (sistema digital), almacenarlas en un
computador y procesarlas según la necesidad.
18
2.1.6 Sensores utilizados para adquirir la posición del árbol de levas y posición de
cigüeñal.
Los sensores que utilizamos para adquirir las señales para su respectivo análisis son los
siguientes CKP y CMP.
El hardware para la adquisición de datos que nosotros utilizamos es una tarjeta (NI USB-
6212) de National Instruments, ya que sus características satisfacen a nuestras necesidades.
19
Figura 2. 6 Configuración para adquisición de señales sensor CKP, CMP.
Fuente. Autores.
Figura 2. 7 Conexión de los sensores con la (DAQ) para adquisición de datos CKP, CMP.
Fuente. (Hurtado, Morocho & Pinguil, 2014).
El procesamiento de señales tanto del sensor CKP como CMP se realizó mediante el software
LabVIEW.
20
LabVIEW (Laboratory Virtual Instrumentation Engineering Workbench), utiliza un lenguaje
de programación de alto nivel tipo gráfico. El lenguaje que usa se llama lenguaje G, donde la
G simboliza que es lenguaje Gráfico.
Los programas que son desarrollados en este software se llaman VI´s (Virtual Instruments),
estos programas se dividen en dos partes, la una es denominada “Panel Frontal”, y otra
“Diagrama de Bloques”.
Panel Frontal: es utilizado para interactuar con el usuario cuando el programa se está
ejecutando. Aquí el usuario pude manipular, controlar y monitorear datos. En esta interfaz
se definen los controles (que son utilizados como entradas, pueden ser botones, marcadores
etc.) e indicadores (que son utilizados como salidas, pueden ser gráficas, etc.).
Diagrama de Bloques: es el programa donde se define su funcionalidad del VI, aquí se
observa la estructura del programa, se colocan íconos que realizan una determinada función
según las necesidades del usuario, y se interconectan, (el código que controla el programa).
Las características que se tomaron en cuenta para utilizar este software son las siguientes:
21
En la figura 2.9 podemos observar la señal del sensor CMP.
El automóvil que utilizamos es el SANTA FE 2.0. Este automóvil posee 60 dientes menos 2
dientes faltantes es igual a 58 dientes en la rueda fónica.
22
Figura 2. 11 Rueda fónica del cigüeñal.
Fuente. Autores.
Estos dos dientes faltantes son detectados por la ECU y son utilizados para determinar el punto
muerto superior y sincronizar el sistema.
Por cada diente se genera un ciclo de una señal alterna es decir que en una vuelta completa de
la rueda fónica se tendrá 58 ciclos de la señal y 2 ciclos sin señal.
Cada ciclo de un motor se realiza en dos giros del cigüeñal por lo que tiene que girar 720 grados
para completar las cuatro fases (admisión, compresión, expansión y escape).
Según los giros que realiza el cigüeñal podemos decir que un diente equivale a seis grados.
El punto muerto superior (PMS) está situado 19 dientes después del diente grande debido a la
velocidad del procesamiento de la ECU.
El punto muerto inferior (PMI) está situado 30 dientes después del diente 19, es decir del punto
muerto superior (PMS).
Para obtener las cotas de reglaje de nuestro motor se procedió a desarmarlo para poder tener
acceso a las válvulas de admisión y escape, para de esta manera identificar los eventos
transitorios del motor, mediante el giro del cigüeñal. El proceso para llegar a obtener esto fue
el siguiente:
23
Primero se observa las mangueras que conectan con el intercooler y se marcan para no conectar
erróneamente cuando se vuelva a colocar el intercooler. Se procede a desconectar las
mangueras y desmontar el intercooler.
De igual forma se desconectan los sockets de los inyectores y las cañerías de combustible para
poder retirar el inyector posteriormente.
Siguiendo con el proceso se debe retirar el perno del seguro de cada inyector, remover el fijador
interno con una varilla o destornillador y por ultimo sacar el inyector manualmente.
24
Figura 2. 15 Ubicación Perno de seguridad del Inyector.
Fuente. Autores.
Procedemos a retirar los pernos de la tapa de válvulas tanto de admisión como de escape.
Se identifica que par de válvulas son las de admisión y las de escape, se coloca las válvulas del
primer cilindro en compresión y el pistón en el punto muerto superior mediante la marca del
cigüeñal con la carcasa.
25
Figura 2. 17 Ubicación de Válvulas de Admisión y Escape.
Fuente. Autores.
Para identificar las cotas de reglaje se empieza a girar el cigüeñal en forma horaria y se observa
cuando las válvulas comienzan su apertura y cierre, además de ir anotando el grado en que
sucede cada evento.
Notamos que al girar el cigüeñal desde que el cilindro 1 está en compresión y en el PMS, la
primera válvula en abrirse es la válvula de escape 32º antes del PMI.
26
Luego continuando con el giro del cigüeñal observamos que la válvula de admisión es la
próxima en abrirse a 5º antes de que el pistón llegue al PMS mientras la válvula de escape
continua abierta.
Seguimos girando el cigüeñal y vemos que la válvula de escape se cierra a 5º después de que
el pistón haya pasado el PMS, por lo que deducimos que el cruce de válvulas es de 10º.
Para finalizar el análisis de las cotas de reglaje en las válvulas observamos que la válvula de
admisión se cierra a los 51º después del PMI.
27
Figura 2. 22 Ubicación de cierre de válvula de Admisión.
Fuente. Autores.
28
En la siguiente figura se indica los eventos transitorios del motor de combustión interna,
representado en cotas de reglaje.
Al momento de que el motor hace la combustión este genera ruido y VIBRACIONES por lo
que el sensor Knock convierte este ruido en voltaje debido a que dicho sensor está compuesto
por un material piezoeléctrico, esté material son químicos que generan voltaje cuando se aplica
presión o alguna vibración.
El sensor Knock siempre va estar transmitiendo señal, debido a que el motor siempre está
vibrando pero al momento de la detonación va haber mayor vibración en el motor por lo que
hay mayor voltaje en la salida del sensor.
29
Figura 2. 26 Comportamiento señal sensor Knock.
Fuente. (Padilla, Pulla, 2014).
En la figura 2.26 se puede visualizar el sensor Knock (Detonación) que utilizamos para la
adquisición de la señal de detonación durante el proceso de la combustión en el motor.
30
2.3.2 Ubicación del sensor Knock (Detonación).
El sensor Knock fue colocado en el block en forma perpendicular al motor mediante un pad
especial, adherido al motor para que no se caiga durante que el motor está en funcionamiento.
Cuando el sensor es único se coloca sobre el cilindro uno o dos y cuando el sensor es doble se
coloca sobre el cilindro uno o dos y tres o cuatro.
En la figura 2.28 se puede observar la ubicación del sensor.
31
2.3.4 Señal de Micrófono (Ruido).
El micrófono se utilizó para analizar el ruido que se genera durante el proceso de combustión
en motores de combustión interna.
En el vehículo se da una serie de fuentes responsables de producir ruido como por ejemplo en
el Motor, Neumáticos, Sistema de escape, etc.
Motor: El ruido que se da en el motor puede definirse como ruido de combustión y ruido
mecánico.
Neumáticos: El ruido que se da en los neumáticos al tener contacto con la calzada se debe a
varios factores como la velocidad a que circula el vehículo, las características que posee el
neumático, superficie del suelo, etc.
Sistema de Escape: el ruido que se genera en el escape es producido por los gases de admisión
y escape que se da durante el proceso de combustión.
El ruido mecánico es producido por el cabeceo del pistón, sistema de inyección, holguras en
los cojinetes, fallas en los engranes, etc.
El micrófono que se utilizó para realizar la adquisición de las señales acústicas que se dan en
el motor al momento de que se realiza la combustión fue el micrófono piezoeléctrico
unidireccional de marca PCBR modelo HT378B2, cuyas características son las siguientes.
MICRÓFONO HT378B02
DESCRIPCIÓN CARACTERÍSTICAS
Rango de frecuencia 6-10000Hz
Sensibilidad 50Mv/Pa
Diámetro nominal del 1/2”
micrófono
Rango de frecuencia (+/- 1bB) 6-10000Hz
Rango de frecuencia (+/- 2bB) 3.5-200000Hz
Rango de temperatura de -400C a 1200 C
operación
Límite de frecuencia más baja 1-2.5Hz
En este apartado se analizaran las normas para la toma de señales de ruido en el motor.
Existen diferentes normas de medición de presión sonora como la SAE, ISO, DIN, nosotros
tomaremos como referencia para nuestro respectivo análisis la norma ISO, específicamente la
norma 3745 la cual especifica métodos para medir los niveles presión acústica sobre una
32
superficie de medición que envuelve una fuente de ruido (maquina o equipo) en cámaras
anecoicas y semianecoicas.
Esta norma específica diferentes métodos para medir los niveles de presión acústica sobre una
superficie de medición que envuelve una fuente de ruido (maquinaria o equipo) en cámaras
anecoicas y semianecoicas. El tamaño máximo de la fuente de ruido depende de los requisitos
especificados con respecto al radio de la esfera hipotética o hemisferio utilizado como
superficie de medición envolvente.
Los métodos especificados en esta norma internacional son válidos para todos los tipos de ruido
(estacionario, no estacionario, etc.). También esta norma establece que la temperatura debe
estar en un rango de 150 C a 300 C (Instituto Ecuatoriano de Normalización, 2014). (Ver Anexo
2).
Para realizar la medición de ruido en el motor el micrófono debe ser colocado según la norma
ISO 3745 en una esfera en el caso de que la medición se realice en una cámara anecoica, o una
semiesfera en caso de una cámara semianecoica. El centro de la esfera o semiesfera debe
coincidir con el centro del motor.
33
Figura 2. 31 Ubicación de micrófono.
Fuente. (Guarango, Lazo, 2014).
De acuerdo a la norma ISO 3745 se dice que el radio de la semiesfera de medición debe ser
igual o mayor a la dimensión de la fuente de ruido, debe ser por lo menos de 1metro.
Las dimensiones del motor SANTA FE 2.0 que utilizamos es de 110cm de largo y 75 cm de
alto.
34
Para la ubicación del centro acústico (centro geométrico), se idealiza al motor como un
paralelepípedo, en el cual el punto central será el centro acústico.
35
ACELERÓMETRO UNIAXIAL ERBESSD
DESCRIPCIÓN CARACTERÍSTICAS
Modelo ACS 3411LN
Marca ACS
Grado de protección IP 67
Sensibilidad 330 mV/ g
Contiene un acelerómetro unidireccional regulador de
voltaje, un amplificador y un filtro que limpia la señal.
Posee una fabricación a prueba de agua.
Ensamble fabricado de acero inoxidable.
Ideal para realizar trabajos al aire libre.
En la figura 2.35 se puede visualizar el acelerómetro uniaxial que utilizamos para la adquisición
de la señal de VIBRACIONES en el motor.
36
2.3.10.1 Norma ISO 10816.
Esta norma habla sobre la vibración mecánica donde se evalúa las VIBRACIONES en una
máquina mediante medidas en partes no rotativas.
Los criterios de vibración de este estándar se aplican a un conjunto de máquinas con potencia
superior a 15 kW y velocidad entre 120 RPM y 15.000 RPM. Los criterios son sólo aplicables
para VIBRACIONES producidas por la propia máquina y no para VIBRACIONES que son
transmitidas a la máquina desde fuentes externas. Las mediciones deben realizarse cuando el
rotor y los descansos principales han alcanzado sus temperaturas estacionarias de trabajo y con
la máquina funcionando bajo condiciones nominales o específicas (por ejemplo de velocidad,
voltaje, flujo, presión y carga), (NORMA ISO 10816-1995, 2014), (Ver Anexo 3).
Para realizar las mediciones hay que ubicar el sensor en una posición y dirección adecuada
para lo cual se tomó en cuenta normas estandarizadas las mismas que se indica en la figura
2.36
37
2.3.12 Adquisición de la señal del acelerómetro Uniaxial (VIBRACIONES).
Primero se lee las señales que envía los sensores CKP y CMP ya que estas señales son las que
permiten identificar las fases de admisión, compresión, explosión y escape. La señal del sensor
CMP es la que da la referencia de inicio del ciclo de combustión, mientras que la señal CKP es
la indica las fases de combustión.
Después de haber adquirido las señales de los sensores, se digitaliza a las señales entre un valor
comprendido de 0 a 1.
38
Una vez que se tiene la señal CKP con los valores entre 0 y 1 se procedió a realizar el código
en matlab para el conteo de dientes. Cada vez que la señal toma un valor de 1 y 0 se cuenta
como un diente en la señal CKP, para la señal CMP se hace de igual forma cuando toma un
valor de 1 y 0 es un ciclo.
Después se necesita saber cuándo debe iniciar el conteo, es aquí donde se implementa la señal
CMP. Cuando la señal CMP conto un ciclo, inmediatamente la señal CKP también comenzara
a contar.
Para el recorte de la señal CKP en las fases que se desea, se creó una variable con el nombre
de (recorte), esta variable tiene valores de 0 donde no se desea visualizar la información y 1
donde si se desea visualizar, en otras palabras es un ventana para visualizar las fases.
A esta ventana le damos el número de diente para que comience a contar y el número de diente
para que termine el conteo.
Figura 2. 40 Ventana para analizar las fases de admisión, compresión, explosión y escape.
Fuente. Autores.
Para analizar la señales de ruido, vibración y detonación en las distintas fases del motor
simplemente se multiplico la ventana con la con cada una de estas señales.
39
En la figura 2. 42 se puede observar un diagrama de flujo de los procesos involucrados para
el desarrollo del algoritmo.
INICION
Digitalizar CMP
Conteo Ciclo
Conteo Ciclo =1
no
si
Cuenta dientes de
CKP
Se crea la
variable
Ventana
Analizar=Knock
Señal=Ventana x Analizar Analizar=Acelerómetro
Analizar=Micrófono
FIN
40
CAPÍTULO 3
Para el análisis de eventos transitorios en el motor de combustión interna se tomó diez muestras
de cada señal como se mencionó anteriormente en condiciones normales del motor, a las
señales adquiridas se les aplicó diferentes tipos de ventanas, para luego calcular la media,
varianza, desviación estándar, mediana, máximo, mínimo, potencia, energía, factor curtosis,
asimetría, rms.
En procesamiento de señales, un ventaneo se aplica cuando nos interesa analizar una señal de
longitud voluntariamente limitada. Por lo que la señal real tiene que ser de tiempo finito (A.V.
Oppenheim, 2008).
41
Si la ventana seleccionada es muy ancha se tendrá una buena resolución en frecuencia pero una
mala resolución en tiempo, una ventana de ancho infinito o ancho igual al tamaño de la señal
es la TF clásica.
Un defecto de la STFT es el no poder dar una alta resolución tanto en tiempo como en la
frecuencia de manera simultánea. La raíz de este problema se basa en el principio de
incertidumbre de Heisenberg, el cual establece que es imposible conocer una representación
exacta tiempo-frecuencia de una señal, es decir no podemos saber qué valor de frecuencia
existe en un instante de tiempo determinado, sólo podemos conocer qué componentes de
frecuencia existen dentro de un intervalo de tiempo dado (Cardona, 2009).
𝐾
𝑉ℎ(𝑘) = 0.5 (1 − 𝑐𝑜𝑠 (2𝜋 )) , 𝑘 = 1, … … … 𝑛
𝑛+1
42
Figura 3. 2 Ventana Hanning tiempo y frecuencia.
Fuente. (Quintero 2014).
Para valores fuera del rango 0 a n-1, estos poseen valor de cero.
La ventana Hamming es muy similar a la de Hanning, pero su respuesta en frecuencia variará.
1 𝑛−(𝑁−1)/2 2
− ( )
𝑣(𝑛) = 𝑒 2 𝜎(𝑁−1)/2 , 𝜎 ≤ 0,5
43
Para seleccionar el tipo de ventana debemos considerar que tipo de señal se va a analizar, a
continuación se presenta la tabla 3 se muestra como seleccionar una ventana según el contenido
de la señal.
Desconocida. Hanning
El usuario tiene beneficios al momento de seleccionar la ventana más idónea para la aplicación
de la STFT a la señal que desea analizar, si tiene conocimiento del comportamiento de las
diferentes ventanas las mismas que se muestran en la tabla 3.
Media Aritmética.
La media aritmética o promedio, muestra el valor central de los datos constituyendo ser la
medida de ubicación que más se utiliza. En general, es calculada sumando los valores de interés
y dividiendo entre el número de valores sumados.
(𝑋1 + 𝑋2 + 𝑋3 + ⋯ + 𝑋𝑛)
𝑋̅ =
𝑛
Varianza.
Constante que representa una medida de dispersión media de una variable aleatoria X, respecto
a su valor medio. La varianza se define como:
2
∑𝑛𝑖=1(𝑋𝑖 − 𝑋̅)2
𝛿 =
𝑛
Desviación estándar.
La desviación estándar (s) mide cuánto se separan los datos y se calcula sacando la raíz
cuadrada de la varianza.
44
∑𝑛 (𝑋𝑖 − 𝑋̅)2
𝛿 = √ 𝑖=1
𝑛
Mediana.
Es el valor que ocupa el lugar central de todos los datos cuando estos están ordenados de menor
a mayor, la mediana solo se puede hallar para variables cuantitativas. La mediana se calcula de
la siguiente manera.
𝑀𝑒 = 𝑋(𝑛+1)/2
Máximo.
Es el valor máximo que toman los datos.
Mínimo.
Es el valor mínimo que toman los datos.
Potencia.
Es una magnitud directamente proporcional al trabajo, e inversamente proporcional al tiempo,
ósea el trabajo que se puede hacer en cada unidad de tiempo.
Energía.
Se define como el área limitada por la gráfica de magnitud al cuadrado de la señal y el eje del
tiempo, se trata siempre de una cantidad positiva.
Factor Curtosis.
El Coeficiente de Curtosis analiza el grado de concentración que presentan los valores
alrededor de la zona central de la distribución.
1
(𝑛) ∗ ∑(𝑋𝑖 − 𝑋𝑚)4 ∗ 𝑛𝑖
𝑔= −3
1 2 2
((𝑛) ∗ ∑(𝑋𝑖 − 𝑋𝑚) ∗ 𝑛𝑖))
Asimetría.
La Asimetría de una distribución hace referencia al grado en que los datos se reparten por
encima y por debajo de la tendencia central.
Interpretación de la gráfica 3.4
45
Figura 3. 5 Asimetría de una muestra.
Fuente. (T Seoane, JLR Martin 2014).
RMS.
La media cuadrática (RMS) es muy útil para calcular la media de variables que toman valores
negativos y positivos. Se suele utilizar cuando el símbolo de la variable no es importante y lo
que interesa es el valor absoluto del elemento. El valor RMS de un conjunto de valores es la
raíz cuadrada de la media aritmética de los cuadrados de los valores originales.
DISEÑO FACTORIAL
Tipo de Ventana. Completa, Rectangular, Hanning, Hamming,
Gaussiana.
Ancho de la 12 grados menos y 12 grados más.
Ventana.
Fase. Admisión, Compresión, Explosión, Escape.
Variable. Vibración, Ruido, Knock.
Tabla 3. 2 Diseño Factorial.
Fuente. Autores.
Para el diseño experimental se varió el ancho de las ventanas de cada evento transitorio del
motor de combustión interna, para después realizar una comparación de resultados frente a
cada tipo de ventana.
Para lo cual aplicamos tres anchos de ventanas diferentes a cada evento de las señales del
Knock, Acelerómetro, Micrófono.
46
En el análisis realizado en la identificación de eventos se obtuvo los siguientes resultados:
En las siguientes figuras podemos visualizar cada evento con un tipo de ventana diferente en
la señal del sensor Knock (Detonación). Para observar todos los eventos con todos los tipos de
ventanas, en el sensor knock, Acelerómetro, Micrófono, (Ver Anexo 3).
47
Carrera de Escape tiene 214 grados.
El ancho de ventana se redujo 12 grados al inicio y al final en cada proceso tiendo los siguientes
grados:
En las siguientes figuras podemos visualizar cada evento con un tipo de ventana diferente en
la señal del Acelerómetro (Vibración). Para observar todos los eventos con todos los tipos de
ventanas, en la señal del Acelerómetro, Knock, Micrófono, (Ver Anexo 4).
48
Carrera de Explosión tiene 148 grados.
El ancho de ventana se incrementó 12 grados al inicio y al final en cada proceso tiendo los
siguientes grados:
En las siguientes figuras podemos visualizar cada evento con un tipo de ventana diferente en
la señal del Micrófono (Ruido). Para observar todos los eventos con todos los tipos de ventanas,
en la señal del Micrófono, Knock, Acelerómetro, (Ver Anexo 5).
49
Carrera de Compresión tiene 134 grados.
50
Cuando la variabilidad que tengo DENTRO de mis tratamientos es igual a la de ENTRE
tratamiento las medias se alinean verticalmente dando como resultado que las medias son
iguales, y por lo contrario la variabilidad que tengo DENTRO de mis tratamientos no son
iguales a la de ENTRE tratamiento las medias no se alinean verticalmente dando como
resultado que las medias son diferentes es aquí donde se da un análisis ANOVA para verificar
cuál de las hipótesis es la correcta.
SEÑALES ADQUIRIDAS
51
Para la señal de Detonación, en la carrera de admisión, compresión, explosión, escape, sacamos
la media, varianza, desviación Estándar, mediana, máximo, mínimo, potencia, energía, factor
curtosis, asimetría y rms, utilizando la ventana Completa, Rectangular, Hanning, Hamming,
Gaussiana.
Todos esos valores son exportados a Excel para realizar el análisis ANOVA.
Para el análisis Anova se utilizó el programa Minitab el cual da información tanto gráfica como
numérica, a cerca de datos estadísticos, de las muestras adquiridas de cada sensor para sacar
las respectivas conclusiones con respecto a la hipótesis.
3.3.1.1 Hipótesis.
Nuestra hipótesis es determinar cuál es el mejor tratamiento para analizar los eventos
transitorios en motor de combustión interna.
52
Hipótesis Nula es cuando todas las medias son iguales, no existe diferencia significativa entre
medias.
Hipótesis Alterna es cuando por lo menos una media es diferente, existe diferencia
significativa entre medias.
También llamado nivel de significación. Es la probabilidad (en tanto por uno) de fallar en
nuestra estimación, esto es, la diferencia entre la certeza (1) y el nivel de confianza (1-α). Por
ejemplo, en una estimación con un nivel de confianza del 95%, el valor α es (100-95)/100 =
0,05.
Completa, Gaussiana,
Hanning, Hamming,
VENTANA 5 Rectangular.
En general, mientras más alto sea R2, mejor se ajusta el modelo a los datos. R2 está siempre
entre 0 y 100%. Se conoce como el coeficiente de determinación o de determinación múltiple
(en regresión múltiple).
Al graficar valores observados por valores ajustados, se ilustran gráficamente los valores de R2
para modelos de regresión.
53
Figura 3. 21 Representación de respuestas observadas
Respuestas ajustadas para dos modelos de regresión.
Figura. Autores.
Residuo.
Los valores residuales indican el grado hasta el cual un modelo representa la variación en los
datos observados. Los puntos de esta gráfica deben generalmente formar una línea recta si los
residuos están normalmente distribuidos.
54
Figura 3. 23 Representación Residuos vs Ajustes.
Fuente. Autores.
Histograma de residuos.
Es una herramienta que muestra las características generales de los datos, incluyendo:
La presencia de largas colas en la gráfica podría indicar asimetría en los datos. Si una o dos
barras están lejos de las demás, esos puntos pueden ser valores atípicos. Debido a que la
apariencia del histograma cambia según el número de intervalos utilizados para agrupar los
datos (Minitab 2014).
Residuos - orden.
Esta es una gráfica de todos los residuos en el orden en el que se tomaron los datos y puede
utilizarse para hallar errores no aleatorios, especialmente de efectos relacionados con el tiempo.
Una correlación positiva es indicada por un conglomerado de residuos con el mismo signo.
Una correlación negativa es indicada por los cambios rápidos en los signos de residuos
consecutivos.
55
Figura 3. 25 Representación Residuos vs Orden.
Fuente. Autores.
Cuatro en uno.
Se seleccionó esta opción para generar una gráfica normal de residuos, un histograma de
residuos, una gráfica de residuos versus ajustes y una gráfica de residuos versus orden, en una
sola ventana de gráfica.
Gráfica de caja muestra en forma gráfica la distribución como la simetría de una muestra.
Valor Atípico (*) es la observación que se encuentra más allá del bigote superior o inferior.
56
Bigote Superior es cuando se extiende hasta el punto máximo de datos dentro de la altura de
la caja.
Cuartiles son medidas de posición que permiten conocer otros puntos característicos de la
distribución.
Línea superior Q3 (Tercer Cuartil), el 75% es el valor que sobrepasa al 75% de los valores
de la distribución.
Línea media Q2 (Mediana), el 50% intermedio de los datos.
Línea inferior Q1 (Primer Cuartil), es el valor mayor que el 25% de los valores de la
distribución.
57
3.3.2 Análisis Anova con Minitab para identificación de la Hipótesis.
Método
Hipótesis nula Todas las medias son iguales
Hipótesis alterna Por lo menos una media es diferente
Nivel de significancia α = 0,05
R-cuad. R-cuad.
S R-cuad. (ajustado) (pred)
0,0009163 97,86% 97,67% 97,36%
58
Figura 3. 30 Representación de caja de moda.
Fuente. Autores.
Para determinar que valores estadísticos nos sirven para nuestra hipótesis procedimos a realizar
una tabla donde se observa los valores de P y R.CUAD., donde el valor de P debe ser menor a
su nivel de significancia (α= 0.05) y R.CUAD debe ser mayor al 70% para considerar que estos
valores son significativos.
En la siguiente tabla se observa los valores de P y R. Cuadrado., para Admisión del sensor
Knock.
59
En el siguiente cuadro se observa los valores de P y R. Cuadrado, para compresión del sensor
knock.
En la siguiente tabla se observa los valores de P y R. Cuadrado para Combustión del sensor
Knock.
En la siguiente tabla se observa los valores de P y R. Cuadrado, para Escape del sensor Knock.
60
ENERGIA 0,000 78,29
CURTOSIS 0,000 40,67
ASIMETRIA 0,000 47,80
RMS 0,000 49,02
Tabla 3. 8 Valores de P y RAIZ CUADRADO para Admisión sensor Knock.
Fuente. Autores.
En la siguiente tabla se observa los valores de P y R. Cuadrado, para Admisión de la señal del
Acelerómetro.
61
En la siguiente tabla se observa los valores de P y R. Cuadrado, para Combustión de la señal
del Acelerómetro.
En la siguiente tabla se observa los valores de P y R. Cuadrado, para Escape de la señal del
Acelerómetro.
En la siguiente tabla se observa los valores de P y R. Cuadrado, para Admisión de la señal del
Micrófono.
62
MINIMO 0,000 75,24
POTENCIA 0,000 80,75
ENERGIA 0,000 98,10
CURTOSIS 0,000 52,31
ASIMETRIA 0,005 27,58
RMS 0,000 81,33
63
En la siguiente tabla se observa los valores de P y R. Cuadrado para Escape de la señal del
Micrófono.
Realizada ya las tablas con sus respectivos valores de P, R. Cuadrado, procedimos a realizar
una selección de los valores de P menor a 0,05 y R. Cuadrado mayor a 70% para realizar el
análisis de la hipótesis.
En la tabla 3.17 se muestra los valores del sensor Knock para determinar qué valor estadístico
sirvió para el análisis de la hipótesis.
En la tabla 3.18 se muestra los valores de la señal del Acelerómetro para determinar qué valor
estadístico servirá para el análisis de la hipótesis.
64
VIBRACIÓN ADMISIÓN COMPRESIÓN COMBUSTIÓN ESCAPE
MEDIA SI SI SI SI
VARIANZA SI SI SI SI
DESVIACION SI SI SI SI
MEDIANA SI SI SI SI
MAXIMO SI SI SI SI
MINIMO SI SI SI SI
POTENCIA SI SI SI SI
ENERGIA SI SI SI SI
CURTOSIS SI SI SI SI
ASIMETRIA SI SI SI SI
RMS SI SI SI SI
Tabla 3. 18 Valores Estadísticos de la señal del Acelerómetro.
Fuente. Autores.
En la tabla 3.19 se muestra los valores de la señal del Micrófono para determinar qué valor
estadístico servirá para el análisis de la hipótesis.
65
La señal del Acelerómetro (VIBRACIONES), nos indica que todos los valores estadísticos son
idóneos para el análisis de la hipótesis debido a que el sensor es normalizado.
La señal de Knock (Detonación), no indica que los valores estadísticos son pocos confiables
debido a que el sensor no es normalizado.
La señal de micrófono (Ruido), tiene una regularidad de datos estadísticos tanto para admisión
y combustión debido a la ubicación del micrófono.
66
Este procedimiento se realizó para la carrera de admisión, compresión, explosión y escape,
para las señales de Detonación, Vibración y Ruido. (Ver Anexo 8).
En la tabla 3.21 se observa las diferencias entre ventanas del sensor Knock, para la Media.
En la tabla 3.22 se observa las diferencias entre ventanas del sensor Knock, para la Mediana.
67
HANNING NO NO NO NO
HAMMING
RECTANGULAR SI SI SI SI
HAMMING
RECTANGULAR SI SI SI SI
HANNING
Tabla 3. 22 Diferencias entre ventanas sensor Knock, para mediana.
Fuente. Autores.
En la tabla 3.23 se observa las diferencias entre ventanas del sensor Knock para la Energía.
En la tabla 3.24 se observa las diferencias entre ventanas de la señal del Acelerómetro, para la
Media.
68
RECTANGULAR NO SI NO NO
COMPLETA
HAMMING SI SI SI SI
GAUSSIANA
HANNING NO NO NO NO
GAUSSIANA
RECTANGULAR SI SI SI SI
GAUSSIANA
HANNING SI SI SI SI
HAMMING
RECTANGULAR SI SI SI SI
HAMMING
RECTANGULAR SI SI SI SI
HANNING
Tabla 3. 24 Diferencias entre ventanas señal del Acelerómetro, para media.
Fuente. Autores.
En la tabla 3.25 se observa las diferencias entre ventanas de la señal del Acelerómetro, para la
Varianza.
69
En la tabla 3.26 se observa las diferencias entre ventanas de la señal del Acelerómetro, para la
Desviación.
En la tabla 3.27 se observa las diferencias entre ventanas de la señal del Micrófono, para la
Varianza.
70
HANNING NO NO NO NO
HAMMING
RECTANGULAR SI SI SI NO
HAMMING
RECTANGULAR SI SI SI NO
HANNING
Tabla 3. 27 Diferencias entre ventanas señal del Micrófono, para varianza.
Fuente. Autores.
En la tabla 3.28 se observa las diferencias entre ventanas de la señal del Micrófono, para la
desviación.
En la tabla 3.29 se observa las diferencias entre ventanas de la señal del Micrófono, para el
mínimo.
71
RECTANGULAR SI SI NO NO
COMPLETA
HANNING NO NO NO NO
GAUSSIANA
RECTANGULAR SI SI SI SI
GAUSSIANA
HANNING NO NO NO NO
HAMMING
RECTANGULAR SI SI SI SI
HAMMING
RECTANGULAR SI SI SI SI
HANNING
Tabla 3. 29 Diferencias entre ventanas señal del Micrófono, para Mínimo.
Fuente. Autores.
En la tabla 3.30 se observa las diferencias entre ventanas de la señal del Micrófono, para la
potencia.
En la tabla 3.31 se observa las diferencias entre ventanas de la señal del Micrófono, para la
energía.
72
ENERGIA ADMISIÓN COMPRESIÓN COMBUSTIÓN ESCAPE
GAUSSIAN SI SI SI SI
COMPLETA
HAMMING SI SI SI SI
COMPLETA
HANNING SI SI SI SI
COMPLETA
RECTANGULAR SI SI SI SI
COMPLETA
HAMMING NO NO NO NO
GAUSSIANA
HANNING NO NO NO NO
GAUSSIANA
RECTANGULAR NO NO NO NO
GAUSSIANA
HANNING NO NO NO NO
HAMMING
RECTANGULAR NO NO NO NO
HAMMING
RECTANGULAR NO NO NO NO
HANNING
Tabla 3. 31 Diferencias entre ventanas señal del Micrófono, para Energía.
Fuente. Autores.
En la tabla 3.32 se observa las diferencias entre ventanas de la señal del Micrófono, para RMS.
73
CAPÍTULO 4
Para Compresión y escape se analiza la señal del acelerómetro ya que es la que tiene mayor
prioridad en comparación de las señales del kanock y micrófono. Y finalmente se analiza la
combustión con la señal del sensor knock ya que es importante el momento de detonación en
dicha fase.
Como su nombre lo indica la representación de cajas y bigotes se hace por medio de una caja
rectangular donde los lados muestran el recorrido de la diferencia entre el primer y tercer cuartil
de una distribución. Esta caja está dividida por una línea que indica la mediana y la relación
entre el primer y tercer cuartil. Las líneas que están fuera de la caja se les conocen como bigotes
estos son los valores mínimos y máximos de la variable.
Para el cálculo de las cajas y bigotes se ordena la distribución, desde el número más pequeño
al más grande para un ejemplo se usa los datos estadísticos del sensor knock en la fase de
admisión para desviación estándar en la ventana completa.
Como se tiene 10 muestras N=10 y para calcular el primer cuartil (Q1) que es el 25% se coge
el primer 25% de las muestras que en este caso es el número 0,0169 y se suma más la siguiente
muestra que es 0,0184 dividido para 2
0,0165 0,0169 0,0184 0,020 0,0321 0,0357 0,0369 0,044 0,046 0,0494
74
Q1 = (0,0169 + 0,0184 ) / 2 = 0,0176
Para el segundo cuartil (Q2) que es el 50% de las muestra se hace de la misma forma de Q1, se
coge la muestra ubicada en el 50% de muestras y se suma más la siguiente y se divide para 2.
Este segundo cuartil va ser la media aritmética.
0,0165 0,0169 0,0184 0,020 0,0321 0,0357 0,0369 0,044 0,046 0,0494
Para el tercer cuartil (Q3) se hace de igual que Q1 y Q2 solo que ahora es el 75% de la muestra
0,0165 0,0169 0,0184 0,020 0,0321 0,0357 0,0369 0,044 0,046 0,0494
75
Para la fase de admisión se analizara las señales del micrófono:
Los valores más grandes de varianza y datos muy dispersos tiene la ventana rectangular. En el
ventaneo completo se encuentra la mediana más alta en comparación de las otras ventanas. Las
ventanas hanning, y gaussiana son las que tienen menor varianza y datos no muy dispersos.
Claramente se puede descartar que el ventaneo rectangular no sirva para el análisis
experimental de varianza en la fase de admisión, ya que los bigotes se traslapan con los valores
de las demás ventanas
La parte inferior de la caja de ventaneo completo es mayor que la parte superior, esto quiere
decir que la varianza comprendida entre el 25% y el 50% de los datos está más dispersa que
entre el 50% y el 75%, al igual que las ventanas gaussiana y hanning
Los datos muy dispersos y valores más grandes de desviación std se encuentra en la ventana
rectangular. En el ventaneo completo se encuentra la mediana más alta en comparación de las
otras ventanas.
Las ventanas hanning, y gaussiana son las que tienen menor desviación std y datos no muy
dispersos. La ventana rectangular no sirve para el análisis en admisión debido a que se traslapan
con las otras ventanas, la ventana completa es la que tiene más desviación std y por lo contrario
la ventana gaussiana y hanning son las que tienen menor desviación std.
76
Figura 4. 3 Representación de cajas de desviación std en Admisión.
Fuente. Autores.
La ventana rectangular no sirve para el análisis ya que se traslapan con las otras ventanas. Las
ventanas gaussiana y hanning son las que tienen menores mínimos, a comparación de la
ventana completa.
77
4.1.5 Potencia en Admisión para el Micrófono.
Los datos muy dispersos y valores más grandes de potencia se encuentran en la ventana
rectangular. En el ventaneo completo se halla la mediana más alta en comparación de las otras
ventanas. Las ventanas hanning, y gaussiana son las que tienen menor potencia y datos no muy
dispersos.
A igual que las anteriores variables estadísticas el ventaneo rectangular no sirve para el análisis
en potencia ya que se traslapan con las otras ventanas.
La ventana que tiene más potencia es la completa y la de menor potencia son las ventanas
hanning y gaussiana.
Los valores más altos de energía y media más alta tienen el ventaneo completo. En las otras
ventanas restantes tienen valores muy bajos de energía.
En energía el ventaneo completo es la que más tiene, y las demás ventanas están con los mismos
valores de energía.
78
4.1.7 RMS en Admisión para el Micrófono.
Los valores más grandes de RMS y datos muy dispersos tiene la ventana rectangular. En el
ventaneo completo se encuentra la mediana más alta en comparación de las otras ventanas. Las
ventanas hanning, y gaussiana son las que tienen menor RMS y datos no muy dispersos.
El ventaneo rectangular no sirve para el análisis en RMS ya que se traslapan con las otras
ventanas. La ventana que tiene más RMS es la completa y la de menor RMS son las ventanas
hanning y gaussiana.
Los valores más grandes de media y mediana se hallan en la ventana rectangular y también se
halla la mediana más alta de la media y mediana. Las otras ventanas son las que tienen menor
media y mediana.
79
Figura 4. 9 Representación de cajas de mediana en Compresión.
Fuente. Autores.
Los valores más grandes de varianza y desviación tiene la ventana hanning y se halla la
mediana más alta tanto en varianza como en desviación std. En el ventaneo completo y
rectangular hay menor varianza y desviación std que las demás ventanas. Los datos no son
muy dispersos en todas las ventanas.
80
Figura 4. 11 Representación de cajas de desviación en Compresión.
Fuente. Autores.
Los valores más grandes de potencia y RMS tienen las ventanas rectangular y completa
también se halla la mediana más alta tanto en potencia como en RMS. En las ventanas sobrantes
hay menor potencia y RMS. Los datos no son muy dispersos en todas las ventanas.
81
Figura 4. 13 Representación de cajas de potencia en Compresión.
Fuente. Autores.
82
Figura 4. 15 Representación de cajas de mínimo en Compresión.
Fuente. Autores.
Los valores más altos de energía y media tiene el ventaneo completo. En las otras ventanas
restantes tienen valores muy bajos de energía. Todas las ventanas no tienen datos muy
dispersos.
83
4.2.7 Factor Curtosis en Compresión para el Acelerómetro.
En la ventana rectangular hay mayor asimetría y se encuentra la media más alta, también tiene
los datos muy dispersos. La ventana completa es la que tienen menor asimetría y datos no muy
dispersos en comparación de la ventana rectangular.
84
4.3 Fase de Combustión.
4.3.1 Fase de Combustión con la señal del micrófono.
Los valores más grandes de varianza y datos muy dispersos tiene la ventana rectangular. En el
ventaneo completo se encuentra la mediana más alta en comparación de las otras ventanas. Las
ventanas hanning, y gaussiana son las que tienen menor varianza y datos no muy dispersos.
85
4.3.1.3 Mínimo en Combustión para el Micrófono.
Los datos muy dispersos y valores más grandes de potencia se encuentran en la ventana
rectangular. En el ventaneo completo se halla la mediana más alta en comparación de las otras
ventanas. Las ventanas hanning, y gaussiana son las que tienen menor potencia y datos no muy
dispersos en comparación de la ventana rectangular.
86
4.3.1.5 Energía en Combustión para el Micrófono.
Los valores más altos de energía y media tienen el ventaneo completo. En las otras ventanas
restantes tienen valores bajos de energía.
Los valores más grandes de RMS tiene la ventana rectangular. En el ventaneo completo se
encuentra la mediana más alta en comparación de las otras ventanas. Las ventanas hanning, y
gaussiana son las que tienen menor RMS y datos no muy dispersos. En cuanto a la dispersión
de datos la ventana completa tiene menor dispersión.
87
4.3.2 Fase de Combustión con la señal del sensor Knock.
En el ventaneo completo la media y mediana son las que tienen más bajos en comparación de
las otras ventanas. La ventana completa se solapa con los datos de la ventana rectangular lo
que es mejor tomar los datos de la ventana completa tanto en media como en mediana.
Las ventanas hanning, hamming y gaussiana son las que tienen mayores medias y medianas
con datos no muy dispersos, estas tres ventanas también se solapan los datos por lo que da igual
coger los datos de cualquiera de las tres ventanas gaussiana, hanning, hamming.
Las cruces que aparecen en el diagrama de cajas de media, son valores atípicos es decir,
observación inusualmente grande o pequeña, se consideran que una observación es un valor
atípico si es diferente de la media de la muestra en más del doble de la desviación estándar
agrupada.
88
4.3.2.2 Energía en Combustión para el Knock.
Los valores más altos de energía y media más alta tienen el ventaneo completo. En las otras
ventanas restantes tienen valores muy bajos de energía.
La de mayor energía es el ventaneo completo, y las ventanas restantes están con los mismos
valores de energía.
Los valores más grandes de media y mediana se hallan en la ventana rectangular y también se
halla la mediana más alta de la media y mediana. Las otras ventanas son las que tienen menor
media y mediana.
89
Figura 4. 29 Representación de cajas de mediada en Escape.
Fuente. Autores.
Los valores más grandes de varianza y desviación tiene la ventana hanning, también se halla
la mediana más alta tanto en varianza como en desviación std. En el ventaneo completo y
rectangular hay menor varianza y desviación std que las demás ventanas. Los datos no son
muy dispersos en todas las ventanas.
90
Figura 4. 31 Representación de cajas de desviación std en Escape.
Fuente. Autores.
Los valores más grandes de potencia y RMS tiene la ventana rectangular y completa también
se halla la mediana más alta tanto en potencia como en RMS. En las ventanas sobrantes hay
menor potencia y RMS que las demás ventanas. Los datos no son muy dispersos en todas las
ventanas.
91
Figura 4. 33 Representación de cajas de potencia en Escape.
Fuente. Autores.
92
Figura 4. 35 Representación de cajas de mínimo en Escape.
Fuente. Autores.
Los valores más altos de energía y media tienen el ventaneo completo. En las otras ventanas
restantes hay valores muy bajos de energía. Todas las ventanas no tienen datos muy dispersos.
93
Figura 4. 37 Representación de cajas de F. curtosis en Escape.
Fuente. Autores.
En la ventana rectangular hay menor asimetría y se encuentra la media más baja, también tiene
los datos muy dispersos. La ventana gaussiana es la que tienen mayor asimetría y datos no muy
dispersos en comparación de la ventana rectangular y completa.
94
CAPÍTULO 5
5. RESULTADOS Y ANÁLISIS.
En este capítulo se realizó la comparación de resultados frente a cada tipo de ventana, para lo
cual se realizó un análisis en Minitab de tipo de ventana, ancho de ventana, fase, variable, de
cada valor estadístico, para identificar que interacción hay entre ellas y concluir que tipo y
ancho de ventana sirve para cada proceso en motores de combustión interna.
5.1.1 Interacción.
En la figura 5.2 podemos observar los efectos principales para Factor de Cresta.
En el tipo de Ventana se puede observar que la ventana Rectangular es diferente a la completa,
las ventanas Hanning, Hamming, Gaussiana son similares entre ellas pero diferentes a la
Completa y Rectangular.
El ancho de ventana los valores del factor de cresta no varían significativamente ya sea el
ventaneo más grande o más pequeño con referencia al ventaneo normal.
95
Para la Fase la admisión se tiene un factor de cresta mayor a la fase de compresión, explosión
y escape, pero son similares entre ellas.
Las variables de Vibración, Ruido y Knock son diferentes, una con respecto a la otra.
En la figura 5.3 se pude observar la interacción entre tipo de ventana por ancho de ventana,
fase y variable. Ancho de ventana por fase y variable. Fase por variable.
En la figura 5.3 se analiza la interacción que existe entre ancho de ventana por fase, debido a
que en las otras no hay interacción entre ellas.
En Ancho de ventana por fase, la fase de compresión y explosión son similares, mientras que
el factor de cresta es menor en explosión cuando el ancho de ventana es pequeña, cuando el
ancho de ventana es grande el factor de cresta es menor en Escape.
96
En la figura 5.4 podemos observar los efectos principales para RMS.
El ancho de ventana los valores de rms no varían significativamente ya sea el ventaneo más
grande o más pequeño con referencia al ventaneo normal.
Para la Fase en compresión se tiene un valor de rms menor a la fase de admisión, explosión y
escape, pero son similares entre ellas.
Las variables de Ruido y Knock son parecidas una con respecto a la otra y diferentes a la
variable de Vibración.
En la figura 5.4 Se pude observar la interacción entre tipo de ventana por ancho de ventana,
fase y variable. Ancho de ventana por fase y variable. Fase por variable.
En la figura 5.4 Se analiza una interacción pequeña que existe entre ancho de ventana por fase,
debido a que en las otras no hay interacción entre ellas.
En Ancho de ventana por fase, el valor de rms es menor en escape cuando el ancho de ventana
es pequeña, cuando el ancho de ventana es grande el valor de rms es menor en Compresión.
97
Figura 5. 5 Representación de Interacción para RMS.
Fuente. Autores.
En la figura 5.6 podemos observar los efectos principales para el Factor de Curtosis.
En el tipo de Ventana se puede observar que la ventana Rectangular es diferente a la completa,
las ventanas Hanning, Hamming, Gaussiana son similares entre ellas pero diferentes a la
Completa y Rectangular.
El ancho de ventana los valores del factor de Curtosis no varían significativamente ya sea el
ventaneo más grande o más pequeño con referencia al ventaneo normal.
Para la Fase en escape se tiene un factor de Curtosis menor a la fase de admisión, explosión y
compresión, también se tiene un factor de Curtosis en compresión similar al de explosión, y
finalmente el factor de Curtosis en admisión es mayor con respecto a los demás fases.
Las variables de Ruido y Knock y Vibración, son diferentes entre ellas.
98
En la figura 5.7 Se pude observar la interacción entre tipo de ventana por ancho de ventana,
fase y variable. Ancho de ventana por fase y variable. Fase por variable.
En la figura 5.7 Se analiza la interacción que existe entre tipo de ventana por fase, tipo de
ventana por variable, debido a que en las otras no hay interacción entre ellas.
En tipo de ventana por fase, en fase de admisión el factor de curtosis es menor en la ventana
rectangular y mayor en la ventana Hanning, Hamming y Gaussiana con respecto a la fase de
explosión.
En tipo de ventana por variable, tenemos una intersección entre vibración y ruido, en la ventana
rectangular el factor de curtosis en la variable de ruido es menor a las demás variables y en las
ventanas Hanning, Hamming, Gaussiana el factor de curtosis es similar en las variables de
ruido, vibración y Knock.
99
Figura 5. 8 Representación de efectos principales para Energía.
Fuente. Autores.
En la figura 5.9 Se pude observar la interacción entre tipo de ventana por ancho de ventana,
fase y variable. Ancho de ventana por fase y variable. Fase por variable.
En la figura 5.9 no se analiza ningún tipo de interacción debido a que ningún efecto de un
factor depende del otro.
100
Para la Fase se tiene una potencia similar en admisión, compresión, explosión y escape.
Las variables de Ruido y Knock son diferentes con respecto a la vibración y similares entre
ellas.
En la figura 5.11 Se pude observar la interacción entre tipo de ventana por ancho de ventana,
fase y variable. Ancho de ventana por fase y variable. Fase por variable.
En la figura 5.11 no se analiza ningún tipo de interacción debido a que ningún efecto de un
factor depende del otro.
101
Las variables de Vibración y Ruido son diferentes mientras que la variable de ruido y knock
son similares.
En la figura 5.13 Se analiza la interacción que existe entre ancho de ventana por fase, debido a
que en las otras no hay interacción entre ellas. En Ancho de ventana por fase admisión,
explosión y escape son similares, mientras que el mínimo es menor en explosión cuando el
ancho de ventana es pequeña, cuando el ancho de ventana es grande el mínimo es menor en
compresión.
102
Figura 5. 14 Representación de efectos principales para Máximo.
Fuente. Autores.
En la figura 5.15 Se analiza la interacción que existe entre ancho de ventana por fase, debido a
que en las otras no hay interacción entre ellas.
En Ancho de ventana por fase admisión, explosión y escape son similares, mientras que el
máximo es menor en explosión cuando el ancho de ventana es pequeña, cuando el ancho de
ventana es grande el máximo es menor en compresión. No hay interacciones fuertes lo que
descarta los valores de máximo para nuestro análisis debido a que cualquier ventana que se
escoja y ancho vamos a tener los mismos resultados.
103
El ancho de ventana los valores de mediana no varían cuando es más grande o más pequeño el
ventaneo con referencia al ventaneo normal.
Para la Fase no hay una variación significativa entre admisión, comprensión, explosión y
escape todos los valores de mediana son muy similares.
Las variables de Vibración y Ruido son diferentes, mientras que la variable de ruido y knock
son similares.
La figura 5.17 Se analiza la interacción que existe entre ancho de ventana por fase, debido a
que en las otras no hay interacción entre ellas.
En Ancho de ventana por fase admisión, explosión y escape son similares, mientras que la
mediana es menor en explosión cuando el ancho de ventana es pequeña, cuando el ancho de
ventana es grande la mediana es menor en compresión. No hay interacciones fuertes lo que
descarta los valores de mediana para nuestro análisis debido a que cualquier ventana que se
escoja y ancho vamos a tener los mismos resultados.
104
En la figura 5.18 podemos observar los efectos principales para desviación estándar.
La ventana Rectangular es similar a la completa, las ventanas Hanning, Hamming, Gaussiana
son similares entre ellas pero diferentes a la Completa y Rectangular.
El ancho de ventana los valores de desviación no varían ya sea el ventaneo más grande o más
pequeño con referencia al ventaneo normal. Para la Fase todas dan el mismo resultado en el
valor estadístico de desviación. Las variables de Vibración y Ruido son diferentes mientras que
la variable de ruido y knock son similares.
En interacción para desviación estándar no es necesario para el análisis debido a que no hay
ninguna interacción en los cuadros de tipo de ventana, ancho de ventana y fase porque ningún
efecto de un factor depende del otro.
105
valor estadístico de varianza. Las variables de Vibración y Ruido son diferentes mientras que
la variable de ruido y knock son similares.
En interacción para varianza no es necesario para el análisis debido a que no hay ninguna
interacción en los cuadros de tipo de ventana, ancho de ventana y fase porque ningún efecto de
un factor depende del otro.
106
Figura 5. 22 Representación de efectos principales para Media.
Fuente. Autores.
En la figura 5.23 no se analiza ningún tipo de interacción debido a que ningún efecto de un
factor depende del otro.
107
Para la Fase, la asimetría tanto en admisión como en compresión son semejantes y diferentes
entre explosión y escape.
En la figura 5.25 Se pude observar la interacción entre tipo de ventana por ancho de ventana,
fase y variable. Ancho de ventana por fase y variable. Fase por variable.
En la figura 5.25 Se analiza la interacción que existe entre tipo de ventana por fase, tipo de
ventana por variable, ancho de ventana por fase y fase por variable, debido a que en las otras
no hay interacción entre ellas.
En tipo de ventana por fase, en fase de admisión la asimetría es menor en la ventana rectangular
y mayor en la ventana Hanning, Hamming y Gaussiana con respecto a la fase de comprensión.
También en fase de escape la asimetría es menor en la ventana rectangular y mayor en la
ventana Hanning, Hamming y Gaussiana con respecto a la fase de explosión.
En tipo de ventana por variable, tenemos una intersección entre knock y ruido, en la ventana
rectangular la asimetría en la variable de knock es menor a las demás variables y en las ventanas
Hanning, Hamming, Gaussiana la asimetría es similar en las variables de ruido, vibración y
Knock.
En Ancho de ventana por fase, la fase de explosión y escape son similares, y existe una
intersección entre compresión y admisión mientras que la asimetría es menor en compresión
cuando el ancho de ventana es pequeña, cuando el ancho de ventana es grande la asimetría es
menor en admisión.
En fase por variable existe una intersección entre vibración y ruido, la asimetría es menor en
ruido en la fase de explosión, y la asimetría es menor en vibración en la fase de escape.
108
Figura 5. 25 Representación de Interacción para Asimetría.
Fuente. Autores.
En la figura 5.26 se puede visualizar las predicciones para analizar qué tipo y ancho de ventana
sirve para cada evento en motores de combustión interna (Ver Anexo 9).
109
En la tabla 5.1 se tiene los resultados de las predicciones para el análisis de cada evento con la
variable de Vibración.
FACTOR DE CRESTA
Tipo de Ajuste Ancho Fase Variable
Ventana Ventana
Gaussiana 1,70541 262 G Admisión Vibración
Rectangular 2,24458 134 G Compresión Vibración
Gaussiana 1,74305 196 G Explosión Vibración
Gaussiana 1,73294 238 G Escape Vibración
RMS
Rectangular 1,93891 262 G Admisión Vibración
Rectangular 1,99522 86 P Compresión Vibración
Rectangular 1,96319 196 G Explosión Vibración
Rectangular 1,96690 190 P Escape Vibración
FACTOR DE CURTOSIS
Rectangular 3,09613 262 G Admisión Vibración
Rectangular 2,44381 86 P Compresión Vibración
Rectangular 3,86943 148 P Explosión Vibración
Rectangular 3,74690 190 P Escape Vibración
ASIMETRIA
Gaussiana 0,185454 262 G Admisión Vibración
Rectangular 0,169355 86 P Compresión Vibración
Gaussiana 0,209041 196 G Explosión Vibración
Gaussiana 0,180579 190 P Escape Vibración
VARIANZA
Hanning 0,50548 262 G Admisión Vibración
Hanning 0,493791 86 P Compresión Vibración
Hanning 0,467118 148 P Explosión Vibración
Hanning 0,471835 190 P Escape Vibración
Tabla 5. 1 Predicciones para análisis de eventos con Vibración.
Fuente. Autores.
En la tabla 5.2 se tiene los resultados de las predicciones para el análisis de cada evento con la
variable de Ruido.
FACTOR DE CRESTA
Tipo de Ajuste Ancho Fase Variable
Ventana Ventana
Gaussiana 4,26497 262 G Admisión Ruido
Gaussiana 3,15945 134 G Compresión Ruido
Rectangular 2,39303 196 G Explosión Ruido
Hanning 3,47348 238 G Escape Ruido
RMS
Rectangular 0,070817 215 P Admisión Ruido
Rectangular 0,129802 86 P Compresión Ruido
110
Rectangular 0,089905 148 P Explosión Ruido
Rectangular 0,091482 190 P Escape Ruido
FACTOR DE CURTOSIS
Gaussiana 6,10436 262 G Admisión Ruido
Gaussiana 3,74624 134 G Compresión Ruido
Gaussiana 3,97101 196 G Explosión Ruido
Hanning 4,24704 190 P Escape Ruido
ASIMETRIA
Gaussiana 1,32263 215 P Admisión Ruido
Gaussiana 1,36636 134 G Compresión Ruido
Rectangular 0,54771 196 G Explosión Ruido
Hanning 0,667385 190 P Escape Ruido
VARIANZA
Rectangular 0,004995 215 P Admisión Ruido
Rectangular 0,010875 134 G Compresión Ruido
Rectangular 0,007896 148 P Explosión Ruido
Rectangular 0,008335 190 P Escape Ruido
Tabla 5. 2 Predicciones para análisis de eventos con Ruido.
Fuente. Autores.
En la tabla 5.3 se tiene los resultados de las predicciones para el análisis de cada evento con la
variable de Knock.
FACTOR DE CRESTA
Tipo de Ajuste Ancho Fase Variable
Ventana Ventana
Hanning 4,91158 262 G Admisión Knock
Rectangular 4,41175 134 G Compresión Knock
Gaussiana 5,09005 196 G Explosión Knock
Rectangular 3,30713 190 P Escape Knock
Tabla 5. 3 Predicciones para análisis de eventos con Knock.
Fuente. Autores.
111
Gaussiana 3,97101 196 G Explosión Ruido
Hanning 4,24704 190 P Escape Ruido
FACTOR CRESTA
Tipo de Ajuste Ancho Fase Variable
Ventana Ventana
Hanning 4,91158 262 G Admisión Knock
Hamming 4,41175 134 G Compresión Knock
Gaussiana 5,09005 196 G Explosión Knock
Rectangular 3,30713 190 P Escape Knock
Tabla 5. 4 Priorización de predicciones para análisis de eventos con Vibración y Ruido.
Fuente. Autores.
Para analizar los eventos transitorios en motores de combustión interna se utilizó un tratamiento
de la señal con las siguientes características.
Se procede a sacar la FFT de la tabla 5.5, para analizar las señales en frecuencia de los eventos
de combustión interna frente a diferentes tipos de ventanas.
Para analizar la Detonación en Compresión se debe utilizar la ventana Hamming, con un ancho
de ventana de 134 grados.
En la figura 5.27 se puede observar la FFT de la señal de Detonación.
112
Figura 5. 27 FFT de Detonación.
Fuente. Autores.
Para el análisis de los picos se relacionó las frecuencias con la velocidad del motor en rpm
(revoluciones por minuto). La velocidad de motor en el que se tomó las muestras es de 830
rpm. Para identificar el pico de 55 Hz de la figura 5.27, se relacionó de la siguiente manera.
1Hz es igual a 60 rpm.
1𝐻𝑧
830𝑟𝑝𝑚 × = 13.8 𝐻𝑧 (𝐹𝑟𝑒𝑐𝑢𝑒𝑛𝑐𝑖𝑎 𝐹𝑢𝑛𝑑𝑎𝑚𝑒𝑛𝑡𝑎𝑙)
60𝑟𝑝𝑚
f = 13.8 Hz
Para los armónicos con 2f, 3f y 4f da como resultado 27.6 Hz, 41.4 Hz y 55.2 Hz
respectivamente, como se puede observar en los armónicos donde hay más alta la detonación
es en el cuarto armónico lo que representa la detonación en el motor.
En la figura 5.28 se puede observar la FFT de la señal de Detonación en Compresión, con una
ventana Hamming y un ancho de ventana de 134 grados.
113
En la figura 5.29 se puede observar la FFT de la señal de Detonación en Compresión, con una
ventana Hanning y un ancho de ventana de 134 grados.
Los resultados de la FFT tanto para la ventana Hamming como Hanning muestran que son
parecidos sus espectros debido a las características que tienen estas ventanas, ambas poseen
picos bien definidos para el análisis de la Detonación en Compresión.
Para realizar el análisis se utilizó los resultados de las tablas 5.4 y 5.5
Para analizar la Vibración en explosión se debe utilizar la ventana Rectangular, con un ancho
de ventana de 148 grados.
114
La velocidad de motor en el que se tomó las muestras es de 830 rpm. Para identificar el pico
de 27 Hz de la figura 5.30, se relacionó de la siguiente manera.
1𝐻𝑧
830𝑟𝑝𝑚 × = 13.8 𝐻𝑧 (𝐹𝑟𝑒𝑐𝑢𝑒𝑛𝑐𝑖𝑎 𝐹𝑢𝑛𝑑𝑎𝑚𝑒𝑛𝑡𝑎𝑙)
60𝑟𝑝𝑚
f = 13.8 Hz
Para los armónicos con 2f y 4f da como resultado 27.6 Hz y 55.2 Hz respectivamente, como se
puede observar en los armónicos donde hay más alta vibración es en el segundo armónico lo
que representa la vibración del motor.
En la figura 5.31 se puede observar la FFT de la señal de Vibración en Explosión, con una
ventana Rectangular y un ancho de ventana de 148 grados.
En la figura 5.32 se puede observar la FFT de la señal de Vibración en Explosión, con ventana
Hamming y un ancho de ventana de 148 grados.
115
Los resultados de la FFT tanto para la ventana Rectangular como Hamming muestran que la
ventana Hamming tiene un pico mayor en 35Hz y la ventana Rectangular tiene picos más
definidos lo que indica que para analizar la Vibración en Explosión se debe utilizar la ventana
Rectangular.
Para analizar el Ruido en Escape se debe utilizar la ventana Hanning, con un ancho de ventana
de 190 grados.
La velocidad de motor en el que se tomó las muestras es de 830 rpm. Para identificar el pico
de 55 Hz de la figura 5.33, se relacionó de la siguiente manera.
1𝐻𝑧
830𝑟𝑝𝑚 × = 13.8 𝐻𝑧 (𝐹𝑟𝑒𝑐𝑢𝑒𝑛𝑐𝑖𝑎 𝐹𝑢𝑛𝑑𝑎𝑚𝑒𝑛𝑡𝑎𝑙)
60𝑟𝑝𝑚
f = 13.8 Hz
Para los armónicos con 2f y 4f, da como resultado 27.6 Hz y 55.2 Hz respectivamente, como
se puede observar en los armónicos donde hay más alto ruido es en el cuarto armónico, lo que
representa el ruido del motor.
116
En la figura 5.34 se puede observar la FFT de la señal de Ruido en Escape, con una ventana
Hanning y un ancho de ventana de 190 grados.
En la figura 5.35 se puede observar la FFT de la señal de Ruido en Escape, con ventana
Gaussiana y un ancho de ventana de 190 grados.
Los resultados de la FFT tanto para la ventana Hanning como Gaussiana muestran que las dos
ventanas tienen picos bien definidos para analizar Ruido en Escape.
Para analizar el ruido en admisión se debe utilizar la ventana Gaussiana, con un ancho de
ventana 262 grados.
117
Figura 5. 36 FFT de Ruido.
Fuente. Autores.
La velocidad de motor en el que se tomó las muestras es de 830 rpm. Para identificar el pico
de 55 Hz de la figura 5.36, se relacionó de la siguiente manera.
1𝐻𝑧
830𝑟𝑝𝑚 × = 13.8 𝐻𝑧 (𝐹𝑟𝑒𝑐𝑢𝑒𝑛𝑐𝑖𝑎 𝐹𝑢𝑛𝑑𝑎𝑚𝑒𝑛𝑡𝑎𝑙)
60𝑟𝑝𝑚
f = 13.8 Hz
Para los armónicos con 2f y 4f, da como resultado 27.6 Hz y 55.2 Hz respectivamente, como
se puede observar en los armónicos donde hay más alto ruido es en el cuarto armónico, lo que
representa el ruido del motor.
En la figura 5.37 se puede observar la FFT de la señal de Ruido en Admisión, con una ventana
Gaussiana y un ancho de ventana de 262 grados.
118
En la figura 5.38 se puede observar la FFT de la señal de Ruido en Admisión, con ventana
Rectangular y un ancho de ventana de 262 grados.
Los resultados de la FFT tanto para la ventana Gaussiana como Rectangular muestran que la
ventana Rectangular tiene un pico mayor en 52Hz y la ventana Gaussiana tiene picos un poco
más definidos, lo que indica que para analizar el Ruido en Admisión se debe utilizar la ventana
Gaussiana.
En la figura 5.40 se puede observar la FFT de la señal de Detonación con falla en la inyección
en el motor.
119
Figura 5. 40 FFT de la señal de Detonación con falla en la inyección.
Fuente. Autores.
Como se puede ver la señal de detonación sin fallas tiene menos armónicos que la señal de
detonación con falla en la inyección, esto indica que el motor no funciona correctamente al
momento de la detonación, la FFT con falla no tiene una frecuencia constante de detonación
porque tiene algunos picos grandes. Mientras que en la FFT sin Falla tiene una frecuencia
constante de detonación que es 55.2 Hz.
En la figura 5.42 se puede observar la FFT de la señal de Detonación con falla en la inyección
en el motor.
120
En la señal de Vibración en condiciones normales como se ve en la figura 5.41 tiene una
vibración constante de 27.6 Hz es donde hay mayor vibración después disminuye la vibración
a 55.2 Hz y por último en 81 Hz, estos picos indican que la vibración es alta y va disminuyendo
a un solo ritmo en la frecuencia. En la señal con falla en la inyección la vibración no es
constante en el motor, vibra aleatoriamente y no hay un solo ritmo en frecuencia.
En la figura 5.43 se puede observar la FFT de la señal de Ruido en condiciones normales del
motor.
En la figura 5.44 se puede observar la FFT de la señal de Detonación con falla en la inyección
en el motor.
En la señal de Ruido en condiciones normales como se ve en la figura 5.43 tiene una amplitud
de 0.06 en 55.2 Hz, pero en la señal de ruido con falla tiene una amplitud de 0.12 en 18 y 21
Hz, lo que indica que la señal sin fallas suena distinta a la señal con fallas en la inyección.
121
CONCLUSIONES Y RECOMENDACIONES
6. CONCLUSIONES Y RECOMENDACIONES.
CONCLUSIONES:
Finalizado este proyecto de tesis podemos concluir lo siguiente:
Realizado el análisis Anova se llegó a establecer que todos los valores estadísticos de
la señal del acelerómetro (VIBRACIONES) son idóneos para el análisis de la hipótesis
debido a que el sensor es normalizado, la señal del micrófono (Ruido) permite analizar
la carrera de admisión y explosión debido a la ubicación del micrófono, y la señal de
Knock (Detonación) nos indicó que los valores estadísticos son pocos confiables
debido a que el sensor no es normalizado.
Para identificar qué tipo de ventana, ancho de ventana y tipo de variable sirve para un
mejor tratamiento de los eventos transitorios del motor, se realizó una Interacción la
cual determinan qué factores dependen de otros llegando a la conclusión de que el
mejor tratamiento para analizar eventos transitorios en motores de combustión interna
es:
122
Para la fase de Compresión es adecuado realizar un análisis con la variable de
Knock, Factor de Cresta y con una ventana Hamming de 134 grados.
Las fallas en el motor se pueden dar por diferentes variables, nosotros provocamos una
falla en la inyección del motor lo cual provoco grandes cambios en el funcionamiento
del mismo, la señal de los sensores fue diferente tanto en condiciones normales como
con falla. Analizando el espectro de la señal de micrófono se pudo diferenciar que el
ruido era diferente para cada caso.
RECOMENDACIONES:
Adentro de un proyecto tan anheloso como fue este, siempre se pretende que haya
mejora continua del mismo, por lo tanto se recomienda a futuros interesados, que el
software hecho en Labview, se puede seguir desarrollando para hacer diferentes
estudios necesarios en los motores de combustión interna, para así lograr encontrar
fallas y alargar la vida útil de dichos motores.
También se recomienda, que todas las señales que se van a ocupar para un análisis en
el motor de combustión interna sean gravadas, así se evita de estar prendiendo y
apagando el motor a cada instante, ya que así se desgasta la vida útil del motor. Con
esta recomendación también se evita de absorber los gases las personas que están
realizando el estudio en el motor, debido a que estos gases expulsados por el motor son
muy dañinos para la salud. En el caso de nosotros, otro motivo por el que se grava las
señales, es debido a que el estudio que se realizó fue en el taller de la Universidad
Politécnica Salesiana y había muchos alumnos que manipulaban el motor en el que se
estaba trabajando para el estudio, motivo por el cual no se encontraba siempre en las
mismas condiciones el motor
.
Con lo que se refiere al software se recomienda llenar todos los espacios en blanco que
se piden tanto el número de grados de inicio y final del ventaneo y el número de vuelta
del cigüeñal una vez que este todo lo necesario se manda a correr el programa.
123
El software que se realizó en Labview puede usarse tanto en motores de diesel como
de gasolina, es decir se podría ampliar el estudio para el análisis en los motores de
gasolina y diesel, para esto es necesario identificar muy bien los sensores de posición
del cigüeñal (CKP) y árbol de levas (CMP), estos sensores son muy importantes ya
que están sincronizados para realizar las fases de combustión del motor. También se
puede usar este software para motores de combustión interna de diferentes fábricas ya
que todos tienen el mismo funcionamiento.
124
REFERENCIAS BIBLIOGRÁFICAS
7. REFERENCIAS.
125
Laplace. (01 de 03 de 2013). Departamento de Fisica Aplicada III Universidad de
Sevilla. Obtenido de Ciclo Diesel:
http://laplace.us.es/wiki/index.php/Ciclo_Diesel
126
ANEXOS
8. ANEXOS.
ANEXO I
Características de DAQ (NI USB-6212).
General
Entrada Analógica
127
General
Canales Diferenciales 8
Rango -10 V - 10 V
Precisión 2.71 mV
Sensibilidad 118 µV
Precisión 0.089 mV
Sensibilidad 5.2 µV
Número de Rangos 4
Muestreo Simultáneo No
Salida Analógica
Número de Canales 2
Resolución 16 bits
Rango -10 V - 10 V
Precisión 3.512 mV
Rango -10 V - 10 V
Precisión 3.512 mV
128
General
E / S Digital
Canales Bidireccionales 32
Temporización Software
Entrada Digital
Salida Digital
Contadores / Temporizadores
Temporizador Watchdog No
Contadores 2
129
General
Operaciones a Búfer Sí
Generación de Pulso Sí
Tamaño 32 bits
Disparo Digital
Especificaciones Físicas
Longitud 16.9 cm
Ancho 9.4 cm
Altura 3.1 cm
130
ANEXO II
NORMA ISO 3745
131
132
133
134
135
136
137
ANEXO III
Norma ISO 10816-1995
VIBRACIÓN MECÁNICA. - EVALUACIÓN DE LA VIBRACIÓN EN UNA
MÁQUINA MEDIANTE MEDIDAS EN PARTES NO ROTATIVAS.
Los criterios de vibración de este estándar se aplican a un conjunto de máquinas con potencia
superior a 15 kW y velocidad entre 120 RPM y 15.000 RPM. Los criterios son sólo aplicables
para VIBRACIONES producidas por la propia máquina y no para VIBRACIONES que son
transmitidas a la máquina desde fuentes externas. El valor eficaz (RMS) de la velocidad de la
vibración se utiliza para determinar la condición de la máquina. Este valor se puede determinar
con casi todos los instrumentos convencionales para la medición de vibración.
Se debe prestar especial atención para asegurar que los sensores estén montados correctamente
y que tales montajes no degraden la precisión de la medición. Los puntos de medida
típicamente son tres, dos puntos ortogonales en la dirección radial en cada caja de descanso y
un punto en la medición axial.
138
Las mediciones deben realizarse cuando el rotor y los descansos principales han alcanzado sus
temperaturas estacionarias de trabajo y con la máquina funcionando bajo condiciones
nominales o específicas (por ejemplo de velocidad, voltaje, flujo, presión y carga).
En máquinas con velocidad o carga variable, las velocidades deben realizarse bajo todas las
condiciones a las que se espera que la máquina trabaje durante períodos prolongados de tiempo.
Los valores máximos medidos, bajo estas condiciones, serán considerados representativos de
la vibración. Si la vibración es superior a lo que el criterio permite y se sospecha de excesiva
vibración de fondo, las mediciones se deben realizar con la máquina detenida para determinar
el grado de influencia de la vibración externa. Si con la máquina detenida excede el 25% de la
vibración medida con la máquina operando, son necesarias acciones correctivas para reducir el
efecto de la vibración de fondo. En algunos casos el efecto de la vibración de fondo se puede
anular por análisis espectral o eliminando las fuentes externas que provocan las
VIBRACIONES de fondo.
Tipo de máquina.
Potencia o altura de eje.
Flexibilidad del soporte.
Grupo 1: Máquinas rotatorias grandes con potencia superior 300 kW. Máquinas
eléctricas con altura de eje H >= 315 mm.
Grupo 2: Máquinas rotatorias medianas con potencia entre 15 y 300 kW. Máquinas
eléctricas con altura de eje 160 =< H =< 315 mm.
Grupo 3: Bombas con impulsor de múltiples álabes y con motor separado (flujo
centrífugo, axial o mixto) con potencia superior a 15 kW.
Grupo 4: Bombas con impulsor de múltiples álabes y con motor integrado (flujo
centrífugo, axial o mixto) con potencia superior a 15 kW.
NOTA: La altura del eje H de una máquina está definida como la distancia medida entre la
línea de centro del eje y el plano basal de la máquina misma. La altura del eje H de una máquina
sin patas o de una máquina con pies levantados o cualquier máquina vertical, se debe tomar
como la altura de eje H de una máquina horizontal en el mismo marco básico. Cuando el soporte
es desconocido, la mitad del diámetro de máquina puede ser utilizada.
139
ANEXO IV
En las siguientes figuras podemos observar todos los eventos con todos los tipos de ventanas,
en el sensor knock, Acelerómetro, Micrófono.
Sensor Knock.
Carrera de Admisión tiene 239 grados.
140
Carrera de Compresión tiene 110 grados.
141
Carrera de Explosión tiene 172 grados.
142
Carrera de Escape tiene 214 grados.
143
Sensor de Vibración (Acelerómetro).
144
Carrera de Compresión tiene 110 grados.
145
Carrera de Explosión tiene 172 grados.
146
Carrera de Escape tiene 214 grados
.
147
Sensor de Ruido (Micrófono).
148
Carrera de Compresión tiene 110 grados.
149
Carrera de Explosión tiene 172 grados.
150
Carrera de Escape tiene 214 grados.
151
ANEXO V
En las siguientes figuras podemos visualizar todos los eventos con todos los tipos de ventanas,
en los sensores de knock, Acelerómetro, Micrófono.
Sensor Knock.
152
Carrera de Compresión tiene 86 grados
153
Carrera de Explosión tiene 148 grados
154
Carrera de Escape tiene 190 grados.
155
Sensor Vibración (Acelerómetro).
156
Carrera de Compresión tiene 86 grados
157
Carrera de Explosión tiene 148 grados
158
Carrera de Escape tiene 190 grados.
159
Sensor de Ruido (Micrófono).
160
Carrera de Compresión tiene 86 grados
161
Carrera de Explosión tiene 148 grados.
162
Carrera de Escape tiene 190 grados.
163
ANEXO VI
En las siguientes figuras podemos visualizar todos los eventos con todos los tipos de ventanas,
en los sensores de knock, Acelerómetro, Micrófono.
Sensor Knock.
164
Carrera de Compresión tiene 134 grados.
165
Carrera de Explosión tiene 196 grados.
166
Carrera de Escape tiene 238 grados.
167
Sensor Vibración (Acelerómetro).
168
Carrera de Compresión tiene 134 grados
.
169
Carrera de Explosión tiene 196 grados.
170
Carrera de Escape tiene 238 grados.
171
Sensor de Ruido (Micrófono).
172
Carrera de Compresión tiene 134 grados.
173
Carrera de Explosión tiene 196 grados.
174
Carrera de Escape tiene 238 grados.
175
ANEXO VII
Análisis Anova para la carrera de admisión, compresión, explosión y escape, para las señales
de Detonación, Vibración y Ruido.
ADMISIÓN (DETONACIÓN).
Método
Hipótesis nula Todas las medias son iguales
Hipótesis alterna Por lo menos una media es diferente
Nivel de significancia α = 0,05
176
Gráfica de caja de Media
-0,010
-0,015
Media
-0,020
-0,025
-0,030
Completa Gaussian Hamming Hanning Rectangular
Ventana
ADMISIÓN (DETONACIÓN).
Método
Hipótesis nula Todas las medias son iguales
Hipótesis alterna Por lo menos una media es diferente
Nivel de significancia α = 0,05
R-cuad. R-cuad.
S R-cuad. (ajustado) (pred)
0,0005918 28,64% 22,30% 11,90%
Medias
177
Gráfica de caja de Varianza
0,0030
0,0025
0,0020
Varianza
0,0015
0,0010
0,0005
0,0000
Completa Gaussian Hamming Hanning Rectangular
Ventana
COMPRESIÓN (RUIDO).
ANOVA unidireccional: Desviación vs. Ventana
Método
Hipótesis nula Todas las medias son iguales
Hipótesis alterna Por lo menos una media es diferente
Nivel de significancia α = 0,05
R-cuad. R-cuad.
S R-cuad. (ajustado) (pred)
0,0133487 70,34% 67,70% 63,38%
178
Gráfica de caja de Desviacion
0,12
0,11
0,10
0,09
Desviacion
0,08
0,07
0,06
0,05
0,04
0,03
Completa Gaussian Hamming Hanning Rectangular
Ventana
COMPRESIÓN (RUIDO)
ANOVA unidireccional: Mediana vs. Ventana
Método
R-cuad. R-cuad.
S R-cuad. (ajustado) (pred)
0,0180468 10,27% 2,29% 0,00%
179
Gráfica de caja de Mediana
0,075
0,050
0,025
Mediana
0,000
-0,025
-0,050
Completa Gaussian Hamming Hanning Rectangular
Ventana
COMBUSTIÓN (VIBRACIÓN).
Método
Hipótesis nula Todas las medias son iguales
Hipótesis alterna Por lo menos una media es diferente
Nivel de significancia α = 0,05
180
Factor Niveles Valores
Ventana 5 Completa. Gaussian. Hamming. Hanning. Rectangular
2,25
2,20
2,15
Maximo
2,10
2,05
2,00
1,95
Completa Gaussian Hamming Hanning Rectangular
Ventana
COMBUSTIÓN (VIBRACIÓN).
Método
Hipótesis nula Todas las medias son iguales
Hipótesis alterna Por lo menos una media es diferente
Nivel de significancia α = 0,05
181
Información del factor
Factor Niveles Valores
Ventana 5 Completa. Gaussian. Hamming. Hanning. Rectangular
R-cuad. R-cuad.
S R-cuad. (ajustado) (pred)
0,0531920 99,53% 99,48% 99,41%
1,6
1,4
1,2
1,0
Minimo
0,8
0,6
0,4
0,2
0,0
ESCAPE (DETONACIÓN).
Método
Hipótesis nula Todas las medias son iguales
Hipótesis alterna Por lo menos una media es diferente
Nivel de significancia α = 0,05
182
Información del factor
Factor Niveles Valores
Ventana 5 Completa. Gaussian. Hamming. Hanning. Rectangular
R-cuad. R-cuad.
S R-cuad. (ajustado) (pred)
0,0012141 95,87% 95,50% 94,90%
-0,012
-0,014
-0,016
-0,018
Media
-0,020
-0,022
-0,024
-0,026
-0,028
Completa Gaussian Hamming Hanning Rectangular
Ventana
183
ESCAPE (DETONACIÓN).
Método
Hipótesis nula Todas las medias son iguales
Hipótesis alterna Por lo menos una media es diferente
Nivel de significancia α = 0,05
R-cuad. R-cuad.
S R-cuad. (ajustado) (pred)
0,0005880 24,31% 17,58% 6,56%
184
Gráfica de caja de Varianza
0,0030
0,0025
0,0020
Varianza
0,0015
0,0010
0,0005
0,0000
Completa Gaussian Hamming Hanning Rectangular
Ventana
185
ANEXO VIII
Análisis Tukey para la carrera de admisión, compresión, explosión y escape, para las señales
de Detonación, Vibración y Ruido.
ADMISIÓN (DETONACIÓN).
Media vs. Ventana.
ADMISIÓN (DETONACIÓN).
Varianza vs. Ventana.
186
COMPRESIÓN (RUIDO).
Desviación vs. Ventana.
COMPRESIÓN (RUIDO).
Mediana vs. Ventana.
187
COMBUSTIÓN (VIBRACIÓN).
Máximo vs. Ventana.
COMBUSTIÓN (VIBRACIÓN).
Mínimo vs. Ventana.
188
ESCAPE (DETONACIÓN).
Media vs. Ventana.
ESCAPE (DETONACIÓN).
Varianza vs. Ventana.
189
ANEXO IX
Predicciones para analizar qué tipo y ancho de ventana sirve para cada evento en motores de
combustión interna.
VIBRACIÓN.
RUIDO.
kNOCK.
190
ANEXO X
Código en Labview para Identificación de Eventos Transitorios en motores de combustión
interna.
191
MANUAL DE USUARIO Y MANTENIMIENTO
Una vez que estén cargadas las señales se escribiré los grados en el número 1, esto indica en
qué grado del cigüeñal comienza el ventaneo, hay que tomar en cuenta que los grados tiene un
rango de 0° a 360°. En el número 2 se escribe que vuelta del cigüeñal se desea para el comienzo
del ventaneo. El número 3 sirve para visualizar el número de diente que se insertó en grados,
en el número 1.
Después se escribiré los grados en el número 4 esto indica en qué grado del cigüeñal termina
el ventaneo, de igual forma tiene un rango 0° a 360°. En el número 5 se escribe que vuelta del
cigüeñal se desea para terminar el ventaneo. El número 6 sirve para visualizar el número de
diente que se insertó en grados, en el número 4.
Esto quiere decir que siempre el número de diente de inicio del ventaneo, tiene que ser menor
al número de diente del final del ventaneo, en caso que se inserte mal los dientes saldrá el
siguiente mensaje.
192
EL número 7 sirve para capturar las señales en el momento que deseemos y el numero 8 sirve
para visualizar en que instante hace la captura de las señales, esto nos sirve para tomar todas
las muestras en un mismo tiempo.
Portada.
Señales.
Se ve todas las señales adquiridas (CMP, CKP, Inyección, Detonación, Vibración y Sonido)
en condiciones normales del motor es decir lo que está cargado en Read From File.
193
S. Cortada
En la ventana de señal cortada se ve todas las señales anteriores pero sincronizadas. Y también
se ve en que diente comienza el ventaneo y en que diente termina el ventaneo.
S. Analizada
En donde dice Channel se introduce el número de canal que se quiere analizar, el canal 0 es la
señal de detonación, el canal 1 es la señal de vibración y el canal 2 es la señal del sonido. En
la primera grafica se muestra la señal total que se analiza ya sea bien detonación, vibración o
sonido. En la segunda grafica se muestra la señal analizada, pero recortada en los límites que
se introducido al principio. Y en la tercera grafica se ve la señal analizada amplificada solo en
ese instante.
194
Características.
Para observar los datos estadísticos de las ventanas (completa, rectangular, hanning, hamming
y gaussian) se da un clik en el botón de TOMAR DATOS, e inmediatamente se guardaran los
datos y parara el programa. El botón que tiene una luz verde sirve para borrar los datos.
Ventaneo.
Aquí se ve la señal que escogimos a analizar pero echo un ventaneo Hanning, Hamming y
Gausian.
195
FFT.
FFT Fallas.
196