Un cadre inspiré des circuits pour les systèmes de fabrication

Résumé

Cet article présente un cadre de modélisation et d'analyse des systèmes de fabrication par analogie avec les circuits électriques. L'approche, appelée modèle de circuit de fabrication, offre une perspective fondée sur la physique qui complète les méthodes traditionnelles. En mettant les éléments de fabrication en correspondance avec des composants de circuit, le cadre permet d'appliquer la théorie des circuits pour quantifier la dynamique du système, dimensionner les tampons et traduire un changement opérationnel en coût. Les concepts clés incluent le facteur de puissance de fabrication (MPF), le dimensionnement des tampons par l'amortissement, et l'analyse fréquentielle du flux de production. Chaque paramètre du modèle se rattache à une mesure ou à une valeur de politique choisie, si bien que le cadre ne comporte aucune constante d'ajustement libre.

1. Introduction

Les systèmes de fabrication sont des réseaux complexes de procédés interconnectés, de flux de matières et d'échanges d'informations. Analyser et optimiser ces systèmes est essentiel pour améliorer l'efficacité, réduire le gaspillage et renforcer la réactivité. Les approches traditionnelles telles que la théorie des files d'attente, la simulation et la modélisation à événements discrets ont apporté des enseignements précieux, mais manquent souvent d'un cadre cohérent, fondé sur la physique, qui capture l'interaction dynamique des composants du système.

Cet article établit une analogie directe entre les systèmes de fabrication et les circuits électriques. Le modèle de circuit de fabrication qui en résulte offre une manière compacte et mesurable de raisonner sur la façon dont une ligne de production répond aux variations de la demande.

2. État de l'art

Les méthodes existantes d'analyse des systèmes de fabrication comprennent :

• La théorie des files d'attente : modélise les files et la congestion mais se concentre sur les moyennes en régime permanent plutôt que sur la réponse transitoire.
• La simulation à événements discrets : puissante pour les interactions complexes, mais dépourvue de cadre général sous forme close et pouvant être coûteuse en calcul.
• Le lean manufacturing : met l'accent sur la réduction du gaspillage et sur le flux, mais manque d'une approche quantitative de la dynamique au niveau du système.
• La dynamique des systèmes : utilise des boucles de rétroaction et des retards pour modéliser le comportement ; le modèle de circuit en est une spécialisation avec des éléments nommés et mesurables.
• La factory physics : fournit des relations quantitatives comme la loi de Little, mais s'arrête avant une analogie de circuit au niveau des composants.

Le modèle de circuit de fabrication s'appuie sur ces fondations avec un cadre qui est intuitif, mathématiquement explicite et ancré dans la mesure.

3. Le modèle de circuit de fabrication

Unités de base tout au long : pièces $\text{P}$ et secondes $\text{s}$.

Tableau 1 : analogies circuit-fabrication

La demande entre comme une source de courant, pas une source de tension : le client fixe un débit, et soit la ligne le tient, soit le tampon et le carnet de commandes absorbent l'écart. Comme la demande est une source de courant, les grandeurs fréquentielles ci-dessous proviennent des fonctions de transfert des équations d'état de la section 3.3, et non d'un triangle d'impédance RLC série : dans ce modèle, $R_{\text{eff}}$ figure dans le terme de raideur, pas dans le terme dissipatif.

3.1 Les trois éléments, ancrés dans la mesure

