Um modelo inspirado em circuitos para sistemas de manufatura

Resumo

Este artigo apresenta um modelo para representar e analisar sistemas de manufatura por analogia com circuitos elétricos. A abordagem, chamada de modelo de circuito de manufatura, oferece uma perspectiva baseada em física que complementa os métodos tradicionais. Ao mapear elementos de manufatura em componentes de circuito, o modelo permite aplicar a teoria de circuitos para quantificar a dinâmica do sistema, dimensionar pulmões e traduzir mudança operacional em custo. Os conceitos centrais incluem o fator de potência de manufatura (MPF), o dimensionamento de pulmões pelo amortecimento e a análise no domínio da frequência do fluxo de produção. Cada parâmetro do modelo remete a uma medição ou a um valor de política escolhido, de modo que o modelo não carrega constantes de ajuste livres.

1. Introdução

Sistemas de manufatura são redes complexas de processos interconectados, fluxos de material e trocas de informação. Analisar e otimizar esses sistemas é essencial para melhorar a eficiência, reduzir o desperdício e aumentar a capacidade de resposta. Abordagens tradicionais como teoria de filas, simulação e modelagem de eventos discretos trouxeram contribuições valiosas, mas muitas vezes carecem de um modelo coeso, baseado em física, que capture a interação dinâmica dos componentes do sistema.

Este artigo traça uma analogia direta entre sistemas de manufatura e circuitos elétricos. O modelo de circuito de manufatura resultante oferece uma forma compacta e mensurável de raciocinar sobre como uma linha de produção responde a mudanças na demanda.

2. Revisão da literatura

Os métodos existentes de análise de sistemas de manufatura incluem:

• Teoria de filas: modela filas de espera e congestionamento, mas foca em médias de regime permanente em vez da resposta transitória.
• Simulação de eventos discretos: poderosa para interações complexas, mas carece de um modelo geral em forma fechada e pode ser cara em termos computacionais.
• Manufatura enxuta: enfatiza a redução de desperdício e o fluxo, mas carece de uma abordagem quantitativa da dinâmica no nível do sistema.
• Dinâmica de sistemas: usa laços de realimentação e atrasos para modelar o comportamento; o modelo de circuito é uma especialização dessa ideia com elementos nomeados e mensuráveis.
• Factory physics: fornece relações quantitativas como a lei de Little, mas para antes de uma analogia de circuito no nível dos componentes.

O modelo de circuito de manufatura se apoia nessas bases com um arcabouço intuitivo, matematicamente explícito e ancorado na medição.

3. O modelo de circuito de manufatura

Unidades de base ao longo do texto: peças $\text{P}$ e segundos $\text{s}$.

Tabela 1: analogias circuito-manufatura

A demanda entra como uma fonte de corrente, não de tensão: o cliente define uma taxa, e ou a linha a atende, ou o pulmão e a carteira de pedidos absorvem a diferença. Como a demanda é uma fonte de corrente, as grandezas no domínio da frequência a seguir vêm das funções de transferência das equações de estado da seção 3.3, e não de um triângulo de impedância RLC série: neste modelo, $R_{\text{eff}}$ está no termo de rigidez, não no dissipativo.

3.1 Os três elementos, ancorados na medição

