A Circuit-Inspired Framework for Manufacturing Systems

Abstract

This paper presents a framework for modeling and analyzing manufacturing systems by drawing analogies with electrical circuits. The approach, termed the Manufacturing Circuit Model, provides a physics-based perspective that complements traditional methods. By mapping manufacturing elements to circuit components, the framework enables the application of circuit theory to quantify system dynamics, size buffers, and translate operational change into cost. Key concepts include the Manufacturing Power Factor (MPF), damping-based buffer sizing, and frequency-domain analysis of production flow. Every parameter in the model traces to a measurement, so the framework carries no free fitting constants.

1. Introduction

Manufacturing systems are complex networks of interconnected processes, material flows, and information exchanges. Analyzing and optimizing these systems is critical for improving efficiency, reducing waste, and enhancing responsiveness. Traditional approaches such as queuing theory, simulation, and discrete-event modeling have provided valuable insights but often lack a cohesive physics-based framework that captures the dynamic interplay of system components.

This paper draws a direct analogy between manufacturing systems and electrical circuits. The resulting Manufacturing Circuit Model gives a compact, measurable way to reason about how a production line responds to changes in demand.

2. Literature Review

Existing manufacturing system analysis methods include:

Queuing theory: models waiting lines and congestion but focuses on steady-state averages rather than transient response.
Discrete-event simulation: powerful for complex interactions but lacks a general closed-form framework and can be computationally intensive.
Lean manufacturing: emphasizes waste reduction and flow but lacks a quantitative approach to system-level dynamics.
System dynamics: uses feedback loops and delays to model behavior; the circuit model is a specialization of this idea with named, measurable elements.
Factory physics: provides quantitative relationships such as Little's Law but stops short of a component-level circuit analogy.

The Manufacturing Circuit Model builds on these foundations with a framework that is intuitive, mathematically explicit, and grounded in measurement.

3. The Manufacturing Circuit Model

Base units throughout: parts $\text{P}$ and seconds $\text{s}$ .

Table 1: Circuit-manufacturing analogies

Circuit element	Manufacturing equivalent	Manufacturing units	Interpretation
Charge $q$	Work in process $Q$	P	units held in a buffer
Current $I$	Throughput	P/s	flow of units through a station
Current source	Demand rate $\lambda$	P/s	exogenous pull from the customer
Resistance $R_{\text{eff}}$	Effective cycle time	s/P	$t_c / (m\,D)$ : cycle time over parallelism and availability
Capacitance $C$	Buffer capacity $Q_{\max}$	P	units the buffer can hold
Node potential $V_C$	Buffer fill fraction	dimensionless	$Q / Q_{\max}$ , ranges 0 to 1
Time constant $\tau_L$	Throughput adjustment lag	s	63% rise time of output after a demand step
Inductance $L$	Production inertia	s²/P	$R_{\text{eff}}\,\tau_L$
Reactance $X$	Out-of-phase response	s/P	$\omega L - 1/(\omega C)$
Impedance $Z$	Total dynamic opposition	s/P	$\sqrt{R_{\text{eff}}^2 + X^2}$

Demand enters as a current source, not a voltage source: the customer sets a rate, and the line either meets it or the buffer and backlog absorb the difference.

3.1 The three elements, grounded

Resistance $R_{\text{eff}} = t_c / (m\,D)$ in s/P. The effective cycle time of a station: nominal cycle time $t_c$ divided by the number of parallel machines $m$ and the duty cycle $D$ (fraction of calendar time the station can run). Its reciprocal $\mu = 1 / R_{\text{eff}}$ in P/s is the station's maximum rate.
Capacitance $C = Q_{\max}$ in P. The physical buffer between two stations. The buffer's potential is the dimensionless fill fraction $V_C = Q / Q_{\max}$ , so $C = Q / V_C$ has units of parts. The capacitor relation is exact: $\dot Q = I_{\text{in}} - I_{\text{out}}$ .
Inductance $L = R_{\text{eff}}\,\tau_L$ in s²/P. A station cannot change its output rate instantly: WIP has to be repositioned, tooling changed, staffing and line balance adjusted. The lag $\tau_L$ in s is the measured time for output to reach 63% of a new level after a sustained demand step. The stored quantity $\tfrac{1}{2} L I^2 = \tfrac{1}{2} R_{\text{eff}}\,\tau_L\,I^2$ has units of parts and represents committed in-process work that will complete even if releases stop.