• Résistance $R_{\text{eff}} = t_c / (m\,D)$ en s/P. Le temps de cycle effectif d'un poste : temps de cycle nominal $t_c$ divisé par le nombre de machines en parallèle $m$ et le taux de marche $D$ (fraction du temps calendaire pendant laquelle le poste peut fonctionner). Son inverse $\mu = 1 / R_{\text{eff}}$ en P/s est le débit maximal du poste.
• Capacité $C = Q_{\max}$ en P. Le tampon physique entre deux postes. Le potentiel du tampon est la fraction de remplissage sans dimension $V_C = Q / Q_{\max}$, donc $C = Q / V_C$ a pour unité la pièce. La relation du condensateur est exacte : $\dot Q = I_{\text{in}} - I_{\text{out}}$.
• Inductance $L = R_{\text{eff}}\,\tau_L$ en s²/P. Un poste ne peut pas changer son débit de sortie instantanément : il faut repositionner les en-cours, changer l'outillage, ajuster les effectifs et l'équilibrage de la ligne. Le retard $\tau_L$ en s est le temps mesuré pour que le débit atteigne 63 % d'un nouveau niveau après un échelon de demande soutenu. La quantité stockée $\tfrac{1}{2} L I^2 = \tfrac{1}{2} R_{\text{eff}}\,\tau_L\,I^2$ a pour unité la pièce et représente le travail engagé en cours qui s'achèvera même si les lancements s'arrêtent.

3.2 Lois de circuit, avec saturation

Loi de débit (une résistance qui sature). Tant que le tampon a de la marge, le débit est fixé par la fraction de remplissage et le temps de cycle effectif :

Unités : sans dimension divisé par s/P donne P/s.

Ce que représente la loi de débit. C'est une politique de lancement, pas la physique de la machine. Un poste laissé à lui-même fonctionne à son plein débit $\mu = 1/R_{\text{eff}}$ dès qu'une pièce attend, ce qui rendrait le débit brut égal à $\mu$ pour tout $Q > 0$, et non proportionnel à $Q$. Ici, $I = V_C / R_{\text{eff}}$ se lit comme une règle d'appel proportionnelle aux en-cours : le poste est autorisé à appeler et à lancer du travail à un débit proportionnel au taux de remplissage de son tampon amont, la forme linéaire du principe CONWIP et kanban selon lequel une étape ne doit pas devancer son signal aval. C'est ce qui rend légitime la forme linéaire en $Q$ : c'est une loi de commande choisie, pas une affirmation sur l'équipement. Le gain d'appel de cette politique est fixé à $1/R_{\text{eff}}$ pour que le débit autorisé atteigne exactement la vraie capacité lorsque le tampon est plein, et le reste de l'article utilise ce réglage, si bien que $\mu = 1/R_{\text{eff}}$ partout et que la clause de saturation ci-dessous se déclenche toujours à la vraie capacité. Faire fonctionner un gain volontairement bridé $g < 1/R_{\text{eff}}$ est possible, mais c'est un choix de conception distinct qui abaisse la capacité effective à $g$ ; il n'est pas considéré ici. Ainsi $\omega_n$, $\zeta$, $\phi$ et $A_Q$ ci-dessous décrivent une ligne commandée, pas une ligne en boucle ouverte.

Quand la demande dépasse la capacité, $\lambda \ge \mu = 1 / R_{\text{eff}}$, la fraction de remplissage se cale à 1, le débit se plafonne à $\mu$, et la demande non satisfaite s'accumule en carnet de commandes $B$ avec $\dot B = \lambda - \mu$. Une approximation lissée à la transition est la forme en conductances en parallèle

avec les trois termes en s/P. Cela remplace la loi linéaire $V = I R$, qui ne ferme pas dimensionnellement et ne capture pas le fait que le débit sature à la capacité.

Conservation du flux (loi des nœuds de Kirchhoff). À chaque nœud, les unités se conservent :

3.3 Dynamique du second ordre

On couple le tampon (condensateur) au retard de débit (bobine). Avec $Q$ le niveau du tampon et $I$ la sortie du poste :

En éliminant $I$, on obtient un système amorti du second ordre en niveau de tampon :

avec fréquence propre et taux d'amortissement

Les deux expressions vérifient les unités. Chaque paramètre est fixé avant l'exécution du modèle : $R_{\text{eff}}$ par une étude de temps de cycle (et servant de gain d'appel, 3.2), $C$ par le tampon, $\tau_L$ par un essai en échelon de sortie. Il n'y a aucune constante libre à ajuster.

4. Concepts clés et applications

4.1 Facteur de puissance de fabrication (MPF)