• Resistência $R_{\text{eff}} = t_c / (m\,D)$ em s/P. O tempo de ciclo efetivo de uma estação: tempo de ciclo nominal $t_c$ dividido pelo número de máquinas em paralelo $m$ e pela fração operacional $D$ (fração do tempo de calendário em que a estação pode operar). Seu inverso $\mu = 1 / R_{\text{eff}}$ em P/s é a taxa máxima da estação.
• Capacitância $C = Q_{\max}$ em P. O pulmão físico entre duas estações. O potencial do pulmão é a fração de preenchimento adimensional $V_C = Q / Q_{\max}$, então $C = Q / V_C$ tem unidade de peças. A relação do capacitor é exata: $\dot Q = I_{\text{in}} - I_{\text{out}}$.
• Indutância $L = R_{\text{eff}}\,\tau_L$ em s²/P. Uma estação não pode mudar sua vazão de saída instantaneamente: o material em processo precisa ser reposicionado, o ferramental trocado, a mão de obra e o balanceamento da linha ajustados. O atraso $\tau_L$ em s é o tempo medido para a vazão atingir 63% de um novo nível após um degrau de demanda sustentado. A quantidade armazenada $\tfrac{1}{2} L I^2 = \tfrac{1}{2} R_{\text{eff}}\,\tau_L\,I^2$ tem unidade de peças e representa o trabalho em processo já comprometido, que se conclui mesmo que os lançamentos parem.

3.2 Leis de circuito, com saturação

Lei da vazão (uma resistência que satura). Enquanto o pulmão tem folga, a vazão é definida pela fração de preenchimento e pelo tempo de ciclo efetivo:

Unidades: adimensional dividido por s/P dá P/s.

O que a lei da vazão representa. É uma política de liberação, não a física da máquina. Uma estação deixada por conta própria opera na vazão plena $\mu = 1/R_{\text{eff}}$ sempre que há uma peça esperando, o que tornaria a saída bruta igual a $\mu$ para qualquer $Q > 0$, e não proporcional a $Q$. Aqui, $I = V_C / R_{\text{eff}}$ é lida como uma regra de puxada proporcional ao material em processo: a estação está autorizada a puxar e liberar trabalho a uma taxa proporcional a quão cheio está seu pulmão de entrada, a forma linear do princípio CONWIP e kanban de que uma etapa não deve ultrapassar seu sinal a jusante. É isso que torna legítima a forma linear em $Q$: é uma lei de controle escolhida, não uma afirmação sobre o equipamento. O ganho de puxada dessa política é fixado em $1/R_{\text{eff}}$ para que a taxa autorizada atinja exatamente a capacidade real quando o pulmão está cheio, e o restante do artigo usa esse ajuste, de modo que $\mu = 1/R_{\text{eff}}$ em todo o texto e a cláusula de saturação abaixo ainda vincula à capacidade real. Operar com um ganho deliberadamente reduzido $g < 1/R_{\text{eff}}$ é possível, mas é uma decisão de projeto separada que baixa a capacidade efetiva para $g$; não é considerada aqui. Assim, $\omega_n$, $\zeta$, $\phi$ e $A_Q$ abaixo descrevem uma linha controlada, não uma em malha aberta.

Quando a demanda excede a capacidade, $\lambda \ge \mu = 1 / R_{\text{eff}}$, a fração de preenchimento trava em 1, a vazão satura em $\mu$, e a demanda não atendida se acumula como carteira de pedidos $B$ com $\dot B = \lambda - \mu$. Uma aproximação suave na transição é a forma de condutâncias em paralelo

com os três termos em s/P. Isso substitui a lei linear $V = I R$, que não fecha dimensionalmente e não captura que a vazão satura na capacidade.

Conservação do fluxo (lei dos nós de Kirchhoff). Em cada nó, as unidades se conservam:

3.3 Dinâmica de segunda ordem

Acopla-se o pulmão (capacitor) ao atraso de vazão (indutor). Com $Q$ o nível do pulmão e $I$ a saída da estação:

Eliminando $I$, obtém-se um sistema amortecido de segunda ordem no nível do pulmão:

com frequência natural e razão de amortecimento

Ambas as expressões conferem em unidades. Cada parâmetro é fixado antes de o modelo rodar: $R_{\text{eff}}$ por um estudo de tempo de ciclo (e usado como ganho de puxada, 3.2), $C$ pelo pulmão, $\tau_L$ por um ensaio de degrau na saída. Não há constantes livres a ajustar.

4. Conceitos centrais e aplicações

4.1 Fator de potência de manufatura (MPF)