3.2 Circuit laws, with saturation

Throughput law (a saturating resistor). While the buffer has headroom, throughput is set by the fill fraction and the effective cycle time:

I = \frac{V_C}{R_{\text{eff}}} = \frac{Q}{Q_{\max}\,R_{\text{eff}}}, \qquad 0 \le V_C \le 1

Units: dimensionless divided by s/P gives P/s. When demand exceeds capacity, $\lambda \ge \mu = 1 / R_{\text{eff}}$ , the fill fraction pins at 1, throughput clamps at $\mu$ , and unmet demand accumulates as backlog $B$ with $\dot B = \lambda - \mu$ . A smooth approximation across the transition is the parallel-conductance form

\frac{1}{I} = \frac{1}{\lambda} + R_{\text{eff}} \quad\Longrightarrow\quad I \to \min(\lambda,\ \mu)

with all three terms in s/P. This replaces the linear $V = I R$ , which does not close dimensionally and does not capture that throughput saturates at capacity.

Conservation of flow (Kirchhoff's current law). At every node, units are conserved:

\sum I_{\text{in}} - \sum I_{\text{out}} = \frac{dQ}{dt}

3.3 Second-order dynamics

Couple the buffer (capacitor) to the throughput lag (inductor). With $Q$ the buffer level and $I$ the station output:

\dot Q = \lambda(t) - I, \qquad \tau_L\,\dot I = \frac{Q}{Q_{\max}\,R_{\text{eff}}} - I

Eliminating $I$ gives a damped second-order system in the buffer level:

\tau_L\,\ddot Q + \dot Q + \frac{Q}{Q_{\max}\,R_{\text{eff}}} = \tau_L\,\dot\lambda + \lambda

with natural frequency and damping ratio

\omega_n = \frac{1}{\sqrt{L\,C}} = \frac{1}{\sqrt{R_{\text{eff}}\,\tau_L\,Q_{\max}}} \;\; [\text{s}^{-1}], \qquad \zeta = \frac{R_{\text{eff}}}{2}\sqrt{\frac{C}{L}} = \frac{1}{2}\sqrt{\frac{Q_{\max}\,R_{\text{eff}}}{\tau_L}} \;\; [\text{--}]

Both expressions unit-check. Every parameter is measured: $R_{\text{eff}}$ from a cycle-time study, $C$ from the buffer, $\tau_L$ from an output step test. There are no free constants.

4. Key Concepts and Applications

4.1 Manufacturing Power Factor (MPF)

At the dominant frequency $\omega$ of demand variation, the station's throughput lags demand by a phase angle $\phi(\omega)$ :

\tan\phi(\omega) = \frac{X}{R_{\text{eff}}} = \omega\,\tau_L - \frac{1}{\omega\,Q_{\max}\,R_{\text{eff}}}

\text{MPF}(\omega) = \cos\phi(\omega) = \frac{R_{\text{eff}}}{\sqrt{R_{\text{eff}}^2 + \left(\omega L - \dfrac{1}{\omega C}\right)^2}}

Both terms in $\tan\phi$ are dimensionless. MPF runs from 0 to 1.

Interpretation. $\sin\phi$ is the share of production effort that builds work in process above the demand-tracking level during one part of the cycle and draws it back down during the rest. Its net contribution to throughput over a full cycle is zero. Unlike reactive power in an electrical circuit, this effort is not returned: the labor, machine time, energy, and carrying cost of that circulating WIP are all spent. MPF is therefore the fraction of production activity that yields net output at the demand frequency, and $1 - \text{MPF}$ is genuine loss, not stored energy.

The lever is $\tau_L$ . MPF improves as $\omega\,\tau_L$ shrinks, that is, as the throughput adjustment lag falls. Changeover reduction (SMED), flexible staffing, and smaller transfer batches all cut $\tau_L$ , which at the same time lowers $L$ , raises the damping ratio $\zeta$ , and flattens $\phi(\omega)$ across the demand band. Chasing $\omega = \omega_n$ is not the goal: see 4.2.

4.2 Buffer sizing by damping

An earlier version of this framework prescribed operating at resonance, $\omega_n = \omega_{\text{takt}}$ . That is wrong. Near $\omega_n$ a lightly damped line amplifies buffer swings (the bullwhip or hunting mode), which is the behavior to suppress, not to seek. The design target is a demand-to-WIP response with no resonant peak, which requires $\zeta \ge 1/\sqrt{2}$ :

\frac{1}{2}\sqrt{\frac{Q_{\max}\,R_{\text{eff}}}{\tau_L}} \ge \frac{1}{\sqrt{2}} \quad\Longrightarrow\quad Q_{\max} \ge \frac{2\,\tau_L}{R_{\text{eff}}}

Units: s divided by s/P gives P. For a fully non-oscillatory (critically damped) response, require $\zeta \ge 1$ :

Q_{\max} \ge \frac{4\,\tau_L}{R_{\text{eff}}}

$\omega_n$ is still useful: it is the frequency of demand variation the line is least able to reject. Keeping $\zeta$ above $1/\sqrt{2}$ ensures fluctuations in that band are not amplified in the buffer.

4.3 Buffer fill and drain timing

These quantities are kinematic (pure flow balance) and hold regardless of the dynamic model. Let $I_{\text{up}}$ and $I_{\text{down}}$ be the upstream and downstream throughputs.

Time to fill if the downstream station stops:

t_{\text{fill}} = \frac{Q_{\max} - Q}{I_{\text{up}}}

Time to drain if the upstream station stops:

t_{\text{drain}} = \frac{Q}{I_{\text{down}}}

Net rate of change of the buffer level:

\frac{dQ}{dt} = I_{\text{up}} - I_{\text{down}}

In subcritical steady state both stations run at the demand rate, $I_{\text{up}} = I_{\text{down}} = \lambda$ , so $dQ/dt = 0$ and the buffer settles at a fill fraction equal to utilization:

V_C^{\ast} = \lambda\,R_{\text{eff}} = \frac{\lambda}{\mu} \quad (\text{valid while } \lambda \le \mu)

Minimum buffer to cover a planned upstream outage of duration $t_{\text{out}}$ without starving downstream:

Q_{\min} = I_{\text{down}}\,t_{\text{out}} = \frac{t_{\text{out}}}{R_{\text{eff,down}}}

Take the larger of $Q_{\min}$ (outage coverage) and the damping requirement $2\,\tau_L / R_{\text{eff}}$ from 4.2.

4.4 Duty cycle and effective resistance

Availability lives inside $R_{\text{eff}}$ . For a station with nominal cycle time $t_c$ in s/P, $m$ parallel machines, and duty cycle

D = \frac{\text{available production time}}{\text{total calendar time}}

the effective resistance and maximum rate are

R_{\text{eff}} = \frac{t_c}{m\,D} \;\; [\text{s/P}], \qquad \mu = \frac{1}{R_{\text{eff}}} = \frac{m\,D}{t_c} \;\; [\text{P/s}], \qquad I_{\max,\text{day}} = 86{,}400\,\mu

Typical values of $D$ : one 8-hour shift, $8/24 = 0.33$ ; two shifts, $16/24 = 0.67$ ; three shifts, $1.0$ ; weekend-only, $16/168 = 0.095$ .

The duty cycle also scales the dynamic response. Since $R_{\text{eff}} \propto 1/D$ ,

\omega_n = \sqrt{\frac{m\,D}{t_c\,\tau_L\,Q_{\max}}} \;\propto\; \sqrt{D}

Lower availability slows the line's response to demand shifts as well as its rate. This replaces the earlier linear claim $\omega_{\text{eff}} = \omega_n\,D$ ; the square-root scaling is what the state equations give.

4.5 Batch processing

Batch processes (ovens, soak tanks, leak tests) hold units for a fixed time $\tau$ in s, then release them. They behave as a pure delay plus a moving average.

Normalized transfer function, from input rate to output rate, with unit DC gain:

H(\omega) = e^{-j\omega\tau/2}\,\operatorname{sinc}\!\left(\frac{\omega\tau}{2}\right) \quad (\text{dimensionless})

The magnitude $\lvert\operatorname{sinc}(\omega\tau/2)\rvert$ shows the batch stage passing low-frequency demand trends and attenuating high-frequency variation, a low-pass filter with its first null at $\omega = 2\pi/\tau$ .

Capacity, with batch size $b$ in P and $N$ parallel processors, kept separate from the frequency response:

\mu_{\text{batch}} = \frac{N\,b}{\tau} \;\; [\text{P/s}]

Staggering for continuous flow. To make a batch stage approximate continuous flow at takt period $T_{\text{takt}}$ , run

n = \operatorname{round}\!\left(\frac{\tau}{T_{\text{takt}}}\right)

overlapping batches, each offset by $\tau / n$ . The buffer feeding the stage must hold at least $n\,b$ units so a full batch is always ready.

Worked example. Oven with $\tau = 84{,}807$ s, $b = 4$ P, $N = 2$ . At a takt period of $T_{\text{takt}} \approx 17{,}300$ s, $n = \operatorname{round}(4.9) = 5$ overlapping batches, feed buffer at least $20$ units, and $\mu_{\text{batch}} = 2 \cdot 4 / 84{,}807 \approx 9.4 \times 10^{-5}$ P/s, about 8 units per day.

Integration.

Duty cycle: limited availability of a batch station lengthens the effective stagger interval and forces schedule adjustments to protect downstream flow.
Frequency domain: a batch stage is a low-pass filter, so it smooths high-frequency demand variation before it reaches downstream stations.
MPF: aligning batch starts with takt and holding $n$ near $\tau / T_{\text{takt}}$ minimizes the WIP that circulates through the stage, which raises system MPF.

4.6 Financial translation

The phase split from 4.1 carries directly into cost. Let $K$ be total production cost over a period (labor, machine time, energy, overhead), and $\phi$ the demand-band phase angle with $\cos\phi = \text{MPF}$ .

Value-adding cost:

K_{\text{VA}} = K\cdot\cos\phi = K\cdot\text{MPF}

Circulating (non-value-adding) cost, $K_{\text{NVA}} = K\,(1 - \text{MPF})$ , splits by the two reactive terms, with $Z = \sqrt{R_{\text{eff}}^2 + X^2}$ and $X = \omega L - 1/(\omega C)$ :

K_L = K\cdot\frac{\omega L}{Z}\,\sin\phi \qquad (\text{inertia: setups, changeovers, ramp effort})

K_C = K\cdot\frac{1}{\omega C\,Z}\,\sin\phi \qquad (\text{inventory: carrying and handling of oscillating WIP})

Each fraction $\omega L / Z$ and $1/(\omega C Z)$ is dimensionless. These costs are spent, not recoverable; the electrical picture of energy returning to the source does not apply here.

Return on a lag-reduction project. A project that lowers $\tau_L$ raises MPF by $\Delta\text{MPF}$ :

\text{ROI} = \frac{K\cdot\Delta\text{MPF}}{\text{project cost}}

Example: a SMED effort costing $250,000 lifts MPF from 0.70 to 0.85 on a line with $4,000,000 annual cost:

\text{ROI} = \frac{\$4{,}000{,}000 \times 0.15}{\$250{,}000} = \frac{\$600{,}000}{\$250{,}000} = 2.4 \quad (240\%\ \text{first-year return})

Industry forms. High fixed-cost lines (semiconductor, automotive): recoverable value is $(\text{full-capacity revenue} \times \text{MPF}) - \text{fixed cost}$ . High variable-cost lines (assembly, food processing): $\text{margin} \times \text{output} \times \text{MPF}$ .

Fit with standard tools. MPF is a dimensionless multiplier on the cost base, so its projected gain slots into net present value, payback, and activity-based costing without new machinery.