En transformant les équations d'état de la section 3.3, avec la demande $\lambda$ comme source de courant et la sortie du poste $I$ comme réponse, on obtient la fonction de transfert demande-vers-débit

un passe-bas du second ordre avec les $\omega_n$ et $\zeta$ de la section 3.3. À une fréquence de demande $\omega$, le débit est en retard sur la demande du déphasage

Le numérateur et le dénominateur ont tous deux pour unité $\text{s}^{-2}$, donc $\tan\phi$ est sans dimension. On définit

Interprétation. On décompose la réponse en débit à la fréquence $\omega$ en une part en phase avec la demande et une part en quadrature. La part en phase, proportionnelle à $\cos\phi$, est le flux qui sert la demande. La part en quadrature, proportionnelle à $\sin\phi$, remplit le tampon sur une moitié du cycle et le vide sur l'autre, avec une contribution nette nulle au débit, tout en consommant de la main-d'œuvre, du temps machine et un coût de portage. Le MPF est la fraction de la réponse en débit qui suit la demande, sur une échelle de 0 à 1 : il vaut 1 quand la sortie est en phase et tombe à 0 en $\omega = \omega_n$, où la réponse est un pur échange de tampon. Pour $\omega \ge \omega_n$, le $\cos\phi$ brut est négatif, ce qui signifie que la sortie évolue à l'inverse de la demande, et le MPF est ramené à 0 ; la grandeur pertinente dans ce régime est le déphasage $\phi$ lui-même, pas son cosinus. Le MPF atteignant 0 en $\omega_n$ est la forme fréquentielle de la mise en garde de la section 4.2 : ne pas laisser la variation dominante de la demande se placer à la fréquence propre.

Le levier, c'est $\tau_L$. Pour $\omega < \omega_n$, réduire le retard d'ajustement du débit $\tau_L$ abaisse $\tan\phi$ et augmente le MPF, et augmente en même temps $\omega_n$ et $\zeta$ (4.2). La réduction des changements de série (SMED), la flexibilité des effectifs et des lots de transfert plus petits réduisent tous $\tau_L$. Le retard a un plancher : quand $\tau_L \to 0$, le modèle se réduit au retard du premier ordre du tampon $\phi \to \arctan(\omega\,Q_{\max}R_{\text{eff}})$, fixé par la constante de temps du tampon $Q_{\max}R_{\text{eff}}$. Un tampon plus grand augmente l'amortissement (4.2) mais augmente aussi ce retard résiduel, si bien que la taille du tampon arbitre entre le suivi de la demande et la stabilité au lieu d'améliorer les deux.

Oscillation des en-cours. La réponse du tampon à la demande découle de $\dot Q = \lambda - I$ :

avec $A_Q$ en secondes. Une composante de demande d'amplitude $\Delta\lambda$ à la fréquence $\omega$ entraîne une oscillation du tampon d'amplitude $\Delta Q = \Delta\lambda\,A_Q(\omega)$. À basse fréquence, $A_Q \to Q_{\max}R_{\text{eff}}$, si bien qu'un changement lent de demande de taille $\Delta\lambda$ décale la fraction de remplissage de $\Delta\lambda/\mu$, ce qui correspond à $V_C^{\ast} = \lambda/\mu$ de la section 4.3. Près de $\omega_n$, $A_Q$ dépasse ce plancher dès que $\zeta < \tfrac{1}{2}\sqrt{1 + \sqrt{2}} \approx 0{,}777$ : c'est l'effet coup de fouet (bullwhip). Le seuil n'est pas $1/\sqrt{2}$, car le zéro en $s = -1/\tau_L$ dans $Q/\lambda$ le relève ; voir 4.2. Respecter la règle de tampon de la section 4.2 maintient $A_Q$ au niveau de son plancher basse fréquence ou en dessous sur toute la bande de demande, et bien au-dessus de $\omega_n$ le tampon filtre l'oscillation, si bien que $A_Q$ passe sous le plancher.

4.2 Dimensionnement des tampons par l'amortissement