Transformando as equações de estado da seção 3.3, com a demanda $\lambda$ como fonte de corrente e a saída da estação $I$ como resposta, obtém-se a função de transferência demanda-para-vazão

um passa-baixas de segunda ordem com os $\omega_n$ e $\zeta$ da seção 3.3. Numa frequência de demanda $\omega$, a vazão atrasa em relação à demanda pelo ângulo de fase

Numerador e denominador têm ambos unidade de $\text{s}^{-2}$, então $\tan\phi$ é adimensional. Define-se

Interpretação. Decompõe-se a resposta de vazão na frequência $\omega$ numa parte em fase com a demanda e numa parte em quadratura. A parte em fase, proporcional a $\cos\phi$, é o fluxo que atende à demanda. A parte em quadratura, proporcional a $\sin\phi$, enche o pulmão em metade do ciclo e o esvazia na outra metade, com contribuição líquida nula à vazão, mas ainda consumindo mão de obra, tempo de máquina e custo de manutenção de estoque. O MPF é a fração da resposta de vazão que acompanha a demanda, numa escala de 0 a 1: vale 1 quando a saída está em fase e cai a 0 em $\omega = \omega_n$, onde a resposta é pura troca de pulmão. Para $\omega \ge \omega_n$, o $\cos\phi$ bruto é negativo, ou seja, a saída se move contra a demanda, e o MPF é limitado a 0; a grandeza relevante nesse regime é a própria defasagem $\phi$, não seu cosseno. O MPF chegando a 0 em $\omega_n$ é a forma no domínio da frequência do alerta da seção 4.2: não deixe a variação dominante da demanda ficar na frequência natural.

A alavanca é $\tau_L$. Para $\omega < \omega_n$, reduzir o atraso de ajuste da vazão $\tau_L$ diminui $\tan\phi$ e aumenta o MPF, e ao mesmo tempo aumenta $\omega_n$ e $\zeta$ (4.2). Redução de setup (SMED), mão de obra flexível e lotes de transferência menores reduzem todos $\tau_L$. O atraso tem um piso: quando $\tau_L \to 0$, o modelo se reduz ao atraso de primeira ordem do pulmão $\phi \to \arctan(\omega\,Q_{\max}R_{\text{eff}})$, fixado pela constante de tempo do pulmão $Q_{\max}R_{\text{eff}}$. Um pulmão maior aumenta o amortecimento (4.2), mas também aumenta esse atraso residual, então o tamanho do pulmão troca acompanhamento da demanda por estabilidade em vez de melhorar os dois.

Oscilação do material em processo. A resposta do pulmão à demanda decorre de $\dot Q = \lambda - I$:

com $A_Q$ em segundos. Uma componente de demanda de amplitude $\Delta\lambda$ na frequência $\omega$ produz uma oscilação do pulmão de amplitude $\Delta Q = \Delta\lambda\,A_Q(\omega)$. Em baixa frequência, $A_Q \to Q_{\max}R_{\text{eff}}$, então uma mudança lenta de demanda de tamanho $\Delta\lambda$ desloca a fração de preenchimento em $\Delta\lambda/\mu$, coincidindo com $V_C^{\ast} = \lambda/\mu$ da seção 4.3. Perto de $\omega_n$, $A_Q$ sobe acima desse piso sempre que $\zeta < \tfrac{1}{2}\sqrt{1 + \sqrt{2}} \approx 0{,}777$: é o efeito chicote (bullwhip). O limiar não é $1/\sqrt{2}$, porque o zero em $s = -1/\tau_L$ em $Q/\lambda$ o eleva; ver 4.2. Atender à regra de pulmão da seção 4.2 mantém $A_Q$ no piso de baixa frequência ou abaixo dele em toda a banda de demanda, e bem acima de $\omega_n$ o pulmão filtra a oscilação, então $A_Q$ cai abaixo do piso.

4.2 Dimensionamento de pulmões pelo amortecimento

