BRPI0803875B1 - método para estimar a posição angular do virabrequim na qual 50% da massa de combustível foi queimada em um cilindro de um motor de combustão interna com ignição espontânea da mistura - Google Patents
método para estimar a posição angular do virabrequim na qual 50% da massa de combustível foi queimada em um cilindro de um motor de combustão interna com ignição espontânea da mistura Download PDFInfo
- Publication number
- BRPI0803875B1 BRPI0803875B1 BRPI0803875-9A BRPI0803875A BRPI0803875B1 BR PI0803875 B1 BRPI0803875 B1 BR PI0803875B1 BR PI0803875 A BRPI0803875 A BR PI0803875A BR PI0803875 B1 BRPI0803875 B1 BR PI0803875B1
- Authority
- BR
- Brazil
- Prior art keywords
- tooth
- angular
- fact
- internal combustion
- delayed
- Prior art date
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1497—With detection of the mechanical response of the engine
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D35/00—Controlling engines, dependent on conditions exterior or interior to engines, not otherwise provided for
- F02D35/02—Controlling engines, dependent on conditions exterior or interior to engines, not otherwise provided for on interior conditions
- F02D35/028—Controlling engines, dependent on conditions exterior or interior to engines, not otherwise provided for on interior conditions by determining the combustion timing or phasing
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/0097—Electrical control of supply of combustible mixture or its constituents using means for generating speed signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/04—Testing internal-combustion engines
- G01M15/06—Testing internal-combustion engines by monitoring positions of pistons or cranks
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/24—Electrical control of supply of combustible mixture or its constituents characterised by the use of digital means
- F02D41/26—Electrical control of supply of combustible mixture or its constituents characterised by the use of digital means using computer, e.g. microprocessor
- F02D41/28—Interface circuits
- F02D2041/286—Interface circuits comprising means for signal processing
- F02D2041/288—Interface circuits comprising means for signal processing for performing a transformation into the frequency domain, e.g. Fourier transformation
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D2200/00—Input parameters for engine control
- F02D2200/02—Input parameters for engine control the parameters being related to the engine
- F02D2200/10—Parameters related to the engine output, e.g. engine torque or engine speed
- F02D2200/1002—Output torque
- F02D2200/1004—Estimation of the output torque
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
- Testing Of Engines (AREA)
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
- Electrical Control Of Ignition Timing (AREA)
Abstract
MÉTODO PARA ESTIMAR A POSIÇÃO ANGULAR DO VIRABREQUIM NA QUAL 50% DA MASSA DE COMBUSTÍVEL FOI QUEIMADA EM UM CILINDRO DE UM MOTOR DE COMBUSTÃO INTERNA COM IGNIÇÃO ESPONTÂNEA DA MISTURA. Trata-se de um método de estimação do ângulo do virabrequim no qual 50% da massa de combustível foi queimada em um cilindro (2) de um motor de combustão interna (1) com ignição espontânea da mistura, provido de um eixo de acionamento (4) acoplado a uma roda fônica (5) apresentando uma série N de dentes (6); o método compreendendo as etapas de: ler a passagem de cada dente (6) da roda fônica (5) diante de um sensor (7); determinar a velocidade angular ((Ômega)i) do eixo de acionamento (4) em cada evento de dente da roda tônica (5); determinar pelo menos um harmônico ((Ômega)n) do sinal de velocidade; determinar um modelo mecânico inverso (H(j(Ômega) = T / (Õmega)) da transmissão; determinar pelo menos um harmônico de torque (Tn) aplicando o modelo mecânico inverso (Hü(Ômega)) = T / (Ômega)) ao harmônico do sinal de velocidade ((Ômega)n); determinar uma função algébrica (Psi) que coloca um índice de combustão (MFB50) em relação com a fase (Arg(Tn)) do n-ésimo harmônico de torque (Tn); e determinar o índice de combustão (MFB50) aplicando (...).
Description
[001] A presente invenção se refere a um método para estimar a posição angular do virabrequim na qual 50% da massa de combustível foi queimada em um cilindro de um motor de combustão interna com ignição espontânea da mistura.
[002] Nos motores de combustão interna com ignição controlada da mistura (ou seja, motores de combustão interna que operam de acordo com o ciclo de "otto" e são alimentados com gasolina, metano, GLP, entre outros), utiliza-se há muitos anos como grandeza de controle o índice de combustão MFB50 (Fração de Massa de Combustível Queimada), que representa a posição angular do virabrequim na qual 50% da massa de combustível já foi queimada dentro de um cilindro. De fato, nos motores de combustão interna com ignição controlada da mistura, o instante de partida da combustão (isto é, o instante de centelhamento da centelha de ignição) é conhecido de forma exata e precisa, sendo escolhido de antemão pela estratégia de controle do motor por meio do controle de avanço de centelha.
[003] Já os motores de combustão interna com ignição espontânea da mistura, (quer dizer, motores de combustão interna operando conforme o ciclo "Diesel" e alimentados com óleo ou similar), o índice de combustão MFB50 não é utilizado como grandeza de controle face às consideráveis dificuldades em se estimar tal índice de combustão MFB50 de maneira eficiente (isto é, com precisão suficiente), eficaz (isto é, com rapidez e sem exigir muita capacidade de cálculo) e com boa relação custo-benefício (ou seja, sem exigir a instalação de componentes adicionais aos normalmente presentes).
[004] Nos dias atuais, utilizam-se sensores de pressão capazes de medir diretamente a pressão dentro dos cilindros para a determinação do índice de combustão MFB50; porém, tais sensores de pressão são extremamente caros e apresentamconfiabilidade limitada com o tempo, sendo assim adequados apenas a testes em laboratórios e impróprios para use extenso nos carros de produção convencionais.
[005] O documento US 5200899 descreve um método para a detecção de detonações (misfares) em motores a combustão interna. O método explora uma medição da velocidade angular do eixo árvore do motor em conjunto com o processamento de um sinal eletrônico híbrido; uma vez que o sinal da velocidade angular é condicionado para minimizar os efeitos de erros aleatórios e de distúrbios externos, sendo que quatro métodos eficientes de computação alternativas podem ser usados para extrair as informações pertinentes a produção de torque pelos cilindros individuais.
[006] O objetivo da presente invenção é o de oferecer um método para estimar com eficácia, eficiência e boa relação custo-beneficio a posição angular do virabrequim na qual 50% da massa de combustível foi queimada em um cilindro de um motor de combustão interna com ignição espontânea da mistura.
[007] De acordo com a presente invenção, propõe-se um método para estimar a posição angular do virabrequim na qual 50% da massa de combustível foi queimada em um cilindro de um motor de combustão interna com ignição espontânea da mistura conforme reivindicado nas reivindicações apensas.
[008] A presente invenção é descrita a seguir, com referência aos desenhos concomitantes, os quais ilustram uma concretização não restritiva da mesma, nos quais: - a figura 1 é uma vista diagramática de um motor de combustão interna provido de uma unidade de controle que implementa o método de estimação objeto da presente invenção; e- a figura 2 é um diagrama de blocos que ilustra a operação da unidade de controle na figura 1 para estimar o índice de combustão MFB50.
[009] Na figura 1, o número 1 indica, como um todo, um motor de combustãointerna alimentado com óleo (isto é, operando conforme o ciclo Diesel que contempla a ignição espontânea da mistura) e montado a bordo de um veículo de estrada (não ilustrado) que é provido de uma linha de transmissão (não ilustrada) para transmitir o torque gerado pelo motor de combustão interna 1 ao solo. O motor de combustão interna 1 compreende quatro cilindros 2, cada um dos quais acomoda um pistão correspondente 3 mecanicamente conectado, por meio de uma biela, a um eixo de acionamento 4 para transmitir a força gerada pela combustão dentro do cilindro 2 ao próprio eixo de transmissão 4.
[0010] Uma roda fônica 5 provida de uma série N (por exemplo, 60) de dentes 6 é encaixada no eixo de acionamento 4 e acoplada a um sensor 7, que é adaptado para detectar o tempo decorrido entre a passagem de dois dentes consecutivos 6. Os dentes 6 da roda fônica 5 são igualmente espaçados, com a exceção de um par de dentes 6 que são dispostos mutuamente a uma distância maior do que os demais para constituir uma singularidade que permite identificar de forma precisa cada um dos dentes 6 e, dessa forma, a posição angular da roda fônica 5 (isto é, do eixo de acionamento 4). O motor 1 compreende ainda uma unidade de controle 8 que é conectada a sensor 7 e adaptada para estimar, para cada cilindro 2, o índice de combustão MFB50 (50% de fração de massa de combustível queimada) correspondente, o qual representa a posição angular do virabrequim na qual 50% da massa de combustível foi queimada dentro do cilindro 2.
[0011] Com referência à figura 2, é descrito o modo usado pela unidade de controle 8 para estimar o índice de combustão MFB50 em cada cilindro 2 utilizando-se a informação fornecida pelo sensor 7 acoplado à roda fônica 5.
[0012] Conforme mencionado anteriormente, a roda fônica 5 é provida de uma série N (por exemplo, 60) de dentes 6, que são igualmente espaçados, com a exceção de um par de dentes 6 que são mutuamente dispostos a uma distancia maior do que os demais para constituir uma singularidade que permite identificar com precisão cada dente 6; dessa forma, cada dente 6 é identificado par um número sequencial característico i compreendido entre 1 e N.
[0013] O sensor 7 detecta o número i do dente atual 6 (i-ésimo dente), o tempo Ti do dente atual 6, isto é, o intervalo do tempo decorrido entre as detecções dos dois eventos angulares que caracterizam o dente atual 6, a velocidade angular bruta ωraw,i do dente atual 6, e o ângulo de acionamento α (isto é, a posição angular do eixo de acionamento 4 que está compreendia entre 0 e 4π radianos). A velocidade angular bruta (úraw,i do i-ésimo dente 6 é fornecida pela seguinte equação:Wraw,i - velocidade angular bruta do i-ésimo dente 6 [rad/s];Δαi - amplitude angular do i-ésimo dente 6 [rad];Ti - tempo do i-ésimo dente 6 [s].
[0014] As grandezas detectadas pelo sensor 7 e listadas acima são transmitidas a um bloco de compensação 9, que corrige cada velocidade angular bruta ωraw,i para compensar os erros provenientes das inevitáveis tolerâncias de fabricação dos dentes 6. Em princípio, todos os dentes 6 tem a mesma amplitude angular, mas na verdade, devido às inevitáveis tolerâncias de fabricação, cada dente 6 apresenta um erro de forma característico, isto é, apresenta uma amplitude angular característica mais ou menos diferente da amplitude angular teórica; como consequência, as velocidades angulares cürawj são afetadas por um ruído intenso que torna praticamente impossível seu uso sem uma compensação apropriada. Com o objetivo de compensar os erros de forma dos dentes 6, é armazenado um vetor (δpi, δp2,..., δpi,..., ÕPN) composto de N valores de compensação δpi, cada um dos quais está associado a um i-ésimo dente 6 correspondente e representa a diferença entre a amplitude angular real do i-ésimo dente 6 e a amplitude teórica do i-ésimo dente 6.
[0015] A velocidade angular compensada correspondente é calculada para cada dente 6 no bloco de compensação 9 utilizando-se a seguinte equação:WI - velocidade angular compensada do i-ésimo dente 6 [rad/s];Δαi - amplitude angular do i-ésimo dente 6 [rad]; δµi - valor de compensação δµi do i-ésimo dente 6 [rad];Ti - tempo do i-ésimo dente 6 [s].
[0016] O modo de calibração utilizado para determinar os valores de compensação δpi do dente 6 é descrito a seguir. O modo de calibração contempla calcular a velocidade angular bruta ωraw,i do i-ésimo dente 6 utilizando a equação supradescrita e calculando a velocidade angular média ωav,i por rotação do i-ésimo dente 6 por meio de um filtro de média móvel de N valores da sequência de velocidades angulares brutas ωraw,i; é importante observar que o filtro de média de N valores da sequência de velocidades angulares brutas ωraw,i introduz um retardo de fase igual a N/2 intervalos de dente nas velocidades angulares médias Em seguida, calcula-se o erro de velocidade retardado N/2 para o i-ésimo dente 6 utilizando-se a seguinte equação:Δωj-N/2 - erro de velocidade retardado do i-N/2-ésimo dente 6 [rad/s];ωav,i - velocidade angular média por rotação do i-ésimo dente 6 [rad/s] apresentando um retardo de fase igual a N/2 intervalos de dente;ωraw,i - N/2 - velocidade angular bruta ωraw,i do i-N/2-ésimo dente 6 [rad/s].
[0017] A sequência de velocidades angulares brutas ωraw,i é retardada propositalmente por um retardo de fase igual a N/2 intervalos de dente para se igualar ao retardo de fase introduzido pelo filtro de média móvel.
[0018] Em seguida, determina-se o erro angular retardado δαi - N/2 para o i-ésimo dente 6 utilizando-se a seguinte equação:δαi - N/2 - erro angular retardado do i-N/2-ésimo dente 6 [rad/s];Δωi - N/2 - erro angular retardado do i-N/2-ésimo dente 6 [rad/s];Ti - N/2 - tempo do i-N/2-ésimo dente 6 [s].
[0019] Finalmente, o vetor de erros angulares retardados δαi - N/2 pode ser utilizado tanto para determinar (no caso da primeira operação de calibração absoluta) quanto para atualizar o vetor de valores de compensação δpi. De acordo com uma concretizaçãopreferida, antes de efetivamente utilizar um vetor de erros angulares retardados verifica- se cada erro angular retardado δθi - N/2 do próprio vetor é inferior, como valor absoluto, a um limiar de aceitabilidade que indica a tolerância máxima possível para a forma do dente 6; somente se cada um dos erros angulares retardados δαi - N/2 do vetor for inferior, como valor absoluto, ao limiar de aceitabilidade e que o vetor de erros angulares ôa, - N/2 poderá ser efetivamente usado. Se um vetor de valores de compensação δpi for atualizado, utiliza-se, de preferência, uma média ponderada entre os valores de compensação anteriores δpi e os erros angulares retardados correspondentes δαi - N/2 que normalmente designam um peso maior aos valores de compensado anteriores δpi (por exemplo, quanto maior o valor, maior a confiabilidade dos valores de compensação anteriores δpi) Por exemplo, podemos utilizar a seguinte equação:δ|Ji,j - i-ésimo valor de compensação atualizado [rad];K - peso adimensional compreendido entre 0 e 1;δαi - i-ésimo erro angular [rad];δ|Jij -1 - i-1-ésimo valor de compensação anterior [rad].
[0020] De preferência, a operação de calibração descrita acima para determinação/atualização do vetor de valores de compensação ÕJJÍ é realizada apenas diante certas circunstâncias, como por exemplo, quando o motor de combustão interna 1 está na fase de corte e a velocidade de rotação do eixo de acionamento 4 está incluída em uma certa faixa de calibração, de modo que os torques inerciais recíprocos indicados sejam aproximadamente equilibrados com o motor de combustão interna 1 na fase de corte.
[0021] As velocidades angulares compensadas coi, calculadas pelo bloco de compensação 9, são fornecidas para um bloco de análise 10, no qual uma análise de frequência do tipo conhecido é realizada para as velocidades angulares compensadas ωi para determinar pelo menos um harmônico do sinal de velocidade Ωn. O valor genérico do harmônico Ωndo n-ésimo sinal de velocidade é sempre um número complexo,caracterizado por seu próprio módulo |Ωn| e sua própria fase Arg(Ωn).
[0022] Os harmônicos Ωndo sinal de velocidade calculados pelo bloco de análise 10 são fornecidos a um bloco de cálculo 11, no qual é implementado um algoritmo baseado no modelo mecânico inverso da transmissão que é simbolicamente indicado pela equação H(jω) = T/Ω e representa, no domínio de frequência ω, a relação entre as transformações de Fourier da velocidade angular Ω e do torque T do motor de combustão interna 1. A função H(jω) é uma função de valor complexo e permite a obtenção do n-ésimo harmônico de torque Tn por meio da equação:Tn - n-ésimo harmônico de torque;Ωn- n-ésimo sinal harmônico de velocidade.
[0023] Assim como o sinal harmônico Ωn, o harmônico do torque Tn é um número complexo caracterizado por seu próprio módulo |Tn| e por sua própria fase Arg(Tn). è importante observar que de modo a ser capaz de obter estimativas precisas, o modelo mecânico inverso da transmissão deve considerar os torques inerciais causados pelas massas em movimento recíproco (isto é, as massas dos pistões 3).
[0024] O módulo |Tn| do harmônico do torque Tn calculado pelo bloco de cálculo 11 é fornecido a um bloco de cálculo 12, no qual é implementada uma função algébrica r que coloca o torque médio indicado Tind em relação com o módulo |Tn| do n-ésimo harmônico do torque Tn de acordo com a seguinte equação:Tind - torque médio indicado Tind;|Tn| - n-ésimo módulo harmônico do torque Tn;P - conjunto de parâmetros de engenharia do motor de combustão interna 1.
[0025] Normalmente, os parâmetros de engenharia P do motor de combustão interna 1 compreendem a velocidade de rotação do eixo de acionamento 4 e a pressão de admissão (isto é, a pressão presente dentro de um tubo de admissão).
[0026] A fase Arg(Tn) dos harmônicos de torque Tn calculado pelo bloco de calculo11 é fornecida a um bloco de cálculo 13, no qual é implementada uma função algébrica ψ que coloca o índice de combustão MFB50 em relação com o a fase Arg(Tn) do n-ésimo harmônico de torque Tn de acordo com a seguinte equação:MFB50 - índice de combustão MFB50;(Tn) - fase do n-ésimo harmônico de torque Tn;P - conjunto de parâmetros de engenharia do motor de combustão interna 1.
[0027] A título de exemplo não restritivo, as equações do modelo mecânico direto da transmissão no domínio de frequência podem ser as seguintes:Lado do motor de torqueTorque inérciaI Torque do cilindro Ji - momento de inércia do motor de combustão interna 1;J2 - momento de inércia do veículo em relação ao eixo de acionamento 4;Keb - rigidez da linha de transmissão;peb - coeficiente de amortecimento da linha de transmissão;j - unidade imaginária;n - ordem do harmônico;ó)e - pulso;cúe - velocidade angular do eixo de acionamento 4;Te - torque gerado pelo motor de combustão interna 1; Tait,n - torque devido às massas em movimento recíproco;meq ■ massas em movimento recíproco;θait - fator de fase do torque inicial (movimento recíproco);Tind - torque médio indicado.
[0028] θ modelo mecânico inverso da transmissão (indicado simbolicamente antes pela equação H(jω) = T / Ω) é obtido a partir das três equações do modelo mecânico direto da transmissão ilustrada substituindo-se a primeira e a segunda equações na terceira equação e de modo a colocar o torque médio indicado Tind,n em relação com a velocidade angular ωe,n do eixo de acionamento 4; a relação entre o torque médio indicado Tind,n e a velocidade angular ωe,n do eixo de acionamento 4 fornece o modelo mecânico inverso da transmissão (a fórmula final do modelo mecânico inverso da transmissão não é ilustrada pois ela não acrescenta informações e é complexa de ser anotada). É importante observar que, na prática, realiza-se uma metodologia experimental para obter experimentalmente por pontos a função de transferência do modelo mecânico inverso da transmissão.
[0029] A título de exemplo não restritivo, a função algébrica ψ que coloca o índice de combustão MFB50 em relação com a fase Arg(Tn) do harmônico Tn do n-ésimo torque poderia ser a seguinte:MFB50 - índice de combustão;O - fase do n-ésimo harmônico de torque Tn;P - pressão de admissão; o - coeficientes.
[0030] Os valores dos coeficientes Q podem ser obtidos experimentalmente minimizando-se 0 erro quadrado médio entre uma serie de valores medidos do índice de combustão MFB50 e uma série correspondente de valores estimados do índice de combustão MFB50.
[0031] A título de exemplo não restritivo, a função algébrica T que coloca 0 torquemédio indicado ting em relação com o módulo |Tn| do harmônico Tn do n-ésimo torque poderia ser a seguinte:Tind - torque médio indicado Tnd;|Tn| - n-ésimo módulo harmônico do torque Tn;P - pressão de admissão;N - velocidade de rotação do eixo de acionamento 4;bi - coeficientes.
[0032] Os valores dos coeficientes bi podem ser obtidos experimentalmente minimizando-se o erro quadrado médio entre uma série de valores medidos do torque médio indicado Tind e uma série correspondente de valores estimados do torque médio indicado Tind.
[0033] O modo de estimação descrito acima do índice de combustão MFB50 em cada cilindro 2 apresenta muitas vantagens por ser eficiente (isto é, permite estimar o índice de combustão MFB50 com precisão considerável), eficaz (isto é, permite estimar o índice de combustão MFB50 com rapidez e sem exigir muita capacidade de cálculo) e ter boa relação custo-beneficio (ou seja, não exige a instalação de componentes adicionais aos normalmente presentes em um motor de combustão interna moderno operando de acordo com o ciclo "Diesel").
[0034] Em virtude da possibilidade de se estimar o índice de combustado MFB50 em cada cilindro 2, é possível controlar por realimentação o processo de combustão, em virtude do que se torna possível reduzir o consumo de combustível e as emissões poluentes, com o mesmo desempenho. Por exemplo, o índice de combustão MFB50 em cada cilindro 2 pode ser utilizado de forma vantajosa para controlar a forca sobre as bielas e o fenômeno de batimento do motor.
[0035] Além disso, o modo de estimação descrito acima do índice de combustão MFB50 em cada cilindro 2 permite obter o torque médio indicado Tind produzido por cada cilindro 2, que pode ser usado de forma eficaz para melhorar ainda mais o controle doprocesso de combustão.
Claims (15)
1. Método de estimação do ângulo do virabrequim no qual 50% da massa de combustível foi queimada em um cilindro (2) de um motor de combustão interna (1) com ignição espontânea da mistura provido de um eixo de acionamento (4) acoplado a uma roda fônica (5) apresentando uma série de dentes (6); o método de estimação compreendendo as etapas de:ler a passagem de cada dente (6) da roda fônica (5) diante de um sensor (7); e determinar a velocidade angular (ωi) do eixo de acionamento (4) em cada evento de dente da roda fônica (5),o método de estimação sendo caracterizado pelo fato de compreender as etapas adicionais de:- determinar, por meio de uma analise de frequência, as velocidades angulares (ω,) do eixo de acionamento (4), pelo menos um harmônico (Ωn) do sinal de velocidade apresentando seu próprio modulo (|Ωn|) e sua própria fase (Arg(Ωn));- determinar um modelo mecânico inverso (H(jω) = T/Ω) da transmissão que representa, no domínio de frequência (ω), a relação entre as transformadas de Fourier da velocidade angular (Ω) e do torque (T) do motor de combustão interna (1);- determinar pelo menos um harmônico de torque (Tn) apresentando seu próprio módulo (|Tn|) e sua própria fase (Arg(Tn)), aplicando o modelo mecânico inverso (H(jω) = T/ Ω) da transmissão ao harmônico do sinal de velocidade (Ωn);- determinar uma primeira função algébrica (ψ) que coloca um índice de combustão (MFB50) em relação com a fase (Arg(Tn)) do n-ésimo harmônico de torque (Tn); e- determinar o índice de combustão (MFB50) aplicando a primeira função algébrica (ψ) ao n-ésimo harmônico de torque (Tn).
2. Método, de acordo com a reivindicação 1, caracterizado pelo fato de que o modelo mecânico inverso (H(jω) = T/Ω) da transmissão leva em consideração os torques inerciais devido às massas em movimento recíproco do motor de combustão interna (1).
3. Método, de acordo com a reivindicação 2, caracterizado pelo fato de que a primeira função algébrica (ψ) depende de alguns parâmetros de engenharia (P) do motor da rnmhi ictãn intorna (1
4. Método, de acordo com a reivindicação 3, caracterizado pelo fato de que os parâmetros de engenharia (P) do motor de combustão interna (1) compreendem a velocidade de rotação do eixo de acionamento (4) e a pressão de admissão.
5. Método, de acordo com qualquer uma das reivindicações 1 a 4, caracterizado por compreender as etapas adicionais de:determinar uma segunda função algébrica (I") que coloca um torque médio indicado (Tind) em relação com o módulo (|Tn|) do n-ésimo harmônico de torque (Tn); edeterminar o torque médio indicado (Tind) aplicando a segunda função algébrica (I") ao n- ésimo harmônico de torque (Tn).
6. Método, de acordo com a reivindicação 5, caracterizado pelo fato de que a segunda função algébrica (I") depende de alguns parâmetros de engenharia (P) do motor de combustão interna (1).
7. Método, de acordo com a reivindicação 6, caracterizado pelo fato de que os parâmetros de engenharia (P) do motor de combustão interna (1) compreendem a velocidade de rotação do eixo de acionamento (4) e a pressão de admissão.
8. Método, de acordo com qualquer uma das reivindicações 1 a 7, caracterizado pelo fato de que a etapa de determinar a velocidade angular (ωi) do eixo de acionamento (4) em cada evento de dente da roda fônica (5) compreende as etapas adicionais de: detectar a amplitude angular (Δαi) para cada dente 6;detectar o tempo (T) de cada dente (6); edeterminar a velocidade angular bruta (ωraw,i) de cada dente (6) de acordo com a amplitude angular correspondente (Δαi) e o tempo (T).
9. Método, de acordo com a reivindicação 8, caracterizado por adicionalmente compreender as etapas de:determinar, para cada dente (6), urn valor de compensacao que representa a diferença entre a amplitude angular real do dente (6) e a amplitude angular teórica do dente (6); e determinar a velocidade angular (ωQ de cada dente (6) corrigindo a velocidade angular bruta por (ωraw,i) meio do valor de compensação correspondente (õpi).
10. Método, de acordo com a reivindicação 9, caracterizado pelo fato de que a vplnridadp snnular G.iA dp rada dpnt-p P dpt-prminada dp arnrdn rnm a cpniiint-pequação:ωi - velocidade angular compensada do i-ésimo dente (6) [rad/s];Δαi - amplitude angular do i-ésimo dente (6) [rad];δpij - valor de compensação δpij do i-ésimo dente (6) [s].
11. Método, de acordo com a reivindicação 9 ou 10, caracterizado pelo fato de que a etapa de determinar o valor de compensação (δpij) de cada dente (6) compreende as etapas adicionais de:calcular a velocidade angular bruta (ωraw,i) do dente (6);calcular uma velocidade angular média (ωav,i) por rotação do dente (6) por meio de um filtro de média móvel de N valores da sequência de velocidades angulares brutas (ωraw,i);calcular um erro de velocidade retardado (Δωi - N/2) do dente (6) como a diferença entre a velocidade angular média (ωav,i) por rotação e a velocidade angular bruta retardada (lüraw,i) de um retardo de fase igual a N/2 intervalos de dente;determinar um erro angular retardado (δαi - N/2) do dente (6) multiplicando 0 erro de velocidade retardado (Δωi - N/2) do dente (6) pelo tempo (T,) do dente (6) retardado por um retardo de fase igual a N/2 intervalos de dente; eutilizar 0 erro angular retardado (δαi - N/2) tanto para determinar como para atualizar 0 valor de compensação (δpij) do dente (6).
12. Método, de acordo com a reivindicação 11, caracterizado pelo fato de que 0 erro de velocidade retardado (Δω, - N/2) do dente (6) é calculado utilizando-se a seguinte equação:Δωi - N/2 - velocidade angular retardada do i-N/2-ésimo dente (6) [rad/s];ωav,i - velocidade angular média por rotação do i-ésimo dente (6) [rad/s] apresentando um retardo de fase igual a N/2 intervalos de dente;ωraw,i-N/2 - velocidade angular bruta ωraw,i do i-N/2-ésimo dente (6) [rad/s].
14. Método, de acordo com a reivindicação 11, 12 ou 13, caracterizado pelo fato de que um vetor de N erros angulares retardados (δαi - N/2) é utilizado tanto para determinar quanto para atualizar um vetor de valores de compensação (δpi) apenas se cada erro angular retardado (δαi - N/2) do próprio vetor for inferior, como valor absoluto, a um limiar de aceitabilidade que indica a tolerância máxima possível sobre a forma do dente (6).
15. Método, de acordo com qualquer uma das reivindicações 9 a 14, caracterizado pelo fato de que os valores de compensação (δpi) dos dentes (6) são determinados ou atualizados somente quando 0 motor de combustão interna (1) está na fase de corte e a velocidade de rotação do eixo de acionamento (4) está compreendida em uma certa faixa de calibração, de modo que os torques inerciais indicados recíprocos sejam aproximadamente equilibrados com 0 motor de combustão interna (1) na fase de corte.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
BR122018074529-0A BR122018074529B1 (pt) | 2007-07-25 | 2008-07-24 | método para determinar a velocidade angular média do virabrequim de um motor a combustão interna |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP07425457.4 | 2007-07-25 | ||
EP07425457A EP2022967B1 (en) | 2007-07-25 | 2007-07-25 | Method for estimating the crank angle at which the 50% mass fraction burn time occurs in an internal combustion engine |
Publications (2)
Publication Number | Publication Date |
---|---|
BRPI0803875A2 BRPI0803875A2 (pt) | 2009-06-16 |
BRPI0803875B1 true BRPI0803875B1 (pt) | 2020-11-17 |
Family
ID=38894116
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
BR122018074529-0A BR122018074529B1 (pt) | 2007-07-25 | 2008-07-24 | método para determinar a velocidade angular média do virabrequim de um motor a combustão interna |
BRPI0803875-9A BRPI0803875B1 (pt) | 2007-07-25 | 2008-07-24 | método para estimar a posição angular do virabrequim na qual 50% da massa de combustível foi queimada em um cilindro de um motor de combustão interna com ignição espontânea da mistura |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
BR122018074529-0A BR122018074529B1 (pt) | 2007-07-25 | 2008-07-24 | método para determinar a velocidade angular média do virabrequim de um motor a combustão interna |
Country Status (6)
Country | Link |
---|---|
US (2) | US7827857B2 (pt) |
EP (2) | EP2022967B1 (pt) |
CN (2) | CN101408140B (pt) |
AT (2) | ATE481561T1 (pt) |
BR (2) | BR122018074529B1 (pt) |
DE (2) | DE602007009330D1 (pt) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ITBO20100562A1 (it) * | 2010-09-16 | 2012-03-17 | Magneti Marelli Spa | Metodo di stima di un indice di combustione e/o della coppia indicata in un cilindro di un motore a combustione interna |
KR101189486B1 (ko) * | 2010-09-30 | 2012-10-12 | 한양대학교 산학협력단 | 엔진의 연소 위상 검출 방법 |
DE102011090077A1 (de) * | 2011-12-29 | 2013-07-04 | Robert Bosch Gmbh | Verfahren zum Bestimmen einer mittleren Drehzahl einer sich drehenden Antriebswelle einer Brennkraftmaschine |
ITBO20120023A1 (it) * | 2012-01-19 | 2013-07-20 | Magneti Marelli Spa | Metodo di stima dello squilibrio di coppia indicata erogata dai cilindri di un motore a combustione interna |
CN105283768A (zh) | 2013-06-12 | 2016-01-27 | 瓦锡兰芬兰有限公司 | 确定发动机的角速度 |
KR102373491B1 (ko) * | 2015-07-15 | 2022-03-11 | 삼성전자주식회사 | 회전체의 회전 인식 방법 및 그 방법을 처리하는 전자 장치 |
ITUB20154998A1 (it) * | 2015-11-03 | 2017-05-03 | Magneti Marelli Spa | Metodo di stima dell'indice mfb50 di combustione e della coppia istantanea generata dai cilindri di un motore a combustione interna |
GB2555436B (en) * | 2016-10-27 | 2019-05-08 | Delphi Tech Ip Ltd | Method to determine tooth error in period of one or more teeth of a rotating toothed wheel |
FR3086387B1 (fr) * | 2018-09-24 | 2020-08-28 | Continental Automotive France | Procede de determination de la position d'un vilebrequin de vehicule automobile |
DE102019207252A1 (de) * | 2018-11-14 | 2020-05-14 | Vitesco Technologies GmbH | Erfassung von zylinderindividuellen Brennverlaufsparameterwerten für einen Verbrennungsmotor |
CN109630289B (zh) * | 2018-12-13 | 2020-07-07 | 清华大学 | 基于燃烧模型重构缸内压力的发动机瞬时转矩估计方法 |
CN109991438B (zh) * | 2019-03-26 | 2021-12-07 | 武汉理工大学 | 一种汽车齿圈角速度测量误差消除方法与装置 |
FR3129983B1 (fr) | 2021-12-08 | 2023-10-27 | Vitesco Technologies | procédé et un dispositif de contrôle moteur |
EP4299889A1 (en) | 2022-06-30 | 2024-01-03 | Marelli Europe S.p.A. | Method to estimate the maximum pressure inside a combustion chamber of a cylinder of an internal combustion engine |
CN115935840B (zh) * | 2022-12-02 | 2023-08-11 | 哈尔滨工程大学 | 一种基于倒拖缸压和扭矩信号的缸内压力重构方法 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3122376A1 (de) * | 1981-06-05 | 1982-12-23 | Robert Bosch Gmbh, 7000 Stuttgart | Vorrichtung zur erfassung der drehzahl von rotierenden teilen |
US5200899A (en) * | 1990-04-20 | 1993-04-06 | Regents Of The University Of Michigan | Method and system for detecting the misfire of an internal combustion engine utilizing angular velocity fluctuations |
US6142026A (en) * | 1994-06-06 | 2000-11-07 | Toyota Jidosha Kabushiki Kaisha | Wheel information estimating apparatus |
DE19540674C2 (de) * | 1995-10-31 | 1999-01-28 | Siemens Ag | Adaptionsverfahren zur Korrektur von Toleranzen eines Geberrades |
US5977765A (en) * | 1997-12-05 | 1999-11-02 | Ford Global Technologies, Inc. | Speed, direction, and acceleration sensor for a rotating shaft having a rotor with teeth having unequal spacing |
FR2793558B1 (fr) * | 1999-05-10 | 2001-07-06 | Renault | Procede et dispositif de mesure du couple d'un moteur a combustion interne |
US6363318B1 (en) * | 2000-06-21 | 2002-03-26 | Cummins Engine Company, Inc. | Method to compensate errors in engine speed measurement |
WO2002071308A1 (en) * | 2001-03-05 | 2002-09-12 | The Ohio State University | Engine control using torque estimation |
US6752009B2 (en) * | 2001-08-03 | 2004-06-22 | General Motors Corporation | Encoded crank position sensor |
ATE376620T1 (de) * | 2002-01-24 | 2007-11-15 | Fraunhofer Ges Forschung | Verfahren zur bestimmung und kompenstion von geometriefehlern eines drehgeberrades |
ITBO20030529A1 (it) * | 2003-09-15 | 2005-03-16 | Magneti Marelli Powertrain Spa | Metodo per rilevare la accelerazione angolare di un |
DE102004048133A1 (de) * | 2004-10-02 | 2006-04-06 | Robert Bosch Gmbh | Verfahren zur Messung der Drehzahl einer Kurbelwelle |
FI121150B (fi) * | 2005-11-30 | 2010-07-30 | Waertsilae Finland Oy | Mäntäpolttomoottorisysteemin laitteisto ja menetelmä tunnistamaan epäyhtenäinen sylinteriteho-osuus |
JP2007198366A (ja) * | 2005-12-26 | 2007-08-09 | Denso Corp | 内燃機関の制御システム及びタイミングロータ |
US7400993B2 (en) * | 2006-10-31 | 2008-07-15 | General Electric Company | Method and system for accuracy of speed and acceleration measurements on rotating machinery with a speed wheel |
US7814780B2 (en) * | 2007-04-09 | 2010-10-19 | Bg Soflex Llc | Engine position tracking for internal combustion engines |
-
2007
- 2007-07-25 DE DE602007009330T patent/DE602007009330D1/de active Active
- 2007-07-25 AT AT09168056T patent/ATE481561T1/de not_active IP Right Cessation
- 2007-07-25 DE DE602007009520T patent/DE602007009520D1/de active Active
- 2007-07-25 EP EP07425457A patent/EP2022967B1/en active Active
- 2007-07-25 AT AT07425457T patent/ATE483103T1/de not_active IP Right Cessation
- 2007-07-25 EP EP09168056A patent/EP2136061B1/en active Active
-
2008
- 2008-07-23 US US12/177,922 patent/US7827857B2/en active Active
- 2008-07-24 BR BR122018074529-0A patent/BR122018074529B1/pt not_active IP Right Cessation
- 2008-07-24 BR BRPI0803875-9A patent/BRPI0803875B1/pt active IP Right Grant
- 2008-07-25 CN CN200810134756XA patent/CN101408140B/zh active Active
- 2008-07-25 CN CN201210393607.1A patent/CN103147867B/zh active Active
-
2010
- 2010-10-05 US US12/897,968 patent/US8028568B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US20090095060A1 (en) | 2009-04-16 |
EP2022967A1 (en) | 2009-02-11 |
EP2136061B1 (en) | 2010-09-15 |
US8028568B2 (en) | 2011-10-04 |
EP2136061A9 (en) | 2010-05-12 |
EP2136061A1 (en) | 2009-12-23 |
DE602007009330D1 (de) | 2010-10-28 |
ATE481561T1 (de) | 2010-10-15 |
BR122018074529B1 (pt) | 2020-11-17 |
EP2022967B1 (en) | 2010-09-29 |
ATE483103T1 (de) | 2010-10-15 |
CN101408140A (zh) | 2009-04-15 |
US7827857B2 (en) | 2010-11-09 |
BRPI0803875A2 (pt) | 2009-06-16 |
US20110029267A1 (en) | 2011-02-03 |
DE602007009520D1 (de) | 2010-11-11 |
CN103147867B (zh) | 2015-10-21 |
CN101408140B (zh) | 2013-06-12 |
CN103147867A (zh) | 2013-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
BRPI0803875B1 (pt) | método para estimar a posição angular do virabrequim na qual 50% da massa de combustível foi queimada em um cilindro de um motor de combustão interna com ignição espontânea da mistura | |
US7623955B1 (en) | Method for estimation of indicated mean effective pressure for individual cylinders from crankshaft acceleration | |
JP4276685B2 (ja) | 内燃機関の異常検出装置および異常検出方法 | |
US10215111B2 (en) | Controller and control method for internal combustion engine | |
JP5026334B2 (ja) | 角速度及び角加速度算出装置、トルク推定装置、燃焼状態推定装置 | |
US8256278B2 (en) | Engine misfire detection systems and methods using discrete fourier transform approximation | |
US8532908B2 (en) | System and method for estimating indicated mean effective pressure of cylinders in an engine | |
JPH0392570A (ja) | 内燃機関の燃料量及び点火時期制御方法および装置 | |
JP6822820B2 (ja) | 内燃機関のシリンダーによって生じさせられるmfb50燃焼指数と瞬時トルクとを推定する方法 | |
EP2431595B1 (en) | Method of estimating a combustion index and/or the indicated torque in a cylinder of a four stroke internal combustion engine with spontaneous mixture ignition | |
TWI418702B (zh) | Egr率推測檢測裝置 | |
JP5131230B2 (ja) | 燃焼判定装置 | |
JP5508631B2 (ja) | 内燃機関の制御装置 | |
BR102013027800B1 (pt) | método para reconhecimento do tipo de combustível efetivamente utilizado em um motor de combustão interna | |
TWI221504B (en) | Engine control device | |
BR112012005586B1 (pt) | método para predição do regime de rotação de um virabraquim de motor em fase final de rotação e uso do referido método para predição do regime de rotação de um virabrequim de motor em fase final de rotação | |
JP2009174483A (ja) | 内燃機関の制御装置 | |
CN116324152A (zh) | 曲柄角传感器控制装置和内燃机控制装置 | |
JP2004218605A (ja) | 内燃機関の失火検出装置 | |
JP5876297B2 (ja) | 車両制御装置 | |
JP2005307945A (ja) | 内燃機関の失火判定装置 | |
ITBO20120023A1 (it) | Metodo di stima dello squilibrio di coppia indicata erogata dai cilindri di un motore a combustione interna | |
Ponti | In-Cylinder Pressure Measurement: Requirements for On-Board Engine Control | |
JPH1123424A (ja) | エンジンの発生トルク推定方法 | |
JP2006152859A (ja) | 内燃機関の失火検出装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
B03A | Publication of a patent application or of a certificate of addition of invention [chapter 3.1 patent gazette] | ||
B06T | Formal requirements before examination [chapter 6.20 patent gazette] | ||
B06F | Objections, documents and/or translations needed after an examination request according [chapter 6.6 patent gazette] | ||
B09A | Decision: intention to grant [chapter 9.1 patent gazette] | ||
B16A | Patent or certificate of addition of invention granted [chapter 16.1 patent gazette] |
Free format text: PRAZO DE VALIDADE: 10 (DEZ) ANOS CONTADOS A PARTIR DE 17/11/2020, OBSERVADAS AS CONDICOES LEGAIS. |