Une version antérieure de ce cadre prescrivait de fonctionner à la résonance, $\omega_n = \omega_{\text{takt}}$. C'est une erreur. Près de $\omega_n$, une ligne peu amortie amplifie les oscillations de tampon (le mode coup de fouet ou de pompage), un comportement à supprimer, pas à rechercher.

L'objectif de conception est une réponse demande-vers-en-cours $A_Q(\omega)$ (4.1) sans pic au-dessus de sa valeur basse fréquence. La fonction de transfert du tampon $Q/\lambda$ porte un zéro en $s = -1/\tau_L$ qui relève son côté haute fréquence, si bien que la condition de réponse plate n'est pas $\zeta \ge 1/\sqrt{2}$. En écrivant $|Q/\lambda|^2$ sous la forme $(x + 4\zeta^2)/\bigl((1-x)^2 + 4\zeta^2 x\bigr)$ avec $x = (\omega/\omega_n)^2$, sa pente en $x = 0$ est $1 + 8\zeta^2 - 16\zeta^4$, qui reste positive (la réponse s'élève au-dessus de sa valeur continue) jusqu'à

Unités : s divisé par s/P donne P. Le seuil simple $\zeta \ge 1/\sqrt{2}$ s'applique quand même, mais à la réponse demande-vers-débit $I/\lambda$, qui est un passe-bas pur du second ordre sans zéro : dimensionner selon celui-ci quand l'objectif est un débit de sortie lisse, et selon la règle en $(1 + \sqrt{2})$ pour dimensionner le tampon contre l'oscillation des en-cours. Pour un retour des en-cours monotone, sans aucun dépassement, exiger $\zeta \ge 1$ :

$\omega_n$ reste utile : c'est la fréquence de variation de la demande que la ligne rejette le moins bien. Maintenir $\zeta$ au-dessus de $\tfrac{1}{2}\sqrt{1 + \sqrt{2}}$ garantit que les fluctuations de tampon dans cette bande ne sont pas amplifiées.

4.3 Temps de remplissage et de vidange des tampons

Ces grandeurs sont cinématiques (pur bilan de flux) et valent quel que soit le modèle dynamique. Soient $I_{\text{up}}$ et $I_{\text{down}}$ les débits amont et aval.

Temps de remplissage si le poste aval s'arrête :

Temps de vidange si le poste amont s'arrête :

Taux net de variation du niveau du tampon :

En régime permanent sous-critique, les deux postes fonctionnent au débit de demande, $I_{\text{up}} = I_{\text{down}} = \lambda$, donc $dQ/dt = 0$ et le tampon se stabilise à une fraction de remplissage égale au taux d'utilisation :

Tampon minimal pour couvrir un arrêt amont planifié de durée $t_{\text{out}}$ sans affamer l'aval :

Prendre le plus grand de $Q_{\min}$ (couverture d'arrêt) et de l'exigence d'amortissement $(1 + \sqrt{2})\,\tau_L / R_{\text{eff}}$ de la section 4.2.

4.4 Taux de marche et résistance effective

La disponibilité vit à l'intérieur de $R_{\text{eff}}$. Pour un poste de temps de cycle nominal $t_c$ en s/P, de $m$ machines en parallèle et de taux de marche

la résistance effective et le débit maximal sont

Valeurs typiques de $D$ : une équipe de 8 heures, $8/24 = 0{,}33$ ; deux équipes, $16/24 = 0{,}67$ ; trois équipes, $1{,}0$ ; week-end seulement, $16/168 = 0{,}095$.

Le taux de marche met aussi à l'échelle la réponse dynamique. Comme $R_{\text{eff}} \propto 1/D$,

Une disponibilité plus faible ralentit la réponse de la ligne aux variations de demande autant que son débit. Cela remplace l'affirmation linéaire antérieure $\omega_{\text{eff}} = \omega_n\,D$ ; la mise à l'échelle en racine carrée est ce que donnent les équations d'état.

4.5 Traitement par lots

Les procédés par lots (fours, bacs de trempe, tests d'étanchéité) retiennent les unités pendant un temps fixe $\tau$ en s, puis les libèrent. Ils se comportent comme un retard pur plus une moyenne glissante.

Fonction de transfert normalisée, du débit d'entrée au débit de sortie, de gain continu unitaire :

Le module $\lvert\operatorname{sinc}(\omega\tau/2)\rvert$ montre l'étape par lots laisser passer les tendances de demande basse fréquence et atténuer la variation haute fréquence, un filtre passe-bas dont le premier zéro est en $\omega = 2\pi/\tau$.

Capacité, avec une taille de lot $b$ en P et $N$ processeurs en parallèle, tenue à l'écart de la réponse fréquentielle :

Décalage pour un flux continu. Pour qu'une étape par lots approche un flux continu à la période takt $T_{\text{takt}}$, faire fonctionner

lots qui se chevauchent, chacun décalé de $\tau / n$. Le tampon qui alimente l'étape doit contenir au moins $n\,b$ unités pour qu'un lot complet soit toujours prêt.

Exemple résolu. Four avec $\tau = 84{,}807$ s, $b = 4$ P, $N = 2$. À une période takt de $T_{\text{takt}} \approx 17{,}300$ s, $n = \operatorname{round}(4.9) = 5$ lots qui se chevauchent, tampon d'alimentation d'au moins $20$ unités, et $\mu_{\text{batch}} = 2 \cdot 4 / 84{,}807 \approx 9.4 \times 10^{-5}$ P/s, soit environ 8 unités par jour.

Intégration.

• Taux de marche : la disponibilité limitée d'un poste par lots allonge l'intervalle de décalage effectif et impose des ajustements d'ordonnancement pour protéger le flux aval.
• Domaine fréquentiel : une étape par lots est un filtre passe-bas, elle lisse donc la variation de demande haute fréquence avant qu'elle n'atteigne les postes aval.
• MPF : aligner les départs de lots sur le takt et maintenir $n$ près de $\tau / T_{\text{takt}}$ minimise les en-cours qui circulent dans l'étape, ce qui réduit sa contribution à l'oscillation de tampon $A_Q$ de la section 4.1.

4.6 Traduction financière

Le coût récupérable de l'inefficacité dynamique est le coût des en-cours qui circulent sans débit net, plus le coût de faire monter et descendre le débit pour produire cette circulation. Les deux découlent des grandeurs de la section 4.1, et les deux se dimensionnent à partir de paramètres mesurés. La version antérieure de cette section répartissait le coût par une identité de facteur de puissance, $K_{\text{VA}} = K\cos\phi$ avec un reste réactif ; cette identité provenait d'un triangle d'impédance RLC série que le modèle ne produit pas, et ses parts ne totalisaient pas l'ensemble. Elle est remplacée ci-dessous.

Prenons une bande de demande de fréquence dominante $\omega$ et d'amplitude $\Delta\lambda$ autour du débit moyen. Le tampon oscille avec une amplitude $\Delta Q = \Delta\lambda\,A_Q(\omega)$. La part inévitable est le décalage statique de remplissage $\Delta\lambda\,Q_{\max}R_{\text{eff}}$ ; le reste,

correspond aux en-cours oscillants qu'une ligne mieux amortie et plus réactive ne porterait pas. Sous la fréquence de coupure, le tampon amplifie l'oscillation et ce terme est positif ; bien au-dessus de $\omega_n$, $A_Q$ passe sous le plancher parce que le tampon filtre l'oscillation, le terme se ramène à zéro, et le modèle prédit aucun surcoût de portage venant des composantes de demande haute fréquence. C'est une prédiction vérifiable.

Coût de portage. Une oscillation autour d'une moyenne fixe n'augmente pas à elle seule les en-cours moyens : la demi-période au-dessus de la moyenne annule celle en dessous. Le coût entre par le stock de sécurité. Pour empêcher le creux de l'oscillation d'affamer le poste aval, il faut relever le niveau moyen du tampon d'environ $\Delta Q_{\text{exc}}$, et cette moyenne relevée est portée tout le temps. Avec $h$ le coût de détention par pièce et par unité de temps (capital, surface au sol, manutention, risque d'obsolescence),

Si la contrainte déterminante est la taille physique du tampon plutôt que le coût de portage, le même $\Delta Q_{\text{exc}}$ entraîne plutôt un coût d'investissement ponctuel pour le tampon plus grand.

Coût de rampe. Produire l'oscillation force le débit lui-même à monter et descendre de $\Delta I = \Delta\lambda\,\lvert I/\lambda(j\omega)\rvert$. Avec $c_{\text{ramp}}$ le coût marginal mesuré par unité d'oscillation de débit (changements de série supplémentaires, effectifs flexibles, accélérations) et $N_{\text{cyc}}$ cycles de demande dans la période,

Retour sur un projet de réduction du retard. Réduire $\tau_L$ augmente $\omega_n$ et $\zeta$ et rapproche $A_Q(\omega)$ de son plancher, réduisant à la fois $\Delta Q_{\text{exc}}$ et $\Delta I$. Avec les indices 0 et 1 pour l'avant et l'après,

Exemple illustratif. Les nombres sont ronds, non mesurés. Une ligne avec $R_{\text{eff}} = 100\ \text{s/P}$, $Q_{\max} = 1000\ \text{P}$ et $\tau_L = 5\times10^{5}\ \text{s}$ a $\omega_n \approx 4.5\times10^{-6}\ \text{s}^{-1}$ et $\zeta \approx 0.22$. Une oscillation de demande mensuelle, $\omega \approx 2.4\times10^{-6}\ \text{s}^{-1}$, d'amplitude $\Delta\lambda = 0.002\ \text{P/s}$, donne $A_Q \approx 2.1\times10^{5}\ \text{s}$ contre un plancher $Q_{\max}R_{\text{eff}} = 1.0\times10^{5}\ \text{s}$, donc $\Delta Q \approx 420\ \text{P}$ dont environ $220\ \text{P}$ évitables. Un programme SMED et de réduction des lots de transfert qui divise $\tau_L$ par deux à $2.5\times10^{5}\ \text{s}$ porte $\omega_n$ à $\approx 6.3\times10^{-6}\ \text{s}^{-1}$ et $\zeta$ à $\approx 0.32$, faisant tomber $A_Q$ à $\approx 1.3\times10^{5}\ \text{s}$ et les en-cours évitables à environ $64\ \text{P}$. L'économie annuelle de portage est $(220 - 64)\,\text{P}\times h \times 8760\ \text{h/an}$ grâce au stock de sécurité réduit ; l'économie de rampe due au plus petit $\Delta I$ s'y ajoute. Diviser la somme par le coût du projet donne le ROI de première année, calculé à partir d'un changement de $\tau_L$ plutôt que d'un gain de MPF supposé.

L'exemple montre aussi le versant diagnostic de la section 4.2. Les deux états ont un $\zeta$ bien en dessous du seuil de $0{,}777$ ($0{,}22$, puis $0{,}32$), et la règle de tampon $Q_{\max} \ge (1 + \sqrt{2})\,\tau_L / R_{\text{eff}}$ réclame $\approx 12{,}000\ \text{P}$ avant le projet et $\approx 6{,}000\ \text{P}$ après, contre un tampon réel de $1{,}000\ \text{P}$. Le modèle signale la ligne comme fortement sous-amortie : diviser $\tau_L$ par deux est une vraie amélioration mais reste très loin du compte, et il faut que le tampon ou le retard bouge bien davantage avant que l'effet coup de fouet ne disparaisse.

Compatibilité avec les outils standard. $K_{\text{WIP}}$ et $K_{\text{ramp}}$ sont des postes de coût ordinaires, si bien que leurs réductions projetées s'intègrent telles quelles dans la valeur actuelle nette, le délai de récupération et la comptabilité par activités. Ce que le modèle apporte, c'est une manière de les dimensionner à partir de $\tau_L$, $R_{\text{eff}}$, $Q_{\max}$ et du spectre de la demande.