Uma versão anterior deste modelo prescrevia operar na ressonância, $\omega_n = \omega_{\text{takt}}$. Isso está errado. Perto de $\omega_n$, uma linha pouco amortecida amplifica as oscilações do pulmão (o modo chicote ou de caça), um comportamento a suprimir, não a buscar.

O objetivo de projeto é uma resposta demanda-para-material-em-processo $A_Q(\omega)$ (4.1) sem pico acima de seu valor de baixa frequência. A função de transferência do pulmão $Q/\lambda$ carrega um zero em $s = -1/\tau_L$ que eleva seu lado de alta frequência, então a condição de resposta plana não é $\zeta \ge 1/\sqrt{2}$. Escrevendo $|Q/\lambda|^2$ na forma $(x + 4\zeta^2)/\bigl((1-x)^2 + 4\zeta^2 x\bigr)$ com $x = (\omega/\omega_n)^2$, sua inclinação em $x = 0$ é $1 + 8\zeta^2 - 16\zeta^4$, que permanece positiva (a resposta sobe acima de seu valor DC) até

Unidades: s dividido por s/P dá P. O limiar simples $\zeta \ge 1/\sqrt{2}$ ainda se aplica, mas à resposta demanda-para-vazão $I/\lambda$, que é um passa-baixas puro de segunda ordem sem zero: dimensione por ele quando o objetivo for uma vazão de saída suave, e pela regra em $(1 + \sqrt{2})$ ao dimensionar o pulmão contra a oscilação de material em processo. Para um assentamento monótono do material em processo, sem qualquer sobressinal, exija $\zeta \ge 1$:

$\omega_n$ continua útil: é a frequência de variação da demanda que a linha rejeita pior. Manter $\zeta$ acima de $\tfrac{1}{2}\sqrt{1 + \sqrt{2}}$ garante que as flutuações do pulmão nessa banda não sejam amplificadas.

4.3 Tempos de enchimento e esvaziamento do pulmão

Essas grandezas são cinemáticas (puro balanço de fluxo) e valem independentemente do modelo dinâmico. Sejam $I_{\text{up}}$ e $I_{\text{down}}$ as vazões a montante e a jusante.

Tempo de enchimento se a estação a jusante para:

Tempo de esvaziamento se a estação a montante para:

Taxa líquida de variação do nível do pulmão:

Em regime permanente subcrítico, as duas estações operam na taxa de demanda, $I_{\text{up}} = I_{\text{down}} = \lambda$, então $dQ/dt = 0$ e o pulmão se estabiliza numa fração de preenchimento igual à utilização:

Pulmão mínimo para cobrir uma parada planejada a montante de duração $t_{\text{out}}$ sem faltar material a jusante:

Tome o maior entre $Q_{\min}$ (cobertura de parada) e o requisito de amortecimento $(1 + \sqrt{2})\,\tau_L / R_{\text{eff}}$ da seção 4.2.

4.4 Fração operacional e resistência efetiva

A disponibilidade vive dentro de $R_{\text{eff}}$. Para uma estação com tempo de ciclo nominal $t_c$ em s/P, $m$ máquinas em paralelo e fração operacional

a resistência efetiva e a taxa máxima são

Valores típicos de $D$: um turno de 8 horas, $8/24 = 0{,}33$; dois turnos, $16/24 = 0{,}67$; três turnos, $1{,}0$; apenas fim de semana, $16/168 = 0{,}095$.

A fração operacional também escala a resposta dinâmica. Como $R_{\text{eff}} \propto 1/D$,

Menor disponibilidade retarda a resposta da linha às mudanças de demanda tanto quanto sua taxa. Isso substitui a afirmação linear anterior $\omega_{\text{eff}} = \omega_n\,D$; a escala em raiz quadrada é o que as equações de estado dão.

4.5 Processamento em lotes

Processos em lote (fornos, tanques de imersão, testes de vazamento) retêm as unidades por um tempo fixo $\tau$ em s e depois as liberam. Comportam-se como um atraso puro somado a uma média móvel.

Função de transferência normalizada, da taxa de entrada para a taxa de saída, com ganho DC unitário:

O módulo $\lvert\operatorname{sinc}(\omega\tau/2)\rvert$ mostra a etapa em lote deixando passar as tendências de demanda de baixa frequência e atenuando a variação de alta frequência, um filtro passa-baixas com o primeiro nulo em $\omega = 2\pi/\tau$.

Capacidade, com tamanho de lote $b$ em P e $N$ processadores em paralelo, mantida separada da resposta em frequência:

Escalonamento para fluxo contínuo. Para que uma etapa em lote se aproxime de fluxo contínuo no período takt $T_{\text{takt}}$, opere

lotes sobrepostos, cada um deslocado de $\tau / n$. O pulmão que alimenta a etapa deve conter pelo menos $n\,b$ unidades para que um lote completo esteja sempre pronto.

Exemplo resolvido. Forno com $\tau = 84{,}807$ s, $b = 4$ P, $N = 2$. Num período takt de $T_{\text{takt}} \approx 17{,}300$ s, $n = \operatorname{round}(4.9) = 5$ lotes sobrepostos, pulmão de alimentação de pelo menos $20$ unidades, e $\mu_{\text{batch}} = 2 \cdot 4 / 84{,}807 \approx 9.4 \times 10^{-5}$ P/s, cerca de 8 unidades por dia.

Integração.

• Fração operacional: a disponibilidade limitada de uma estação em lote alonga o intervalo de escalonamento efetivo e força ajustes de programação para proteger o fluxo a jusante.
• Domínio da frequência: uma etapa em lote é um filtro passa-baixas, então suaviza a variação de demanda de alta frequência antes que ela chegue às estações a jusante.
• MPF: alinhar os inícios de lote ao takt e manter $n$ próximo de $\tau / T_{\text{takt}}$ minimiza o material em processo que circula pela etapa, reduzindo sua contribuição à oscilação de pulmão $A_Q$ da seção 4.1.

4.6 Tradução financeira

O custo recuperável da ineficiência dinâmica é o custo do material em processo que circula sem vazão líquida, mais o custo de fazer a vazão subir e descer para produzir essa circulação. Ambos vêm das grandezas da seção 4.1, e ambos são dimensionados a partir de parâmetros medidos. A versão anterior desta seção repartia o custo por uma identidade de fator de potência, $K_{\text{VA}} = K\cos\phi$ com um resto reativo; essa identidade vinha de um triângulo de impedância RLC série que o modelo não produz, e suas partes não somavam o todo. Ela é substituída a seguir.

Tome uma banda de demanda com frequência dominante $\omega$ e amplitude $\Delta\lambda$ em torno da taxa média. O pulmão oscila com amplitude $\Delta Q = \Delta\lambda\,A_Q(\omega)$. A parte inevitável é o deslocamento estático de preenchimento $\Delta\lambda\,Q_{\max}R_{\text{eff}}$; o restante,

é material em processo oscilante que uma linha mais reativa e mais bem amortecida não carregaria. Abaixo da frequência de corte, o pulmão amplifica a oscilação e esse termo é positivo; bem acima de $\omega_n$, $A_Q$ cai abaixo do piso porque o pulmão filtra a oscilação, o termo se anula, e o modelo prevê nenhum custo extra de manutenção de estoque vindo das componentes de demanda de alta frequência. Essa é uma previsão verificável.

Custo de manutenção de estoque. Uma oscilação em torno de uma média fixa não aumenta, por si só, o material em processo médio: a meia-onda acima da média cancela a meia-onda abaixo. O custo entra pelo estoque de segurança. Para impedir que o vale da oscilação deixe a estação a jusante sem material, o nível médio do pulmão precisa ser elevado em cerca de $\Delta Q_{\text{exc}}$, e essa média elevada é carregada o tempo todo. Com $h$ o custo de posse por peça por unidade de tempo (capital, área de piso, movimentação, risco de obsolescência),

Se a restrição determinante for o tamanho físico do pulmão em vez do custo de posse, o mesmo $\Delta Q_{\text{exc}}$ gera um custo de capital único para o pulmão maior.

Custo de rampa. Produzir a oscilação força a própria vazão a subir e descer em $\Delta I = \Delta\lambda\,\lvert I/\lambda(j\omega)\rvert$. Com $c_{\text{ramp}}$ o custo incremental medido por unidade de oscilação de vazão (setups extras, mão de obra flexível, aceleração de pedidos) e $N_{\text{cyc}}$ ciclos de demanda no período,

Retorno de um projeto de redução de atraso. Reduzir $\tau_L$ aumenta $\omega_n$ e $\zeta$ e puxa $A_Q(\omega)$ na direção de seu piso, encolhendo tanto $\Delta Q_{\text{exc}}$ quanto $\Delta I$. Com os índices 0 e 1 para antes e depois,

Exemplo ilustrativo. Os números são redondos, não medidos. Uma linha com $R_{\text{eff}} = 100\ \text{s/P}$, $Q_{\max} = 1000\ \text{P}$ e $\tau_L = 5\times10^{5}\ \text{s}$ tem $\omega_n \approx 4.5\times10^{-6}\ \text{s}^{-1}$ e $\zeta \approx 0.22$. Uma oscilação mensal de demanda, $\omega \approx 2.4\times10^{-6}\ \text{s}^{-1}$, amplitude $\Delta\lambda = 0.002\ \text{P/s}$, dá $A_Q \approx 2.1\times10^{5}\ \text{s}$ contra um piso $Q_{\max}R_{\text{eff}} = 1.0\times10^{5}\ \text{s}$, então $\Delta Q \approx 420\ \text{P}$ com cerca de $220\ \text{P}$ evitáveis. Um programa de SMED e redução de lotes de transferência que corta $\tau_L$ pela metade, para $2.5\times10^{5}\ \text{s}$, eleva $\omega_n$ para $\approx 6.3\times10^{-6}\ \text{s}^{-1}$ e $\zeta$ para $\approx 0.32$, baixando $A_Q$ para $\approx 1.3\times10^{5}\ \text{s}$ e o material em processo evitável para cerca de $64\ \text{P}$. A economia anual de posse é $(220 - 64)\,\text{P}\times h \times 8760\ \text{h/ano}$ vinda do estoque de segurança reduzido; a economia de rampa do menor $\Delta I$ soma-se a ela. Dividir a soma pelo custo do projeto dá o ROI de primeiro ano, calculado a partir de uma mudança de $\tau_L$ em vez de um ganho de MPF suposto.

O exemplo também mostra o lado diagnóstico da seção 4.2. Os dois estados têm $\zeta$ bem abaixo do limiar de $0{,}777$ ($0{,}22$ e depois $0{,}32$), e a regra de pulmão $Q_{\max} \ge (1 + \sqrt{2})\,\tau_L / R_{\text{eff}}$ pede $\approx 12{,}000\ \text{P}$ antes do projeto e $\approx 6{,}000\ \text{P}$ depois, contra um pulmão real de $1{,}000\ \text{P}$. O modelo está sinalizando a linha como fortemente subamortecida: cortar $\tau_L$ pela metade é uma melhoria real, mas longe de suficiente, e ou o pulmão ou o atraso precisa se mover muito mais antes de o efeito chicote desaparecer.

Encaixe com ferramentas padrão. $K_{\text{WIP}}$ e $K_{\text{ramp}}$ são linhas de custo comuns, então suas reduções projetadas entram sem alteração no valor presente líquido, no payback e no custeio baseado em atividades. O que o modelo acrescenta é uma forma de dimensioná-las a partir de $\tau_L$, $R_{\text{eff}}$, $Q_{\max}$ e do espectro da demanda.